Análisis de fragmentación del bosque

Este estudio se basó en el análisis del uso de la tierra, la caracterización de la cobertura vegetal por medio de un análisis estadístico de los parches de vegetación y una investigación sobre la red de drenaje mediante un tratamiento morfométrico y el establecimiento de la relación entre la densidad del sistema hídrico y los parches de bosque. Para el análisis de fragmentación, se delimitaron cuatro sectores en la subcuenca (microcuencas), a saber: Parque Nacional Juan Castro Blanco -en adelante, PNJCB-, Pueblo Nuevo, Tapezco y Las Brisas, con el fin de comparar resultados entre diferentes zonas, demarcadas por su ubicación en la parte alta o baja, así como por la intensidad de las actividades.

A partir de la fotointerpretación de imágenes de la Misión Costa Rica Airborne Research and Technology Aplications del Programa de Investigaciones Aerotransportadas y Sensores Remotos y la Administración Nacional de Aeronáutica y el Espacio de Estados Unidos (^{PRIAS y NASA 2005}), a escala 1:5000, se obtuvo una capa preliminar de usos de la tierra, cuyos límites y contenido fueron corroborados y actualizados en campo para obtener un mapa de uso de la tierra en el 2012. Este fue el insumo principal para los análisis posteriores.

De las once categorías iniciales de uso de la tierra, se seleccionaron cuatro para determinar la fragmentación y la formulación de enlaces de vegetación: bosque, vegetación ribereña, bosque en regeneración y pastos arbolados. Estos tipos de cobertura fueron escogidos porque favorecen la conectividad entre parches de bosque y contribuyen a diversificar la estructura de la matriz y el paisaje en general (^{Herrera y García 2009}). Además de estas coberturas, también se consideraron las cercas vivas y las cortinas rompevientos, porque, según sea su disposición en el espacio, ofrecen una serie de beneficios para la estructura y configuración del paisaje, así como para la protección de fuentes de agua (^{Harvey et ál. 2008}). A partir de una nueva capa, se realizó con dichas coberturas un análisis estadístico mediante el uso la extensión Patch Analyst, del software ArcGIS 10.0. El estudio se realizó en cada sector de la subcuenca con fines comparativos. Las métricas empleadas fueron de área, forma y contagio (^{MacGarigal y Marks 1995}).

Las métricas de área correspondieron a: 1) tamaño promedio de fragmentos, que es el área promedio de los fragmentos de un mismo tipo contenidos en un paisaje; 2) número de fragmentos, que se refiere a la cantidad de parches correspondientes a un tipo de cobertura específico. La métrica de forma consistió en el índice que mide la complejidad de la configuración externa de los parches, comparándolos con una forma estándar correspondiente a un círculo. Cuando la forma se asemeja a este, el índice es menor y crece conforme el parche se deforma. Finalmente, la métrica de contagio consistió en la distancia promedio al vecino más cercano, que constituye la distancia de un parche a otro del mismo tipo basada en la distancia de borde a borde.

Para conocer en términos numéricos la continuidad del bosque en la subcuenca, se determinó el índice de Vogelmann, el cual indica qué tan continuo es el bosque tomando como parámetros el área y el perímetro. En este caso, los valores altos implican mayor continuidad del bosque, mientras que los valores bajos reflejan mayor fragmentación y discontinuidad de los parches (^{Vogelmann 1995}). El índice se calculó mediante la ecuación 1:

Donde:

FCI = índice de continuidad del bosque (Forest continuity index, por sus siglas en inglés).

A = área total de parches de bosque del paisaje.

P = perímetro total de parches de bosque del paisaje.

En la caracterización de la red hídrica se determinaron los órdenes de cauces. Para establecer la forma y la densidad del drenaje se usó la ley de Horton (^{Chorley 1995}). Además, para relacionar la presencia de parches de bosque con la red de drenaje de la subcuenca y estimar así el grado de cobertura de los cauces, se determinó la densidad de parches por área, según el número de parches y el área de estos, así como la densidad de drenaje.

Formulación de enlaces de vegetación

Los análisis de conectividad se realizaron en tres etapas. La primera consistió en un análisis de los vacíos de vegetación a lo largo de la red hídrica, identificando su localización y los usos de la tierra en estas zonas, para determinar sitios que requieren mayor atención por los posibles impactos de los usos de la tierra en el recurso hídrico (^{Njue et ál. 2016}; ^{Uriarte et ál. 2011}). En la segunda se realizó un análisis del dinamismo en el cambio de uso de la tierra y se determinaron zonas críticas de cambio en la subcuenca, ya que los usos de la tierra y los cambios que estos presentan pueden tener implicaciones al considerar la gestión del agua (^{Uriarte et ál. 2011}). Finalmente, en la tercera etapa se analizó en forma integrada la conectividad en la subcuenca, a partir de los resultados previos.

