SciELO - Scientific Electronic Library Online

 
vol.43 issue2Spatial information gaps in continental terrestrial mammal richness from ColombiaEffect of the matrix-edge-forest interior gradient on the phyllostomid bats assemblage in sub-Andean forest fragments author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Caldasia

Print version ISSN 0366-5232On-line version ISSN 2357-3759

Caldasia vol.43 no.2 Bogotá July/Dec. 2021  Epub Apr 07, 2022

https://doi.org/10.15446/caldasia.v43n2.84845 

Ecología

Effect of landscape modification on primate assemblages of the Magdalena River Valley, Colombia

Efecto de la modificación del paisaje sobre el ensamblaje de primates en el valle del río Magdalena, Colombia

Néstor Roncando-Duque1  * 
http://orcid.org/0000-0001-8575-8272

1 Facultad de Ciencias Agropecuarias - Universidad Nacional de Colombia. nroncanciod@gmail.com


ABSTRACT

Population and community dynamics are affected by local landscape disturbances, and the answer of each species to these changes could differentially affect species survival. Natural land cover in the Magdalena River valley, Colombia, has been reduced by almost 70 % and fragmented. There inhabit isolated and reduced populations of Alouatta seniculus, Ateles hybridus, Cebus versicolor, and Sagui-nus leucopus (Primates). To understand the effect of fragmentation over these species, the relationship between altitude, landscape and class metrics, and plant structure and diversity with the population density of these four primates in 20 localities were assessed. Additionally, the relationship between the density of A. hybridus, assumed as a better competitor, and the density of the other primates was assessed. Mean shape index was the most related variable and affected positively the density of A. hybridus, C. versicolor, and S. leucopus and negatively that of A. seniculus. Weight class index was the second more important variable, the four primates increased their density where more mature is the plant cover. Diameter at breast height was positively related with the density of A. hybridus and C. versicolor, negatively related for A. seniculus and not related to S. leucopus density. Tree density was related negatively with the density of A. seniculus and positively with S. leucopus. A. hybridus was related to the density of the other primates, positively to A. seniculus and negatively to the other two species.

Keywords: New World monkeys; relative abundance; habitat transformation

RESUMEN

Las dinámicas poblacionales y de las comunidades son afectadas por los disturbios en los paisajes locales, y las respuestas de cada especie a esos cambios podrían afectar diferencialmente su sobrevivencia. La cobertura vegetal natural en el valle del río Magdalena en Colombia ha sido reducida casi en un 70 % y está fragmentada. En este paisaje habitan en poblaciones reducidas y aisladas Alouatta seniculus, Ateles hybridus, Cebus versicolor y Saguinus leucopus. Para evaluar el efecto de la fragmentación sobre estas especies se evaluó la relación de la altitud, métricas de paisaje y clase, y estructura y diversidad de plantas con la densidad poblacional de los cuatro primates en 20 localidades. Adicionalmente, se evaluó la relación entre la densidad de A. hybridus, asumiéndolo como mejor competidor, con la de los otros primates. El índice medio de forma fue la variable más relacionada con la densidad y afectó positivamente a A. hybridus, C. versicolor y S. leucopus y negativamente a A. seniculus. El índice ponderado de coberturas fue la segunda, los cuatro primates incrementan la densidad cuando más madura es la cobertura vegetal. El diámetro a la altura del pecho estuvo relacionado positivamente con la densidad de A. hybridus y C. versicolor, negativamente con A. seniculus y no tuvo efecto para S. leucopus. La densidad de árboles estuvo relacionada negativamente con la densidad de A. seniculus y positivamente con la de S. leucopus. La densidad de A. hybridus resultó relacionada con la de los otros primates, positivamente con A. seniculus y negativamente con las otras dos especies.

Palabras clave: Primates del nuevo mundo; abundancia relativa; transformación de hábitat

INTRODUCTION

Anthropogenic disturbances such as logging and forest clearing constitute major threats to biodiversity (Michalski and Peres 2005). When a disturbance changes environmental conditions, the natural community can suffer profound, immediate, and cumulative effects (Allesina et al. 2006). However, as different species have different ecological attributes-their scale of movement, life history, longevity, and what constitutes its habitat-that determine their niche and their capacity to survive in a modified landscape. Thus, some become more abundant or expand their range, while others may drop and even become locally extinct (Boyle and Smith 2010, Chapman et al. 2010). The loss of a species or a change in its abundance, particularly for those that interact with many other species, can have a considerable effect on ecological processes (Lindenmayer and Fischer 2013). It is, therefore, important to consider multiple anthropogenic factors when assessing the synergistic effects of human activities, in order to understand how species, respond in fragmented landscapes (Benchimol and Peres 2013).

The west margin of the Magdalena River valley and east slope of the Cordillera Central have a history around 150 years of intensive transformation, and the natural land cover has been reduced by almost 70 %, transforming it into a highly fragmented landscape (IDEAM 2009). The floral association in the forest remnants have been impacted at different levels and different ways, affecting successional processes.

Four primates are found in this region: Red howler monkeys (Alouatta seniculus Linnaeus, 1766), brown spider monkeys (Ateles hybridus Geoffroy, 1829), white fronted capuchins (Cebus versicolor Pucheran, 1845) and white footed tamarins (Saguinus leucopus GUnther, 1877). A. seniculus is an atelid that weights around 7 kg that inhabits a broad range of forest covers from transitional vegetation until riparian and mature forest up to 3200 masl. It is a facultative leaf-eater monkey when the fruit is scarce (Link et al. c2021c). A. hybridus is an atelid primate too around that weights 7.5 kg. It is dependent of mature forest, but some groups could be found in fragmented forests. It is highly fru-givorous and feed largely on ripe fleshy fruits, which comprise approximately 80 % of its diet (Link et al. c2020b). C. versicolor is a cebid primate that weights around 2.5 kg that inhabits tropical lowland forests and seasonally flooded forests. It has an omnivorous diet, eating predominantly ripe fruits and invertebrates, but also other vegetative parts of plants, vertebrate prey and eggs (Link et al. c2021d). S. leucopus is a small (460 g) callitrichid that inhabits from transitional to mature, mainly humid, and subhumid forests. Its diet is based on fruits and is supplemented with invertebrates, bark, and flowers (Roncancio et al. 2012). and S. leucopus are endemic to Colombia. C. versicolor is Endangered (Link et al. c2021d) and S. leucopus is Vulnerable (Link et al. c2021a). A. hybridus is categorized as Critically Endangered (Link et al. c2020b). A. seniculus is categorized as Least Concern even though its populations are decreasing (Link et al. c2021c). The main threat to these taxa is the reduction and fragmentation of their habitat resulting from the combined effects of selective extraction, hunting and trade, and the lack of no-take protected areas as National Natural Parks (Londoño and Gómez-Posada 2010, Roncancio et al. 2011, Defler et al. 2013).

