Preprint
Article

Spatial and Temporal Dynamics of the Vegetation Cover from the Bijagual Massif, Boyacá, Colombia, during the 1986–2021 Period

Altmetrics

Downloads

112

Views

57

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

28 December 2023

Posted:

29 December 2023

You are already at the latest version

Alerts
Abstract
Landscape changes based on spectral responses allow showing plant cover changes through diversity, composition, and ecological connectivity. The spatial and temporal vegetation dynamics of the Bijagual Massif from 1986 to 2021 were analyzed as a measure of ecological integrity, conservation, and territory. The covers identified were high open forest (Hof), dense grassland of non-wooded mainland (Dgnm), a mosaic of pastures and crops (Mpc), lagoons (Lag), and bare and degraded lands (Bdl). The Bijagual Massif has 8,574.1 ha. In 1986, Dgnm occupied 42.6% of the total area, followed by Mpc (32.8%) and Hof (24.5%); by 2000, Mpc and Hof increased (43.7 and 28.1%, respectively), while Dgnm decreased (28%); by 2021, Dgnm was restricted to the northeastern zone and continued to decrease (25.2%), Mpc occupied 52.9%, Hof 21.7% and Bdl 0.1%. Of the three fractions of the connectivity probability index, only dPCintra and dPCflux contribute to ecological connectivity. Hof and Dgnm show patches with biota habitat quality and availability. Between 1986 and 2021, Dgnm lost 1,489 ha (41%) and Hof 239.5 ha (11%). Mpc replaced various covers (1,722.2 ha; 38%) in 2021. Bijagual has a valuable biodiversity potential limited by Mpc. Territorial planning and sustainable agroecological and ecotourism proposals are required due to the context of the ecosystems.
Keywords: 
Subject: Environmental and Earth Sciences  -   Ecology

1. Introduction

High mountain ecosystems in Colombia are characterized by a degree of endemicity and richness of taxonomic entities [1,2] due to the outcome of landscape and topographic evolution, as well as climatic modeling and ecological interaction [3,4]. For example, plant physiognomy, functionality, and diversity have allowed the identification of ecological associations related to different altitudinal belts, such as high Andean forest and paramo [5,6,7,8]. In fact, these ecological associations have been considered strategic [9] because they provide ecosystem services such as food and climate, air and soil regulation, and water capacity filtration and concentration [10,11]. However, for example, in the Andean forest, human transformations have shaped the native forest landscapes into a matrix of forest patches, cultivated and grazing areas, and socio-ecosystem interactions or relationships [12,13].
Due to the socio-ecosystem relationship that occurred in this region since the 19th century [14,15,16], the vegetation cover changes over time can be identified and, hence, produce accurate information on the distribution and rates of vegetation cover chances, particularly regarding forest and paramo, which are essential for conservation strategies [17,18]. Change detection at the landscape scale can be identified by analyzing the spatiotemporal effect of vegetation reflectance [19,20,21]. Spectral variations model the vegetation reflectance using combinations of satellite image bands [17,21]. Furthermore, because Landsat satellite images can be obtained from 1970, the historical changes in plant covers can be identified over time [20,21].
In Colombia, some Andean forest and paramo areas are delimited by administrative units called complexes [22], which include the Bijagual Massif under the Tota-Bijagual-Mamapacha complex. The land use and forest deforestation history of Bijagual has been related to intensive and extensive agricultural and livestock activities [23,24,25]. However, the magnitude of human activities under the spatiotemporal effect of vegetation reflectance and, therefore, the historical, spatial, and temporal dynamics of fragmentation and ecological connectivity are unknown [26,27].
Accordingly, this work studied the spatial and temporal dynamics of the Bijagual Massif vegetation in Colombia to measure the ecological integrity of the natural area replaced by excessive human activities [25]. The supervised classification method with satellite images for a 35-year period from 1986 to 2021 was employed to obtain diversity, composition, and ecological connectivity information through landscape-scale information and coverage monitoring to suggest natural conservation strategies and territorial planning development. The results and their analyses are explained in three sections: 1. Land covers - field verification, 2. Classification accuracy, and 3. Multitemporal dimension.

2. Results

2.1. Land covers – field verification

Field verification allowed the identification of five cover types spectrally assigned to the Corine Land Cover categories (CLC; IDEAM [28]), as follows: 1. High open forest (Hof): includes the high Andean forest formation dominated by trees > 10 m tall; 2. Dense grasslands of non-wooded mainland (Dgnm): belongs to the paramo ecosystem characterized by scrubs, pajonal-chuscal, frailejonal, and sphagnal type communities; 3. Mosaic of pastures and crops (Mpc): lands whose use-occupation are pastures and crops, particularly potatoes; 4. Lagoons, lakes, and natural swamps (Lag): natural water surfaces that, for Bijagual, correspond to the La Calderona and El Pato lagoons; 5. Bare and degraded lands (Bdl): deforested area (Figure 1).
The vegetation of the Bijagual Massif belongs to Andean mountain formations that characterize the high mountain ecosystems in Colombia. According to Gil-Leguizamón et al. [29,30], the strip of high Andean forest (Hof) shows a wide species turnover at short distances, with trees and shrubs of encenillos (Weinmannia tomentosa L. f., W. rollottii Killip, W. fagaroides Kunth, W. balbisiana Kunth, and W. reticulata Ruiz & Pav.), gaque (Clusia alata Planch. & Triana, C. elliptica Kunth, and C. multiflora Kunth), susque (Brunellia propinqua Kunth, and B. cf. comocladifolia Bonpl.), garrocho (Viburnum tinoides L. f., and V. triphyllum Benth.), uva camarona (Macleania rupestris (Kunth) A.C. Sm.), tunos, charne, and siete cueros (Miconia cataractae Triana, M. cundinamarcensis Wurdack, M. jahnii Pittier, M. ligustrina (Sm.) Triana, M. theizans (Bonpl.) Cogn., Axinaea scutigera Triana, Bucquetia glutinosa (L. f.) DC., and Tibouchina grossa (L. f.) Cogn.), and manos de oso (Oreopanax bogotensis Cuatrec., and O. mutisianus (Kunth) Decne. & Planch.), among other species identified in field verification visits and corroborated with the collections of the UPTC Herbarium of Universidad Pedagógica y Tecnológica de Colombia.
In the dense grassland of non-wooded mainland or paramo coverage (Dgnm), there are associations of scrubs, pajonales-chuscales, frailejonales, and sphagnals included in the Linochilo phylicoidis-Arcytophyllion nitidi Rangel & Ariza 2000 nom. mut. nov. (all.) alliance described by Pinto Zárate [31] for the Eastern Cordillera of Colombia. Shrub communities include species such as Bejaria resinosa Mutis ex. L. f., Escallonia myrtilloides L. f., Miconia cleefii L. Uribe, Geissanthus andinus Mez, Ageratina theaefolia (Benth.) R.M. King & H. Rob., Diplostephium phylicoides (Kunth) Wedd., Pentacalia pulchella (Kunth) Cuatrec., Monnina salicifolia Ruiz & Pav., Hesperomeles obtusifolia (Pers.) Lindl., and Symplocos theiformis (L. f.) Oken., among others. The pajonal-chuscal formations are comprised of Agrostis perennans (Walter) Tuck., Anthoxanthum odoratum L., Calamagrostis effusa (Kunth) Steud., and Chusquea tessellata Munro., among others. Frailejonal formations include Espeletia argentea Bonpl., E. boyacensis Cuatrec., and E. murilloi Cuatrec., and the sphagnals formation comprises a complex of Sphagnum species [29,32].
Endemic and endangered species are found in the Bijagual Massif, including forest formation species, such as Greigia stenolepis L.B. Sm., Hieronyma rufa P. Franco, Tillandsia pallescens Betancur & García Nestor, Symplocos venulosa Cuatrec., Diplostephium oblongifolium Cuatrec., and Dunalia trianaei Dammer [30]. Furthermore, Espeletia cayetana (Cuatrec.) Cuatrec., E. murilloi Cuatrec., and Puya goudotiana Mez. [32] were found in the paramo ecosystem. In addition, species widely distributed in the high Andean forests and paramos in Boyacá were recorded, including angelitos (Monochaetum myrtoideum Naudin, M. uribei Wurdack, and M. meridense Naudin), chites (Hypericum laricifolium Juss., H. mexicanum L., H. myricariifolium Hieron., and H. lycopodioides Triana & Planch.), and cortadera (Cortaderia nitida (Kunth) Pilg.) [33,34].
The species composition in the Bijagual Massif is related to a high biodiversity potential [29]. Nevertheless, vegetation is limited by the pastures and crops cover (Mpc) mosaics extending from south to north in the study area. This cover type is dominated by potato (Solanum tuberosum L.), mora (Rubus urticifolius Poir.), and beans (Phaseolus vulgaris L.), as well as grazing cattle and sheep, considered potential threats to the ecological integrity of the verified ecosystems [35].

