Food Web Ecosystem Research Articles

Genetic structure of invasive ctenophore Mnemiopsis leidyi populations in temperate northern European waters is consistent with a southern North Sea overwintering refuge

The ctenophore Mnemiopsis leidyi is native to the Atlantic coasts of the Americas and is considered a threat to biodiversity and food webs in recipient ecosystems. Most of these invasive ctenophores perish during the cold winter in temperate northern European waters (North and Baltic Seas) but spread again every summer/autumn. We collected ctenophores from the southern North Sea, inner Danish waters, and western Baltic Sea in 2017 and 2018, and sequenced genomic DNA obtained from tissue samples. We obtained sequences of internal transcribed spacer (ITS) 1 and 2, and the intervening 5.8S fragment and partial fragments of the flanking 18S rRNA and 28S rRNA genes from nuclear DNA. Excluding the rRNA regions and focussing on only the two ITS sequences, we compared with ITS sequences archived in databases from specimens in its native (Northwest Atlantic) and other invaded habitats (Caspian Sea, Mediterranean Sea, central-western Atlantic Ocean). The comparison revealed no significant difference in genetic variation among temperate northern European sampling locations but showed variation when compared with native and other invaded habitats. No significant genetic difference was identified between specimens collected in the temperate northern European waters for two consecutive years. The very low interregional and interannual genetic variations in M. leidyi specimens in the Northeast Atlantic are consistent with the presence of one single panmictic population in northern European waters and with the hypothesis, also suggested in previous studies, that annual dispersal from the southeastern North Sea into the western Baltic Sea can be traced back to the same population, which most likely originates every spring from the Dutch Wadden Sea and the English Channel.

Sexual reproduction during diatom bloom.

Phytoplankton supports food webs in all aquatic ecosystems. Ecological studies highlighted the links between environmental variables and species successions in situ. However, the role of life cycle characteristics on phytoplankton community dynamics remains poorly characterized. In diatoms, sexual reproduction creates new genetic combinations and prevents excessive cell size miniaturization. It has been extensively studied in vitro but seldom in the natural environment. Here, analyzing metatranscriptomic data in the light of the expression patterns previously characterized in vitro, we identified a synchronized and transient sexual reproduction event during a bloom of the toxic diatom species Pseudo-nitzschia australis. Despite the complexity of environmental conditions encountered in situ, sexual reproduction appeared to be the strongest differential gene expression signal that occurred during the bloom. The potential link between environmental conditions and the initiation of sexual reproduction remain to be determined, but sexual reproduction probably had a major impact on the bloom dynamic.

Inhibitory Effects of Microalgal Biofilm on Vibrio parahaemolyticus

Microalgae are primary producers, the first trophic level in the food webs of aquatic ecosystems. Therefore, they have many potential applications in aquaculture. This study aimed to create autotrophic biofilms from microalgae strains belonging to the Nostocaceae, Protosiphonaceae, and Scenedesmaceae, isolated in Vietnam and to select a microalgal biofilm capable of inhibiting Vibrio parahaemolyticus. The microalgal biofilms were created by attaching the microalgal strains (cyanobacterium and green algae) to substrates (coconut fibre, sugarcane bagasse and cheesecloth). Additionally, the inhibitory effect of the microalgal biofilm on V. parahaemolyticus was evaluated based on the reduction proportion of V. parahaemolyticus cell numbers in an environment with the microalgal biofilm (co-culture method) compared to those in an environment without microalgal biofilm. The results showed that coconut fibre and cheesecloth had higher coverage density and biomass content of microalgae after 72 hours of biofilm formation than the sugarcane bagasse. We also found that microalgal biofilms formed by the combination of a cyanobacterial strain NK1111 with a green microalgal strain MBN1 or MHN122 reduced 76.2% - 99.8% of V. parahaemolyticus in 96 - 120 hours of co-culture. Meanwhile, microalgal biofilm formed by the strain NK1111 and a green microalgal strain MHB231 inhibited 91.8% of V. parahaemolyticus bacterium as early as 24 hours and inhibited 99.4% of this bacteria at 120 hours of co-culture. Therefore, the microalgal biofilms consisting of the strain NK1111 and the strain MHB231 on cheesecloth or coconut fibre substrate are highly appreciated for their potential application in shrimp farms in Vietnam.