The responses of these primates to disturbances, and their specific tolerance thresholds to the resulting landscape configuration are, however, not well understood (French and Smith 2005, Ewers and Didham 2007) due to the lack, among other things, of a landscape model grounded on a robust mathematical approach (Reckhow 1990, Wade 2001, Fahrig 2003, Arroyo-Rodríguez et al. 2013). In this study, I estimated the relationship of different traits such as size, shape, isolation distance and land cover composition, plant structure, and plant diversity, with the population density of these four primates. I also estimated the relationship between the population densities of these primates under a theoretical framework of compensation density (Peres and Dolman 2000). I assumed that A. hybridus is a better competitor than other species due to its large size and specialized frugivorous diet (Defler 2010). My hypothesis was that the changes in the population density of primates could be explained by changes in size, shape, isolation distance, composition of plant covers and plant structure and diversity, although, differentially in type and intensity and affected by the interspecific interactions of the assembly.

MATERIALS AND METHODS

Study area

I conducted fieldwork in 20 localities (local landscape), randomly selected with respect to the primate presence. In other words, the presence of the primates was not a criterion to select the localities (WGS 84 - decimal degrees: Amaní -74.88611W, 5.56083N; Berlin -74.57354W, 6.52191N; Cerro Las Campanas -74.45416W, 8.23381N; Cerros de Corcovado -74.46314W, 8.46191N; Charca de Guarinocito -74.73472W, 5.33972N; EL Pajuil -74.69466W, 5.88033N; Falan -74.97000W, 5.11889N; Fresno -75.01475W, 5.20222N; Guamo 2 -74.69180W, 5.86057N; Guamo1 -74.70102W, 5.86893N; La Florida -74.19829W, 8.57968N; La Pedrera -74.95722, 5.57611; La Primavera -74.96667W, 5.51667N; Las Margaritas -74.59761W, 6.56341N; Mariquita -74.94606, 5.22742N; Playas -74.92331W, 6.26989N; Porce III -75.18909W, 6.83280N; Río Manso -74.77361W, 5.67556N; Venecia -74.83861W, 5.64278N; Yalí -74.87281, 6.67851N). These areas are located in the middle Magdalena River valley in Colombia, distributed from the north of the department of Tolima to the south of the department of Bolivar, extending through the eastern parts of the departments of Caldas and Antioquia (Fig. 1). The altitude of the locations ranged from 190 to almost 1400 m (Table 1), with mean temperatures of 19 to 27 °C and mean annual precipitation ranging from 1800 to 4200 mm (Hijmans et al. 2005). The region is approximately 48 000 km2, comprising crops (23 %), pasture (22 %) and natural vegetation (32 %) (the remaining 23 % is for other land use). Natural vegetation includes dense forest fragments of less than 0.1 to 300 and 6000 km2 (mean = 60 km2 ±1640), with shapes that tend to be irregular (Mean Shape Index = 1.8) and medium connectivity (interspersion juxtaposition index = 52, range: 0 < IJI ≤100) (Statistics estimated for author from IDEAM 2009).

Figure 1 Location of the sampling localities along the Magdalena River valley, Colombia 

Table 1 Local features, landscape and class metrics selected as explanatory variables. 

Primate population density

Sampling was done between 2005 and 2011. To estimate primate population densities, I used distance sampling along linear transects (Buckland et al. 2001) and censuses in two localities where I could assume that complete counts were possible. With this method, transects can be resampled to increase the sample size to calculate encounter rate, and individual recognition is not necessary. The number and length of transects varied depending on the extension, shape and topography of each site (Table 2). Transects were distributed throughout each site to obtain a representative density estimation. In two localities, Guarinocito and Guamo 1, population density was estimated using censuses, because these localities are fragments without structural connectivity of forest, of 0.21 and 0.1 km2, respectively, and are elongated. I therefore, assumed that it was possible to evaluate the entire area of interest and count all its primates individually. Population densities of A. seniculus, A. hybridus and C. versicolor in Florida, A. seniculus in Pedrera and Rio Manso, and A. hybridus and C. versicolor in Margaritas, were estimated using the detection function obtained in similar localities of this study, because the number of detections were too low (<10) in these localities to enable a robust analysis in Distance. Some population densities were taken from secondary information (Santamaria et al. 2007, Rojas-V. 2009). I assumed that the results are comparable given that they used de same methodology.

Table 2 Sampling design and sampling effort to estimate the population density of the four primate species using distance sampling with line transects. Guarinocito and Guamo 1 are not in this table because no distance sampling line transects were used there. 

1 Rojas-V. (2009); Santamaria et al. (2007)

Local features and landscape and class metrics of each site

Initially, I defined a local landscape for each locality through a circle that was the smallest perimeter containing the transects used to estimate primate population densities (Fig. 2). For Guarinocito and Guamo 1, the entire fragments were contained in the smallest perimeter. This circle was overlaid with a raster of land cover (IDEAM 2009), a raster of mean annual temperature and mean annual precipitation from BIOCLIM (Hijmans et al. 2005), and with a Digital Elevation Model (USGS 2000). For each locality I measured three local features, two landscape metrics and twelve class metrics (in the class where the population density was estimated). Local features included: mean altitude, mean annual temperature, and mean annual precipitation. Landscape metrics included: Shannon index (SI) and a weighted class index (WCI). I developed the WCI to integrate an index that allowed to measure the land cover composition and structure, not only the land cover structure as SI (Fahrig et al. 2011). To estimate the WCI, I qualified class types as follows: Dense forest = 8, Open forest = 7, Fragmented forest = 6, Scrubland = 5, Grassland (natural) = 4, Crops = 3, Pasture = 2, and wetland, open soil and urbanized area = 1. Then, I multiplied the value of each class type by its relative size in the local landscape. Finally, I added the results of the multiplication to obtain the index values. In this way, the value for a local landscape with predominantly dense forest, was close to eight, while for one with more pasture and crops, the values were between two and three. Class metrics, for the class in the local landscape where the population density estimation was done, were as follows: mean core area (MCA), total core area (TCA), total core area index (TCAI), core area density (CAD), interspersion juxtaposition index (IJI), area weighted mean shape index (AWMSI), mean shape index (MSI), mean patch fractal dimension (MPFD), area weighted mean patch fractal dimension (AWMPFD), total edge (TE), mean patch size (MPS) and core area (CA). All analyses were done using ArcGis 10.5, Patch Grid tool.

