Seasonal and Interspeciﬁc Variation in Frugivory by a Mixed Resident-Migrant Overwintering Songbird Community

: Many temperate passerine bird species switch from diets of mostly invertebrates in the spring and summer to diets that include fruit and seeds in the fall and winter. However, relatively few studies have quantiﬁed diet composition or the extent of seasonal shifts during the non-breeding period, particularly among species and across communities with both residents and migrants. We measured carbon and nitrogen stable isotope values in food items (fruits, C 3 and C 4 seeds, and insects from various trophic levels and plant communities) and in multiple tissues (feathers and plasma/whole blood) from 11 species of songbirds wintering in the southeastern U.S. We combined these diet and tissue values with empirically derived discrimination factors and used concentration-dependent mixing models to quantify seasonal diet shifts. We also validated mixing model results with data from fecal samples. Diets in this bird community, as delineated N and C isotopic space, diverged in the fall and winter relative to the summer as consumption of fruits and seeds increased. Across this songbird community, estimated contributions of fruit to plasma/whole blood increased from 16.2 ± 7.5% in the fall (mean ± SD; range: 4–26%) to 21.7 ± 10.3% (range: 9–37%) in the winter, while contributions of seeds increased from 29.4 ± 2.6% (range: 28–32%) in the fall to 36.6 ± 4.8% (range: 32–42%) in the winter. Fecal data showed qualitatively similar trends to mixing models, but consistently estimated higher contributions of fruit. Our work indicates that fruits and seeds constitute substantial sources of sustenance for non-breeding songbirds, there is considerable separation of resource use among species in the fall and winter, and fecal estimates of contributions to songbird tissues should be interpreted cautiously.


Introduction
Animal diets are constrained by the availability of different food items, so many animals seasonally shift their diets as resource availability changes [1][2][3]. Because diet is a key interface between an organism and its environment [4,5], investigating the direction and scope of dietary shifts is often essential for understanding both the influence of the environment on individual condition and the influence of species on ecological communities [6][7][8]. More specifically, understanding diet shifts can help illuminate aspects of physiological function such as immune responses [9] or fat accumulation [10] and also aspects of community structure such as niche partitioning [11] and predation pressure [2]. However, despite the potential value of studying diet shifts, diets have not often been described in both qualitative and quantitative terms for multiple species that compose a bird community.
A relatively well-established example of a highly predictable diet shift is the seasonal transition by many species of passerine birds in temperate regions from invertebrate prey in the spring and summer to alternate resources (e.g., fruits and seeds) in the fall and winter [12][13][14]. This shift is likely driven by reduced biomass of invertebrates during fall and winter and a concurrent increase in availability of fruits and seeds [15][16][17], which are generally conspicuous, energy rich, and easily "captured" [14,18,19]. The most common methods for studying this phenomenon have included both direct and indirect observations of fruit consumption by birds and analysis of the contents of fecal samples. These methods have been used to successfully identify the species of birds that engage in frugivory, the kinds of fruit they consume, and the timing of fruit consumption [13,16,17,20,21]. Thus, prior work has made considerable progress in characterizing networks of bird-plant interactions and other community-level patterns of frugivory involving songbirds.
Methodological constraints have limited the investigation of seasonal frugivory by songbirds in a manner that is both quantitative and readily integrated with organismal biology. For example, whereas the impact of frugivory on fruits is highly visible, direct observations and fecal analysis can only provide information on what is ingested, not what is assimilated by songbirds [7,18]. Moreover, direct observations can be affected by bias in detectability of different food items and effort [20], and fecal analysis is complicated by variation in digestibility, nutrient content, and retention time of food items [22][23][24]. Many of these issues can be addressed through the use of stable isotope diet reconstruction methods, which provide information on both assimilated nutrients and establish a common currency that can be used to calculate relative contributions of food items to consumers (e.g., [19,25,26]). This is particularly the case for studies using Bayesian mixing models, which are able to produce robust estimates of contributions by different food items to animal tissues through the propagation of the underlying uncertainty in various isotopic measurements [27,28]. For example, Gagnon and Hobson [29] were able to document a surprisingly high reliance on fruits by American Redstarts (Setophaga ruticilla) and blackand-white Warblers (Mniotilta varia), and surprisingly low reliance on fruits by American Robins (Turdus migratorius) and Cedar Waxwings (Bombycilla cedrorum) during both spring and fall migration. Meanwhile, Gómez et al. [30] estimated that fruits contributed approximately 14% of the diet of Gray-cheeked Thrushes (Catharus minimus) during spring migration and Ruhl et al. [31] estimated that Rubus fruits contributed 9-13% of the diets of Scarlet Tanagers (Piranga olivacea) during the post-fledgling period. Such quantitative estimates of assimilated diet and seasonal frugivory are still relatively rare, in large part due to the challenge of obtaining high-quality information on potential food sources and diet-tissue discrimination factors [32,33]. Furthermore, few studies document variation in resource use within and among species, despite its known effect on both the structure and stability of ecological communities [34,35] and potential insights into behavioral variation among individuals [36][37][38].
We used measurements of stable isotope ( 13 C and 15 N) enrichment of tissues and three types of food resource (arthropods, fruits, and seeds) to estimate resource use by a songbird community composed of eleven resident and migratory species during fall and winter at a site in South Carolina, USA. Our primary objective was to provide estimates of fruit and seed use by this non-breeding bird community, with those estimates made more organismally relevant by focusing on assimilated nutrients and made more robust by providing contemporaneous measurements of the isotope values of potential food items in our study area. These estimates also allowed us to directly compare estimates of diet from the stable isotope mixing models with those from traditional fecal analyses and assess the scope of differences between these methods. As a secondary objective, we characterized changes in the composition and variability of songbird diets between seasons, where changes in resource availability might drive changes in diet. In particular, we expected fruit and seed contributions to increase as insect abundance decreased from summer to fall to winter, and we expected diet variability to increase as resource use diverged across the songbird community over time. Finally, our database of isotope values for a diverse array of prey items may be particularly useful for future estimates of consumer diet composition.