Reindeer carcasses modulate vegetation composition and greenness in High-Arctic tundra

Vertebrate carrion is an integral part of foodwebs in ecosystems and can impact biodiversity at the local as well as the landscape scale. However, very little knowledge currently exists about the ecological role of carrion in the Arctic ecosystems. We conducted a ground survey on the cover of five plant functional groups at paired reindeer carcass and control sites and analysed the relationship between cover and carcass presence in the Arctic tundra of Svalbard. Vegetation indices from Red–Green–Blue (RGB) imagery captured by drones complemented this, assessing plant productivity in terms of “spectral greening” and modelling the relationship between vegetation index values and carcass distance. We show that graminoids capitalised most from carcass presence, whereas bryophytes and lichen showed decreases in cover. Woody plant and forb covers were not significantly impacted by carcass presence. The Red Green Blue Vegetation Index decreased locally at fresh carcasses (i.e., <1 year old) but showed an increase at more established carcass sites (i.e., >1 year). We show that carcasses have differential impacts on the plant functional groups of Svalbard's tundra and induce a local “green-up” through secondary succession within 2 m of the carcass. Given their non-random distribution, carcasses may contribute to vegetation heterogeneity at landscape scales. This is relevant for understanding how climate change-induced reindeer mortalities will impact tundra plant community composition in the future.

Enhancing yield and economic benefits through sustainable pest management in Okra cultivation.

Okra (Abelmoschus esculentus) is a prominent vegetable crop in Asia, confronting persistent threats from pests such as leafhoppers, whiteflies, and shoot and fruit borers. Conventional chemical control methods, despite their adverse ecological effects, remain the primary approach for pest management. Indiscriminate chemical use has led to reduced biodiversity among natural predators and the disruption of food webs in ecosystems. To address these challenges, this study assessed the efficacy of integrated (IM) and biointensive (BM) pest management modules in comparison to conventional chemical methods (CM) for mitigating insect damage to okra leaves and fruits, and subsequently, their impact on okra yield. Our result revealed that the BM exhibited the least effectiveness but outperformed untreated control plots significantly. In contrast, both IM and CM significantly reduced damage from sap-sucking insects and borer pests. Notably, plots treated with the chemical module found decreased populations of natural enemies. The IM demonstrated the lowest fruit infestation rate (5.06%), yielding the highest crop production (8.97 t ha-1), along with the maximum net return (Indian Rupees: 44,245) and incremental cost-benefit ratio (3.31). Thus, the study suggested that the implementation of integrated pest management practices can result in higher okra yields and greater economic benefits. These findings shed light on the potential of sustainable agricultural practices as a safer and more economically viable alternative to chemical-intensive pest control in okra cultivation.

Effects of Aquatic Physicochemical Parameters Variation on the Phytoplankton Abundance and Diversity in the Babon River, (Semarang, Central Java, Indonesia)

The Babon River is an integral component of Semarang City's drainage system, which is heavily impacted by the influx of waste, leading to elevated water temperatures caused by global warming and decreased pH levels due to acid rain. The two main environmental factors that affect phytoplankton growth and productivity are temperature and nutrient availability. Nutrients can change the balance of the natural food webs of aquatic ecosystems. Meanwhile, metals are toxic and can accumulate in biota tissues. This study examined how nutrients, pH, temperature, and heavy metals in water affected the diversity and abundance of phytoplankton. Water samples and phytoplankton were collected from 7 observation stations representing the upstream, middle, and downstream rivers. Observations at each station were made 3 times in the months of April-May, June-July, and August-September. Temperature and pH were measured in situ with the HORIBA Water Checker. Analysis of nitrate, phosphate, and heavy metals was carried out with an AAS. The results showed that the dominant phytoplankton in the Babon Rivers consisted of Bacillariophyceae, Chlorophyceae, and Cyanophyceae. There is no diversity against all independent variables (temperature, pH, nitrate, phosphate, Cd, Total Cr, Pb concentrations) and dependent variables (number of genera, abundance, and diversity index of phytoplankton) due to the sampling station's location. However, sampling time caused diversity in nitrate, phosphate, Cd, Total Cr, and Pb concentrations, the sampling location effect was more dominant than sampling time, so the results showed that all independent variables did not contribute significantly to the dependent variables.

Hydrodynamic and trophic variations reshape macroinvertebrate food webs in urban ecosystems