La determinación de la red de conectividad se basó en el análisis de vacíos de vegetación a lo largo de los ríos, según la ley Forestal 7575 de Costa Rica, artículo 33, que establece las áreas de protección respectivas en los alrededores de nacientes y márgenes de ríos, quebradas y arroyos. La ley determina una franja de quince metros en zona rural, de diez metros en zona urbana si el terreno es plano y de cincuenta metros si es quebrado (^{República de Costa Rica 1996}). Se consideran terrenos quebrados u ondulados aquellos cuya pendiente supera el 15%, según el Decreto 23214-MAG-MIRENEM (^{Presidencia de la República de Costa Rica 1994}).

De acuerdo con estas áreas de protección, se determinó un buffer a lo largo de la red hídrica y alrededor de las nacientes de agua (cien metros). Para esto se elaboró un modelo digital de elevación a partir de curvas de nivel cada cincuenta metros, para determinar el grado de la pendiente y así definir el ancho del buffer, según la variación de la topografía a través de la red hídrica. Este proceso se realizó mediante la herramienta Buffer del módulo ArcMap de ArcGis 10.0. A la capa buffer generada se sobrepuso la capa de coberturas vegetales seleccionadas para el análisis de la fragmentación. El resultado fue un mapa con las zonas en donde existen vacíos de vegetación y que no permiten la continuidad del bosque en los ríos, es decir, dicho mapa posibilita identificar el grado de interrupción de la vegetación ribereña. En estos vacíos se determinó el tipo de cobertura presente en cada uno, para posteriormente evaluar las zonas donde sería más viable restaurar la cobertura vegetal.

El análisis del dinamismo en el cambio de uso de la tierra se realizó en tres periodos: 2005, 2008 y 2012, con ello se determinaron las zonas de mayor y menor cambio a través del tiempo. Se identificaron zonas críticas que requieren acciones prioritarias para mejorar la conectividad del bosque. Con el fin de realizar comparaciones estadísticas entre un año y otro, así como para analizar el cambio temporal del uso de la tierra en el paisaje, se utilizó el índice de diversidad de Shannon -en adelante, SDI (por sus siglas en inglés)-, el cual mide la heterogeneidad a partir de la diversidad de fragmentos mediante la ecuación 2:

SDI toma valores entre cero e infinito. Será igual a cero cuando una categoría ocupe todo el espacio y aumentará a medida que aumente heterogeneidad (^{Vila et ál. 2006}).

La propuesta de enlaces tuvo como eje principal la red hídrica y la continuidad del bosque a lo largo de ella, pero complementariamente se consideró la presencia de parches boscosos y las nacientes, ya que estas constituyen puntos terminales de la red de conectividad propuesta. Este análisis se realizó mediante la sobreposición de la capa correspondiente al buffer con las capas de nacientes y terrenos adquiridos por las ASADAS o alguna ONG para hacer reforestación.

Con el fin de proponer una red de conectividad integrada, se tomó como punto de partida las zonas de vegetación en la red hídrica definidas por el buffer. De allí, se partió hacia las fincas en regeneración y las nacientes que contaran con alguna cobertura vegetal o que requirieran ser cubiertas. Los criterios utilizados para trazar las rutas fueron: el uso de la tierra entre las zonas a conectar, así como la menor distancia entre estos puntos.

Los enlaces se diseñaron en dos etapas: la primera consistió en el análisis visual del mapa integrado de conectividad en la subcuenca y la sectorización definida, mientras que en la segunda etapa se realizó la sobreposición de las capas de buffer, terrenos adquiridos en regeneración, nacientes y coberturas vegetales con el software ArcMap de ArcGIS 10.0, a partir de lo cual se trazaron posibles rutas de conectividad. La combinación de estos criterios se logró mediante la herramienta Cost Path del paquete ArcGIS 10.0, que determina la ruta de menor coste entre un punto de destino y un origen, en relación con las unidades de coste. Los enlaces partieron de la red hídrica, teniendo como destino las coberturas seleccionadas y como unidades de coste el uso de la tierra. Se dio prioridad a los parches de bosque y bosque en regeneración como unidades de menor coste para diseñar los enlaces. El resultado de este proceso fue una red de conectividad que integra toda la subcuenca.

Fragmentación del bosque en la subcuenca del río Tapezco

La colonización agrícola del cantón Zarcero se realizó mayormente por unidades familiares que cultivaron su porción de tierra y comercializaron sus productos, lo que da indicios del patrón de uso de la tierra desde la década de los cuarenta. La riqueza de este cantón se basó en los bosques maderables, la tierra cultivable y la abundancia del recurso hídrico; sin embargo, la explotación maderera tuvo como consecuencia la deforestación en la zona desde su fundación en 1915 (^{Murillo Vargas 2009}).

