Phytoplankton growth and stoichiometric responses to warming,nutrient addition and grazing depend on lake productivity and cell size

Global change involves shifts in multiple environmental factors that act in concert to shape ecological systems in ways that depend on local biotic and abiotic conditions. Little is known about the effects of combined global change stressors on phytoplankton communities, and particularly how these are mediated by distinct community properties such as productivity, grazing pressure and size distribution. Here, we tested for the effects of warming and eutrophication on phytoplankton net growth rate and C:N:P stoichiometry in two phytoplankton cell size fractions (<30 µm and >30 µm) in the presence and absence of grazing in microcosm experiments. Because effects may also depend on lake productivity, we used phytoplankton communities from three Dutch lakes spanning a trophic gradient. We measured the response of each community to multifactorial combinations of temperature, nutrient, and grazing treatments and found that nutrients elevated net growth rates and reduced carbon:nutrient ratios of all three phytoplankton communities. Warming effects on growth and stoichiometry depended on nutrient supply and lake productivity, with enhanced growth in the most productive community dominated by cyanobacteria, and strongest stoichiometric responses in the most oligotrophic community at ambient nutrient levels. Grazing effects were also most evident in the most oligotrophic community, with reduced net growth rates and phytoplankton C:P stoichiometry that suggests consumer‐driven nutrient recycling. Our experiments indicate that stoichiometric responses to warming and interactions with nutrient addition and grazing are not universal but depend on lake productivity and cell size distribution.     

Temperature‐dependent body size effects determine population responses to climate warming

Current understanding of animal population responses to rising temperatures is based on the assumption that biological rates such as metabolism, which governs fundamental ecological processes, scale independently with body size and temperature, despite empirical evidence for interactive effects. Here, we investigate the consequences of interactive temperature‐ and size scaling of vital rates for the dynamics of populations experiencing warming using a stage‐structured consumer‐resource model. We show that interactive scaling alters population and stage‐specific responses to rising temperatures, such that warming can induce shifts in population regulation and stage‐structure, influence community structure and govern population responses to mortality. Analysing experimental data for 20 fish species, we found size–temperature interactions in intraspecific scaling of metabolic rate to be common. Given the evidence for size–temperature interactions and the ubiquity of size structure in animal populations, we argue that accounting for size‐specific temperature effects is pivotal for understanding how warming affects animal populations and communities.     

Do marine phytoplankton follow Bergmann's rule sensu lato?

Global warming has revitalized interest in the relationship between body size and temperature, proposed by Bergmann's rule 150 years ago, one of the oldest manifestations of a ‘biogeography of traits’. We review biogeographic evidence, results from clonal cultures and recent micro‐ and mesocosm experiments with naturally mixed phytoplankton communities regarding the response of phytoplankton body size to temperature, either as a single factor or in combination with other factors such as grazing, nutrient limitation, and ocean acidification. Where possible, we also focus on the comparison between intraspecific size shifts and size shifts resulting from changes in species composition. Taken together, biogeographic evidence, community‐level experiments and single‐species experiments indicate that phytoplankton average cell sizes tend to become smaller in warmer waters, although temperature is not necessarily the proximate environmental factor driving size shifts. Indirect effects via nutrient supply and grazing are important and often dominate. In a substantial proportion of field studies, resource availability is seen as the only factor of relevance. Interspecific size effects are greater than intraspecific effects. Direct temperature effects tend to be exacerbated by indirect ones, if warming leads to intensified nutrient limitation or copepod grazing while ocean acidification tends to counteract the temperature effect on cell size in non‐calcifying phytoplankton. We discuss the implications of the temperature‐related size trends in a global‐warming context, based on known functional traits associated with phytoplankton size. These are a higher affinity for nutrients of smaller cells, highest maximal growth rates of moderately small phytoplankton (ca. 10² µm³), size‐related sensitivities for different types of grazers, and impacts on sinking rates. For a phytoplankton community increasingly dominated by smaller algae we predict that: (i) a higher proportion of primary production will be respired within the microbial food web; (ii) a smaller share of primary production will be channeled to the classic phytoplankton – crustacean zooplankton – fish food chain, thus leading to decreased ecological efficiency from a fish‐production point of view; (iii) a smaller share of primary production will be exported through sedimentation, thus leading to decreased efficiency of the biological carbon pump.     

