Reconciling complexity with stability in naturally assembling food webs

Newcomb-Benford number law and ecological processes

The Newcomb-Benford number law has found applications in the natural and social sciences for decades, with limited ecological attention. The aim of this communication is to highlight a statistical correspondence between the first significant digit frequencies of the Benford probability distribution and ecological systems in a balanced state of dynamic equilibrium. Analytical methods, including multidimensional Euclidean distance, Cohen-W effect size, Kossovsky sum of squared deviations (SSD), and Pearson residuals are presented to facilitate the identification of this canonical representation across multiple levels of ecological organization and scale. Case studies reveal novel applications of the Benford distribution for detecting impending state transitions in marginally stable systems, as well as temporal and spatial divergence of ecological information through the measurement of Kullback-Leibler relative entropy. Widespread documentation of the leading digit phenomenon is expected as ecologists revisit empirical datasets and formalize sampling protocols for its detection. The conversion of randomly collected sets of arithmetic data into logarithmic probabilities of first significant digits presents unique opportunities to advance our understanding of ecological processes related with stability, complexity, and maturity.

Soil biodiversity and function under global change

Soil organisms represent the most abundant and diverse organisms on the planet and support almost every ecosystem function we know, and thus impact our daily lives. Some of these impacts have been well-documented, such as the role of soil organisms in regulating soil fertility and carbon sequestration; processes that have direct implications for essential ecosystem services including food security and climate change mitigation. Moreover, soil biodiversity also plays a critical role in supporting other aspects from One Health-the combined health of humans, animals, and the environment-to the conservation of historic structures such as monuments. Unfortunately, soil biodiversity is also highly vulnerable to a growing number of stressors associated with global environmental change. Understanding how and when soil biodiversity supports these functions, and how it will adapt to changing environmental conditions, is crucial for conserving soils and maintaining soil processes for future generations. In this Essay, we discuss the fundamental importance of soil biodiversity for supporting multiple ecosystem services and One Health, and further highlight essential knowledge gaps that need to be addressed to conserve soil biodiversity for the next generations.

Impacts of trophic interactions on carbon accrual in soils

The transformation and stabilization of soil organic carbon (SOC) are important processes of global carbon (C) cycling, with implications for climate change. Much attention has been given to microbial anabolic processes driving SOC accrual. These are referred to as the soil microbial carbon pump (MCP), which emphasizes the contribution of microbial metabolism and necromass to the stable soil C pool. However, we still lack a fundamental understanding of how trophic interactions between soil fauna and microbiota modulate microbial necromass production and, consequently, SOC formation. Here, we provide an ecological perspective on the impacts of trophic interactions on modulating necromass formation and C accrual in soils. We discuss the mechanisms of trophic interactions in the context of food web ecology, with a focus on trophic control of microbial population densities and their influences on soil microbiota assembly. We foresee that integrating trophic interactions into the soil MCP framework can provide a more comprehensive basis for guiding future research efforts to elucidate the mechanisms modulating microbial necromass and SOC formation in terrestrial ecosystems. This perspective offers an ecological foundation for leveraging the use of biological interventions to enhance SOC accrual, providing valuable insights for sustainable C management strategies.

Diversity, connectivity and negative interactions define robust microbiome networks across land, stream, and sea

In this era of rapid global change, factors influencing the stability of ecosystems and their functions have come into the spotlight. For decades the relationship between stability and complexity has been investigated in modeled and empirical systems, yet results remain largely context dependent. To overcome this we leverage a multiscale inventory of fungi and bacteria ranging from single sites along an environmental gradient, to habitats inclusive of land, sea and stream, to an entire watershed. We use networks to assess the relationship between microbiome complexity and robustness and identify fundamental principles of stability. We demonstrate that while some facets of complexity are positively associated with robustness, others are not. Beyond positive biodiversity x robustness relationships we find that the number of "gatekeeper" species or those that are highly connected and central within their networks, and the proportion of predicted negative interactions are universal indicators of robust microbiomes. With the potential promise of microbiome engineering to address global challenges ranging from human to ecosystem health we identify properties of microbiomes for future experimental studies that may enhance their stability. We emphasize that features beyond biodiversity and additional characteristics beyond stability such as adaptability should be considered in these efforts.

Seasonal Shifts in Trophic Interaction Strength Drive Stability of Natural Food Webs

It remains challenging to understand why natural food webs are remarkably stable despite highly variable environmental factors and population densities. We investigated the dynamics in the structure and stability of Lake Constance's pelagic food web using 7 years of high-frequency observations of biomasses and production, leading to 59 seasonally resolved quantitative food web descriptions. We assessed the dynamics in asymptotic food web stability through maximum loop weight, which revealed underlying stability mechanisms. Maximum loop weight showed a recurrent seasonal pattern with a consistently high stability despite pronounced dynamics in biomasses, fluxes and productivity. This stability resulted from seasonal rewiring of the food web, driven by energetic constraints within loops and their embedding into food web structure. The stabilising restructuring emerged from counter-acting effects of metabolic activity and competitiveness/susceptibility to predation within a diverse grazer community on loop weight. This underscores the role of functional diversity in promoting food web stability.

Differential impacts of small hydropower plants on macroinvertebrate communities upstream and downstream under ecological flow

Hydropower dams influence freshwater biodiversity by altering river flow patterns and habitat conditions. With the global surge in small hydropower plants (SHPs), their impacts on aquatic ecosystems have become increasingly recognized. However, most previous studies did not consider the recently implemented ecological flows. Consequently, the effects of SHPs under ecological flow conditions on aquatic organisms, such as macroinvertebrate communities, remain unclear. We surveyed 15 SHPs in the Oujiang region, establishing sampling sites upstream of the intake dams (S1), in dam-induced reservoirs (S2), in dewatered sections downstream of the dams with ecological flows (S3), and in sections with restored natural flow (S4). By comparing macroinvertebrate community composition, diversity, functional feeding groups, and network structures in these areas, we assessed the ecological response of macroinvertebrates to SHPs under ecological flows. Our research found that SHPs significantly impact macroinvertebrate communities. Specifically, at site S2, stagnant water species replaced those typically found in flowing conditions, resulting in a marked difference in species composition between S2 and other sites. Compared to S1 and S4, diversity indices at S2 and S3 were lower, with filterers and collectors dominating the functional feeding groups at S2 and S3. Co-occurrence network analysis revealed that network complexity at S2 and S3 was lower than at S1 and S4. Additionally, S3 was less affected by SHPs than S2, underscoring the importance of ecological flow replenishment. Overall, our research confirmed the remarkable influence of SHPs on S2 macroinvertebrate community, and emphasized the importance of maintaining sufficient ecological flow to the downstream aquatic organism of S3 reach. We suggest a comprehensive assessment of the potential environmental impacts of SHPs, particularly the negative effects caused by insufficient ecological flow, to ensure the sustainable development of ecosystems.