Zarcero es considerado un municipio fundamentalmente agropecuario, donde el 90% de las actividades económicas a las que se dedica la población son de este tipo. La diversidad en el uso de la tierra y su división han propiciado un patrón de tenencia fragmentada. Las fincas de hortalizas van desde 0,5 ha hasta 2 ha de terreno, siendo generalmente fincas familiares. Por su parte, las fincas ganaderas oscilan entre las 8 y 12 ha en Las Brisas y Palmira. El patrón paisajístico característico es de pequeñas unidades agrícolas, principalmente hortícolas y fincas dedicadas a la producción de leche. La matriz está compuesta por pastos y cultivos, con parches de bosque dispersos, y usos urbanos generalmente a lo largo de las vías de comunicación (figura 2).

Fuente: vuelo LightHawk 2011, Instituto Nectandra.

Figura 2 a) campos de pastos y cultivos con escasa presencia de bosque en Pueblo Nuevo. b) parches de bosque a lo largo de los ríos y concentraciones urbanas siguiendo las carreteras en Tapezco. 

En la figura 3 y la tabla 1 se presenta la heterogeneidad de usos en la subcuenca, donde se observa un total de once categorías con el predominio del área de pastos (26,5%), seguido de la cobertura de bosque (26%), concentrado fundamentalmente en la parte alta de la subcuenca, y cultivos (19%). Los pastos arbolados representan un 14% y, por su parte, el bosque en regeneración tan solo un 5,2%.

Datos: comprobación de campo 2012, fotointerpretación imágenes de PRIAS y NASA 2005; ^{Ortiz-Malavassi 2008}.

Figura 3 Usos de la tierra, subcuenca del río Tapezco, 2012 

Tabla 1 Usos de la tierra en la subcuenca del río Tapezco Zarcero (2012) 

Datos: derivados de fotointerpretación de imágenes de PRIAS y NASA 2005.

La heterogeneidad del uso de la tierra en el paisaje origina una diversidad espacial del mosaico y crea contrastes entre parches de diferentes tipos. El patrón espacial que origina la heterogeneidad influye en diversos procesos y flujos de materia y energía que ocurren en el paisaje, tal como el flujo hídrico (^{Farina 2006}; ^{Uriarte et ál. 2011}). Por su parte, en virtud de la diversidad de usos y su disposición espacial en el paisaje, el nivel de contraste de vecindad entre cada uso tendrá mayor o menor significatividad para la continuidad del bosque y los procesos ecológicos asociados (^{Herrera 2011}).

Según un análisis por sector, en comparación con los demás sectores, el PNJCB presenta el mayor porcentaje de área cubierta por bosque (42,43%), seguido de pastos con árboles dispersos (17,55%), lo que resulta beneficioso para la conectividad de la parte alta de la subcuenca, adyacente a un área protegida que tiene entre sus objetivos la defensa del recurso hídrico en la región. Los sectores Pueblo Nuevo, Tapezco y Brisas, en la parte baja de la subcuenca, presentan un predominio de cultivos y pastos con presencia de parches de bosque aislado (figura 4) o con formas alargadas que siguen la red hídrica.

Fuente: Vuelo LightHawk 2011, Instituto Nectandra.

Figura 4 Parche boscoso de gran tamaño en medio de pastos y cultivos en el sector de Pueblo Nuevo. 

La cobertura boscosa está concentrada principalmente en el sector PNJCB, así como en las márgenes de los ríos en los sectores de Pueblo Nuevo, PNJCB, Tapezco y, más escasamente, en Las Brisas. Los pastos, por su parte, presentan un patrón de distribución disperso en la mayoría de la subcuenca a excepción del sector PNJCB. Los cultivos, de igual forma, se encuentran distribuidos en toda la subcuenca, pero tienden a formar pequeñas concentraciones.

En el sector de Tapezco, que corresponde a la parte baja de la subcuenca, se observa un claro predominio de cultivos y pastos que ocupan casi el 60% del área. De forma similar, el sector de Las Brisas está cubierto en un 70% por cultivos y pastos, y apenas en un 8% por bosque, lo que lo convierte en el sector más degradado en términos de recursos para la conectividad. En casos donde el paisaje natural se encuentra muy intervenido por el ser humano, las cercas vivas y rompevientos representan un recurso viable, en el corto plazo, para iniciar alguna estrategia de conectividad entre parches de bosque, ya que se ha demostrado que aquellas, al igual que los árboles dispersos en zonas de pastoreo, contribuyen con la conectividad en el paisaje, además de brindar aportes para la conservación del suelo, el agua y la diversidad biológica (^{Harvey et ál. 2008}).