Figure 2 Example of the local landscapes delimitation to estimate the local features, landscape and class metrics. A = Playas, B = Berlín. 

Plant structure and diversity

In the forest where the primate densities were estimated of each locality, six or seven rectangular plots of 0.20 ha (50 m x 4 m) were randomly placed. I stopped sampling when the standard deviation of the diameter breast height (DBH) was stable. I took samples for species identification, and determinations were carried out in the herbaria of the Antioquia -HUA-, and Caldas -FAUC- universities. Mean DBH and height, tree density, and the Inverse Simpson Index (ISI) were estimated for each locality as explanatory variables.

Data Analyses

Population density. I analysed data using Distance 7.1 software (Thomas et al. 2010). Analyses were performed separately for each locality. Detection functions were selected according to the fit between the frequency distribution of detection distances and theoretical models (Key functions and adjustment terms). The models tested were: Half-normal (Cosine, Hermite polynomial), Uniform (Cosine, simple polynomial), and Hazard rate (Cosine, simple polynomial). I chose the best fitting model according to the Akaike Information Criterion. Population density variance (CV) was empirically calculated as the sum of the sampling variances of encounter rate, estimates of detection probability, and group size. I used the expected group size estimated from the regression between group size and detection probability to estimate the density (Buckland et al. 2001). In Guarinocito and Guamo 1 where I used censuses, the density was estimated divided the number of individuals over the size of the fragments. Distance sampling is one of the most often used methods to obtain precise and unbiased estimations with a given sampling effort, so I assumed that the results obtained with distance sampling and the censuses are comparable (Cassey and McArdle 1999, Norvell et al. 2003).

Relationship of the explanatory variables in primate density. Because the vegetation was sampled in only fifteen out of the 20 localities, I analysed the relationship of the density of primates and the local features, landscape, and class metrics of the localities independent of the plant structure and diversity. First, I performed analyses of association between the explanatory variables in order to eliminate collinearity (Autocorrelation) (Cush-man et al. 2008). I used the Spearman coefficient of correlation to evaluate this association (data did not distribute normally, Shapiro-Wilk test: W = 0.90, P = 0.04). When two variables were correlated (Spearman rank correlation coefficient: P-rs < 0.05), one of them was eliminated (Sokal and Rohlf 1995), considering the precision of the measurement. However, the main criterion was mathematical. Finally, four out of the 17 local features, landscape and class metrics were considered (altitude, interspersion juxtaposition index, mean shape index and weighted cover index) and, likewise, only the density of the Variegated spider monkey out the four primates and the four plant structure and diversity variables were considered as explanatory variables.

Given that the population density is a variable expressed in area, I assumed that it has Poisson probability distribution and, thus, the relationships were evaluated using Poisson regressions. This was done using a Bayesian approach and using an uninformative prior distribution to the precision of the explanatory variable effects, and an uninformative prior distribution to the intercept (alpha) and slopes (beta). I assumed that the posterior distributions of the intercept and the slope of each variable had a normal distribution (McCarthy 2007, Pfeiffer et al. 2008). I did not consider the interaction (multiplicative effects) between the explanatory variables. In order to select the best model, the Deviance Information Criterion (DIC) was used. Initially, I ran the model with all of the possible explanatory variables and estimated the DIC. Later, I extracted one variable, ran the model and estimated de DIC again. If the DIC was higher, I returned the variable and extracted another. If the DIC resulted lower, I removed definitely that variable and extracted another, until I got the model with the lowest DIC. Estimation of the intercept and slopes was done using Markov chains with 100 000 iterations, and taking account from the 10 001 iteration to the final estimation. The DIC was estimated with 100 000 additional iterations. The analyses were performed using OpenBugs 3.2.2 software (Lunn et al. 2000).

RESULTS

Population densities

Saguinus leucopus was found in 20 localities with densities of 6.0 to 149 ind/km2. A. hybridus was found in seven localities with densities of 0.62 to 58 ind/km2. C. versicolor was found in six localities with densities of 1.04 to 103 ind/km2, and A. seniculus was found in nine localities with densities of 0.21 to 75 ind/km2 (Table 3).

Table 3 Population density (individuals /km2) of the four primate species recorded in the surveys. D= Population density, Cl= Confidence interval 95%, CV= Coefficient of variation. 

1 Rojas-V. (2009); 2Santamaría et al. (2007)

Local features, landscape, and class metrics

The localities were at elevations between 191 m and 1385 m. Corcovado, Campanas and Margaritas had the largest area of dense forest while Guarinocito and Fresno were the localities with largest areas of crops and pastures (Table 1; WCI). The localities with the more aggregated forests were Fresno, Río Manso, Guamo 1, Guamo 2, and Corcovado (One fragment with all pixels together), while Margaritas had the most disaggregated forest patches (Table 1; IJI). The localities with the most regular shaped forests were Guamo 1, Guamo 2 and Pajuil, while the forests of Pedrera and Primavera were the most elongated (Table 1; MSI).

Regarding plant structure and diversity, Guarinocito showed the highest mean DBH following by Corcovado and Río Manso. Playas had the smallest mean DBH. Río Manso, Amaní and Berlín had the tallest trees. The highest tree densities were found in Primavera and Berlín, and the lowest in Porce III (Table 4). The localities with the highest tree diversity were Río Manso and Berlin, while Porce III, Playas and Guarinocito showed the lowest tree diversity values.

Table 4 Plant structure and diversity for each locality. DBH = Diameter breast height 

1 Rojas-V. (2009); 2Santamaria et al. (2007)

Relationship of the explanatory variables in primate densities

Higher population densities of S. leucopus were recorded in localities with elongated, denser, and taller forests, in localities at lower elevations and with higher structural connectivity. These results also show that S. leucopus densities decrease when the density of A. hybridus increase (Table 5). Population densities of C. versicolor and A. hybridus tended to be higher when there was less connectivity. Regarding with plant structure and diversity, these results suggest that C. versicolor densities are positively affected only by tree DBH, while A. hybridus populations tend to increase in forests with higher DBH values, lower tree-height, and higher plant diversity (Table 5).

Table 5 Relationship between the explanatory variables and the density of the studied species. WCI = Weighted core index; IJI = Interspersion juxtaposition index; MSI = Mean shape index; TD = Tree density; DBH = Diameter at breast height; SI = Shannon index. 

**Main effect; *Secondary effect; Blank = no significant effect

Population densities of A. seniculus were higher with denser forests, and lower in more elongated forests. Contrary to the pattern found for S. leucopus and C. versicolor, the population densities of A. seniculus increase when the density of A. hybridus increase and decrease when the lower DBH and tree density increase (Table 5).