Stability of Ecological Systems: A Theoretical Review

The stability of ecological systems is a fundamental concept in ecology, which offers profound insights into species coexistence, biodiversity, and community persistence. In this article, we provide a systematic and comprehensive review on the theoretical frameworks for analyzing the stability of ecological systems. Notably, we survey various stability notions, including linear stability, sign stability, diagonal stability, D-stability, total stability, sector stability, and structural stability. For each of these stability notions, we examine necessary or sufficient conditions for achieving such stability and demonstrate the intricate interplay of these conditions on the network structures of ecological systems. We further discuss the stability of ecological systems with higher-order interactions.

The disruption of birds' double mutualistic interactions in novel ecosystems

Non-native trees disrupt ecological processes vital to native plant communities. We studied how forests dominated by Acacia dealbata and Eucalyptus globulus affect the role of birds as dual pollinators and seed dispersers in a region heavily impacted by these two non-native species. We compared bird-plant interactions in the native and in the two non-native forest types. We constructed a multilayer regional network for each forest type and evaluated differences in network dissimilarity between networks. We also calculated the bird's importance in connecting processes and variables associated with module diversity. To determine how the networks react to changes in species richness, we did a simulation of species richness gradient and link percentage for each forest type. The number of birds acting both as pollinators and seed dispersers was higher in native than in non-native forests. However, birds in non-native forests still play a crucial role in maintaining the ecological services provided to native plant communities. However, the eucalyptus network exhibited a concerning simplification, forcing bird species to fully exploit the few remaining resources, leaving little room for structural adjustments and limiting the ecosystem's ability to withstand further species loss. These findings highlight how non-native trees may trigger cascading effects across trophic levels.

Spartina alterniflora invasion reduces soil microbial diversity and weakens soil microbial inter-species relationships in coastal wetlands

Soil microorganisms play a crucial role in the plant invasion process, acting as both drivers of and responders to plant invasion. However, the effects of plant invasion on the complexity and stability of co-occurrence networks of soil microbial communities remain unclear. Here, we investigated how the invasion of Spartina alterniflora affected the diversity, composition, and co-occurrence networks of soil bacterial and fungal communities in the Yellow River Delta, China. Compared to the native plant (Suaeda salsa), S. alterniflora invasion decreased the α-diversity of soil bacterial communities but did not affect that of fungal communities. The β-diversity of soil bacterial and fungal communities under S. salsa and S. alterniflora habitats also differed dramatically. S. alterniflora invasion increased the relative abundance of the copiotrophic phylum Bacteroidota, whereas decreased the relative abundances of the oligotrophic phyla Acidobacteriota and Gemmatimonadota. Additionally, the relative abundance of Chytridiomycota, known for its role in degrading recalcitrant organic matter, increased substantially within the soil fungal community. Functional predictions revealed that S. alterniflora invasion increased the relative abundance of certain soil bacteria involved in carbon and nitrogen cycling, including aerobic chemoheterotrophy, nitrate reduction, and nitrate respiration. More importantly, S. alterniflora invasion reduced the complexity and stability of both soil bacterial and fungal community networks. The shifts in soil microbial community structure and diversity were mainly induced by soil available nutrients and soil salinity. Overall, our study highlights the profound impacts of S. alterniflora invasion on soil microbial communities, which could further indicate the modification of ecosystem functioning by invasive species.

Putting the Asymmetric Response Concept to the test: Modeling multiple stressor exposure and release in a stream food web

Communities in stream ecosystems often respond asymmetrically to increase and release of stressors, as indicated by slow and incomplete recovery. The Asymmetric Response Concept (ARC) posits that this is due to a shift in the relative importance of three mechanisms: tolerance, dispersal, and biotic interactions. In complex natural communities, these mechanisms may produce alternative outcomes through poorly understood indirect effects. To understand how the three mechanisms respond to different temporal stressor scenarios, we studied multiple scenarios using a stream food web model. We asked the following questions: Do groups of species decline as expected on the basis of individual tolerance rankings derived from laboratory experiments when they are embedded in a complex dynamic food web? Does the response of ecosystem function match that of communities? To address these questions, we aggregated data on individual tolerances at the level of functional groups and studied how single and multiple stressors affect food web dynamics and nutrient cycling. Multiple stressor scenarios involved different intensities of salt and temperature increase. Functional groups exhibited a different relative tolerance ranking between the laboratory and dynamic food web contexts. Salt as a single stressor had only minor and transient effects at low level but led to the loss of one or more functional groups at high level. In contrast, high temperature, alone or in combination with salt, caused the loss of functional groups at all tested levels. Patterns often differed between the response of communities and ecosystem function. We discuss our findings with respect to the ARC.

A successional shift enhances stability in ant symbiont communities

Throughout succession, communities undergo structural shifts, which can alter the relative abundances of species and how they interact. It is frequently asserted that these alterations beget stability, i.e. that succession selects for communities better able to resist perturbations. Yet, whether and how alterations of network structure affect stability during succession in complex communities is rarely studied in natural ecosystems. Here, we explore how network attributes influence stability of different successional stages of a natural network: symbiotic arthropod communities forming food webs inside red wood ant nests. We determined the abundance of 16 functional groups within the symbiont community across 51 host nests in the beginning and end stages of succession. Nest age was the main driver of the compositional shifts: symbiont communities in old nests contained more even species abundance distributions and a greater proportion of specialists. Based on the abundance data, we reconstructed interaction matrices and food webs of the symbiont community for each nest. We showed that the enhanced community evenness in old nests leads to an augmented food web stability in all but the largest symbiont communities. Overall, this study demonstrates that succession begets stability in a natural ecological network by making the community more even.