2.2. Classification accuracy

During the field verification, 276 control points were generated. Producer and user accuracy fluctuated between 98 and 100% of the covers. The Kappa coefficient established a thematic reliability of 95% (gray diagonal boxes; Table 1).

2.3. Multitemporal dimension

The Bijagual Massif has a total landscape area (TA) of 8,574.1 ha. In 1986, Dgnm (paramo) occupied 42.6% of the total area, followed by Mpc (32.8%) and Hof (24.5%). In 2000, the spatial composition was maintained; however, the productive systems and the forest increased (43.7% and 28.1%, respectively), while the paramo decreased (28%). In 2021, Bijagual showed paramo extensions restricted to the northeast in the municipality of Viracacha and isolated patches in the center and south (2,160 ha; 25.2% of the total area). Mpc occupies 52.9% (4,534.1 ha) and Hof 21.7% (1,861.8 ha; Table 2; Figure 1).
Between 1986 and 2000, Dgnm lost 1,248 ha (34%), while Mpc and Hof increased 936.6 ha (25%) and 309.2 ha (13%), respectively. From 2000 to 2021, Dgnm and Hof lost area (241 ha, equivalent to 10% and 548.7 ha or 23%, respectively). Furthermore, Mpc increased (785.6 ha or 17%). In 35 years, paramo and forest areas decreased (1,489 ha or 41% and 239.5 ha or 11%, respectively), while pastures and crops increased (1,722.2 ha or 38%) and occupied more than 50% of the land extension of the Bijagual Massif (Table 2; Figure 1).

2.3.1. Spatiotemporal changes 1986-2000

From 2,101.3 ha Hof registered in the year 1986, 130.6 ha changed to Dgnm (6.2%) and 92.1 ha to Mpc (4.4%) until 2000. Likewise, Dgnm was replaced in 1,182.8 ha for Mpc (32.4%) and 342.5 ha for Hof (9.4%). Similar to Hof, the coverage of Mpc changed by 191.4 ha to Hof (6.8%), 146.8 ha to Dgnm (5.2%), and 2,473.6 ha (88%) remained unchanged (Table 3). Between 1986 and 2000, of the total area of 8,574.1 ha, 6,484.3 ha maintained the same coverage (75.6%), and 2,089.9 ha changed to another cover type (24.4%).

2.3.2. Spatiotemporal changes 2000-2021

From the 2,410.5 ha Hof registered, 447.3 ha changed to Dgnm (18.6%) and 521.9 ha to Mpc (21.7%). Further, from the 2,401 ha Dgnm recorded, 257.9 ha changed to Hof (10.7%) and 740.6 ha to Mpc (30.8%). In addition, 3,748.5 ha of Mpc changed during the 2000-2021 period to another cover type; 161.6 ha changed to Hof (4.3%), 310.3 ha to Dgnm (8.3%), and solely 5.1 ha to Bdl (0.1%) (Table 3). Between 2000 and 2021, only 6,128.6 ha maintained its cover (71.5%), and 2,445.5 ha changed from one to another (28.5%).

2.3.3. Diversity

In 1986, 76 patches (TNP) corresponded to natural or transformed areas, with a density (PD) of 0.88 patches/ha. In 2000, 75 patches (PD: 0.87 patches/ha) were recorded, and in 2021, 396 patches (PD: 4.6 patches/ha) were found.

2.3.4. Composition

From the total number of patches in 1986, 25 belong to Mpc (PD: 0.3 patches/ha), but these decreased in 2000 (21 patches/PD: 0.24 patches/ha) and increased by 2021 (74/PD: 0.86). For Mpc in 1986, the largest patch index (LPI) was 18.4%, rising to 22.7% in 2000 and 48.7% in 2021. In Dgnm, the number of patches between 1986 and 2000 increased (from 25/PD: 0.3 to 31/PD: 0.36); however, the number of patches was higher in 2021 (168/PD: 1.9). Conversely, the LPI values of this coverage decreased in 35 years (1986: 21%; 2000: 8.1%; 2021: 6.7%). The abundance of Hof patches decreased between 1986 and 2000 (from 23/PD: 0.25 to 20/PD: 0.23) but increased in 2021 (150/PD: 1.7). The LPI index for Hof shows a decrease over time (1986: 8.5%; 2000: 8.2%; 2021: 6.5%; Figure 2).

2.3.5. Connectivity

Of the three fractions of the dPCk probability index, only dPCintra and dPCflux contribute to connectivity, while the contribution of the dPCconnector fraction is limited (Figure 3-A).
In 1986 and 2000, the Hof and Dgnm covers showed patches with habitat quality and availability for biota (connectivity within the patch) and act as dispersal flows (connectivity between patches); however, these ecological characteristics decreased drastically in 2021 for the two covers. Both Hof and Dgnm patches provide different connectivity values, evidenced in the dispersion of the standard deviation data concerning the average of each fraction by cover type and time (Figure 3-A).
Habitat availability contributions and flow between patches in 1986 and 2021 were low (dPCintra <10 and dPCflux <5; Figure 3-B). In the spatiotemporal dynamics, changes in the connectivity contributions of each patch were identified for the dPCintra and dPCflux fractions. In this way, in 1986, Hof showed two patches that contribute to habitat quality and availability. In 2000, one patch contributed to habitat quality and availability, two contributed to habitat quality and ecological flow, and the other two to ecological flow (Figure 3B). In 2021, only one patch contributed to habitat quality and availability and five to ecological flow (Figure 3B,C). For its part, Dgnm in 1986 showed one patch with contributions to habitat quality, habitat availability, and ecological flow (or dispersion), another with contributions to habitat quality and availability, and solely one patch with contributions to ecological flow. In 2000, only one patch registered contributions to habitat quality and availability (Figure 3B), and in 2021, one patch contributed to habitat quality, habitat availability, and ecological flow, and another solely to ecological flow (Figure 3B,C).

3. Discussion

3.1. Land cover – field verification