El patrón de distribución de uso de la tierra refleja las estrategias tradicionales de conservación, en las que se protege las zonas altas de las cuencas y las zonas bajas se dedican mayoritariamente a la explotación de recursos. Esto genera diferentes tipos de matriz, más homogéneas en la parte superior y con más cantidad de hábitat boscoso disponible. Las partes bajas presentan mayor heterogeneidad de usos, usualmente con producción intensiva. Sin embargo, se debe considerar que no solo la cantidad de área boscosa resulta importante, sino que también lo es la configuración espacial de los fragmentos remanentes en la matriz. En el caso de Las Brisas y Tapezco, la vegetación ribereña remanente y las cercas vivas mejoran la estructura en un paisaje agrícola degradado, lo cual es asimismo de interés para la conservación, tal como señalan ^{Harvey et ál. (2006)} al estudiar los paisajes de pastos arbolados y fauna asociada en Nicaragua.

De acuerdo con los usos mayoritarios en la subcuenca, se define una matriz compuesta principalmente por pastos y cultivos, a esta se suman los restantes usos (pastos con árboles, bosque en regeneración, vegetación ribereña, rompevientos, uso urbano, invernaderos y plantación forestal), que corresponden al 30% del área, lo que genera una matriz muy diversa. Desde la aproximación que promueve las matrices agroecológicamente diversas como alternativa para la conservación en paisajes antropizados, este tipo de matriz resulta ventajosa, ya que, como apuntan ^{Gurrutxaga y Lozano (2006)} y ^{Herrera (2011)}, las matrices con estructura y configuración heterogénea favorecen, hasta cierto punto, la diversidad de hábitats potenciales, así como los flujos de materia y energía.

El proceso de fragmentación en un paisaje se manifiesta por medio de diferentes efectos, como el aumento en el número de parches, la disminución de su tamaño, el incremento del aislamiento entre aquellos y su forma (^{Fahrig 2003}). Dichos efectos introducen cambios en la configuración del patrón espacial del paisaje, que ocasionan a su vez cambios en la funcionalidad ecológica de este. Estos cambios ecológicos pueden tener efectos perdurables en la flora, fauna, suelo y en los recursos hídricos (^{Bennett 2004}). El cálculo de métricas de vegetación permite determinar y analizar el grado de fragmentación en la subcuenca y su comparación entre sectores.

Las métricas de coberturas vegetales se presentan en la tabla 2, donde se observa que en el sector Las Brisas hay 51 parches de bosque con un tamaño promedio de casi 3 ha. Este es el sector que posee los parches de menor tamaño, lo que se explica por la gran diversidad de usos del suelo que existe allí. En los sectores de Pueblo Nuevo y Tapezco el tamaño promedio es similar: de 2 y 1,9 ha, respectivamente; sin embargo, la cantidad de parches en Pueblo Nuevo es 86, mientras que en Tapezco es 30, esto refleja que la fragmentación es mucho mayor en el sector de Pueblo Nuevo.

Tabla 2 Métricas de coberturas vegetales por sector, subcuenca del río Tapezco 

Datos: derivados de fotointerpretación de PRIAS y NASA 2005, procesados mediante la herramienta Patch Analyst de ArcGIS 10.

En el sector del PNJCB existen 57 parches de bosque, siendo el segundo con mayor número de parches, por lo que se podría considerar uno de los más fragmentados. Sin embargo, presenta un parche de 291 ha, que tanto por su área como por su forma dendrítica hace que este sector sea el de mayor continuidad de bosque en la subcuenca, aunque el número de parches sea mayor que en otros sectores. Sin este fragmento, el tamaño promedio de los parches se reduciría notablemente, por lo que su presencia hace que el sector se convierta, en términos numéricos, en el menos fragmentado.

El mayor porcentaje de bosque en regeneración se encuentra en el sector del PNJCB, esto obedece a que las organizaciones comunitarias han adquirido progresivamente terrenos para reforestación en el área de amortiguamiento del PNJCB, con lo que se ha creado un efecto buffer. Aquí es importante resaltar que si bien el mantenimiento de parches grandes y complejos resulta primordial, también lo es considerar la complejidad estructural y la configuración de usos que rodea dichos parches, para disminuir el efecto borde (figura 5) y de contraste estructural, es decir, enfocarse en la calidad de la matriz circundante, cuando no es posible dedicar completamente los terrenos a la regeneración de la vegetación (^{Williams et ál. 2006}).

Fuente: Vuelo LightHawk 2011, Instituto Nectandra.

Figura 5 Parches de bosque continuos con efecto borde pronunciado. 

Aspectos como el uso de la tierra, las formas del terreno y la red fluvial condicionan la configuración externa y distribución de los parches, con lo que se generan patrones particulares (^{Meddens et ál. 2008}). La cobertura boscosa presenta, en la mayor parte del área de estudio, una configuración que sigue la sinuosidad de la red de drenaje. En el sector del PNJCB, además de esto, el bosque se presenta concentrado en la parte más alta de la subcuenca. Esto obedece, aparte de la cercanía al parque nacional, a las pendientes más pronunciadas en esa zona, lo que restringe los usos agrícolas y pecuarios.