Study Site
We worked at the Savannah River National Environmental Research Park, located in Aiken and Barnwell counties, South Carolina, USA (33 • 18 N, 81 • 37 W). This 780 km 2 site occurs in the Sandhill and upper Coastal Plain physiographic provinces. Forested areas cover >80% of the study site and are primarily comprised of managed stands of longleaf (Pinus palustris) and loblolly pine (P. taeda), with riparian forest and large patches of secondgrowth and hardwoods where fruit is seasonally abundant [39]. Detailed descriptions of the site, its history and management are provided by Kilgo and Blake [40].

Study Design
We collected feathers, blood, and fecal samples from passerines captured in the fall and winter of 1998 and 1999 and used two methods to evaluate dietary contributions to tissues: (1) a comparison of the isotopic values of tail feathers to those of plasma and blood, and (2) measuring the remains of arthropods, fruits, and seeds in fecal samples collected at the same time as the tissue samples. Flight feathers are grown during a~3 week interval [41] on the breeding grounds predominantly between late June and early September (hereafter "summer") for the species included in this study [42][43][44][45][46][47][48][49][50][51][52], although several species extend molt beyond this period and there may be considerable variation among individuals ( Figure 1). The isotopic value of feathers reflects an integration of dietary inputs during their growth, after which the tissue is inert [53,54]. In contrast, isotopic values of plasma and whole blood are constantly changing and represent the integration of dietary inputs during the hours-to-days (plasma) or weeks (blood) prior to sampling [19,25,55,56]. Gut passage rates depend on substrate, but fecal samples generally represent consumption no more than several hours prior to sample collection [22,57].
To parameterize our mixing model with isotopic signatures of the resources that birds were consuming locally, we collected arthropods, fruits, and seeds from the same site and at the same time of year as bird captures. We sampled fruit species known to be important to birds at the study site [39]. We collected invertebrates by sweeping vegetation with a net and by visually searching leaves and litter at randomly determined locations within all bird capture sites. We included only those arthropod taxa commonly contributing to the diets of our study species, identified via literature review or that we observed in feces of captured birds. We included only taxa of seeds that we observed birds eating on or near our study site. We did not include seeds found in fruit because the fruit-eating species we captured typically defecate those seeds intact (S. Pearson and D. Levey, pers. obs.).

Defining Fall and Winter
For seasonal comparisons, we used the change-point analysis approach and software developed by Taylor [58] to identify changes between a warmer than average to a colder than average ambient temperature, based on cumulative sum graphs of locally measured daytime high temperatures throughout study periods. For the 1998-1999 study period, the daily high temperature shifted from relatively mild fall conditions (average ± SD high temperature: 26.9 ± 4.6 • C) to cooler winter conditions (17.5 ± 6.4 • C) on 8 December 1998. For the 1999-2000 study period, the transition from fall conditions (27.0 ± 5.8 • C) to winter conditions (14.6 ± 6.1 • C) occurred on 12 November 1999. For subsequent analyses, we use these transition dates to define fall (before 8 December in 1998 and before 12 November in 1999) and winter (on or after 8 December in 1998, and on or after 12 November in 1999).