DISCUSSION

These results show that the population densities of the four studied species are very variable among different localities, from local extinction, for example to A. hybridus in all localities in Caldas, to some of the highest population densities estimated for the taxon or congeneric species (Crockett et al. 1987, Chapman and Balcomb 1998, Gómez-Posada et al. 2007, 2009, Roncancio and Gómez-Posada 2009, Gómez-Posada et al. 2010, Williams-Guillén et al. 2013, De Luna and Link 2018). This indicates that different zones in the Magdalena River valley have been affected in different ways. Differences in the local occupation between species indicate that S. leucopus in this landscape could resists better the degradation, reduction and fragmentation of its habitat and synergic threats such as trade, given that, in spite of the large variation in its densities, it was found in all 20 localities. The other three species, in contrast, show less tolerance to changes in their habitats, and face a local extinction probability close to 60 % even though there are remnant forests, based in proportion of the occupied localities to each primate species.

Results obtained for C. versicolor and A. hybridus are consistent with the patterns found for Cebus capucinus (Linnaeus, 1758) and Ateles geoffroyi (Kuhl, 1820) in Nicaragua and with Cebus libidinosus (Spix, 1823) in Bolivia (Pyritz et al. 2010) where they have suffered local extinctions. Usually, A. seniculus and other howlers are recognized as being able to tolerate relatively higher levels of habitat fragmentation, as in many areas it is the only primate taxon to persist and in some it can be found at very high densities (Crockett et al. 1987, Chapman and Balcomb 1998, Gómez-Posada et al. 2007, 2009, Roncancio and Gómez-Posada 2009, Gómez-Posada et al. 2010, Williams-Guillén et al. 2013, De Luna and Link 2018). Red howlers have been locally extirpated in a number of localities in the Magdalena valley, and their densities are very variable in the localities where they are still present. A study of Alouatta guariba (Humboldt, 1812) carried out in southern Brazil showed a similar situation to that found for red howler populations in this study (Silva and Bicca-Marques 2013). The population depletion of howler monkeys could become driven by synergic threats and the degradation of the ecological conditions in the forest fragments.

Different features, landscape and class metrics are related differently the population densities of the four studied primate species, even though forest shape had the most significant effect on population density for all species. This situation was also found in Mexico, where the occupancy of all primates was negatively affected by shape irregularity (Arroyo-Rodríguez et al. 2008). In this study the effect of the shape was not the same for the four primates. While the population density of A. seniculus tended to decrease in elongated forests, the densities of the others three primates tended to increase (Lenz et al. 2014). High densities, however, do not always mean high abundance or better conditions for the population. High densities in forest remnants have been explained as a crowding effect (Chiarello and de Melo 2001, Peres 2001, Baranga 2004, Martins 2005, Rode et al. 2006). If the crowding is the first effect that suffer these primates, probably, the howler monkeys become affected later than the other species given that it could use some types of forest land cover that the other species do not use and in this way the population is diluted in a more extended area (Crockett 1998, Gómez-Posada and Londoño 2012). On the other hand, depending on the moment of the fragmentation process and the regeneration dynamics, early fragmentation, reversal isolation by regeneration in transitional crops, the smaller or more mobile species could or must use the elongated fragments to satisfy their ecological requirements, increasing their relative abundance, while the howler just need a small home range in the core areas (Matthews 2009, Spehar et al. 2010, Palma et al. 2011, Alba-Mejia et al. 2013). In this sense the fragment shape is an important factor that determines population size in fragmented landscapes, thus, habitat restoration efforts may be more effective if they focus on saving the forest fragments with less elongated shapes and keep a minimal width in the corridors (Ewers and Didham 2007).

Forest cover (WCI) was the second most important factor, when WCI was higher the population densities of all four species tend to be higher. Similarly, forest land cover was one of the main predictors of species persistence in isolated forests in a continental-scale analysis of neotropical primates (Benchimol and Peres 2013) and in two primate-assemblages in French Guiana and Ecuadorian Amazonia (Youlatos 2004). Forest land cover could be a determinant factor in the survival of these species (Sorensen and Fedigan 2000). Additionally, the general structure of the land cover landscape, associated with the matrix has been recognized as a factor that determines the presence of some primates in a fragment (da Silva et al. 2015).

Regarding plant structure, DBH was the main determining factor but, while the howler densities tend to be lower in forest with higher densities of large-trunked trees, densities of spider monkeys and capuchins tends to higher. In contrast, Alouatta palliata (Gray, 1849) has been shown to prefer fragments with a greater density of large trees (DBH > 60 cm) (Arroyo-Rodríguez et al. 2008). The equitability in the plant diversity was important for A. hybridus only, highlighting the high ecological requirements to this species and determine the necessity of focal intervention as specific plant enrichment to habitat restoration (Boyle and Smith 2010).

There is broad evidence that some types of agro-ecosystems in a fragmented landscape could be useful as temporary or additional habitat for different taxa (Estrada et al. 2006, Guzmán et al. 2016). However, I found that the primates in the Magdalena River valley are concentrated in zones with considerable areas of mature forest (WCI), including S. leucopus, even though tamarins and marmosets have preferences for regenerating forest (Snowdon and Soini 1988, Rylands 1996, Garber 1998, Defler 2010). This pattern is probably related to the ecological stability and resilience provided by mature forests. Thus, the biodiversity conservation at a landscape scale requires considerations regarding a minimum proportion of dense forest in the mosaic.

The densities of A. hybridus were significantly related to those of the other three primates, but, while howlers tend to be found in higher densities in forests with higher densities of Ateles, the other two species tend to be less abundant. The reason may be related to hunting pressure Howler monkeys come second to spider monkeys as preferred game by hunters (Freese et al. 1982, Rodríguez-Luna et al. 1996, Sorensen and Fedigan 2000, da Silva et al. 2005). So, it could be possible that the populations of A. seniculus will not start to decline by hunt until the density of A. hybridus is impacted by overhunting.

Ateles hybridus, by size, could be a better competitor than Cebus versicolor and Saguinus leucopus. The three species are overlapping in some niche dimensions (ecological requirements), so, by a compensation density process (MacArthur et al. 1972), the abundance of A. hybridus could regulate the density of the smaller primates. As a result, when the density of the bigger primates decreases, the density of the smaller primates increases (Peres and Dolman 2000, Hilário and Ferrari 2015), as has been observed in fragments where there were no spider monkeys, howlers or capuchins, and the density of the tamarins was of over 100 ind/km2. Thus, the high densities of S. leucopus reflect a very disrupted assemblage of primates, and I must not interpret them as a good indicator for the whole system.

