Marine biodiversity and ecosystem functioning: what's known and what's next?

Marine eDNA sampling from submerged surfaces with paint rollers

Environmental DNA (eDNA) analyses of species present in marine environments is the most effective biological diversity measurement tool currently available. eDNA sampling methods are an intrinsically important part of the eDNA biodiversity analysis process. Identification and development of eDNA sampling methods that are as rapid, affordable, versatile and practical as possible will improve rates of detection of marine species. Optimal outcomes of eDNA biodiversity surveys come from studies employing high levels of sampling replication, so any methods that make sampling faster and cheaper will improve scientific outcomes. eDNA sampling methods that can be applied more widely will also enable sampling from a greater range of marine surface micro-habitats, resulting in detection of a wider range of organisms. In this study, we compared diversity detection by several methods for sampling eDNA from submerged marine surfaces: polyurethane foam, nylon swabs, microfibre paint rollers, and sediment scoops. All of the methods produced a diverse range of species identifications, with >250 multicellular species represented by eDNA at the study site. We found that widely-available small paint rollers were an effective, readily available and affordable method for sampling eDNA from underwater marine surfaces. This approach enables the sampling of marine eDNA using extended poles, or potentially by remotely operated vehicles, where surface sampling by hand is impractical.

Biogeographic and seasonal differences in consumer pressure underlie strong predation in the tropics

Biotic interactions play a critical role in shaping patterns of global biodiversity. While several macroecological studies provide evidence for stronger predation in tropical regions compared with higher latitudes, results are variable even within the tropics, and the drivers of this variability are not well understood. We conducted two complementary standardized experiments on communities of sessile marine invertebrate prey and their associated predators to test for spatial and seasonal differences in predation across the tropical Atlantic and Pacific coastlines of Panama. We further tested the prediction that higher predator diversity contributes to stronger impacts of predation, using both direct observations of predators and data from extensive reef surveys. Our results revealed substantially higher predation rates and stronger effects of predators on prey in the Pacific than in the Atlantic, demonstrating striking variation within tropical regions. While regional predator diversity was high in the Atlantic, functional diversity at local scales was markedly low. Peak predation strength in the Pacific occurred during the wet, non-upwelling season when ocean temperatures were warmer and predator communities were more functionally diverse. Our results highlight the importance of regional biotic and abiotic drivers that shape interaction strength and the maintenance of tropical communities, which are experiencing rapid environmental change.

Transient species driving ecosystem multifunctionality: Insights from competitive interactions between rocky intertidal mussels

Anthropogenic biodiversity loss poses a significant threat to ecosystem functioning worldwide. Numerically dominant and locally rare (i.e., transient) species are key components of biodiversity, but their contribution to multiple ecosystem functions (i.e., multifunctionality) has been seldomly assessed in marine ecosystems. To fill this gap, here we analyze the effects of a dominant and a transient species on ecosystem multifunctionality. In an observational study conducted along ca. 200 km of the southeastern Pacific coast, the purple mussel Perumytilus purpuratus numerically dominated the mid-intertidal and the dwarf mussel Semimytilus patagonicus exhibited low abundances but higher recruitment rates. In laboratory experiments, the relative abundances of both species were manipulated to simulate the replacement of P. purpuratus by S. patagonicus and five proxies for ecosystem functions-rates of clearance, oxygen consumption, total biodeposit, organic biodeposit, and excretion-were analyzed. This replacement had a positive, linear, and significant effect on the combined ecosystem functions, particularly oxygen consumption and excretion rates. Accordingly, S. patagonicus could well drive ecosystem functioning given favorable environmental conditions for its recovery from rarity. Our study highlights therefore the key role of transient species for ecosystem performance. Improving our understanding of these dynamics is crucial for effective ecosystem conservation, especially in the current scenario of biological extinctions and invasions.

A meta-evaluation of the quality of reporting and execution in ecological meta-analyses

Quantitatively summarizing results from a collection of primary studies with meta-analysis can help answer ecological questions and identify knowledge gaps. The accuracy of the answers depends on the quality of the meta-analysis. We reviewed the literature assessing the quality of ecological meta-analyses to evaluate current practices and highlight areas that need improvement. From each of the 18 review papers that evaluated the quality of meta-analyses, we calculated the percentage of meta-analyses that met criteria related to specific steps taken in the meta-analysis process (i.e., execution) and the clarity with which those steps were articulated (i.e., reporting). We also re-evaluated all the meta-analyses available from Pappalardo et al. [1] to extract new information on ten additional criteria and to assess how the meta-analyses recognized and addressed non-independence. In general, we observed better performance for criteria related to reporting than for criteria related to execution; however, there was a wide variation among criteria and meta-analyses. Meta-analyses had low compliance with regard to correcting for phylogenetic non-independence, exploring temporal trends in effect sizes, and conducting a multifactorial analysis of moderators (i.e., explanatory variables). In addition, although most meta-analyses included multiple effect sizes per study, only 66% acknowledged some type of non-independence. The types of non-independence reported were most often related to the design of the original experiment (e.g., the use of a shared control) than to other sources (e.g., phylogeny). We suggest that providing specific training and encouraging authors to follow the PRISMA EcoEvo checklist recently developed by O'Dea et al. [2] can improve the quality of ecological meta-analyses.

Sponge functional roles in a changing world

Sponges are ecologically important benthic organisms with many important functional roles. However, despite increasing global interest in the functions that sponges perform, there has been limited focus on how such functions will be impacted by different anthropogenic stressors. In this review, we describe the progress that has been made in our understanding of the functional roles of sponges over the last 15 years and consider the impacts of anthropogenic stressors on these roles. We split sponge functional roles into interactions with the water column and associations with other organisms. We found evidence for an increasing focus on functional roles among sponge-focused research articles, with our understanding of sponge-mediated nutrient cycling increasing substantially in recent years. From the information available, many anthropogenic stressors have the potential to negatively impact sponge pumping, and therefore have the potential to cause ecosystem level impacts. While our understanding of the importance of sponges has increased in the last 15 years, much more experimental work is required to fully understand how sponges will contribute to reef ecosystem function in future changing oceans.

Evenness, biodiversity, and ecosystem function of intertidal communities along the Italian coasts: Experimental short-term response to ambient and extreme air temperatures