Dado que la mayoría de parches siguen la forma de la red hídrica, el índice de forma para el bosque en la subcuenca es 2,26, por lo que se considera que los parches tienen forma irregular. Por razones ecológicas las formas irregulares y alargadas hacen que los parches sean más vulnerables a perturbaciones por efecto borde (^{Santos y Tellería 2006}); no obstante, esta configuración, desde la perspectiva del paisaje, resulta beneficiosa para la conectividad, debido a que se constituyen hábitats lineales a lo largo de los cauces, a partir de los cuales se puede establecer toda una red de conectividad hacia los parches más grandes y alejados del sistema de drenaje.

El aislamiento entre parches boscosos es otra consecuencia del proceso de fragmentación: a medida que avanza el proceso, la distancia entre parches aumenta y estos son separados por los nuevos usos, lo que interrumpe su continuidad (^{Guariguata y Kattan 2002}). En términos de distancia respecto del vecino más cercano, el sector del PNJCB es el que presenta mejores condiciones; tanto el promedio como la mediana son los menores valores para todas las coberturas con respecto a los demás sectores; por su parte, los sectores de Pueblo Nuevo y Tapezco presentan mayores promedios de distancia y mayores rangos, lo que significa que estos sectores son los que presentan menor continuidad del bosque.

El índice de continuidad del bosque de Vogelmann -en adelante, FCI (por sus siglas en inglés)- confirma lo descrito anteriormente: el sector del PNJCB, con un FCI de 4,02, es el que presenta mayor continuidad del bosque; por el contrario, el sector de Las Brisas es el menos continuo, con un FCI de 2,74. Por su parte, para los sectores de Pueblo Nuevo y Tapezco el FCI arrojó valores de 3,29 y 3,10, respectivamente. En estos dos sectores la intensidad y configuración de los usos agropecuarios ha impedido un patrón continuo del bosque, presentando una matriz con menor cantidad de hábitats remanentes y más distanciados. De acuerdo con ^{Fahrig (2003)}, dicho entorno circundante genera mayor vulnerabilidad para los parches relicto.

Respecto de la caracterización del sistema hidrográfico, se identificó que la subcuenca del río Tapezco presenta una forma alargada e irregular, donde el sistema de drenaje, según su forma, se cataloga como dendrítico. En cuanto al número de orden de cauces y según escala de trabajo, se considera la subcuenca de orden cuatro y el tipo de actividades socioeconómicas desarrolladas en la subcuenca hacen que esta se defina como agropecuaria, principalmente para la producción de hortalizas, así como para el desarrollo de la producción de leche. Estas actividades se extienden, en varios casos, hasta las márgenes de los ríos donde ha sido eliminada la vegetación de ribera, lo que aumenta el riesgo de contaminación por vertido de agroquímicos y efluentes de las explotaciones agropecuarias. Este riesgo aumenta en zonas donde la cobertura vegetal es escasa o inexistente, ya que la presencia de esta contribuiría a mejorar el flujo y calidad del agua (^{Allan 2004}; ^{Bruijnzeel 2004}; ^{Carvalho-Santos, Honrado y Hein 2014}).

La estrecha relación que la vegetación de ribera tiene con el estado del medio fluvial hace que las alteraciones de estos sistemas susciten graves problemas para el mantenimiento y conservación de ambos. La vegetación ribereña modera una gran cantidad de procesos ecológicos en los cursos fluviales, influye en los regímenes de temperatura y luz, y produce detritos orgánicos, como restos leñosos y de hojas, encauzamiento del agua y de los sedimentos, así como la estructuración del hábitat físico en diferentes escalas, proporcionando sustrato y hábitat para la actividad biológica (^{Ward et ál. 2002}). Debido a esta importante interacción, el análisis de la red de drenaje que comprendió la relación entre la densidad de parches, según número de parches y su área, además de la densidad de drenaje, permitió tener una noción clara de la cobertura de bosque existente con respecto a la red hídrica, como un aporte de base para futuros estudios a profundidad en materia de calidad de agua. Como resultado, se encontró que los parches de bosque obedecen generalmente a un patrón de distribución a lo largo de la red hídrica, lo que, además de los aportes mencionados para la conectividad estructural y ecológica, resulta muy beneficioso en términos de protección de los cauces. Estos cálculos se realizaron para cada sector de la subcuenca (tabla 3) y se representan en la figura 6.

Tabla 3 Relación parches de bosque-densidad de drenaje por sector, subcuenca del río Tapezco (2012) 

Datos: derivados de fotointerpretación de PRIAS y NASA 2005.

Datos: derivados de fotointerpretación de PRIAS y NASA 2005.

Figura 6 Relación distribución y densidad de parches de bosque-densidad de drenaje, subcuenca del río Tapezco (2012). 