I therefore suggest the use of A. hybridus as a landscape species to play the role as umbrella, keystone and maybe flagship surrogate to guide the priorities in the management of the humid and sub-humid forests of the Magdalena River valley (Sanderson et al. 2002). A. hybridus was the most affected species and by more factors (local vulnerability) (Rimbach et al. 2013) and it is a keystone species in this system (Christianou and Ebenman 2005). Thus, consolidating sustainable conservation scenarios for this primate, will satisfy the requirements for the other species, as a good representative of the ecological integrity of the whole system in a landscape approach.

The landscape approach to the conservation management usually demands the mosaic design that determines the minimal proportion and type of the core areas, their shape and aggregation and landscape composition. The landscape species approach based in the requirements or in the landscape structure effects allow us to plan and monitor the management actions in more cost-effective way.

However, landscape structure needs to be measured at the appropriate scale (Arroyo-Rodriguez and Fahrig 2014). The majority of the studies, including this, quantified landscape predictors within a single spatial scale, potentially missing significant primate-landscape responses (Galán-Acedo et al. 2019). Thus, studies on landscape effects on primates need to use a multiscale approach to ensure that landscape-species associations are correctly evaluated (Arroyo-Rodriguez and Fahrig 2014, Jackson and Fahrig 2015, Galán-Acedo et al. 2019). This study was done as patch-landscape analysis, and the scales were adjusted to each local landscape for a relative criteria based in the inferential area to primate density estimation. I think that standard measures of scales or direct scales are not applicable since the size of the inferential areas was highly variable, but a multiscale stratified analyses, with enough sampling sites, could be addressed in regions with similar land use patterns such as Caldas western region, Antioquia southeast and Serranía de San Lucas, that allow the use of standard scales, based on ecological criteria, to determine the optimal scale for each species to each region (Jackson and Fahrig 2012, 2015).

ACKNOWLEDGEMENTS

I mainly acknowledge to Wildlife Conservation Society Colombian program by the technical and economical support. Additionally, I acknowledge to European Association of Zoo and Aquarium -EAZA- and Colombia Association of Zoos and Aquarium -ACOPAZOA, Corporacion autonoma regional de Caldas -Corpocaldas, Environmental Ministry and Margot Mash Biodiversity Foundation whom supported with resources in different moments of this data collection. This work would not have been possible without the participation of almost 30 students a recent graduate biologist of the different universities of Colombia and Bolivia. Special gratitude to Diana Carolina Guzman Caro and Olena Kovalenko and the anonymous editor by the English revision and the relevant advises about this paper.

LITERATURE CITED

Alba-Mejia L, Caillaud D, Montenegro OL, Sánchez-Palomino P, Crofoot MC. 2013. Spatiotemporal interactions among three neighboring groups of free-ranging whitefooted tamarins (Saguinus leucopus) in Colombia. Int. J. Primatol. 34:1281- 1297. doi: https://doi.org/10.1007/s10764-013-9740-6Links ]

Allesina S, Bodini A, Bondavalli C. 2006. Secondary extinctions in ecological networks: bottlenecks unveiled. Ecol. Model. 194(1-3): 150-161. doi: https://doi.org/10.1016/j.ecolmodel.2005.10.016Links ]

Arroyo-Rodríguez V, Mandujano S, Benítez-Malvido J. 2008. Landscape attributes affecting patch occupancy by howler monkeys (Alouatta palliata mexicana) at Los Tuxtlas, Mexico. Am. J. Primatol. 70(1): 69-77. doi: https://doi.org/10.1002/ajp.20458Links ]

Arroyo-Rodríguez V, Cuesta-del Moral E, Mandujano S, Chapman CA, Reyna-Hurtado R, Fahrig L. 2013. Assessing habitat fragmentation effects on primates: The importance of evaluating questions at the correct scale. In: Marsh LK, Chapman CA, editors. Primates in Fragments. New York: Springer. p. 13-28. [ Links ]

Arroyo-Rodríguez V, Fahrig L. 2014. Why is a landscape perspective important in studies of primates? Am. J. Primatol. 76(10): 901-909. doi: https://doi.org/10.1002/ajp.22282Links ]

Baranga D. 2004. Forest fragmentation and primates' survival status in non-reserved forests of the "Kampala area", Uganda. Afr. J. Ecol. 42(1): 70-77. doi: https://doi.org/10.1111/j.1365-2028.2004.00464.xLinks ]

Benchimol M, Peres CA. 2013. Anthropogenic modulators of species-area relationships in Neotropical primates: a continental-scale analysis of fragmented forest landscapes. Divers. Distrib. 19(11): 1339-1352. doi: https://doi.org/10.1111/ddi.12111Links ]

Boyle SA, Smith AT. 2010. Can landscape and species characteristics predict primate presence in forest fragments in the Brazilian Amazon? Biol. Conserv. 143(5): 1134-1143. doi: https://doi.org/10.1016/j.biocon.2010.02.008Links ]

Buckland ST, Anderson DR, Burnham KP, Laake JL, Borchers DL, Thomas LJ. 2001. An Introduction to Distance Sampling: Estimating Abundance of Biological Populations. Oxford, UK: Oxford University Press. [ Links ]

Cassey P, McArdle BH. 1999. An assessment of distance sampling techniques for estimating animal abundance. Environmetrics. 10(3): 261-278. doi: https://doi.org/10.1002/(SICI)1099-095X(199905/06)10:3<261::AID-ENV351>3.0.CO;2-OLinks ]

Chapman CA, Balcomb SR. 1998. Population characteristics of howlers: ecological conditions or group history. Int. J. Primatol. 19: 385-403. doi: https://doi.org/10.1023/A:1020352220487Links ]

Chapman CA, Struhsaker TT, Skorupa JP, Snaith TV, Rothman JM. 2010. Understanding long-term primate community dynamics: implications of forest change. Ecol. Appl. 20(1): 179-191. doi: https://doi.org/10.1890/09-0128.1Links ]

Chiarello AG, de Melo FR. 2001. Primate population densities and sizes in Atlantic Forest remnants of northern Espírito Santo, Brazil. Int. J. Primatol. 22: 379-396. doi: https://doi.org/10.1023/A:1010751527749Links ]

Christianou M, Ebenman B. 2005. Keystone species and vulnerable species in ecological communities: strong or weak interactors? J. Theoret. Biol. 235(1): 95-103. doi: https://doi.org/10.1016/j.jtbi.2004.12.022Links ]

Crockett CM. 1998. Conservation biology of the genus Alouatta. Int. J. Primatol. 19: 549-578. doi: https://doi.org/10.1023/A:1020316607284Links ]

Crockett CM, Eisenberg JF, Smuts BB, Cheney DL, Seyfarth RM, Wrangham RW, Struhsaker TT, editors. 1987. Primate Societies. Chicago: University of Chicago Press. [ Links ]

Cushman SA, McGarigal K, Neel MC. 2008. Parsimony in landscape metrics: strength, universality, and consistency. Ecol. Indic. 8(5): 691-703. doi: https://doi.org/10.1016/j.ecolind.2007.12.002Links ]

Galán-Acedo C, Arroyo-Rodríguez V, Cudney-Valenzuela SJ, Fahrig L. 2019. A global assessment of primate responses to landscape structure. Biol. Rev. 94(5)1605-1618. doi: https://doi.org/10.1111/brv.12517Links ]

da Silva MNF, Shepard GH, Douglas WY. 2005. Conservation implications of primate hunting practices among the Matsigenka of Manu National Park. Neotrop. Primates 13(2): 31-36. doi: https://doi.org/10.1896/1413-4705.13.2.31Links ]

da Silva LG, Ribeiro MC, Hasui É, da Costa CA, da Cunha RGT. 2015. Patch size, functional isolation, visibility and matrix permeability influences Neotropical primate occurrence within highly fragmented landscapes. PloS One. 10: e0114025. doi: https://doi.org/10.1371/journal.pone.0114025Links ]

de Luna AG, Link A. 2018. Distribution, population density and conservation of the critically endangered brown spider monkey (Ateles hybridus) and other primates of the inter-Andean forests of Colombia. Biodiv. Conserv. 27: 3469-35H. doi: https://doi.org/10.1007/s10531-018-1611-1Links ]

Defler TR. 2010. Historia Natural de Los Primates Colombianos. Bogotá: Universidad Nacional de Colombia, Facultad de Ciencias. Departamento de Biología. [ Links ]

Defler TR, Stevenson PR, Bueno ML, Guzmán D, editors. 2013. Primates Colombianos en Peligro de Extinción. Bogotá: Asociación Primatológica Colombiana. [ Links ]

Estrada A, Saenz J, Harvey C, Naranjo E, Muñoz D, Rosales-Meda M. 2006. Primates in agroecosystems: conservation value of some agricultural practices in Mesoamerican landscapes. In: Estrada A, Garber PA, Pavelka MSM, Luecke L, editors. New Perspectives in the Study of Mesoamerican Primates. New York: Springer . p. 437-470. [ Links ]

Ewers RM, Didham RK. 2007. Habitat fragmentation: panchreston or paradigm? Trends Ecol. Evol. 22(10): 511-511. doi: https://doi.org/10.1016/j.tree.2007.06.004Links ]

Fahrig L. 2003. Effects of habitat fragmentation on biodiversity. Ann. Rev. Ecol. Evol. Syst. 34: 487-515. doi: https://doi.org/10.1146/annurev.ecolsys.34.011802.132419Links ]

Fahrig L, Baudry J, Brotons L, Burel FG, Crist TO, Fuller RJ, Sirami C, Siriwardena GM, Martin JL. 2011. Functional landscape heterogeneity and animal biodiversity in agricultural landscapes. Ecol. Lett. 14(2): 101-112. doi: https://doi.org/10.1111/j.1461-0248.2010.01559.xLinks ]

Freese CH, Heltne PG, Napoleon CR, Whitesides G. 1982. Patterns and determinants of monkey densities in Peru and Bolivia, with notes on distributions. Int. J. Primatol. 3: 53-90. doi: https://doi.org/10.1007/BF02693490Links ]

French AR, Smith TB. 2005. Importance of body size in determining dominance hierarchies among diverse tropical frugivores. Biotropica. 37(1): 96-101. doi: https://doi.org/10.1111/j.1744-7429.2005.04051.xLinks ]

Garber PA. 1998. Within- and between-site variability in moustached tamarin (Saguinus mystax) positional behavior during food procurement. In: Strasser E, Fleagle JG, Rosenberger AL, McHenry HM, editors. Primate Locomotion. New York: Springer . p. 61-78. [ Links ]

Gómez-Posada C, Londoño JM. 2012. Alouatta seniculus: density, home range and group structure in a bamboo forest fragment in the Colombian Andes. Folia Primatol. 83: 56-65. doi: https://doi.org/10.1159/000339803Links ]

Gómez-Posada C, Martínez J, Giraldo P. Kattan GH. 2007. Density, habitat use, and ranging patterns of red howler monkeys in a Colombian Andean forest. Neotrop. Primates 14(1): 2-10. doi: https://doi.org/10.1896/044.014.0102Links ]

Gómez-Posada C, Álvarez Z, Giraldo-Chavarriaga P. 2009. Densidad y estatus poblacional de monos aulladores rojos en un guadual, fragmento aislado, La Tebaida, Quindío, Colombia. Univ. Sci. 14(1): 8-15. doi: https://doi.org/10.11144/javeriana.SC14-1.dyepLinks ]

Gómez-Posada C, Roncancio N, Hincapié P, Betancourt A. 2010. Density and group composition in three populations of red howler monkey (Alouatta seniculus) in Valle and Cauca, Colombia. Bol. Cient. Mus. Hist. Nat. Univ. Caldas. 14: 79-91. [ Links ]

Guzmán A, Link A, Castillo JA, Botero JE. 2016. Agroecosystems and primate conservation: shade coffee as potential habitat for the conservation of Andean night monkeys in the northern Andes. Agr. Ecosyst. Environ. 215: 57-67. doi: https://doi.org/10.1016/j.agee.2015.09.002Links ]

Jackson HB, Fahrig L. 2012. What size is a biologically relevant landscape? Landsc. Ecol. 27(7): 929-941. doi: https://doi.org/10.1007/s10980-012-9757-9Links ]

Jackson HB, Fahrig L. 2015. Are ecologists conducting research at the optimal scale? Glob. Ecol. Biogeogr. 24(1): 52-63. doi: https://doi.org/10.1111/geb.12233Links ]

Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. 2005. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25(15): 1965-1978. doi: https://doi.org/10.1002/joc.1276Links ]

Hilário RR, Ferrari SF. 2015. Dense understory and absence of capuchin monkeys (Sapajus xanthosternos) predict higher density of common marmosets (Callithrix jacchus) in the Brazilian Northeast. Am. J. Primatol. 77(4): 425-433. doi: https://doi.org/10.1002/ajp.22358Links ]

[IDEAM] Instituto de Hidrología, Meteorología y Estudios Ambientales. 2009. Leyenda Nacional de Coberturas de la Tierra. Metodología CORINE Land Cover Adaptada para Colombia Escala 1:100.000. Bogotá, D. C., Colombia: Instituto de Hidrología, Meteorología y Estudios Ambientales. [ Links ]

Lenz BB, Jack KM, Spironello WR. 2014. Edge effects in the primate community of the biological dynamics of forest fragments project, Amazonas, Brazil. Am. J. Phys. Anthropol. 155(3): 436-446. doi: https://doi.org/10.1002/ajpa.22590Links ]

Lindenmayer DB, Fischer J. 2013. Habitat Fragmentation and Landscape Change: An Ecological and Conservation Synthesis. Washington, DC: Island Press. [ Links ]

Link A, Guzmán-Caro DC, Roncancio N, Mittermeier RA, Rodríguez V. c2021a. Saguinus leucopus (amended version of 2020 assessment). The IUCN Red List of Threatened Species 2021: e.T19819A192550769. [Last accessed: 20 Abr 2021] doi: [Last accessed: 20 Abr 2021] doi: https://dx.doi.org/10.2305/IUCN.UK.2021-1.RLTS.T19819A192550769Links ]

Link A, Urbani B, Stevenson PR. Mittermeier RA. c2020b. Ateles hybridus. The IUCN Red List of Threatened Species web site. [Last accessed: 11 Oct 2020] [Last accessed: 11 Oct 2020] https://dx.doi.org/10.2305/IUCN.UK.2020-2.RLTS.T39961A17929680.en . [ Links ]

Link A, Palacios E, Cortés-Ortiz L, Stevenson PR, Cornejo FM, Mittermeier RA, Shanee S, de la Torre S, Boubli JP, Guzmán-Caro DC, Moscoso P, Urbani B. c2021c. Alouatta seniculus (amended version of 2020 assessment). The IUCN Red List of Threatened Species 2021: e.T70547436A190432393. [Last accessed: 20 Apr 2021] [Last accessed: 20 Apr 2021] https://dx.doi.org/10.2305/IUCN.UK.2021-1.RLTS.T70547436A190432393.enLinks ]

Link A, Boubli JP, Lynch-Alfaro JW. c2021d. Cebus versicolor. The IUCN Red List of Threatened Species 2021: e.T39952A81282279. [Last accessed: 20 Apr 2021] [Last accessed: 20 Apr 2021] https://dx.doi.org/10.2305/IUCN.UK.2021-1.RLTS.T39952A81282279.enLinks ]

Londoño JM, Gómez-Posada C. 2010. Status of the red howler monkey population in the Barbas Canion montane fragment, central mountains, Colombia. Acta Biol. Colomb. 15: 25-36. [ Links ]

Lunn DJ, Thomas A, Best N, Spiegelhalter D. 2000. A Bayesian modelling framework: concepts, structure, and extensibility. Stat. Comput. 10: 325-337. doi: https://doi.org/10.1023/A:1008929526011Links ]

MacArthur RH, Diamond JM, Karr JR. 1972. Density compensation in island faunas. Ecology 53(2): 330-342. doi: https://doi.org/10.2307/1934090Links ]

Martins MM. 2005. Density of primates in four semi-deciduous forest fragments of São Paulo, Brazil. Biodiv. Conserv. 14: 2321-2329. doi: https://doi.org/10.1007/s10531-004-1666-zLinks ]

Matthews LJ. 2009. Activity patterns, home range size, and intergroup encounters in Cebus albifrons support existing models of capuchin socioecology. Int. J. Primatol. 30: 709-728. doi: https://doi.org/10.1007/s10764-009-9370-1Links ]

McCarthy MA. 2007. Bayesian Methods for Ecology. Cambridge, UK: Cambridge University Press. [ Links ]

Michalski F, Peres CA. 2005. Anthropogenic determinants of primate and carnivore local extinctions in a fragmented forest landscape of southern Amazonia. Biol. Conserv. 124(3): 383-396. doi: https://doi.org/10.1016/j.biocon.2005.01.045Links ]

Norvell RE, Howe FP, Parrish JR. 2003. A seven-year comparison of relative-abundance and distance-sampling methods. Auk. 120(4): 1013 -1028. doi: https://doi.org/10.1093/auk/120.4.1013Links ]

Palma AC, Vélez A, Gómez-Posada C, López H, Zárate DA, Stevenson PR. 2011. Use of space, activity patterns, and foraging behavior of red howler monkeys (Alouatta seniculus) in an Andean forest fragment in Colombia. Am. J. Primatol. 73(10): 1062-1071. doi: https://doi.org/10.1002/ajp.20973Links ]

Peres CA. 2001. Synergistic effects of subsistence hunting and habitat fragmentation on Amazonian forest vertebrates. Conserv. Biol. 15(6): 1490-1505. doi: https://doi.org/10.1046/j.1523-1739.2001.01089.xLinks ]

Peres CA, Dolman PM. 2000. Density compensation in Neotropical primate communities: evidence from 56 hunted and nonhunted Amazonian forests of varying productivity. Oecologia. 122: 175-189. doi: https://doi.org/10.1007/PL00008845Links ]

Pfeiffer DU, Robinson TP, Stevenson M, Stevens KB, Rogers DJ, Clements AC. 2008. Spatial Analysis in Epidemiology. New York: Oxford University Press. [ Links ]

Pyritz LW, Büntge ABS, Herzog SK, Kessler M. 2010. Effects of habitat structure and fragmentation on diversity and abundance of primates in tropical deciduous forests in Bolivia. Int. J. Primatol. 31: 796-812. doi: https://doi.org/10.1007/s10764-010-9429-zLinks ]

Reckhow KH. 1990. Bayesian inference in non-replicated ecological studies. Ecology 71(6): 2053-2059. doi: https://doi.org/10.2307/1938619Links ]

Rimbach R, Link A, Heistermann M, Gómez-Posada C, Galvis N, Heymann EW. 2013. Effects of logging, hunting, and forest fragment size on physiological stress levels of two sympatric ateline primates in Colombia. Conserv. Physiol. 1(1): cot031. doi: https://doi.org/10.1093/conphys/cot031Links ]

Rode KD, Chapman CA, McDowell LR, Stickler C. 2006. Nutritional correlates of population density across habitats and logging intensities in redtail monkeys (Cercopithecus ascanius). Biotropica. 38(5): 625-634. doi: https://doi.org/10.1111/j.1744-7429.2006.00183.xLinks ]

Rodríguez-Luna, Cortés-Ortiz EL, Miller PH, Ellis S. 1996. Population and Habitat Viability Assessment for the Mantled Howler Monkey (Alouatta palliata mexicana). Apple Valley: Minnesota: IUCN/SSC Conservation Breeding Specialist Group. [ Links ]

