Invasive species impacts on ecosystem structure and function: A comparison of the Bay of Quinte, Canada, and Oneida Lake, USA, before and after zebra mussel invasion

As invasion rates of exotic species increase, an ecosystem level understanding of their impacts is imperative for predicting future spread and consequences. We have previously shown that network analyses are powerful tools for understanding the effects of exotic species perturbation on ecosystems. We now use the network analysis approach to compare how the same perturbation affects another ecosystem of similar trophic status. We compared food web characteristics of the Bay of Quinte, Lake Ontario (Canada), to previous research on Oneida Lake, New York (USA) before and after zebra mussel (Dreissena polymorpha) invasion. We used ecological network analysis (ENA) to rigorously quantify ecosystem function through an analysis of direct and indirect food web transfers. We used a social network analysis method, cohesion analysis (CA), to assess ecosystem structure by organizing food web members into subgroups of strongly interacting predators and prey. Together, ENA and CA allowed us to understand how food web structure and function respond simultaneously to perturbation. In general, zebra mussel effects on the Bay of Quinte, when compared to Oneida Lake, were similar in direction, but greater in magnitude. Both systems underwent functional changes involving focused flow through a small number of taxa and increased use of benthic sources of production; additionally, both systems structurally changed with subgroup membership changing considerably (33% in Oneida Lake) or being disrupted entirely (in the Bay of Quinte). However, the response of total ecosystem activity (as measured by carbon flow) differed between both systems, with increasing activity in the Bay of Quinte, and decreasing activity in Oneida Lake. Thus, these analyses revealed parallel effects of zebra mussel invasion in ecosystems of similar trophic status, yet they also suggested that important differences may exist. As exotic species continue to disrupt the structure and function of our native ecosystems, food web network analyses will be useful for understanding their far-reaching effects.     

Evidence for the dominance of indirect effects in 50 trophic ecosystem networks

Indirect effects are powerful influences in ecosystems that may maintain species diversity and alter apparent relationships between species in surprising ways. Here, we applied network environ analysis to 50 empirically-based trophic ecosystem models to test the hypothesis that indirect flows dominate direct flows in ecosystem networks. Further, we used Monte Carlo based perturbations to investigate the robustness of these results to potential error in the underlying data. To explain our findings, we further investigated the importance of the microbial food web in recycling energy-matter using components of the Finn Cycling Index and analysis of environ centrality. We found that indirect flows dominate direct flows in 37/50 (74.0%) models. This increases to 31/35 (88.5%) models when we consider only models that have cycling structure and a representation of the microbial food web. The uncertainty analysis reveals that there is less error in the I/D values than the ±5% error introduced into the models, suggesting the results are robust to uncertainty. Our results show that the microbial food web mediates a substantial percentage of cycling in some systems (median = 30.2%), but its role is highly variable in these models, in agreement with the literature. Our results, combined with previous work, strongly suggest that indirect effects are dominant components of activity in ecosystems.     

Invasive species impacts on ecosystem structure and function: A comparison of Oneida Lake, New York, USA, before and after zebra mussel invasion

Exotic species invasion is widely considered to affect ecosystem structure and function. Yet, few contemporary approaches can assess the effects of exotic species invasion at such an inclusive level. Our research presents one of the first attempts to examine the effects of an exotic species at the ecosystem level in a quantifiable manner. We used ecological network analysis (ENA) and a social network analysis (SNA) method called cohesion analysis to examine the effect of zebra mussel (Dreissena polymorpha) invasion on the Oneida Lake, New York, USA, food web. We used ENA to quantify ecosystem function through an analysis of food web carbon transfer that explicitly incorporated flow over all food web paths (direct and indirect). The cohesion analysis assessed ecosystem structure through an organization of food web members into subgroups of strongly interacting predators and prey. Our analysis detected effects of zebra mussel invasion throughout the entire Oneida Lake food web, including changes in trophic flow efficiency (i.e., carbon flow among trophic levels) and alterations of food web organization (i.e., paths of carbon flow) and ecosystem activity (i.e., total carbon flow). ENA indicated that zebra mussels altered food web function by shunting carbon from pelagic to benthic pathways, increasing dissipative flow loss, and decreasing ecosystem activity. SNA revealed the strength of zebra mussel perturbation as evidenced by a reorganization of food web subgroup structure, with a decrease in importance of pelagic pathways, a concomitant rise of benthic pathways, and a reorganization of interactions between top predator fish. Together, these analyses allowed for a holistic understanding of the effects of zebra mussel invasion on the Oneida Lake food web.     