Multitrophic Higher-Order Interactions Modulate Species Persistence

AbstractEcologists increasingly recognize that interactions between two species can be affected by the density of a third species. How these higher-order interactions (HOIs) affect species persistence remains poorly understood. To explore the effect of HOIs stemming from multiple trophic layers on a plant community composition, we experimentally built a mesocosm with three plants and three pollinator species arranged in a fully nested and modified network structure. We estimated pairwise interactions among plants and between plants and pollinators, as well as HOIs initiated by a plant or a pollinator affecting plant species pairs. Using a structuralist approach, we evaluated the consequences of the statistically supported HOIs on the persistence probability of each of the three competing plant species and their combinations. HOIs substantially redistribute the strength and sign of pairwise interactions between plant species, promoting the opportunities for multispecies communities to persist compared with a non-HOI scenario. However, the physical elimination of a plant-pollinator link in the modified network structure promotes changes in per capita pairwise interactions and HOIs, resulting in a single-species community. Our study provides empirical evidence of the joint importance of HOIs and network structure in determining species persistence within diverse communities.

Trophic tug-of-war: Coexistence mechanisms within and across trophic levels

Ecological communities encompass rich diversity across multiple trophic levels. While modern coexistence theory has been widely applied to understand community assembly, its traditional formalism only allows assembly within a single trophic level. Here, using an expanded definition of niche and fitness differences applicable to multitrophic communities, we study how diversity within and across trophic levels affects species coexistence. If each trophic level is analysed separately, both lower- and higher trophic levels are governed by the same coexistence mechanisms. In contrast, if the multitrophic community is analysed as a whole, different trophic levels are governed by different coexistence mechanisms: coexistence at lower trophic levels is predominantly limited by fitness differences, whereas coexistence at higher trophic levels is predominantly limited by niche differences. This dichotomy in coexistence mechanisms is supported by theoretical derivations, simulations of phenomenological and trait-based models, and a case study of a primeval forest ecosystem. Our work provides a general and testable prediction of coexistence mechanism operating in multitrophic communities.

Unlocking ensemble ecosystem modelling for large and complex networks

The potential effects of conservation actions on threatened species can be predicted using ensemble ecosystem models by forecasting populations with and without intervention. These model ensembles commonly assume stable coexistence of species in the absence of available data. However, existing ensemble-generation methods become computationally inefficient as the size of the ecosystem network increases, preventing larger networks from being studied. We present a novel sequential Monte Carlo sampling approach for ensemble generation that is orders of magnitude faster than existing approaches. We demonstrate that the methods produce equivalent parameter inferences, model predictions, and tightly constrained parameter combinations using a novel sensitivity analysis method. For one case study, we demonstrate a speed-up from 108 days to 6 hours, while maintaining equivalent ensembles. Additionally, we demonstrate how to identify the parameter combinations that strongly drive feasibility and stability, drawing ecological insight from the ensembles. Now, for the first time, larger and more realistic networks can be practically simulated and analysed.

Systems theory, thermodynamics and life: Integrated thinking across ecology, organization and biological evolution

In this paper we explore the relevance and integration of system theory and thermodynamics in terms of the Earth system. It is proposed that together, these fields explain the evolution, organization, functionality and directionality of life on Earth. We begin by summarizing historical and current thinking on the definition of life itself. We then investigate the evidence for a single unit of life. Given that any definition of life and its levels of organization are intertwined, we explore how the Earth system is structured and functions from an energetic perspective, by outlining relevant thermodynamic theory relating to molecular, metabolic, cellular, individual, population, species, ecosystem and biome organization. We next investigate the fundamental relationships between systems theory and thermodynamics in terms of the Earth system, examining the key characteristics of self-assembly, self-organization (including autonomy), emergence, non-linearity, feedback and sub-optimality. Finally, we examine the relevance of systems theory and thermodynamics with reference to two specific aspects: the tempo and directionality of evolution and the directional and predictable process of ecological succession. We discuss the importance of the entropic drive in understanding altruism, multicellularity, mutualistic and antagonistic relationships and how maximum entropy production theory may explain patterns thought to evidence the intermediate disturbance hypothesis.

Shorebirds-driven trophic cascade helps restore coastal wetland multifunctionality

Ecosystem restoration has traditionally focused on re-establishing vegetation and other foundation species at basal trophic levels, with mixed outcomes. Here, we show that threatened shorebirds could be important to restoring coastal wetland multifunctionality. We carried out surveys and manipulative field experiments in a region along the Yellow Sea affected by the invasive cordgrass Spartina alterniflora. We found that planting native plants alone failed to restore wetland multifunctionality in a field restoration experiment. Shorebird exclusion weakened wetland multifunctionality, whereas mimicking higher predation before shorebird population declines by excluding their key prey - crab grazers - enhanced wetland multifunctionality. The mechanism underlying these effects is a simple trophic cascade, whereby shorebirds control crab grazers that otherwise suppress native vegetation recovery and destabilize sediments (via bioturbation). Our findings suggest that harnessing the top-down effects of shorebirds - through habitat conservation, rewilding, or temporary simulation of consumptive or non-consumptive effects - should be explored as a nature-based solution to restoring the multifunctionality of degraded coastal wetlands.

Competitive hierarchies in bryozoan assemblages mitigate network instability by keeping short and long feedback loops weak

Competitive hierarchies in diverse ecological communities have long been thought to lead to instability and prevent coexistence. However, system stability has never been tested, and the relation between hierarchy and instability has never been explained in complex competition networks parameterised with data from direct observation. Here we test model stability of 30 multispecies bryozoan assemblages, using estimates of energy loss from observed interference competition to parameterise both the inter- and intraspecific interactions in the competition networks. We find that all competition networks are unstable. However, instability is mitigated considerably by asymmetries in the energy loss rates brought about by hierarchies of strong and weak competitors. This asymmetric organisation results in asymmetries in the interaction strengths, which reduces instability by keeping the weight of short (positive) and longer (positive and negative) feedback loops low. Our results support the idea that interference competition leads to instability and exclusion but demonstrate that this is not because of, but despite, competitive hierarchy.