El sector del PNJCB presenta menor densidad de parches por número y mayor densidad de parches por área, es decir, pocos parches de bosque, pero de gran tamaño y con mayor densidad de drenaje, debido a que corresponde a la parte alta de la subcuenca, donde la red es más intrincada. El sector de Las Brisas presenta un mayor número de parches y de menor tamaño respecto de la red hídrica, lo que expresa una menor continuidad y mayor fragmentación de las coberturas vegetales, debido a la presión de usos de la tierra en ese sector.

Se ha comprobado que la cobertura boscosa ofrece una serie de servicios hidrológicos en la escala más amplia del paisaje y en la red hidrográfica (^{Carvalho-Santos, Honrado y Hein 2014}), ya que la presencia no solamente de parches extensos de bosque, sino también de vegetación ribereña, árboles dispersos en pastos y cercas vivas, mediante el sistema de raíces y el dosel, contribuye a los procesos de infiltración de agua en el suelo y disminución de escorrentía, lo que la intercepción de humedad por las copas de los árboles (^{Ilstedt et ál. 2007}). La calidad del agua se ve favorecida, en términos de carga de sedimentos, debido también al mejoramiento de la estabilidad del suelo (^{Calder 2002}; ^{Lele 2009}).

De acuerdo con los resultados obtenidos de las métricas de fragmentación para cada sector, se tiene que el sector del PNJCB es el que presenta mejores condiciones en cuanto a continuidad, densidad y superficie de bosque. Por el contrario, al sector de Las Brisas corresponde al sector más afectado en términos de distancia entre fragmentos y tamaño más reducido, así como en su continuidad. Esto permite identificar una marcada diferencia entre la parte alta y baja de la subcuenca, en cuanto a la configuración espacial de los parches de bosque.

En las partes bajas el bosque se localiza principalmente en las orillas de los cauces, mientras que en las partes altas los parches, además de estar a lo largo de la red hídrica, también se presentan fuera de esta y con mayores dimensiones. Esto obedece a que en las partes bajas la topografía es menos abrupta y existen más y mejores vías de comunicación para el desarrollo de las actividades económicas, por lo que hay mayor presión en el uso de la tierra y son menores las oportunidades para adquirir terrenos que se dediquen a la regeneración.

En las partes altas la presión sobre el uso de la tierra es menor, debido a condiciones topográficas y a que los esfuerzos de conservación se han concentrado en esta sección de la subcuenca para dar continuidad a los ecosistemas boscosos del PNJCB.

Respecto de estos patrones de diferenciación de usos y coberturas ^{Cotler y Pineda López (2005)}, señalan que en las cuencas hidrográficas la identificación y el reconocimiento de las interconexiones existentes entre la parte alta y baja de aquellas resultan fundamentales, considerando las externalidades generadas por las actividades productivas de forma integrada en todo el sistema hidrográfico.

Formulación de enlaces de vegetación

Los diferentes usos de la tierra, su dinamismo en cuanto a los cambios temporales y su disposición en el espacio juegan un papel determinante en diferentes aspectos relativos a las propiedades y el funcionamiento de los sistemas fluviales (^{Njue et ál. 2016}; ^{Uriarte et ál. 2011}). Siguiendo estas consideraciones y en el caso de la subcuenca del río Tapezco, se encontró que la vegetación ribereña ha sufrido cambios pocos significativos en el periodo 2005-2012, presentando una leve mejoría en cuanto al porcentaje de vegetación faltante (tabla 4). Las áreas con más grandes vacíos de vegetación se ubican en los afluentes de los ríos principales (Tapezco, San Rafael, Catarata y Jilguero), que corresponden a zonas de mayor pendiente y donde se llevan a cabo actividades agrícolas y de pastoreo.

Tabla 4 Porcentaje de vacíos de vegetación, subcuenca del río Tapezco (2005, 2008 y 2012) 

Datos: derivados comprobación en campo, de fotointerpretación de PRIAS y NASA 2005.

Un análisis más detallado de dichos vacíos se realizó para cada sector de la subcuenca y se determinó el área de cobertura vegetal faltante a lo largo de la red hídrica y las nacientes. Los resultados se muestran en la tabla 5 y figura 7.

Tabla 5 Vegetación ribereña faltante por sector, según Ley forestal 7575, subcuenca del río Tapezco Zarcero, 2012 

Datos: derivados de comprobación en campo y fotointerpretación de PRIAS y NASA 2005.

Datos: derivados de fotointerpretación de PRIAS y NASA 2005, procesados mediante la herramienta Buffer de ArcGIS 10.0.

Figura 7 Vacíos de vegetación, según Ley forestal 7575, subcuenca del río Tapezco (2012). 

