Effect of free-ranging cattle on mammalian diversity: an Austral Yungas case study
←
→
Page content transcription
If your browser does not render page correctly, please read the page content below
Effect of free-ranging cattle on mammalian diversity: an Austral Yungas case study GRIET AN ERICA CUYCKENS, NOELIA VIVIANA GONZALEZ BAFFA TRASCI P A B L O G A S T Ó N P E R O V I C and L U C I O R I C A R D O M A L I Z I A Abstract Extensive cattle ranging is an important economic Introduction activity in mountains, with diverse effects on native mam- mal communities. The effects of cattle Bos taurus can be negative, positive or neutral, mostly depending on the stock- F ree-ranging cattle affect native fauna, species interac- tions and ecological communities across geographical re- gions in a variety of ways (Elliott & Barrett, ; Moser & ing rate. We examined the effect of cattle on the diversity and abundance of native mammalian species in the Witmer, ; Hettinger, ; Pia et al., ; Shepherd & Austral Yungas region of Argentina, considering environ- Ditgen, ; Chaikina & Ruckstuhl, ). In forest mental variables, land protection status, and human influ- ecosystems in Australia, the presence of free-ranging cattle ence. Using , trap-nights from camera-trap negatively affects vegetation regeneration by compacting stations over years (–), we calculated a relative the soil through trampling and diminishing the abundance abundance index using camera events and used generalized of seedlings by grazing (Eldridge et al., ). In South linear models to estimate the effect of cattle on small mam- American forests, the effect of cattle Bos taurus is hetero- mals, large herbivores, species of conservation concern and geneous and studies have been restricted to a few forest eco- felids. Cattle had different effects on each group of native systems and ecological variables (Mazzini et al., ). mammals. We observed a lower abundance of large native The intensive use of the forest understorey by cattle can herbivores and the absence of small mammals in areas cause a reduction in plant biomass, potentially reducing the with high cattle abundance. The tapir Tapirus terrestris, jag- complexity of the understorey (Loeser et al., ) and in- uar Panthera onca and white-lipped peccary Tayassu pecari creasing the density of shrubs resistant to browsing are rare in the Yungas and therefore potentially vulnerable (Vandenberghe et al., ), thus changing forest structure to extinction there. Conservation of small felids and low cat- and composition. These changes in vegetation could have tle abundance could be compatible, but felids are threatened both direct and indirect effects on native mammal biodiver- by other anthropogenic influences. Native mammalian di- sity. By affecting the understorey of forests and reducing ref- versity and richness were related to land protection status. uge and food availability, cattle have the potential to affect The entire ecoregion is potentially suitable for cattle, sug- small forest mammals negatively (Tabeni et al., ). Cattle gesting the potential for further threats, and that cattle may also affect large native herbivores by reducing food should be excluded from strictly protected areas. To ensure diversity and competing for pastures (Madhusudan, extensive cattle ranging is compatible with wildlife conser- ). By altering small mammal abundance (i.e. prey for vation in areas where exclusion is not possible, we recom- carnivores), cattle also cause cascading effects on higher mend improved management of cattle and moderate trophic levels (Pia et al., ). stocking rates. The influence of cattle is not straightforward (Hettinger, ) and is mainly determined by the stocking rate Keywords Bos taurus, cattle, empty forest, human foot- (Schieltz & Rubenstein, ). Mammalian species exhibit print, livestock, national parks, niche-based distribution different responses linked to their specific life traits modelling, threatened species (Suraci et al., ). Forest specialist mammal species decline Supplementary material for this article is available at rapidly when forest cover decreases and are unlikely to be doi.org/./S found in secondary forests, whereas non-specialists survive in human-modified habitats (Newbold et al., ), and some species such as the white-eared opossum Didelphis al- biventris and Molina’s hog-nosed skunk Conepatus chinga GRIET AN ERICA CUYCKENS (Corresponding author, orcid.org/0000-0002-9019- could benefit from cattle-modified habitat (Di Bitetti et al., 8916, grietcuyckens@yahoo.com), NOELIA VIVIANA GONZALEZ BAFFA TRASCI ( orcid.org/0000-0001-5666-8614) and LUCIO RICARDO MALIZIA Instituto de ). Ecorregiones Andinas–Universidad Nacional de Jujuy, Consejo Nacional de To ensure the protection of biodiversity, a total exclusion Investigaciones Científicas y Tecnológicas and Centro de Estudios Territoriales Ambientales y Sociales, Alberdi 47, 4600 San Salvador de Jujuy, of cattle and strict controls against hunting are imposed Argentina by national parks in Argentina (IUCN category II pro- PABLO GASTÓN PEROVIC Administración de Parques Nacionales, Salta, Argentina tected areas). Other categories (VI) allow traditional cattle Received December . Revision requested July . raising. In regions where cattle are managed extensively Accepted October . and with little veterinary control, the management of a This is an Open Access article, distributed under the terms of the Creative Commons-Attribution-ShareAlike licence (https://creativecommons.org/licenses/by-sa/4.0/), which permits re-use, distribution, reproduction, transformation, and adaptation in any medium and for any purpose, provided the original work is properly cited and any transformation/adaptation is distributed under the same Creative Commons licence. Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
878 G. A. E. Cuyckens et al. site plays a fundamental role in the local stocking rate, and Species groups therefore in fauna conservation (Schieltz & Rubenstein, ). Of the recorded species, we selected four groups that can Since the s, the range of domestic cattle in Argentina potentially be influenced by cattle. By reducing refuge and has expanded into areas marginal for agriculture, mostly food availability for small mammals (Tabeni et al., ), forested areas (Guevara et al., ) such as the Austral we expected a negative influence of cattle on the presence Yungas in the north-west. In the Argentine Yungas, cattle and relative abundance of small mammal species. In raising entails releasing cattle into the forests, without nu- this group we included unidentified small mammals tritional supplements and with little to no management (i.e. # kg), the agouti Dasyprocta sp. and the tapeti (Quiroga et al., ). This extensive ranging of livestock, Sylvilagus brasiliensis. Cattle may also affect large native combined with selective logging and firewood extraction, herbivores by reducing food diversity and competing for has adverse effects on the structure of native forests pastures (Madhusudan, ). In the large herbivore species (Campanello et al., ; Blundo et al., ). However, group we included the red brocket deer Mazama americana the effect of cattle has not been considered as an explanatory and gray brocket deer Mazama gouazoubira. We expected a variable when examining patterns of mammalian diversity, lower relative abundance of both species with higher cattle nor is it known how cattle affect fauna in one of the most abundance, and absence of the species at a certain, un- biodiverse ecosystems in Argentina. known, threshold of cattle abundance. We predict similar We aimed to assess the effect of cattle on the native mam- effects on these two species, and therefore we pooled data malian community of the Austral Yungas and any potential for a more robust statistical analysis. For species of conser- interaction with altitude, latitude and land protection status. vation concern we included threatened species based on We considered four species groups: small mammals, large national and/or international standards (Ministerio de herbivores, species of conservation concern and the felid Ambiente y Desarrollo Sostenible & Sociedad Argentina community. We also examined the potential of cattle to in- para el Estudio de los Mamíferos, ; IUCN, ): the habit the Yungas ecoregion on a regional scale. We expected lowland tapir Tapirus terrestris, white-lipped peccary Tayassu () lower mammalian richness and diversity at higher alti- pecari and jaguar Panthera onca. Based on the potential tudes, increasing latitude and in areas of lower protection influence of cattle on the relative abundance of small mam- status, () lower abundance or absence of the four species mals and cascading effects on higher trophic levels (Pia groups in areas with greater cattle abundance and anthropo- et al., ), we expect a negative influence of cattle on the genic influence, and () that cattle could potentially inhabit presence and relative abundance of small and medium- the entire Yungas ecoregion. sized felids. This group comprised the medium-sized ocelot Leopardus pardalis, and five small felids: jaguarundi Herpailurus yagouaroundi, Pampas cat Leopardus colocolo, Geoffroy’s cat Leopardus geoffroyi, oncilla Leopardus Study area tigrinus and margay Leopardus wiedii (Table ). The study area is the Austral Yungas of Argentina (sensu Brown & Pacheco, ) on the eastern slope of the Andes. Camera-trap survey This ecoregion is characterized by subtropical cloudy montane forests and has an altitudinal gradient of vege- We placed camera stations: in national parks, in tation physiognomy and species composition (Brown et al., provincial reserves, in private protected areas, in ). The Yungas is considered a vulnerable ecosystem private lands without protection, in state properties (Olson & Dinerstein, ) of high conservation value without management, and nine in Indigenous territories. (Malizia et al., ) because of the high faunal diversity The trapping period was January – September (Narosky & Yzurieta, ; Ojeda, ). Cattle were intro- (, effective trap-nights). We recorded geographical duced into the region c. years ago (Brown & Grau, ), coordinates and altitude at each point using a GPS. but their distribution is limited by the terrain and therefore Camera-trap stations consisted of one camera, generally cattle density is highly variable. located along a trail, road or river bank, to optimize the capture of larger mammals, which prefer to walk along lin- ear features (Harmsen et al., ). In some cases, to ensure Methods cameras were not interfered with or stolen, we placed them in the forest interior. The mean distance between We surveyed using camera traps across latitudinal and nearest neighbouring cameras was . ± SD . km. Cam- altitudinal gradients in areas with forest cover (Fig. ). eras were programmed to obtain a set of three photographs, To measure the success of our method, we developed a po- with a -minute delay between successive sets, and were tential species list for the Yungas forests (Table ). active continuously. Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
Effect of free-ranging cattle on mammalian diversity 879 FIG. 1 Study area and location of camera traps in the Austral Yungas of Argentina. Size of circles indicates the number of native mammalian species recorded per camera trap. Environmental variables brooks, rills and runnels). All variables were projected to the WGS datum and were at a spatial resolution of We used bioclimatic variables (Karger et al., ). From arc-seconds or resampled at this pixel size, equivalent altitude, we derived slope and roughness, using QGIS .. to c. km. (QGIS, ), which we also used for all other geographical analyses. We obtained the human influence index from WCS & CIESIN (); this index incorporates human Mammalian diversity population pressure, human land use, infrastructure and human access. Values range from (no influence) to We identified all large and medium-sized animals recorded (maximum influence). We obtained the coupled evapo- by the camera traps to species level, noting whether they transpiration and gross primary production (hereafter re- were native or exotic, and calculated diversity indices for ferred to as primary production) from the National the native species. Records were considered independent Tibetan Plateau Data Center (Zhang et al., , ; when they were at least h apart. From the number of events Gan et al., ), and we generated Euclidean distances and effort (i.e. the number of effective camera-trap nights from cameras to the nearest water line (rivers, streams, of each camera), we calculated the relative abundance Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
880 G. A. E. Cuyckens et al. TABLE 1 Species of large and medium-sized native and exotic mammals that could potentially occur in the Austral Yungas (Ministerio de Ambiente y Desarrollo Sostenible & Sociedad Argentina para el Estudio de los Mamíferos, ; Fig. ), with the species group for those species included in our analysis (i.e. small mammals, large herbivores, species of conservation concern and felids), national and IUCN Red List status for native species, and whether or not recorded in our camera-trap study during –. The table does not include uniden- tified small mammal species. National Native/ Species (by family) Species group status1 Red List status2 exotic Recorded Didelphidae White-eared oposum Didelphis albiventris LC LC Native Yes Thick-tailed oposum Lutreolina massoia NT LC Native Yes Mymercophagidae Southern tamandua Tamandua tetradactyla NT LC Native Yes Giant anteater Myrmecophaga tridactyla VU VU Native Yes Dasypodidae Armadillo Dasypus mazzai DD DD (as Dasypus Native No yepesi) Nine-banded long-nosed armadillo Dasypus LC LC Native No novemcinctus Six-banded armadillo Euphractus sexcinctus LC LC Native Yes Cebidae Brown capuchin monkey Sapajus cay VU LC Native Yes Canidae Crab-eating fox Cerdocyon thous LC LC Native Yes Pampas fox Lycalopex gymnocercus LC LC Native Yes Domestic dog Canis lupus Exotic Yes Felidae Jaguarundi Herpailurus yagouaroundi Felid LC LC Native Yes Pampas cat Leopardus colocolo Felid VU NT Native Yes Geoffroy’s cat Leopardus geoffroyi Felid LC LC Native Yes Oncilla Leopardus tigrinus Felid VU VU Native Yes Ocelot Leopardus pardalis Felid VU LC Native Yes Margay Leopardus wiedii Felid VU NT Native Yes Cougar Puma concolor LC LC Native Yes Jaguar Panthera onca Conservation CR NT Native Yes concern Domestic cat Felis catus Exotic No Mustelidae Molina’s hog-nosed skunk Conepatus chinga LC LC Native Yes Neotropical river otter Lontra longicaudis NT NT Native No Tayra Eira barbara NT LC Native Yes Little grison Galictis cuja LC LC Native Yes Procyonidae Crab-eating raccoon Procyon cancrivorus LC LC Native Yes South American coati Nasua nasua LC LC Native Yes Tapiridae Lowland tapir Tapirus terrestris Conservation VU VU Native Yes concern Tayassuidae Collared peccary Pecari tajacu VU LC Native Yes White-lipped peccary Tayassu pecari Conservation EN VU Native Yes concern Cervidae Red brocket deer Mazama americana Herbivore VU DD Native Yes Gray brocket deer Mazama gouazoubira Herbivore LC LC Native Yes Leporidae Tapeti Sylvilagus brasiliensis Small mammal LC EN Native Yes Common hare Lepus europaeus Exotic No Sciuridae Yungas squirrel Sciurus ignitus Native Yes Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
Effect of free-ranging cattle on mammalian diversity 881 TABLE 1 (Cont.) Species (by family) Species group National Red List status2 Native/ Recorded status1 exotic Erithizontidae Bicolored-spined porcupine Coendou bicolor VU LC Native No Prehensile-tailed porcupine Coendou prehensilis VU LC Native No Hydrochaeridae Capybara Hydrochoerus hydrochaeris LC LC Native No Dasyproctidae Yungas agouti Dasyprocta sp. Small mammal LC LC Native Yes Myocastoridae Coypu Myocastor coypus LC LC Native No Suidae Pig Sus scrofa Exotic Yes Bovidae Cattle Bos taurus Exotic Yes Equidae Horse Equus caballus Exotic Yes Donkey Equus asinus Exotic Yes Ovidae Sheep Ovis aries Exotic Yes Capridae Goat Capra sp. Exotic Yes Ministerio de Ambiente y Desarrollo Sostenible & Sociedad Argentina para el Estudio de los Mamíferos (). IUCN (). DD, Deficient Data; LC, Least Concern; NT, Near Threatened; EN, Endangered; CR Critically Endangered. index (RAI) as: abundance, primary production and the human influence index; for species of conservation concern: cattle abun- RAIi = ntot/daystot × 100 (1) dance, primary production, human influence index, land where ntot is the number of independent events of the ith protection status and distance to water lines; for the felid species and daystot is the total number of effective trap- community: cattle relative abundance index and the nights, using the package camtrapR (Bengsen et al., ; human impact index. For presence/absence models, we Mandujano & Pérez-Solano, ) in R .. (R Core used the negative binomial error and Gaussian distributions Team, ). We calculated two measures of diversity, for abundance, using the package MASS (Venables & Ripley, using the package vegan (Oksanen et al., ) in R: species ) in R. We checked for homogeneity by plotting resi- richness (S) and the Shannon–Weaver index (H). The latter duals vs fitted values, for normality using quantile-quantile was calculated as: plots, and for independence by plotting residuals vs each ex- planatory variable. Because we had several combinations of S H=− RAI i × ln(RAI i ) (2) variables and therefore multiple models, we used single- i=1 term deletions to obtain the most parsimonious model, The value of the relative abundance index increases with in- using Akaike’s information criterion (AIC; Burnham & creasing richness and evenness in the abundance of species Anderson, ). in the community. We developed generalized linear models (GLM) to examine the effect of cattle and land protection status on A niche-based model for cattle S and H. As biodiversity follows global patterns, with an in- crease in species richness toward the tropics and a decline in Species distribution models examine the potential influence species richness with increasing elevation (Pianka, ; of environmental variables on species presence. MaxEnt Lomolino, ; Hillebrand, ), we included altitude finds the distribution of maximum entropy (i.e. the largest and latitude as factors in the models. For each of the four spread in a geographical dataset of species presences), sub- species groups we examined the potential influence of ex- ject to the constraint that the projected value of each variable planatory variables on presence/absence of species and on is close to its empirical average (Phillips et al., ). This relative abundance index; for small mammals: cattle abun- information can then be extrapolated to non-sampled dance and primary production; for large herbivores: cattle areas (Phillips et al., ). Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
882 G. A. E. Cuyckens et al. FIG. 2 Mean native species richness recorded by camera traps set in areas of varying protection status. Boxes represent the interquartile range, box widths are proportional to the square root of the sample sizes, the horizontal solid line indicates the FIG. 3 Habitat suitability for cattle Bos taurus (from the median, and the whiskers indicate the % CI. niche-based model) in relation to the human influence index. The study area for modelling was the Yungas ecoregion recorded by (%) cameras. Gray brocket deer had a with a -km buffer. We generated , random points relative abundance index of –, and red brocket deer to extract values for the environmental variables, and tested –. Among the species of conservation concern, the for correlations with the Pearson test, selecting only those lowland tapir had the highest relative abundance index (re- with a correlation # .. These were: mean diurnal tem- corded by cameras), the white-lipped peccary had the perature range, temperature seasonality, mean temperature lowest relative abundance index (three cameras) and the jag- of wettest quarter, precipitation of driest quarter, precipita- uar was recorded by cameras. The ocelot was the most tion of warmest quarter, roughness, human influence index, abundant felid species, followed by the margay, jaguarundi, and distance to nearest water lines. oncilla, Pampas cat and Geoffroy’s cat. The maximum num- We used presence records of cattle obtained from the ber of small and medium-sized felid species detected by one camera traps. We ran the species distribution model camera was five. times and used the average outcome. To measure the general Native species richness was almost % higher in nation- performance, we used the area under the receiver operating al parks and private protected areas than in provincial characteristic curve (AUC), which measures the probability protected areas, with intermediate values for Indigenous that a randomly chosen presence point will rank above a territories and private lands (Fig. ). Native species richness randomly chosen background point (AUC . = random; and diversity decreased with altitude, without interactions values closer to indicate better discrimination power; with cattle relative abundance index or latitude (Fig. , Bellamy et al., ). We converted the model into a binary Table ). model by applying the minimum presence logistic threshold Presence, but not relative abundance, of small mammals (.) and categorized habitat suitability for cattle into five and large native herbivores was influenced by cattle relative classes (very low: , .; low: .–.; medium: .– abundance index. For small mammals to be present, the .; high: .–.; very high: ˃ .) for visualization as relative abundance index of cattle had to be , , and large a map. herbivores were absent at a cattle relative abundance index of and present at . Thus, to assure the presence of large herbivores, the number of independent cattle records Results should not be greater than five times the number of effective trap-nights (for example, in a survey of trap-nights the We recorded native species, % of the species that could number of independent cattle events should be , ). potentially occur in the region (Table ), and seven exotic Lowland tapir relative abundance increased with distance species, with of cameras (%) recording at least from water lines and was not influenced by cattle relative one native species and cameras (%) at least one exotic abundance index. Jaguar relative abundance and presence species. The range of species richness per camera was – were not significantly associated with any of the measured (mean .) for native and – (mean .) for exotic species. variables. Data for the white-lipped peccary could not be Bos taurus had the highest relative abundance index and was analysed because there were only three records. The also the most frequently recorded species. Other domestic human influence index negatively affected felid richness, species recorded were horses, pigs, goats and dogs. but not relative abundance, with an abrupt decrease in rich- Domestic cats were not recorded. ness at a human influence index of $ . Primary produc- Small mammals were recorded by cameras (%). tion did not influence the presence or relative abundance of Brocket deer (i.e. M. gouazoubira or M. americana) were small mammals or the tapir (Table ). Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
Effect of free-ranging cattle on mammalian diversity 883 TABLE 2 Generalized linear models for native species richness absent where the cattle relative abundance index was (number of species) and diversity (Shannon−Weaver index), above and , respectively. Felids were negatively affected small mammal presence/absence, deer Mazama sp. presence/ab- by the human influence index but not by cattle relative sence, lowland tapir Tapirus terrestris relative abundance index abundance, suggesting that felids are influenced not by cat- and felid species richness in the Austral Yungas. Only significant models are presented, with t and P values for potentially influential tle directly but by activities related to their presence. As cat- variables. tle have the potential to inhabit most of the Austral Yungas and their presence is associated with hunting, exotic species t P (domestic dogs and cats), and selective logging (Perovic, Species richness ), we recommend cattle should be reduced to abun- Latitude 0.036 0.846 dances that allow coexistence with wildlife in all areas Altitude −0.002 0.0002** with forest cover, and excluded from strictly protected Protection status −0.850 2.75 × 10−08** Cattle relative abundance index −0.005 0.415 areas. Strictly protected areas (private or state-managed) Species diversity are the only management regime ensuring long-term Latitude 0.060 0.197 fauna conservation in the Yungas. Altitude −4.96 × 10−04 1.54 × 10−06** We recorded % of the native species that could poten- Land protection status −0.123 0.001* tially occur in the Austral Yungas and therefore we consider Cattle relative abundance index −0.002 0.239 our methodology successful. We did not record water- Small mammals associated species such as the capybara Hydrochoerus hy- Cattle relative abundance index −2.100 0.041* drochaeris, nutria Myocastor coypus and neotropical otter Primary productivity −0.078 0.938 Deer Lontra longicaudis, the latter being rare in the Yungas Cattle relative abundance index 0.023 0.021* (Albanesi et al., ). As we did not place cameras in trees, Human influence index 1.859 0.067 our methodology was not suitable for detecting arboreal Tapir relative abundance index species such as the bicolored-spined porcupine Coendou Cattle relative abundance index −0.039 0.781 bicolor and prehensile-tailed porcupine Coendou prehensi- Primary productivity 0.937 0.352 lis, which are also rare. We did, however, record arboreal Distance to rivers 2.070 0.042* species such as the capuchin monkey Sapajus cay and Felid species richness Cattle relative abundance index −0.814 0.418 Bolivian squirrel Sciurus ignitus, but on the ground. Only Human influence index −2.870 0.005* two of the six cingulate species (armadillos) were recorded, suggesting they have a naturally low relative abundance or *P , .; **P , .. are difficult to record with camera traps. A national-scale as- sessment indicated the need to protect these six species We obtained a habitat suitability model for cattle with (Abba et al., ). AUC = .. Values . . are considered to indicate mod- The niche-based distribution model for cattle shows their els with a good general performance (Phillips & Dudík, potential to occupy almost the entire Yungas ecoregion ex- ). Cattle encounter suitable habitat in most (%) of cept for the mountain peaks, probably because of the low the Yungas, with higher habitat suitability towards lower la- winter temperatures, low carrying capacity and the difficulty titudes and lower suitability towards higher altitudes of access for people. The positive association between habi- (Fig. ). Habitat suitability for cattle increases with the tat suitability for cattle and the human influence index is a human influence index, peaking at and then decreasing result of the association of cattle with people, but areas with abruptly (Fig. ). Lower suitability corresponds to areas al- a human influence index ˃ are no longer suitable for cat- ready transformed into croplands, human settlements or tle. The raster layer of suitable habitat for cattle could be roads (Fig. ). used in areas not surveyed to estimate cattle presence locally, and could serve as a tool to analyse the effects of cattle in the Austral Yungas. Discussion Environmental variables, which are influenced by lati- tude and altitude, affect the diversity and composition of As far as we are aware, this is the first study to use a large native biodiversity. As predicted, we found a decrease in camera-trap data set for the Yungas ecoregion over an ex- species richness and diversity with an increase in elevation, tended period ( years) and the first to examine the factors in accordance with global patterns (Lomolino, ) and affecting the native mammalian community. In contrast to with a previous study in the northern Yungas (Di Bitetti our expectation, general mammalian species richness and et al., ). We found high mountain areas in the Yungas diversity were not directly related to the cattle relative abun- to be naturally poorer in native and exotic species than for- dance index, but the latter affected two of the four groups of ests at lower elevations. The northern Austral Yungas has species studied. Small mammals and brocket deer were higher species richness than the central and southern Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
884 G. A. E. Cuyckens et al. FIG. 4 Presence records of cattle with (a) predicted habitat suitability (see text for details), and (b) the human influence index (see text for details) in the Austral Yungas of Argentina. ), depending on their management. National Parks in Argentina are legally required to exclude cattle, although this regulation is not always implemented. They are gener- ally larger and have stricter controls than provincial and pri- vate protected areas and have historically been established in areas with low potential for economic development (Margules & Pressey, ; Rodrigues et al., ), and hence could have a higher intrinsic value for conservation. Thus, various anthropogenic factors seem to be implicated in the persistence of native large mammals in these areas. FIG. 5 Native species richness per camera trap in relation to Indigenous territories, with intermediate cattle abundance, elevation during an -year camera-trap study in the Austral may contribute to conservation and offer a complementary Yungas of Argentina. institutional model to state-run protected areas (Johnson & Nelson, ). Provincial protected areas, with higher cattle Austral Yungas and is the southern limit of the distributions relative abundance index, had the lowest diversity indices, of several mammal species (Sapajus cay, Tapirus terrestris, similar to unprotected areas (private or state property). In Leopardus wiedii), emphasizing the role of this part of the low-income countries, nature conservation is not neces- ecoregion in the conservation of these species. In contrast sarily a priority and so-called paper parks (i.e. protected to general patterns (Brown & Lomolino, ) and to areas that only exist on paper and do not achieve conserva- those of small forest mammals in the Yungas (Ojeda et al., tion goals), are common (Rodríguez & Rodríguez-Clark, ), we found that latitude did not influence native spe- ). Nevertheless, even paper parks matter (Rodríguez & cies richness and diversity, indicating a latitudinally homo- Rodríguez-Clark, ) if they still have forest cover. To re- genized native mammal community. verse the failure in the achievement of conservation objec- Land protection status was the most important variable tives of such paper parks in Argentina and to achieve in explaining native mammalian biodiversity. We recorded conservation goals, we recommend lowering cattle relative the highest values of native species richness and diversity in abundance index to # (the highest limit compatible with national parks, highlighting the importance of strictly pro- the presence of both deer and small mammals) and estab- tected areas and the complementary role of small private lishing adaptive management plans that include stricter protected lands (Johnson & Nelson, ; Kamal et al., controls than at present. Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
Effect of free-ranging cattle on mammalian diversity 885 Our findings indicate that high cattle abundance (relative cattle and human activities such as habitat destruction and abundance ˃ ) is incompatible with the presence of small hunting for food (peccaries) and in retaliation (jaguars), are mammals. Trampling and browsing by cattle reduce the het- probably having negative, synergistic effects on species of erogeneity of the forest understorey for this group (Smith conservation concern (Romero-Muñoz et al., ). These et al., ; Hayward et al., ). We could not identify species depend on protected forests with extensive cover small mammals to species level and the responses of indi- and protection against hunting. vidual species to cattle may vary (Schieltz & Rubenstein, This is the first study based on an extensive camera-trap ). Nevertheless, when present, small mammals had high survey to provide evidence that cattle affect the assemblage relative abundance indices, so a prey base for medium-sized of mammals in the Austral Yungas, both directly and indir- and small felids remains available. ectly. We have provided guidelines for cattle abundance that Small felids are indirectly affected by cattle, as indicated should be implemented in protected areas where cattle rais- by the human influence index. This index includes human ing is allowed. However, we cannot provide guidelines for population pressure, human land use and infrastructure and cattle abundance compatible with species of conservation human access. Differential logging (i.e. logging of only par- concern (jaguar and white-lipped peccary), and their low ticular tree species and individuals of a certain size, resulting abundances indicate their high risk of extinction in this in impoverished species richness and affecting forest struc- region. Following Mazzini et al. (), we used directly ture) could also affect forest specialists such as the oncilla measured cattle relative abundance at the local scale and and margay. Therefore, we suggest direct hunting and asso- complemented this with an indirect method at the regional ciated activities, such as the presence of dogs (Perovic, ) scale (distribution modelling). Our work therefore provides and habitat transformation, are underlying explanations for both a method for future assessments of cattle impacts and the influence of humans on felids, not the presence of cattle an indicator of potential cattle abundance in unsurveyed directly. We did not record domestic cats, so these may not areas in the Austral Yungas. yet be a threat to small felids in the Yungas. We found that the Vulnerable oncilla had the lowest relative abundance Acknowledgements We thank all field assistants and landowners who allowed us to conduct fieldwork, the Tinkunaku community, and was the rarest of the felid group and therefore its con- park rangers of the Secretariat of Environment and Sustainable servation status should be monitored. Development of Salta province, the National Park Administration of In agreement with Nanni (), we found brocket deer Argentina; IDEA WILD for providing equipment; Leonidas only in areas with low cattle relative abundance. The human Lizarraga for help downloading data; Agustin Abba for assistance influence index did not affect these large herbivores, sug- with identification of ungulates; Tadeu de Oliveira for help with iden- tification of small felids; two anonymous reviewers for advice; and gesting a direct effect of cattle on these two species. In agree- Martin Fisher and the teachers and classmates of the Conservation ment with Mazzini et al. (), we suggest biological inter- Leadership Programme Writing for Conservation Course in Colombia. actions such as competition and dietary overlap between cattle and native herbivores are the cause of this effect. Author contributions Study conception: GAEC, LRM; study de- Species of conservation concern had low relative abun- sign, fieldwork: GAEC, PGP; data input: GAEC, NGBT; data analysis: GAEC, writing: all authors. dance indices. Lowland tapir relative abundance was nega- tively influenced by distance to linear watercourses. Like the Conflicts of interest None. Malayan tapir Tapirus indicus and Baird’s tapir Tapirus baardii, which depend on proximity to water (Dudgeon, Ethical standards This research abided by the Oryx guidelines on ; Reyna-Hurtado et al., ), lowland tapirs use ethical standards. water banks to browse and mate, and enter the water to take refuge from predators (Brooks et al., ). Lowland tapir relative abundance was not negatively influenced by References cattle relative abundance index. Tapirs probably differ in feeding habits from cattle, browsing more on seedlings A B B A , A.M., T O G N E L L I , M.F., S E I T Z , V.P., B E N D E R , J.B. & V I Z C A Í N O , and fruits. Prey remains available for jaguars, but the ab- S.F. () Distribution of extant xenarthrans (Mammalia: Xenarthra) in Argentina using species distribution models. undance of wild prey is decreasing, which could provoke Mammalia, , –. jaguars to start predating on cattle, with consequent A L B A N E S I , S.A., J AY AT , J.P., A L B E R T I , P. & B R O W N , A.D. () New escalations in human–predator conflict (Perovic, ; record of river otter (Lontra longicaudis Olfers, ) in the extreme Cuyckens et al., ). The species of conservation concern south of Yungas Northwestern Argentina. IUCN/Species Survival considered here (lowland tapir, white-lipped peccary and Commission Otter Specialist Group Bulletin, , –. B E L L A M Y , C., S C OT T , C. & A LT R I N G H A M , J. () Multiscale, jaguar) are large mammals, and body size is an indicator presence-only habitat suitability models: fine-resolution maps for of extinction risk (Cardillo, ) as it determines suscepti- eight bat species. Journal of Applied Ecology, , –. bility to hunting pressure and habitat selectivity. Our study B E N G S E N , A.J., L E U N G , L.K.-P., L A P I D G E , S.J. & G O R D O N , I.J. () was in areas that still have forest cover. Hence, the effects of Using a general index approach to analyze camera-trap abundance Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
886 G. A. E. Cuyckens et al. indices: improving camera-trap abundance indices. The Journal of production in the Monte Desert, Argentina. Journal of Arid Wildlife Management, , –. Environments, , –. B L U N D O , C., G A S P A R R I , N.I., M A L I Z I A , A., C L A R K , M., G AT T I , G., H A R M S E N , B.J., F O S T E R , R.J., S I LV E R , S., O S T R O , L. & D O N C A S T E R , C A M P A N E L LO , P.I. et al. () Relationships among phenology, C.P. () Differential use of trails by forest mammals and climate and biomass across subtropical forests in Argentina. the implications for camera-trap studies: a case study from Belize. Journal of Tropical Ecology, , –. Biotropica, , –. B R O O K S , D.M., B O D M E R , R.E. & M AT O L A , S. () Tapirs – Status, H A Y WA R D , B., H E S K E , E.J. & P A I N T E R , C.W. () Effects of livestock Survey and Conservation Action Plan. Tapir Specialist Group, grazing on small mammals at a desert Cienaga. The Journal of IUCN, Gland, Switzerland and Cambridge, UK. Wildlife Management, , –. B R O W N , A.D. & G R A U , H.R. () La Naturaleza y el Hombre H E T T I N G E R , N. () Defining and evaluating exotic species: issues en las Selvas de Montaña. Proyecto GTZ – Desarrollo for Yellowstone Park policy. Western North American Naturalist, Agroforestal en Comunidades Rurales del Noroeste Argentino, , –. Salta, Argentina. H I L L E B R A N D , H. () On the generality of the latitudinal diversity B R O W N , J.H. & L O M O L I N O , M.V. () Biogeography, nd edition. gradient. The American Naturalist, , –. Sinauer Asociates, Inc., Sunderland, USA. IUCN () The IUCN Red List of Threatened Species. Version -. B R O W N , A.D., G R A U , H.R., M A L I Z I A , L.R. & G R A U , A. () iucnredlist.org [accessed May ]. Argentina. In Bosques nublados del Neotrópico (eds M. Kapelle & J O H N S O N , K.A. & N E L S O N , K.C. () Common property and A.D. Brown), pp. –. INBio, San José, Costa Rica. conservation: the potential for effective communal forest B R O W N , A.D. & P AC H E C O , S. () Propuesta de actualización del management within a national park in Mexico. Human Ecology, mapa ecorregional de la Argentina. In La Situación Ambiental , –. Argentina (eds A. Brown, U. Martínez Ortíz, M. Acerbi & K A M A L , S., G R O D Z I Ń S K A -J U R C Z A K , M. & B R O W N , G. () J. Corcuera), p. . Fundación Vida Silvestre, Buenos Aires, Conservation on private land: a review of global strategies with a Argentina. proposed classification system. Journal of Environmental Planning B U R N H A M , K.P. & A N D E R S O N , D.R. () Model Selection and and Management, , –. Multimodel Interference: A practical Information-Theoretic K A R G E R , D.N., C O N R A D , O., B Ö H N E R , J., K AW O H L , T., K R E F T , H., Approach. nd edition. Springer, New York, USA. S O R I A -A U Z A , R.W. et al. () Climatologies at high resolution for C A M P A N E L LO , P.I., G E N O V E VA G AT T I , M., A R E S , A., M O N T T I , L. & the Earth’s land surface areas. Scientific Data, , . G O L D S T E I N , G. () Tree regeneration and microclimate in a L O E S E R , M.R., M E Z U L I S , S.D., S I S K , T.D. & T H E I M E R , T.C. () liana and bamboo-dominated semideciduous Atlantic Forest. Vegetation cover and forb responses to cattle exclusion: implications Forest Ecology and Management, , –. for pronghorn. Rangeland Ecology & Management, , –. C A R D I L LO , M. () Multiple causes of high extinction risk in large L O M O L I N O , M.V. () Elevation gradients of species-density: mammal species. Science, , –. historical and prospective views. Global Ecology and Biogeography, C H A I K I N A , N.A. & R U C K S T U H L , K.E. () The effect of cattle , –. grazing on native ungulates: The good, the bad, and the ugly. M A D H U S U D A N , M.D. () Recovery of wild large herbivores Rangelands, , –. following livestock decline in a tropical Indian wildlife reserve: C U Y C K E N S , G., F A L K E , F. & P E T R AC C A , L. () Jaguar Panthera onca livestock decline and wild herbivore recovery. Journal of Applied in its southernmost range: use of a corridor between Bolivia and Ecology, , –. Argentina. Endangered Species Research, , –. M A L I Z I A , L.R., P A C H E C O , S., B L U N D O , C. & B R O W N , A.D. () D I B I T E T T I , M.S., D I B L A N C O , Y.E. & J I M É N E Z P É R E Z , J.I. () Caracterización altitudinal, uso y conservación de las Yungas Cambios Asociados a la Presencia de Ganado Vacuno en los Subtropicales de Argentina. Ecosistemas, , –. Ensambles de Mamíferos de tres Ambientes de la Reserva Natural M A N D U J A N O , S. & P É R E Z -S O L A N O , L.A. (eds) () Fototrampeo en Iberá, Corrientes, Argentina. Asociación Argentina de Ecología, R: Organización y Análisis de Datos. Instituto de Ecología, Xalapa, Buenos Aires, Argentina. Mexico. D I B I T E T T I , M.S., A L B A N E S I , S.A., F O G U E T , M.J., D E A N G E LO , C. & M A R G U L E S , C.R. & P R E S S E Y , R.L. () Systematic conservation B R O W N , A.D. () The effect of anthropic pressures and elevation planning. Nature, , –. on the large and medium-sized terrestrial mammals of the M A Z Z I N I , F., R E LVA , M.A. & M A L I Z I A , L.R. () Impacts of subtropical mountain forests (Yungas) of NW Argentina. domestic cattle on forest and woody ecosystems in southern South Mammalian Biology, , –. America. Plant Ecology, , –. D U D G E O N , D. () Tropical Stream Ecology. Elsevier/Academic M I N I S T E R I O D E A M B I E N T E Y D E S A R R O L LO S O S T E N I B L E & S O C I E D A D Press, Amsterdam, The Netherlands. A R G E N T I N A P A R A E L E S T U D I O D E LO S M A M Í F E R O S () E L D R I D G E , D.J., P O O R E , A.G.B., R U I Z -C O L M E N E R O , M., L E T N I C , M. & Categorización de los Mamíferos de Argentina. cma.sarem.org.ar/es/ S O L I V E R E S , S. () Ecosystem structure, function, and especies-nativas [accessed September ]. composition in rangelands are negatively affected by livestock M O S E R , B.W. & W I T M E R , G.W. () The effects of elk and cattle grazing. Ecological Applications, , –. foraging on the vegetation, birds, and small mammals of the Bridge E L L I O T T , H.W. & B A R R E T T , R.H. () Dietary overlap among axis, Creek Wildlife Area, Oregon. International Biodeterioration, , –. fallow, and black-tailed deer and cattle. Journal of Range N A N N I , A.S. () Dissimilar responses of the gray brocket deer Management, , . (Mazama gouazoubira), crab-eating fox (Cerdocyon thous) and G A N , R., Z H A N G , Y., S H I , H., Y A N G , Y., E A M U S , D., C H E N G , L. et al. pampas fox (Lycalopex gymnocercus) to livestock frequency in () Use of satellite leaf area index estimating evapotranspiration subtropical forests of NW Argentina. Mammalian Biology, , –. and gross assimilation for Australian ecosystems: Coupled estimates N A R O S K Y , T. & Y Z U R I E T A , D. () Guía para la Identificación de las of ET and GPP. Ecohydrology: Ecosystems, Land and Water Process Aves de Argentina y Uruguay. Asociación Ornitológica del Plata, Interactions, Ecohydrogeomorphology, , e. Buenos Aires, Argentina. G U E VA R A , J.C., G R Ü N WA L D T , E.G., E S T E V E Z , O.R., B I S I G AT O , A.J., N E W B O L D , T., H U D S O N , L.N., P H I L L I P S , H.R.P., H I L L , S.L.L., C O N T U , B L A N C O , L.J., B I U R R U N , F.N. et al. () Range and livestock S., L Y S E N KO , I. et al. () A global model of the response of Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
Effect of free-ranging cattle on mammalian diversity 887 tropical and sub-tropical forest biodiversity to anthropogenic R O D R I G U E S , A.S.L., A N D E L M A N , S.J., B A K A R R , M.I., B O I T A N I , L., pressures. Proceedings of the Royal Society B: Biological Sciences, B R O O K S , T.M., C O W L I N G , R.M. et al. () Effectiveness of the , . global protected area network in representing species diversity. O J E D A , R.A. () Biodiversidad y conservación de mamíferos de la Nature, , –. interfase tropical-templada de Argentina. In Biodiversidad y Uso de R O D R Í G U E Z , J.P. & R O D R Í G U E Z -C L A R K , K.M. () Even ‘paper la Tierra: Conceptos y Ejemplos de Latinoamérica (eds parks’ are important. Trends in Ecology & Evolution, , . S.D. Matteucci, O.T. Solbrig, J. Morello & G. Halffter), pp. –. R O M E R O -M U Ñ O Z , A., B E N Í T E Z -L Ó P E Z , A., Z U R E L L , D., B A U M A N N , EUDEBA-UNESCO Colección CEA, Buenos Aires, Argentina. M., C A M I N O , M., D E C A R R E , J. et al. () Increasing synergistic O J E D A , R.A., B A R Q U E Z , R.M., S TA D L E R , J. & B R A N D L , R. () effects of habitat destruction and hunting on mammals over three Decline of mammal species diversity along the Yungas forest of decades in the Gran Chaco. Ecography, , –. Argentina. Biotropica, , –. S C H I E LT Z , J.M. & R U B E N S T E I N , D.I. () Evidence based review: O K S A N E N , J., B L A N C H E T , F.G., F R I E N D LY , M., K I N D T , R., L E G E N D R E , positive versus negative effects of livestock grazing on wildlife. What P.R., M C G L I N N , D. et al. () vegan: Community Ecology Package. do we really know? Environmental Research Letters, , . cran.r-project.org/web/packages/vegan/index.html [accessed S H E P H E R D , J.D. & D I T G E N , R.S. () Human use and small mammal December ]. communities of Araucaria frests in Neuquén, Argentina. O L S O N , D.M. & D I N E R S T E I N , E. () The Global : priority Mastozoología Neotropical, , –. ecoregions for global conservation. Annals of the Missouri Botanical S M I T H , G.T., A R N O L D , G.W., S A R R E , S., A B E N S P E R G -T R A U N , M. & Garden, , –. S T E V E N , D.E. () The effect of habitat fragmentation and P E R OV I C , P.G. () Ecología de la comunidad de félidos en las selvas livestock grazing on animal in remnants of gimlet Eucalyptus nubladas del Noroeste argentino. PhD thesis. Universidad Nacional salubris woodland in the Western Australian Wheatbelt. II. Lizards. de Córdoba, Córdoba, Argentina. The Journal of Applied Ecology, , –. P H I L L I P S , S.J., A N D E R S O N , R.P. & S C H A P I R E , R.E. () Maximum S U R AC I , J.P., G AY N O R , K.M., A L L E N , M.L., A L E X A N D E R , P., entropy modeling of species geographic distributions. Ecological B R A S H A R E S , J.S., C E N D E J A S -Z A R E L L I , S. et al. () Disturbance Modelling, , –. type and species life history predict mammal responses to humans. P H I L L I P S , S.J. & D U D Í K , M. (). Modeling of species distributions Global Change Biology, , –. with Maxent: new extensions and a comprehensive evaluation. T A B E N I , S., S P I R I T O , F. & O J E D A , R.A. () Conservation of small Ecography, , – and medium-sized mammals following native woodland regrowth: P H I L L I P S , S.J., D U D Í K , M. & S C H A P I R E , R.E. () Maxent software a case study in a long-term UNESCO Biosphere Reserve, Argentina. for modeling species niches and distributions. Version ... Journal of Arid Environments, , –. biodiversityinformatics.amnh.org/open_source/maxent [accessed V A N D E N B E R G H E , C., S M I T , C., P O H L , M., B U T T L E R , A. & May ]. F R E L É C H O U X , F. () Does the strength of facilitation by nurse P I A , M.V., L Ó P E Z , M.S. & N O VA R O , A.J. () Effects of livestock on shrubs depend on grazing resistance of tree saplings? Basic and the feeding ecology of endemic culpeo foxes (Pseudalopex culpaeus Applied Ecology, , –. smithersi) in central Argentina. Revista Chilena de Historia Natural, V E N A B L E S , W.N. & R I P L E Y , B.D. () Modern Applied Statistics with , –. S. Springer, New York, USA. P I A N K A , E.R. () Latitudinal gradients in species diversity: a review WCS & CIESIN (W I L D L I F E C O N S E R VAT I O N S O C I E T Y & C E N T E R of concepts. The American Naturalist, , –. F O R I N T E R N AT I O N A L E A R T H S C I E N C E I N F O R M AT I O N N E T W O R K ) QGIS () QGIS Geographic Information System. Open Source () Last of the Wild Project, Version , (LWP-): Global Geospatial Foundation Project. qgis.org/en/site/ [accessed May Human Footprint Dataset (Geographic). NASA Socioeconomic Data ]. and Applications Center (SEDAC), Palisades, USA. doi.org/./ Q U I R O G A , P.A., F E R N Á N D E Z , H.R., S I R O M B R A , M.D. & D O M Í N G U E Z , HMHF [accessed December ]. E. () Riparian forests and cattle management problems in Y O U N G , T.P., P A L M E R , T.M. & G A D D , M.E. () Competition and Andean subtropical streams: in the search of water quality compensation among cattle, zebras, and elephants in a semi-arid sustainability. Lilloa, , –. savanna in Laikipia, Kenya. Biological Conservation, , –. R C O R E T E A M (). R: A Language and Environment for Statistical Z H A N G , M., L I U , Z. & T E N G , L. () Seasonal habitat selection of the Computing. R Foundation for Statistical Computing, Vienna, red deer (Cervus elaphus alxaicus) in the Helan Mountains, China. Austria. Zoologia, , –. R E Y N A -H U R T A D O , R., S I M A -P A N T Í , D., A N D R A D E , M., P A D I L L A , A., Z H A N G , Y., K O N G , D., G A N , R., C H I E W , F.H.S., M C V I C A R , T.R., R E T A N A -G U I A S C O N , O. & S A N C H E Z -P I N Z Ó N , K. () Tapir Z H A N G , Q. & Y A N G , Y. () Coupled estimation of m and -day population patterns under the disappearance of free-standing water. resolution global evapotranspiration and gross primary production Therya, , –. in –. Remote Sensing of Environment, , –. Oryx, 2022, 56(6), 877–887 © The Author(s), 2022. Published by Cambridge University Press on behalf of Fauna & Flora International doi:10.1017/S0030605321001538 https://doi.org/10.1017/S0030605321001538 Published online by Cambridge University Press
You can also read