Temperature, productivity, and habitat characteristics collectively drive lake food web structure

While many efforts have been devoted to understand variations in food web structure among terrestrial and aquatic ecosystems, the environmental factors influencing food web structure at large spatial scales remain hardly explored. Here, we compiled biodiversity inventories to infer food web structure of 67 French lakes using an allometric niche-based model and tested how environmental variables (temperature, productivity, and habitat) influence them. By applying a multivariate analysis on 20 metrics of food web topology, we found that food web structural variations are represented by two distinct complementary and independent structural descriptors. The first is related to the overall trophic diversity, whereas the second is related to the vertical structure. Interestingly, the trophic diversity descriptor was mostly explained by habitat size (26.7% of total deviance explained) and habitat complexity (20.1%) followed by productivity (dissolved organic carbon: 16.4%; nitrate: 9.1%) and thermal variations (10.7%). Regarding the vertical structure descriptor, it was mostly explained by water thermal seasonality (39.0% of total deviance explained) and habitat depth (31.9%) followed by habitat complexity (8.5%) and size (5.5%) as well as annual mean temperature (5.6%). Overall, we found that temperature, productivity, and habitat characteristics collectively shape lake food web structure. We also found that intermediate levels of productivity, high levels of temperature (mean and seasonality), as well as large habitats are associated with the largest and most complex food webs. Our findings, therefore, highlight the importance of focusing on these three components especially in the context of global change, as significant structural changes in aquatic food webs could be expected under increased temperature, pollution, and habitat alterations.

Disentangling the effects of phosphorus loading on food web stability in a large shallow lake

Excessive nutrient loads reduce ecosystem resilience, resulting in fundamental changes in ecosystem structure and function when exceeding a certain threshold. However, quantitative analysis of the processes by which nutrient loading affects ecosystem resilience requires further exploration. Food web stability is at the heart of ecosystem resilience. In this study, we simulated the dynamics of the food web under different phosphorus loads for Lake Baiyangdian using the PCLake model and calculated the food web stability. Our results showed that there was a good correspondence between the food web stability and ecosystem state response to phosphorus loads. This relationship confirmed that food web stability could be regarded as a signal for the state transition in a real lake ecosystem. Moreover, our estimates suggested that food web stability was influenced only by several functional groups and their interaction strength. Diatoms and zooplankton were the key functional groups that affected food web stability. Phosphorus loads alter the distribution of functional group biomass, which in turn affects energy delivery and, ultimately, the stability of the food web. Corresponding to functional groups, the interactions among zooplankton, diatoms and detritus had the greatest impact, and the interaction strength of the three was positively correlated with food web stability. Overall, our study explained that food-web stability was critical to characterize ecosystem resilience response to external disturbances and can be turned into a scientific tool for lake ecosystem management.

Impacts of extreme climatic events on trophic network complexity and multidimensional stability

Untangling the relationship between network complexity and ecological stability under climate change is an arduous challenge for theoretical and empirical ecology. Even more so, when considering extreme climatic events. Here, we studied the effects of extreme climatic events (heatwaves) on the complexity of realistic freshwater ecosystems using topological and quantitative trophic network metrics. Next, we linked changes in network complexity with the investigation of four stability components (temporal stability, resistance, resilience, and recovery) of community's functional, compositional, and energy flux stability. We found reduction in topological network complexity to be correlated with reduction of functional and compositional resistance. However, temperature-driven increase in link-weighted network complexity increased functional and energy flux recovery and resilience, but at the cost of increased compositional instability. Overall, we propose an overarching approach to elucidate the effects of climate change on multidimensional stability through the lens of network complexity, providing helpful insights for preserving ecosystems stability under climate change.

Stabilization through self-coupling in networks of small-world and scale-free topology

Mechanisms that ensure the stability of dynamical systems are of vital importance, in particular in our globalized and increasingly interconnected world. The so-called connectivity-stability dilemma denotes the theoretical finding that increased connectivity between the components of a large dynamical system drastically reduces its stability. This result has promoted controversies within ecology and other fields of biology, especially, because organisms as well as ecosystems constitute systems that are both highly connected and stable. Hence, it has been a major challenge to find ways to stabilize complex systems while preserving high connectivity at the same time. Investigating the stability of networks that exhibit small-world or scale-free topology is of particular interest, since these topologies have been found in many different types of real-world networks. Here, we use an approach to stabilize recurrent networks of small-world and scale-free topology by increasing the average self-coupling strength of the units of a network. For both topologies, we find that there is a sharp transition from instability to asymptotic stability. Then, most importantly, we find that the average self-coupling strength needed to stabilize a system increases much slower than its size. It appears that the qualitative shape of this relationship is the same for small-world and scale-free networks, while scale-free networks can require higher magnitudes of self-coupling. We further explore the stabilization of networks with Kronecker-Leskovec topology. Finally, we argue that our findings, in particular the stabilization of large recurrent networks through small increases in the unit self-regulation, are of practical importance for the stabilization of diverse types of complex systems.

Life-history omnivory in the fairy shrimp Branchinecta orientalis (Branchiopoda: Anostraca)

Very little is known about the feeding of naupliar and juvenile life stages of omnivorous fairy shrimps (Crustacea: Anostraca). Here, we aim to reveal whether the fairy shrimp Branchinecta orientalis is an ontogenetic omnivore and at which age and ontogenetic stage they gain the ability to feed on zooplankton. We assess how food uptake rates change with age until reaching maturity by providing algae (pico- and nanoplanktonic unicellular algae) and zooplankton (rotifers and copepod nauplii) as food in individual experiments. We found that the fairy shrimp B. orientalis started to feed on both types of algal prey immediately after hatching. Nanoplanktonic algae likely represented the most important food source until reaching maturity. Moreover, fairy shrimps started to feed on zooplankton already when they were 7 days old. Slow-moving rotifers gradually gained importance in the fairy shrimp diet with time. Our results reveal an ontogenetic change in the prey spectrum of fairy shrimp. The systematic shift towards omnivory likely affects both phyto- and zooplankton community composition, possibly contributing to temporal changes in food web dynamics in fairy shrimp habitats, and temporary ponds, which may warrant more detailed investigations in future studies.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10750-022-05132-z.

