Disentangling the influences of mean body size and size structure on ecosystem functioning: an example of nutrient recycling by a non‐native crayfish

Body size is a fundamental functional trait that can be used to forecast individuals' responses to environmental change and their contribution to ecosystem functioning. However, information on the mean and variation of size distributions often confound one another when relating body size to aggregate functioning. Given that size‐based metrics are used as indicators of ecosystem status, it is important to identify the specific aspects of size distributions that mediate ecosystem functioning. Our goal was to simultaneously account for the mean, variance, and shape of size distributions when relating body size to aggregate ecosystem functioning. We take advantage of habitat‐specific differences in size distributions to estimate nutrient recycling by a non‐native crayfish using mean‐field and variance‐incorporating approaches. Crayfishes often substantially influence ecosystem functioning through their omnivorous role in aquatic food webs. As predicted from Jensen's inequality, considering only the mean body size of crayfish overestimated aggregate effects on ecosystem functioning. This bias declined with mean body size such that mean‐field and variance‐incorporating estimates of ecosystem functioning were similar for samples at mean body sizes >7.5 g. At low mean body size, mean‐field bias in ecosystem functioning mismatch predictions from Jensen's inequality, likely because of the increasing skewness of the size distribution. Our findings support the prediction that variance around the mean can alter the relationship between body size and ecosystem functioning, especially at low mean body size. However, methods to account for mean‐field bias performed poorly in samples with highly skewed distributions, indicating that changes in the shape of the distribution, in addition to the variance, may confound mean‐based estimates of ecosystem functioning. Given that many biological functions scale allometrically, explicitly defining and experimentally or statistically isolating the effects of the mean, variance, and shape of size distributions is necessary to begin generalizing relationships between animal body size and ecosystem functioning. Functional traits provide a mechanistic way to scale up the effects of individuals to aggregate measures of ecosystem functioning. However, trait‐based approaches in ecosystem functioning research often consider only the mean of trait distributions, ignoring substantial inter‐ and intraspecific variation. Here we examine a non‐n