Biodiversity can promote ecosystem functioning in both terrestrial and marine environments, emphasizing the necessity of biodiversity conservation in order to preserve critical ecosystem functions and associated services. However, the role of biodiversity in buffering ecosystem functioning under extreme events caused by climate change remains a major scientific issue, especially for intertidal systems experiencing stressors from both terrestrial and marine drivers. We performed a regional-scale field experiment along the Italian coast to investigate the response of unmanipulated intertidal communities (by using a natural biodiversity gradient) to low tide aerial exposure to both ambient and short-term extreme temperatures. We specifically investigated the relationship between Biodiversity and Ecosystem Functioning (BEF) using different biodiversity indexes (species richness, functional diversity and evenness) and the response of the intertidal communities' ecosystem functioning (community respiration rates). Furthermore, we investigated which other environmental variables could influence the BEF relationship. We show that evenness explained a greater variation in intertidal community ecosystem functioning under both temperature conditions. Species richness (the most often used diversity metric in BEF research) was unrelated to ecosystem functioning, while functional diversity was significantly related to respiration under ambient but not extreme temperatures. We highlight the importance of the short-term thermal history of the communities (measured as body temperature) in the BEF relationship as it was consistently identified as the best predictor or response under both temperature conditions. However, Chlorophyll a in seawater and variation in sea surface temperature also contributed to the BEF relationship under ambient but not under extreme conditions, showing that short-duration climate-driven events can overcome local physiological adaptations. Our findings support the importance of the BEF relationship in intertidal communities, implying that systems with more diverse and homogeneous communities may be able to mitigate the effects of extreme temperatures.

Foundational biodiversity effects propagate through coastal food webs via multiple pathways

Relatively few studies have attempted to resolve the pathways through which the effects of biodiversity on ecosystem functioning cascade from one trophic level to another. Here, we manipulated the richness of habitat-forming seaweeds in a western Atlantic estuary to explore how changes in foundation species diversity affect the structure and functioning of the benthic consumer communities that they support. Structural equation modeling revealed that macroalgal richness enhanced invertebrate abundance, biomass, and diversity, both directly by changing the quality and palatability of the foundational substrate and indirectly by increasing the total biomass of available habitat. Consumer responses were largely driven by a single foundational seaweed, although stronger complementarity among macroalgae was observed for invertebrate richness. These findings with diverse foundational phyla extend earlier inferences from terrestrial grasslands by showing that biodiversity effects can simultaneously propagate through multiple independent pathways to maintain animal foodwebs. Our work also highlights the potential ramifications of human-induced changes in marine ecosystems.

Benthic ecosystem functioning under climate change: modelling the bioturbation potential for benthic key species in the southern North Sea

Climate change affects the marine environment on many levels with profound consequences for numerous biological, chemical, and physical processes. Benthic bioturbation is one of the most relevant and significant processes for benthic-pelagic coupling and biogeochemical fluxes in marine sediments, such as the uptake, transport, and remineralisation of organic carbon. However, only little is known about how climate change affects the distribution and intensity of benthic bioturbation of a shallow temperate shelf sea system such as the southern North Sea. In this study, we modelled and projected changes in bioturbation potential (BP_p) under a continuous global warming scenario for seven southern North Sea key bioturbators: Abra alba, Amphiura filiformis, Callianassa subterranea, Echinocardium cordatum, Goniada maculata, Nephtys hombergii, and Nucula nitidosa. Spatial changes in species bioturbation intensity are simulated for the years 2050 and 2099 based on one species distribution model per species driven by bottom temperature and salinity changes using the IPCC SRES scenario A1B. Local mean bottom temperature was projected to increase between 0.15 and 5.4 °C, while mean bottom salinity was projected to moderately decrease by 1.7. Our results show that the considered benthic species are strongly influenced by the temperature increase. Although the total BP remained rather constant in the southern North Sea, the BP_p for four out of seven species was projected to increase, mainly due to a simultaneous northward range expansion, while the BP_p in the core area of the southern North Sea declined for the same species. Bioturbation of the most important species, Amphiura filiformis and Echinocardium cordatum, showed no substantial change in the spatial distribution, but over time. The BP_p of E. cordatum remained almost constant until 2099, while the BP_p of A. filiformis decreased by 41%. The northward expansion of some species and the decline of most species in the south led to a change of relative contribution to bioturbation in the southern North Sea. These results indicate that some of the selected key bioturbators in the southern North Sea might partly compensate the decrease in bioturbation by others. But especially in the depositional areas where bioturbation plays a specifically important role for ecosystem functioning, bioturbation potential declined until 2099, which might affect the biochemical cycling in sediments of some areas of the southern North Sea.

Stressors Increase the Impacts of Coastal Macrofauna Biodiversity Loss on Ecosystem Multifunctionality

There is substantial evidence that biodiversity underpins ecosystem functioning, but it is unclear how these relationships change with multiple stressors in complex real-world settings. Coastal zones are affected by numerous stressors (for example, sediment input and nutrient runoff from land) and the cumulative effects of these stressors may result in pronounced and unexpected changes in the functioning of ecosystems. To investigate the cumulative effects of turbidity and elevated nutrients on coastal biodiversity-ecosystem functioning relationships, we performed a large-scale field experiment manipulating in situ sediment porewater ammonium concentrations and measured multiple ecosystem functions related to carbon fixation and mineralisation in 15 estuaries with varying levels of turbidity. The results indicated that the benthic macrofauna diversity (species richness, abundance, and functional richness) declined with increased porewater ammonium concentrations and there were clear thresholds in light levels at the seafloor in relation to the biodiversity-ecosystem function relationships. Multifunctionality indices (an integrated index of all measured functions) in moderately turbid and turbid estuaries (daily mean seafloor PAR < 420 µmol m−2 s−1) decreased with the loss of macrofauna biodiversity. Functioning in low-turbidity estuaries (daily mean PAR > 420 µmol m−2 s−1) however remained relatively constant, suggesting that they were more resilient against the nutrient-induced biodiversity loss. Our results demonstrate that ecosystems already stressed by stressors that alter functional performance (turbidity) may be more prone to loss of overall functioning if biodiversity is reduced by another stressor (nutrient enrichment), highlighting the potential snowballing effects of cumulative change.

Accounting for Uncertainties in Biodiversity Estimations: A New Methodology and Its Application to the Mesopelagic Sound Scattering Layer of the High Arctic

Measures of biological diversity (biodiversity) are important for monitoring the state of ecosystems. Several indices and methods are used to describe biodiversity from field observations. Marine faunal biodiversity is often quantified based on analysis of samples collected using a trawl during research surveys. To monitor spatial and temporal changes in biodiversity between surveys, samples are generally collected from a series of stations. Inference regarding changes in biodiversity must account for uncertainties in the estimation of the values for the different biodiversity indices used. Estimation for a single station is affected by spatial-temporal variation in the species composition in the area and by uncertainty due to the finite sample size taken by the trawl. Therefore, variation between stations needs to be accounted for when estimating uncertainty for values of different indices during a survey as total or as mean for the survey. Herein, we present a method based on nested bootstrapping that accounts for uncertainties in the estimation of various indices and which can be used to infer changes in biodiversity. Application of this methodology is demonstrated using data collected in the mesopelagic sound scattering layer in the high Arctic.

Metagenomic analysis reveals the different characteristics of microbial communities inside and outside the karst tiankeng

Background

Karst tiankengs serve as a reservoir of biodiversity in the degraded karst landscape areas. However, the microbial diversity of karst tiankengs is poorly understood. Here, we investigated the composition and function of the microbial community in a karst tiankeng.

Results

We found that habitat differences inside and outside the karst tiankeng changed the composition and structure of the soil microbial communities, and the dominant phyla were Proteobacteria, Actinobacteria, Chloroflexi and Acidobacteria. The Shannon–Wiener diversity of microbial communities inside and outside the tiankeng was significantly different, and it was higher inside the tiankeng (IT). Venn and LEfSe analysis found that the soil microbial communities inside the tiankeng had 640 more endemic species and 39 more biomarker microbial clades than those identified outside of the tiankeng (OT)..Functional prediction indicated that soil microorganisms in outside the tiankeng had a high potential for carbohydrate metabolism, translation and amino acid metabolism. There were biomarker pathways associated with several of human diseases at both IT and OT sites. Except for auxiliary activities (AA), other CAZy classes had higher abundance at IT sites, which can readily convert litter and fix carbon and nitrogen, thereby supporting the development of underground forests. The differences in microbial communities were mainly related to the soil water content and soil total nitrogen.

Conclusions

Our results provide a metagenomic overview of the karst tiankeng system and provide new insights into habitat conservation and biodiversity restoration in the area.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-022-02513-1.

Managing human-mediated range shifts: understanding spatial, temporal and genetic variation in marine non-native species

The use of molecular tools to manage natural resources is increasingly common. However, DNA-based methods are seldom used to understand the spatial and temporal dynamics of species' range shifts. This is important when managing range shifting species such as non-native species (NNS), which can have negative impacts on biotic communities. Here, we investigated the ascidian NNS Ciona robusta, Clavelina lepadiformis, Microcosmus squamiger and Styela plicata using a combined methodological approach. We first conducted non-molecular biodiversity surveys for these NNS along the South African coastline, and compared the results with historical surveys. We detected no consistent change in range size across species, with some displaying range stability and others showing range shifts. We then sequenced a section of cytochrome c oxidase subunit I (COI) from tissue samples and found genetic differences along the coastline but no change over recent times. Finally, we found that environmental DNA metabarcoding data showed broad congruence with both the biodiversity survey and the COI datasets, but failed to capture the complete incidence of all NNS. Overall, we demonstrated how a combined methodological approach can effectively detect spatial and temporal variation in genetic composition and range size, which is key for managing both thriving NNS and threatened species.

This article is part of the theme issue ‘Species’ ranges in the face of changing environments (part I)’.

Biodiversity influences the effects of oil disturbance on coastal ecosystems

Biodiversity can enhance the response of ecosystems to disturbance. However, whether diversity can reduce the ecological effect of human-induced novel and extreme disturbances is unclear. In April 2010, the Deepwater Horizon (DwH) platform exploded, allowing an uncontrolled release of crude oil into the northern Gulf of Mexico. Initial surveys following the spill found that ecological impacts on coastal ecosystems varied greatly across habitat-type and trophic group; however, to date, few studies have tested the influence of local biodiversity on these responses. We used a meta-analytic approach to synthesize the results of 5 mesocosm studies that included 10 independent oil experiments and 5 independent oil + dispersant experiments. We tested whether biodiversity increased the resistance and/or resilience of coastal ecosystems to oil disturbance and whether a biodiversity effect depended on the type of diversity present (taxonomic or genetic) and/or the response type measured (population, community, or ecosystem level). We found that diversity can influence the effects of oiling, but the direction and magnitude of this diversity effect varied. Diversity reduced the negative impact of oiling for within-trophic-level responses and tended to be stronger for taxonomic than genetic diversity. Further, diversity effects were largely driven by the presence of highly resistant or quick to recover taxa and genotypes, consistent with the insurance hypothesis. However, we found no effect of diversity on the response to the combination of oil and dispersant exposure. We conclude that areas of low biodiversity may be particularly vulnerable to future oil disturbances and provide insight into the benefit of incorporating multiple measures of diversity in restoration projects and management decisions.

Underwater Hyperspectral Imaging (UHI): A Review of Systems and Applications for Proximal Seafloor Ecosystem Studies

Marine ecosystem monitoring requires observations of its attributes at different spatial and temporal scales that traditional sampling methods (e.g., RGB imaging, sediment cores) struggle to efficiently provide. Proximal optical sensing methods can fill this observational gap by providing observations of, and tracking changes in, the functional features of marine ecosystems non-invasively. Underwater hyperspectral imaging (UHI) employed in proximity to the seafloor has shown a further potential to monitor pigmentation in benthic and sympagic phototrophic organisms at small spatial scales (mm–cm) and for the identification of minerals and taxa through their finely resolved spectral signatures. Despite the increasing number of studies applying UHI, a review of its applications, capabilities, and challenges for seafloor ecosystem research is overdue. In this review, we first detail how the limited band availability inherent to standard underwater cameras has led to a data analysis “bottleneck” in seafloor ecosystem research, in part due to the widespread implementation of underwater imaging platforms (e.g., remotely operated vehicles, time-lapse stations, towed cameras) that can acquire large image datasets. We discuss how hyperspectral technology brings unique opportunities to address the known limitations of RGB cameras for surveying marine environments. The review concludes by comparing how different studies harness the capacities of hyperspectral imaging, the types of methods required to validate observations, and the current challenges for accurate and replicable UHI research.

Stronger predation intensity and impact on prey communities in the tropics

The hypothesis that biotic interactions strengthen toward lower latitudes provides a framework for linking community-scale processes with the macroecological scales that define our biosphere. Despite the importance of this hypothesis for understanding community assembly and ecosystem functioning, the extent to which interaction strength varies across latitude and the effects of this variation on natural communities remain unresolved. Predation in particular is central to ecological and evolutionary dynamics across the globe, yet very few studies explore both community-scale causes and outcomes of predation across latitude. Here we expand beyond prior studies to examine two important components of predation strength: intensity of predation (including multiple dimensions of the predator guild) and impact on prey community biomass and structure, providing one of the most comprehensive examinations of predator-prey interactions across latitude. Using standardized experiments, we tested the hypothesis that predation intensity and impact on prey communities were stronger at lower latitudes. We further assessed prey recruitment to evaluate the potential for this process to mediate predation effects. We used sessile marine invertebrate communities and their fish predators in nearshore environments as a model system, with experiments conducted at 12 sites in four regions spanning the tropics to the subarctic. Our results show clear support for an increase in both predation intensity and impact at lower relative to higher latitudes. The predator guild was more diverse at low latitudes, with higher predation rates, longer interaction durations, and larger predator body sizes, suggesting stronger predation intensity in the tropics. Predation also reduced prey biomass and altered prey composition at low latitudes, with no effects at high latitudes. Although recruitment rates were up to three orders of magnitude higher in the tropics than the subarctic, prey replacement through this process was insufficient to dampen completely the strong impacts of predators in the tropics. Our study provides a novel perspective on the biotic interaction hypothesis, suggesting that multiple components of the predator community likely contribute to predation intensity at low latitudes, with important consequences for the structure of prey communities.

Towards a systematic method for assessing the impact of chemical pollution on ecosystem services of water systems

Chemical pollution impinges on the quality of water systems and the ecosystem services (ESs) they provide. Expression of ESs in monetary units has become an essential tool for sustainable ecosystem management. However, the impact of chemical pollution on ESs is rarely quantified, and ES valuation often focuses on individual services without considering the total services provided by the ecosystem. The purpose of the study was to develop a stepwise approach to quantify the impact of sediment pollution on the total ES value provided by water systems. Thereby, we calculated the total ES value loss as a function of the multi-substance potentially affected fraction of species at the HC50 level (msPAF(HC50)). The function is a combination of relationships between, subsequently: the msPAF(HC50), diversity, productivity and total ES value. Regardless of the inherent differences between terrestrial and aquatic ecosystems, an increase of diversity generally corresponded to an increase in productivity with curvilinear or linear effects. A positive correlation between productivity and total values of ESs of biomes was observed. The combined relationships showed that 1% msPAF(HC50) corresponded to on average 0.5% (0.05-1.40%) of total ES value loss. The ES loss due to polluted sediments in the Waal-Meuse river estuary (the Netherlands) and Flemish waterways (Belgium) was estimated to be 0.3-5 and 0.6-10 thousand 2007$/ha/yr, respectively. Our study presents a novel methodology to assess the impact of chemical exposure on diversity, productivity, and total value that ecosystems provide. With sufficient monitoring data, our generic methodology can be applied for any chemical and region of interest and help water managers make informed decisions on cost-effective measures to remedy pollution. Acknowledging that the ES loss estimates as a function of PAF(HC50) are crude, we explicitly discuss the uncertainties in each step for further development and application of the methodology.

Specialized metabolic functions of keystone taxa sustain soil microbiome stability

BACKGROUND

The relationship between biodiversity and soil microbiome stability remains poorly understood. Here, we investigated the impacts of bacterial phylogenetic diversity on the functional traits and the stability of the soil microbiome. Communities differing in phylogenetic diversity were generated by inoculating serially diluted soil suspensions into sterilized soil, and the stability of the microbiome was assessed by detecting community variations under various pH levels. The taxonomic features and potential functional traits were detected by DNA sequencing.

RESULTS

We found that bacterial communities with higher phylogenetic diversity tended to be more stable, implying that microbiomes with higher biodiversity are more resistant to perturbation. Functional gene co-occurrence network and machine learning classification analyses identified specialized metabolic functions, especially "nitrogen metabolism" and "phosphonate and phosphinate metabolism," as keystone functions. Further taxonomic annotation found that keystone functions are carried out by specific bacterial taxa, including Nitrospira and Gemmatimonas, among others.

CONCLUSIONS

This study provides new insights into our understanding of the relationships between soil microbiome biodiversity and ecosystem stability and highlights specialized metabolic functions embedded in keystone taxa that may be essential for soil microbiome stability. Video abstract.

Linking multiple facets of biodiversity and ecosystem functions in a coastal reef habitat

Reef-building species play key roles in promoting local species richness and regulating ecosystem functions like biogeochemical fluxes. We evaluated the functioning of a habitat engineered by the reef-building polychaete Sabellaria alveolata, by measuring oxygen and nutrient fluxes in the reef structures and in the soft-sediments nearby. Then, we investigated the relative importance of temperature, the engineer S. alveolata, and different facets of macrofauna diversity (taxonomic, functional diversity and identity), on the reef biogeochemical fluxes using multiple linear regressions and effect sizes. The reef fluxes were more intense than the soft-sediment fluxes and mainly driven by the engineer biomass and abundance, stressing the importance of these biogenic structures. Higher water temperatures and an intermediate level of associated macrofauna functional dispersion weighted only by abundance (i.e. intermediate biological trait variability) maximized the reef's global biogeochemical functioning. Ultimately, the physical degradation of the reefs could lead to lower levels of functioning.

Does the use of biological traits predict a smooth landscape of ecosystem functioning?

The biodiversity crisis has increased interest in understanding the role of biodiversity for ecosystem functioning. Functional traits are often used to infer ecosystem functions to increase our understanding of these relationships over larger spatial scales. The links between specific traits and ecosystem functioning are, however, not always well established. We investigated how the choice of analyzing either individual species, selected modalities, or trait combinations affected the spatial patterns observed on a sandflat and how this was related to the natural variability in ecosystem functioning. A large dataset of 400 benthic macrofauna samples was used to explore distribution patterns. We hypothesized that (1) if multiple species (redundancy) represent a trait combination or a modality their spatial patterns would be smoothed out, and (2) the lost spatial variability within a trait combination or modality, due to the smoothing effect, would potentially affect their utility for predicting ecosystem functioning (tested on a dataset of 24 samples). We predicted that species would show heterogeneous small spatial patterns, while modalities and trait combinations would show larger and more homogeneous patterns because they would represent a collection of many distributions. If modalities and trait combinations are better predictors of ecosystem functioning than species, then the smoother spatial patterns of modalities and trait combinations would result in a more homogeneous landscape of ecosystem function and the number of species exhibiting specific traits would provide functional redundancy. Our results showed some smoothing of spatial patterns progressing from species through modalities to trait combinations, but generally spatial patterns reflected a few dominant key species. Moreover, some individual modalities and species explained more or equal proportions of the variance in the ecosystem functioning than the combined traits. The findings thus suggest that only some spatial variability is lost when species are combined into modalities and trait combinations and that a homogeneous landscape of ecosystem function is not likely.

Large seasonal and spatial variation in nano- and microphytoplankton diversity along a Baltic Sea-North Sea salinity gradient

Aquatic phytoplankton experience large fluctuations in environmental conditions during seasonal succession and across salinity gradients, but the impact of this variation on their diversity is poorly understood. We examined spatio-temporal variation in nano- and microphytoplankton (> 2 µm) community structure using almost two decades of light-microscope based monitoring data. The dataset encompasses 19 stations that span a salinity gradient from 2.8 to 35 along the Swedish coastline. Spatially, both regional and local phytoplankton diversity increased with broad-scale salinity variation. Diatoms dominated at high salinity and the proportion of cyanobacteria increased with decreasing salinity. Temporally, cell abundance peaked in winter-spring at high salinity but in summer at low salinity. This was likely due to large filamentous cyanobacteria blooms that occur in summer in low salinity areas, but which are absent in higher salinities. In contrast, phytoplankton local diversity peaked in spring at low salinity but in fall and winter at high salinity. Whilst differences in seasonal variation in cell abundance were reasonably well-explained by variation in salinity and nutrient availability, variation in local-scale phytoplankton diversity was poorly predicted by environmental variables. Overall, we provide insights into the causes of spatio-temporal variation in coastal phytoplankton community structure while also identifying knowledge gaps.

Long-term spatio-temporal changes of the muddy fine sand benthic community of the Bay of Seine (eastern English Channel)

In the English Channel, the eastern Bay of Seine is exposed to numerous anthropogenic disturbances, in particular major changes in sediment dynamics, which are expected to greatly impact benthic communities. To assess the long-term effects of these stressors on the muddy fine sand benthic community, an original long-term monitoring program has been implemented since 1988. It is based on the sampling of a network of 60 stations during seven surveys over 28 years from 1988 to 2016. We investigate changes of species density, species composition and species diversity at different scales (α-diversity, β-diversity and γ-diversity). Contrary to results obtained in many coastal areas, our results showed a long-term persistence of the community in terms of species composition and structure although a general shift towards muddy sediment has resulted in increased colonisation by species associated with muddy habitats and a decrease in spatial beta diversity.

Functional diversity of habitat formers declines scale-dependently across an environmental stress gradient

Marine habitat formers such as seaweeds and corals are lynchpins of coastal ecosystems, but their functional diversity and how it varies with scale and context remains poorly studied. Here, we investigate the functional diversity of seaweed assemblages across the rocky intertidal stress gradient at large (zones) and small (quadrat) scales. We quantified complementary metrics of emergent group richness, functional richness (functional space occupied) and functional dispersion (trait complementarity of dominant species). With increasing shore height, under species loss and turnover, responses of functional diversity were scale- and metric-dependent. At the large scale, functional richness contracted while—notwithstanding a decline in redundancy—emergent group richness and functional dispersion were both invariant. At the small scale, all measures declined, with the strongest responses evident for functional and emergent group richness. Comparisons of observed versus expected values based on null models revealed that functional richness and dispersion were greater than expected in the low shore but converged with expected values higher on the shore. These results show that functional diversity of assemblages of marine habitat formers can be especially responsive to environmental stress gradients at small scales and for richness measures. Furthermore, niche-based processes at the small—neighbourhood—scale can favour co-occurrence of functionally distinctive species under low, but not high, stress, magnifying differences in functional diversity across environmental gradients. As assemblages of marine habitat formers face accelerating environmental change, further studies examining multiple aspects of functional diversity are needed to elucidate patterns, processes, and ecosystem consequences of community (dis-)assembly across diverse groups.

How community adaptation affects biodiversity-ecosystem functioning relationships

Evidence is growing that evolutionary dynamics can impact biodiversity-ecosystem functioning (BEF) relationships. However the nature of such impacts remains poorly understood. Here we use a modelling approach to compare random communities, with no trait evolutionary fine-tuning, and co-adapted communities, where traits have co-evolved, in terms of emerging biodiversity-productivity, biodiversity-stability and biodiversity-invasion relationships. Community adaptation impacted most BEF relationships, sometimes inverting the slope of the relationship compared to random communities. Biodiversity-productivity relationships were generally less positive among co-adapted communities, with reduced contribution of sampling effects. The effect of community-adaptation, though modest regarding invasion resistance, was striking regarding invasion tolerance: co-adapted communities could remain very tolerant to invasions even at high diversity. BEF relationships are thus contingent on the history of ecosystems and their degree of community adaptation. Short-term experiments and observations following recent changes may not be safely extrapolated into the future, once eco-evolutionary feedbacks have taken place.

Biodiversity increases multitrophic energy use efficiency, flow and storage in grasslands

The continuing loss of global biodiversity has raised questions about the risk that species extinctions pose for the functioning of natural ecosystems and the services that they provide for human wellbeing. There is consensus that, on single trophic levels, biodiversity sustains functions; however, to understand the full range of biodiversity effects, a holistic and multitrophic perspective is needed. Here, we apply methods from ecosystem ecology that quantify the structure and dynamics of the trophic network using ecosystem energetics to data from a large grassland biodiversity experiment. We show that higher plant diversity leads to more energy stored, greater energy flow and higher community-energy-use efficiency across the entire trophic network. These effects of biodiversity on energy dynamics were not restricted to only plants but were also expressed by other trophic groups and, to a similar degree, in aboveground and belowground parts of the ecosystem, even though plants are by far the dominating group in the system. The positive effects of biodiversity on one trophic level were not counteracted by the negative effects on adjacent levels. Trophic levels jointly increased the performance of the community, indicating ecosystem-wide multitrophic complementarity, which is potentially an important prerequisite for the provisioning of ecosystem services.

Low redundancy and complementarity shape ecosystem functioning in a low-diversity ecosystem

Ecosystem functioning is positively linked to biodiversity on land and in the sea. In high-diversity systems (e.g. coral reefs), species coexist by sharing resources and providing similar functions at different temporal or spatial scales. How species combine to deliver the ecological function they provide is pivotal for maintaining the structure, functioning and resilience of some ecosystems, but the significance of this is rarely examined in low-diversity systems such as estuaries. We tested whether an ecological function is shaped by biodiversity in a low-diversity ecosystem by measuring the consumption of carrion by estuarine scavengers. Carrion (e.g. decaying animal flesh) is opportunistically fed on by a large number of species across numerous ecosystems. Estuaries were chosen as the model system because carrion consumption is a pivotal ecological function in coastal seascapes, and estuaries are thought to support diverse scavenger assemblages, which are modified by changes in water quality and the urbanization of estuarine shorelines. We used baited underwater video arrays to record scavengers and measure the rate at which carrion was consumed by fish in 39 estuaries across 1,000 km of coastline in eastern Australia. Carrion consumption was positively correlated with the abundance of only one species, yellowfin bream Acanthopagrus australis, which consumed 58% of all deployed carrion. The consumption of carrion by yellowfin bream was greatest in urban estuaries with moderately hardened shorelines (20%-60%) and relatively large subtidal rock bars (>0.1 km² ). Our findings demonstrate that an ecological function can be maintained across estuarine seascapes despite both limited redundancy (i.e. dominated by one species) and complementarity (i.e. there is no spatial context where the function is delivered significantly when yellowfin bream are not present) in the functional traits of animal assemblages. The continued functioning of estuaries, and other low-diversity ecosystems, might therefore not be tightly linked to biodiversity, and we suggest that the preservation of functionally dominant species that maintain functions in these systems could help to improve conservation outcomes for coastal seascapes.

Predator identity dominates non-consumptive effects in a disease-impacted rocky shore food web

Predicting the effects of predator diversity loss on food webs is challenging, because predators can both consume and induce behavioral responses in their prey (i.e., non-consumptive effects or NCEs). Studies manipulating predator diversity and investigating NCEs are rare, especially in marine systems. Recently, a severe outbreak of sea star wasting syndrome (SSWS) on the west coast of North America resulted in unprecedented declines of the sea star Pisaster ochraceus. We investigated the consequences of Pisaster loss on an abundant grazer, the black turban snail Tegula funebralis, through NCEs. We combined a laboratory experiment and field surveys to examine the importance of identity vs. diversity in a predator assemblage (Pisaster, crabs, and octopuses) on Tegula behavior, feeding, and growth. Laboratory and field results indicated that predator identity, not diversity, drives Tegula behavior and causes NCEs. Mesocosm treatments with Pisaster caused greater NCEs on Tegula than assemblages without Pisaster. Tegula's distribution in the field, which is driven primarily by anti-predator behavior, was strongly associated only with Pisaster abundance, and not with the abundance of crabs, octopuses, and other predatory sea stars (Leptasterias spp.). We conclude that Pisaster primarily drives Tegula vertical distribution and may be having strong NCEs on Tegula on northern California rocky shores. Furthermore, predator diversity in northern California does not provide functional redundancy, in terms of NCEs on Tegula, to buffer the system from Pisaster loss. Thus, predator-induced vertical distributions and grazing suppression may not be maintained in areas where Pisaster populations are reduced or slow to recover from SSWS.

Realistic scenarios of environmental disturbance lead to functionally important changes in benthic species-environment interactions

Changes in community structure concurrent with environmental forcing often form a precursor to changes in species diversity, and can have substantive consequences for ecosystem functioning. Here, we assess the effects of altered levels of evenness that are representative of different levels of eutrophication and changes in salinity associated with altered precipitation patterns, on the mediation of nutrient release by sediment-dwelling invertebrate communities. We find that an adjustment towards a more even distribution of species corresponds with an increase in sediment particle reworking that, in general, translates to increased levels of nutrient release. This response, however, is dependent on the functional role of each species in the community and is influenced by concomitant changes in salinity, especially when salinity extends beyond the range typically experienced by the community. Overall, our findings highlight the dynamic nature of species contributions to functioning and reinforce the importance of understanding when, and how, the mechanistic basis of species-environment interactions are modified as the influence of abiotic and biotic factors flex under periods of directional forcing.

Ecological and functional consequences of coastal ocean acidification: Perspectives from the Baltic-Skagerrak System

Ocean temperatures are rising; species are shifting poleward, and pH is falling (ocean acidification, OA). We summarise current understanding of OA in the brackish Baltic-Skagerrak System, focussing on the direct, indirect and interactive effects of OA with other anthropogenic drivers on marine biogeochemistry, organisms and ecosystems. Substantial recent advances reveal a pattern of stronger responses (positive or negative) of species than ecosystems, more positive responses at lower trophic levels and strong indirect interactions in food-webs. Common emergent themes were as follows: OA drives planktonic systems toward the microbial loop, reducing energy transfer to zooplankton and fish; and nutrient/food availability ameliorates negative impacts of OA. We identify several key areas for further research, notably the need for OA-relevant biogeochemical and ecosystem models, and understanding the ecological and evolutionary capacity of Baltic-Skagerrak ecosystems to respond to OA and other anthropogenic drivers.

Biodiversity-ecosystem functioning relationships in fish communities: biomass is related to evenness and the environment, not to species richness

The relationship between biodiversity and ecosystem functioning (BEF) is a topic of considerable interest to scientists and managers because a better understanding of its underlying mechanisms may help us mitigate the consequences of biodiversity loss on ecosystems. Our current knowledge of BEF relies heavily on theoretical and experimental studies, typically conducted on a narrow range of spatio-temporal scales, environmental conditions, and trophic levels. Hence, whether a relationship holds in the natural environment is poorly understood, especially in exploited marine ecosystems. Using large-scale observations of marine fish communities, we applied a structural equation modelling framework to investigate the existence and significance of BEF relationships across northwestern European seas. We find that ecosystem functioning, here represented by spatial patterns in total fish biomass, is unrelated to species richness-the most commonly used diversity metric in BEF studies. Instead, community evenness, differences in species composition, and abiotic variables are significant drivers. In particular, we find that high fish biomass is associated with fish assemblages dominated by a few generalist species of a high trophic level, who are able to exploit both the benthic and pelagic energy pathway. Our study provides a better understanding of the mechanisms behind marine ecosystem functioning and allows for the integration of biodiversity into management considerations.

Biochemical and morphological responses in Ucides cordatus (Crustacea, Decapoda) as indicators of contamination status in mangroves and port areas from northern Brazil

The aims of this study were to analyze biochemical and morphological responses (glutathione S-transferase activity and branchial lesions) in Ucides cordatus (crabs) and to verify how the species is responding to environmental contamination in a port (potentially contaminated area) and mangrove (reference area; Amazon Coast, Maranhão, Brazil). Adult males were captured bimonthly for a period of 1 year. Higher GST activity (1.03 ± 0.07 μmol min^-1 mg protein^-1) was observed in crabs in the port when compared with those in the reference area (p < 0.05). The greatest number of branchial lesions (serious alterations) was recorded in crabs only in the port area. The GST activity increased until serious lesions appeared; after this limit, GST activity decreased dramatically to very low levels, thus resulting in irreversible lesions (lamella collapse). The mathematical model based on the two parameters evaluated in U. cordatus showed that the port area experienced substantial contamination impact, while the mangroves (reference area) presented moderate environmental quality.

The effects of functional diversity on biomass production, variability, and resilience of ecosystem functions in a tritrophic system

Diverse communities can adjust their trait composition to altered environmental conditions, which may strongly influence their dynamics. Previous studies of trait-based models mainly considered only one or two trophic levels, whereas most natural system are at least tritrophic. Therefore, we investigated how the addition of trait variation to each trophic level influences population and community dynamics in a tritrophic model. Examining the phase relationships between species of adjacent trophic levels informs about the strength of top-down or bottom-up control in non-steady-state situations. Phase relationships within a trophic level highlight compensatory dynamical patterns between functionally different species, which are responsible for dampening the community temporal variability. Furthermore, even without trait variation, our tritrophic model always exhibits regions with two alternative states with either weak or strong nutrient exploitation, and correspondingly low or high biomass production at the top level. However, adding trait variation increased the basin of attraction of the high-production state, and decreased the likelihood of a critical transition from the high- to the low-production state with no apparent early warning signals. Hence, our study shows that trait variation enhances resource use efficiency, production, stability, and resilience of entire food webs.

Biodiversity and ecosystem functioning in naturally assembled communities

Approximately 25 years ago, ecologists became increasingly interested in the question of whether ongoing biodiversity loss matters for the functioning of ecosystems. As such, a new ecological subfield on Biodiversity and Ecosystem Functioning (BEF) was born. This subfield was initially dominated by theoretical studies and by experiments in which biodiversity was manipulated, and responses of ecosystem functions such as biomass production, decomposition rates, carbon sequestration, trophic interactions and pollination were assessed. More recently, an increasing number of studies have investigated BEF relationships in non-manipulated ecosystems, but reviews synthesizing our knowledge on the importance of real-world biodiversity are still largely missing. I performed a systematic review in order to assess how biodiversity drives ecosystem functioning in both terrestrial and aquatic, naturally assembled communities, and on how important biodiversity is compared to other factors, including other aspects of community composition and abiotic conditions. The outcomes of 258 published studies, which reported 726 BEF relationships, revealed that in many cases, biodiversity promotes average biomass production and its temporal stability, and pollination success. For decomposition rates and ecosystem multifunctionality, positive effects of biodiversity outnumbered negative effects, but neutral relationships were even more common. Similarly, negative effects of prey biodiversity on pathogen and herbivore damage outnumbered positive effects, but were less common than neutral relationships. Finally, there was no evidence that biodiversity is related to soil carbon storage. Most BEF studies focused on the effects of taxonomic diversity, however, metrics of functional diversity were generally stronger predictors of ecosystem functioning. Furthermore, in most studies, abiotic factors and functional composition (e.g. the presence of a certain functional group) were stronger drivers of ecosystem functioning than biodiversity per se. While experiments suggest that positive biodiversity effects become stronger at larger spatial scales, in naturally assembled communities this idea is too poorly studied to draw general conclusions. In summary, a high biodiversity in naturally assembled communities positively drives various ecosystem functions. At the same time, the strength and direction of these effects vary highly among studies, and factors other than biodiversity can be even more important in driving ecosystem functioning. Thus, to promote those ecosystem functions that underpin human well-being, conservation should not only promote biodiversity per se, but also the abiotic conditions favouring species with suitable trait combinations.

Foundation species patch configuration mediates salt marsh biodiversity, stability and multifunctionality

Foundation species enhance biodiversity and multifunctionality across many systems; however, whether foundation species patch configuration mediates their ecological effects is unknown. In a 6-month field experiment, we test which attributes of foundation species patch configuration - i.e. patch size, total patch area, perimeter, area-perimeter ratio, or connectivity - control biodiversity, stability and multifunctionality by adding a standardised density of mussel foundation species in patches of 1, 5, 10, 30, 60, 90 or 180 individuals to a southeastern US salt marsh. Over 67% of response variables increased with clustering of mussels, responses that were driven by increases in area-perimeter ratio (33%), decreases in perimeter (29%), or increases in patch size (5%), suggesting sensitivity to external stressors and/or dependence on foundation species-derived niche availability and segregation. Thus, mussel configuration - by controlling the relative distribution of multidimensional patch interior and edge niche space - critically modulates this foundation species' effects on ecosystem structure, stability and function.

Long-term nitrogen fertilization decreases bacterial diversity and favors the growth of Actinobacteria and Proteobacteria in agro-ecosystems across the globe

Long-term elevated nitrogen (N) input from anthropogenic sources may cause soil acidification and decrease crop yield, yet the response of the belowground microbial community to long-term N input alone or in combination with phosphorus (P) and potassium (K) is poorly understood. We explored the effect of long-term N and NPK fertilization on soil bacterial diversity and community composition using meta-analysis of a global dataset. Nitrogen fertilization decreased soil pH, and increased soil organic carbon (C) and available N contents. Bacterial taxonomic diversity was decreased by N fertilization alone, but was increased by NPK fertilization. The effect of N fertilization on bacterial diversity varied with soil texture and water management, but was independent of crop type or N application rate. Changes in bacterial diversity were positively related to both soil pH and organic C content under N fertilization alone, but only to soil organic C under NPK fertilization. Microbial biomass C decreased with decreasing bacterial diversity under long-term N fertilization. Nitrogen fertilization increased the relative abundance of Proteobacteria and Actinobacteria, but reduced the abundance of Acidobacteria, consistent with the general life history strategy theory for bacteria. The positive correlation between N application rate and the relative abundance of Actinobacteria indicates that increased N availability favored the growth of Actinobacteria. This first global analysis of long-term N and NPK fertilization that differentially affects bacterial diversity and community composition provides a reference for nutrient management strategies for maintaining belowground microbial diversity in agro-ecosystems worldwide.

Selective and context-dependent effects of chemical stress across trophic levels at the basis of marine food webs

Human activities increasingly impact the functioning of marine food webs, but anthropogenic stressors are seldom included in ecological study designs. Diet quality, as distinct from just diet quantity, has moreover rarely been highlighted in food web studies in a stress context. We measured the effects of metal and pesticide stress (copper and atrazine) on the contribution of a benthic intertidal diatom community to two processes that are key to the functioning of intertidal systems: biomass (diet quantity) and lipid (diet quality) production. We then examined if stressors affected diatom functioning by selectively targeting the species contributing most to functioning (selective stress effects) or by changing the species' functional contribution (context-dependent effects). Finally, we tested if stress-induced changes in diet quality altered the energy flow to the diatoms' main grazers (harpacticoid copepods). Diatom diet quantity was reduced by metal stress but not by low pesticide levels due to the presence of an atrazine-tolerant, mixotrophic species. Selective effects of the pesticide reduced diatom diet quality by 60% and 75% at low and high pesticide levels respectively, by shifting diatom community structure from dominance by lipid-rich species toward dominance by an atrazine-tolerant, but lipid-poor, species. Context-dependent effects did not affect individual diatom lipid content at low levels of both stressors, but caused diatoms to lose 40% of their lipids at high copper stress. Stress-induced changes in diet quality predicted the energy flow from the diatoms to their copepod consumers, which lost half of their lipids when feeding on diatoms grown under low and high pesticide and high metal stress. Selective pesticide effects were a more important threat for trophic energy transfer than context-dependent effects of both stressors, with shifts in diatom community structure affecting the energy flow to their copepod grazers at stress levels where no changes in diatom lipid content were detected.

Increasing densities of an invasive polychaete enhance bioturbation with variable effects on solute fluxes

Bioturbation is a key process affecting nutrient cycling in soft sediments. The invasive polychaete genus Marenzelleria spp. has established successfully throughout the Baltic Sea increasing species and functional diversity with possible density-dependent effects on bioturbation and associated solute fluxes. We tested the effects of increasing density of M. arctia, M. viridis and M. neglecta on bioturbation and solute fluxes in a laboratory experiment. Benthic communities in intact sediment cores were manipulated by adding increasing numbers of Marenzelleria spp. The results showed that Marenzelleria spp. in general enhanced all bioturbation metrics, but the effects on solute fluxes varied depending on the solute, on the density and species identity of Marenzelleria, and on the species and functional composition of the surrounding community. M. viridis and M. neglecta were more important in predicting variation in phosphate and silicate fluxes, whereas M. arctia had a larger effect on nitrogen cycling. The complex direct and indirect pathways indicate the importance of considering the whole community and not just species in isolation in the experimental studies. Including these interactions provides a way forward regarding our understanding of the complex ecosystem effects of invasive species.

Bioturbation activity of three macrofaunal species and the presence of meiofauna affect the abundance and composition of benthic bacterial communities

Given concerns of increasing rates of species extinctions, the relationship between biodiversity and ecosystem functioning has become a major research focus over the past two decades. Many studies have shown that biodiversity per se (e.g. species richness) or species-specific traits may be good predictors of changes in ecosystem function. Although numerous studies on this subject have focused on terrestrial systems, few have evaluated benthic marine systems. We used the Limecola balthica community as a model to test whether the number or identity of three well-studied macrofaunal species influence the sediment bacterial compartment, which drives important biogeochemical processes and influence ecosystem functioning. We also investigated the poorly known role of meiofauna in the interactions between macrofauna and bacteria. Eight combinations of 0-3 species were maintained in microcosms for 34 days in the presence or absence of meiofauna. The abundance and composition of the bacterial community, defined by the relative percentage of cells with a high (HNA) vs low (LNA) nucleic acid content, were measured. Species identity of macrofauna was a better indicator of changes in the microbial compartment than was species richness per se. In particular, the gallery-diffuser behaviour of the polychaete Alitta virens likely induced strong changes in sediment physical and geochemical properties with a major impact on the bacterial compartment. Moreover, the presence of meiofauna modulated the influence of macrofauna on bacterial communities. This study provides evidence that species identity provides greater explanatory power than species richness to predict changes in the bacterial compartment. We propose that multi-compartment approaches to describe interactions amongst different size classes of organisms and their ecological roles should be further developed to improve our understanding of benthic ecosystem functioning.

Predator and prey biodiversity relationship and its consequences on marine ecosystem functioning-interplay between nanoflagellates and bacterioplankton

The importance of biodiversity effects on ecosystem functioning across trophic levels, especially via predatory-prey interactions, is receiving increased recognition. However, this topic has rarely been explored for marine microbes, even though microbial biodiversity contributes significantly to marine ecosystem function and energy flows. Here we examined diversity and biomass of bacteria (prey) and nanoflagellates (predators), as well as their effects on trophic transfer efficiency in the East China Sea. Specifically, we investigated: (i) predator diversity effects on prey biomass and trophic transfer efficiency (using the biomass ratio of predator/prey as a proxy), (ii) prey diversity effects on predator biomass and trophic transfer efficiency, and (iii) the relationship between predator and prey diversity. We found higher prey diversity enhanced both diversity and biomass of predators, as well as trophic transfer efficiency, which may arise from more balanced diet and/or enhanced niche complementarity owing to higher prey diversity. By contrast, no clear effect was detected for predator diversity on prey biomass and transfer efficiency. Notably, we found prey diversity effects on predator-prey interactions; whereas, we found no significant diversity effect on biomass within the same trophic level. Our findings highlight the importance of considering multi-trophic biodiversity effects on ecosystem functioning in natural ecosystems.

Energy Flux: The Link between Multitrophic Biodiversity and Ecosystem Functioning

Relating biodiversity to ecosystem functioning in natural communities has become a paramount challenge as links between trophic complexity and multiple ecosystem functions become increasingly apparent. Yet, there is still no generalised approach to address such complexity in biodiversity-ecosystem functioning (BEF) studies. Energy flux dynamics in ecological networks provide the theoretical underpinning of multitrophic BEF relationships. Accordingly, we propose the quantification of energy fluxes in food webs as a powerful, universal tool for understanding ecosystem functioning in multitrophic systems spanning different ecological scales. Although the concept of energy flux in food webs is not novel, its application to BEF research remains virtually untapped, providing a framework to foster new discoveries into the determinants of ecosystem functioning in complex systems.