El sector de Pueblo Nuevo presenta el mayor porcentaje de vegetación ribereña faltante, debido a que existen zonas con pendientes mayores al 15% que deberían tener una banda ancha de cobertura vegetal según la Ley forestal y que están desprovistas de ella, por lo que dichas zonas son ocupadas por cultivos, pastos con árboles o solamente pastos.

En el sector de Las Brisas, los vacíos de vegetación se concentran, de igual forma, en las zonas de alta pendiente. Este sector se caracteriza por una mayor presión en el uso de la tierra, por lo cual, la mayor parte de los vacíos están ocupados por cultivos e incluso algunas riberas están ocupadas por algún tipo de infraestructura urbana. En la parte baja de la subcuenca, el sector de Tapezco presenta junto con el PNJCB, en la parte alta, el menor porcentaje de vegetación ribereña faltante, esto obedece a que casi la totalidad del bosque en este sector se localiza a lo largo de la red hídrica. En el sector del PNJCB se encuentra la mayor área de bosque en la subcuenca, lo cual corresponde también con el menor porcentaje de vegetación ribereña faltante, con apenas un 34,48%.

Con relación a los usos de la tierra, la subcuenca del río Tapezco se caracteriza por la gran intensidad en la ocupación de ella, pero a la vez las actividades que se realizan cambian rápidamente con el paso del tiempo. Se trata del cultivo de hortalizas, que pueden ser sustituidos por otros cultivos o por el pastoreo en periodos cortos. Esta práctica resulta atractiva para los propietarios, ya que con ella estos diversifican la producción y tienen más posibilidades de obtener ingresos.

El análisis del cambio de uso de la tierra permitió, a partir de las capas de uso del 2005, 2008 y 2012, determinar las zonas que presentan mayor dinamismo en dicho cambio y cómo este afecta las coberturas boscosas en términos de fragmentación, así como la relación de estos cambios con la red de drenaje. Esto contribuyó a la identificación de zonas prioritarias para establecer enlaces de vegetación, las cuales se representan en la figura 8.

Datos: derivados de fotointerpretación de PRIAS y NASA 2005.

Figura 8 Zonas críticas por dinamismo en el uso de la tierra, subcuenca del río Tapezco (2012). 

El sector de Tapezco presenta una disminución importante de la vegetación de ribera, pasando de 11% del área total del sector en el 2005 a apenas 4% en el 2012. Algunas áreas que estaban ocupadas por bosque fueron reemplazadas por pastos. Aquí se presenta una situación similar a la del sector de Las Brisas, ya que ambos espacios constituyen la parte media-baja de la subcuenca y se encuentran muy degradados desde el 2005, punto de partida del análisis. La vegetación de ribera experimentó una disminución significativa en el periodo estudiado por lo que se presentan considerables vacíos.

Además del análisis del dinamismo en el uso de la tierra, el SDI permitió también tener una noción más acertada de la variación de los usos a través del tiempo. Se comprobó que la heterogeneidad aumentó de un valor de 1.714 en el 2005 a 1.766 en el 2012. Este aumento en la heterogeneidad del paisaje obedece a que la tierra históricamente ha sido dividida en pequeñas parcelas, que permiten mayor cambio en periodos cortos. Además de esto, la intensidad del uso de la tierra es cada vez mayor, lo que responde a las exigencias del mercado interno. La heterogeneidad de usos en el paisaje refleja la complejidad espacial en el mosaico y genera impactos no solo en el funcionamiento de los procesos ecológicos, sino también en los recursos hídricos. Esto también fue abordado en el estudio de ^{Uriarte et ál. (2011)}, donde los autores encontraron que los cambios de uso y cobertura entre periodos de medición, más que en la variación espacial, influyeron en la variación de la calidad del agua, según el indicador analizado. El tipo de cobertura presentó diferentes implicaciones en la calidad del agua, así como también varió la escala espacial de afectación según el indicador estudiado. En la subcuenca del río Tapezco es necesario contar con este tipo de análisis de mayor profundidad, para establecer los efectos producidos por cambios de uso y coberturas en la calidad de las aguas, considerando también diversas escalas.

Análisis integrado para la formulación de los enlaces

Los elementos que componen el paisaje presentan relaciones estructurales. Este vínculo estructural se analiza en términos de conectividad y está determinado por la disposición espacial de elementos. Cada sector de la subcuenca presenta un patrón de distribución distinto de dichos elementos y por esta razón se estudiaron por separado.

En el sector del PNJCB los vacíos de vegetación están ocupados mayormente por pastos arbolados y pastos. La mayoría de los vacíos se localiza en los afluentes de los cauces principales; sin embargo, se observa que el río San Rafael está descubierto en un 65% de su cauce, por lo que requeriría una mayor atención a la hora de realizarse los enlaces de vegetación. La presencia de pastos arbolados es importante, porque casi todos se encuentran en medio de parches de bosque, por lo que son un elemento clave para la conectividad entre dichos parches. En términos de cobertura vegetal, a lo largo de la red hídrica este sector presenta poca conectividad; sin embargo, en términos de conectividad entre parches la situación es más favorable, debido al tamaño y forma de estos.

