Title : Temporal changes in the individual size distribution decouple long-term trends in 1 abundance, biomass, and energy use of North American breeding bird communities

A core objective of contemporary biodiversity science is to understand long-term trends in the structure and function of ecological communities. Different currencies of ecological function – specifically, total abundance, total standing biomass, and total metabolic flux – are naturally linked, but may become decoupled if the underlying size structure of a system changes. Here, we seek to establish how changes in community size composition modulate long-term relationships between different currencies of ecological function for North American birds. We used species’ traits and allometric scaling to estimate individual size measurements and basal 41 metabolic rate for birds observed in the North American Breeding Bird Survey. We compared 42 the long-term trajectories for community-wide standing biomass and energy use to the long-term 43 trends driven by changes in individual abundance alone. Finally, we used dissimilarity metrics to 44 evaluate the link between changes in species and size composition and changes in the 45 relationship between abundance- and size-driven dynamics.


Introduction 63
Understanding the interrelated dynamics of size-and abundance-based dimensions of 64 biodiversity is key to understanding biodiversity change in the Anthropocene. Total abundance -65 i.e. the total number of individual organisms present in a system -and size-based currencies -66 such as the total biomass or total metabolic flux ("energy use") of a system -are intertwined, but 67 not necessarily equivalent, measures of biological function. Abundance is more closely tied to 68 species-level population dynamics, while size-based metrics more directly reflect assemblage- Specifically, at the community scale, changes in the relationship between size and abundance can 82 signal important shifts in community structure and functional composition. To the extent that 83 size is a proxy for other functional traits, the dynamics of the community-level size structure 84 (individual size distribution, ISD) over time may reflect processes related to niche structure 85 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint consistency in the size structure may maintain the relationship between size-and -abundance 92 based currencies, even as species composition, total abundance, and total biomass and total 93 energy use fluctuate over time, which can reflect consistency in the niche structure over time 94 (Holling 1992). 95 It is important to improve our understanding of these dynamics for terrestrial animal 96 communities in particular. In contrast to terrestrial trees and aquatic systems (Kerr and Dickie Macroecological-scale synthesis on the interrelated dynamics of the ISD, total abundance, and 107 community function for terrestrial animals has been constrained by 1) a lack of community-level 108 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made either 1) how the size structures for these systems behave over time or 2) the extent to which 122 changes in community size structure decouple the community-level dynamics of abundance, 123 biomass, and energy use in these systems. 124 Here, we begin to address this gap by exploring how temporal changes in the size structure 125 modulate the relationship between total abundance, energy, and biomass for communities of 126 North American breeding birds. We used allometric scaling to estimate community size and 127 abundance data for the North American Breeding Bird Survey, and evaluated how changes in 128 total abundance, biomass, and energy use have co-varied from 1988-2018. Specifically, we 129 examined: 1) How often do these currencies change together vs. exhibit decoupled dynamics?; 2) 130 What are the dominant directions and magnitudes of the overall change over time and degree of 131 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made routes that had been sampled in at least 27 of 30 years in this window (n = 739) and compared 150 these results to a more strict selection of routes that were sampled in every year (n = 199). 151 Results for this more stringent subset of routes were qualitatively the same as for the more 152 inclusive selection of routes (Appendix S1). We take the route to be the "community" scale 153 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint distributions at a spatial and temporal scale that would otherwise be impossible (Thibault et al. 177

2011). 178
Comparing abundance-and size-based currencies 179 Comparing trends across different currencies is a nontrivial statistical problem. Because different 180 currencies vary widely in their scale of measure (e.g. abundance in the hundreds of individuals; 181 total biomass in the thousands of grams), it is challenging to interpret differences in magnitude of 182 slope across different currencies. Transformation and scaling using common approaches (such as 183 a square-root transformation, or rescaling each currency to a mean of 0 and a standard deviation 184 of 1) destroys information about the degree of variability within each currency that is necessary 185 in order to make comparisons between currencies for the same timeseries. 186 Therefore, rather than attempting to compare slopes across currencies or to transform different 187 currencies to a common scale, we used a simple null model to generate the expected dynamics in 188 biomass and energy use if the individual size distribution had remained constant over time, but 189 allowed abundance to vary consistent with observed dynamics. In effect, we generated the 190 expected dynamics of biomass and energy use if only abundance drove changes in those 191 currencies over time. For each route, we characterized the "observed" timeseries of total biomass 192 and total energy use by simulating size measurements for all individuals observed in each time 193 step and summing across individuals, using the method described above. We then simulated 194 timeseries for "abundance-driven" dynamics of biomass and energy use incorporating observed 195 changes in community-wide abundance over time, but under a scenario of consistent species (and 196 therefore approximate size) composition over time. For each community, we characterized the 197 timeseries-wide probability of an individual drawn at random from the community belonging to 198 a particular species (ܲሺ‫ݏ‬ ሻ ) as each species' mean relative abundance taken across all timesteps: 199 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint is the abundance of species