The Bijagual Massif corresponds to one of the six high mountain natural areas in Boyacá with broad floristic diversity [22,34]. The field verification process and the consultation of herbarium biological collections [30,32] allowed estimating a richness of up to 327 species for the Hof cover and up to 120 species for Dgnm. The results of the spatiotemporal dynamics showed that this diversity is vulnerable to the effects of the pasture and crop matrix since it has caused modifications in the vegetation [36]. According to Gil-Leguizamón et al. [29,30], the presence and dominance of some species of Chusquea Kunth in forest interiors and of Cenchrus L., Brachiaria (Trin.) Griseb., and Holcus lanatus L. in paramos indicate that the composition and physiognomy of these covers have been modified, enhancing the loss of richness and forming more heterogeneous plant communities. Even so, the research carried out by Gil-Leguizamón et al. [29,30] and Carrillo et al. [32] relate some species of ecological importance, including trees, shrubs, rosettes, and grasses that characterize the physiognomy of the Andean forest, such as C. multiflora, C. elliptica, C. alata, W. rollottii, B. comocladifolia, V. triphyllum, and E. myrtilloides, and then the species Espeletia murilloi Cuatrec., Paepalanthus columbiensis Ruhland, Puya goudotiana Mez, Hypericum lycopodioides Triana & Planch., Blechnum auratum (Fée) R. M. Tryon & Stolze, Calamagrostis effusa (Kunth) Steud., and Rhynchospora ruiziana Boeckeler typical of paramos. These records agree with the altitudinal delimitation of the vegetation proposed by Cuatrecasas [37].
Research is required for these groups to allow inferring biogeographic relationships, particularly of the Andean forest, which, through distribution records, are not only indicators that corroborate land cover but also allow interpreting the influence of the environmental gradients in the distribution patterns of the flora, its endemisms, and affinities with other biogeographic units, as has been done for the non-vascular flora of Bijagual [38].
Besides generating thematic reliability in land cover layers, field verifications allow identifying species that, due to their phenological, reproductive, and dispersal characteristics, can become a potential for ecological restoration processes in Bijagaul [35,39]. Within these, Myrcianthes rhopaloides (Kunth) McVaugh, Pentacalia pulchella (Kunth) Cuatrec., Diplostephium floribundum (Benth.) Wedd., Vallea stipularis L. f., Weinmannia fagaroides Kunth, Macleania rupestris (Kunth) A.C. Sm., Hypericum lycopodioides Triana & Planch., Bucquetia glutinosa (L. f.) DC and some species of Viburnum L. and Cestrum L. are highlighted.

3.2. Classification accuracy

From the 276 control points obtained during the field verification and the overlap of these points with the land cover raster layer (spectral response obtained from the supervised classification for the 2021 image), a thematic reliability of 95% (Kappa 0.95) was obtained. According to cartographic quality standards, the data are associated with a "very good spatial concordance" [40]. Of these 276 points, nine did not coincide with their respective classification cover due to the similarity of the digital value of a pixel in the image associated with another cover. This was evidenced for Hof with six points corresponding to forest cover in the field verification and that in the supervised classification algorithm were assigned to Dgnm (4 points) and Mpc (2 points). A similar case was obtained for Lag with only two points and Dgnm with one (Table 1). The remaining 267 control points maintain a concordance between verification and classification, which allowed producer and user accuracy to be higher than 98%, and the errors of omission and commission were not higher than 2% (Table 1). For this reason, the pre- and post-processing images were considered reliable.

3.3. Multitemporal dimension

Changes in plant cover due to changes in land use affect the properties and functioning of ecosystems [41,42]. The vegetation and natural covers in the Colombian Andean relief have been modified as a result of human colonization processes. For 35 years in Bijagual, extensive and small areas of forest and paramo were substituted for agricultural and livestock activities, which currently affect habitat loss for biodiversity [23,43]. The landscape structure [44], species dynamics [5], and the ecological network conformation of the natural space of the Tota-Bijagual-Mamapacha complex [22], currently threatened by human activities, have caused [45] a decrease in forest area and an increase in the natural distance of patches, and, therefore, the loss of this ecological continuity corridor [45,46,47] (Figure 1; Table 2).
Forest fragmentation has been more evident since 2021 due to the increase of Dgnm and Hof in the number of patches, but also due to area reduction and patch isolation. For 35 years, 1,923.4 and 614 ha of Dgnm and Hof, respectively, were converted to Mpc (Table 3), supported by the increase in the number of patches, the area, and the ones that become connected (Figure 1 and Figure 2).
Areas around the massif, which have been proposed as transition zones, were found to be affected by crops and livestock activities long before the 1980s. According to Rodríguez-Eraso et al. [48], the loss of natural area in the 80s and 90s was centered in the eastern Cordillera in Boyacá and Cundinamarca, caused by deforestation patterns of the high Andean forest and paramo due to food demand. In addition, Sánchez-Cuervo et al. [24], in the departments of Santander and Boyacá, pointed out that potato cultivation increased from 380 to 482 km2 between 2006 and 2008, reducing forest covers.
Potato cultivation has a high per capita consumption in Colombia, and, therefore, expansion has occurred under native forest areas, increasing the continuous degradation of Andean forests and the paramos. Furthermore, another factor influencing forest cover is climate change, which has been enhancing during the last few years and has affected species distribution [49,50]. Since the Andean high mountain ecosystems are seasonal, the increase in temperature would prolong the intensity of the dry periods, and the paramo would reduce the capacity to capture and regulate water [36]. In the Colombian eastern mountain range, only 36.4% of the total paramo area (322,925 ha) and 25.3% of the Andean montane forest (427,675 ha that include the sub-Andean and Andean forest belts) are protected. However, of these areas, only 41 and 45% maintain their original extension [46]. The remaining areas have experienced land use changes with subsequent biodiversity loss [36].
The paramo cover (Dgnm) decreased in area between 1986 and 2021 (1,489 ha), while the forest expanded its area between 1986 and 2000 (309.2 ha) and decreased between 2000 and 2021 (548.7 ha). The forest area that decreased in the second period is larger than the area that expanded in the first (with a difference of 239.5 ha) (Table 2). Instead, the pasture and crop cover increased by 1,722.2 ha in 35 years. Therefore, it is necessary to conserve the Dgnm and Hof areas that still exist in the municipalities of Viracachá, Cienaga, and Tibana (4,034.9 ha [15,18,47]) and include structure, functionality, and ecological diversity attributes in the Bijagual extension [6,35,39], particularly in the municipality of Ramiriquí, which is the most affected by the expansion of the agricultural and livestock frontier (Figure 1).
Spectral image response outcomes show that natural succession has occurred around the paramo and forest areas (Figure 1; Table 3) [23,51,52]. Nevertheless, this process seems to be slow. Lequerica et al. [53] and Prado-Castillo et al. [35] suggest that succession is slow in areas previously used for cultivation due to degradation, and the compaction of soil produced by agriculture can limit the development of propagules and modifies the germinable or viable seed bank [54]. Hence, natural succession needs to be assisted by ecological restoration processes to improve the increase of paramo vegetation.
The increase in the total number of patches (TNP 76 to 396 between 1986 and 2021), density (PD), and percentage of patch occupation (LPI) are evidence of the intensive and extensive land use and the changes from one coverage to another. While in 1986, one patch was registered for each hectare, in 2021, four or five patches (PD) were recorded for each hectare. In 2021, Mpc coverage presented the lowest number of patches compared to Dgnm and Hof (NP 74, 168, and 150, respectively; Figure 2). However, Mpc presents patches with a larger occupation area compared to Dgnm and Hof (LPI 48.7, 6.7, and 6.5%, respectively; Figure 2). In 35 years, the LPI metric describes greater spatiotemporal occupancy of the Mpc patches and lower occupancy of the Dgnm and Hof patches (Figure 2).
The cover and topological position of the Bijagual landscape per patch allow the identification of different contributions to ecological connectivity. Each patch can act as a space in which there is connectivity, i.e., a larger area with a suitable habitat condition or quality that may have higher connectivity [55,56]. The functionality of some patches is evident, even more over those that show higher habitat availability (Figure 3 B,C) since these can act as refuge, feeding and reproduction areas or are habitat units that produce or receive dispersal flows of individuals to other units or patches [56]. In fact, these areas can be analyzed in the metapopulation dynamic to establish possible unknown connections with the biotic complex of Tota to the north and Mamapacha to the south, which, together, comprise a biological corridor.
The spatial configuration of the Hof and Dgnm patches has contributed to the ecological connectivity in Bijagual from 1986 to 2021, evidenced by the contributions of the dPCintra and dPCflux fractions; however, during this period, also natural area loss has occurred, affecting the ecological diversity [29,57,58]. Ecological connectivity in Bijagual decreased noticeably from 2000 to 2021 (Figure 3A). The massif may experience a fragmentation process in the near future [15,45,46], and the current landscape of this natural area shows a loss in habitat quality and availability [45,55,56].
As a result of the distancing and isolation of the Dgnm and Hof patches, the continuity of the ecological corridor has been interrupted (Figure 1). This, in turn, poses difficulties for the movement of species and genetic exchange [56,59] and is evident in the reduced contribution of the dPCconnector fraction (Figure 3A) that determines limited dispersion fluxes. However, the Bijagual Massif still has paramo and Andean forest patches that can act as critical habitat areas for maintaining ecological connectivity (Figure 3C), where habitat loss or deterioration would have an even more negative impact on landscape connectivity. These areas must be prioritized and considered reference ecosystems for ecological restoration processes.
This research states that the Bijagual Massif requires the inclusion of sustainable management practices and territorial planning with agroecological, forestry, environmental education, ecotourism, biogeographical, and climate research proposals, according to the natural environment considering the strategic and vulnerable ecosystem context where the State interacts through its territorial entities (government and environment), the peasant communities (through citizen participation actions of the associated communities), agricultural companies, and the academy as knowledge generators.