Roncancio Duque, N, Gómez-Posada C. 2009. Sympatric Alouatta seniculus and Cebus capucinus in an Andean forest fragment in Colombia: a survey of population density. Neotrop. Primates 16(29): 51-56. doi: https://doi.org/10.1896/044.016.0201Links ]

Roncancio N, Rojas W, Defler T. 2011. Densidad poblacional de Saguinus leucopus en remanentes de bosque con diferentes características físicas y biológicas. Mastozool. Neotrop. 18: 105-117. [ Links ]

Roncancio-Duque NJ, Rojas-Díaz V, Ríos-Franco C, Gómez-Posada C, Gutiérrez-Chacón C, Giraldo P, Velasco J, Franco P. 2012. Plan de conservación y manejo del tití gris (Saguinus leucopus). Pereira, Colombia: Carder - Corpocaldas [ Links ]

Rojas-V W. 2009. Densidad Poblacional y Tamaño de Grupo de Titi Gris (Saguinus leucopus Günther 1876) en los Municipios de Mariquita Fresno y Falan (Tolima). [Unpublished report]. [Bogotá, Colombia]: Fundación Biodiversa. [ Links ]

Rylands AB. 1996. Habitat and the evolution of social and reproductive behavior in Callitrichidae. Am. J. Primatol. 38(1): 5-18. doi: https://doi.org/10.1002/(SICI)1098-2345(1996)38:1<5::AID-AJP2>3.0.CO;2-2Links ]

Santamaría M, Devenís C, Franco AM, Sánchez JD, Mendoza H, Pedraza C, Forero F, Medina S. 2007. Densidad y Ecología Básica del Tití Gris Saguinus leucopus en el Área de Influencia del Embalse Amaní, Norcasia, Caldas. [Unpublished report]. [Bogotá, Colombia]: Instituto Alexander Von Humboldt - ISAGEN. [ Links ]

Sanderson EW, Redford KH, Vedder A, Coppolillo PB, Ward SE. 2002. A conceptual model for conservation planning based on landscape species requirements. Landsc. Urban Plan. 58(1): 41-56. doi: https://doi.org/10.1016/S0169-2046(01)00231-6Links ]

Silva FE, Bicca-Marques JC. 2013. Do patch size and interpatch distance influence the distribution of brown howler monkeys (Alouatta guariba clamitans) in a fragmented landscape in south Brazil? In: Marsh LK, Chapman CA, editors. Primates in Fragments. New York: Springer . p. 137-145. [ Links ]

Snowdon CT, Soini P. 1988. The tamarins, genus Saguinus. In: Mittermeier RA, Rylands AB, Coimbra-F AF, da Fonseca GAB, editors. Ecology and Behaviour of Neotropical Primates. Vol. 2. Washington, DC: World Wildlife Fund. p. 223-298. [ Links ]

Sokal RR, Rohlf FJ. 1995. Biometry: The Principals and Practice of Statistics in Biological Research. New York: W. H. Freeman Co. [ Links ]

Sorensen TC, Fedigan LM. 2000. Distribution of three monkey species along a gradient of regenerating tropical dry forest. Biol. Conserv. 92(2): 227-240. doi: https://doi.org/10.1016/S0006-3207(99)00068-3Links ]

Spehar SN, Link A, Di Fiore A. 2010. Male and female range use in a group of white-bellied spider monkeys (Ateles belzebuth) in Yasunı National Park, Ecuador. Am. J. Primatol. 72(2): 129-141. doi: https://doi.org/10.1002/ajp.20763Links ]

Thomas L, Buckland ST, Rexstad EA, Laake JL, Strindberg S, Hedley SL, Bishop JRB, Marques TA, Burnham KP. 2010. Distance software: design and analysis of distance sampling surveys for estimating population size. J. Appl. Ecol. 47(1): 5-14. doi: https://doi.org/10.1111/j.1365-2664.2009.01737.xLinks ]

USGS. 2000. Shuttle radar topography mission, 3 Arc Second, Unfilled Unfinished 2.0, Global Land Cover Facility. College Park, Maryland: University of Maryland. [ Links ]

Wade PR. 2001. Bayesian methods in conservation biology. Conserv. Biol. 14(5): 1308-1316. doi: https://doi.org/10.1046/j.1523-1739.2000.99415.xLinks ]

Williams-Guillén K, Hagell S, Otterstrom S, Spehar S. Gómez C. 2013. Primate populations in fragmented tropical dry forest landscapes in southwestern Nicaragua. In: Marsh L, Chapman C, editors. Primates in Fragments. New York: Springer . p. 105-120. [ Links ]

Youlatos D. 2004. Multivariate analysis of organismal and habitat parameters in two neotropical primate communities. Am. J. Phys. Anthropol. 123(2): 181-194. doi: https://doi.org/10.1002/ajpa.10307 Links ]

Citation: Roncancio-Duque N. 2021. Effect of landscape modification on primate assemblages of the Magdalena River Valley, Colombia. Caldasia 43(2):261-273. doi: https://doi.org/10.15446/caldasia.v43n2.84845

SUPPLEMENTARY MATERIAL

Table 1 Model summary of the GLMMs for the effect of the gradient sampling station (matrix, edge, 75m, 150m) on bat richness and bat abundance. GLMM was fitted using a Poisson error structure and a log link function with fragment as a random factor. Significant p-values (<0.05) are given in bold. SE= Standard Error. 

Table 2 Post-hoc multiple comparisons of means using Tukey’s all-pairwise comparisons to determine differences in bat richness and bat abundance across gradient sampling stations. Significant post-hoc comparisons are given in bold (p-values <0.05). SE= Standard Error. 

Table 3 Post-hoc multiple comparisons of means using Wilcoxon pairwise comparisons to determine differences in the inverse of Simpson index across gradient sampling stations and between fragments. Significant post-hoc comparisons are given in bold (p-values <0.05). 

Table 4 Results of multivariate generalized linear model analysis testing for the effect of the gradient sampling station (matrix, edge, 75m, 150m) and fragment (Cachalú, Pericos, San Benito, and Tapias) on assemblage composition and species abundance for each trophic guilds. Significant p-values for predictor variables are given in bold and only species with p-values <0.05 in adjusted post-hoc univariate tests are shown. 

Figure 1 Species accumulation curves for captured individuals in a. matrix, b. edge, c. 75m, d. 150m, and e. total sampling, in four Sub-Andean forest fragments in Encino (Santander, Colombia) 

Received: February 05, 2020; Accepted: May 04, 2021; Published: May 05, 2021

* Author for correspondence.

CONFLICT OF INTEREST

The authors declares that they do not have conflict of interest.

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