Comparing network analysis methodologies for consumer–resource relations at species and ecosystems scales

This research compares two existing methodologies, mixed trophic impact analysis and utility analysis, which use network analysis to evaluate the direct, pair-wise, and indirect, holistic, ecological relations between ecosystem compartments. The two approaches have many similarities, but differ in some key assumptions which affect both the final results and interpretations. Here, we briefly introduce both methodologies through a series of two simple examples; a 3-compartment competition model and a 3-compartment food chain model, and then apply the methodologies to a 15-compartment ecosystem model of the Chesapeake Bay. This example demonstrates how implementing the various conceptual and methodological assumptions lead to differing results. Notably, the overall number of positive relations is greatly affected by the treatment of the self-interactions and the handling of detritus compartments lead to a distinction between ecological or trophic relations. We recommend slight changes to both methodologies, not necessarily in order to bring them completely together, but because each has some points which are stronger and better defensible.     

The concepts of emergent and collective properties in individual-based models—Summary and outlook of the Bornhöved case studies

Ecology requires the conceptual and technical ability to analyse complex and dynamic systems consisting of a high and variable number of components and relations. These components are part of a variable interaction structure in a spatially heterogeneous context. The components of ecological interaction networks can give rise to self-organised, and scale-dependent interaction patterns and processes, which are the underlying causes of the overall ecological systems states.The individual-based modelling approach provides a widely applicable simulation framework based on a ‘hierarchy theory’ view of ecological systems.Here, we summarise and generalise the theoretical implications of the modelling studies presented in this volume in the field of terrestrial and aquatic, animal and plant ecology. The case studies cover a representative profile of processes related to ecological applications, such as food web interactions, population dynamics, dispersal, energy physiology, nutrient allocation and mutual impact of morphological and physiological development. The generic approach applied in this context allows a hierarchical representation of ecological systems and their components. Model results are obtained as self-organised structural relation networks and as aggregated quantitative states. In order to address different model characteristics we distinguish collective and emergent properties. Collective properties are those that are attributed equally to different organisation levels of the system. Emergent properties result from the activities of lower level entities on a higher organisation level, while not being present on the lower level. They can be subdivided into aggregational and connective properties. Emergent properties that are aggregational are those which emerge as a result of an aggregation procedure by an observer on the higher level which does not make sense or is not applicable on lower levels. Emergent properties that are connective, however, are based on an interaction network of lower level entities, which brings about the specific system characteristic.This classification of model results will allow to generalise the achievements and potential of the individual-based modelling approach in ecology.     

The effects of aggregation on the performance of the inverse method and indicators of network analysis

Food webs are usually aggregated into a manageable size for their interpretation and analysis. The aggregation of food web components in trophic or other guilds is often at the choice of the modeler as there is little guidance in the literature as to what biases might be introduced by aggregation decisions. We examined the impacts of the choice of the a priori model on the subsequent estimation of missing flows using the inverse method and on the indices derived from ecological network analysis of both inverse method-derived flows and on the actual values of flows, using the fully determined Sylt-Rømø Bight food web model. We used the inverse method, with the least squares minimization goal function, to estimate ‘missing’ values in the food web flows on 14 aggregation schemes varying in number of compartments and in methods of aggregation. The resultant flows were compared to known values; the performance of the inverse method improved with increasing number of compartments and with aggregation based on both habitat and feeding habits rather than diet similarity. Comparison of network analysis indices of inverse method-derived flows with that of actual flows and the original value for the unaggregated food web showed that the use of both the inverse method and the aggregation scheme affected indices derived from ecological network analysis. The inverse method tended to underestimate the size and complexity of food webs, while an aggregation scheme explained as much variability in some network indices as the difference between inverse-derived and actual flows. However, topological network indices tended to be most robust to both the method of determining flows and to the inverse method. These results suggest that a goal function other than minimization of flows should be used when applying the inverse method to food web models. Comparison of food web models should be done with extreme care when different methodologies are used to estimate unknown flows and to aggregate system components. However, we propose that indices such as relative ascendency and relative redundancy are most valuable for comparing ecosystem models constructed using different methodologies for determining missing flows or for aggregating system components.     