En el sector de Pueblo Nuevo existen pocos parches de bosque en regeneración que puedan ser utilizados como potenciales conectores entre parches boscosos; como alternativa a esto, existe un área importante de pastos arbolados, que eventualmente podrían ser utilizados para el establecimiento de enlaces de vegetación. La vegetación arbórea a lo largo de la red hídrica en este sector presenta poca conectividad y aunque existen varios parches de bosque, estos presentan un patrón de distribución muy disperso, por lo que no se ve favorecida la conectividad.

Los sectores de Tapezco y Las Brisas se analizaron en conjunto, porque presentan características similares. La cobertura vegetal a lo largo de los cauces, aunque presenta cierta continuidad, es escasa en cuanto a área y existen varios afluentes desprovistos de cobertura boscosa, cuya área está ocupada por cultivos y pastos. Los escasos parches de bosque, presentan muy baja conectividad espacial, a pesar de que su distribución corresponde en su mayoría con la red hídrica. Las cercas vivas y las cortinas rompevientos son elementos del paisaje que podrían mejorar la conectividad en este sector. Los resultados del análisis integrado se presentan en la figura 9.

Datos: derivados de fotointerpretación de PRIAS y NASA 2005, procesados mediante la herramienta Cost Path de ArcGIS 10.0.

Figura 9 Propuesta de red de conectividad en la subcuenca del río Tapezco. 

Conclusiones

El grado de fragmentación que presenta la subcuenca del río Tapezco es producto de un intenso proceso histórico de ocupación y uso de la tierra, que inició desde la colonización del cantón de Zarcero. El uso actual de la tierra continúa reflejando gran diversidad en su empleo, presentando gran heterogeneidad en la configuración espacial del paisaje. Los usos de pastos y cultivos ocupan casi un 46% del total de área de la subcuenca, lo que da lugar a una matriz agropecuaria con una distribución muy dispersa de dichas coberturas. Este patrón resulta poco beneficioso para la continuidad del bosque, ya que aumenta la distancia entre los parches y, generalmente, las coberturas existentes entre ellos no facilitan la conectividad.

La fragmentación es desigual en la subcuenca, a partir de la cual se distinguen dos áreas diferenciadas: a) la parte alta, donde la cercanía al PNJCB influye de forma determinante en las políticas y acciones de conservación, por lo que existe una mejor conectividad y continuidad de la cobertura boscosa, es decir, una menor fragmentación; b) la parte baja de la subcuenca, debido a una mayor ocupación humana, variedad y dinamismo en el uso de la tierra, presenta una mayor fragmentación, con parches de bosque ubicados casi exclusivamente a orilla de los ríos.

Esta diferenciación se manifiesta también en la relación de densidad de parches de bosque y la densidad de drenaje por área en la subcuenca, donde una relación positiva entre estos aspectos se encontró solamente en el sector del PNJCB, contrario a lo hallado en el sector de Las Brisas, en la parte baja de la subcuenca, donde la cantidad de parches es mucho mayor y cubre menos área que en el resto de los sectores, lo cual evidenció una mayor fragmentación del bosque.

La cobertura boscosa presenta, en la mayor parte del área de estudio, una configuración que sigue el patrón de la red hídrica, con excepción del sector del PNJCB, donde, además de esto, el bosque se presenta concentrado en la parte más alta de la subcuenca, debido a que corresponde al área protegida PNJCB.

Según el análisis de cambio en el uso de la tierra, se constató que la subcuenca presenta un alto dinamismo en la transformación de los usos. Esto se comprobó mediante la determinación de las zonas críticas en los cuatro sectores que ameritan mayor atención para la mejora de la conectividad.

De manera consecuente con los resultados del análisis de fragmentación y conectividad, así como del dinamismo en el cambio de usos, los enlaces de vegetación propuestos en los sectores de la parte alta de la subcuenca se determinaron de una forma más localizada, en la mayoría de los casos, uniendo parches de bosque a través de pastos arbolados y siguiendo los recorridos más cortos entre ellos. Por su parte, en la subcuenca baja, la poca presencia de parches de bosque y las numerosas áreas de pastos y cultivos impidieron la formulación de enlaces localizados. Los enlaces propuestos en estos sectores se orientaron a mejorar la continuidad de la vegetación ribereña.

De manera complementaria a los enlaces propuestos, las barreras rompevientos y las cercas vivas fueron otro tipo de estructuras vegetales lineales de importancia en el paisaje estudiado, las cuales pueden contribuir a la conectividad estructural en un paisaje agropecuario de uso muy intensivo como el estudiado.