Urbanization is modifying aquatic ecosystems, with hydrodynamic and trophic variations altering biotic assemblages in rapidly expanding cities worldwide. Despite the fundamental bioenergetic role of food webs within these assemblages, their responding mechanism to the hydrodynamic and trophic variations remains largely unknown. Here we show that hydrodynamic and trophic loss, coupled with the weakening of cascade controls by key trophic guilds, leads to a significant decline in the structure, function and stability of macroinvertebrate food webs. Utilizing the allometric diet breadth model and biomass balance model, we established representative food webs for macroinvertebrate groups under varying hydrodynamic and trophic stresses. We found that such losses have reduced ∼75% trophic guild richness, ∼85% biomass flux, and ∼80% biomass storage. These reductions promote trophic guild specialization, further destabilizing food web, eroding interactive strength asymmetry, and diminishing the control of trophic guilds. Furthermore, macroinvertebrate food webs show divergent stability responses under similar stress levels, mainly driven by differences in the cascade controls exerted by key trophic guilds. Our results underscore the critical role of hydrodynamic and trophic variations in shaping urban aquatic ecosystems and highlight the significance of both external environmental revitalization and internal food web dynamics enhancement in restoring the ecological stability in urban settings.

Temporal variations of food web in a marine bay ecosystem based on LIM-MCMC model

Temporal variations of food web in a marine bay ecosystem based on LIM-MCMC model

Determining key mammalian species and food web robustness across different land cover vegetation using network analysis

Abstract. Risdiyanto I, Santosa Y, Santoso N, Sunkar A. 2024. Determining key mammalian species and food web robustness across different land cover vegetation using network analysis. Biodiversitas 25: 3175-3190. The ecosystem's food web depicts the intricate energy flow and complex interactions among its diverse organisms. Organisms are not confined to single trophic levels or food chains; they establish multiple food connections with other organisms within the ecosystem. This study aims to identify critical mammalian species and assess the robustness of various ecosystem types, considering the diversity of species presence, through food web analysis. Computational graphic-based network analysis is employed to achieve this goal. Ecosystem types, categorized by their land cover, include plantation/agricultural forests, Bushes and shrubss, forests, and mixed landscapes. Network centrality metrics such as degree, closeness, betweenness, and eigenvector centrality are utilized to evaluate the presence of key species and ecosystem robustness. The relative contribution of mammalian species as connectors and regulators of energy flow in the food web ranges from 8-23% of all nodes involved. Key mammalian species are classified into ecosystem stability and sustainability keys. In Bushes and shrubss ecosystems, key species predominantly consist of mammalian predator species that are crucial for maintaining ecosystem stability through population control. Conversely, in other ecosystems, key species are primarily connectors, ensuring the sustained energy flow. The most resilient food webs are observed in Bushes and shrubss and mixed ecosystems due to their higher biomass growth rates and abundant presence of mammalian species. Utilizing food web analysis can significantly contribute to species and ecosystem conservation efforts by offering a comprehensive understanding of interspecies interactions, food web structures, and ecosystem dynamics.

Turbidity effects on prey consumption and survival of larval European smelt (Osmerus eperlanus)

The anadromous European smelt (Osmerus eperlanus) plays a key role in food webs of many riverine ecosystems in Europe. However, population sizes in several German rivers (e.g. Elbe or Weser rivers) have diminished drastically over the past decade. Turbidity has been considered one of the stressors affecting the successful recruitment of European smelt, as their early life stages may be particularly sensitive to changes in the abiotic environment. In this study, we investigated whether prey consumption and survival of European smelt larvae would be negatively affected by an acute exposure to elevated turbidity. We reared the larvae in the laboratory and exposed them in four separate trials (18 to 26 days post hatch, 9.5 ± 0.8 mm standard length, mean ± SD) to six turbidity levels (0–500 NTU, nephelometric turbidity units). We found that prey uptake increased at low turbidity levels and decreased at high turbidity levels, with an optimum between 100 and 200 NTU. Survival started to decrease at turbidity levels above 300 NTU. In addition, we conducted a systematic literature analysis in which we found that prey consumption of larval and juvenile fishes had been tested across a wide range of turbidity levels, mostly using pelagic (e.g. planktonic) prey items, with more studies focusing on perciform fishes and juvenile rather than larval life stages. Our empirical findings contribute to establishing thresholds for optimal larval European smelt performance under increased turbidity and provide valuable information for developing mechanistic models that assess potential consequences for European smelt recruitment dynamics.