Bird Capture and Sample Collection
We captured a total of 653 birds of 47 species in mist nets placed in upland (n = 3) and bottomland (n = 1) areas of second growth and in bottomland hardwood stands (n = 3). Nets were active from 22 September 1998 to 19 February 1999 and from 24 August 1999 to 5 February 2000. Approximately 51% of these captures occurred in upland second growth, 24% in bottomland hardwood stands and 17% in bottomland second growth. Yellow-Rumped Warblers (Setophaga coronata) and Ruby-Crowned Kinglets (Regulus calendula) represented 24% and 15% of captures, respectively. We were able to obtain at least 5 blood samples from the following species, which were therefore included in subsequent analyses: residents = Eastern Towhee (Pipilio erythrophthalmus), Tufted Titmouse (Baeolophus bicolor), and Northern Cardinal (Cardinalis cardinalis); migrants = Gray Catbird (Dumetella carolinensis; rarely nests in the region and is common during migration) and Swainson's Thrush (Catharus ustulatus); overwintering migrants = American Robin (which nests in the region but not in our study area), Dark-Eyed Junco (Junco hyemalis), Hermit Thrush (Catharus guttatus), Ruby-Crowned Kinglet, White-Throated Sparrow (Zonotrichia albicollis), and Yellow-Rumped Warbler. Some species showed preferences for different habitats, with the majority of Yellow-Rumped Warblers (91%) and Ruby-Crowned Kinglets (70%) captured in upland second growth habitats and a majority of White-Throated Sparrows captured in bottomland hardwood stands (65%). However, none of our focal species changed habitat types between seasons, meaning that our subsequent analyses comparing seasonal diet shifts are unlikely to be driven by seasonal shifts in habitat use. Migratory movements by individuals have the potential to complicate stable isotope-based diet estimation by introducing the signatures of unmeasured, non-local diet items. However, we do not believe that this strongly influenced our diet estimates for two reasons: (1) for all but one species we used plasma, which has a rapid turnover rate (half-life no more than 24 h [19,25,55,56]); and (2) we found relatively high recapture rates during winter among resident (8-23%) and overwintering migrant (3-25%) species, with some species, such as the Hermit Thrush, known to hold overwintering territories in the region (e.g., [59]). Short turnover times likely result in local diet integration, while local residency implies consistent use of local resources. The one species for which we conducted stable isotope analysis on whole blood (Ruby-Crowned Kinglet, see below) was also a species with a high recapture rate and local winter residency, which may compensate for the integration of diet over a longer time.
Upon capture, birds were immediately placed in clean cloth bags. For all individuals, we plucked a center tail feather and, when possible, collected a blood sample from the brachial vein in a non-heparinized micro-hematocrit capillary tube. For most species, we were able to obtain a relatively large blood sample (~120 µL), which we then centrifuged to obtain the plasma for analysis. However, we were only able to obtain a small volume of blood from Ruby-Crowned Kinglets (~60 µL), and so analyzed samples from that species as whole blood. We transferred blood and plasma samples from capillary tubes to cryovials and stored them at −70 • C until analysis. Tail feathers were stored individually in plastic bags.

Fecal Analysis
When available, fecal samples were scraped from the inside of the cloth bags that held birds after capture and stored in 90% ethanol. We estimated percentages of fruit, invertebrates and seeds by placing each sample into a Petri dish and separating the sample into a relatively uniform layer of material. We then examined each sample under a dissecting microscope against a 2 × 2 mm grid background to count the number of grid cells obscured by arthropods, fruits, seeds, and other items in 4-5 randomly selected views, from which we derived an estimate of the percent of cells at least partially obscured by each dietary item. We identified many of the arthropod parts to order, following Moorman et al. [60].

Isotope Sample Preparation and Analysis
Prior to analysis in the winter of 2007/2008, we freeze-dried blood samples and prepared tail feathers by washing them in ether to remove oils [61] and cutting them into small pieces (~1 mm2). Fruit pulp (hereafter "fruit"), arthropods, and seeds were placed in a drying oven at 55 • C for seven days, and then ground into a powder using a mortar and pestle. We then loaded aliquots of samples into tin capsules for isotopic analysis.
The University of California, Davis Stable Isotope Facility performed all isotopic analyses. Samples were combusted in an elemental analyzer coupled to a continuous flow isotope ratio mass spectrometer (ANCA-CFIRMS; Europa Scientific Crewe, England; ANCA combustion unit and 20-20 mass spectrometer). The ratio of C and N stable isotopes in the CO 2 and N 2 gases produced is expressed in delta (δ) notation: [(Rsample/Rstandard) −1] × 1000, where δ is the isotope ratio of the sample relative to a standard (C = Vienna Peedee Belemnite (VPDB) limestone formation; N = atmospheric N (AIR)). Isotope ratios were derived by comparison of sample with reference gas injections in each analytical cycle. Raw delta values were normalized to international values δ 15 N AIR and δ 13 C VPDB using combusted standard samples (ammonium sulfate + sucrose; δ 13 C = −23.83‰, δ 15 N = +1.33‰) inserted at intervals in each batch of samples. Masses of C and N in each sample were derived by integrating the ion beam intensities and comparing the integrals to those of standard samples with known concentrations of C and N. Carbon and N mass were converted to percent by dividing the mass of N and C in the sample by the sample's mass and multiplying by 100.

Comparing Community-Level Diet Divergence between Summer, Fall, and Winter
We plotted the average tissue location in isotope space (defined by plotting mean δ 13 C against mean δ 15 N) to graphically assess the degree of bird community diet divergence in summer, fall and winter. Differences among species within and among seasons in δ 13 C identify the relative contribution to tissue synthesis from food resources originating from primary producers with different photosynthetic pathways. Higher values indicate greater contributions of C 4 pathways and lower values indicate greater contributions of C 3 pathways. Differences in δ 15 N identify the contribution to tissue synthesis from trophically distinct food resources or groups of food resources (e.g., resources from primary producers, low-trophic arthropods, and high-trophic arthropods). Values of δ 15 N generally increase with trophic level. We compared the dispersion of isotope values by calculating the average distance to the centroid for each season.
To distinguish high-trophic-level arthropods (e.g., spiders) from low-trophic arthropods (e.g., larval lepidoptera), we applied a k-means cluster analysis to the entire arthropod dataset, which identified two groups of arthropods separated at a δ 15 N of 3%; we classified those above and below this cutoff as high-and low-trophic arthropods, respectively. There was a gap in arthropod δ 13 C values from −21‰ to −18‰, so we arbitrarily used these values to distinguish between arthropods from C 3 and C 4 environments. We calculated δ 13 C values for C 3 and C 4 arthropods independent of trophic level. Likewise, we calculated δ 15 N values for high-and low-trophic arthropods independent of C 3 or C 4 status. For example, low-trophic C 4 arthropods were assigned the average δ 13 C value of all C 4 arthropods, regardless of trophic status. This method maximized the use of our arthropod samples and produced more consistent δ 13 C and δ 15 N values among dietary endpoints.
We used a similar approach to calculate seed δ 15 N endpoints. Our small sample of C 3 seeds resulted in an exceptionally low standard deviation in δ 13 C values, so we scaled the standard deviation to produce a coefficient of variation equivalent to that in C 4 seeds. For fruits and seeds, we weighted by sample size the contribution from a given taxon when calculating average δ 15 N and δ 13 C values, but used the uncertainty associated with the overall averages. For arthropods, each sample represented a different taxon (although often within the same order), so we elected to calculate unweighted averages for δ 15 N and δ 13 C.

Assessing Distinctiveness of Dietary Endpoints
Only isotopically distinct diet sources can be distinguished using mixing models. Consequently, we conducted a multivariate analysis of variance (MANOVA; PROC GLM; SAS Version 9.1, SAS Institute, Cary, NC, USA) to determine whether 15 N and 13 C varied among diet items. The MANOVA indicated that diets differed (see results). We then used Scheffe's multiple contrasts with follow-up univariate ANOVAs to establish how δ 15 N and δ 13 C differed among dietary items. Because PROC GLM cannot accommodate missing data, only diet items with measures of both δ 15 N and δ 13 C were used in the MANOVA analysis; we conducted follow-up ANOVAs and pair-wise contrasts in PROC MIXED (SAS Version 9.1, SAS Institute, Cary, NC, USA) using maximum likelihood estimation so that statistical comparisons are based on the full δ 15 N and δ 13 C datasets of dietary endpoints (i.e., those used in the mixing model analysis).

