Social network analysis and valid Markov chain Monte Carlo tests of null models

Analyses of animal social networks derived from group-based associations often rely on randomisation methods developed in ecology (Manly, Ecology 76:1109–1115, 1995) and made available to the animal behaviour community through implementation of a pair-wise swapping algorithm by Bejder et al. (Anim Behav 56:719–725, 1998). We report a correctable flaw in this method and point the reader to a wider literature on the subject of null models in the ecology literature. We illustrate the importance of correcting the method using a toy network and use it to make a preliminary analysis of a network of associations among eagle rays.

Ecosystem network analysis indicators are generally robust to parameter uncertainty in a phosphorus model of Lake Sidney Lanier, USA

Understanding how data uncertainty influences ecosystem analysis is critical as we move toward ecosystem-based management. Here, we investigate how 18 Ecological Network Analysis (ENA) indicators that characterize ecosystem growth, development, and condition are affected by uncertainty in an ecosystem model of Lake Sidney Lanier (USA). We applied ENA to 122 plausible parameterizations of the ecosystem developed by Borrett and Osidele (2007, Ecological Modelling 200, 371-387), and then used the coefficient of variation (CV) to compare system indicator variability. We considered Total System Throughput (TST) as a measure of the underlying model uncertainty and tested three hypotheses. First, we hypothesized that non-ratio indicators whose calculation includes the TST would be at least as variable as TST if not more variable. Second, we postulated that indicators calculated as ratios, with TST in the numerator and denominator would tend to be less variable than TST because its influence will cancel. Last, we expected the Average Mutual Information (AMI) to be less variable than TST because it is a bounded function. Our work shows that the 18 indicators grouped into four categories. The first group has significantly larger CVs than the CV for TST. In this group, model uncertainty is amplified rendering these three indicators less useful. The second group of four indicators shows no significant difference in variability with respect to TST. Finally, there are two groups whose CV values are significantly lower than that for TST. The least variable group includes the ratio-based indicators and Average Mutual Information. Due to their low variability, we conclude that these indicators are the most robust to the parameter uncertainty and most useful for ecosystem assessment and comparative ecosystem analysis. In summary, this work suggests that we can be as certain, or more certain, in most of the selected ENA indicators as we are in the parameters of the model analyzed.     

Using a fisheries ecosystem model with a water quality model to explore trophic and habitat impacts on a fisheries stock: A case study of the blue crab population in the Chesapeake Bay

Recent calls for the development of ecosystem-based fisheries management compel the development of resource management tools and linkages between existing fisheries management tools and other resource tools to enable assessment and management of multiple impacts on fisheries resources. In this paper, we describe the use of the Chesapeake Bay Fisheries Ecosystem Model (CBFEM), developed using the Ecopath with Ecosim (EwE) software, and the Chesapeake Bay Water Quality Model (WQM) to demonstrate how linkages between available modeling tools can be used to inform ecosystem-based natural resource management. The CBFEM was developed to provide strategic ecosystem information in support of fisheries management. The WQM was developed to assess impacts on water quality. The CBFEM was indirectly coupled with the WQM to assess the effects of water quality and submerged aquatic vegetation (SAV) on blue crabs. The output from two WQM scenarios (1985-1994), a baseline scenario representing actual nutrient inputs and another with reduced inputs based on a tributary management strategy, was incorporated into the CBFEM. The results suggested that blue crab biomass could be enhanced under management strategies (reduced nutrient input) when the effective search rate of blue crab young-of-the-year's (YOY's) predators or the vulnerability of blue crab YOY to its predators was adjusted by SAV. Such model linkages are important for incorporating physical and biological components of ecosystems in order to explore ecosystem-based fisheries management options.     

Response of balanced network models to large-scale perturbation: Implications for evaluating the role of small pelagics in the Gulf of Maine

