- 1School of Biological Sciences, University of Utah, Salt Lake City, UT, United States
- 2Center for Conservation Biology, Department of Biology, Stanford University, Stanford, CA, United States
- 3Department of Molecular Biology and Genetics, Koç University, Sarıyer, Istanbul, Türkiye
- 4KuzeyDoğa Society, Kars, Türkiye
- 5Cancer Science Institute of Singapore, National University of Singapore, Singapore, Singapore
- 6Department of Geology, Geophysics College of Mines and Earth Sciences, University of Utah, Salt Lake City, UT, United States
- 7300m Norte de la Escuela de Buena Vista, Barva, Heredia, Costa Rica
- 8Department of Geological Sciences, Doerr School of Sustainability, Stanford University, Stanford, CA, United States
- 9Department of Biology, University of New Mexico, Albuquerque, NM, United States
How human modification of native habitats changes the feeding patterns and nutritional ecology of tropical birds is critical to conserving avian biodiversity, but tropical bird diets are laborious to investigate using the traditional methods of diet analysis. Stable isotope analysis provides a cost-effective and efficient proxy to identify general foraging patterns, especially when dietary shifts spanning multiple trophic levels have occurred due to ecosystem disturbance or transformation. To characterize the diets of forest bird species that persist in tropical agricultural countryside, we compared feather carbon (δ13C) and nitrogen (δ15N) isotope values of four species caught and radio-tracked in a 270 hectare forest reserve, smaller forest remnants (including mature forest, secondary forest, and riparian strips), and coffee plantations in mid-elevation (ca. 800–1,400 m) southern Costa Rica. Bird habitat choice had a significant effect on diet composition as revealed by δ13C and δ15N values. Three of the four species studied showed evidence of significantly reduced consumption of invertebrates in coffee plantations, with the isotope values of two species (Tangara icterocephala and Mionectes oleaginous) indicating, by comparison, nearly a doubling of invertebrate consumption in forest remnants. Our results suggest that coffee plantations are deficient in invertebrates preferred by forest generalist birds that forage in both native forest remnants and coffee plantations. In this region, typical of mountainous American tropics, small forest remnants and a larger forest reserve provide critical dietary resources for native forest birds that utilize the agricultural countryside.
Introduction
Forests worldwide are being reduced to biologically-impoverished remnants (Laurance and Bierregaard, 1997) embedded in agricultural countryside (Daily et al., 2001; Şekercioğlu et al., 2019) – croplands, pasture, gardens, open second growth, and a scattering of forest fragments, riparian strips, and remnant trees. Even though human-dominated areas are the preferred habitat of fewer than 1% of the world’s ~ 11,000 avian species (Şekercioğlu et al., 2004), about a third make some use of such habitats (Şekercioğlu et al., 2007), including many temperate-tropical migratory birds (Blount et al., 2021) that are increasingly threatened with extinction (Horns and Şekercioğlu, 2018). Tropical agricultural countryside can even favor generalist species by providing them with additional food resources (Clough et al., 2011; Fahrig et al., 2011). Given the high rate of tropical deforestation, the extent and ecological qualities of deforested countryside will determine whether some forest species can persist in the absence of large, intact forest patches.
Nitrogen is an essential component of protein synthesis but is found in limited quantities in most ecosystems (White, 1993). Compared to tropical fruits, which average ~ 5% dry mass as protein (Bosque and Pacheco, 2000), invertebrates have substantially higher protein content (~60% of dry mass; Bell, 1990) and are key sources of nitrogen for most tropical bird species.
A limitation of protein-rich food sources could have large-scale effects on avian communities because the relative level of protein intake has important impacts on bird growth and reproduction (Izhaki, 1998; Gill et al., 2019). Proteins are critical in many cellular processes, and the growth, development and consequently the fitness of organisms may be limited by a deficiency of nitrogen in their dietary resources. The underappreciated research frameworks of ecological stoichiometry and ionomics address how such changes in the environmental nutritional supply affect the ecophysiology, behavior, health, and fitness of individuals, influencing their ecological interactions, population functioning and conservation biology (Kaspari, 2021; Filipiak and Filipiak, 2022).
Stoichiometric mismatch between metabolic demands and nutritional supply can be especially critical for juveniles (Filipiak and Filipiak, 2022). For example, even for bird species that predominantly eat fruits as adults, invertebrates are preferentially caught by adults and fed to offspring during critical periods of growth and development (Izhaki, 1998), enabling chicks to rapidly fledge in an environment full of nest predators, with some exceptions in granivorous birds (Díaz, 1996). Some tropical forest birds time their breeding with increased insect availability (Hau et al., 2000). Declines in protein sources can thus lead to reduced growth rates, increased chick mortality, declines in population recruitment, and other limitations on the growth rate and overall size of a bird population.
We used stable isotope analysis (SIA) to test the hypothesis that the breeding season diets of native forest birds that continue to persist in and adjacent to coffee plantations in southern Costa Rica are deficient in nitrogen (protein) in coffee plantations when compared to their diets in native forest remnants. Because birds have higher energy and protein needs during reproduction (Klasing, 1998) and songbirds tend to consume more arthropods during this period (Izhaki, 1998), we conducted our study to coincide with the main breeding season of our study species between March and June. By examining the carbon (δ13C) and nitrogen (δ15N) isotope values of feathers (n = 170) collected from four resident bird species, we tested the hypothesis that individual birds foraging in southern Costa Rican human-modified habitats dominated by coffee plantations consume fewer invertebrates in comparison to individuals of the same species foraging in native forest vegetation. We also collected and analyzed samples of common sources of food (fruits and invertebrates) consumed by our focal bird species (Şekercioğlu et al., 2007). We sampled birds and dietary items from two primary habitats, partially shaded coffee plantations and forest remnants that include primary forest, secondary forest, and riparian strips; hereafter referred to as “coffee” and “non-coffee” habitats, respectively. We also used radiotelemetry data to assess the overall habitat use of each individual bird and to understand the level of use of coffee versus non-coffee habitats. A total of 49 individuals from three species were intensively radio-tracked and their habitat use patterns were documented in detail (Şekercioğlu et al., 2007), which allowed for a unique coupling of isotopically derived dietary information with information on individual movement patterns.
Based on a priori knowledge of the birds’ habitat and dietary preferences, we hypothesized that there would be significant isotopic differences among bird species related to their degree of insectivory, which would mirror those recorded from field observations (Şekercioğlu et al., 2007). Coffee plantations in our study area have depauperate invertebrate communities (Goehring et al., 2002; Ricketts et al., 2002), so we also predicted that birds that spend a greater amount of time in these areas should consume fewer invertebrates than do the birds spending more time in forest remnants. The unique coupling of extensive radio-tracking data with isotopically derived dietary information on resource use and assimilation provides a rare perspective on the potential effects of landscape modification on the nutrient uptake of a tropical bird community. This is one of the few ornithological studies that combines radio tracking and SIA to address a critical conservation issue for tropical avian species.
Materials and methods
Study site
Our study area in southern Costa Rica (Figure 1), near the Las Cruces Biological Station, is a formerly forested region now dominated by coffee plantations (~50% of land cover), cattle pastures (~20%), and other human-modified environments (~10%). In these agricultural areas, large quantities of chemicals in the form of fungicides, pesticides, and fertilizers are applied on coffee plantations, and the regulations restricting the use of pesticides are not adequately enforced (Stiles and Skutch, 1989). Thus, the conversion of forests to agricultural countryside may not only reduce the abundance of invertebrates via decreases in vegetation structure and botanical diversity, but the elimination of invertebrates with pesticides can lead to the poisoning of birds and other non-target species while directly decreasing the amount of protein available to the resident bird species. Previous studies at our study site found significant decreases in the numbers of invertebrates such as Blatteria, Coleoptera, Diptera, Lepidoptera, and Orthoptera in coffee plantations in comparison to the forest remnants (Goehring et al., 2002; Ricketts et al., 2002). These insect orders are major components of local bird diets (Şekercioğlu, 2002) and such reductions could result in nitrogen limitation for birds that spend a considerable amount of their time foraging in coffee plantations.
Figure 1. Map of the study area. Dark gray represents arboreal cover (forest fragments, riparian strips, and remnant trees), gray represents second growth with few trees, and white represents human-dominated habitats (pastures, coffee plantations, other agriculture, and towns).
Study species
Catharus aurantiirostris (Orange-billed Nightingale-thrush, hereafter Catharus) is a species with low forest dependence and prefers forest edge, second growth, thickets, plantations, and gardens (Stiles and Skutch, 1989). Tangara icterocephala (Silver-throated Tanager, hereafter Tangara) has intermediate forest dependence and is a common species of forest canopy, forest edge and agricultural habitats in the study area. Tangara was the most common seed disperser visiting Miconia trees (a major food source analyzed in this study) isolated in the pastures surrounding the intact Las Cruces forest (Luck and Daily, 2003). Turdus assimilis (White-throated Thrush, hereafter Turdus) is most often found in moist forests, but also enters bordering thickets, riparian corridors, second growth and isolated fruiting trees, especially outside the breeding season (Stiles and Skutch, 1989). Finally, Mionectes oleaginous (Ochre-bellied Flycatcher, hereafter Mionectes) is most often found in the lower levels of humid forest, secondary growth, and clearings (Stiles and Skutch, 1989). Mionectes was the only species that was not radio-tracked. All four species are regularly found in forest and open countryside, feed on fruits and invertebrates, and are likely to be important dispersers of forest plants into non-forested areas. Catharus is mostly insectivorous, whereas the others are mostly frugivorous (Stiles and Skutch, 1989; Şekercioğlu et al., 2007).
We analyzed a subset of feathers from 170 individuals of these four bird species (Table 1). The feathers were part of a collection of contour and tail feathers from over 6,000 individuals of 211 Costa Rican bird species obtained during a long-term bird banding and radio telemetry study centered around the Las Cruces Biological Station in the Coto Brus canton of southern Costa Rica (Şekercioğlu et al., 2002, 2007, 2015, 2019).
Table 1. Mean (±SD) δ13C and δ15N values and the sample sizes of tail feathers of four bird species and potential food sources collected in coffee and non-coffee habitats.
Radio tracking
For three of these four species, we collected extensive radiotelemetry data on 49 individuals (Şekercioğlu et al., 2007). We used radiotelemetry data to assess the habitat use of each individual bird, especially the use of forest habitats (the mostly intact primary forest of the Las Cruces Forest Reserve, small forest remnants, and riparian strips) versus agricultural habitats (mostly coffee plantations). Various studies have shown that radio tags have minimal to no impact on bird physiology and behavior (Kenward, 2001; Wells et al., 2003), and the radio tags in this study had no observable effects on our study species (Şekercioğlu et al., 2007).
From the average radio-tracked individual, we collected ≥55 Global Positioning System (GPS) points during more than 45 h of radio-tracking distributed across 11 days. To assess the level of confidence we could place on bird capture locations for understanding habitat use, we analyzed the fit between 97 capture site data points, each of which describes the location where an individual bird was caught, and each bird’s radio tracking locations, which give detailed information on habitat use by the same individual bird. For each bird, we calculated the ratio of the number of telemetry data points collected in coffee plantations to all telemetry data points. To assess whether the capture location of the bird could predict the habitat where the bird spent most of its time according to radiotelemetry, the vegetation cover of the study area was classified as “coffee” versus “non-coffee,” each of which had ~50% cover based on extensive vegetation surveys (Şekercioğlu et al., 2007).
If an individual bird spent more than 50% of its time in coffee plantations, we interpreted it as a preference for that habitat type and vice versa for non-coffee areas. We divided the number of times capture data and radio-tracking data assigned an individual to the same habitat (coffee or non-coffee) by the total number of radio-tracked birds. A z-test was performed to see if this was significantly different than the 50% value that would be expected by chance alone. We expected the fit between the capture locations and the corresponding telemetry data points to be better than that expected by chance alone. As the coffee and non-coffee areas each represent ~50% of the study site, the expected agreement rate due to chance alone was ~50%.
For the purposes of the isotope mixing models and statistical analyses, percent agricultural use was calculated for each bird by dividing the number of telemetry points from agricultural habitats by the total number of telemetry points. An index value of 34% or less was arbitrarily taken to indicate light use of agricultural countryside, between 34 and 67% moderate use, and above 67%, heavy use. For birds that were not radio-tracked, capture data was used instead, and 0% was assigned to a bird captured in a non-coffee site, while 100% was assigned to a bird captured in a coffee site. We also calculated the percent of each bird’s home range that was covered by trees and used linear regression to analyze the relationship between percent tree cover and the isotopic composition of each bird’s feathers.
We assumed that the location of each bird during the time of feather formation (and reflected in its isotopic signatures) corresponds to the location in which the bird was caught or to the locations in which the bird spent most of its time as revealed by radio tracking studies. In addition, we assumed that a bird’s diet during the period of feather formation is representative of its typical diet. This might not be true because feather formation requires high amounts of protein and thus, during this time, birds may be preferentially consuming more protein (Earle and Clarke, 1991; Wolf et al., 2003). Nevertheless, even if an individual does modify its foraging behavior to consume more protein during the molting/re-growth period, a bird impoverished in protein post-formation was likely even more nutrient impoverished during feather formation.
Stable isotopes as avian dietary proxies
SIA has yielded important insights into the life history, ecology, migration, and evolution of birds (Hobson and Clark, 1992; Chamberlain et al., 1997; Hobson, 1999a; Rubenstein et al., 2002; Smith et al., 2002; Hobson and Bairlein, 2003; Pain et al., 2004). Specifically, carbon (δ13C) and nitrogen (δ15N) isotope analyses of bird tissues have been used to investigate trophic relationships (Hobson, 1999b; Kelly, 2000; Herrera et al., 2003) and to characterize avian diet composition (Mizutani et al., 1990; Thompson et al., 1995; Thompson and Furness, 1995; Bearhop et al., 1999; Hocking and Reimchen, 2002). The use of SIA to analyze diet has several advantages. Traditionally, the examination of regurgitates, feces, and pellets, and other laborious means have been necessary to study avian diets (Rosenberg and Cooper, 1990; Şekercioğlu et al., 2002). These techniques provide a snapshot of dietary information and typically represent food ingested over a few days prior to capture. These techniques are also highly invasive, may create intolerable levels of stress for study subjects, and can even lead to mortality (Poulin et al., 1994). Isotopic analysis of feathers, easily collected during any study that involves mist netting (Smith et al., 2002), removes the need for the laborious and invasive processes involved in obtaining regurgitates and feces, and provides direct insights into the diet of an individual bird recorded in its feathers during the period of feather growth (Mizutani et al., 1990).
Finally, quantitative estimates of dietary sources can sometimes be made with SIA by applying trophic discrimination factors to consumer tissue isotope values to account for isotopic sorting during metabolic and biosynthetic processes, and then using mixing models designed to calculate the relative proportion of dietary sources in consumer diets (Phillips and Koch, 2002; Phillips and Gregg, 2003; Parnell, 2019). In most cases, consumers are enriched in the rare heavy isotope (13C or 15N) relative to their diets by ~1–2‰ for δ13C and 3–5‰ for δ15N for each increase in trophic level (Kelly, 2000; Caut et al., 2009).
Fruit and invertebrate samples
We conducted SIA of commonly consumed fruit and invertebrate samples from the sites where the birds were captured. We sampled the fruits of Miconia trinervia (family Melastomataceae; hereafter Miconia), a regular food resource of all the study species and a plant that is commonly found in natural and human-dominated habitats in the study area. We also collected the fruits of Cecropia spp. (family Cecropiaceae; hereafter Cecropia), also fruits commonly consumed by our study species. For invertebrates, we sampled predaceous arachnids (spiders) and grazing Orthoptera (grasshoppers), which represent two trophic levels and are two of the four invertebrate groups (others include Coleoptera and Formicidae) that are most frequently encountered in the diet samples of birds from our study area we collected via regurgitation (Şekercioğlu et al., 2002). We collected and analyzed 4–25 samples of each dietary item from coffee and non-coffee areas to compare the isotopic differences between these treatments.
Stable isotope analyses
Bird tail feathers were prepared for isotope analysis by washing in a 2:1 chloroform/methanol mixture to remove surface contaminants, followed by drying for at least 8 h in a fume hood to remove solvents (Kelly, 2000). Fruits and invertebrates were freeze-dried and homogenized using a mortar and pestle. Dried samples of feathers (0.5 mg) and potential food items (invertebrates: ~0.5 mg; fruits: ~5.0 mg) were sealed in tin capsules, and δ13C or δ15N values were measured using a Costech elemental analyzer interfaced with a Finnegan Delta Plus gas source mass spectrometer in the Department of Geological and Environmental Sciences at Stanford University (Palo Alto, CA). Results are expressed as δ values, δ13C or δ15 N = 1,000 [(Rsample/Rstandard) − 1], where Rsample and Rstandard are the 13C/12C or 15N/14N ratios of the sample and standard, respectively; units are expressed as parts per thousand or per mil (‰). The standards are Vienna-Pee Dee Belemnite limestone (V-PDB) for carbon and atmospheric N2 for nitrogen. Repeated within-run measurements of a gelatin standard yielded an average standard deviation of <0.2‰ for both δ13C and δ15N values. Duplicate isotopic measurements were performed on ~20% of all unknown samples and yielded an absolute difference of 0.2‰ for both δ13C and δ15N values. We also measured the weight percent carbon and nitrogen concentrations of all samples, which are presented as nitrogen concentrations and C/N ratios in Table 1.
An analysis of variance (ANOVA) with a Tukey-HSD post hoc comparison test was used to test for significance among groups (bird species or prey types), which was assessed at an α-value of 0.05. In addition, we used a Bayesian concentration-dependent mixing model, Stable Isotope Mixing Models in R (simmr; Parnell, 2019), to estimate the relative contribution of four prey types to the diet of each bird species. Based on the results of a controlled feeding experiment on an omnivorous passerine (Pearson et al., 2003), we varied δ13C and δ15N trophic discrimination factors (TDF) among prey types to account for the differences in assimilation due to prey protein (nitrogen) content. For δ13C, we applied mean (±SD) TDFs of 3.0 ± 0.5‰ for invertebrates (spiders and grasshoppers) and 2.0 ± 0.5‰ for fruits (Cecropia and Miconia). For δ15N, we applied mean (±SD) TDFs of 3.5 ± 0.5‰ for invertebrates (spiders and grasshoppers) and 3.0 ± 0.5‰ for fruits (Cecropia and Miconia). The model also included elemental concentration data and we used mean weight percent [C] and [N] for each prey type reported in Table 1. Gelman diagnostics for all model runs produced values of one, indicating model convergence. We used a two-tailed t-test to assess the significant differences in prey contributions between habitats within species, and significance was assessed using an α-value of 0.05. The dietary percentages presented in Table 2 represent the mean source contributions (SD in parentheses) of all individuals within a species versus habitat type based on telemetry data or capture location (coffee versus non-coffee).
Table 2. Mean (±SD) contributions of potential prey to the diets of four bird species collected in coffee and non-coffee habitats.
Results
Habitat use
The agreement between 97 capture site data points and telemetry data points was 68%. This was significantly higher (z = 3.70; p = 0.0002) than would be expected based on chance alone, meaning that the capture location was a good predictor of the habitat use of an individual bird over time. Hence, capture location data were used in data analyses for birds that were not radio-tracked. We also calculated the average percent (±SD) agricultural use of radio-tracked birds to determine the extent of use of agricultural countryside. Percent agricultural use was 14% (±5.9%) for Mionectes, 22% (±5.0%) for Turdus, 33% (±4.9%) for Catharus, and 34% (±6.9%) for Tangara, indicating that these species made relatively low use of the agricultural countryside. Percent agricultural use for Mionectes was calculated using only the capture locations since this species was not radio-tracked.
Foraging observations
Each species differed from the other two in its resource use (Figure 2), based on the proportion of foraging observations on the resources that each species was observed to exploit (all χ2 > 120; all p < 0.001; Figure 2). Catharus (Figure 2), although observed to forage on 14 plant taxa, foraged predominantly on invertebrates, either directly or on substrates (i.e., soil, rotting fruit) where invertebrates are commonly found. Turdus (Figure 2) and Tangara (Figure 2) were more frugivorous and foraged on 71 and 45 plant taxa, respectively. Fruits of Cecropia peltata, Ficus spp., and Syzygium malacanse were favorite dietary items for both species.
Figure 2. The relative percentage of foraging observations on common food resources (each type comprising >5% of all total foraging observations for three of the four bird species examined in this study). No foraging data were collected for Mionectes. Data are from Şekercioğlu et al., 2007.
Isotopic composition of dietary sources
There were significant differences between δ15N values, δ13C values, and weight percent carbon:nitrogen (C/N) concentrations among potential food sources (Tukey HSD, p < 0.01; Table 1; Figure 3). Spiders collected within coffee habitats had significantly higher δ15N values than any other food source (Figure 3). Miconia fruits collected in non-coffee habitats had significantly lower δ15N values than any other food source. δ15N and δ13C values of Cecropia fruits were not statistically distinguishable from those of grasshoppers in either habitat. For both habitats, spiders had significantly higher values than did grasshoppers, but only grasshoppers collected in non-coffee habitats had significantly higher δ15N values than non-coffee Miconia fruits. δ15N values of grasshoppers were indistinguishable from δ15N values of Miconia fruits in coffee plantations. Lastly, spiders collected in coffee habitat had significantly higher δ15N values than spiders from non-coffee habitats.
Figure 3. Mean δ15N (A) and δ13C (B) values of the four bird species and common dietary resources collected in coffee versus forested non-coffee habitats. Error bars denote one standard deviation from the mean. Refer to Table 1 for the sample sizes. Lowercase letters signify results of a one-way ANOVA comparing mean values among prey types and uppercase letters signify ANOVA results for the comparisons of a single prey type between habitats.
Like δ15N results, spiders collected within coffee habitats also had significantly higher δ13C values than any other food source (Figure 3), with the only exception being spiders collected from non-coffee habitats. The δ13C composition of grasshoppers collected in coffee or non-coffee habitats was indistinguishable from that of Cecropia fruits. The mean δ13C value of Cecropia fruits from either habitat was indistinguishable from that of Miconia fruits. However, Miconia fruits collected from non-coffee habitats had significantly lower δ13C values than grasshoppers from either habitat. Lastly, within-type comparisons of food sources collected in coffee versus non-coffee habitats showed no significant differences in δ13C values except for spiders, which (like δ15N) had higher δ13C values in coffee versus non-coffee habitats (Figure 3B).
Grasshoppers had significantly higher variance in δ13C and δ15N values than any other food source type (Table 1; Figure 3). There were no differences in the concentrations of carbon [C] and nitrogen [N] between samples of particular food types collected in coffee versus non-coffee habitats (Table 1). Invertebrates (spiders and grasshoppers) had significantly higher [N] than fruits, resulting in significant differences in C/N ratios between these two general types of food (Table 1).
Isotopic composition of bird feathers
We measured the δ13C and δ15N values of feathers collected from 63 Catharus, 38 Turdus, 34 Tangara, and 36 Mionectes individuals (Table 1). At the species level, sex, wingspan, and body mass did not have significant effects on δ13C and δ15N values (ANOVA, p > 0.05). For intraspecific comparisons (uppercase letters, Figure 3), the location of capture had a significant effect on feather δ15N values for Catharus, Mionectes, and Tangara, where individuals captured in coffee habitats had significantly higher δ15N values than individuals from non-coffee habitats (Figure 3; ANOVA). For δ13C, intraspecific comparisons showed significant effects of habitat for only one species. Turdus individuals captured in coffee plantations had significantly higher δ13C values than those caught in non-coffee habitats (Figure 3).
Interspecific ANOVA comparisons (lower case letters, Figure 3) show significant differences in mean δ15N but not mean δ13C values among species. Catharus and Mionectes have significantly higher δ15N values than do Tangara. Furthermore, Turdus captured in non-coffee habitats have significantly lower δ15N values than Catharus captured in either habitat type. Likewise, Tangara captured in non-coffee habitats have lower δ15N values than Turdus caught in either habitat type.
Results of the δ13C and δ15N concentration-dependent mixing model are presented in Table 2. Consumption of invertebrates was higher in forested versus coffee habitats for three of the four focal bird species (Figure 4). Specifically, the mean contribution of spiders to the diets of Catharus, Tangara, and Mionectes was higher in forested versus coffee habitats. The mean contribution of grasshoppers to the diet of Mionectes was higher in forested versus coffee habitats. Conversely, the consumption of Miconia fruits was higher in coffee versus forested habitats for all four bird species, but the consumption of Cecropia fruits by Tangara was higher in forested versus coffee habitats. These results mostly confirm our hypotheses. All species except Turdus consume substantially more invertebrates in forested non-coffee habitats than they do in coffee plantations, with Tangara and Mionectes consuming nearly double the invertebrates in forest remnants than they do in coffee plantations (Figure 4). Patterns among species are similar; overall, Mionectes has the highest invertebrate consumption and Turdus and Tangara have the lowest (Figure 4).
Figure 4. Mean diet proportions of potential prey types consumed by the four bird species in coffee (C) and forested non-coffee (NC) habitats.
Of the 170 individual birds sampled for isotopes, 19 Catharus, 16 Tangara, and 14 Turdus were radio-tracked. The percentage of trees in their home ranges had no significant effect (r2 < 0.10; p > 0.05 for both isotope systems) on feather δ13N (Figure 5A) or δ13C values (Figure 5B) for any of the species that were radio-tracked.
Figure 5. Effect of percent tree cover in radio-tracked birds’ home ranges on the nitrogen (A) and carbon (B) isotopic composition of bird feathers. Regression statistics for the entire bird community (rather than single species) are provided because bivariate analyses of isotopic composition versus percent tree cover for each species were similar (r2 < 0.10; p > 0.05 for both δ13C and δ15N).
Discussion
Our results indicate that three of the four bird species examined in this study fed on fewer invertebrates in coffee plantations than in forested habitats comprised mostly of forest remnants as well as riparian strips and a larger (270 ha) forest reserve. In the mosaic of habitat types in our study area, extensive coffee plantations and other human-modified habitats have a significant negative effect on the diet composition of the four frugivore-insectivore bird species we examined in this study. Individual birds that were captured in coffee plantations, or known (via radiotelemetry data) to spend a considerable portion of time in coffee plantations, may be able to make up for some of their resulting deficiency in invertebrates because their mobility allows them to move among coffee plantations and the surrounding forest remnants, including primary forests, small forest fragments, secondary forest, and/or riparian strips adjacent to coffee-dominated habitat. This explanation is supported by our radio-tracking data, which showed that Tangara and Turdus were highly mobile and remained in coffee plantations for relatively short periods of time in comparison to the amount of time spent in forest remnants (Şekercioğlu et al., 2007). Moreover, we were able to capture a greater number of Tangara and Turdus individuals in non-coffee habitats versus coffee plantations in comparison to Catharus. Of the species radio tracked, Tangara and Turdus have significantly larger home ranges than Catharus. Average home range estimates for Tangara range from ~16–22 ha and between ~30–70 ha for Turdus (Şekercioğlu et al., 2007), but only ~1.7–3.6 ha for Catharus, depending on whether minimum convex polygon versus kernel densities were used to calculate home range size. This suggests that highly mobile birds that can rapidly move between forest remnants in agricultural countryside may have greater access to nitrogen-rich foods (invertebrates), but less mobile species such as Catharus may be stuck between paying a high energetic cost by constantly moving between forest remnants or getting less invertebrate protein in coffee plantations.
Catharus has relatively low forest dependence compared to most other forest birds in the region (Stiles, 1985), and its natural preference for forest gaps, edges, and second growth has resulted in its successful adaptation to coffee plantations (Şekercioğlu et al., 2007). Most other insectivorous birds in the region are more forest dependent (e.g., Lophotriccus pileatus, Platyrinchus mystaceus, Thamnophilus bridgesi) and generally avoid the nonforest matrix dominated by coffee plantations (Şekercioğlu, 2002; Şekercioğlu et al., 2015, 2019). Future tracking and dietary studies should focus on these species, although capturing sufficient numbers of more forest-dependent species in coffee plantations will remain a challenge.
While our isotopic results show differences in the percentages of invertebrates consumed among species, all the species examined are generalist omnivores and likely have broader dietary tolerances than do the more specialist species. In comparison to many tropical forest insectivore specialists that are particularly sensitive to habitat disturbance and fragmentation and are rarely found outside forested areas (Şekercioğlu and Sodhi, 2007; Sherry et al., 2020), the relatively mobile and omnivorous bird species we examined here are more capable of using the agricultural countryside effectively to obtain the protein they require. Nonetheless, our results show that even these more adaptable species consume substantially fewer invertebrates in coffee plantations than they do in forest remnants.
Intraspecific comparisons showed that capture habitat sometimes had a significant effect on δ15N values. This effect appears to be mainly driven by the differences in spider δ15N values (Figure 3) and not in the differences in the proportion of dietary items consumed in different habitats. We do not know why spiders captured in coffee plantations have significantly higher δ15N values than do the spiders in non-coffee habitats. However, it is unlikely to be due to the influence of anthropogenic nitrogenous inputs (e.g., fertilizers, atmospheric deposition) in agricultural areas because other dietary items (grasshoppers, Cecropia, and Miconia) at lower trophic levels collected in coffee plantations do not have significantly higher δ15N values than those collected in forested habitats. Spiders are generalist insectivores (Sanders and Platner, 2007) and could forage at multiple trophic levels, so the observed differences in mean δ15N values could relate to the differences in their diet composition in coffee versus forested habitats.
Our results indicate significant differences in diet composition among species that agree with these birds’ foraging preferences based on our observational evidence (Şekercioğlu et al., 2007). These results mostly confirm our hypotheses. All species except Turdus consume substantially more invertebrates in forested non-coffee habitats than they do in coffee plantations, with Tangara and Mionectes consuming nearly double the invertebrates in forest remnants than they do in coffee plantations (Figure 4). Patterns among species are similar; overall, Mionectes has the highest invertebrate consumption and Turdus and Tangara have the lowest (Figure 4).
One limitation of our study was that we were only able to analyze four food sources that likely represent a small subset of possible resources consumed by these bird species, especially in southern Costa Rica’s very diverse tropical habitats. Although Cecropia and Miconia fruits are consumed by these species regularly, and spiders and grasshoppers make up two of the top four invertebrate types eaten by these species in our study area (Şekercioğlu et al., 2002), these generalist bird species likely feed on a wide range of fruits and invertebrates in this diverse landscape mosaic. When dealing with omnivorous species inhabiting diverse ecosystems, it is sometimes difficult to isotopically characterize all potential food sources. Thus, we designed our sampling protocol to capture potential items representative of three trophic levels, including primary producers (fruits), primary consumers (grasshoppers), and secondary consumers (spiders) to specifically evaluate the relative proportion of nitrogen-rich (invertebrates) versus nitrogen-poor (fruits) resources in individual bird diets.
The second potential source of error in our analyses is associated with the precision of our habitat use classifications for individual birds that were not radio-tracked, especially in relation to the timing of feather formation in the species examined. While our statistical analysis suggests that the location of capture is a dependable indicator of individual habitat use, the high mobility of some species may bias our coffee versus non-coffee classification scheme if all individuals (within a species) are utilizing habitats in equal proportion to the average habitat use reflected in the radio-tracking data (Şekercioğlu et al., 2007). A comparison of habitat use data derived from radiotelemetry, however, suggests a considerable amount of variance in the types of habitats utilized and the home range sizes of individual birds. For instance, some Catharus and Turdus individuals spend a lot of time in agricultural areas even though average habitat use at the species level suggests low use of this habitat type in comparison to forested habitats (Figure 5). Furthermore, just how persistent individual movement patterns are over the course of feather growth remains unresolved. For individuals that were radio-tracked, each bird was followed for an average of 45 h distributed across an average of 11 days. We must consider the possibility that the temporal persistence of individual habitat preferences is shorter than the average duration of feather growth. For the four resident tropical forest bird species examined, feather molt and regrowth occur from June to September, which coincides with the wet season in southern Costa Rica. The radiotelemetry data and the dietary samples for isotopic analysis were collected during the months of March–June. If the individual habitat preferences, as recorded in the radio-tracking data, are more ephemeral than the duration of time required to regrow tail feathers during the molting period, than it is possible that the telemetric and isotopic data do not correlate in time (as seen in Figure 5).
Conclusion
Our results indicate that coffee plantations in southern Costa Rica may be deficient in the invertebrates preferred by generalist birds that use both forest remnants and coffee plantations. Forest remnants surrounding coffee plantations in our study area provide important dietary resources for native forest birds that utilize the agricultural countryside. At least three of the four species examined in this study may be obtaining most of their food from forest remnants and spending time in coffee plantations only while moving between remnants. This interpretation is supported by our extensive radio-tracking data (Şekercioğlu et al., 2007), which showed that most of these same forest species have persisted in agricultural landscapes by being highly mobile foragers that feed on dozens of different native plant and invertebrate species scattered in small forest remnants. Landscape context is increasingly shown to be of paramount importance for biodiversity conservation in human-dominated, agricultural landscapes (Kati and Şekercioğlu, 2006; Concepción et al., 2008; Şekercioğlu, 2009; Martin et al., 2019; Meier et al., 2022) and our findings support this. Conservation strategies should focus on increasing the connectivity of forest remnants – intact forest, secondary growth, riparian strips – in tropical agricultural countryside, and it is urgent to prioritize the conservation and regeneration of forest remnants in increasingly human-dominated agricultural areas that continue to replace the world’s most biodiverse tropical forests.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was reviewed and approved by Stanford University IACUC.
Author contributions
ÇŞ developed the idea for the paper, conceptualized and designed the study, with support from SN, PC, GD, and PE. GD, PE, and ÇŞ acquired funding and resources for the field work. ÇŞ and FB conducted the field work and collected the radio tracking data, feather and diet samples. ÇŞ analyzed the radio tracking and habitat use data. MF conducted isotopic analyses on the diet and feather samples, with help from SN and PC. PC provided the mass spectrometer and other resources for isotopic analyses. ÇŞ, MF, and SN and led the writing of the manuscript. PC, TC, GD, and PE reviewed the manuscript and provided intellectual insights. All authors contributed to the article and approved the submitted version.
Funding
We are grateful to H. Bing, P. Bing, T. Brokaw, Center for Latin American Studies, W. Loewenstern, Koret Foundation, Moore Family Foundation, National Geographic Society (grants #7730-04 and #8411-08), Wildlife Conservation Society, Winslow Foundation, and the Carnegie Institution for Science for providing financial support.
Acknowledgments
We thank the Costa Rican government (MINAE) and the Organization for Tropical Studies for allowing us to work at the Las Cruces Biological Research Station.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Bearhop, S., Thompson, D. R., Waldron, S., Russell, I. C., Alexander, G., and Furness, R. W. (1999). Stable isotopes indicate the extent of freshwater feeding by cormorants (Phalacrocorax carbo) shot at island fisheries in England. J. Appl. Ecol. 36, 75–84. doi: 10.1046/j.1365-2664.1999.00378.x
Bell, G. P. (1990). Birds and mammals on an insect diet: a primer on diet composition analysis in relation to ecological energetics. Stud. Avian Biol. 13, 416–422.
Blount, D., Horns, J. J., Kittelberger, K., Neate-Clegg, M. H. C., and Şekercioğlu, Ç. H. (2021). Avian use of agricultural areas as migration stopover sites: a review of crop management practices and ecological correlates. Front. Ecol. Evol. 9:e664764. doi: 10.3389/fevo.2021.650641
Bosque, C., and Pacheco, M. A. (2000). Dietary nitrogen as a limiting nutrient in frugivorous birds. Rev. Chil. Hist. Nat. 73, 441–450. doi: 10.4067/S0716-078X2000000300007
Caut, S., Angulo, E., and Courchamp, F. (2009). Variation in discrimination factors (Δ15N and Δ13C): the effect of diet isotopic values and applications for diet reconstruction. J. Appl. Ecol. 46, 443–453. doi: 10.1111/j.1365-2664.2009.01620.x
Chamberlain, C. P., Blum, J. D., Holmes, R. T., Feng, X., Sherry, T. W., and Graves, G. R. (1997). The use of isotope tracers for identifying populations of migratory birds. Oecologia 109, 132–141. doi: 10.1007/s004420050067
Clough, Y., Barkmann, J., Juhrbandt, J., Kessler, M., Wanger, T. C., Anshary, A., et al. (2011). Combining high biodiversity with high yields in tropical agroforests. Proc. Natl. Acad. Sci. 108, 8311–8316. doi: 10.1073/pnas.1016799108
Concepción, E. D., Díaz, M., and Baquero, R. A. (2008). Effects of landscape complexity on the ecological effectiveness of Agri-environment schemes. Landsc. Ecol. 23, 135–148. doi: 10.1007/s10980-007-9150-2
Daily, G. C., Ehrlich, P. R., and Sanchez-Azofeifa, G. A. (2001). Countryside biogeography: use of human-dominated habitats by the avifauna of southern Costa Rica. Ecol. Appl. 11, 1–13. doi: 10.1890/1051-0761(2001)011[0001:CBUOHD]2.0.CO;2
Díaz, M. (1996). Food choice by seed-eating birds in relation to seed chemistry. Comp. Biochem. Physiol. A Physiol. 113, 239–246. doi: 10.1016/0300-9629(95)02093-4
Earle, K. E., and Clarke, N. R. (1991). The nutrition of the budgerigar (Melopsittacus-undulatus). J. Nutr. 121, S186–S192. doi: 10.1093/jn/121.suppl_11.S186
Fahrig, L., Baudry, J., Brotons, L., Burel, F. G., Crist, T. O., Fuller, R. J., et al. (2011). Functional landscape heterogeneity and animal biodiversity in agricultural landscapes. Ecol. Lett. 14, 101–112. doi: 10.1111/j.1461-0248.2010.01559.x
Filipiak, M., and Filipiak, Z. M. (2022). Application of ionomics and ecological stoichiometry in conservation biology: nutrient demand and supply in a changing environment. Biol. Conserv. 272:e109622:109622. doi: 10.1016/j.biocon.2022.109622
Gill, F. B., Prum, R. O., and Robinson, S. K. (2019). Ornithology, 4th (New York: W.H. Freeman, Macmillan Learning).
Goehring, D. M., Daily, G. C., and Şekercioğlu, Ç. H. (2002). Distribution of ground-dwelling arthropods in tropical countryside habitats. J. Insect Conserv. 6, 83–91. doi: 10.1023/A:1020905307244
Hau, M., Wikelski, M., and Wingfield, J. C. (2000). Visual and nutritional food cues fine-tune timing of reproduction in a neotropical rainforest bird. J. Exp. Zool. 286, 494–504. doi: 10.1002/(SICI)1097-010X(20000401)286:5<494::AID-JEZ7>3.0.CO;2-3
Herrera, L. G., Hobson, K. A., Rodríguez, M., and Hernandez, P. (2003). Trophic partitioning in tropical rainforest birds: insights from stable isotope analysis. Oecologia 136, 439–444. doi: 10.1007/s00442-003-1293-5
Hobson, K. A. (1999a). Tracing origins and migrations of wildlife using stable isotopes: a review. Oecologia 120, 314–326. doi: 10.1007/s004420050865
Hobson, K. A. (1999b). Stable-carbon and nitrogen isotope ratios of songbird feathers grown in two terrestrial biomes: implications for evaluating trophic relationships and breeding origins. Condor 101, 799–805. doi: 10.2307/1370067
Hobson, K. A., and Bairlein, F. (2003). Isotopic fractionation and turnover in captive garden warblers (Sylvia borin): implications for delineating dietary and migratory associations in wild passerines. Can. J. Zool. 81, 1630–1635. doi: 10.1139/z03-140
Hobson, K. A., and Clark, R. G. (1992). Assessing avian diets using stable isotopes II: factors influencing diet-tissue fractionation. Condor 94, 189–197. doi: 10.2307/1368808
Hocking, M. D., and Reimchen, T. E. (2002). Salmon-derived nitrogen in terrestrial invertebrates from coniferous forests of the Pacific northwest. BMC Ecol. 2:4. doi: 10.1186/1472-6785-2-4
Horns, J. J., and Şekercioğlu, Ç. H. (2018). Conservation of migratory species. Curr. Biol. 28, R980–R983. doi: 10.1016/j.cub.2018.06.032
Izhaki, I. (1998). Essential amino acid composition of fleshy fruits versus maintenance requirements of passerine birds. J. Chem. Ecol. 24, 1333–1345. doi: 10.1023/A:1021274716062
Kaspari, M. (2021). The invisible hand of the periodic table: how micronutrients shape ecology. Annu. Rev. Ecol. Evol. Syst. 52, 199–219. doi: 10.1146/annurev-ecolsys-012021-090118
Kati, V., and Şekercioğlu, Ç. H. (2006). Diversity, ecological structure, and conservation of the land bird community of Dadia reserve, Greece. Divers. Distrib. 12, 620–629. doi: 10.1111/j.1366-9516.2006.00288.x
Kelly, J. F. (2000). Stable isotopes of carbon and nitrogen in the study of avian and mammalian trophic ecology. Can. J. Zool. 78, 1–27. doi: 10.1139/z99-165
Laurance, W. F., and Bierregaard, R. O. (1997). Tropical Forest Remnants: Ecology, Management, and Conservation of Fragmented Communities (Chicago: University of Chicago Press).
Luck, G. W., and Daily, G. C. (2003). Tropical countryside bird assemblages: richness, composition, and foraging differ by landscape context. Ecol. Appl. 13, 235–247. doi: 10.1890/1051-0761(2003)013[0235:TCBARC]2.0.CO;2
Martin, E. A., Dainese, M., Clough, Y., Báldi, A., Bommarco, R., Gagic, V., et al. (2019). The interplay of landscape composition and configuration: new pathways to manage functional biodiversity and agroecosystem services across Europe. Ecol. Lett. 22, 1083–1094. doi: 10.1111/ele.13265
Meier, E. S., Lüscher, G., and Knop, E. (2022). Disentangling direct and indirect drivers of farmland biodiversity at landscape scale. Ecol. Lett. 25, 2422–2434. doi: 10.1111/ele.14104
Mizutani, H., Fukuda, M., and Wada, E. (1990). Carbon isotope ratio of feathers reveals feeding behavior of cormorants. Auk 107, 400–403. doi: 10.2307/4087626
Pain, D. J., Green, R. E., Gieing, B., Kozulin, A., Poluda, A., Ottosson, U., et al. (2004). Using stable isotopes to investigate migratory connectivity of the globally threatened aquatic warbler Acrocephalus paludicola. Oecologia 138, 168–174. doi: 10.1007/s00442-003-1416-z
Parnell, A. (2019). Simmr: A Stable Isotope Mixing Model. R package version 0.4.6.9000. Available at: https://CRAN.R-project.org/package=simmr (Accessed January 11, 2023).
Pearson, S. F., Levey, D. J., Greenberg, C. H., and Martinez del Rio, C. (2003). Effects of elemental composition on the incorporation of dietary nitrogen and carbon isotopic signatures in an omnivorous songbird. Oecologia 135, 516–523. doi: 10.1007/s00442-003-1221-8
Phillips, D. L., and Gregg, J. W. (2003). Source partitioning using stable isotopes: coping with too many sources. Oecologia 136, 261–269. doi: 10.1007/s00442-003-1218-3
Phillips, D. L., and Koch, P. L. (2002). Incorporating concentration dependence in stable isotopes: a critique. Oecologia 130, 114–125. doi: 10.1007/s004420100786
Poulin, B., Lefebvre, G., and McNeil, R. (1994). Effect and efficiency of tartar emetic in determining the diet of tropical land birds. Condor 96, 98–104. doi: 10.2307/1369067
Ricketts, T. H., Daily, G. C., and Ehrlich, P. R. (2002). Does butterfly diversity predict moth diversity? Testing a popular indicator taxon at local scales. Biol. Conserv. 103, 361–370. doi: 10.1016/S0006-3207(01)00147-1
Rosenberg, K. V., and Cooper, R. J. (1990). Approaches to avian diet analysis. Stud. Avian Biol. 13, 80–90.
Rubenstein, D. R., Chamberlain, C. P., Holmes, R. T., Ayres, M. P., Waldbauer, J. R., Graves, G. R., et al. (2002). Linking breeding and wintering ranges of a migratory songbird using stable isotopes. Science 295, 1062–1065. doi: 10.1126/science.1067124
Sanders, D., and Platner, C. (2007). Intraguild interactions between spiders and ants and top-down control in a grassland food web. Oecologia 150, 611–624. doi: 10.1007/s00442-006-0538-5
Şekercioğlu, Ç. H. (2002). Forest fragmentation hits insectivorous birds hard. Dir. Sci. 1, 62–64. doi: 10.1100/tsw.2002.190
Şekercioğlu, Ç. H. (2009). Tropical conservation: riparian corridors connect fragmented forest populations. Curr. Biol. 19, R210–R213. doi: 10.1016/j.cub.2009.01.006
Şekercioğlu, Ç. H., Daily, G. C., and Ehrlich, P. R. (2004). Ecosystem consequences of bird declines. Proc. Natl. Acad. Sci. 101, 18042–18047. doi: 10.1073/pnas.0408049101
Şekercioğlu, Ç. H., Ehrlich, P. R., Daily, G. C., Aygen, D., Goehring, D., and Sandi, R. F. (2002). Disappearance of insectivorous birds from tropical forest fragments. Proc. Natl. Acad. Sci. 99, 263–267. doi: 10.1073/pnas.012616199
Şekercioğlu, Ç. H., Loarie, S., Oviedo-Brenes, F., Daily, G. C., and Ehrlich, P. R. (2007). Persistence of forest birds in tropical countryside. Conserv. Biol. 21, 482–494. doi: 10.1111/j.1523-1739.2007.00655.x
Şekercioğlu, Ç. H., Loarie, S., Oviedo-Brenes, F., Mendenhall, C. D., Daily, G. C., and Ehrlich, P. R. (2015). Tropical countryside riparian corridors provide critical habitat and connectivity for seed-dispersing forest birds in a fragmented landscape. J. Ornithol. 156, 343–353. doi: 10.1007/s10336-015-1299-x
Şekercioğlu, Ç. H., Mendenhall, C. D., Oviedo-Brenes, F., Horns, J. J., Ehrlich, P. R., and Daily, G. C. (2019). Long-term declines in bird populations in tropical agricultural countryside. Proc. Natl. Acad. Sci. 116, 9903–9912. doi: 10.1073/pnas.1802732116
Şekercioğlu, Ç. H., and Sodhi, N. S. (2007). Conservation biology: predicting birds’ responses to forest fragmentation. Curr. Biol. 17, R838–R840. doi: 10.1016/j.cub.2007.07.037
Sherry, T. W., Kent, C. M., Sánchez, N. V., and Şekercioğlu, Ç. H. (2020). Insectivorous birds in the Neotropics: ecological radiations, specialization, and coexistence in species-rich communities. Auk 137:ukaa049. doi: 10.1093/auk/ukaa049
Smith, T. B., Marra, P. P., Webster, M. S., Lovette, I., Gibbs, H. L., Holmes, R. T., et al. (2002). A call for feather sampling. Auk 20, 218–221. doi: 10.2307/4090162
Stiles, F. G. (1985). “Conservation of forest birds in Costa Rica: problems and perspectives” in Conservation of Tropical Forest Birds. eds. A. W. Diamond and T. E. Lovejoy (Cambridge, UK: International Council for Bird Preservation), 141–170.
Stiles, F. G., and Skutch, A. F. (1989). A Guide to the Birds of Costa Rica (New York: Cornell University Press).
Thompson, D. R., and Furness, R. W. (1995). Stable isotope of carbon and nitrogen in feathers indicate seasonal dietary shifts in northern fulmars. Auk 112, 493–498. doi: 10.2307/4088739
Thompson, D. R., Furness, R. W., and Lewis, S. A. (1995). Diets and long-term changes in δ15N and δ13C values in northern fulmars (Fulmarus glacialis) from two Northeast Atlantic colonies. Mar. Ecol. Prog. Ser. 125, 3–11. doi: 10.3354/meps125003
Wells, K. M. S., Washburn, B. E., Millspaugh, J. J., Ryan, M. R., and Hubbard, M. W. (2003). Effects of radio-transmitters on fecal glucocorticoid levels in captive dickcissels. Condor 105, 805–810. doi: 10.1093/condor/105.4.805
White, T. C. R. (1993). The Inadequate Environment: Nitrogen and the Abundance of Animals (Heidelberg, Springer Berlin).
Keywords: ornithology, avian ecology, conservation biology, tropical biology, ecosystem services, deforestation, neotropics, stable isotope ecology
Citation: Şekercioğlu ÇH, Fullwood MJ, Cerling TE, Brenes FO, Daily GC, Ehrlich PR, Chamberlain P and Newsome SD (2023) Using stable isotopes to measure the dietary responses of Costa Rican forest birds to agricultural countryside. Front. Ecol. Evol. 11:1086616. doi: 10.3389/fevo.2023.1086616
Edited by:
Todd Jason McWhorter, University of Adelaide, AustraliaReviewed by:
Michał Filipiak, Jagiellonian University, PolandMario Diaz, Spanish National Research Council (CSIC), Spain
Copyright © 2023 Şekercioğlu, Fullwood, Cerling, Brenes, Daily, Ehrlich, Chamberlain and Newsome. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Çağan H. Şekercioğlu, c.s@utah.edu