Dietary Discrimination Values
We used experimentally derived results in Pearson et al. [56] to estimate arthropod and fruit diet discrimination values (∆) relating diet-tissue discrimination values to the elemental concentration of the diet items (Table S2) [56]. These discrimination values are within the range reported by other authors (e.g., [25,61,62]) and avoid potential bias associated with assuming constant values [63]. For plasma and whole blood, we calculated ∆N and ∆C for fruit based on the average N and C concentrations of all fruit samples, and for arthropods based on the average N and C concentrations of all arthropod samples (Table 1). We calculated the uncertainty (standard deviations) associated with ∆N and ∆C by summing in quadrature the standard deviation of tissue-specific regression residuals from the published equations [56] and the standard deviation of discrimination estimates resulting from uncertainty in the elemental concentrations of the dietary endpoints. That is, the change in ∆N or ∆C when using N and C concentrations one standard deviation above or below their average represented the standard deviation of the discrimination value due to uncertainty in elemental concentrations. Similarly, we calculated diet-plasma δ 15 N and δ 13 C discrimination values and their uncertainty for the seed endpoint based on the average C and N concentrations of all seed samples (Table 1). Although not used in any mixing models, estimated ∆ values for whole blood from the Pearson et al. [56] equations agreed generally with published estimates of discrimination in blood cells of birds on known corn diets [64]. The tissue and diet specific discrimination values used in our mixing models are provided in Table S2. Table 1. Species-specific probabilities that dietary endpoints differ among seasons. Probabilities above 0.5 indicate support for the one-sided comparisons (e.g., winter fruit > fall fruit) whereas probabilities below 0.5 indicate support in the other direction (e.g., fall fruit > winter fruit). Dashes indicate that contributions of those dietary endpoints were not estimated for that species-season combination. Abbreviations: Arth = arthropod; H-C 3 = high-trophic C 3 arthropods; L-C 3 = low-trophic C 3 arthropods, etc. Winter samples were pooled across years for all species except the Yellow-Rumped Warbler, for which there was sufficient sample size to estimate the two years separately (winter 1 refers to 1998/1999 and winter 2 refers to 1999/2000).

Species
Comparison

Mixing Model Diet Estimation and Comparison of Fall and Winter Resource Use
We used concentration-dependent Bayesian mixing models to estimate the isotopic contributions of fruit, arthropods, and seeds to blood and plasma collected from birds in the fall and winter (sample sizes in Table S3). For this analysis, we used the SIAR program (version 4.2) [65], operated within the R statistical computing environment (version 3.1.2; R development core team, Vienna, Austria).
We related δ 15 N and δ 13 C values of plasma (or blood for Ruby-Crowned Kinglets) to δ 15 N and δ 13 C values of fruits, arthropods and seed species that plausibly occurred in the diet of each species, using N and C elemental concentrations from each food source and empirically derived discrimination factors from Yellow-Rumped Warblers (see above). The dietary endpoints used in mixing models (see above for numbers) varied among species based primarily on whether their typical diet included fruits or certain seeds: the potential diets of American Robin, Gray Catbird, Hermit Thrush, Ruby-Crowned Kinglet, Swainson's Thrush, and Yellow-Rumped Warbler included endpoints 1-3; Tufted Titmouse included endpoints 1-3 and 6; Eastern Towhee and Northern Cardinal included endpoints 1-7; and Dark-Eyed Juncos and White-Throated Sparrows included endpoints 2-7. We estimated average contributions of food items to songbird diets at the population level (≥5 individual samples) for each possible combination of species and season. Based on our experience in captive studies using known mixtures of diet items, this approach trades off the ability to estimate individual resource use for less variation and bias in estimates for species. The default settings for uniform priors, burn-in iterations, subsequent iterations, and thinning interval successfully resulted in probability distributions for estimated dietary proportions based on 10,000 iterations.
To estimate the probability that the contribution of a particular dietary endpoint exceeded that of another within or between seasons, we compared the mixing model output on an iteration-by-iteration basis (A. Parnell, pers. comm.). For example, to find the probability that the contribution of fruit in the winter diet of Ruby-Crowned Kinglets exceeded that of fruit in their fall diet, we divided the number of iterations in which the proportion of winter fruit was greater than the proportion of fall fruit (i.e., 8399) by the total number of iterations (i.e., 10,000), which in this case produced a probability of 0.8399. A probability of 0.5 indicates no difference in dietary contributions. Probabilities are symmetrical at approximately 0.5, so that a probability of 0.75 represents the same magnitude of difference as a probability of 0.25. In such a comparison, probabilities above 0.5 indicate support for the one-sided comparisons (e.g., winter fruit > fall fruit) whereas probabilities below 0.5 indicate support in the other direction (e.g., fall fruit > winter fruit). We restricted our quantitative comparisons between seasons to species with at least five samples in each season (Northern Cardinal, Ruby-Crowned Kinglet, White-Throated Sparrow, and Yellow-Rumped Warbler).

Seasonal Changes in Diet Dispersion
Feather (i.e., summer) δ 15 N and δ 13 C values were relatively clumped compared to more dispersed fall and winter plasma/blood δ 15 N and δ 13 C values, when plotted in isotope space (Figure 2). This increased dispersion was largely the result of a clear divergence in δ 13 C values of plasma/blood among the songbird community, where the range of 13 C values increased from 2.3‰ in the summer to 7.6‰ in the winter. The average distance to the centroid was 1.1 ± 0.8 (Mean ± SD) in the summer, 1.6 ± 0.6 in the fall, and 3.0 ± 0.7 in the winter. Specific changes in the winter included increased δ 13 C values in Dark-Eyed Junco, White-Throated Sparrow, and Northern Cardinal plasma, likely resulting from an increased contribution of C 4 dietary items. On average, δ 15 N values decreased in the fall and winter, suggesting lower trophic contributions to these tissues at that time of year ( Figure 2). However, this trend in δ 15 N was inconsistent among species with, for example, Ruby-Crowned Kinglet blood being relatively stable across time and seasons whereas that for Dark-Eyed Junco, White-Throated Sparrow, and Northern Cardinal plasma decreasing especially in winter (Figure 2).

Variability and Seasonal Patterns of Resource Use
For birds known to eat at least some fruit, results of the mixing model indicate that median fruit contribution to plasma/blood ranged from 4 to 26% in the fall (mean ± SD: 16.2 ± 7.5%; n = 7 species) and 9-37% in the winter (21.7 ± 10.3%; n = 5 species, Figure 3, Table S4). There was a similar pattern for seed-consuming species, with median seed contributions ranging from 28 to 32% in the fall (29.4 ± 2.6%; n = 3 species) and from 32 to 42% in the winter (36.6 ± 4.8%; n = 3 species, Figure 3, Table S4). Thus, (1) both fruits and seeds constitute substantial proportions of the fall and winter diets of this songbird community and (2) those contributions are generally larger and more variable in the winter than the fall, although this trend is more pronounced for seeds than fruit. These community-level patterns were mirrored by similar changes in the two overwintering migrants with adequate sample sizes for inter-seasonal comparisons. For both the Yellow-Rumped Warbler and Ruby-Crowned Kinglet, fruit resources contributed more to tissues during the winter than during the fall (Probabilities ≥ 0.77; Table 1). Specifically, median fruit contribution increased from 10.2% in the fall to 36.8% (1998/1999) or 19.2% (1999/2000) in the winter for Yellow-Rumped Warblers and median fruit contribution increased from 4.3% in the fall to 18.2% in the winter for Ruby-Crowned Kinglets (Figure 3, Table S4). Interestingly, this was not the case among resident Northern Cardinals, for whom the median fruit contribution to plasma decreased from 25.8% in the fall to 16.0% in the winter (Table 1, probability = 0.25). However, cardinals did use more C 4 seeds in the winter (17.8%) than in the fall (6.0%; Table 1, probability = 0.86). Meanwhile, the one overwintering migrant granivore with sufficient data for inter-seasonal comparisons, the White-Throated Sparrow, showed no difference between the fall and winter in median seed contribution (Table 1, probability = 0.57), but did increase their C 4 seed usage from 11.9% in the fall to 17.2% in the winter (Table 1, probability = 0.65).

Comparison of Fecal and Mixing Model Estimates of Songbird Fruit Consumption
Fecal samples from all our focal species except Dark-Eyed Junco and White-Throated Sparrow contained at least some fruit tissue. American Robin, Hermit Thrush, Gray Catbird, and Yellow-Rumped Warbler, all species previously described as being fruit-eaters during the fall and winter months, consistently contained at least 50% fruit pulp and skin in their feces (Figure 3). Fruit was also consistently consumed by Northern Cardinal and Tufted Titmouse (Figure 3), neither of which are noted frugivores. Mixing model estimates of fruit contributions to songbird diets showed a general positive relationship with fecal estimates of fruit consumption, although the average amount of fruit observed in the feces was consistently and substantially greater than mixing model estimates ( Figure 4). Moreover, both fecal and mixing model estimates of fruit contributions to songbird diets showed extreme variability regardless of sample size, with between-species coefficients of variation of 62.3% and 53.8%, respectively, and within-species coefficients of variation ranging from 7.8 to 100%. Details on the taxonomic identity of fruits found in songbird fecal samples are provided in the footnote to Table S1.

Fruit and Seed Contributions to Assimilated Songbird Diets
Consistent with prior data on fruit production [39], fruit disappearance [66,67], and songbird resource tracking [68][69][70] at our study site, results from our concentrationdependent mixing models indicated that fruit contributed substantially to songbird tissue synthesis during the fall and winter. For fruit-eating species at our site, fruit resources contributed as much as 26% to tissue synthesis in the fall (Northern Cardinal) and 37% in the winter (Yellow-Rumped Warbler). These higher estimates are generally greater than contributions estimated in other studies on related species that used similar stable isotope mixing models. For example, the greatest contribution of fruits among three species in the post-fledging period was only 13% in southern Indiana [31]. Similarly, the median contribution to fruit was only 14% in Gray-cheeked thrushes during spring migration in Colombia [30] and 15-22% in overwintering Wood Thrushes (Hylocichla mustelina) in Belize [71], although the estimated contribution increased to~35% during spring migration in the latter case. The lower estimates of fruit use by songbirds in these studies may be the result of differences in fruit availability between study sites, as the post-fledging period likely precedes peak fruit abundance in temperate areas and fruit resources are less pulsed in tropical areas [72]. Alternately, differences in nutritional requirements between seasons (e.g., higher energetic requirements during migration and overwintering) may explain the high fruit consumption present in our bird community during winter in a temperate area and high fruit consumption during migration in Belize [71]. However, the most likely explanation for the differences across studies is simply differences between focal species. Such interspecific differences in fruit consumption were readily apparent within our study; although we found large contributions of fruit for most bird species, the minimum estimated contribution was 4% in the fall (Ruby-Crowned Kinglet) and 9% in the winter (Hermit Thrush). This interspecific variation in contribution of fruits may be weak evidence for partitioning of fruit and other food resources among songbirds at our study site. For example, Eastern Towhee and Tufted Titmouse showed similar overall reliance on arthropods and seeds (Figure 3a), but Eastern Towhee showed markedly greater use of resources from C 4 plant communities ( Figure 3b).
As with fruits, mixing model estimates of songbird reliance on seeds matched our expectations based on prior research, with contributions as high as 32% in the fall (White-Throated Sparrow) and 42% in the winter (Dark-Eyed Junco). There are relatively few studies that have used mixing models to quantify seed use by songbirds, but seed use in our study was comparable to the reliance on tree mast by overwintering Rusty Blackbirds (Euphagus carolinus) in Georgia [73]. The relatively high seed contributions to tissues in our study were similarly comparable to the reliance of overwintering Black-capped Chickadees (Poecile atricapillus) on sunflower seed subsidies in Alaska [74], although, interestingly, the contribution of native seeds in that study was much lower. Overall, seed use appears to be much more consistent across studies and species than fruit use. However, despite this difference, our isotopic reconstructions of diet indicate that both fruits and seeds make substantial contributions to the tissues of individuals in this songbird community.

Comparing Estimates Derived from Stable Isotope Modeling to Those from Feces
Fecal analysis is a well-established method for analyzing the composition of animal diets and has been used to make estimates of fruit use by songbird populations [15,16]. However, few studies have contrasted fecal estimates of fruit use by songbirds with estimates derived using other methods such as stable isotope analysis. In our study, fecal data independently supported the importance of fruit to birds in the fall and winter. However, fecal data suggested a much greater contribution of fruit to songbird diets than did stable isotope data in all cases (Figures 3 and 4). Overall, fecal data estimated an average fruit contribution that was 214 ± 186% (mean ± SD) greater than values derived from mixing models.
There are several potential explanations for this consistent difference. For one, fecal data may overestimate the occurrence of fruit in the diet because birds were always caught within 2.6 m of the ground and usually in the vicinity of fruit-producing shrubs (e.g., Morella cerifera and Ilex opaca), which were widely scattered throughout our study site. As a result, we likely captured birds when they were near the ground foraging on fruit and not, for example, when in the canopy foraging on arthropods. Thus, considering the short time window of diet represented by fecal samples [22,57], estimates of reliance on fruit from feces may differ from estimates integrated over longer periods, with a bias towards fruit. Alternately, it is well established that songbirds are generally inefficient at digesting and absorbing fruits, particularly in comparison to arthropods [24,75,76], meaning that fruit tissue may be overrepresented in fecal samples. Mixing model estimates of resource use, by contrast, represent the contributions of assimilated nutrients and will therefore more accurately account for the contributions of arthropods that are more thoroughly digested and less apparent in feces. The underrepresentation of arthropods in fecal samples due to more efficient digestion is supported by the difficulty of determining the taxonomic identity of arthropod parts, which were highly fragmented and restricted to limited hard parts (e.g., mandibles, head capsules, elytra). Finally, isotopic routing could result in underor over-representation of dietary items to tissues because stable isotopes from various macronutrients in the diet may be routed differentially into tissues [62,77]. Specifically, if the carbon from dietary protein (primarily from arthropods) is preferentially routed to songbird plasma and whole blood, then mixing model estimates based on those tissues will tend to exaggerate the contribution of arthropods and undervalue the contributions of fruits. It is important to note that none of these explanations are exclusive-it is entirely possible that all of them contributed to the observed discrepancies between fecal and mixing model estimates of resource use. Considering both the scope of the disparity between these methods and the difficulty of identifying the source of the difference, it is clear that stable isotope mixing model and fecal estimates of diet provide very different results and, due to assimilation and processing differences, are not readily substitutable for one another. Accordingly, the choice of method will depend on the context of the study, and studies focused on organism-level patterns and consequences of frugivory, where assimilated nutrients are particularly important, should ideally use stable isotope analysis and interpret fecal estimates of diet cautiously.