4. Materials and Methods

4.1. Study area

The Bijagual Massif is located in the Andes eastern mountain range over the Department of Boyacá (Colombia) between 5°26′0.0″-5°16′0.0″ N and 73°14′0.0″-73°22′0.0″ W [29]. The vegetation is composed of paramo (between 2,990 and 3,460 m a.s.l.) and high Andean forest (between 2,682 and 3,268 m a.s.l.) in an extension of 8,574.1 ha (Figure 4). Its borders delimit to the north with Lake Tota and the paramos of Toquilla, Sarna, Suse, and Los Curíes and to the south with the Mamapacha Massif. Together, they shape the Tota-Bijagual-Mampacha paramo complex [19]. Precipitation values range between 1,000 and 2,500 mm/year, and temperatures range between 13 and 17°C [32].

4.2. Data source

Satellite data interpreted with landscape metrics (i.e., diversity, composition, and configuration) were used to analyze the spatial and temporal dynamics of the Bijagual Massif vegetation. The vegetation cover thematic maps were edited according to the method proposed by The Nature Conservancy [60] and Posada [61] at a scale of 1:100,000. Three images from the Landsat 5 Thematic Mapper (TM; year 1986), Landsat 7 Enhanced Thematic Mapper (ETM+; year 2000), and Landsat 8 Operational Land Imager (OLI; year 2021) sensors with a spatial resolution of 30 m, Path/ Row 7/56, and cloudiness < 10% (source: Earth Explorer of the USGS - U.S. Geological Survey) were employed. Band combinations in each image were made according to the specifications of each sensor to differentiate covers. For the 1986 and 2000 images, the combination was Red: 4, Green: 5, and Blue: 3, and for 2021, Red: 5, Green: 6, and Blue: 2. Pre- and post-processing was performed with the programs ERDAS Image ver. 2016 and ArcGIS ver. 10.7 (academic licenses) described below.

4.3. Geometric correction

One hundred and fifty control points (in images 1986 and 2021) were generated from a third-order polynomial model and a georeferenced image (image from 2000; [17,40]); in the 1986 image, 21 automatic points and 129 manual points were projected (standard error-RMSE 0.59; Std. Dev: 0.218); in 2021, 10 automatic and 128 manual points (RMSE 0.55; Erros Std. Dev: 0.208) were projected. These were considered reliable because their values were < 1.0 [40]. Subsequently, the study area was cut from the Bijagual polygon.

4.4. General contrasts (spatial enhancement)

The images were smoothed in their spatial contrasts similar to the digital Level of each pixel with respect to neighboring pixels using low-pass filters (Kernel 3X3) to reduce spatial variability and visualize less sharp, more blurred, and homogeneous profiles.

4.5. Land covers (thematic categorization) – field verification

A statistical digital classification was performed for each image to group the pixels by cover, i.e., going from a continuous measurement (digital levels of each pixel) scale to a categorical one (vegetation coverage in the Corine Land Cover system – IDEAM [28]) using the supervised classification method and the selection and delimitation of training areas (grouping of pixels of the same category).
The supervised classification method was supported with field verification activities, which involved obtaining geographic coordinates (control points in the WGS84 system) and floristic and vegetation characterizations in each natural cover type (published in Gil-Leguizamón et al. [29]). In addition, for endemic and threatened species, Herbarium collections and online catalogs were consulted [30,32].

4.6. Classification accuracy

The reliability of the training areas (spectral separability of the categories) was statistically evaluated through an error matrix that contrasted the real classes (the truth of the terrain) and the cartographic units (cover). Additionally, accuracy measurements were used for each cover to determine the reliability of the supervised classification, including:

4.6.1. Producer accuracy

The omission error measure indicates the probability that a sample point is correctly classified. It is calculated by dividing the number of pixels correctly classified into a category by the total number of pixels in that category according to Equation (1).
P P % = ( X i i / X + i ) 100
Where PP% is producer accuracy in %, Xii is the value of the column diagonal, and X+i is the total marginals of column i.

4.6.2. User accuracy

The probability that a sampling point corresponds to the category assigned in the field (commission error). It is calculated by dividing the total number of correct pixels in a category by the total number of pixels that, in fact, belong to this category according to Equation (2).
P U % = ( X i i / X i + ) 100
Where PU% is user accuracy in %, Xii is the diagonal value of row I, and Xi+ is the total marginals of this row.

4.6.3. Kappa coefficient

It measures the difference between what is observed between classes and what is evidenced in the classification according to Equation 3. The veracity corresponds to values close to or equal to 1.
K = ( O E ) / ( 1 E )
Where O is the correct proportion of the classes identified in the image, and E is the correct observed percent change estimate.

4.7. Multitemporal dimension (landscape metrics)

Statistics to identify the spatial and temporal dynamics of the coverages for 1986, 2000, and 2021 were obtained according to the spatial pattern theory [62]. This incorporates the landscape as a mosaic of patches or discrete covers, which explain the distribution of geographic objects, their patterns, and processes in time [62,63]. The statistics used quantified the landscape with three parameters: diversity and composition with Fragstats [62] and connectivity with Conefor [55,56,64], as follows:

4.7.1. Diversity

It was analyzed through the total landscape area (TA, ha), the total number of patches (TNP), patch richness (PR), and patch richness density (PRD).

4.7.2. Composition

Includes the area of each cover or class (CA, ha), the number of patches (NP), patch density (PD), landscape percentage (PLAND, %), and the largest patch index (LPI, %) to analyze the composition.

4.7.3. Connectivity

The topo-ecological index of connectivity probability (dPCk) was used and interpreted according to three fractions that quantified how the patches (as habitat or corridor) contribute to the global connectivity of the Bijagual landscape, according to Equation 4, to obtain connectivity information.
d P C k = d P C i n t r a k + d P C f l u x k + d P C c o n n e c t o r k
Where k is a landscape element (patch), intra is the contribution of patch k in terms of the habitat area (or quality) available within it (intrapatch connectivity), and flux corresponds to the dispersion flux (weighted by the area) received or originated through the connections of patch k with the rest of the patches (it reflects how well-connected k is with the rest of the existing habitat in the landscape). Connector evaluates the contribution of k as a connecting element or bridge patch between the rest of the patches (estimates to what extent k facilitates dispersive flows and which do not have their origin or destination in k but are enhanced and pass through k) [56,65]. Furthermore, an intermediate dispersal distance (0.5 probability) and an approximate average distance between patches of 50 m were assumed. The results are shown as bar and scatter plots obtained with the ggplot2 package [66] included in the open-source RStudio Desktop version software of 2003 for Windows 10/11.