Dynamic trophic shifts in bacterial and eukaryotic communities during the first 30 years of microbial succession following retreat of an Antarctic glacier

We examined microbial succession along a glacier forefront in the Antarctic Peninsula representing ∼30 years of deglaciation to contrast bacterial and eukaryotic successional dynamics and abiotic drivers of community assembly using sequencing and soil properties. Microbial communities changed most rapidly early along the chronosequence, and co-occurrence network analysis showed the most complex topology at the earliest stage. Initial microbial communities were dominated by microorganisms derived from the glacial environment, whereas later stages hosted a mixed community of taxa associated with soils. Eukaryotes became increasingly dominated by Cercozoa, particularly Vampyrellidae, indicating a previously unappreciated role for cercozoan predators during early stages of primary succession. Chlorophytes and Charophytes (rather than cyanobacteria) were the dominant primary producers and there was a spatio-temporal sequence in which major groups became abundant succeeding from simple ice Chlorophytes to Ochrophytes and Bryophytes. Time since deglaciation and pH were the main abiotic drivers structuring both bacterial and eukaryotic communities. Determinism was the dominant assembly mechanism for Bacteria, while the balance between stochastic/deterministic processes in eukaryotes varied along the distance from the glacier front. This study provides new insights into the unexpected dynamic changes and interactions across multiple trophic groups during primary succession in a rapidly changing polar ecosystem.

Characterization of soil microbial community activity and structure for reducing available Cd by rice straw biochar and Bacillus cereus RC-1

The combination of biochar and specific bacteria has been widely applied to remediate Cadmium-contaminated soil. But little is known about how such composites affect the dynamic distribution of metal fractions. This process is accompanied by the alternations of soil properties and microbial community structures. Composite of rice straw biochar and Bacillus cereus RC-1 were applied to investigate its impacts on Cd alleviation and soil microbial diversity and structure. The bacterial/biochar composite treatment decreased the fraction of HOAc-extractable Cd by 38.82%, and increased residual Cd by 23.95% compared to the untreated control. Moreover, compared with the untreated control, the composite treatment significantly increased the soil pH by about 1.5 units, and the activities of catalase, urease and invertase enzymes were increased by 42.39%, 30.50% and 31.20%, respectively. Composite treatment increased soil bacterial and fungal alpha diversity, the relative abundance of Bacillus, Streptomyces, Arthrobacter, and Aspergillus species were also increased. Mantel test and correlation analysis indicated that the effects associated with fungal communities in influencing soil properties were lower than that those of bacterial communities by different treatment. Aggregated boosted tree (ABT) models analysis showed that soil chemical proprieties (as determined by SOM, CEC, AN, etc.,) contributed over 50% of the changes in bacterial and fungal communities by the composite treatment. The co-occurrence network results showed that all treatments enhanced the correlation between OUT groups and improved the possible relationships in the bacterial and fungal communities, especially the interrelationships between bacteria and fungi after the Cd fractions stabilized. These findings provide a new insight to optimal strategies for the remediation of Cd-contaminated soil.

Linking animal behavior to ecosystem change in disturbed environments

Environmental disturbances often cause individuals to change their behavior. The behavioral responses can induce a chain of reactions through the network of species interactions, via consumptive and trait mediated connections. Given that species interactions define ecosystem structure and functioning, changes to these interactions often have ecological repercussions. Here, we explore the transmission of behavioral responses through the network of species interactions, and how the responses influence ecological conditions. We describe the underlying mechanisms and the ultimate impact that the behavioral responses can have on ecosystem structure and functioning, including biodiversity and ecosystems stability and services. We explain why behavioral responses of some species have a larger impact than that of others on ecosystems, and why research should focus on these species and their interactions. With the work, we synthesize existing theory and empirical evidence to provide a conceptual framework that links behavior responses to altered species interactions, community dynamics, and ecosystem processes. Considering that species interactions link biodiversity to ecosystem functioning, a deeper understanding of behavioral responses and their causes and consequences can improve our knowledge of the mechanisms and pathways through which human activities alter ecosystems. This knowledge can improve our ability to predict the effects of ongoing disturbances on communities and ecosystems and decide on the interventions needed to mitigate negative effects.

Bridging Theories for Ecosystem Stability Through Structural Sensitivity Analysis of Ecological Models in Equilibrium

Ecologists are challenged by the need to bridge and synthesize different approaches and theories to obtain a coherent understanding of ecosystems in a changing world. Both food web theory and regime shift theory shine light on mechanisms that confer stability to ecosystems, but from different angles. Empirical food web models are developed to analyze how equilibria in real multi-trophic ecosystems are shaped by species interactions, and often include linear functional response terms for simple estimation of interaction strengths from observations. Models of regime shifts focus on qualitative changes of equilibrium points in a slowly changing environment, and typically include non-linear functional response terms. Currently, it is unclear how the stability of an empirical food web model, expressed as the rate of system recovery after a small perturbation, relates to the vulnerability of the ecosystem to collapse. Here, we conduct structural sensitivity analyses of classical consumer-resource models in equilibrium along an environmental gradient. Specifically, we change non-proportional interaction terms into proportional ones, while maintaining the equilibrium biomass densities and material flux rates, to analyze how alternative model formulations shape the stability properties of the equilibria. The results reveal no consistent relationship between the stability of the original models and the proportionalized versions, even though they describe the same biomass values and material flows. We use these findings to critically discuss whether stability analysis of observed equilibria by empirical food web models can provide insight into regime shift dynamics, and highlight the challenge of bridging alternative modelling approaches in ecology and beyond.

Complexity-stability trade-off in empirical microbial ecosystems