An allometric model of home range formation explains the structuring of animal communities exploiting heterogeneous resources

Understanding and predicting the composition and spatial structure of communities is a central challenge in ecology. An important structural property of animal communities is the distribution of individual home ranges. Home range formation is controlled by resource heterogeneity, the physiology and behaviour of individual animals, and their intra‐ and interspecific interactions. However, a quantitative mechanistic understanding of how home range formation influences community composition is still lacking. To explore the link between home range formation and community composition in heterogeneous landscapes we combine allometric relationships for physiological properties with an algorithm that selects optimal home ranges given locomotion costs, resource depletion and competition in a spatially‐explicit individual‐based modelling framework. From a spatial distribution of resources and an input distribution of animal body mass, our model predicts the size and location of individual home ranges as well as the individual size distribution (ISD) in an animal community. For a broad range of body mass input distributions, including empirical body mass distributions of North American and Australian mammals, our model predictions agree with independent data on the body mass scaling of home range size and individual abundance in terrestrial mammals. Model predictions are also robust against variation in habitat productivity and landscape heterogeneity. The combination of allometric relationships for locomotion costs and resource needs with resource competition in an optimal foraging framework enables us to scale from individual properties to the structure of animal communities in heterogeneous landscapes. The proposed spatially‐explicit modelling concept not only allows for detailed investigation of landscape effects on animal communities, but also provides novel insights into the mechanisms by which resource competition in space shapes animal communities.     

Energetic determinants of abundance in winter landbird communities

There is increasing evidence that individual energetics constrain macroecological patterns. Here we model total abundance within winter landbird communities as a function of (1) energy supply, as measured by ecosystem net primary productivity, and (2) energy use of individuals, as influenced by body mass and ambient temperature. Using data from the North American Christmas Bird Count, we find that total abundance increases with productivity to the 0.61 power, and decreases with body mass and environmental temperature as predicted by metabolic theory when individuals are below their thermoneutral zone. We note a negative relationship between ambient temperature and average body mass, and suggest that this community‐level pattern, reminiscent of Bergmann's Rule, is related to a tendency for small species to be less‐abundant or absent from cold locations. Results from this study emphasize the importance of individual‐level metabolism for understanding large‐scale ecological patterns.     

Isometric size-scaling of metabolic rate and the size abundance distribution of phytoplankton

The relationship between phytoplankton cell size and abundance has long been known to follow regular, predictable patterns in near steady-state ecosystems, but its origin has remained elusive. To explore the linkage between the size-scaling of metabolic rate and the size abundance distribution of natural phytoplankton communities, we determined simultaneously phytoplankton carbon fixation rates and cell abundance across a cell volume range of over six orders of magnitude in tropical and subtropical waters of the Atlantic Ocean. We found an approximately isometric relationship between carbon fixation rate and cell size (mean slope value: 1.16; range: 1.03-1.32), negating the idea that Kleiber's law is applicable to unicellular autotrophic protists. On the basis of the scaling of individual resource use with cell size, we predicted a reciprocal relationship between the size-scalings of phytoplankton metabolic rate and abundance. This prediction was confirmed by the observed slopes of the relationship between phytoplankton abundance and cell size, which have a mean value of -1.15 (range: -1.29 to -0.97), indicating that the size abundance distribution largely results from the size-scaling of metabolic rate. Our results imply that the total energy processed by carbon fixation is constant along the phytoplankton size spectrum in near steady-state marine ecosystems.     

High predictability of direct competition between marine diatoms under different temperatures and nutrient states

The distribution of marine phytoplankton will shift alongside changes in marine environments, leading to altered species frequencies and community composition. An understanding of the response of mixed populations to abiotic changes is required to adequately predict how environmental change may affect the future composition of phytoplankton communities. This study investigated the growth and competitive ability of two marine diatoms, Phaeodactylum tricornutum and Thalassiosira pseudonana, along a temperature gradient (9–35°C) spanning the thermal niches of both species under both high‐nitrogen nutrient‐replete and low‐nitrogen nutrient‐limited conditions. Across this temperature gradient, the competitive outcome under both nutrient conditions at any assay temperature, and the critical temperature at which competitive advantage shifted from one species to the other, was well predicted by the temperature dependencies of the growth rates of the two species measured in monocultures. The temperature at which the competitive advantage switched from P. tricornutum to T. pseudonana increased from 18.8°C under replete conditions to 25.3°C under nutrient‐limited conditions. Thus, P. tricornutum was a better competitor over a wider temperature range in a low N environment. Being able to determine the competitive outcomes from physiological responses of single species to environmental changes has the potential to significantly improve the predictive power of phytoplankton spatial distribution and community composition models.     

Integrating elements and energy through the metabolic dependencies of gross growth efficiency and the threshold elemental ratio

Metabolic theory proposes that individual growth is governed through the mass‐ and temperature‐dependence of metabolism, and ecological stoichiometry posits that growth is maximized at consumer‐specific optima of resource elemental composition. A given consumer's optimum, the threshold elemental ratio (TER), is proportional to the ratio of its maximum elemental gross growth efficiencies (GGEs). GGE is defined by the ratio of metabolism‐dependent processes such that GGEs should be independent of body mass and temperature. Understanding the metabolic‐dependencies of GGEs and TERs may open the path towards a theoretical framework integrating the flow of energy and chemical elements through ecosystems. However, the mass and temperature scaling of GGEs and TERs have not been broadly evaluated. Here, we use data from 95 published studies to evaluate these metabolic‐dependencies for C, N and P from unicells to vertebrates. We show that maximum GGEs commonly decline as power functions of asymptotic body mass and exponential functions of temperature. The rates of change in maximum GGEs with mass and temperature are relatively slow, however, suggesting that metabolism may not causally influence maximum GGEs. We additionally derived the theoretical expectation that the TER for C:P should not vary with body mass and this was supported empirically. A strong linear relationship between carbon and nitrogen GGEs further suggests that variation in the TER for C:N should be due to variation in consumer C:N. In general we show that GGEs may scale with metabolic rate, but it is unclear if there is a causal link between metabolism and GGEs. Further integrating stoichiometry and metabolism will provide better understanding of the processes governing the flow of energy and elements from organisms to ecosystems.     

Global warming may disproportionately affect larger adults in a predatory coral reef fish

Global warming is expected to reduce body sizes of ectothermic animals. Although the underlying mechanisms of size reductions remain poorly understood, effects appear stronger at latitudinal extremes (poles and tropics) and in aquatic rather than terrestrial systems. To shed light on this phenomenon, we examined the size dependence of critical thermal maxima (CTmax) and aerobic metabolism in a commercially important tropical reef fish, the leopard coral grouper (Plectropomus leopardus) following acclimation to current‐day (28.5 °C) vs. projected end‐of‐century (33 °C) summer temperatures for the northern Great Barrier Reef (GBR). CTmax declined from 38.3 to 37.5 °C with increasing body mass in adult fish (0.45–2.82 kg), indicating that larger individuals are more thermally sensitive than smaller conspecifics. This may be explained by a restricted capacity for large fish to increase mass‐specific maximum metabolic rate (MMR) at 33 °C compared with 28.5 °C. Indeed, temperature influenced the relationship between metabolism and body mass (0.02–2.38 kg), whereby the scaling exponent for MMR increased from 0.74 ± 0.02 at 28.5 °C to 0.79 ± 0.01 at 33 °C, and the corresponding exponents for standard metabolic rate (SMR) were 0.75 ± 0.04 and 0.80 ± 0.03. The increase in metabolic scaling exponents at higher temperatures suggests that energy budgets may be disproportionately impacted in larger fish and contribute to reduced maximum adult size. Such climate‐induced reductions in body size would have important ramifications for fisheries productivity, but are also likely to have knock‐on effects for trophodynamics and functioning of ecosystems.     

Food web framework for size-structured populations

We synthesise traditional unstructured food webs, allometric body size scaling, trait-based modelling, and physiologically structured modelling to provide a novel and ecologically relevant tool for size-structured food webs. The framework allows food web models to include ontogenetic growth and life-history omnivory at the individual level by resolving the population structure of each species as a size-spectrum. Each species is characterised by the trait ‘size at maturation’, and all model parameters are made species independent through scaling with individual body size and size at maturation. Parameter values are determined from cross-species analysis of fish communities as life-history omnivory is widespread in aquatic systems, but may be reparameterised for other systems. An ensemble of food webs is generated and the resulting communities are analysed at four levels of organisation: community level, species level, trait level, and individual level. The model may be solved analytically by assuming that the community spectrum follows a power law. The analytical solution provides a baseline expectation of the results of complex food web simulations, and agrees well with the predictions of the full model on biomass distribution as a function of individual size, biomass distribution as a function of size at maturation, and relation between predator-prey mass ratio of preferred and eaten food. The full model additionally predicts the diversity distribution as a function of size at maturation.     

Linking community size structure and ecosystem functioning using metabolic theory

Understanding how biogeochemical cycles relate to the structure of ecological communities is a central research question in ecology. Here we approach this problem by focusing on body size, which is an easily measured species trait that has a pervasive influence on multiple aspects of community structure and ecosystem functioning. We test the predictions of a model derived from metabolic theory using data on ecosystem metabolism and community size structure. These data were collected as part of an aquatic mesocosm experiment that was designed to simulate future environmental warming. Our analyses demonstrate significant linkages between community size structure and ecosystem functioning, and the effects of warming on these links. Specifically, we show that carbon fluxes were significantly influenced by seasonal variation in temperature, and yielded activation energies remarkably similar to those predicted based on the temperature dependencies of individual-level photosynthesis and respiration. We also show that community size structure significantly influenced fluxes of ecosystem respiration and gross primary production, particularly at the annual time-scale. Assessing size structure and the factors that control it, both empirically and theoretically, therefore promises to aid in understanding links between individual organisms and biogeochemical cycles, and in predicting the responses of key ecosystem functions to future environmental change.     

Zooming in on size distribution patterns underlying species coexistence in Baltic Sea phytoplankton

Scale is a key to determining which processes drive community structure. We analyse size distributions of phytoplankton to determine time scales at which we can observe either fixed environmental characteristics underlying communities structure or competition‐driven size distributions. Using multiple statistical tests, we characterise size distributions of phytoplankton from 20‐year time series in two sites of the Baltic Sea. At large temporal scales (5–20 years), size distributions are unimodal, indicating that fundamental barriers to existence are here subtler than in other systems. Frequency distributions of the average size of the species weighted by biovolume are multimodal over large time scales, although this is the product of often unimodal short‐term (<1 year) patterns. Our study represents a much‐needed structured, high‐resolution analysis of phytoplankton size distributions, revealing that short‐term analyses are necessary to determine if, and how, competition shapes them. Our results provide a stepping‐stone on which to further investigate the intricacies of competition and coexistence.     

Metabarcoding and metabolome analyses of copepod grazing reveal feeding preference and linkage to metabolite classes in dynamic microbial plankton communities

In order to characterize copepod feeding in relation to microbial plankton community dynamics, we combined metabarcoding and metabolome analyses during a 22‐day seawater mesocosm experiment. Nutrient amendment of mesocosms promoted the development of haptophyte (Phaeocystis pouchetii)‐ and diatom (Skeletonema marinoi)‐dominated plankton communities in mesocosms, in which Calanus sp. copepods were incubated for 24 h in flow‐through chambers to allow access to prey particles (<500 μm). Copepods and mesocosm water sampled six times spanning the experiment were analysed using metabarcoding, while intracellular metabolite profiles of mesocosm plankton communities were generated for all experimental days. Taxon‐specific metabarcoding ratios (ratio of consumed prey to available prey in the surrounding seawater) revealed diverse and dynamic copepod feeding selection, with positive selection on large diatoms, heterotrophic nanoflagellates and fungi, while smaller phytoplankton, including P. pouchetii, were passively consumed or even negatively selected according to our indicator. Our analysis of the relationship between Calanus grazing ratios and intracellular metabolite profiles indicates the importance of carbohydrates and lipids in plankton succession and copepod–prey interactions. This molecular characterization of Calanus sp. grazing therefore provides new evidence for selective feeding in mixed plankton assemblages and corroborates previous findings that copepod grazing may be coupled to the developmental and metabolic stage of the entire prey community rather than to individual prey abundances.     

Cell size trade-offs govern light exploitation strategies in marine phytoplankton