Effects of multiple stressors on freshwater food webs: Evidence from a mesocosm experiment

Natural and anthropogenic pressures exert influence on ecosystem structure and function by affecting the physiology and behavior of organisms, as well as the trophic interactions within assemblages. Therefore, understanding how multiple stressors affect aquatic ecosystems can improve our ability to manage and protect these ecosystems and contribute to understanding fundamental ecological principles. Here, we conducted a mesocosm experiment to ascertain the individual and combined effects of multiple stressors on trophic interactions within species in freshwater ecosystems. Furthermore, we investigated how species respond to such changes by adapting their food resources. To mimic a realistic food web, we selected fish and shrimp as top predators, gastropods, zooplankton and zoobenthos as intermediate consumers, with producers (macrophytes, periphyton and phytoplankton) and detritus as basal resources. Twelve different treatments included a control, nutrient loading only, herbicide exposure only, and a combination of nutrient loading and herbicide exposure, each replicated under ambient temperature, constant warming and multiple heat waves to simulate environmental stressors. Our results demonstrated that antagonistic interactions between environmental stressors were widespread in trophic interactions, with a more pronounced and less intense impact observed for the high trophic level species. The responses of freshwater communities to environmental stressors are complex, involving direct effects on individual species as well as indirect effects through species interactions. Moreover, our results confirmed that the combinations of stressors, but not individual stressors, led to a shift to herbivory in top predators, indicating that multiple stressors can be more detrimental to organisms than individual stressors alone. These findings elucidate how changes in the resource utilization of species induced by environmental stressors can potentially influence species interactions and the structural dynamics of food webs in freshwater ecosystems.

Overview of Predation by Birds, Cephalopods, Fish and Marine Mammals on Marine Benthic Amphipods

With about 8000 marine benthic species, the amphipod crustaceans form one of the richest animal groups of the worldwide Ocean. They have colonized a wide range of soft- and hard-bottom natural and artificial habitats extending from the intertidal to hadal zones. Moreover, they show a broad size spectrum, with numerous giant species exceeding 20 cm in length and some species smaller than 2 mm. When biofouling artificial hard surfaces, some tube-building species can form very dense populations comprising up to 100,000 individuals per square meter. Amphipods are important prey for fish and mammals. Along with cephalopod juveniles, they are also included in the trophic diet of shorebirds that consume amphipods mostly during the low tide on tidal flats. They display diel migration, which reinforces the predation by demersal fish in the suprabenthic zone just above the sea bed, as well as by pelagic fish in the water column. Despite their importance in terms of biodiversity and trophic transfer, no general overview is available on the role of benthic amphipods in marine ecosystem food webs. Various methods, including laboratory and field experiments, as well as the analysis of stomach contents and DNA extraction, have been used to identify the prey/predator trophic links. Based on an extensive literature review, this study discusses the role of marine benthic amphipods as potential food for higher trophic levels in natural and artificial hard-bottom communities created via the construction of offshore wind farms.

Size-dependence of food intake and mortality interact with temperature and seasonality to drive diversity in fish life histories.

Understanding how growth and reproduction will adapt to changing environmental conditions is a fundamental question in evolutionary ecology, but predicting the responses of specific taxa is challenging. Analyses of the physiological effects of climate change upon life history evolution rarely consider alternative hypothesized mechanisms, such as size-dependent foraging and the risk of predation, simultaneously shaping optimal growth patterns. To test for interactions between these mechanisms, we embedded a state-dependent energetic model in an ecosystem size-spectrum to ask whether prey availability (foraging) and risk of predation experienced by individual fish can explain observed diversity in life histories of fishes. We found that asymptotic growth emerged from size-based foraging and reproductive and mortality patterns in the context of ecosystem food web interactions. While more productive ecosystems led to larger body sizes, the effects of temperature on metabolic costs had only small effects on size. To validate our model, we ran it for abiotic scenarios corresponding to the ecological lifestyles of three tuna species, considering environments that included seasonal variation in temperature. We successfully predicted realistic patterns of growth, reproduction, and mortality of all three tuna species. We found that individuals grew larger when environmental conditions varied seasonally, and spawning was restricted to part of the year (corresponding to their migration from temperate to tropical waters). Growing larger was advantageous because foraging and spawning opportunities were seasonally constrained. This mechanism could explain the evolution of gigantism in temperate tunas. Our approach addresses variation in food availability and individual risk as well as metabolic processes and offers a promising approach to understand fish life-history responses to changing ocean conditions.

