nueva página del texto (beta)
Español (pdf)
Artículo en XML
Referencias del artículo
Traducción automática
Enviar artículo por email
Citado por SciELO 
Similares en
SciELO 
Permalink
INTRODUCCIÓN
El territorio analizado se caracteriza por la heterogeneidad de su relieve, un régimen pluvial que supera los 2500 mm anuales, así como bruscos cambios en altitud y pendientes abruptas (Quesada-Román, 2016). El factor condicionante para las mayores inundaciones y procesos de ladera en el territorio son los ciclones tropicales del Mar Caribe y el océano Pacífico; entre los eventos más recordados por sus cuantiosos daños y muertes están los relacionados a los huracanes Joan en 1988 y César durante 1996, así como las tormentas tropicales Alma en 2008 y Nate en 2017 (Alfaro, Quesada-Román y Solano, 2010; LA RED, 2016).
El área de estudio se encuentra ubicada en su mayoría en la provincia de San José, en el cantón de Pérez Zeledón, aunque también abarca pequeñas porciones de la provincia de Puntarenas, Cartago y Limón. Dicho territorio se localiza entre las coordenadas geográficas: 9° 12’ 46” - 9° 35’ 03” Latitud Norte y 83° 28’ 47” - 83° 49’ 07” Longitud Oeste (Figura 1).
La cuenca alta del río General se encuentra en la zona centro-sur de Costa Rica en la vertiente pacífica, y está integrada por las subcuencas Buenavista, Chirripó Pacífico, Peñas Blancas y San Pedro. Todas ellas confluyen en el río más grande del país, el Térraba. Las elevaciones del área de estudio varían entre los 320 y 3820 msnm, desde el Cerro Chirripó hasta la llanura aluvial del río General.
Una de las características de este territorio es su variada litología, el sector septentrional, en particular la Cordillera de Talamanca, la cual está integrada de rocas volcánicas extrusivas (brechas, andesitas y depósitos de flujos piroclásticos) e intrusivas (granitoides con facies gabroides) del Mioceno. El piedemonte es resultado de la sobreposición de secuencias de fanglomerados que definen una morfología de abanico (Plio-Pleistoceno). El cauce del río General ocupa una llanura aluvial amplia y de fondo plano (Cuaternario). La Fila Brunqueña define el límite sur de la zona de estudio y se trata de un conjunto montañoso menor, de composición sedimentaria (lutitas y areniscas) del Oligoceno-Mioceno (Denyer y Alvarado, 2007).
La precipitación, en este territorio, condiciona la evolución del relieve y es de mayo a noviembre cuando su influencia es mayor. En este sentido, los datos promedio registrados en las estaciones meteorológicas cercanas (Repunta, Coopeagri y Chirripó), son de 2500 mm anuales, no obstante, en las montañas con altitudes mayores a los 1000 msnm, con facilidad superan los 4000 mm de lluvia. El histograma que resulta de estos registros es el típico de las zonas tropicales húmedas, con lluvias intensas relacionadas con sistemas de baja presión y ciclones tropicales entre agosto, septiembre y octubre (IMN, 2009).
Por otro lado, la afectación ciclónica es importante en esta zona, de manera directa se ha visto perjudicada por el huracán Joan en 1988 y César en 1996, las tormentas tropicales Alma en 2008 y Nate en 2017, eventos que generaron cuantiosos daños, e incluso muertes (Quesada-Román y Zamorano-Orozco, 2018b). Estos episodios extraordinarios de precipitación son los responsables de las inundaciones y de numerosos movimientos en masa. A esto se suma el carácter sismogenerador asociado a las fallas activas de Pangolín, Buenavista y División, todas ellas con historial sísmico reciente y a las que se les atribuye el devastador sismo de Buenavista en 1983 (Mora y Peraldo, 2011).
MATERIALES Y MÉTODOS
El trabajo tiene como punto de partida el análisis morfométrico, con el objetivo de establecer la intensidad de procesos exógenos y su vínculo con la dinámica del relieve, a partir del cálculo y análisis matemático de las formas terrestres (Goudie, 2004). Además, se define como un importante componente del análisis del terreno y el modelado de la superficie terrestre (Huggett, 2007). Este tipo de mapas son muy útiles en estudios de evolución de laderas, evaluación de erosión, así como de riesgos (Peña-Monné, 1997).
Los métodos morfométricos utilizados siguen los criterios de Simonov (1985), Lugo (1988) y Zamorano (1990). Los parámetros elegidos en este análisis son: altimetría, inclinación del terreno, densidad de la disección, profundidad de la disección, energía del relieve, erosión potencial y erosión total. La interpretación integral de este tipo de cartografía, permitió entender los mecanismos de ocurrencia de inundaciones y procesos de ladera, y, del mismo modo, zonificar los terrenos que favorecen su existencia.
Para establecer la susceptibilidad a procesos de ladera e inundaciones se sobrepusieron, a través de un Sistema de Información Geográfica (ArcGIS 10.3), los siguientes mapas morfométricos: densidad de la disección, profundidad de la disección, energía del relieve, erosión potencial y erosión total. En cada uno de ellos los datos se agruparon en tres rangos (altos, medio y bajos), proceder que permitió el análisis del cruce entre variables. De esta manera, en la medida que los valores bajos coinciden en uno o varios sectores de la cartografía síntesis, se hacen evidentes los terrenos subhorizontales o deprimidos y, por tanto, es posible deducir la dinámica asociada, inundaciones en este caso.
Por otro lado, los polígonos que representan valores altos ponen de manifiesto los terrenos montañosos, y entre más cruces de información coincidan, la morfología es más heterogénea, las laderas son más abruptas y la dinámica que se deduce para estas áreas es la gravitacional. Los datos intermedios, manifiestan una aparente estabilidad del terreno (Figura 2). Este tipo de información permite conocer la espacialidad de los procesos moderadores a partir de una zonificación precisa, lo que facilita una planeación objetiva del trabajo campo y la elaboración de la cartografía geomorfológica detallada.
Figura 2 Procedimiento utilizado para obtener la zonificación de procesos de ladera e inundaciones en la cuenca alta del río General por métodos indirectos.
Los territorios en donde coincidían hasta dos indicadores morfométricos se consideraron áreas potenciales; las zonas donde llegaban a sobreponerse tres índices se clasificaron de ocurrencia frecuente y los espacios catalogados como de máxima ocurrencia fueron aquellos donde se hacían presentes cuatro o cinco variables. Estos criterios se aplicaron tanto para procesos de ladera como a inundaciones.
Por último, fue posible comparar las áreas obtenidas como de ocurrencia de procesos gravitacionales e inundaciones con los eventos registrados en la base de datos de desastres DesInventar para Costa Rica entre 1970 y 2015 (LA RED, 2016).
ANÁLISIS MORFOMÉTRICO
Altimetría
Este documento simplifica la información topográfica, muestra la espacialidad y características de los contrastes altitudinales. Cuando esta información se sobrepone con el mapa de inclinación del terreno y se analiza junto con información geológica, es posible establecer unidades geomorfológicas; cada una de ellas define territorios similares en cuanto a génesis, morfología y dinámica. De este proceso resultaron las siguientes regiones morfológicas: zona de montaña (Cordillera de Talamanca), sierras menores (Fila Brunqueña), llanura aluvial, rampa acumulativa proluvial y estribaciones montañosas cubiertas por detritos proluviales (Figura 3).
La zona de montaña se presenta intensamente modelada, y en su expresión ha tenido que ver la existencia de potentes cortezas de intemperismo sobre laderas inclinadas y fracturadas, lo que favorece la existencia de procesos fluvio-gravitacionales. Las sierras menores sirven de límite a la planicie aluvial del río Pacuar y General, el carácter sedimentario, la falta de consolidación de los materiales, la existencia de laderas con fuerte inclinación, el régimen pluvial y cambios de uso de la tierra, los que convierten en susceptibles a estos terrenos a desarrollar procesos de ladera. Las estribaciones montañosas presentan una cima plana, y definen una morfología de mesetas con cubiertas aluviales sobre la base de las laderas. La rampa acumulativa proluvial es la transición entre el sector montañoso y la parte baja, su morfología está definida por la presencia de abanicos aluviales sobrepuestos y continuos a manera de glacis. Por último, en la llanura aluvial se reconocen el área de inundación estacional, la más baja, y la extraordinaria con expresión en gradería por la existencia de terrazas.
Inclinación del terreno
Esta variable tiene como objetivo reconocer terrenos potenciales a desarrollar procesos de ladera e inundaciones, al estar vinculados con la gravedad, la inclinación del terreno es de importancia para su desarrollo. De esta manera se pueden reconocer los rangos críticos y su concordancia con áreas estables e inestables. A su vez, se busca detectar la conexión entre la pendiente y la predisposición a procesos erosivos y acumulativos (Figura 4).
La zona de montaña presenta laderas con inclinaciones >26º y se caracterizan por ser superficies inestables, sobre todo en las cabeceras de ríos Chirripó Pacífico y Buenavista, en donde la pendiente en interacción con la precipitación, litología y estructuras disyuntivas, favorecen y magnifican los procesos de ladera. Por otro lado, los valores >15º son característicos de la llanura aluvial y la rampa acumulativa proluvial, en ellas, además de existir una dinámica de deposición sobre la erosiva, existen inundaciones en periodos tanto en la temporada de lluvias como en los periodos extraordinarios. Con respecto a las sierras menores, su inclinación es heterogénea y responde a una litología variada que favorece la incisión de valles.
Regiones morfológicas
El análisis de relieve en general se hace a partir de territorios comunes en cuanto a origen, dinámica, morfología, evolución y edad relativa. En este marco la zona de montaña (Cordillera de Talamanca) se considera una región escarpada en donde la altitud varía de 640 a 3820 msnm; su valor máximo corresponde a la cima del Cerro Chirripó (el punto más alto del país).
La zona de montaña se orienta de NW-SE, ocupa el sector centro-norte de la zona de estudio, y a nivel regional representan una pequeña porción del sector SE de la Cordillera de Talamanca, uno de los sistemas montañosos más importantes de Costa Rica. La morfología característica es de cimas agudas y vertientes heterogéneas en longitud, geometría y orientación. Otro de sus rasgos es la presencia de marcados desniveles que definen laderas escarpadas y profundos valles excavados en granitos miocénicos, sustrato que predomina en toda la estructura.
Por último, hay que mencionar que, por encima de los 3500 msnm, existen rasgos morfológicos de origen glaciar desarrollados durante el Pleistoceno y que continuaron hasta 12 000-10 000 A.P. (Veas, Quesada-Román, Hidalgo y Alfaro, 2018). En estos terrenos el relieve característico es de circos y valles glaciares, morrenas (de fondo, laterales, frontales y colgadas), superficies pulidas, estriadas, aborregadas o con morfología de lomo de ballena y tarns (lagos) (Quesada-Román y Zamorano-Orozco, 2018a).
Las sierras menores (Fila Brunqueña), de origen sedimentario, tienen su comienzo en periodos de sedimentación marino-costeros durante el Oligoceno-Mioceno (Denyer y Alvarado, 2007). Esta estructura es de poca altitud (320 a 960 msnm), no obstante, existen desniveles profundos y escarpados que forman parte de numerosos valles fluviales; de todos ellos, destacan el Pacuar y el General. Los ríos que tienen desarrollo en esta unidad presentan un patrón dendrítico muy denso, lo que hace evidente un proceso erosivo intenso que llega a magnificarse como respuesta al cambio del uso de la tierra.
La llanura aluvial representa la región más deprimida de la zona de estudio y se dispone entre las cotas de 320-1920 msnm e inicia en la zona de montaña, como un lecho amplio que integra los valles intermontanos, conforme el cauce se desplaza río abajo incrementa su amplitud hasta adoptar una morfología subhorizontal en donde la inclinación es mínima. Estas condiciones persisten en la porción sur de la zona de estudio y es la principal característica de los ríos Buenavista, Chirripó Pacífico, Pacuar y General.
La llanura aluvial se desarrolla en diferentes sustratos, tomando este hecho se ha dividido en tres segmentos; el primero, se ubica por encima de la cota de los 960 msnm, es estrecho y se emplaza sobre granitos con facies gabroides del Mioceno que se intercalan con tobas, areniscas, brechas e ignimbritas. El segundo sector se localiza entre los 600 y 960 msnm y en él predominan los detríticos, materiales que se asocian con extensos abanicos proluviales de edad Plio-Pleistoceno que sobreyacen a sedimentos continentales y de transición marina-costera del Cuaternario. Hay que mencionar que la susceptibilidad del material a la erosión fluvial favorece el desarrollo de las planicies fluviales en la horizontal (amplitud). La última región asociada a las planicies aluviales se dispone entre 320-600 msnm, en este sector afloran areniscas, conglomerados y lutitas de miocénicas (Denyer y Alvarado, 2007). Si bien el sustrato es competente a la erosión fluvial, su papel es secundario, ya que el desarrollo y disposición de las planicies aluviales está controlado por el trazo de los sistemas disyuntivos y la disposición de los ejes neotectónicos de fracturamiento.
La rampa acumulativa proluvial es una superficie poco inclinada, inicia a una altitud de 1600 msnm y su frontera más distal termina a los 320 msnm. Se ha desarrollado a partir de la coalición de abanicos aluviales de grandes dimensiones, en la mayoría de los casos no es posible reconocer de manera individual los límites de cada estructura, la orientación preferencial de todo el conjunto es NW-SE. Los detritos que la constituyen son relativamente recientes y poco consolidados; por un lado, se tienen conglomerados y areniscas fluvio-continentales del Plio-Pleistoceno como el sustrato principal en el sector NW, mientras que en el SE son sedimentos continentales y de transición marina-costera (Denyer y Alvarado, 2007). La existencia de una morfología conservada en cuatro de las grandes estructuras confirma, para ellas, una edad reciente.
Las estribaciones montañosas parcialmente cubiertas por detritos proluviales se presentan a manera de alveolos independientes y se localizan al centro-sur del área de estudio, se orientan de NW a SE y se constituyen de detritos del Mioceno y del Plio-Pleistoceno (Denyer y Alvarado, 2007). La altitud varía de 320 a 640 msnm; no obstante, algunas de ellas son > 700 msnm. Esta unidad tiene su origen en la sobreposición de capas de detritos vinculados con antiguos abanicos aluviales; esta estructura es el basamento donde sobreyace la rampa acumulativa proluvial. La morfología que predomina a lo largo de la frontera con la llanura aluvial es abrupta, a manera de frente montañoso, mientras, en sentido opuesto, hacia el norte, el contacto con la unidad III ocurre en forma de cubierta o manto, lo que favorece la existencia de una inclinación del terreno de 7° (Figura 5).
Densidad de la disección
A través de esta cartografía se hacen evidentes los territorios con mayor intensidad erosiva fluvial, y para ello se toma en cuenta la concentración de cauces en un área predeterminada. El primer paso en la elaboración de este mapa consistió en marcar a detalle la red fluvial sobre un mapa topográfico (1:25 000), y después se midió la longitud de los cauces con un curvímetro convencional utilizando una malla de captura equivalente a 1 km2. Los valores obtenidos en cada superficie fueron interpolados usando el método IDW (Inverse Distance Weighting) de ArcGIS 10.3. De esta manera se obtuvieron siete mapas, se presenta el que refleja la mejor distribución espacial de acuerdo con las características del territorio (Figura 6).
Las valores altos (7 a 11 km/km2) se localizan en la zona de montaña, en este territorio existe una morfología heterogénea en cuanto a longitud, orientación, geometría e inclinación del terreno. Estas características morfológicas, sumadas a una importante precipitación, favorecen la existencia de densas redes fluviales que definen patrones detríticos muy densos y valles angostos y profundos. Los valores, que varían entre 5 y 7 km/km2, tienen una amplia distribución en las laderas montañosas medias y bajas, así como en las rampas periféricas a la planicie fluvial e incluso en ella, siguiendo la inclinación regional del terreno. Los valores <5 km/km2 hacen evidentes terrenos ligeramente inclinados, o subhorizontales, en donde la dinámica que predomina es acumulativa sobre la erosiva. Estos índices también son característicos de cimas interfluviales amplias, en donde la morfología tiene características ya mencionadas.
Profundidad de la disección
Por medio de esta variable se conocen las zonas donde la erosión ha incidido con mayor intensidad en el relieve (valles profundos), el corte vertical de los ríos está en función del tipo de roca, la existencia de fallas, la inclinación del terreno y de la precipitación (Lugo, 1988).
La información se obtuvo al medir la distancia entre el talweg y la ruptura de pendiente más más cercana al borde inferior del valle. La lectura se realizó utilizando la base topográfica y la malla de muestreo anterior. La interpolación se realizó por el método Natural Neighbor (ArcGIS 10) y el resultado se presenta en la Figura 7.
En general, la profundidad de la disección está condicionada por la competencia del material a la erosión, la inclinación del terreno, la presencia de estructuras de debilidad (fallas), rupturas de pendiente y una precipitación importante.
Los valles más profundos (> 60 m) han sido excavados en sustratos granitoides y gabroides de la zona de montaña, en donde la fuerte inclinación del terreno y la existencia de densas redes de fallas ha favorecido la erosión vertical de los ríos. Esta dinámica es característica de cabeceras fluviales y de terrenos cercanos que integran la vertiente pacífica, en ambos casos el área que ocupan es poco significativa (Figura 6).
Los valles con profundidades entre 20-60 m se presentan en la mayor parte del relieve montañoso y se excluyen rampas acumulativas y la llanura aluvial, que, por sus características morfológicas, no favorecen el desarrollo de cauces profundos. Los valores <20 m son característicos de terrenos ligeramente inclinados o subhorizontales, como las rampas acumulativas, las llanuras de inundación y los terrenos cumbrales. En todos ellos la dinámica que predomina es deposicional.
Energía del relieve
Este mapa expresa la intensidad relativa de la actividad endógena en relación con la exógena; por tanto, valores altos de energía pueden asociarse a zonas de mayor actividad tectónica como acumulación de materiales volcánicos o movimientos tectónicos de ascenso asociado a la actividad de las fallas, de esta manera, es posible conocer el potencial erosivo al que puede estar sujeto el relieve (Lugo, 1988). Los datos se obtienen de la resta del valor altitudinal mínimo al máximo, en una superficie de un km2 a escala 1: 25,000. De esta manera, se obtiene una matriz de puntos que fue interpolada por el método Kriging (ArcGIS 10.3).
Los valores <60 m se localizan en terrenos donde los contrastes altitudinales no son representativos, como en la llanura aluvial, las rampas acumulativas y sectores de las sierras menores, en estos relieves la dinámica acumulativa esta sobre la erosiva. El rango de energía del relieve de mayor distribución es el de 60-400 m, se presenta como una zona de transición, entre las planicies aluviales y las cimas más altas (Figura 8). La dinámica que caracteriza este territorio es mixta (deposición-erosión), responde a la existencia de un relieve de contrastes altitudinales, con vertientes heterogéneas en longitud, geometría y, alternancia de afloramientos rocosos y detríticos.
Los desniveles de 400 a 900 m se restringen a la zona cumbral y a las laderas altas de la vertiente pacífica, sector donde se localizan las cabeceras más importantes de los ríos de la región. Al interior de estas formas existen numerosos saltos de cabecera activos (erosión remontante), que han desarrollado en conjunto una morfología de desniveles a manera de gradería. Las precipitaciones intensas son comunes en este sector (ciclónicas) y a ellas se asocia la existencia de valles profundos y laderas escarpadas.
Erosión potencial
El mapa que se presenta muestra el estado actual de la disección, así como los terrenos susceptibles a ser erosionados. Su elaboración inició con el cálculo de la inclinación de la ladera más larga en un área de 1 km2, y multiplicar este valor por el índice de densidad de la disección correspondiente (Lugo, 1988). Los valores resultantes fueron georreferenciados con el fin de interpolarlos mediante el método Natural Neighbor (ArcGIS 10.3).
Los valores que se obtuvieron fueron de 0 a 7.2 km/km2, y con el fin de facilitar su análisis se establecieron tres rangos: 4.4-7.2 km/km2, 1.5-4.3 km/km2 y 0-1.4 km/km2, y se analizan tomando en cuenta las unidades morfológicas (Figura 9).
Los valores <1.4 km/km2 son representativos de la llanura aluvial y en los abanicos aluviales y al SE, en algunos sectores de la rampa acumulativa proluvial. Estos relieves se caracterizan por presentar una inclinación que no favorece el desarrollo de cauces y, por lo tanto, la dinámica fluvio-gravitacional; no obstante, estos terrenos están constituidos por materiales competentes al transporte (detritos de compactación media o baja).
Los valores entre 1.5 y 4.3 km/km2, además de mostrar una dinámica mixta entre la deposición y la erosión, se localizan en buena medida en las laderas que descienden de la zona de montaña hacia la rampa acumulativa proluvial y en algunas divisorias de aguas en la montaña. Este territorio es heterogéneo en morfología, inclinación, alternancia de afloramientos rocosos y detríticos. Estas características, sumadas a la existencia de precipitaciones intensas, en particular las ciclónicas, hacen a estos terrenos vulnerables a los procesos fluvio-gravitacional.
Los valores entre 4.4 y 7.2 km/km2 están presentes en los terrenos más abruptos, en las cimas de las montañas del NW y en las laderas altas que integran la vertiente pacífica. Las condiciones morfológicas y de precipitación son similares a las mencionadas en el párrafo anterior, excepto la pendiente que, al ser mayor en esta zona, favorece la intensidad y magnitud de procesos de ladera y acarreo fluvial hacia las partes bajas.
Erosión total
Por medio de este mapa se pueden inferir aquellos terrenos que han tenido más tiempo expuestos al intemperismo y a la posterior actividad erosiva fluvial (gravitacional o incluso glaciar). El punto de partida es la expresión, la densidad de curvas de nivel y la medición de la longitud de cada una de ellas usando la malla inicial de muestreo de 1 km2, y a través de la herramienta Line Density, los datos obtenidos se interpolaron través de Natural Neighbor (ArcGIS 10.3).
El mapa de erosión total hace evidentes los terrenos en donde los procesos erosivos han sido persistentes e intensos, de igual manera, es posible reconocer los relieves recientes asociados a la deposición (llanuras aluviales y rampas acumulativas). Los datos que se obtuvieron fueron entre 0 y 45 km/km2 y con ellos se integraron los siguientes rangos: 0-13 km/km2, 14-22 km/km2 y 23-45 km/km2 (Figura 10).
Los índices <11 km/km2 se localizan en relieves de formación reciente y de origen exógeno acumulativo (abanicos y llanuras aluviales). Los valores entre 12 y 25 km/km2 a manera de una franja de transición entre la zona de montaña y los relieves acumulativos. El sustrato donde se emplazan estos valores es de ígneo intrusivo (granitoides), cuya morfología se ha suavizado por la existencia y desarrollo de cortezas de intemperismo. El último rango (26-45 km/km2) se presenta en la zona cumbral y en las laderas altas de la vertiente pacífica y caribeña, en donde predominan los terrenos rocosos abruptos, lo que dio como resultado una morfología variada. Hay que mencionar que las cubiertas detríticas (cortezas de intemperismo) no son representativas, debido a que fueron removidas por procesos fluviales, gravitacionales y glaciares.
RESULTADOS Y DISCUSIÓN
Zonificación de áreas susceptibles a procesos de ladera e inundaciones
Los métodos morfométricos son sencillos de elaborar, económicos y con resultados en poco tiempo, característica que cobra relevancia cuando se tiene como objetivo identificar territorios susceptibles a desarrollar peligros geomorfológicos, como la erosión, la colmatación, la remoción en masa y las inundaciones.
La zonificación que se presenta resulta de la integración espacial de los valores morfométricos (densidad y profundidad de la disección, energía del relieve, erosión potencial y total). El proceso de elaboración consiste en sobreponer los índices numéricos más bajos y altos de cada mapa morfométrico para obtener uno nuevo usando la herramienta Union (ArcGIS 10.3). El mapa que resultó se analizó con la topografía a fin de revisar los resultados con el relieve, en ese proceso fue determinante el conocimiento adquirido en los trabajos de campo y la interpretación de fotografías aéreas (Figura 11).
El análisis espacial de los datos morfométricos, en el marco de las regiones morfológicas, permitió elaborar la zonificación de procesos de ladera e inundaciones; este método se enriquece con el trabajo de campo e información geomorfológica.
Asociación con los procesos de ladera
La zonificación de los procesos de ladera está muy bien definida, son característicos del relieve montañoso de la Cordillera de Talamanca, unidades que están presentes en las cabeceras de los principales afluentes del río General; todas ellas tienen un arreglo que tiende a ser concéntrico. Las regiones asociadas con procesos de ladera contabilizan un área de 913 km2, lo que representa un 58.7% del territorio total del área de estudio (1560.5 km2); estas zonas se localizan en laderas de valle tanto de la Cordillera de Talamanca como de la Fila Brunqueña. Cabe resaltar que, a partir de un análisis estadístico de frecuencias, se logró establecer que la única variable que se repitió en más ocasiones fue la erosión total.
Los territorios en donde coincidían hasta dos valores morfométricos altos se consideraron áreas potenciales; las zonas donde llegaban a sobreponerse tres índices se clasificaron de ocurrencia frecuente y los espacios catalogados como de máxima ocurrencia, fueron aquellos donde se hacían presentes cuatro o cinco variables (Quesada-Román et al., 2018).
La zona de máxima ocurrencia coincide con los terrenos cumbrales y la porción más alta e inclinada de las laderas de montaña. Los valles característicos de estos territorios son profundos con morfología de cañón; en conjunto integran las cuencas de captación de los ríos más importantes de la región (Buenavista y Chirripó Pacífico), y se extienden a lo largo de 184.3 km2, lo que equivale al 20% de los terrenos asociados a procesos de ladera.
La litología de este sector es de granito y rocas volcánicas del Mioceno (en ambos casos), fracturadas por dos sistemas disyuntivos orientados de NW-SE y de N-S. Hay que mencionar que las condiciones de humedad y temperatura que existen en la zona de estudio han dejado su huella, ambos sustratos se presentan alterados y en la mayoría han desarrollado cortezas de intemperismo de espesores variables que afloran de manera irregular o aislada.
Las condiciones naturales de este territorio, considerado como el más susceptible a desarrollar procesos de ladera, son evidenciadas por el análisis numérico; en estas zonas coinciden los valores de cuatro o cinco variantes morfométricas. Por último, hay que tomar en cuenta que la lluvia (estacional o extraordinaria), los sismos, los cambios de uso de la tierra, la construcción de obras de infraestructura y vivienda, pueden catalizar la dinámica de laderas.
La zona de ocurrencia frecuente se presenta en dos sectores; el primero, se dispone bordeando la zona anterior y si bien presenta las mismas características morfológicas, litológicas, estructurales y climáticas; la segunda, en las Sierras menores, en frentes montañosos que dan al río Pedregoso. La diferencia radica en que en estas áreas coinciden tres de las cinco variantes morfométricas. Por su arreglo espacial pueden considerarse como una superficie de transición integrada por laderas de montaña medias. Representan un área de 191.1 km2, equivalente al 21.2% de los terrenos asociados a procesos de ladera. La mayor parte de la Fila Brunqueña está considerada en esta categoría, a excepción de su extremo SE. Las condiciones naturales se mantienen; es decir, son similares a las ya descritas, solo que el granito está ausente.
Los territorios potenciales se asocian también con relieves montañosos y, en particular, con las laderas bajas de la Cordillera de Talamanca y el sector SE de la Fila Brunqueña. Estas regiones se extienden a lo largo de 537.6 km2, lo que equivale al 58.6% de los terrenos asociados a procesos de ladera. Si bien existen todas las condiciones morfológicas, geológicas y estructurales para que se desarrollen procesos de ladera, para su ocurrencia hace falta un agente catalizador que inicie el proceso (natural o antrópico). Por último, hay que mencionar que en estas áreas coincide de una a dos variantes morfométricas.
Los terrenos donde no coincidió ningún índice morfométrico son casos aislados y en esta fase del análisis, fueron considerados como zonas de estabilidad aparente. Su extensión espacial es de 88 km2, lo que representa un 5.6% del área de estudio.
Asociación con inundaciones
Este aspecto fue interpretado a partir del análisis de la espacialidad de los valores morfométricos bajos, mismos que se disponen en el centro-sur de la zona de estudio. Este sector incluye la porción más deprimida (llanura aluvial) y territorios adyacentes (abanicos aluviales). Los relieves involucrados son homogéneos en morfología, inclinación y litología. Estas características determinan su dinámica; cuando llegan a presentarse lluvias intensas, los volúmenes de agua se distribuyen por toda la porción baja de manera turbulenta en la mayoría de los casos.
Las zonas asociadas con inundaciones tienen una extensión de 555 km2, lo que representa un 35.5% del total de la cuenca alta del río General. A partir de un análisis estadístico de frecuencias se estableció que las variables que se repitieron en más ocasiones fueron la profundidad de la disección y erosión total.
Las zonas en donde coinciden uno o dos valores morfométricos bajos se consideraron como áreas potenciales; las regiones donde llegaban a sobreponerse tres variantes se clasificaron como de ocurrencia frecuente y los sitios catalogados como de máxima ocurrencia fueron aquellos donde se hacían presentes cuatro o cinco variables.
Los territorios a presentar inundaciones con una máxima ocurrencia ocupan la mayor parte de la rampa acumulativa que está formada por amplios abanicos aluviales (Quesada-Román, 2017). Estos terrenos suman un área de 93.3 km2, lo que equivale a un 16.8% de las zonas asociadas a inundaciones en la cuenca alta del río General. La morfología de estas estructuras permite la circulación del agua en sus sectores convexos (porción central), mientras en los extremos o bordes laterales la geometría se invierte y tiende a ser cóncava; por tanto, el agua se concentra y retarda su circulación, lo que favorece la inundación. Esta es la situación de las llanuras aluviales de los ríos Pedregoso y General.
Las zonas de ocurrencia frecuente se localizan al norte de la rampa acumulativa a manera de una franja continua e irregular. Esta disposición refleja un control morfológico de los abanicos sobre esta área; es decir, la geometría ligeramente convexa de estas estructuras funciona como obstáculo y retarda la circulación de las corrientes fluviales. Esta condición cobrará importancia en la temporada de lluvias extraordinarias, donde este factor favorecerá la inundación. Estas regiones se extienden a lo largo de 69.3 km2, lo que equivale al 12.4% de los terrenos asociados a inundaciones.
El área potencial se refiere a aquellas zonas susceptibles a inundación desde un punto de vista geomorfológico; sin embargo, debido a que se encuentran ligeramente más altas en relación con la zona de máxima ocurrencia, deben existir o coincidir con condiciones extraordinarias para favorecer la dinámica, por ejemplo, temporada de lluvias intensas y ciclón o la obstrucción de un río por un proceso de ladera y rompimiento repentino de la represa natural. Representan un área de 392.4 km2, equivalente al 70.6% de los terrenos asociados a inundaciones.
VALIDACIÓN DE LA ZONIFICACIÓN MORFOMÉTRICA CON LOS PROCESOS DE LADERA E INUNDACIONES OCURRIDOS ENTRE 1970-2015
Los valores fueron revisados en campo para delimitar aún mejor las áreas de afectación real, con las que la morfometría podría malinterpretar la información. Al contrastar las áreas resultantes con estudios de los impactos en el área de estudio realizados por Peraldo (2004) y los impactos reportados en la base de datos DesInventar (LA RED, 2016), se confirma que las zonas de mayor susceptibilidad, tanto en procesos de ladera como inundaciones, coinciden con eventos reales suscitados (Tabla 1).
Tabla 1 Área en km2 de cada tipo de susceptibilidad por peligro geomorfológico.
| Tipo de susceptibilidad/Peligro geomorfológico | Procesos de ladera | Inundaciones | Total |
| Máxima ocurrencia | 184.3 | 93.3 | 277.6 |
| Ocurrencia frecuente | 191.1 | 69.3 | 260.4 |
| Área potencial | 541.6 | 392.4 | 934 |
| Estabilidad aparente | 88 | ||
| Total | 917 | 555 | 1560 |
La base de datos DesInventar está estructurada por la fecha de registro, el tipo de evento, la localización a la mayor escala de detalle (municipio o incluso distrito), el lugar, la fuente, las observaciones del evento. Por otro lado, registra el número de muertos, desaparecidos, heridos, víctimas, afectados, evacuados, relocalizados, casas destruidas, casas afectadas, hectáreas de cultivos y plantaciones forestales afectadas, rutas de comunicación afectadas, centros de educación afectados, centros de salud afectados, cabezas de ganado afectadas, pérdidas en moneda local, pérdidas económicas en dólares, otras pérdidas, sectores públicos afectados (transporte, comunicaciones, instalaciones de organizaciones de ayuda, acueductos, alcantarillado, educación, energía, industria, salud, entre otros), duración del evento, magnitud del evento, tipo de causa y observaciones de la causa (LA RED, 2016). Para Costa Rica, los reportes son realizados mediante la utilización de notas periodísticas, donde destacan las siguientes fuentes: periódico La Nación, informes de la Comisión Nacional de Prevención de Riesgos y Atención de Emergencias, entre otros. Los datos abarcan de 1970 a 2015 (LA RED, 2016).
El mapa de zonificación de procesos de ladera e inundaciones fue contrapuesto con la base de datos de desastres DesInventar. Los poblados como Buenavista, Palmital, Rivas, Guadalupe, San Gerardo y Herradura son afectados por movimientos en masa que han sido reportados e incluso cartografiados en estudios previos, y que, a su vez, se exhiben en el mapa resultado de este análisis. Extensos sectores de la Carretera Interamericana entre Siberia y San Isidro con frecuencia presentan procesos de ladera como desprendimientos.
Las comunidades, como Rivas, Peñas Blancas, Hermosa, General Viejo, Repunta y el centro urbano de San Isidro, han sido afectadas por efectos directos e indirectos de los ciclones tropicales con fuertes inundaciones, y dichos poblados resultan dentro del análisis zonas de alta susceptibilidad (Campos-Durán y Quesada-Román, 2017). Amplias zonas relacionadas con la llanura aluvial del río General y sus poblados aledaños se ven afectadas de manera frecuente por las inundaciones y desbordamientos que afectan áreas de cultivo e infraestructura pública, lo que causa importantes pérdidas económicas en el municipio de Pérez Zeledón.
Del total de reportes para el área de estudio (209), 128 se relacionan con inundaciones y 81 con procesos de ladera. Estos datos fueron ubicados con la mayor precisión posible según cada reporte en el territorio, después fueron sobrepuestos en el mapa de zonificación de la susceptibilidad a procesos de ladera e inundaciones (Figura 12).
Figura 12 El mapa de zonificación e impactos de procesos de ladera e inundaciones cotejado con los impactos suscitados por estos procesos entre 1970 y 2015, según la base de datos DesInventar
La ubicación de los reportes de DesInventar en el área de estudio se relacionan con el resultado de la zonificación de los procesos de ladera e inundación, y se logran calcular los porcentajes de validación según el número de datos que coinciden con las áreas susceptibles determinadas. El análisis determina que el 85.15% de los reportes coinciden con las zonas propensas a inundación mientras que los procesos de ladera tuvieron una relación del 76.54% con los sitios propensos a estos peligros geomorfológicos.
CONCLUSIONES
Mediante el análisis morfométrico se analizaron siete variables: altimetría, inclinación del terreno, densidad de la disección, profundidad de la disección, energía del relieve, erosión potencial y erosión total. A partir de la altimetría se lograron delimitar las regiones morfológicas de la cuenca alta del río General: I) zona de montaña (Cordillera de Talamanca); II) llanura aluvial; III) rampa acumulativa (abanicos aluviales); IV) sierras menores y V) estribaciones montañosas cubiertas por detritos aluviales.
La zonificación de procesos de ladera e inundaciones resultó de la integración espacial de las siete variables morfométricas. Se lograron definir para cada peligro geomorfológico tres categorías de susceptibilidad: máxima ocurrencia, ocurrencia frecuente y área potencial. El mapa de zonificación de procesos de ladera e inundaciones fue corroborado con la base de datos de desastres DesInventar. El 85.15% de los reportes coinciden con las zonas propensas a inundación mientras que el 76.54% con las áreas susceptibles a procesos de ladera. Este método puede ser replicado con facilidad dado que solamente se utilizan las curvas de nivel y ríos, las mediciones son sencillas y pueden ser realizadas a partir de un sistema de información geográfica en su totalidad. Esta condición agiliza los tiempos de digitalización de cauces fluviales, el procesamiento de la información, la sistematización del procedimiento y la practicidad tanto para hacer distintas mediciones como para la toma de datos puntuales de las isohipsas y fluvios.
Por medio de este procedimiento se pueden hacer valoraciones de base para el ordenamiento territorial de cuencas hidrográficas y unidades político-administrativas a escalas locales y regionales. La geomorfometría es una herramienta práctica para la mitigación del riesgo a desastres mediante la determinación de áreas susceptibles a inundaciones y procesos de ladera en general. No obstante, requiere una corrección con las formas del relieve para precisar las áreas a afectarse a partir de una cartografía de mayor precisión (escalas 1:25 000 o superiores), lo que solamente puede ser realizado por medio de un mapa geomorfológico a detalle.