Marine phytoplankton show complex community structures and biogeographic distributions, the net results of physiological and ecological trade-offs of species responses to fluctuating, heterogeneous environments. We analysed photosynthesis, responses to variable light and macromolecular allocations across a size panel of marine centric diatoms. The diatoms have strong capacities to withstand and exploit fluctuating light, when compared with picophytoplankton. Within marine diatoms, small species show larger effective cross-sections for photochemistry, and fast metabolic repair of photosystem II after photoinactivation. In contrast, large diatoms show lower susceptibility to photoinactivation, and therefore incur lower costs to endure short-term exposures to high light, especially under conditions that limit metabolic rates. This size scaling of key photophysiological parameters thus helps explain the relative competitive advantages of larger versus smaller species under different environmental regimes, with implications for the function of the biogenic carbon pump. These results provide a mechanistic framework to explain and predict shifts in marine phytoplankton community size structure with changes in surface irradiance and mixed layer depth.     

Elevated temperature alters proteomic responses of individual organisms within a biofilm community

Microbial communities that underpin global biogeochemical cycles will likely be influenced by elevated temperature associated with environmental change. Here, we test an approach to measure how elevated temperature impacts the physiology of individual microbial groups in a community context, using a model microbial-based ecosystem. The study is the first application of tandem mass tag (TMT)-based proteomics to a microbial community. We accurately, precisely and reproducibly quantified thousands of proteins in biofilms growing at 40, 43 and 46 °C. Elevated temperature led to upregulation of proteins involved in amino-acid metabolism at the level of individual organisms and the entire community. Proteins from related organisms differed in their relative abundance and functional responses to temperature. Elevated temperature repressed carbon fixation proteins from two Leptospirillum genotypes, whereas carbon fixation proteins were significantly upregulated at higher temperature by a third member of this genus. Leptospirillum group III bacteria may have been subject to viral stress at elevated temperature, which could lead to greater carbon turnover in the microbial food web through the release of viral lysate. Overall, these findings highlight the utility of proteomics-enabled community-based physiology studies, and provide a methodological framework for possible extension to additional mixed culture and environmental sample analyses.     

Metabolic constraints on the body size scaling of extreme population densities

Pest outbreaks, harmful algal blooms and population collapses are extreme events with critical consequences for ecosystems. Therefore, understanding the ecological mechanisms underlying these extreme events is crucial. We evaluated theoretical predictions on the size scaling and variance of extreme population abundance by combining (i) the generalized extreme value (GEV) theory and (ii) the resource-limited metabolic restriction hypothesis for population abundance. Using the phytoplankton data from the L4 station in the English Channel, we showed a negative size scaling of the expected value of maximal density, whose confidence interval included the predicted metabolic scaling (α = −1) supporting theoretical predictions. The role of resources and temperature in the distribution of the size–abundance pattern and residuals was well characterized by the GEV distribution. This comprehensive modelling framework will allow to elucidate community structure and fluctuations and provide unbiased return times estimates, thereby improving the prediction accuracy of the timing of the population outbreaks.     

Ecological emergence of thermal clines in body size