A non-lethal stable isotope analysis of valued freshwater predatory fish using blood and fin tissues as alternatives to muscle tissue.

Stable isotope analysis (SIA) is widely used to study trophic ecology and food webs in aquatic ecosystems. In the case of fish, muscle tissue is generally preferred for SIA, and the method is lethal in most cases. We tested whether blood and fin clips can be used as non-lethal alternatives to muscle tissue for examining the isotopic composition of two freshwater predatory fish, European catfish (Silurus glanis) and Northern pike (Esox lucius), species of high value for many freshwater systems as well as invasive species in many others. Blood samples from the caudal vein, anal fin clips, and dorsal muscle obtained by biopsy punch were collected from four catfish and pike populations (14-18 individuals per population). Subsequently, these samples were analyzed for δ13C and δ15N. The effects of alternative tissues, study site, and fish body mass on the isotopic offset were investigated. Both species showed a correlation between the isotopic offset and the tissue type, as well as the study site, but no significant relationship with the body mass. The isotopic offsets between tissues were used to calculate the conversion equations. The results demonstrated that both blood and fin clips are suitable and less invasive alternative to muscle in SIA studies focused on European catfish and Northern pike. Blood provided better correspondence to muscle isotope values. However, our results clearly demonstrated that isotopic offsets between tissues vary significantly among populations of the same species. Therefore, obtaining a muscle biopsy from several individuals in any population is advisable to gain initial insights and establish a possible population-specific inter-tissue conversion.

Country-Wide Ecological Health Assessment Methodology for Air Toxics: Bridging Gaps in Ecosystem Impact Understanding and Policy Foundations.

Amid the growing concerns about air toxics from pollution sources, much emphasis has been placed on their impacts on human health. However, there has been limited research conducted to assess the cumulative country-wide impact of air toxics on both terrestrial and aquatic ecosystems, as well as the complex interactions within food webs. Traditional approaches, including those of the United States Environmental Protection Agency (US EPA), lack versatility in addressing diverse emission sources and their distinct ecological repercussions. This study addresses these gaps by introducing the Ecological Health Assessment Methodology (EHAM), a novel approach that transcends traditional methods by enabling both comprehensive country-wide and detailed regional ecological risk assessments across terrestrial and aquatic ecosystems. EHAM also advances the field by developing new food-chain multipliers (magnification factors) for localized ecosystem food web models. Employing traditional ecological multimedia risk assessment of toxics' fate and transport techniques as its foundation, this study extends US EPA methodologies to a broader range of emission sources. The quantification of risk estimation employs the quotient method, which yields an ecological screening quotient (ESQ). Utilizing Kuwait as a case study for the application of this methodology, this study's findings for data from 2017 indicate a substantial ecological risk in Kuwait's coastal zone, with cumulative ESQ values reaching as high as 3.12 × 103 for carnivorous shorebirds, contrasted by negligible risks in the inland and production zones, where ESQ values for all groups are consistently below 1.0. By analyzing the toxicity reference value (TRV) against the expected daily exposure of receptors to air toxics, the proposed methodology provides valuable insights into the potential ecological risks and their subsequent impacts on ecological populations. The present contribution aims to deepen the understanding of the ecological health implications of air toxics and lay the foundation for informed, ecology-driven policymaking, underscoring the need for measures to mitigate these impacts.

Study on Sex Ratio of Lampreys Based on Simulated Ecosystem-Food Web Model

Study on Sex Ratio of Lampreys Based on Simulated Ecosystem-Food Web Model

Increased contribution of parasites in microbial eukaryotic communities of different Aegean Sea coastal systems.

Protistan communities have a major contribution to biochemical processes and food webs in coastal ecosystems. However, related studies are scarce and usually limited in specific groups and/or sites. The present study examined the spatial structure of the entire protistan community in seven different gulfs and three different depths in a regional Mediterranean Sea, aiming to define taxa that are important for differences detected in the marine microbial network across the different gulfs studied as well as their trophic interactions. Protistan community structure analysis was based on the diversity of the V2-V3 hypervariable region of the 18S rRNA gene. Operational taxonomic units (OTUs) were identified using a 97% sequence identity threshold and were characterized based on their taxonomy, trophic role, abundance and niche specialization level. The differentially abundant, between gulfs, OTUs were considered for all depths and interactions amongst them were calculated, with statistic and network analysis. It was shown that Dinophyceae, Bacillariophyta and Syndiniales were the most abundant groups, prevalent in all sites and depths. Gulfs separation was more striking at surface corroborating with changes in environmental factors, while it was less pronounced in higher depths. The study of differentially abundant, between gulfs, OTUs revealed that the strongest biotic interactions in all depths occurred between parasite species (mainly Syndiniales) and other trophic groups. Most of these species were generalists but not abundant highlighting the importance of rare species in protistan community assemblage. Overall this study revealed the emergence of parasites as important contributors in protistan network regulation regardless of depth.

Conversion of methane to organic acids is a widely found trait among gammaproteobacterial methanotrophs of freshwater lake and pond ecosystems.

Aerobic gammaproteobacterial methanotrophs (gMOB) are key organisms controlling methane fluxes at the oxic-anoxic interfaces of freshwater ecosystems. Under hypoxic environments, gMOB may shift their aerobic metabolism to fermentation, resulting in the production of extracellular organic acids. We recently isolated a gMOB strain representing the Methylobacter spp. of boreal lake water columns (i.e., Methylobacter sp. S3L5C) and demonstrated that it converts methane to organic acids (acetate, formate, malate, and propionate) under hypoxic conditions. Annotation for putative genes encoding organic acid production within the isolate's genome and in environmental metagenome-assembled genomes (MAGs) representing Methylobacter spp. suggests that the potential for methane conversion into organic acids is widely found among Methylobacter spp. of freshwater ecosystems. However, it is not known yet whether the capability to convert methane to organic acids is restricted to Methylobacter spp. or ubiquitously present among other freshwater gMOB genera. Therefore, we isolated representatives of two additional gMOB genera from the boreal lake water columns, i.e., Methylomonas paludis S2AM and Methylovulum psychrotolerans S1L, and demonstrated similar bioconversion capacities. These genera could convert methane to organic acids, including acetate, formate, succinate, and malate. Additionally, S2AM produced lactate. Furthermore, we detected genes encoding organic acid production within their genomes and in MAGs representing Methylomonas spp. and Methylovulum spp. of lake and pond ecosystems. Altogether, our results demonstrate that methane conversion to various organic acids is a widely found trait among lake and pond gMOB, highlighting their role as pivotal mediators of methane carbon into microbial food webs of freshwater lake and pond ecosystems. IMPORTANCE Aerobic gammaproteobacterial methanotrophic bacteria (gMOB) play an important role in reducing methane emissions from freshwater ecosystems. In hypoxic conditions prevalent near oxic-anoxic interfaces, gMOB potentially shift their metabolism to fermentation, resulting in the conversion of methane to extracellular organic acids, which would serve as substrates for non-methanotrophic microbes. We intended to assess the prevalence of fermentation traits among freshwater gMOB. Therefore, we isolated two strains representing relevant freshwater gMOB genera, i.e., Methylovulum and Methylomonas, from boreal lakes, experimentally showed that they convert methane to organic acids and demonstrated via metagenomics that the fermentation potential is widely dispersed among lake and pond representatives of these genera. Combined with our recent study showing coherent results from another relevant freshwater gMOB genus, i.e., Methylobacter, we conclude that the conversion of methane to organic acids is a widely found trait among freshwater gMOB, highlighting their role as pivotal mediators of methane carbon into microbial food webs.

Predation has a significant impact on the complexity and stability of microbial food webs in subalpine lakes.

As an important part of microbial food webs, protists transfer organic carbon and nutrients to higher trophic levels in aquatic ecosystems. Protist predation often influences the abundance and composition of bacterial communities. However, we still do not understand whether and how predation affects the complexity and stability of microbial food webs. This study assessed the seasonal dynamic characteristics and driving factors of microbial food webs in terms of complexity and stability. Our findings have implications for future surveys to reveal the effects of climate and environmental changes.

Energy flow across soil food webs of different ecosystems: Food webs with complex structures support higher energy flux

Energy flow across soil food webs of different ecosystems: Food webs with complex structures support higher energy flux