May's stability theory, which holds that large ecosystems can be stable up to a critical level of complexity, a product of the number of resident species and the intensity of their interactions, has been a central paradigm in theoretical ecology. So far, however, empirically demonstrating this theory in real ecological systems has been a long-standing challenge with inconsistent results. Especially, it is unknown whether this theory is pertinent in the rich and complex communities of natural microbiomes, mainly due to the challenge of reliably reconstructing such large ecological interaction networks. Here we introduce a computational framework for estimating an ecosystem's complexity without relying on a priori knowledge of its underlying interaction network. By applying this method to human-associated microbial communities from different body sites and sponge-associated microbial communities from different geographical locations, we found that in both cases the communities display a pronounced trade-off between the number of species and their effective connectance. These results suggest that natural microbiomes are shaped by stability constraints, which limit their complexity.

From diversity to complexity: Microbial networks in soils

Network analysis has been used for many years in ecological research to analyze organismal associations, for example in food webs, plant-plant or plant-animal interactions. Although network analysis is widely applied in microbial ecology, only recently has it entered the realms of soil microbial ecology, shown by a rapid rise in studies applying co-occurrence analysis to soil microbial communities. While this application offers great potential for deeper insights into the ecological structure of soil microbial ecosystems, it also brings new challenges related to the specific characteristics of soil datasets and the type of ecological questions that can be addressed. In this Perspectives Paper we assess the challenges of applying network analysis to soil microbial ecology due to the small-scale heterogeneity of the soil environment and the nature of soil microbial datasets. We review the different approaches of network construction that are commonly applied to soil microbial datasets and discuss their features and limitations. Using a test dataset of microbial communities from two depths of a forest soil, we demonstrate how different experimental designs and network constructing algorithms affect the structure of the resulting networks, and how this in turn may influence ecological conclusions. We will also reveal how assumptions of the construction method, methods of preparing the dataset, and definitions of thresholds affect the network structure. Finally, we discuss the particular questions in soil microbial ecology that can be approached by analyzing and interpreting specific network properties. Targeting these network properties in a meaningful way will allow applying this technique not in merely descriptive, but in hypothesis-driven research. Analysing microbial networks in soils opens a window to a better understanding of the complexity of microbial communities. However, this approach is unfortunately often used to draw conclusions which are far beyond the scientific evidence it can provide, which has damaged its reputation for soil microbial analysis. In this Perspectives Paper, we would like to sharpen the view for the real potential of microbial co-occurrence analysis in soils, and at the same time raise awareness regarding its limitations and the many ways how it can be misused or misinterpreted.

Soil microbial co-occurrence networks become less connected with soil development in a high Arctic glacier foreland succession

Classically, ecologists have considered that biota becomes more integrated and interdependent with ecosystem development in primary successional environments. However, recent work on soil microbial communities suggests that there may in fact be no change in network integration over successional time series. Here, we performed a test of this principle by identifying network-level topological features of the soil microbial co-occurrence networks in the primary successional foreland environment of the retreating high-Arctic glacier of Midtre Lovénbreen, Svalbard. Soil was sampled at sites along the foreland of inferred ages 10-90 years since deglaciation. DNA was extracted and amplicon sequenced for 16 s rRNA genes for bacteria and ITS1 region for fungi. Despite the chronologically-related soil pH decline and organic C/N accumulation, analysis on network-level topological features showed network integration did not change with inferred chronological ages, whereas network integration declined with decreasing pH and increasing total organic carbon (TOC) - both factors that can be viewed as an indicator of soil development. We also found that bacteria played a greater role in the network structure than fungi, with all keystone species in the microbial co-occurrence network being bacteria species. Both number and relative abundance of the keystone species were significantly higher when soil pH increased or TOC decreased. It appears that in the more extreme and less productive conditions of early primary succession, integration between members of soil biota into consortia may play a greater role in niche adaptation and survival. Our finding also emphasizes that ecosystem development is not simply a product of time but is influenced by locally heterogeneous factors.

Percolation in networks with local homeostatic plasticity

Percolation is a process that impairs network connectedness by deactivating links or nodes. This process features a phase transition that resembles paradigmatic critical transitions in epidemic spreading, biological networks, traffic and transportation systems. Some biological systems, such as networks of neural cells, actively respond to percolation-like damage, which enables these structures to maintain their function after degradation and aging. Here we study percolation in networks that actively respond to link damage by adopting a mechanism resembling synaptic scaling in neurons. We explain critical transitions in such active networks and show that these structures are more resilient to damage as they are able to maintain a stronger connectedness and ability to spread information. Moreover, we uncover the role of local rescaling strategies in biological networks and indicate a possibility of designing smart infrastructures with improved robustness to perturbations.

The long-term changes in food web structure and ecosystem functioning of a shallow lake: Implications for the lake management

The food web structure (FWS) and ecosystem functioning (EF) of lakes worldwide are impacted by multiple disturbances. The historical evolution of the FWS and EF are not well understood due to the lack of sufficient long-term records of biotic variations. This study reconstructed the food web models in the 1950s, 1980s, 1990s, 2000s, and 2010s for Baiyangdian Lake (BYDL), the largest shallow lake in northern China, using the Linear Inverse Modeling (LIM) and investigated EF in different periods. Our results confirmed that the FWS and EF of BYDL have undergone significant changes. The biomass of phytoplankton continuously increased, and the primary productivity of phytoplankton began to replace the primary productivity of submerged macrophytes in the 2000s and became the largest energy flow in the food web. Changes in the energy flow of primary producers are transmitted to high-trophic functional groups, which affects the diet composition of fish. Based on the ecological network analysis indices and food web stability indicators, it was concluded that after a turning point in the 1990s, the ecosystem showed initial stability and then gradually became unstable. Water level fluctuations and nutrient enrichment may be the key driving factors for changes in ecosystem state. Therefore, to maintain a good state of the ecosystem, we recommend implementing comprehensive management measures of hydrological management, nutrient-loading reduction, and biomanipulation. Furthermore, this study extended LIM to lake ecosystems, which may provide a new method for lake ecological environment management.

Soil Rehabilitation Promotes Resilient Microbiome with Enriched Keystone Taxa than Agricultural Infestation in Barren Soils on the Loess Plateau