The unprecedented rate of global warming requires a better understanding of how ecosystems will respond. Organisms often have smaller body sizes under warmer climates (Bergmann's rule and the temperature‐size rule), and body size is a major determinant of life histories, demography, population size, nutrient turnover rate, and food‐web structure. Therefore, by altering body sizes in whole communities, current warming can potentially disrupt ecosystem function and services. However, the underlying drivers of warming‐induced body downsizing remain far from clear. Here, we show that thermal clines in body size are predicted from universal laws of ecology and metabolism, so that size‐dependent selection from competition (both intra and interspecific) and predation favors smaller individuals under warmer conditions. We validate this prediction using 4.1 × 10⁶ individual body size measurements from French river fish spanning 29 years and 52 species. Our results suggest that warming‐induced body downsizing is an emergent property of size‐structured food webs, and highlight the need to consider trophic interactions when predicting biosphere reorganizations under global warming.     

Increasing importance of small phytoplankton in a warmer ocean

The macroecological relationships among marine phytoplankton total cell density, community size structure and temperature have lacked a theoretical explanation. The tiniest members of this planktonic group comprise cyanobacteria and eukaryotic algae smaller than 2 μm in diameter, collectively known as picophytoplankton. We combine here two ecological rules, the temperature–size relationship with the allometric size‐scaling of population abundance to explain a remarkably consistent pattern of increasing picophytoplankton biomass with temperature over the ?0.6 to 22 °C range in a merged dataset obtained in the eastern and western temperate North Atlantic Ocean across a diverse range of environmental conditions. Our results show that temperature alone was able to explain 73% of the variance in the relative contribution of small cells to total phytoplankton biomass regardless of differences in trophic status or inorganic nutrient loading. Our analysis predicts a gradual shift toward smaller primary producers in a warmer ocean. Because the fate of photosynthesized organic carbon largely depends on phytoplankton size, we anticipate future alterations in the functioning of oceanic ecosystems.     

Litter elemental stoichiometry and biomass densities of forest soil invertebrates

To maintain constant chemical composition, i.e. elemental homeostasis, organisms have to consume resources of sufficient quality to meet their own specific stoichiometric demand. Therefore, concentrations of elements indicate resource quality, and rare elements in the environment may act as limiting factors for individual organisms scaling up to constrain population densities. We investigated how the biomass densities of invertebrate populations of temperate forest soil communities depend on 1) the stoichiometry of the basal litter according to ecological stoichiometry concepts and 2) the population average body mass as predicted by metabolic theory. We used a large data set on biomass densities of 4959 populations across 48 forests in three regions of Germany. Following various ecological stoichiometry hypotheses, we tested for effects of the carbon‐to‐element ratios of 10 elements. Additionally, we included the abiotic litter characteristics habitat size (represented by litter depth), litter diversity and pH, as well as forest type as an indicator for human management. Across 12 species groups, we found that the biomass densities scaled significantly with population‐averaged body masses thus supporting metabolic theory. Additionally, 10 of these allometric scaling relationships exhibited interactions with stoichiometric and abiotic co‐variables. The four most frequent co‐variables were 1) forest type, 2) the carbon‐to‐phosphorus ratio (C:P), 3) the carbon‐to‐sodium ratio (C:Na), and the carbon‐to‐nitrogen ratio (C:N). Hence, our analyses support the sodium shortage hypothesis for microbi‐detritivores, the structural elements hypothesis for some predator groups (concerning N), and the secondary productivity hypothesis (concerning P) across all trophic groups in our data. In contrast, the ecosystem size hypothesis was only supported for some meso‐ and macrofauna detritivores. Our study is thus providing a comprehensive analysis how the elemental stoichiometry of the litter as the basal resource constrain population densities across multiple trophic levels of soil communities.     

Spatial and body‐size dependent response of marine pelagic communities to projected global climate change

Temperature, oxygen, and food availability directly affect marine life. Climate models project a global warming of the ocean's surface (~+3 °C), a de‐oxygenation of the ocean's interior (~?3%) and a decrease in total marine net primary production (~?8%) under the ‘business as usual’ climate change scenario (RCP8.5). We estimated the effects of these changes on biological communities using a coupled biogeochemical (PISCES) – ecosystems (APECOSM) model forced by the physical outputs of the last generation of the IPSL‐CM Earth System Model. The APECOSM model is a size‐structured bio‐energetic model that simulates the 3D dynamical distributions of three interactive pelagic communities (epipelagic, mesopelagic, and migratory) under the effects of multiple environmental factors. The PISCES‐APECOSM model ran from 1850 to 2100 under historical forcing followed by RCP8.5. Our RCP8.5 simulation highlights significant changes in the spatial distribution, biomass, and maximum body‐size of the simulated pelagic communities. Biomass and maximum body‐size increase at high latitude over the course of the century, reflecting the capacity of marine organisms to respond to new suitable environment. At low‐ and midlatitude, biomass and maximum body‐size strongly decrease. In those regions, large organisms cannot maintain their high metabolic needs because of limited and declining food availability. This resource reduction enhances the competition and modifies the biomass distribution among and within the three communities: the proportion of small organisms increases in the three communities and the migrant community that initially comprised a higher proportion of small organisms is favored. The greater resilience of small body‐size organisms resides in their capacity to fulfill their metabolic needs under reduced energy supply and is further favored by the release of predation pressure due to the decline of large organisms. These results suggest that small body‐size organisms might be more resilient to climate change than large ones.