Author Contributions

Conceptualization, P.A.G.-L., J.F.P-O., D.S.-M., AES-T., J.D.M.-G. and M.E.M.-P.; methodology, P.A.G.-L., J.D.M.-G. and M.E.M.-P.; software, P.A.G.-L.; validation, P.A.G.-L., J.F.P.-O., D.S.-M., A.E.S.-T., J.D.M.-G. and M.E.M.-P.; formal analysis, P.A.G.-L., J.F.P.-O., D.S.-M., A.E.S.-T., J.D.M.-G. and M.E.M.-P.; investigation, P.A.G.-L., J.D.M.-G. and M.E.M.-P.; resources, P.A.G.-L. and M.E.M.-P.; data curation, P.A.G.-L.; writing—original draft preparation, P.A.G.-L. and M.E.M.-P.; writing—review and editing, P.A.G.-L., J.F.P.-O., D.S.-M., A.E.S.-T., J.D.M.-G. and M.E.M.-P.; visualization, P.A.G.-L., J.F.P.-O., D.S.-M., A.E.S.-T., J.D.M.-G. and M.E.M.-P.; supervision, P.A.G.-L. and M.E.M.-P.; project administration, P.A.G.-L. and M.E.M.-P.; funding acquisition, P.A.G.-L. and M.E.M.-P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the collaboration agreement DHS No. 5211071, biodiversity and conservation objects of the strategic ecosystems of Sucuncuca and Bijagual (Province of Lengupá) Boyacá [biodiversidad y objetos de conservación de los ecosistemas estratégicos de Sucuncuca y Bijagual (provincia de Lengupá) Boyacá], established between Universidad Pedagógica y Tecnológica de Colombia and the company Ecopetrol S. A.

Data Availability Statement

All data supporting the reported results are included in the paper. Biological records registered in each coverage are available in SibColombia: https://ipt.biodiversidad.co/cr-sib/pdf.do?r=0724_bijagual_20180824&n=1656C6909E7.

Acknowledgments