Long-term trends 210
For each route, we evaluated the "observed" 30-year trend in biomass (or energy use) and 211 compared this to the trend derived from the "abundance-driven" null model using generalized 212 linear models with a Gamma family and log link (appropriate for strictly-positive response 213 variables such as biomass or total energy use). We fit four model formulas to characterize 1) the 214 trend in biomass (or energy use) over time and 2) whether this trend deviates from the trend 215 expected given only changes in individual abundance. These models correspond to qualitatively 216 different "syndromes" of change: 217 1. biomass ~ year * dynamics or energy use ~ year * dynamics in which "dynamics" refers to 218 being either the "observed" or "abundance-driven" (null model) dynamics. This model fits a 219 slope and intercept for the observed trend in biomass or energy use over time, and a 220 separate slope and intercept for the trend drawn from the abundance-driven, or null model, 221 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint dynamics. We refer to this model as describing a syndrome of "Decoupled trends" between 222 abundance-driven and observed dynamics. 223 2. biomass ~ year + dynamics or energy use ~ year + dynamics. This model fits a separate 224 intercept, but not slope, for the abundance-driven and observed dynamics. This model was 225 never selected as the best-performing description of community dynamics. 226 3. biomass ~ year or energy use ~ year. This model fits a temporal trend, but does not fit 227 separate trends for the observed and abundance-driven dynamics. We refer to this syndrome 228 as "Coupled trends" between abundance-driven and observed dynamics. (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint Relating change in community structure to decoupling between abundance and size-based 241 dynamics 242 We used dissimilarity metrics to explore the extent to which change in community species or 243 size composition caused decoupling between long-term trends in individual abundance and total 244 biomass and energy use. These dissimlilarity metrics are most readily interpretable when making 245 pairwise comparisons (as opposed to repeated comparisons over a timeseries). We therefore 246 made comparisons between the first and last five-year intervals in each timeseries, resulting in a 247 We calculated three metrics to explore how changes in community composition and size 253 structure translate into decoupling between abundance-driven and observed dynamics for 254 biomass and energy use. First, we evaluated the change in average community-wide body size, 255 calculated as the absolute log ratio of mean body size in the last five years relative to the mean 256 body size in the first five years: 257 are the mean body size of all individuals observed in the first and last 5 258 years, respectively. Large changes in average body size are, by definition, expected to translate 259 into decoupling between observed and abundance-driven dynamics. 260 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint Second, we calculated measures of turnover in the size structure and in species composition. We 261 calculated turnover in the ISD using a measure inspired by an overlap measure that has 262 previously been applied to species body size distributions in mammalian communities (Read et 263 al. 2018). We characterized each "begin" or "end" ISD as a smooth probability density function 264 is the probability density from the density smooth for the 275 second ISD at that evaluation point. We subtracted this quantity from 1 to obtain a measure of 276 turnover between two ISDs. 277 To evaluate turnover in species composition between the five-year time periods, we calculated 278 Bray-Curtis dissimilarity between the two communities using the R package vegan (Pinheiro et 279 al. 2020). 280 We tested whether routes whose dynamics were best-described by each "syndrome" of change -281 i.e. "Decoupled trends", "Coupled trends", or "No directional change" -differed in 1) the 282 magnitude of change in mean body size; 2) turnover in the ISD over time; or 3) species 283 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint compositional turnover (Bray-Curtis dissimilarity) over time. For change in mean body size, we 284 fit an ordinary linear model of the form absolute_log_ratio ~ syndrome. We used the absolute 285 log ratio so as to focus on the magnitude, rather than the direction, of change in body size (see 286 also Supp and Ernest (2014) for the use of the absolute log ratio to examine the magnitudes of 287 differences between values). We compared this model to an intercept-only null model of the 288 form abs(log_ratio) ~ 1. Because our metrics for turnover in the ISD and species composition 289 are bounded from 0-1, we analyzed these metrics using binomial generalized linear models of the 290 form ISD_turnover ~ syndrome and Bray_Curtis_dissimilarity ~ syndrome, and again compared 291 these models to intercept-only null models. In instances where the model fit with a term for 292 syndrome outperformed the intercept-only model, we calculated model estimates and contrasts 293 using the R package emmeans (Lenth 2021). 294

295
Of the 739 routes in this analysis, approximately 70% (501/739 for biomass, and 509/739 for 296 energy use) exhibited a significant temporal trend in either abundance or in biomass/energy use 297 that resulted in the route being classified as exhibiting either "Decoupled trends" or "Coupled 298 trends" (Table 1). All results were qualitatively the same using a subset of 199 routes with 299 complete temporal sampling over time (Appendix S1). Trends driven by individual abundance, 300 as reflected by the dynamics of a simple null model, were strongly dominated by declines (67% 301 declines and 33% increases for abundance-driven dynamics in biomass, and 70% decreases and 302 30% increases for abundance-driven dynamics in energy use; Figure 2; Table 2). However, for 303 biomass, the long-term temporal trends were evenly balanced between increases and decreases 304 (49% increasing and 51% decreasing; Figure 2; Table 2). For energy use, there was a greater 305 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint representation of decreasing trends than for biomass, but still less so than for strictly abundance-306 driven dynamics (65% decreasing and 35% increasing trends; Figure 2; Table 2). 307 These divergent aggregate outcomes in individual abundance, energy use, and especially 308 biomass occurred due to decoupling in the long-term trends for these different currencies. For a 309 substantial minority of routes (20% of all routes for biomass, and 7% of all routes for energy 310 use), long-term dynamics were best-described as a syndrome of "Decoupled trends" (that is, with 311 a different slope for biomass or energy use-driven dynamics than for the "null", individual 312 abundance-driven, trend) (Table 1). When this decoupling occurred, it was dominated by 313 instances in which the slope for individual abundance-driven dynamics was more negative than 314 that for biomass or energy use (Figure 3). 315 Decoupling between the long-term trajectories of individual abundance and energy use or 316 biomass is, by definition, indicative of some degree of change in the ISD over time. Routes 317 whose dynamics for biomass were best-described as syndromes of decoupled trends over time 318 had a higher absolute log ratio of mean mass (i.e. greater magnitude of change, either increasing 319 or decreasing, in mean mass over time) than routes with coupled or no directional trends ( Figure  320 4; Appendix S2 Tables S1-S3). However, there was not a detectable difference in the degree of 321 temporal turnover in the ISD overall (Figure 4; Appendix S2 Table S4), or in species 322 composition (Figure 4; Appendix S2 Table S5), compared between routes that exhibited different 323 syndromes of change. 324

Discussion 325
Abundance, biomass, and energy use are nonequivalent currencies 326 Simultaneously examining multiple currencies of community-level abundance revealed 327 qualitatively different continent-wide patterns in the long-term trends for abundance in terms of 328 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint individuals, biomass, and energy use. While long-term trends in individual abundance were 329 dominated by decreases, long-term trends in biomass were evenly split between increases and 330 decreases, and trends in energy use were again dominated by declines (Figure 2) biomass, and energy use naturally co-vary to some extent, shifts in the community size structure 343 can and do produce qualitatively different trends for these different currencies. These may reflect 344 contrasting long-term changes in different types of community processes -for example, shifts in 345 habitat structure that affect the optimal body sizes for organisms in a system, but do not result in 346 overall changes in resource availability (e.g. White et al. (2004)). Moreover, extrapolating the 347 long-term trend from one currency to another may elide underlying changes in the community 348 that complicate these dynamics. To appropriately monitor different dimensions of biodiversity 349 change, it is therefore important to focus on the specific currency most closely aligned with the 350 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. these findings reflect turnover in species composition broadly favoring larger-bodied species. A 364 clear next step for this work is to identify the proximate and ultimate drivers of these shifts -for 365 example, by identifying which groups or species are responsible for these trends, how these 366 shifts vary over regions or habitat types, and if and how they are linked to underlying changes in 367 habitat quality or anthropogenic disturbances. An equally important counterpoint to this work 368 will be integrating potential shifts in intraspecific body size, and particularly declines in body 369 size associated with rising temperatures, with the interspecific dynamics documented here 370 (Youngflesh et al. 2022). In principle, declining body size over time could partially offset the 371 interspecific size shifts observed here. Given that interspecific size shifts occur over a 372 dramatically larger range of body sizes than do intraspecific size shifts, we anticipate that, in 373 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint most instances, the community-wide dynamics would remain qualitatively despite intraspecific 374 size change. While there are not currently robust estimates of intraspecific size changes at the 375 scale of the full Breeding Bird Survey dataset, this is an excellent opportunity for a focused study 376 in a well-documented system to elucidate the net effects of inter and intraspecific size change on 377 community-level properties. Finally, we note that these increases in body size do not generally 378 appear great enough to decouple the long-term trends in energy use from total abundance (Figure  379 3). Energy use scales nonlinearly with body size with an exponent less than 1, which means that 380 community-wide increases in mean body size result in smaller increases in total energy use than 381 in total biomass. 382

Complex relationships between compositional change and community-level properties 383
The decoupling between the long-term trends for biomass, abundance, and energy use 384 demonstrated in many of the communities studied here is symptomatic of a directional shift in 385 the size structure -in these instances, generally favoring larger bodied species. However, 386 examining the community-wide dynamics of turnover in species composition and the overall size 387 structure reveals that the relationship between changes in community structure and changes in 388 the scaling between different currencies of community-wide abundance is considerably more 389 nuanced than simple directional shifts in mean size. Routes that exhibit a statistically detectable 390 decoupling between total biomass and total abundance show large changes in average body size 391 compared to routes for which biomass and abundance either change more nearly in concert with 392 each other or do not show temporal trends (Figure 4; Appendix 2 Tables S1-S3). This aligns 393 naturally with mathematical intuition given the intrinsic relationship between average body size, 394 total abundance, and total biomass. However, these routes are not extraordinary in terms of their 395 overall degree of temporal turnover in either the size structure or in species composition. Rather, 396 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint the levels of turnover in overall community structure are comparable between routes that show 397 decoupling between abundance and biomass, statistically indistinguishable trends, or no temporal 398 trends in either currency (Figure 4; Appendix 2 Tables S4-S5).  399 For many communities, therefore, there has been appreciable change in the species and size 400 composition that simply does not manifest in a shift in the overall community-wide mean body 401 size or mean metabolic rate sufficient to decouple the dynamics of biomass, abundance, and (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint across North America, we find that changes in species and size composition produce 420 qualitatively different aggregate patterns in the long-term trends of abundance, biomass, and 421 energy use, highlighting the nuanced relationship between these related, but decidedly 422 nonequivalent, currencies and reflecting widespread changes in community size structure that 423 may signal substantive changes in functional composition. Simultaneously, the complex 424 relationship between turnover in community species and size composition, and the scaling 425 between different currencies of community-level abundance, highlights opportunities for 426 synergies between recent computational and statistical advances, case studies grounded in 427 empiricism and natural history, and future macroecological-scale synthesis to realize the full 428 potential of this conceptual space. 429 430 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint Data availability: All data and code supporting this manuscript are available online on GitHub. 535 For the purposes of double-blind review, we have uploaded a copy of these analyses to a 536 temporary GitHub repository at https://github.com/bbssizeshifts/BBSsims. Upon manuscript 537 acceptance, these will be archived in perpetuity on Zenodo. 538 539 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint Tables  540   Table 1 31-32% of routes are best described as syndromes of "No directional change" (intercept-only 549 models). For the remaining routes, in most instances, the dynamics of biomass and energy use 550 exhibit a temporal trend, but with no detectable difference in the temporal trends for abundance-551 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint driven and observed dynamics ("Coupled trends"). However, for a substantial minority of routes 552 (20% overall for biomass, or 30% of routes with a temporal trend; 7% overall for energy use, or 553 10% of routes with a temporal trend), there is a detectable deviation between the trends expected 554 due only to changes in abundance and the observed dynamics ("Decoupled trends"). 555 556 557 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint The proportion of trends that are increasing (specifically, for which the ratio of the last fitted 560 value to the first fitted value > 1) for abundance-driven and observed dynamics, for routes 561 exhibiting temporal trends ("coupled trends" or "decoupled trends") in total biomass and total 562 energy use (for biomass, n = 501; for energy use, n = 509). Trends that are not increasing are 563 decreasing. 564 Trends in abundance-driven dynamics are dominated by declines (67% of routes for total 565 biomass, and 70% of routes for total energy). Observed dynamics for biomass differ qualitatively 566 from the abundance-driven dynamics. Specifically, observed trends in biomass are evenly 567 divided between increases and decreases (49% increasing). Observed trends in energy use more 568 closely mirror abundance-driven trends (65% declines). 569 570 571 . CC-BY 4.0 International license available under a (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made The copyright holder for this preprint this version posted November 9, 2022. ; https://doi.org/10.1101/2022.11.08.515659 doi: bioRxiv preprint decoupled the dynamics of biomass from those that would occur due only to changes in 581 abundance. The slope for abundance-driven dynamics is significantly more negative than for the 582 observed dynamics (interaction term p = 0.0013). B. Underlying changes in the ISD. The 583 individual size distributions for the first 5 years (solid lines) and last 5 years (dashed lines) of the 584 timeseries. The x-axis is body size (as mass in grams; note log scale) and the y-axis is probability 585 density from a Gaussian mixture model fit to a vector of simulated individual masses for all 586 individuals observed in the years in questions, standardized to sum to 1. For the abundance-587 driven (blue) scenario, individuals' species identities (which determine their body size estimates) 588 are re-assigned at random weighted by each species' mean relative abundance throughout the 589 timeseries, resulting in a consistent individual size distribution over time. For the observed (gold) 590 scenario, individuals' body sizes are estimated based actual species abundances at each time 591 step. For this route, species composition has shifted over time and produced different ISDs for 592 the "begin" and "end" time periods. Specifically, the "end" ISD has peaks at larger body sizes 593 (ca. 90g and 500g) not present in the "begin" ISD. This redistribution of density towards larger 594