Exploring the response of an ecosystem, and subsequent tradeoffs among its biological community, to human perturbations remains a key challenge for the implementation of an ecosystem approaches to fisheries (EAF). To address this and related issues, we developed two network (or energy budget) models, Ecopath and Econetwrk, for the Gulf of Maine ecosystem. These models included 31 network “nodes” or biomass state variables across a broad range of trophic levels, with the present emphasis to particularly elucidate the role of small pelagics. After initial network balancing, various perturbation scenarios were evaluated to explore how potential changes to different fish, fisheries and lower trophic levels can affect model outputs. Categorically across all scenarios and interpretations thereof, there was minimal change at the second trophic levels and most of the “rebalancing” after a perturbation occurred via alteration of the diet matrix. Yet the model results from perturbations to a balanced energy budget fall into one of three categories. First, some model results were intuitive and in obvious agreement with established ecological and fishing theory. Second, some model results were counter-intuitive upon initial observation, seemingly contradictory to known ecological and fishing theory; but upon further examination the results were explainable given the constraints of an equilibrium energy budget. Finally, some results were counter-intuitive and difficult to reconcile with theory or further examination of equilibrium constraints. A detailed accounting of biomass flows for example scenarios explores some of the non-intuitive results more rigorously. Collectively these results imply a need to carefully track biomass flows and results of any given perturbation and to critically evaluate the conditions under which a new equilibrium is obtained for these types of models, which has implications for dynamic simulations based off of them. Given these caveats, the role of small pelagics as a prominent component of this ecosystem remains a robust conclusion. We discuss how one might use this approach in the context of further developing an EAF, recognizing that a more holistic, integrated perspective will be required as we continue to evaluate tradeoffs among marine biological communities.     

Different response of bacterial community to the changes of nutrients and pollutants in sediments from an urban river network

• Bacterial community varied spatially in sediments from the urban river network. • Key environmental factors shaping bacterial community were detected by RDA. • Bacterial co-occurrence networks changed at different levels of nutrient and metal. • Potential indicator species were selected to predict pollution risk in sediment. Microbial communities in sediment are an important indicator linking to environmental pollution in urban river systems. However, how the diversity and structure of bacterial communities in sediments from an urban river network respond to different environmental factors has not been well studied. The goal of this study was to understand the patterns of bacterial communities in sediments from a highly dense urbanized river network in the lower Yangtze River Delta by Illumina MiSeq sequencing. The correlations between bacterial communities, the environmental gradient and geographical distance were analyzed by redundancy analysis (RDA) and network methods. The diversity and richness of bacterial community in sediments increased from upstream to downstream consistently with the accumulation of nutrient in the urban river network. Bacterial community composition and structure showed obvious spatial changes, leading to two distinct groups, which were significantly related to the characteristics of nutrient and heavy metal in sediments. Humic substance, available nitrogen, available phosphorus, Zn, Cu, Hg and As were selected as the key environmental factors shaping the bacterial community in sediments based on RDA. The co-occurrence patterns of bacterial networks showed that positive interaction between bacterial communities increased but the connectivity among bacterial genera and stability of sediment ecosystem reduced under a higher content of nutrient and heavy metal in average. The sensitive and ubiquitous taxa with an overproportional response to key environmental factors were detected as indicator species, which provided a novel method for the prediction of the pollution risk of sediment in an urban river network.     

Which forcing factors fit? Using ecosystem models to investigate the relative influence of fishing and changes in primary productivity on the dynamics of marine ecosystems

Fishing mortality and primary production (or proxy for) were used to drive the dynamics of fish assemblages in 9 trophodynamic models of contrasting marine ecosystems. Historical trends in abundance were reconstructed by fitting model predictions to observations from stock assessments and fisheries independent survey data. The model fitting exercise derives values for otherwise unknown parameters that specify the relative strength of trophic interactions and, in some instances, a time series anomaly for changes in primary production. We measured how much better or worse were model predictions when bottom-up forcing by primary production were added to top-down forcing by fishing. Searching for cross system patterns, the relative contribution of fishing and changes in primary production, mediated through trophic interactions, are evaluated for the ecosystems as a whole and for selected similar species in different ecosystems. The analysis provides a simple qualitative way to explain which forcing factors have most influence on modeled dynamics. Both fishing and primary production forcing were required to obtain the best model fits to data. Fishing effects more strongly influenced 6 of 9 of the ecosystems, but primary production was more often found to be the main factor influencing the selected pelagic and demersal fish stock trends. Examination of sensitivity to ecological and model parameters suggests that the results are the product of complex food-web interactions rather than simple deterministic responses of the models.