Dynamic environ analysis of compartmental systems: A computational approach

Ecosystems are often modeled as stocks of matter or energy connected by flows. Network environ analysis (NEA) is a set of mathematical methods for using powers of matrices to trace energy and material flows through such models. NEA has revealed several interesting properties of flow–storage networks, including dominance of indirect effects and the tendency for networks to create mutually positive interactions between species. However, the applicability of NEA is greatly limited by the fact that it can only be applied to models at constant steady states. In this paper, we present a new, computationally oriented approach to environ analysis called dynamic environ approximation (DEA). As a test of DEA, we use it to compute compartment throughflow in two implementations of a model of energy flow through an oyster reef ecosystem. We use a newly derived equation to compute model throughflow and compare its output to that of DEA. We find that DEA approximates the exact results given by this equation quite closely – in this particular case, with a mean Euclidean error ranging between 0.0008 and 0.21 – which gives a sense of how closely it reproduces other NEA-related quantities that cannot be exactly computed and discuss how to reduce this error. An application to calculating indirect flows in ecosystems is also discussed and dominance of indirect effects in a nonlinear model is demonstrated.     

Analysis of ecosystem structure and function: extended path and flow analysis of a steady-state oyster reef model

An analysis of the extended path and flow structure of a six compartment steady-state oyster reef model was conducted. The extended path and flow structure were analyzed in the context of a refined canonical path classification system based on the systems theory methods of environ and network unfolding analyses. A computer implementation of an operational path classification system facilitated investigation of a finite portion (path length ≤17 arcs) of the direct and indirect path structure of the oyster reef model. Important results of the path structure analysis include: (1) few simple paths and large numbers of compound paths enumerated; (2) dominance of path numbers by subsequent passage terminal cycle paths; (3) structural evidence in support of feedback control in ecosystems; (4) results provide evidence by analogy to support the hypothesis of network homogenization first described using the systems analysis methods of environ analysis and network unfolding; (5) constancy of the pattern of origin–destination path counts with increasing path length; (6) importance of nonliving compartments in the extended path structure of ecosystems. Simultaneous path and flow analysis of the oyster reef model assessed the flow contributions of the fundamental path categories for this model using a modification of a path-based network unfolding method. First passage paths contribute most of the flow; however, multiple passage cyclic paths also provide a large (22%) flow contribution. Because of cycling in the system, the numerous long paths in the extended path structure of this ecosystem model are significant in its function as represented by the flows. These results provide microscopic evidence for the macroscopic results of environ analysis that implicate cycling as a key ecosystem attribute in the mechanisms of holistic system determination. The principles enunciated here for a model with a low cycling index (11%) carry over to, and would be even more significant for, models with high cycling indexes. These results also serve to form a link between the extended structure of food webs and their functioning as represented by energy-matter flows. The present analysis demonstrates that extended path structure, and the component articulation from which it is generated, have significant consequences for ecosystem function.     

Environ indicator sensitivity to flux uncertainty in a phosphorus model of Lake Sidney Lanier, USA

Effective environmental impact assessment and management requires improved understanding of the organization and transformation of ecosystems in which independent agents are linked through an intricate network of energy, matter, and informational interactions. While advances have been made, we still lack a complete understanding of the processes that create, constrain, and sustain ecosystems. Network environ analysis (NEA) provides one approach for building novel ecosystem insights, but it is model dependent. As ecological modeling is an imprecise art, often complicated by inadequate empirical data, the utility of NEA may be limited by model uncertainty. Here, we investigate the sensitivity of NEA indicators of ecosystem growth and development to flow and storage uncertainty in a phosphorus model of Lake Sidney Lanier, USA. The indicators are total system throughflow (TST), total system storage (TSS), total boundary input (Boundary), Finn cycling index (FCI), ratio of indirect-to-direct flows (Indirect/Direct), indirect flow index (IFI), network aggradation (AGG), network homogenization (HMG), and network amplification (AMP). Our results make two primary contributions. First, they demonstrate that five of the indicators – FCI, Indirect/Direct, IFI, AGG and HMG – are relatively robust to the flow and storage uncertainty in the Lake Lanier model. This stability lets us draw robust conclusions about the Lake Lanier ecosystem organization (e.g., phosphorus flux in the lake is dominated by internal processes) in spite of uncertainties in the model. Second, we show that the majority of the indicators co-vary and that most of their common variation could be mapped onto two latent factors, which we interpret as (1) system integration and (2) boundary influences.     

Network mutualism: Positive community-level relations in ecosystems

Empirically observable energy and matter transfers in ecosystems create network structures commonly called food webs. The relation or interaction type associated with each link between pair-wise objects can be classified as (+, −) or (−, +) depending on the net gain or loss experienced by each object. If objects are not adjacent in the food web, then their observed direct interaction is neutralism (0, 0). From this perspective, a zero-sum balance exists between the number of positive and negative relations in the ecosystem. However, community-level relations arise from observable direct and unobservable indirect pathways within a food web, giving rise to indirectly mediated relations, mutualism (+, +) and competition (−, −). Determination of community-level relations requires a systemic or holistic approach. Utility measures from environ analysis in the broader frame of ecological network analysis (ENA) provide such a methodology to investigate the relations resulting from all observed and indirect transfers. This research demonstrates the methodology and shows three important results from the analysis. First, all objects in ecological networks are related either through their input and output environs, and therefore all objects interact with and influence the others in the web: there are no null community-level relations. Second, the community-level relations can and do differ from direct relations: what you see is not always what you get. Third, due to the web of trophic and non-trophic interactions, community-level relations usually have a greater occurrence of mutualism than competition making them more positive than the direct relations that produced them: this is the property called network mutualism.     

Structural changes in lake functioning induced from nutrient loading and climate variability

Climate variability is increasingly recognized as an important regulatory factor, capable of influencing the structural properties of aquatic ecosystems. Lakes appear to be particularly sensitive to the ecological impacts of climate variability, and several long time series have shown a close coupling between climate, lake thermal properties and individual organism physiology, population abundance, community structure, and food web dynamics. Thus, understanding the complex interplay among meteorological forcing, hydrological variability, and ecosystem functioning is essential for improving the credibility of model-based water resources/fisheries management. Our objective herein is to examine the relative importance of the ecological mechanisms underlying plankton seasonal variability in Lake Washington, Washington State (USA), over a 35-year period (1964–1998). Our analysis is founded upon an intermediate complexity plankton model that is used to reproduce the limiting nutrient (phosphate)–phytoplankton–zooplankton–detritus (particulate phosphorus) dynamics in the lake. Model parameterization is based on a Bayesian calibration scheme that offers insights into the degree of information the data contain about model inputs and allows obtaining predictions along with uncertainty bounds for modeled output variables. The model accurately reproduces the key seasonal planktonic patterns in Lake Washington and provides realistic estimates of predictive uncertainty for water quality variables of environmental management interest. A principal component analysis of the annual estimates of the underlying ecological processes highlighted the significant role of the phosphorus recycling stemming from the zooplankton excretion on the planktonic food web variability. We also identified a moderately significant signature of the local climatic conditions (air temperature) on phytoplankton growth (r = 0.41), herbivorous grazing (r = 0.38), and detritus mineralization (r = 0.39). Our study seeks linkages with the conceptual food web model proposed by Hampton et al. [Hampton, S.E., Scheuerell, M.D., Schindler, D.E., 2006b. Coalescence in the Lake Washington story: interaction strengths in a planktonic food web. Limnol. Oceanogr. 51, 2042–2051.] to emphasize the “bottom-up” control of the Lake Washington plankton phenology. The posterior predictive distributions of the plankton model are also used to assess the exceedance frequency and confidence of compliance with total phosphorus (15 μg L⁻¹) and chlorophyll a (4 μg L⁻¹) threshold levels during the summer-stratified period in Lake Washington. Finally, we conclude by underscoring the importance of explicitly acknowledging the uncertainty in ecological forecasts to the management of freshwater ecosystems under a changing global environment.     

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.     

Nutrient enrichment reduces constraints on material flows in a detritus-based food web

Most aquatic and terrestrial ecosystems are experiencing increased nutrient availability, which is affecting their structure and function. By altering community composition and productivity of consumers, enrichment can indirectly cause changes in the pathways and magnitude of material flows in food webs. These changes, in turn, have major consequences for material storage and cycling in the ecosystem. Understanding mechanisms and predicting consequences of nutrient-induced changes in material flows requires a quantitative food web approach that combines information on consumer energetics and consumer-resource stoichiometry. We examined effects of a whole-system experimental nutrient enrichment on the trophic basis of production and the magnitude and pathways of carbon (C), nitrogen (N), and phosphorus (P) flows in a detritus-based stream food web. We compared the response of the treated stream to an adjacent reference stream throughout the study. Dietary composition and elemental flows varied considerably among invertebrate functional feeding groups. During nutrient enrichment, increased flows of leaf litter and amorphous detritus to shredders and gatherers accounted for most of the altered flows of C from basal resources to consumers. Nutrient enrichment had little effect on patterns of material flows but had large positive effects on the magnitude of C, N, and P flows to consumers (mean increase of 97% for all elements). Nutrient-specific food webs revealed similar flows of N and P to multiple functional groups despite an order of magnitude difference among groups in consumption of C. Secondary production was more strongly related to consumption of nutrients than C, and increased material flows were positively related to the degree of consumer-resource C:P and C:N imbalances. Nutrient enrichment resulted in an increased proportion of detrital C inputs consumed by primary consumers (from -15% to 35%) and a decreased proportion of invertebrate prey consumed by predators (from -80% to 55%). Our results demonstrate that nutrient enrichment of detritus-based systems may reduce stoichiometric constraints on material flows, increase the contribution of consumers to C, N, and P cycling, alter the proportion of C inputs metabolized by consumers, and potentially lead to reduced ecosystem-level storage of C.     

Do microbial processes regulate the stability of a coral atoll's enclosed pelagic ecosystem?

Complex marine ecosystems contain multiple feedback cycles that can cause unexpected responses to perturbations. To better predict these responses, complicated models are increasingly being developed to enable the study of feedback cycles. However, the sparseness of ecological data often limits the direct empirical parameterization of all model parameters. Here we use a Bayesian inverse analysis approach to synthesize empirical data and ecological theory derived from published studies of a coral atoll's enclosed pelagic ecosystem (Takapoto Atoll, French Polynesia). We then use the estimates of flux magnitudes to parameterize probabilistic compartment models with two forms of heterotrophic consumption: (1) “bottom-up” donor-controlled heterotrophic consumption and (2) “top-down” mass-action heterotrophic consumption. We explore how the flux magnitudes affect the ecosystem's stability properties of resilience, reactivity, and resistance under both assumptions for heterotrophic consumption. The models suggest that the microbial uptake of dissolved organic carbon (DOC) regulates the long term rate of return to steady state following a temporary or pulse perturbation (resilience), and the cycling of carbon between abiotic pools and heterotrophic compartments regulates the short-term response (reactivity). In the bottom-up process model, the sensitivity of steady state masses following a sustained or press perturbation (resistance) is highest for the DOC pool following a sustained change to the microbial uptake rate of DOC. Further, a change in the microbial uptake of DOC propagates through the ecosystem and affects the steady state values of zooplankton. The analysis suggests that the food web is highly dependent on the recycling between the abiotic and biotic carbon pools, particularly as mediated by the microbial consumption of DOC, and this recycling determines how the ecosystem responds to perturbations.     

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.     

Weak indirect effects inherent to nitrogen biogeochemical cycling within anthropogenic ecosystems: A network environ analysis

Indirect effects are assumed to be the major causes of the complexity and stability of ecological networks. The complexity of urban-rural complexes (URCs) could also be attributed to the indirect effects associated with human activities. No studies, however, have quantified the strength of indirect effects in relation to urban biogeochemistry. A network environ analysis (NEA) was used for this study to investigate and compare indirect effects in relation to the nitrogen (N) cycling networks of 22 natural ecosystems and five URCs. Results show that indirect effects were proven to be weak for URC N cycling networks (accounting for only ∼2% of the overall effects measured in natural ecosystems). The weak indirect effects found provide a counterexample for the hypothesis that indirect effects are in fact the dominant components of biogeochemical networks. It also implies that human activity in itself does not always raise the complexity of ecological processes as previously suggested. Weak indirect effects also lead to perturbation fragility for URC N cycles (where the decay rate is greater in comparison to natural ecosystems by a factor of 13). In order to improve the robustness and efficiency of URC biogeochemical cycling, a knockout analysis was carried out. By comparing results after removing single interactions between natural ecosystems and URCs it was found that the loss of indirect effects require cooperative strategies to optimize N cycling networks within URCs.     

Disentangling trophic relationships in a High Arctic tundra ecosystem through food web modeling

Determining the manner in which food webs will respond to environmental changes is difficult because the relative importance of top-down vs. bottom-up forces in controlling ecosystems is still debated. This is especially true in the Arctic tundra where, despite relatively simple food webs, it is still unclear which forces dominate in this ecosystem. Our primary goal was to assess the extent to which a tundra food web was dominated by plant-herbivore or predator-prey interactions. Based on a 17-year (1993-2009) study of terrestrial wildlife on Bylot Island, Nunavut, Canada, we developed trophic mass balance models to address this question. Snow Geese were the dominant herbivores in this ecosystem, followed by two sympatric lemming species (brown and collared lemmings). Arctic foxes, weasels, and several species of birds of prey were the dominant predators. Results of our trophic models encompassing 19 functional groups showed that <10% of the annual primary production was consumed by herbivores in most years despite the presence of a large Snow Goose colony, but that 20-100% of the annual herbivore production was consumed by predators. The impact of herbivores on vegetation has also weakened over time, probably due to an increase in primary production. The impact of predators was highest on lemmings, intermediate on passerines, and lowest on geese and shorebirds, but it varied with lemming abundance. Predation of collared lemmings exceeded production in most years and may explain why this species remained at low density. In contrast, the predation rate on brown lemmings varied with prey density and may have contributed to the high-amplitude, periodic fluctuations in the abundance of this species. Our analysis provided little evidence that herbivores are limited by primary production on Bylot Island. In contrast, we measured strong predator-prey interactions, which supports the hypothesis that this food web is primarily controlled by top-down forces. The presence of allochthonous resources subsidizing top predators and the absence of large herbivores may partly explain the predominant role of predation in this low-productivity ecosystem.     

Quantity and quality: unifying food web and ecosystem perspectives on the role of resource subsidies in freshwaters

Although the study of resource subsidies has emerged as a key topic in both ecosystem and food web ecology, the dialogue over their role has been limited by separate approaches that emphasize either subsidy quantity or quality. Considering quantity and quality together may provide a simple, but previously unexplored, framework for identifying the mechanisms that govern the importance of subsidies for recipient food webs and ecosystems. Using a literature review of > 90 studies of open-water metabolism in lakes and streams, we show that high-flux, low-quality subsidies can drive freshwater ecosystem dynamics. Because most of these ecosystems are net heterotrophic, allochthonous inputs must subsidize respiration. Second, using a literature review of subsidy quality and use, we demonstrate that animals select for high-quality food resources in proportions greater than would be predicted based on food quantity, and regardless of allochthonous or autochthonous origin. This finding suggests that low-flux, high-quality subsidies may be selected for by animals, and in turn may disproportionately affect food web and ecosystem processes (e.g., animal production, trophic energy or organic matter flow, trophic cascades). We then synthesize and review approaches that evaluate the role of subsidies and explicitly merge ecosystem and food web perspectives by placing food web measurements in the context of ecosystem budgets, by comparing trophic and ecosystem production and fluxes, and by constructing flow food webs. These tools can and should be used to address future questions about subsidies, such as the relative importance of subsidies to different trophic levels and how subsidies may maintain or disrupt ecosystem stability and food web interactions.     

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.     

Carbon-based balanced trophic structure and flows in the offshore Lake Ontario food web before (1987-1991) and after (2001-2005) invasion-induced ecosystem change

Replicate mass-balanced solutions to Ecopath models describing carbon-based trophic structures and flows were developed for the Lake Ontario offshore food web before and after invasion-induced disruption. The food webs link two pathways of energy and matter flow: the grazing chain (phytoplankton-zooplankton-fish) and the microbial loop (bacteria-protozoans) and include 19 species-groups and three detrital groups. Mass-balance was achieved by using constrained optimization techniques to randomly vary initial estimates of biomass and diet composition. After the invasion, production declined for all trophic levels and species-groups except Chinook salmon. The trophic level (TL) increased for smelt, adult sculpin, adult alewife and Chinook salmon. Changes to ecotrophic efficiencies indicate a reduction in phytoplankton grazing, increased predation pressure on Mysis, adult smelt and alewife and decreased predation pressure on protozoans. Specific resource to consumer TTE changed; increasing for protozoans (8.0-11.5%), Mysis (0.6-1.0%), and Chinook salmon (1.0-2.3%) and other salmonines (0.4-0.5%) and decreasing for zooplankton (20.2-15.1%), prey-fish (9.7-8.8%), and benthos (1.7-0.6%). Direct trophic influences of recent invasive species were low. The synchrony of the decline in PP and species-group production indicates strong bottom-up influence. Mass balance required an increase of two to threefold in lower trophic level biomass and production, confirming a previously observed paradoxical deficit in lower trophic level production. Analysis of food web changes suggest hypotheses that may apply to other similar large pelagic systems including, (1) as pelagic primary productivity declines, overgrazing of zooplankton results in an increase in protozoan production and a loss of trophic transfer efficiency, (2) habitat and food web changes increased Mysis predation on Diporeia and contributed to their recent decline, and (3) production of Chinook salmon, the primary piscivore, was uncoupled from pelagic production processes. This study demonstrates the value of food web models to better understand the impact of invasive species and to develop novel hypotheses concerning trophic influences.     

Integrative trophic network assessments of a lentic ecosystem by key ecological approaches of water chemistry, trophic guilds, and ecosystem health assessments along with an ECOPATH model

The objective of this study was to describe the trophic structure and energy flow in a lentic ecosystem in South Korea. Physicochemical water conditions were evaluated along with the reservoir ecosystem health using a multimetric IBI model. Nutrient analyses of the reservoir showed a nutrient rich and hypereutrophic system. Guild analysis revealed that tolerant and omnivorous species dominated the ecosystem. Tolerant fish, as a proportion of the number of individuals, were associated (R² > 0.90, p < 0.01) with TN and TP, the key indicators of trophic state in lentic ecosystems. The mean Reservoir Ecosystem Health Assessment (REHA) score was 19.3 during the study, which was judged as in ‘fair to poor’ condition. A trophic analysis of the reservoir estimated by the ECOPATH model shows that most activity in terms of energy flow occurred in the lower part of the trophic web, where there was intensive use of primary producers as a food source. Consequently, of the 10 consumer groups, nine fell within trophic levels <2.8. Trophic levels (TL) estimated from the weighted average of prey trophic levels varied from 1.0 for phytoplankton, macrophytes, and detritus to 3.25 for the top predator, Pseudobagrus fulvidraco. Our integrated approach to trophic network analysis may provide a key tool for determining the effects of nutrient influx on energy flow pathways in lentic ecosystems.     

Implications of network particle tracking (NPT) for ecological model interpretation

Network particle tracking (NPT), building on the foundation of network environ analysis (NEA), is a new development in the definition of coherence relations within and between connected systems. This paper evaluates three ecosystem models in a comparison of throughflow- and storage-based NEA and NPT. Compartments in models with high indirect effects and Finn cycling showed low correlation of NEA storage and throughflow with particle repeat visits and numbers of particles in compartments at steady state. Conversely, the correlation between NEA and NPT results was high with two models having lower indirect effects and Finn cycling. Analysis of ecological orientors associated with NEA showed NPT to fully support conventional NEA results when the common conditions of donor control and steady state are satisfied. Particle trajectories are recorded in the new concept of a particle “passport”. Ability to track and record particle in-system histories enables views of multiple scales and opens the possibility of making pathway-dependent modeling decisions. NPT may also enable modeling of time, allowing integration of Newtonian, organismal and stochastic modeling perspectives in a single comprehensive analysis.     

Dealing with microtopography of an ombrotrophic bog for simulating ecosystem-level CO2 exchanges

Peatlands contain approximately 25% of the global soil carbon (C), despite covering only 3% of the earth's land surface. In order to evaluate the role of peatlands in global C cycling, models of ecosystem biogeochemistry are required, but peatland ecosystems present a number of unique challenges, particularly how to deal with the large variability that occurs at scales of one to several metres. In models, spatial variability is considered either explicitly for each individual unit and the outputs averaged, referred to as flux upscaling, or implicitly by weighting model parameters by the fractional occurrence of the individual units, referred to as parameter upscaling. The advantage of parameter upscaling is that it is much more computationally efficient: a requirement for hemispheric scale simulations. In this study we determined the differences between modelling a raised bog peatland with hummock-hollow microtopography using flux and parameter upscaling. We used the McGill Wetland Model (MWM), a process-based ecosystem C model for peatlands, configured for hummocks and hollows separately and then a weighted mixture of both. The simulated output based on flux and parameter upscaling was compared with eddy-covariance tower measurements. We found that net ecosystem production (NEP) for hollows was much larger than that for hummocks because total ecosystem respiration (TER) for hummocks was greater while gross primary production (GPP) did not differ significantly between the two topographic features. However, despite differences in components of NEP between hummocks and hollows, there was no statistically significant difference between the NEP based on flux and parameter upscaling using the MWM. Both flux and parameter upscaling show equivalent capability to capture the magnitude, direction, seasonality and inter-annual variability. The root-mean-square-errors (RMSE) are 0.66, 0.45, and 0.49 g C m⁻² day⁻¹, respectively for GPP, TER and NEP based on the flux upscaling, while 0.67, 0.44, and 0.48 g C m⁻² day⁻¹, respectively based on the parameter upscaling. The degree of agreement (d*) is 0.96, 0.97, and 0.88, respectively for GPP, TER and NEP based on the flux upscaling, while 0.96, 0.97, and 0.89, respectively based on the parameter upscaling. This result suggests that differences in processes caused by peatland microtopography scale linearly, which means an ecosystem-level model set-up (i.e. parameter upscaling scheme), is sufficient to simulate the C cycling.     

Integrated model projections of climate change impacts on a North American lake

Climate change is likely to impact terrestrial and aquatic ecosystems via numerous physical and biological mechanisms. This study outlines a framework for projecting potential impacts of climate change on lakes using linked environmental models. Impacts of climate drivers on catchment hydrology and thermal balance in Onondaga Lake (New York State) are simulated using mechanistic models HSPF and UFILS4. Outputs from these models are fed into a lake ecosystem model, developed in AQUATOX. Watershed simulations project increases in the magnitude of peak flows and consequent increases in catchment nutrient export as the magnitude of extreme precipitation events increases. This occurs concurrently with a decrease in annual stream discharge as a result of increased evapotranspiration. Simulated lake water temperatures increase by as much as 5 °C during the 2040-2069 time period, accompanied by a prolonging of the duration of summer stratification. Projected changes include shifts in the timing of nutrient cycling between lake sediments and water column. Plankton taxa projected to thrive under climate change include green algae and Bosmina longirostris. Responses for species at higher trophic levels are mixed. Benthic macroinvertebrates may either prosper (zebra mussels) or decline (chironomids), while fish (e.g., gizzard shad) exhibit high seasonal variability without any clear trend.     

Elucidation of ecosystem attributes of an oligotrophic lake in Hokkaido,Japan, using Ecopath with Ecosim (EwE)

The fishing practices in the oligotrophic Lake Toya, Hokkaido, Japan, have profound implications in the ecosystem sustainability. The status of the sockeye salmon (Oncorhynchus nerka) population has become a serious concern among the lake managers and policy makers during the last decades. While the decline of the sockeye salmon population has been well documented in Lake Toya, there is considerable uncertainty with regards to the impact on the broader system dynamics. In this study, our objective is to address this knowledge gap by undertaking a synthesis of the Lake Toya food web using the mass-balance modeling software Ecopath with Ecosim (EwE). Our primary research question is to examine the repercussions of the declining sockeye salmon population on the trophic dynamics of the lake. Namely, we assess if there are any competing species that might have benefited from the decrease of sockeye salmon standing biomass and to what extent do these changes propagate through the Lake Toya food web? Our analysis pinpoints the critical role of the Japanese smelt (Hypomesus transpacificus nipponensis) in the system, which demonstrates a wide range of effects on several functional groups at both higher and lower trophic levels in Lake Toya. In particular, being a substantial portion of the masu salmon (Oncorhynchus masou) and adult sockeye salmon diets, the Japanese smelt has a positive impact on the top predators of the system. Amphipods, insects, and shrimp strongly benefit from the autochthonous and allochthonous organic matter in the system, while the tight coupling between phytoplankton and zooplankton seems to be particularly critical for the integrity of the Lake Toya food web. Whereas the values of the different ecosystem attributes (e.g., primary production/biomass, biomass/total throughput, system omnivory index, amount of recycled throughput, Finn's cycling index) provide evidence that Lake Toya is an immature system, we note that the internal redundancy and the system overhead estimates suggest that the lake possesses substantial reserves to overcome external perturbations. We also examined the effects of a variety of fishing policies on the biomass of masu salmon and adult sockeye salmon, which verify the belief that the adult sockeye population is quite fragile with high likelihood to collapse. Our analysis also predicts that sockeye will not rebound unless the fishing pressure exerted is substantially reduced (>50% of the reference levels used). Masu salmon seems to benefit under all the scenarios examined indicating that the intensity of the current fishing activities is significantly lower than its biomass accumulation rate in the system.     

From individuals to ecosystem function: toward an integration of evolutionary and ecosystem ecology

An important goal in ecology is developing general theory on how the species composition of ecosystems is related to ecosystem properties and functions. Progress on this front is limited partly because of the need to identify mechanisms controlling functions that are common to a wide range of ecosystem types. We propose that one general mechanism, rooted in the evolutionary ecology of all species, is adaptive foraging behavior in response to predation risk. To support our claim, we present two kinds of empirical evidence from plant-based and detritus-based food chains of terrestrial and aquatic ecosystems. The first kind comes from experiments that explicitly trace how adaptive foraging influences ecosystem properties and functions. The second kind comes from a synthesis of studies that individually examine complementary components of particular ecosystems that together provide an integrated perspective on the link between adaptive foraging and ecosystem function. We show that the indirect effects of predators on plant diversity, plant productivity, nutrient cycling, trophic transfer efficiencies, and energy flux caused by consumer foraging shifts in response to risk are qualitatively different from effects caused by reductions in prey density due to direct predation. We argue that a perspective of ecosystem function that considers effects of consumer behavior in response to predation risk will broaden our capacity to explain the range of outcomes and contingencies in trophic control of ecosystems. This perspective also provides an operational way to integrate evolutionary and ecosystem ecology, which is an important challenge in ecology.