Drylands provide crucial ecosystem and economic services across the globe. In barren drylands, keystone taxa drive microbial structure and functioning in soil environments. In the current study, the Chinese Loess plateau’s agricultural (AL) and twenty-year-old rehabilitated lands (RL) provided a unique opportunity to investigate land-use-mediated effects on barren soil keystone bacterial and fungal taxa. Therefore, soils from eighteen sites were collected for metagenomic sequencing of bacteria specific 16S rRNA and fungi specific ITS2 regions, respectively, and to conduct molecular ecological networks and construct microbial OTU-based correlation matrices. In RL soils we found a more complex bacterial network represented by a higher number of nodes and links, with a link percentage of 77%, and a lower number of nodes and links for OTU-based fungal networks compared to the AL soils. A higher number of keystone taxa was observed in the RL (66) than in the AL (49) soils, and microbial network connectivity was positively influenced by soil total nitrogen and microbial biomass carbon contents. Our results indicate that plant restoration and the reduced human interventions in RL soils could guide the development of a better-connected microbial network and ensure sufficient nutrient circulation in barren soils on the Loess plateau.

Biogeochemical extremes and compound events in the ocean

The ocean is warming, losing oxygen and being acidified, primarily as a result of anthropogenic carbon emissions. With ocean warming, acidification and deoxygenation projected to increase for decades, extreme events, such as marine heatwaves, will intensify, occur more often, persist for longer periods of time and extend over larger regions. Nevertheless, our understanding of oceanic extreme events that are associated with warming, low oxygen concentrations or high acidity, as well as their impacts on marine ecosystems, remains limited. Compound events-that is, multiple extreme events that occur simultaneously or in close sequence-are of particular concern, as their individual effects may interact synergistically. Here we assess patterns and trends in open ocean extremes based on the existing literature as well as global and regional model simulations. Furthermore, we discuss the potential impacts of individual and compound extremes on marine organisms and ecosystems. We propose a pathway to improve the understanding of extreme events and the capacity of marine life to respond to them. The conditions exhibited by present extreme events may be a harbinger of what may become normal in the future. As a consequence, pursuing this research effort may also help us to better understand the responses of marine organisms and ecosystems to future climate change.

Substrate type and plant phenolics influence epiphytic bacterial assembly during short-term succession

In natural ecosystems, large amounts of epiphytic bacteria live on the surfaces of submerged plants or non-biological substrates. Although it contributes greatly to host plant health or ecological functions in waters, little is known about the temporal dynamics and assembly mechanisms of epiphytic bacteria. To test whether host plant chemistry leads to divergent community dynamics, we investigated the fine scale temporal community successions of both epiphytic bacteria and the bacterioplankton of the surrounding water in two submerged plants and one non-biological artificial substance. We first observed differentiated epiphytic or surrounding water bacterial communities for different substrates in small spaces (approximately 1 m × 1 m). Selection played dominant roles in affecting the assembly of epiphytic bacteria in the high-phenolic plant Hydrilla verticillata, while for the artificial substance and the low-phenolic plant Vallisneria natans, drift and dispersal drove the assembly of both epiphytic bacteria and bacterioplankton. The higher selection may also contribute to higher turnover rates in both bacterioplankton and epiphytic communities of H. verticillata, with the latter changing drastically in approximately one week. Epiphytic bacteria in H. verticillata developed more complex networks with a higher proportion of positive links, suggesting that more intense interactions such as mutualism or facilitation may exist within epiphytic bacterial communities of the high-phenolic plant. Our results also implied that for the submerged macrophytes used in biological purification, the dynamics of epiphytic biofilm in the purification-related functional capacities might also be considered.

Species richness and food-web structure jointly drive community biomass and its temporal stability in fish communities

Biodiversity-ecosystem functioning and food-web complexity-stability relationships are central to ecology. However, they remain largely untested in natural contexts. Here, we estimated the links among environmental conditions, richness, food-web structure, annual biomass and its temporal stability using a standardised monitoring dataset of 99 stream fish communities spanning from 1995 to 2018. We first revealed that both richness and average trophic level are positively related to annual biomass, with effects of similar strength. Second, we found that community stability is fostered by mean trophic level, while contrary to expectation, it is decreased by species richness. Finally, we found that environmental conditions affect both biomass and its stability mainly via effects on richness and network structure. Strikingly, the effect of species richness on community stability was mediated by population stability rather than synchrony, which contrasts with results from single trophic communities. We discuss the hypothesis that it could be a characteristic of multi-trophic communities.

Phenological traits foster persistence of mutualistic networks by promoting facilitation

Morphological and phenological traits are key determinants of the structure of mutualistic networks. Both traits create forbidden links, but phenological traits can also decouple interaction in time. While such difference likely affects the indirect effects among species and consequently network persistence, it remains overlooked. Here, using a dynamic model, we show that networks structured by phenology favour facilitation over competition within guilds of pollinators and plants, thereby increasing network persistence, while the contrary holds for networks structured by morphology. We further show that such buffering of competition by phenological traits mostly beneficiate to specialists, the most vulnerable species otherwise, which propagate the most positive effects within guilds and promote nestedness. Our results indicate that beyond trophic mismatch, phenological shifts such as those induced by climate change are likely to affect indirect effects within mutualistic assemblages, with consequences for biodiversity.

Energetic constraints imposed on trophic interaction strengths enhance resilience in empirical and model food webs

Food web stability and resilience are at the heart of understanding the structure and functioning of ecosystems. Previous studies show that models of empirical food webs are substantially more stable than random ones, due to a few strong interactions embedded in a majority of weak interactions. Analyses of trophic interaction loops show that in empirical food webs the patterns of the interaction strengths prevent the occurrence of destabilizing heavy loops and thereby enhances resilience. Yet, it is still unexplored which biological mechanisms cause these patterns that enhance food web resilience. We quantified food web resilience using the real part of the maximum eigenvalue of the Jacobian matrix of the food web from a seagrass bed in the Yellow River Delta (YRD) wetland, that could be parametrized by the empirical data of the food web. We found that the empirically based Jacobian matrix of the YRD food web indicated a much higher resilience than random matrices with the same element values but arranged in random ways. Investigating the trophic interaction loops revealed that the high resilience was due to a negative correlation between the negative and positive interaction strengths (per capita top-down and bottom-up effects, respectively) within positive feedback loops with three species. The negative correlation showed that when the negative interaction strengths were strong the positive was weak, and vice versa. Our invented reformulation of loop weight in terms of biomasses and specific production rates showed that energetic properties of the trophic groups in the loop and mass-balance constraints, for example, the food uptake has to balance all losses, created the negative correlation between the interaction strengths. This result could be generalized using a dynamic intraguild predation model, which delivered the same pattern for a wide range of model parameters. Our results shed light on how energetic constraints at the trophic group and food web level create a pattern of interaction strengths within trophic interaction loops that enhances food web resilience.