The authors thank Universidad Pedagógica y Tecnológica de Colombia (UPTC), the UPTC Herbarium, the Doctoral Program in Biological and Environmental Sciences (DCBA - UPTC), and the Doctoral excellence scholarships, second Bicentennial cohort – MINCIENCIAS [Becas de excelencia doctoral del Bicentenario segunda cohorte - MINCIENCIAS]. We also acknowledge Ecopetrol-UPTC [DHS 5211071], SisBio, CORPOCHIVOR, and the communities of the municipalities of Viracachá, Ciénega, Ramiriquí, and Tibaná. Further, many thanks to Semillero Phytos.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alzate, F.; Jiménez-Montoya, J.; Sarrazola, H. Sonsón, un nuevo complejo de paramos en Colombia. In Flora del Oriente Antioqueño. Biodiversidad, ecología y estrategias de conservación, Quijano-Abril, M.A., Ed.; Fondo Editorial Universidad Católica de Oriente: Bogotá D.C., Colombia, 2016; pp. 19–24. [Google Scholar]
  2. Jiménez-Rivillas, C.; García, J.J.; Quijano-Abril, M.A.; Daza, J.M.; Morrone, J.J. A new biogeographical regionalisation of the Paramo biogeographic province. Australian Systematic Botany 2018, 31, 296–310. [Google Scholar] [CrossRef]
  3. Cuesta, F., Muriel, P., Llambí, L.D., Halloy, S., Aguirre, N., Beck, S., Carilla, J., Meneses, R.I., Cuello, S., Grau, A., Gámez, L.E, Irazábal, J., Jácome, J., Jaramillo, J., Ramírez, L., Samaniego, N., Suárez-Duque, D., Thompson, N., Tupayachi, A., Viñas, P., … Gosling, W.D. Latitudinal and altitudinal patterns of plant community diversity on mountain summits across the tropical Andes. Ecography 2017, 40, 001–014.
  4. Fadrique, B., Báez, S., Duque, Á., Malizia, A., Blundo, C., Carilla, J., Osinaga-Acosta, O., Malizia, L., Silman, M., Farfán-Ríos, W., Malhi, Y., Young, K.R., Cuesta, C. F., Homeier, J., Peralvo, M., Pinto, E., Jadan, O., Aguirre, N., Aguirre, Z., Feeley, K.J. Widespread but heterogeneous responses of Andean forests to climate change. Nature 2018, 564, 207–212. [CrossRef]
  5. Llambí, L.D.; Cuesta, F. La diversidad de los paramos andinos en el espacio y en el tiempo. In Avances en investigación para la conservación de los paramos andinos, Cuesta, F.J., Sevink, L., Llambí, B., Posner J., Eds.; CONDESAN: Quito, Ecuador, 2014; pp. 8–39. [Google Scholar]
  6. Suárez Duque, D. Formación de un corredor de hábitat de un bosque montano alto en un mosaico de paramo en el norte del Ecuador. Ecología Aplicada 2008, 7, 9–15. [Google Scholar] [CrossRef]
  7. Rangel-Ch. O. La riqueza de las plantas con flores de Colombia. Caldasia 2015, 37, 279–307. [CrossRef]
  8. Abud, M.; Torres, A.M. Caracterización florística de un bosque alto andino en el Parque Nacional Natural Puracé, Cauca, Colombia. Boletín Científico del Museo de Historia Natural 2016, 20, 27–39. [Google Scholar] [CrossRef]
  9. Llambí, L.D., Soto, A., Celleri, R., De Bievre, B., Ochoa, B., Borja, P. Paramos Andinos: Ecología, Hidrología y Suelos de Paramos. Monsalve Moreno: Bogotá, Colombia 2012; pp. 283.
  10. Rincón-Ruíz, A.; Echeverry-Duque, M.; Piñeros, A.M.; Tapia, C.H.; David, A.; Arias-Arévalo, P.; Zuluaga, P.A. Valoración integral de la biodiversidad y los servicios ecosistémicos: aspectos conceptuales y metodológicos. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt (IAvH): Bogotá, D.C., Colombia 2014; pp. 151.
  11. Nieto, M.; Cardona, L.; Agudelo, C. Análisis de servicios ecosistémicos. Provisión y regulación hídrica. In Hojas de ruta. Guías para el estudio socioecológico de la alta montaña en Colombia, Ungar, P., Ed.; Instituto de Investigación de Recursos Biológicos Alexander von Humboldt, Bogotá, Colombia 2015; pp. 1–95.
  12. Salas-Zapata, W.; Ríos-Osorio, L.V. Álvarez, J. Bases conceptuales para una clasificación de los sistemas socioecológicos de la investigación en sostenibilidad. Revista Lasallista de Investigación 2012, 8, 136–142. [Google Scholar]
  13. Uribe-Castro, H. De ecosistema a socioecosistema diseñado como territorio del capital agroindustrial y del Estado-nación moderno en el valle geográfico del río Cauca, Colombia. Revista Colombiana de Sociología 2014, 37, 121–157. [Google Scholar]
  14. Porter-Bolland, L., et al. Community managed forest and forest protected areas: An assessment of their conservation effectiveness across the tropics. Forest Ecology and Management. 268, pp. 6–17.
  15. Arellano-P.H.; Rangel-Ch, J. O. Fragmentación y estado de conservación en algunos paramos de Colombia. In Colombia, diversidad biótica X, cambio global (Natural) y climático (antrópico) en el paramo colombiano, Rangel, J.O., Ed.; Universidad Nacional de Colombia: Bogotá, Colombia 2010; pp. 479–542.
  16. Moncada-Rasmussen, D. Análisis espacio-temporal del cambio en los bosques de roble (Quercus humboldtii Bonpl.) y su relación con la alfarería en Aguabuena (Ráquirá–Boyacá). Colombia Forestal 2010, 13, 275–298. [Google Scholar] [CrossRef]
  17. Chuvieco, E. Teledetección ambiental, la observación de la Tierra desde el espacio. Editorial Ariel S.A.: Barcelona, Spain 2002; pp. 586.
  18. Chapin, F.S.; Matson, P.A.; Vitousek, P.M. Principles of Terrestrial Ecosystem Ecology. 2nd ed.; Springer Science+Business Media: New York, USA 2011; pp. 392. Principles of Terrestrial Ecosystem Ecology.
  19. McDonald, A.J.; Gemmell, F.M.; Lewis, P.E. Investigation of the utility of spectral vegetation indices for determining Information on Coniferous Forests. Remote Sensing of Environment 1998, 66, 250–272. [Google Scholar] [CrossRef]
  20. Hantson, S.; Chuvieco, E.; Pons, X.; Domingo, C.; Cea, C.; More, G.; Cristobal, J.; Peces, J.; Tejeiro, J. Cadena de pre-procesamiento estándar para las imágenes Landsat del Plan Nacional de Teledetección. Revista de Teledetección 2011, 36, 51–61. [Google Scholar]
  21. Coops, N.C.; Tooke, T.R. Introduction to remote sensing. In Learning Landscape Ecology. A Practical Guide to Concepts and Techniques, Gergel, S.E., Turner M.G., Eds.; Springer Verlag: New York, NY, USA, 2017; pp. 3–20. [Google Scholar]
  22. Morales, M.; Otero, J.; van der Hammen, T.; Torres, A.; Cadena, C.; … Cárdenas, L. Atlas de paramos de Colombia. Instituto de Investigación de Recursos Biológicos Alexander von Humboldt: Bogotá, Colombia 2007; pp. 208.
  23. Etter, A.; McAlpine, C.; Wilson, K.; Phinn, S.; Possingham, H. Regional patterns of agricultural land use and deforestation in Colombia. Agriculture, Ecosystem & Environment 2006, 114, 369–386. [Google Scholar]
  24. Sánchez-Cuervo, A.M.; Aide, M.; Clark, M.; Etter, A. Land cover change in Colombia: surprising forest recovery trends between 2001 and 2010. PLoS ONE 2012, 7, 1–14. [Google Scholar] [CrossRef] [PubMed]
  25. López, A.; Caicedo, E.; Mojica, J.A.; López Bermúdez, D.A. Paramos. Altas temperaturas y depredación aceleran su muerte. Available online: https://www.eltiempo.com/vida/medio-ambiente/paramos-asi-los-afecta-el-cambio-climatico-mineria-agricultura-y-conflicto-791673 (accessed on 13 September 2023).
  26. Pérez, M.E.; García, M. Aplicaciones de la teledetección en degradación de suelos. Boletín de la Asociación de Geógrafos Españoles 2013, 61, 285–308. [Google Scholar]
  27. Saura, S. Métodos y herramientas para el análisis de la conectividad del paisaje y su integración en los planes de conservación. In Avances en el análisis espacial de datos ecológicos: aspectos metodológicos y aplicados. De la Cruz, M., Maestre, F.m Eds.; ECESPA-Asociación Española de Ecología Terrestre: Mósteles, España: Mósteles, 2013; pp. 2–46.
  28. IDEAM. Leyenda nacional de coberturas de la Tierra Corine Land Cover adaptada para Colombia, escala 1:100.000. Instituto de Hidrología, Meteorología y Estudios Ambientales. Editorial Scripto Ltda.: Bogotá, Colombia 2010; pp. 72.
  29. Gil-Leguizamón, P.A.; Morales-Puentes, M.E.; Jácome, J. Estructura del bosque altoandino y paramo en el Macizo de Bijagual, Boyacá, Colombia. Revista de Biología Tropical 2020, 68, 765–776. [Google Scholar] [CrossRef]
  30. Gil-Leguizamón, P.A.; Morales-Puentes, M.E.; Carrillo-Fajardo, M.Y. Species richness of the vascular plants of the Bijagual high Andean forest, Colombia. Acta Botánica Mexicana,. 2021, 128, e1869. [Google Scholar] [CrossRef]
  31. Pinto Zárate, J. La vegetación de la región de vida paramuna colombiana: con énfasis en la cordillera Oriental y los Andes Centrales. Doctor Degree, Universidad Nacional de Colombia, Bogotá, Colombia 2023.
  32. Carrillo, M.Y.; Morales-Puentes, M.E.; Gil-Leguizamón, P.A. Catálogo de la flora (angiospermas) del paramo de Bijagual, Boyacá, Colombia. Boletín Científico del Centro de Museos 2017, 21, 15–29. [Google Scholar] [CrossRef]
  33. Gil-Leguizamón, P.A.; Morales-Puentes, M.E.; Díaz-Pérez, C.N. Evaluación poblacional y grado de amenaza de Chalybea macrocarpa (Melastomataceae) especie endémica de Colombia. Acta biológica colombiana 2014, 19, 261–270. [Google Scholar]
  34. Olaya-Angarita, J.A.; Díaz-Pérez, C.N.; Morales-Puentes, M.E. Composición y estructura de la transición bosque-paramo en el corredor Guantiva-La Rusia (Colombia). Revista de Biología Tropical 2019, 67, 755–768. [Google Scholar] [CrossRef]
  35. Prado-Castillo, L.F.; Cárdenas-Tamayo, T.A.; Carrillo-Fajardo. M.Y. Aproximación conceptual para la restauración de ecosistemas andinos. In Restauración de ecosistemas de montaña: cultura y ecología desde el paramo y el piedemonte llanero, Prado-Castillo, L.F., Gil-Leguizamón, P.A., Sabogal-González, A., Morales-Puentes, M.E., Eds.; Editorial: Tunja, Colombia 2018; pp. 17–30.
  36. Caicedo, E. Donde había frailejones, bromelias y mortiños, hoy hay cebolla, papa y vacas. El Tiempo. Available online: https://www.eltiempo.com/vida/medio-ambiente/paramos-asi-es-como-agricultura-y-ganaderia-han-transformado-los-ecosistemas-792130 (accessed on 8 September 2023).
  37. [Cuatrecasas, J. Aspectos de la vegetación natural de Colombia. Revista de la Academia Colombiana de Ciencias 1958, 10, 221–268. [Google Scholar]
  38. Gil-Novoa, J.E.; Morales-Puentes, M.E.; Mercado-Gómez, J.D. Biogeographic origins and floristic affinities of the bryophytes of the Tota-Bijagual-Mamapacha complex in the Páramos of Boyacá, Colombia. Phytotaxa 2020, 477, 171–193. [Google Scholar] [CrossRef]
  39. Gil-Leguizamón, P.A.; Bravo-Pedraza, W.J.; Hernández-Velandia, D.R.; Moreno-Mancilla, O.F.; Morales-Alba, A.F.; Reyes Camargo, J.E.; Ovalle-Pacheco, A.L.; Muños-Avila, J.A.; Vergara, P.; Ortiz, L.; Useche, D.S. Restauración ecológica del bosque altoandino en la vereda Montaya (Ventaquemada, Boyacá): Avances. In Restauración de ecosistemas de montaña: cultura y ecología desde el paramo y el piedemonte llanero, Prado-Castillo, L.F., Gil-Leguizamón, P.A., Sabogal-González, A., Morales-Puentes, M.E., Eds.; Editorial: Tunja, Colombia, 2018; pp. 139–178. [Google Scholar]
  40. Posada, E. Manual de prácticas de percepción remota parte 2. Centro de Investigación y Desarrollo en Información Geográfica (CIAF). Instituto Geográfico Agustín Codazzi (IGAC): Bogotá, Colombia 2012; pp. 154.
  41. Lambin, E.; Turner, B.; Geist, H.; Agbola, S.; Angelsen, A.; … Xu, J. The causes of Land-use and Land-cover change: moving beyond the myths. Global Environmental Change 2001, 11, 261–269. [Google Scholar] [CrossRef]
  42. Hansen, M.; Stehman, S.; Potatov, P.; Loveland, T.; Townshend, J.; Defries, R.; Pittman, K.; Arunarwati, B.; Stolle, F.; Steininger, M.; Carroll, M.; Dimiceli, Ch. Humid tropical forest clearing from 2000 to 2005 quantified by using multitemporal and multiresolution remotely sensed data. Proceedings of the National Academy of Sciences of the United States of America. 2008, 105, 9439–9444. [Google Scholar] [CrossRef]
  43. Etter, A.; McAlpine, C.; Possingham, H. Historical patterns and drivers of landscape change in Colombia since 1500: A regionalized spatial approach. Annals of the Association of American Geographers 2008, 98, 2–23. [Google Scholar] [CrossRef]
  44. Cortés-Duque, J.; Sarmiento, C. Visión socioecosistémica de los paramos y la alta montaña colombiana: memorias del proceso de definición de criterios para la delimitación de paramos; Instituto de Investigación de Recursos Biológicos Alexander von Humboldt: Bogotá, D.C. Colombia 2013; pp. 235.
  45. Bogaert, J.; Ceulemans, R.; Salvador-Van Eysenrode, D. Decision Tree Algorithm for Detection of Spatial Processes in Landscape Transformation. Environmental Management 2004, 33, 62–73. [Google Scholar] [CrossRef]
  46. Armenteras, D.; Gast, F.; Villarreal, H. Andean forest fragmentation and the representativeness of protected natural areas in the eastern Andes, Colombia. Biological Conservation 2003, 113, 245–256. [Google Scholar] [CrossRef]
  47. Fazey, L.; Fischer, J.; Lindenmayer, D. What do conservation biologists publish? Biological Conservation 2005, 124, 63–73. [Google Scholar] [CrossRef]
  48. Rodríguez-Eraso, N.; Armenteras-Pascual, D.; Alumbreros, J. Land use and land cover change in the Colombian Andes: dynamics and future scenarios. Journal of Land Use Science 2013, 8, 154–174. [Google Scholar] [CrossRef]
  49. Buytaert, W.; Sevink, J.; Cuesta, F.; Cuesta, F.; LLambi, L.; Posner, J. Cambio climático: la nueva amenaza para los paramos. Avances en investigación para la conservación de los paramos andinos. CONDESAN: Lima, Perú 2014; pp. 525.
  50. Morocho, C.C., Chuncho, G. Paramos del Ecuador, importancia y afectaciones: Una revisión. Bosques Latitud Cero 2019, 9, 71–83.
  51. Arroyo-Rodríguez, V.; Melo, F.P.L.; Martínez-Ramos, M.; Bongers, F.; Chazdon, R.L.; Meave, J.A.; Norden, N.; Santos, B.A.; Leal, I.R.; Tabarelli, M. Multiple successional pathways in human-modified tropical landscapes: new insights from forest succession, forest fragmentation and landscape ecology research. Biological Reviews Cambridge Philosophical Society 2015, 92, 326–340. [Google Scholar] [CrossRef] [PubMed]
  52. McDonald, T.; Jonson, J.; Dixon, K.W. National standards for the practice of ecological restoration in Australia. Restoration Ecology 2016, 24, S6–S32. [Google Scholar] [CrossRef]
  53. Lequerica, M.; Bernal, M.; Stevenson, P. Evidencia de direccionalidad del proceso de sucesión temprana del bosque altoandino. Colombia Forestal 2017, 20, 63–84. [Google Scholar] [CrossRef]
  54. Van Andel, J.; Aronson, J. Restoration ecology: the new frontier; Wiley Blackwell: Oxford, UK, 2012; p. 400. [Google Scholar]
  55. Pascual-Hortal, L.; Saura, S. Comparison and development of new graph-based landscape connectivity indices: towards the priorization of habitat patches and corridors for conservation. Landscape Ecology 2006, 21, 959–967. [Google Scholar] [CrossRef]
  56. Saura, S.; Rubio, L. A common currency for the different ways in which patches and links can contribute to habitat availability and connectivity in the landscape. Ecography 2010, 33, 523–537. [Google Scholar] [CrossRef]
  57. Forero-Gómez, Y.K., Gil-Leguizamón, P.A., Morales-Puentes, M.E. Structural connectivity between the Páramos of Guacheneque and Los Cristales, Rabanal-Río Bogotá complex, Colombia. Revista de Teledetección 2020, 57, 65–77. [CrossRef]
  58. Anzoategui, L.A.; Gil-Leguizamón, P.A.; Sanabria-Marin, R. Agricultural frontier and multi-temporality of vegetation cover in moorland of the Cortadera Regional Natural Park (Boyacá, Colombia). BOSQUE 2023, 44, 159–170. [Google Scholar] [CrossRef]
  59. Valdés, A. Modelos de paisaje y análisis de fragmentación: de la biogeografía de islas a la aproximación de paisaje continuo. Ecosistemas 2011, 20, 11–20. [Google Scholar]
  60. The Nature Conservancy. Estudio multitemporal del nor-oriente del Ecuador. Memoria técnica estudio multitemporal de la cobertura vegetal y uso del suelo en los años 1990-2008 y proyección 2030. Geoplades: Ecuador 2009; pp. 31.
  61. Posada, E. Manual de prácticas de percepción remota parte 1. Centro de Investigación y Desarrollo en Información Geográfica (CIAF). Instituto Geográfico Agustín Codazzi (IGAC): Bogotá, Colombia 2012; pp. 154.
  62. McGarigal, K.; Cushman, S.; Neel, M.; Ene, E. FRAGSTATS: spatial pattern analysis program for categorical maps. Computer software program produced by the authors at the University of Massachusetts, Amherst. Available online: https://www.umass.edu/landeco/research/fragstats/fragstats.html (accessed on 31 January 2012).
  63. Magle, S.; Theobald, D.; Crooks, K.A. Comparison of metrics predicting landscape connectivity for a highly interactive species along an urban gradient in Colorado, USA. Landscape Ecology 2009, 24, 267–280. [Google Scholar] [CrossRef]
  64. Saura, S.; Pascual-Hortal, L. A new habitat availability index to integrate connectivity in landscape conservation planning: Comparison with existing indices and application to a case study. Landscape and Urban Planning 2007, 83, 91–103. [Google Scholar] [CrossRef]
  65. Pascual-Hortal, L.; Saura, S. Integrating landscape connectivity in broad-scale forest planning through a new graph-based habitat availability methodology: application to capercaillie (Tetrao urogallus) in Catalonia (NE Spain). European Journal of Forest Research volume 2008, 127, 23–31. [Google Scholar] [CrossRef]
  66. Wickham, H.; Chang, W.; Henry, L.; Lin Pedersen, T.; Takahashi, K.; Wilke, C.; Woo, K.; Yutani, H.; Dunnington, D. ggplot2. Available online: https://ggplot2.tidyverse.org/ (accessed on 5 October 2023).
Figure 1. Composition and spatial configuration of the vegetation cover in the Bijagual Massif; A: Cover in 1986, B: Cover in 2000, and C: Cover in 2021; Dgnm: Dense grasslands of non-wooded mainland (paramo); Hof: High open forest (high Andean forest); Mpc: Mosaic of pastures and crops; Lag: Lagoons; Bdl: Bare and degraded lands (deforested areas).
Figure 1. Composition and spatial configuration of the vegetation cover in the Bijagual Massif; A: Cover in 1986, B: Cover in 2000, and C: Cover in 2021; Dgnm: Dense grasslands of non-wooded mainland (paramo); Hof: High open forest (high Andean forest); Mpc: Mosaic of pastures and crops; Lag: Lagoons; Bdl: Bare and degraded lands (deforested areas).
Preprints 94813 g001
Figure 2. Temporal trend between the number of patches (NP in gray bars) and the largest patch index (LPI in black dots) for the mosaic of pastures and crops (Mpc), dense grassland of non-wooded mainland, i.e., paramo (Dgnm), and high Andean forest (Hof). T1-3: Temporality (years 1986, 2000, and 2021).
Figure 2. Temporal trend between the number of patches (NP in gray bars) and the largest patch index (LPI in black dots) for the mosaic of pastures and crops (Mpc), dense grassland of non-wooded mainland, i.e., paramo (Dgnm), and high Andean forest (Hof). T1-3: Temporality (years 1986, 2000, and 2021).
Preprints 94813 g002
Figure 3. Ecological connectivity in the Bijagual landscape. A. Average contribution and standard deviation of the dPCintra, dPCflux, and dPCconnector fractions by cover and year; B. Spatiotemporal contribution of the Hof and Dgnm patches to habitat quality and availability (dPCintra) and ecological flow (dPCflux); C. Hof and Dgnm patches with the highest contribution to ecological connectivity by 2021. Dgnm: Dense grassland of non-wooded mainland (paramo); Hof: High open forest (high Andean forest).
Figure 3. Ecological connectivity in the Bijagual landscape. A. Average contribution and standard deviation of the dPCintra, dPCflux, and dPCconnector fractions by cover and year; B. Spatiotemporal contribution of the Hof and Dgnm patches to habitat quality and availability (dPCintra) and ecological flow (dPCflux); C. Hof and Dgnm patches with the highest contribution to ecological connectivity by 2021. Dgnm: Dense grassland of non-wooded mainland (paramo); Hof: High open forest (high Andean forest).
Preprints 94813 g003
Figure 4. Study area. Geographic location of the Bijagual Massif in Colombia. A. Location of the Boyacá department in Colombia and South America; B. Location of the Bijagual Massif in the Boyacá Department; C. Bijagual Massif.
Figure 4. Study area. Geographic location of the Bijagual Massif in Colombia. A. Location of the Boyacá department in Colombia and South America; B. Location of the Bijagual Massif in the Boyacá Department; C. Bijagual Massif.
Preprints 94813 g004
Table 1. Error matrix, thematic accuracy for 2021. The gray cells are correct assignments; the rest are leak assignments. Dgnm: Dense grassland of non-wooded mainland; Hof: High open forest; Mpc: Mosaic of pastures and crops; Lag: Lagoons; Bdl: Bare and degraded lands.
Table 1. Error matrix, thematic accuracy for 2021. The gray cells are correct assignments; the rest are leak assignments. Dgnm: Dense grassland of non-wooded mainland; Hof: High open forest; Mpc: Mosaic of pastures and crops; Lag: Lagoons; Bdl: Bare and degraded lands.
Classification Dgnm Hof Mpc Lag Bdl Total row User accuracy (%) Commission error (%)
Dgnm 49 4 0 0 0 53 98.55 1.45
Hof 1 140 0 2 0 143 98.91 1.09
Mpc 0 2 71 0 0 73 99.28 0.72
Lag 0 0 0 6 0 6 100 0
Bdl 0 0 0 0 1 1 100 0
Total column 50 146 71 8 1 276
Producer accuracy (%) 99.64 97.83 100 99.28 100 Overall classification accuracy (Kappa) = 0.95
Error of omission (%) 0.36 2.17 0 0.72 0
The residuals in the paramo, forest, pasture, and lagoon columns indicate real coverage that is not included in the map (error of omission, when a pixel has a specific coverage on the terrain and is not assigned to that class on the map); the residuals of the rows were associated with map coverages that do not conform to reality (error of commission, elements that, although they do not belong to a particular class, are found in it).
Table 2. Area occupied per vegetation cover for the years 1986, 2000, and 2021 in Bijagual. Mpc: Mosaic of pastures and crops; Dgnm: Dense grassland of non-wooded mainland; Hof: High open forest; Lag: Lagoons; Bdl: Bare and degraded lands.
Table 2. Area occupied per vegetation cover for the years 1986, 2000, and 2021 in Bijagual. Mpc: Mosaic of pastures and crops; Dgnm: Dense grassland of non-wooded mainland; Hof: High open forest; Lag: Lagoons; Bdl: Bare and degraded lands.
Cover Area (ha) Change from 1986 to 2000 (ha) Change from 2000 to 2021 (ha) Change from 1986 to 2021 (ha)
1986 % 2000 % 2021 %
Mpc 2,811.8 32.8 3,748.5 43.7 4,534.1 52.9 936.6 785.6 1,722.2
Dgnm 3,648.9 42.6 2,400.9 28.0 2,160.0 25.2 -1,248.0 -241.0 -1,489.0
Hof 2,101.3 24.5 2,410.5 28.1 1,861.8 21.7 309.2 -548.7 -239.5
Lag 12.1 0.1 14.1 0.2 13.1 0.2 2.0 -1.0 1.0
Bdl 0.0 0.0 0.0 0.0 5.1 0.1 0.0 0.9 0.9
Total 8,574.1 100.0 8,574.1 100.0 8,574.1 100.0
Table 3. Transition matrix of the vegetation cover in Bijagual from 1986 to 2000 and from 2000 to 2021. Gray boxes correspond to areas (ha) that are maintained between periods. Mpc: Mosaic of pastures and crops; Dgnm: Dense grassland of non-wooded mainland; Hof: High open forest; Lag: Lagoons; Bdl: Bare and degraded lands.
Table 3. Transition matrix of the vegetation cover in Bijagual from 1986 to 2000 and from 2000 to 2021. Gray boxes correspond to areas (ha) that are maintained between periods. Mpc: Mosaic of pastures and crops; Dgnm: Dense grassland of non-wooded mainland; Hof: High open forest; Lag: Lagoons; Bdl: Bare and degraded lands.
Change of orientation or use between 1986 and 2000
Coverage Hof (ha) Dgnm (ha) Mpc (ha) Lag (ha) Bdl (ha) Total
(year 1986)
Hof 1,875.8 130.6 92.1 2.9 0 2,101.3
Dgnm 342.5 2,123.6 1,182.8 0 0 3,648.9
Mpc 191.4 146.8 2,473.6 0 0 2,811.8
Lag 0.9 0 0 11.2 0 12.1
Total (year 2000) 2,410.5 2,401.0 3,748.5 14.1 0.0 8,574.1
Change of orientation or use between 2000 and 2021
Coverage Hof (ha) Dgnm (ha) Mpc (ha) Lag (ha) Bdl (ha) Total
(year 2000)
Hof 1,441.3 447.3 521.9 0 0 2,410.5
Dgnm 257.9 1,402.5 740.6 0 0 2,401.0
Mpc 161.6 310.3 3,271.6 0 5.1 3,748.5
Lag 1.0 0 0 13.1 0 14.1
Total (year 2021) 1,861.8 2,160.0 4,534.1 13.1 5.1 8,574.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated