Predicting abundance–body size relationships in functional and taxonomic subsets of food webs

Maximum likelihood outperforms binning methods for detecting differences in abundance size spectra across environmental gradients

Individual body size distributions (ISD) within communities are remarkably consistent across habitats and spatiotemporal scales and can be represented by size spectra, which are described by a power law. The focus of size spectra analysis is to estimate the exponent (λ ) of the power law. A common application of size spectra studies is to detect anthropogenic pressures. Many methods have been proposed for estimatingλ most of which involve binning the data, counting the abundance within bins, and then fitting an ordinary least squares regression in log-log space. However, recent work has shown that binning procedures return biased estimates ofλ compared to procedures that directly estimateλ using maximum likelihood estimation (MLE). While it is clear that MLE produces less biased estimates of site-specific λ's, it is less clear how this bias affects the ability to test for changes in λ across space and time, a common question in the ecological literature. Here, we used simulation to compare the ability of two normalised binning methods (equal logarithmic and log2 bins) and MLE to (1) recapture known values ofλ , and (2) recapture parameters in a linear regression measuring the change inλ across a hypothetical environmental gradient. We also compared the methods using two previously published body size datasets across a natural temperature gradient and an anthropogenic pollution gradient. Maximum likelihood methods always performed better than common binning methods, which demonstrated consistent bias depending on the simulated values ofλ . This bias carried over to the regressions, which were more accurate whenλ was estimated using MLE compared to the binning procedures. Additionally, the variance in estimates using MLE methods is markedly reduced when compared to binning methods. The error induced by binning methods can be of similar magnitudes as the variation previously published in experimental and observational studies, bringing into question the effect sizes of previously published results. However, while the methods produced different regression slope estimates, they were in qualitative agreement on the sign of those slopes (i.e. all negative or all positive). Our results provide further support for the direct estimation ofλ and its relative variation across environmental gradients using MLE over the more common methods of binning.

Coexistence, Energy, and Trophic Cascade in a Three-Level Food Chain Integrating Body Sizes

Predation is a biological interaction that influences demographic patterns by modifying community structure. In the current ecological crisis, there is a need to better understand the conditions of coexistence between predators, prey and their resources. The body size is considered a key feature to explain community-scale phenomena, energetic, and evolutionary constraints. This raises the question of how species body size directly or indirectly affects the demographic patterns that enable coexistence. Considering the above, we conducted a theoretical study that implements a Rosenzweig-MacArthur type model, which represents a three-level chain that integrates body sizes and includes a Holling type I functional response. In this model, we characterize coexistence through body size-dependent net reproductive rates. Our results suggest that the body sizes of consumer species strongly affect the size-density relations and energy requirements. We obtain the negative relationship between body size and density of intermediate consumers and discuss the energy equivalence rule. Furthermore, larger predators have a more significant impact on the intensity of the trophic cascade than smaller predators. Finally, we discuss potential extensions and applications of our modeling approach.

Seasonal and ontogenetic variation of whiting diet in the Eastern English Channel and the Southern North Sea

An accurate description of trophic interactions is crucial to understand ecosystem functioning and sustainably manage marine ecosystems exploitation. Carbon and nitrogen stable isotopes were coupled with stomach content analyses to investigate whiting (Merlangius merlangus, Linnaeus, 1758) feeding behavior in the Eastern English Channel and Southern North Sea. Whiting juveniles and adults were sampled in autumn and winter to investigate both ontogenetic and seasonal changes. In addition, queen scallops (Aequipecten opercularis) samples were collected along with fish to be used as isotopic benthic baseline. Results indicated an ontogenetic diet change from crustaceans to fish and cephalopods. In autumn, δ¹⁵N values generally increased with fish size while in winter, a decrease of δ¹⁵N values with fish size was observed, as a potential result of spatial variation in baseline δ¹⁵N values. In winter, a nutrient-poor period, an increase in feeding intensity was observed, especially on the copepod Temora longicornis. This study provides further insights into whiting trophic ecology in relation to ontogenetic and seasonal variations, and it confirms the importance of combining several trophic analysis methods to understand ecosystem functioning.