A novel approach to quantifying trophic interaction strengths and impact of invasive species in food webs

Measuring ecological and economic impacts of invasive species is necessary for managing invaded food webs. Based on abundance, biomass and diet data of autochthonous and allochthonous fish species, we proposed a novel approach to quantifying trophic interaction strengths in terms of number of individuals and biomass that each species subtract to the others in the food web. This allowed to estimate the economic loss associated to the impact of an invasive species on commercial fish stocks, as well as the resilience of invaded food webs to further perturbations. As case study, we measured the impact of the invasive bass Micropterus salmoides in two lake communities differing in food web complexity and species richness, as well as the biotic resistance of autochthonous and allochthonous fish species against the invader. Resistance to the invader was higher, while its ecological and economic impact was lower, in the more complex and species-rich food web. The percid Perca fluviatilis and the whitefish Coregonus lavaretus were the two species that most limited the invader, representing meaningful targets for conservation biological control strategies. In both food webs, the limiting effect of allochthonous species against M. salmoides was higher than the effect of autochthonous ones. Simulations predicted that the eradication of the invader would increase food web resilience, while that an increase in fish diversity would preserve resilience also at high abundances of M. salmoides. Our results support the conservation of biodiverse food webs as a way to mitigate the impact of bass invasion in lake ecosystems. Notably, the proposed approach could be applied to any habitat and animal species whenever biomass and diet data can be obtained.

Changes in Soil Features and Phytomass during Vegetation Succession in Sandy Areas

This research was conducted on an area of inland sands characterised by various degrees of overgrowth by vegetation and soil stabilisation. This landscape’s origin is not natural but is connected to human industrial activities dating from early medieval times, which created a powerful centre for mining and metallurgy. This study aims to identify the changes in the above- and belowground phytomass in the initial stages of succession and their influence on the chemical properties and morphology of the soil. It was found that Salix arenaria dominated in primary phytomass production in all plots tested. The amounts of this species found in each community were as follows: 8.55 kg/400 m2 (algae–mosses), 188.97 kg/400 m2 (sand grassland–willow), 123.44 kg/400 m2 (pine–willow–mosses), 14.63 kg/400 m2 (sand grassland–mosses–willow), and 196.55 kg/400 m2 (willow–pine–sand grassland). A notable share of Koeleria glauca was found in the phytomass production of Plots IV (45.73 kg) and V (86.16 kg). Basic soil properties (pH, Corg, Nt), available plant elements (P), and plant nutrients (Ca, Mg, K, P, Fe) beneath the dominant plant species were examined. Soil acidity (pH) varied greatly, ranging from acidic (pH = 3.2) to weakly acidic (pH = 6.3). The content of organic carbon (Corg) in individual plots beneath the dominant species in the humus horizon ranged from 0.28% to 1.42%. The maximum contents of organic carbon and total nitrogen were found in organic (O) and organic-humus (OA) horizons. The highest Pavail content was found in organic and organic-humus horizons, ranging from 10.41 to 65.23 mg/kg, and in mineral horizons under K. glauca (24.10 mg/kg) and Salix acutifola (25.11 mg/kg). The soil features and phytomass were varied differently across individual sites, representing different stages of succession.

Can biomass distribution across trophic levels predict trophic cascades?

The biomass distribution across trophic levels (biomass pyramid) and cascading responses to perturbations (trophic cascades) are archetypal representatives of the interconnected set of static and dynamical properties of food chains. A vast literature has explored their respective ecological drivers, sometimes generating correlations between them. Here we instead reveal a fundamental connection: both pyramids and cascades reflect the dynamical sensitivity of the food chain to changes in species intrinsic rates. We deduce a direct relationship between cascades and pyramids, modulated by what we call trophic dissipation - a synthetic concept that encodes the contribution of top-down propagation of consumer losses in the biomass pyramid. Predictable across-ecosystem patterns emerge when systems are in similar regimes of trophic dissipation. Data from 31 aquatic mesocosm experiments demonstrate how our approach can reveal the causal mechanisms linking trophic cascades and biomass distributions, thus providing a road map to deduce reliable predictions from empirical patterns.

Geographical variation of multiplex ecological networks in marine intertidal communities

Understanding the drivers of geographical variation in species distributions, and the resulting community structure, constitutes one of the grandest challenges in ecology. Geographical patterns of species richness and composition have been relatively well studied. Less is known about how the entire set of trophic and non-trophic ecological interactions, and the complex networks that they create by gluing species together in complex communities, change across geographical extents. Here, we compiled data of species composition and three types of ecological interactions occurring between species in rocky intertidal communities across a large spatial extent (~970 km of shoreline) of central Chile, and analyzed the geographical variability in these multiplex networks (i.e., comprising several interaction types) of ecological interactions. We calculated nine network summary statistics common across interaction types, and additional network attributes specific to each of the different types of interactions. We then investigated potential environmental drivers of this multivariate network organization. These included variation in sea surface temperature and coastal upwelling, the main drivers of productivity in nearshore waters. Our results suggest that structural properties of multiplex ecological networks are affected by local species richness and modulated by factors influencing productivity and environmental predictability. Our results show that non-trophic negative interactions are more sensitive to spatially structured temporal environmental variation than feeding relationships, with non-trophic positive interactions being the least labile to it. We also show that environmental effects are partly mediated through changes in species richness and partly through direct influences on species interactions, probably associated to changes in environmental predictability and to bottom-up nutrient availability. Our findings highlight the need for a comprehensive picture of ecological interactions and their geographical variability if we are to predict potential effects of environmental changes on ecological communities.