Seasonal Changes in Diet: The Importance of Fruits and Seeds
Seasonal diet shifts have been documented in species ranging from water snakes [78] to geese [79], with many songbirds in the northern hemisphere exhibiting a seasonal change in fruit consumption in the fall as fruits become more available and arthropods less available [12][13][14]. Seasonal variation in seed consumption by songbirds has not been quantified to the same extent as frugivory, but there is evidence for seasonal changes in seed production and abundance that could drive similar shifts [80][81][82]. Consistent with our expectations for both fruits and seeds, estimated contributions of these resources to songbird tissues were generally higher in the winter than the fall across our focal community. These community-level patterns were mostly mirrored in species for which we were able to make population-level seasonal comparisons: two out of three species used greater proportions of fruits in the winter than fall and both seed-consuming species used greater proportions of seeds in the winter than fall. Thus, our results largely corroborate previous findings about broad-scale temporal patterns of resource use by songbirds.
Our results suggest interesting seasonal trends in the variability of resource use across this songbird community. At the community level (Figure 2), songbird diets showed a marked increase in dispersion in isotope space from the summer to the fall to the winter, a finding supported by the greater range of fruit and seed contributions estimated by mixing models in the winter than the fall. These results indicate that, although this community generally shifted away from a reliance on arthropod resources in the fall and winter, that shift was not consistent in direction or magnitude. Such divergence in diets is consistent with winter food limitation: as food becomes scarcer or as consumer abundance increases, resource use should diverge to reduce the potential for competition [83][84][85]. Our results provide evidence that some of this divergence comes at the level of species (Table 1, Figure 3). For example, both Yellow-Rumped Warbler and Ruby-Crowned Kinglet increased their reliance on fruit in the winter compared to the fall, whereas Northern Cardinal decreased fruit usage in the winter. Meanwhile, both Northern Cardinal and White-Throated Sparrow increased their reliance on C 4 seeds in the winter compared to the fall. This partitioning may simply result from species-specific differences in foraging and digestive capabilities [4,18], although it is interesting that frugivory increased in the two overwintering migrants and decreased in the one resident. For migrating birds, food limitation during winter can negatively affect body condition and reduce subsequent breeding success [79,86,87], and so there may be more incentive to take advantage of predictably available fruit resources that minimize foraging time and effort [39,88].
In addition to between-species variation, we also found evidence for within-species variation in diet shifts, as levels of frugivory by Yellow-Rumped Warblers differed between the two winters in which we collected plasma samples. It is possible that this disparity reflects interannual variation in fruit availability, which has been documented at our study site [39,68]. However, variation in resource use at the community and population levels may also be the result of distinct preferences among individuals, which have been observed to vary seasonally in some cases [89,90]. Thus, the high degree of variation in our estimates of diet, as well as variation among years may indicate the influence of individual preferences and possibly resource partitioning among individuals. In either case, our results demonstrate a great deal of flexibility in resource use by our focal songbird community, which add to the complexity of describing seasonal diet shifts.

Conclusions
Diet is a key interface between the organism and its environment, and so quantifying diet and dietary shifts can provide valuable insight into animal-environment relationships. Here, we used a stable isotope analysis of tissues and a range of diet items to document a notable reliance on non-arthropod resources by a mixed migrant and resident songbird community during the non-breeding season. Not only do these estimates add to the relatively few studies that quantify the contributions of fruit and seeds to songbird tissues, they also provide a useful step towards identifying critical winter food resources, which has important implications for conservation and management of migratory birds. Furthermore, by characterizing the variability of resource use across this songbird community and in different seasons, our results suggest that there is considerable separation of diet items both between species and within species. Finally, our study demonstrates that, although there may be considerable uncertainty associated with estimates of diet composition, there is a consistent and substantial offset between estimates derived from stable isotope mixing models and those derived from fecal analysis. This finding highlights an important consideration for those wishing to compare historical accounts of songbird diets as well as those designing dietary studies.

Data Availability Statement:
A database containing measurements of the stable isotope values of diet items and songbird tissues as well as the composition of fecal samples is archived in the Dryad data repository and available at the following doi: https://doi.org/10.5061/dryad.xwdbrv1dc.