Ecological network analysis: network construction

Ecological network analysis (ENA) is a systems-oriented methodology to analyze within system interactions used to identify holistic properties that are otherwise not evident from the direct observations. Like any analysis technique, the accuracy of the results is as good as the data available, but the additional challenge is that the data need to characterize an entire ecosystem's flows and storages. Thus, data requirements are substantial. As a result, there have, in fact, not been a significant number of network models constructed and development of the network analysis methodology has progressed largely within the purview of a few established models. In this paper, we outline the steps for one approach to construct network models. Lastly, we also provide a brief overview of the algorithmic methods used to construct food web typologies when empirical data are not available. It is our aim that such an effort aids other researchers to consider the construction of such models as well as encourages further refinement of this procedure.     

Impact of Arenicola marina (Polychaeta) on the microbial assemblages and meiobenthos in a marine intertidal flat

The benthic microbial food web can be responsible for a large proportion of benthic carbon cycling yet there are few data on the trophic interactions between this food web and macrobenthos. A large-scale field experiment was conducted to investigate effects of eliminating the polychaete Arenicola marina on benthic microbes (prokaryotes, heterotrophic and autotrophic protists) and metazoan meiofauna in a marine intertidal flat of the North Sea, Germany. Over a period of 2 years, quantity and composition of micro- and meiobenthos from unmanipulated sites were compared to those from sites deplete of lugworms. These grazer treatments were cross-classified with different sediment characteristics (low- and mid-intertidal areas). Lugworm removal resulted in an initial increase in abundance of prokaryotes and nanoflagellates, which became less pronounced in the second year. Ciliates were not affected quantitatively, but in the absence of lugworms, diversity and the proportion of carnivorous forms increased. Meiobenthos (nematodes, ostracods and copepods) were affected only moderately. The observed changes are probably due to a combination of release from grazing/predation pressure, changes in the species composition of higher trophic levels (namely large polychaetes) and altered environmental conditions (such as depth of the oxygenated layer and sediment grain size). Spatial differences between sites of different tidal exposure/grain size appeared to be as large as temporal differences during the 2 years following the manipulation of the system. We conclude that in intertidal sediments, indirect effects due to habitat transformation are as important as direct biological interactions (grazing pressure and competition) for the dynamics of the benthic microbial food web.     

The modified niche model: Including detritus in simple structural food web models

Food webs are constructed as structural directed graphs that describe “who eats whom,” but it is common to interpret them as energy flow diagrams where predation represents an energy transfer from the prey to the predator. It is the aim of this work to demonstrate that food webs are incomplete as energy flow diagrams if they ignore passive flows to detritus (dead organic material). While many ecologists do include detritus in conceptual and mathematical models, the detrital omission is still commonly found. Often detritus is either ignored or treated as an unlimited energy source, yet all organisms contribute to the detritus pool, which can be an energy source for other species in the system. This feedback loop is of high importance, since it increases the number of pathways available for energy flows, revealing the significance of indirect effects, and making the functional role of the top predators less clear. In this work we propose the modified niche model by adding a detritus compartment to the niche model. We demonstrate the effect of structural loops that result from feeding on detritus, by comparing empirical data sets to five different assembly models: (1) cascade, (2) constant connectance, (3) niche, (4) modified niche (original in this work), and (5) cyber-ecosystem. Of these models, only the last two explicitly include detritus. We show that when passive flows to detritus are included in the food web structure, the structure becomes more robust to the removal of individual nodes or connections. In addition, we show that food web models that include the detritus feedback loop perform better with respect to several structural network metrics.     

Comparative network analysis toward characterization of systemic organization for human–environmental sustainability

A preliminary study in comparative ecological network analysis was conducted to identify key assumptions and methodological challenges, test initial hypotheses and explore systemic and network structural characteristics for environmentally sustainable ecosystems. A nitrogen network for the U.S. beef supply chain – a small sub-network of the industrial food system analyzed as a pilot study – was constructed and compared to four non-human carbon and nitrogen trophic networks for the Chesapeake Bay and the Florida Everglades. These non-human food webs served as sustainable reference systems. Contrary to the main original hypothesis, the “window of vitality” and the number of network roles did not clearly differentiate between a human sub-network and the more complete non-human networks. The effective trophic level of humans (a partial estimate of trophic level based on the single food source of beef) was much higher (8.1) than any non-human species (maximum of 4.88). Network connectance, entropy, total dependency coefficients, trophic efficiencies and the ascendency to capacity ratio also indicated differences that serve as hypotheses for future tests on more comprehensive human food webs. The study elucidated important issues related to (1) the steady state assumption, which is more problematic for industrial human systems, (2) the absence or dearth of data on contributions of dead humans and human wastes to feed other species in an integrated food web, (3) the ambiguity of defining some industrial compartments as living versus non-living, and (4) challenges with constructing compartments and trophic transfers in industrial versus non-human food webs. The two main novel results are (1) the progress made toward adapting ecological network analysis (ENA) methodology for analysis of human food networks in industrial cultures and (2) characterizing the critical aspects of comparative ENA for understanding potential causes of the problems, and providing avenues for solutions, for environmental sustainability. Based on this work, construction and comparative network analysis of a more comprehensive industrial human food network seems warranted and likely to provide valuable insights for modifying structures of industrial food networks to be more like natural networks and more sustainable.     

Parent-absent signalling of need and its consequences for sibling competition in the barn swallow

In species with parental care, competition among siblings for access to limited parental resources is common. Sibling competition can be mediated by begging behaviour, a suite of different visual and acoustic displays by which offspring solicit parental care. These are mostly addressed to the parents upon food provisioning, but can also be performed in the absence of the attending parents. This so-called parent-absent begging (PAB) may function as an intrabrood communication signal and potentially affect intrabrood competition dynamics for access to food. We investigated the role of PAB in moulding sibling interactions and its effect on food intake among altricial barn swallow (Hirundo rustica) nestlings, both under normal and experimentally reduced food intake. Frequency of PAB increased after food deprivation. Nestlings that had performed PAB increased their begging intensity upon the subsequent parental feeding visit, while siblings reduced their own begging level, but only when they had not been food-deprived. As a consequence, nestlings which had performed PAB before parental arrival had larger chances of receiving food. However, nestlings did not benefit from displaying PAB when competing with food-deprived siblings. Our findings show that PAB reliably reflects need of food, indicating that a nestling will vigorously compete for the subsequent food item. By eavesdropping siblings' PAB displays, nestlings may optimally balance the costs of scrambling competition, the direct fitness gains of being fed and the indirect fitness costs of subtracting food to needy kin. However, large asymmetries in satiation between competitors may lead individual offspring to monopolize parental resources, irrespective of PAB displays.     

Nutzen und Grenzen von Modellökosystemen zur Erfassung der Wirkung von Schadstoffen auf benthische Lebensgemeinschaften

The aim of this short article is to deliver insight into the prospects and constraints of model ecosystems in the assessment of effects of contaminants on benthic communities, especially the meiofauna. Numerous substances accumulate in sediments and thus a complex contamination exists which is able to impair benthic organisms. Besides direct effects, indirect effects are of importance and should be considered. Changes in competition and predator-prey-relationships lead to a variety of effects in the benthic food web that cannot be explained by direct toxicity of contaminants. In this article two microcosm studies, one with the heavy metal cadmium and one with the veterinary pharmaceutical Ivermectin, are presented which both showed direct as well as indirect effects. A significant change in abundance and composition of the meiofauna community has been observed. The prediction of field-effects of a certain substance as well as the interpretation of field data are the aims of model ecosystems. Both studies proofed that the investigation of meiofauna in microcosms is a promising tool to achieve these aims. Furthermore, microcosms with meiofauna might be able to assist the effect-directed analysis (EDA) in a risk assessment by identifying the effect-causing contaminants of complex environmental samples.     

Are network indices robust indicators of food web functioning? A Monte Carlo approach

Indices based on network theory are often used to describe food web functioning. These indices take as input food web flows that are estimated based on merging of (scarce) data with linear inverse methods (LIMs). Due to under sampling, most food webs are highly uncertain and can only be quantified within a specific uncertainty range. The linear inverse method (LIM) can estimate food web flows using a variety of techniques, e.g. the parsimonious or minimum norm (MN) solution, which selects one food web, based on a quadratic minimization technique or the Monte Carlo solution where a finitely many random solutions are generated which are then averaged. We use the Monte Carlo approach (MCA) to estimate the values of several indices from four published food webs, the Gulf of Riga for the autumn, summer and spring seasons, and the Takapoto atoll system. We first show that network indices are much better constrained than the uncertain food webs from which they are calculated. Therefore, even in the face of food web uncertainty, they are robust estimators of food web functioning. We then use the MCA-derived network indices to generate cumulative density functions for each index. These serve to compute the probabilities of the MN indices estimates being an extreme solution as compared to the median values. Our findings show that 82% of the MN solutions are smaller than the MCA solutions, and 63% of the network indices are significantly under-estimated.     