Disturbance regulates the density-body-mass relationship of soil fauna

Theory on the density-body-mass (DBM) relationship predicts that the density of animal species decreases by the power of -0.75 per unit increase in their body mass, or by the power of -1 when taxa across trophic levels are studied. This relationship is, however, largely debated, as the slope often deviates from the theoretical predictions. Here, we tested the ability of the DBM relationship to reflect changes in the structure of communities subjected to an anthropogenic disturbance. The slope would become less steep if smaller animals were more impacted by the disturbance than the larger ones, whereas the slope would become steeper if larger animals were more affected than the smaller ones. We tested the changes in the DBM relationship by sampling soil fauna, i.e., nematodes, Collembola, and larger arthropods, from a semiarid grassland before and after spraying diesel fuel as disturbance. We applied three different treatments: a control, a light disturbance, and an intense disturbance. We found that the slopes of the DBM relationships before the disturbance were around -1 as predicted by theory. The slope became more positive (i.e., less steep) just after the disturbance, especially after the intense disturbance as smaller fauna suffered the most and early colonizers had larger body mass. Interestingly, we observed that the slopes converged back to -1 by 2 months post-disturbance. Our findings show that the response of soil fauna communities to anthropogenic disturbances could explain the large variation in observed slopes of the DBM relationships. We experimentally demonstrate that an animal community, when disturbed, shows a temporal pattern of DBM relationships ranging from deviations from the predicted slope to convergence to the predicted slope with time. We recommend that deviations in the DBM relationships after disturbances can provide insights in the trajectory of community recovery, and hence could be used for biomonitoring.

Shift in a large river fish assemblage: body-size and trophic structure dynamics

As the intensity and speed of environmental change increase at both local and global scales it is imperative that we gain a better understanding of the ecological implications of community shifts. While there has been substantial progress toward understanding the drivers and subsequent responses of community change (e.g. lake trophic state), the ecological impacts of food web changes are far less understood. We analyzed Wabash River fish assemblage data collected from 1974-2008, to evaluate temporal variation in body-size structure and functional group composition. Two parameters derived from annual community size-spectra were our major response variables: (1) the regression slope is an index of ecological efficiency and predator-prey biomass ratios, and (2) spectral elevation (regression midpoint height) is a proxy for food web capacity. We detected a large assemblage shift, over at least a seven year period, defined by dramatic changes in abundance (measured as catch-per-unit-effort) of the dominant functional feeding groups among two time periods; from an assemblage dominated by planktivore-omnivores to benthic invertivores. There was a concurrent increase in ecological efficiency (slopes increased over time) following the shift associated with an increase in large-bodied low trophic level fish. Food web capacity remained relatively stable with no clear temporal trends. Thus, increased ecological efficiency occurred simultaneous to a compensatory response that shifted biomass among functional feeding groups.

Warming alters community size structure and ecosystem functioning

Global warming can affect all levels of biological complexity, though we currently understand least about its potential impact on communities and ecosystems. At the ecosystem level, warming has the capacity to alter the structure of communities and the rates of key ecosystem processes they mediate. Here we assessed the effects of a 4°C rise in temperature on the size structure and taxonomic composition of benthic communities in aquatic mesocosms, and the rates of detrital decomposition they mediated. Warming had no effect on biodiversity, but altered community size structure in two ways. In spring, warmer systems exhibited steeper size spectra driven by declines in total community biomass and the proportion of large organisms. By contrast, in autumn, warmer systems had shallower size spectra driven by elevated total community biomass and a greater proportion of large organisms. Community-level shifts were mirrored by changes in decomposition rates. Temperature-corrected microbial and macrofaunal decomposition rates reflected the shifts in community structure and were strongly correlated with biomass across mesocosms. Our study demonstrates that the 4°C rise in temperature expected by the end of the century has the potential to alter the structure and functioning of aquatic ecosystems profoundly, as well as the intimate linkages between these levels of ecological organization.

Application of nitrogen stable isotope analysis in size-based marine food web and macroecological research

Interacting human and environmental pressures influence the structure and dynamics of marine food webs. To describe and predict the effects of these pressures, theoretical advances need to be supported by a capacity to validate the underlying models and assumptions. Here, we review recent applications of nitrogen stable isotope analysis in marine food web and macroecological research, with a focus on work that has paralleled a resurgence of interest in the development and application of size-based models. Nitrogen stable isotope data have been used to estimate intra- and inter-specific variation in trophic level, predator-prey size ratios, transfer efficiency, food chain length, relationships between predator and prey species diversity and the dynamics of energy use. Many of these estimates have contributed to the development, testing and parameterisation of food web and ecosystem models, some of which have been used to establish baselines for assessing the scale of human impacts. The interpretation of results depends on assumed fractionation but, when supported by sensitivity analyses and experimental validation, nitrogen stable isotope data provide valuable insights into the structuring of marine communities and ecosystems.

Predicting the abundance of European stream macroinvertebrates using biological attributes

Is there a relationship between the abundance of organisms and particular biological attributes? To assess this old, yet still acutely debated key question of ecology, we have used large databases on 312 stream macroinvertebrate genera (from 27 orders) that describe (1) invertebrate abundance at 527 least human-impacted European stream sites, (2) 11 biological traits (size, life-history, food, among others) described in 61 biological trait categories (BTCs; e.g. small, intermediate or large size) and (3) 14 attributes indicating specialization (AISs; e.g. species richness, size and food diversity). We applied interactive procedures to obtain models (for BTCs, AISs and a mixture of both descriptions) explaining as much as possible of the abundance variability of the genera with the lowest number of significant and ecologically meaningful attributes and assessed the predictive power of these models (in crosswise validations) by comparing predicted and observed abundances. Mean European invertebrate abundance increased with BTC affinities favouring viability in stream systems (e.g. attachment to the stream bottom to resist the flow, aquatic passive dispersal with the flow, exploitation of abundant food sources) and decreased with BTC affinities disfavouring this viability (e.g. drag force increase associated with larger body size, flow exposure associated with aerial respiration). Abundance consistently decreased with specialization of the genera (e.g. low species richness, oddity of their overall BTC profile from an "average" European genus). The model including a mixture of a few BTCs and AISs had the greatest predictive power: it predicted 35% of the observed abundance (ln-transformed) variability of the genera; these predictions were marginally affected by taxonomy (using orders as categorical variables). We conclude that a better appreciation of the influence of the examined taxonomic diversity, number and type of biological attributes, environmental system and spatial scale could enable abundance predictions using different sets of biological attributes for different taxonomic groups and systems.

Measurement of body size and abundance in tests of macroecological and food web theory

1. Mean body mass (W) and mean numerical (N) or biomass (B) abundance are frequently used as variables to describe populations and species in macroecological and food web studies. 2. We investigate how the use of mean W and mean N or B, rather than other measures of W and/or accounting for the properties of all individuals, can affect the outcome of tests of macroecological and food web theory. 3. Theoretical and empirical analyses demonstrate that mean W, W at maximum biomass (W(mb)), W when energy requirements are greatest (W(me)) and the W when a species uses the greatest proportion of the energy available to all species in a W class (W(mpe)) are not consistently related. 4. For a population at equilibrium, relationships between mean W and W(me) depend on the slope b of the relationship between trophic level and W. For marine fishes, data show that b varies widely among species and thus mean W is an unreliable indicator of the role of a species in the food web. 5. Two different approaches, 'cross-species' and 'all individuals' have been used to estimate slopes of abundance-body mass relationships and to test the energetic equivalence hypothesis and related theory. The approaches, based on relationships between (1) log(10) mean W and log(10) mean N or B, and (2) log(10) W and log(10) N or B of all individuals binned into log(10) W classes (size spectra), give different slopes and confidence intervals with the same data. 6. Our results show that the 'all individuals' approach has the potential to provide more powerful tests of the energetic equivalence hypothesis and role of energy availability in determining slopes, but new theory and empirical analysis are needed to explain distributions of species relative abundance at W. 7. Biases introduced when working with mean W in macroecological and food web studies are greatest when species have indeterminate growth, when relationships between W and trophic level are strong and when the range of species'W is narrow.