Model selection analysis of temporal variation in benefit for an ant-tended treehopper

Recent studies of mutualism have emphasized both that the net benefit to participants depends on the ecological context and that the density-dependent pattern of benefit is key to understanding the population dynamics of mutualism. Indeed, changes in the ecological context are likely to drive changes in both the magnitude of benefit and the density-dependent pattern of benefit. Despite the close linkage between these two areas of research, however, few studies have addressed the factors underlying variation in the density-dependent pattern of benefit. Here I use model selection to evaluate how variation in the benefits of a mutualism drives temporal variation in the density-dependent pattern of net benefit for the ant-tended treehopper Publilia concava. In the interaction between ants and treehoppers in the genus Publilia, ants collect the sugary excretions of treehoppers as a food resource, and treehoppers benefit both directly (e.g., by feeding facilitation) and indirectly (e.g., by predator protection). Results presented here show that temporal changes in the relative magnitude of direct and indirect benefit components of ant tending, especially the effectiveness of predator protection by ants, qualitatively change the overall pattern of density-dependent benefit between years with maximum benefit shifting from treehoppers in small to large aggregations. These results emphasize the need for empirical studies that evaluate the long-term dynamics of mutualism and theoretical studies that consider the population dynamics consequences of variation in the density-dependent pattern of benefit.     

Adaptations in a hierarchical food web of southeastern Lake Michigan

Two issues in ecological network theory are: (1) how to construct an ecological network model and (2) how do entire networks (as opposed to individual species) adapt to changing conditions? We present a novel method for constructing an ecological network model for the food web of southeastern Lake Michigan (USA) and we identify changes in key system properties that are large relative to their uncertainty as this ecological network adapts from one time point to a second time point in response to multiple perturbations. To construct our food web for southeastern Lake Michigan, we followed the list of seven recommendations outlined in Cohen et al. [Cohen, J.E., et al., 1993. Improving food webs. Ecology 74, 252–258] for improving food webs. We explored two inter-related extensions of hierarchical system theory with our food web; the first one was that subsystems react to perturbations independently in the short-term and the second one was that a system's properties change at a slower rate than its subsystems’ properties. We used Shannon's equations to provide quantitative versions of the basic food web properties: number of prey, number of predators, number of feeding links, and connectance (or density). We then compared these properties between the two time-periods by developing distributions of each property for each time period that took uncertainty about the property into account. We compared these distributions, and concluded that non-overlapping distributions indicated changes in these properties that were large relative to their uncertainty. Two subsystems were identified within our food web system structure (p < 0.001). One subsystem had more non-overlapping distributions in food web properties between Time 1 and Time 2 than the other subsystem. The overall system had all overlapping distributions in food web properties between Time 1 and Time 2. These results supported both extensions of hierarchical systems theory. Interestingly, the subsystem with more non-overlapping distributions in food web properties was the subsystem that contained primarily benthic taxa, contrary to expectations that the identified major perturbations (lower phosphorous inputs and invasive species) would more greatly affect the subsystem containing primarily pelagic taxa. Future food-web research should employ rigorous statistical analysis and incorporate uncertainty in food web properties for a better understanding of how ecological networks adapt.     

Chemical communication and coevolution in an ant–plant mutualism

Protective ant–plant interactions provide valuable model systems to study mutualisms. Here, we summarise our recent research on chemical and physiological adaptations that contribute to the stabilisation of the mutualism between Mesoamerican Acacia host plants and their Pseudomyrmex ant inhabitants against exploiters, that is, species using host-derived rewards without rendering a service. Acacia hosts produce food bodies (FBs) and extrafloral nectar (EFN). Both types of reward are chemically adapted to their specific function as ant food and protected from different exploiters. FBs contained higher amounts of specific proteins than the leaves from which they originate. EFN possessed amino acids making it attractive for the mutualist ants and an invertase making its carbohydrate composition nutritionally suitable for the mutualists but unattractive for generalists. Moreover, pathogenesis-related proteins such as glucanases, chitinases and peroxidases were found in EFN, which likely serve as protection from microorganisms. Digestive adaptations were found that make workers of the ant mutualists dependent on the host-derived food sources, a mechanism that likely counteracts the evolution of cheaters. The ants also possessed a high diversity of bacterial associates, several of which appeared involved in nitrogen fixation, thus contributing to the nutrition of these ‘vegetarian’ ants. By contrast, a non-defending ant species that parasitises the host plants appeared physiologically less adapted to the host-derived food rewards; this species, thus, likely is competitively inferior when colony growth is limited by plant-derived rewards. In summary, several physiological adaptations of both host plants and ants stabilise the Acacia–Pseudomyrmex mutualism against exploitation.     

Lethal and nonlethal predator effects on an herbivore guild mediated by system productivity

Indirect effects propagated through intervening species in a food web have important effects on community properties. Traditionally, these indirect effects have been conceptualized as mediated through density changes of the intervening species, but it is becoming increasingly apparent that those mediated through trait (phenotypic) responses also can be very important. Because density- and trait-mediated indirect effects have different properties, it is critical that we understand the mechanisms of transmission in order to predict how they will interact, and when or where they will be important. In this study, we examined the mechanisms and consequences of the lethal (density-mediated) and nonlethal (trait-mediated) effects of a larval odonate predator on a guild of four herbivore species (a larval anuran and three species of snails) and their resources. We also manipulated system productivity in order to explore the effects of environmental context on the transmission of these two types of indirect effects. We show that trait-mediated effects arising from the predator can be very strong relative to density-mediated effects on both the competing herbivores and the species composition and production of their resources. A number of these indirect effects are shown to be contingent on productivity of the system. We further present evidence that trait- and density-mediated indirect effects originating from a predator may be transmitted independently through different routes in a food web, particularly when spatial responses of the transmitting prey are involved. Finally, effects on prey growth due to trait responses to the predator varied from negative to positive in predictable ways as a function of time and indirect effects on the larger food web. These results indicate the important role that trait-mediated indirect effects can play in trophic cascades and keystone predator interactions, and we discuss how the mechanisms involved can be incorporated in theory.     

Simple predator–prey interactions control dynamics in a plankton food web model

A plankton food web model is analysed using interaction parameter values appropriate to the upper mixed layer of the high latitude oceans. The dynamics of this four-variable system are analysed in terms of the dynamics of much simpler two-variable predator–prey subsystems. Thus, the food web's robust, periodic, four-dimensional dynamics are explained by means of two-dimensional spirals and limit cycles. These dynamical subsystems are coupled by means of an omnivore that transfers control of the dynamics between the two predator–prey subsystems. The food web may substantially decouple the predator–prey subsystems so that the oscillating phytoplankton/zooplankton blooms exhibit population collapses when bacterial ‘breathers’ briefly dominate after growing dramatically from low background levels. This regular bloom/breather behaviour becomes benignly chaotic when the system is mildly forced by the annual cycle of the sun's irradiance.     

An inverse model for the analysis of the Venice lagoon food web

A steady-state model of the Venice lagoon food web was constructed, based on a comprehensive set of data, which were collected in the years 2001-2005. Energy flows were estimated by means of an inverse methodology of constrained optimization based on the Minimum Norm criterion, i.e. on the minimization of both the sum of squares of the residuals and of the sum of squares of energy flows. The solution was constrained by a set inequalities, which were derived from general eco-physiological knowledge and site specific data on energy flows. The trophic network was represented by thirty-two nodes, including single-species compartments for the species of high economical or ecological relevance. Mass balance equations were weighted, in order to obtain meaningful results in presence of large differences, up to 5 orders of magnitude, among biomasses. A perturbation technique was applied, with the purpose of reducing the risk of finding solutions heavily affected by the set of constraints and of obtaining a more robust representation of the energy flows. The main patterns of energy flow are consistent with those obtained in previous attempts at modelling the Venice lagoon food web. Micro- and macro-phytobenthos account for the largest fraction of the primary production. Energy is then transferred towards higher trophic levels by means of two main pathways: the recycling of dead biomass through the detritus compartment and the direct consumption by grazers. The first pathway is the most important and accounts for approximately two/thirds of the energy transferred to the second trophic level.