Distinct responses to warming within picoplankton communities across an environmental gradient

Picophytoplankton are a ubiquitous component of marine plankton communities and are expected to be favored by global increases in seawater temperature and stratification associated with climate change. Eukaryotic and prokaryotic picophytoplankton have distinct ecology, and global models predict that the two groups will respond differently to future climate scenarios. At a nearshore observatory on the Northeast US Shelf, however, decades of year-round monitoring have shown these two groups to be highly synchronized in their responses to environmental variability. To reconcile the differences between regional and global predictions for picophytoplankton dynamics, we here investigate the picophytoplankton community across the continental shelf gradient from the nearshore observatory to the continental slope. We analyze flow cytometry data from 22 research cruises, comparing the response of picoeukaryote and Synechococcus communities to environmental variability across time and space. We find that the mechanisms controlling picophytoplankton abundance differ across taxa, season, and distance from shore. Like the prokaryote, Synechococcus, picoeukaryote division rates are limited nearshore by low temperatures in winter and spring, and higher temperatures offshore lead to an earlier spring bloom. Unlike Synechococcus, picoeukaryote concentration in summer decreases dramatically in offshore surface waters and exhibits deeper subsurface maxima. The offshore picoeukaryote community appears to be nutrient limited in the summer and subject to much greater loss rates than Synechococcus. This work both produces and demonstrates the necessity of taxon- and site-specific knowledge for accurately predicting the responses of picophytoplankton to ongoing environmental change.

Predator switching strength controls stability in diamond-shaped food web models

In food web models that include more than one prey type for a single predator, it is common for the predator's functional response to include some form of switching-preferential consumption of more abundant prey types. Predator switching promotes coexistence among competing prey types and increases diversity in the prey community. Here, we show how the dynamics of a diamond-shaped food web model of a marine plankton community are sensitive to a parameter that sets the strength of predator switching. Stronger switching destabilizes the model's coexistence equilibrium and leads to the appearance of limit cycles. Stronger switching also increases the evenness of the asymptotic prey community and promotes synchrony in the dynamics of disparate prey types. Given the dependence of model behavior on the strength of predator switching, it is important that modelers carefully consider the parameterization of functional responses that include switching.

The Contributions of Maternal Age Heterogeneity to Variance in Lifetime Reproductive Output

AbstractVariance among individuals in fitness components reflects both genuine heterogeneity between individuals and stochasticity in events experienced along the life cycle. Maternal age represents a form of heterogeneity that affects both the mean and the variance of lifetime reproductive output (LRO). Here, we quantify the relative contribution of maternal age heterogeneity to the variance in LRO using individual-level laboratory data on the rotifer Brachionus manjavacas to parameterize a multistate age × maternal age matrix model. In B. manjavacas, advanced maternal age has large negative effects on offspring survival and fertility. We used multistate Markov chains with rewards to quantify the contributions to variance in LRO of heterogeneity and of the stochasticity inherent in the outcomes of probabilistic transitions and reproductive events. Under laboratory conditions, maternal age heterogeneity contributes 26% of the variance in LRO. The contribution changes when mortality and fertility are reduced to mimic more ecologically relevant environments. Over the parameter space where populations are near stationarity, maternal age heterogeneity contributes an average of 3% of the variance. Thus, the contributions of maternal age heterogeneity and individual stochasticity can be expected to depend strongly on environmental conditions; over most of the parameter space, the variance in LRO is dominated by stochasticity.

Seasonal environmental variability drives microdiversity within a coastal Synechococcus population

Marine microbes often show a high degree of physiological or ecological diversity below the species level. This microdiversity raises questions about the processes that drive diversification and permit coexistence of diverse yet closely related marine microbes, especially given the theoretical efficiency of competitive exclusion. Here, we provide insight with an 8‐year time series of diversity within Synechococcus, a widespread and important marine picophytoplankter. The population of Synechococcus on the Northeast U.S. Shelf is comprised of six main types, each of which displays a distinct and consistent seasonal pattern. With compositional data analysis, we show that these patterns can be reproduced with a simple model that couples differential responses to temperature and light with the seasonal cycle of the physical environment. These observations support the hypothesis that temporal variability in environmental factors can maintain microdiversity in marine microbial populations. We also identify how seasonal diversity patterns directly determine overarching Synechococcus population abundance features.

Political economy of renewable resource federalism

The authority to manage natural capital often follows political boundaries rather than ecological. This mismatch can lead to unsustainable outcomes, as spillovers from one management area to the next may create adverse incentives for local decision making, even within a single country. At the same time, one-size-fits-all approaches of federal (centralized) authority can fail to respond to state (decentralized) heterogeneity and can result in inefficient economic or detrimental ecological outcomes. Here we utilize a spatially explicit coupled natural-human system model of a fishery to illuminate trade-offs posed by the choice between federal vs. state control of renewable resources. We solve for the dynamics of fishing effort and fish stocks that result from different approaches to federal management that vary in terms of flexibility. Adapting numerical methods from engineering, we also solve for the open-loop Nash equilibrium characterizing state management outcomes, where each state anticipates and responds to the choices of the others. We consider traditional federalism questions (state vs. federal management) as well as more contemporary questions about the economic and ecological impacts of shifting regulatory authority from one level to another. The key mechanisms behind the trade-offs include whether differences in local conditions are driven by biological or economic mechanisms; degree of flexibility embedded in the federal management; the spatial and temporal distribution of economic returns across states; and the status-quo management type. While simple rules-of-thumb are elusive, our analysis reveals the complex political economy dimensions of renewable resource federalism.

Dynamics and functional diversity of the smallest phytoplankton on the Northeast US Shelf

Picophytoplankton are the most abundant primary producers in the ocean. Knowledge of their community dynamics is key to understanding their role in marine food webs and global biogeochemical cycles. To this end, we analyzed a 16-y time series of observations of a phytoplankton community at a nearshore site on the Northeast US Shelf. We used a size-structured population model to estimate in situ division rates for the picoeukaryote assemblage and compared the dynamics with those of the picocyanobacteria Synechococcus at the same location. We found that the picoeukaryotes divide at roughly twice the rate of the more abundant Synechococcus and are subject to greater loss rates (likely from viral lysis and zooplankton grazing). We describe the dynamics of these groups across short and long timescales and conclude that, despite their taxonomic differences, their populations respond similarly to changes in the biotic and abiotic environment. Both groups appear to be temperature limited in the spring and light limited in the fall and to experience greater mortality during the day than at night. Compared with Synechococcus, the picoeukaryotes are subject to greater top-down control and contribute more to the region's primary productivity than their standing stocks suggest.

Seasons of Syn

Synechococcus is a widespread and important marine primary producer. Time series provide critical information for identifying and understanding the factors that determine abundance patterns. Here, we present the results of analysis of a 16-yr hourly time series of Synechococcus at the Martha's Vineyard Coastal Observatory, obtained with an automated, in situ flow cytometer. We focus on understanding seasonal abundance patterns by examining relationships between cell division rate, loss rate, cellular properties (e.g., cell volume, phycoerythrin fluorescence), and environmental variables (e.g., temperature, light). We find that the drivers of cell division vary with season; cells are temperature-limited in winter and spring, but light-limited in the fall. Losses to the population also vary with season. Our results lead to testable hypotheses about Synechococcus ecophysiology and a working framework for understanding the seasonal controls of Synechococcus cell abundance in a temperate coastal system.

Advancing an interdisciplinary framework to study seed dispersal ecology

Although dispersal is generally viewed as a crucial determinant for the fitness of any organism, our understanding of its role in the persistence and spread of plant populations remains incomplete. Generalizing and predicting dispersal processes are challenging due to context dependence of seed dispersal, environmental heterogeneity and interdependent processes occurring over multiple spatial and temporal scales. Current population models often use simple phenomenological descriptions of dispersal processes, limiting their ability to examine the role of population persistence and spread, especially under global change. To move seed dispersal ecology forward, we need to evaluate the impact of any single seed dispersal event within the full spatial and temporal context of a plant's life history and environmental variability that ultimately influences a population's ability to persist and spread. In this perspective, we provide guidance on integrating empirical and theoretical approaches that account for the context dependency of seed dispersal to improve our ability to generalize and predict the consequences of dispersal, and its anthropogenic alteration, across systems. We synthesize suitable theoretical frameworks for this work and discuss concepts, approaches and available data from diverse subdisciplines to help operationalize concepts, highlight recent breakthroughs across research areas and discuss ongoing challenges and open questions. We address knowledge gaps in the movement ecology of seeds and the integration of dispersal and demography that could benefit from such a synthesis. With an interdisciplinary perspective, we will be able to better understand how global change will impact seed dispersal processes, and potential cascading effects on plant population persistence, spread and biodiversity.

Intraguild predation enables coexistence of competing phytoplankton in a well-mixed water column

Resource competition theory predicts that when two species compete for a single, finite resource, the better competitor should exclude the other. However, in some cases, weaker competitors can persist through intraguild predation, that is, by eating their stronger competitor. Mixotrophs, species that meet their carbon demand by combining photosynthesis and phagotrophic heterotrophy, may function as intraguild predators when they consume the phototrophs with which they compete for light. Thus, theory predicts that mixotrophy may allow for coexistence of two species on a single limiting resource. We tested this prediction by developing a new mathematical model for a unicellular mixotroph and phytoplankter that compete for light, and comparing the model's predictions with a laboratory experimental system. We find that, like other intraguild predators, mixotrophs can persist when an ecosystem is sufficiently productive (i.e., the supply of the limiting resource, light, is relatively high), or when species interactions are strong (i.e., attack rates and conversion efficiencies are high). Both our mathematical and laboratory models show that, depending upon the environment and species traits, a variety of equilibrium outcomes, ranging from competitive exclusion to coexistence, are possible.

Physiological and ecological drivers of early spring blooms of a coastal phytoplankter

Climate affects the timing and magnitude of phytoplankton blooms that fuel marine food webs and influence global biogeochemical cycles. Changes in bloom timing have been detected in some cases, but the underlying mechanisms remain elusive, contributing to uncertainty in long-term predictions of climate change impacts. Here we describe a 13-year hourly time series from the New England shelf of data on the coastal phytoplankter Synechococcus, during which the timing of its spring bloom varied by 4 weeks. We show that multiyear trends are due to temperature-induced changes in cell division rate, with earlier blooms driven by warmer spring water temperatures. Synechococcus loss rates shift in tandem with division rates, suggesting a balance between growth and loss that has persisted despite phenological shifts and environmental change.

Sex difference and Allee effects shape the dynamics of sex-structured invasions

The rate at which a population grows and spreads can depend on individual behaviour and interactions with others. In many species with two sexes, males and females differ in key life-history traits (e.g. growth, survival and dispersal), which can scale up to affect population rates of growth and spread. In sexually reproducing species, the mechanics of locating mates and reproducing successfully introduce further complications for predicting the invasion speed (spread rate), as both can change nonlinearly with density. Most models of population spread are based on one sex, or include limited aspects of sex differences. Here we ask whether and how the dynamics of finding mates interact with sex-specific life-history traits to influence the rate of population spread. We present a hybrid approach for modelling invasions of populations with two sexes that links individual-level mating behaviour (in an individual-based model) to population-level dynamics (in an integrodifference equation model). We find that limiting the amount of time during which individuals can search for mates causes a demographic Allee effect which can slow, delay, or even prevent an invasion. Furthermore, any sex-based asymmetries in life history or behaviour (skewed sex ratio, sex-biased dispersal, and sex-specific mating behaviours) amplify these effects. In contrast, allowing individuals to mate more than once ameliorates these effects, enabling polygynandrous populations to invade under conditions where monogamously mating populations would fail to establish. We show that details of individuals' mating behaviour can impact the rate of population spread. Based on our results, we propose a stricter definition of a mate-finding Allee effect, which is not met by the commonly used minimum mating function. Our modelling approach, which links individual- and population-level dynamics in a single model, may be useful for exploring other aspects of individual behaviour that are thought to impact the rate of population spread.

Generational Spreading Speed and the Dynamics of Population Range Expansion

Some of the most fundamental quantities in population ecology describe the growth and spread of populations. Population dynamics are often characterized by the annual rate of increase, λ, or the generational rate of increase, R0. Analyses involving R0 have deepened our understanding of disease dynamics and life-history complexities beyond that afforded by analysis of annual growth alone. While range expansion is quantified by the annual spreading speed, a spatial analog of λ, an R0-like expression for the rate of spread is missing. Using integrodifference models, we derive the appropriate generational spreading speed for populations with complex (stage-structured) life histories. The resulting measure, relevant to locations near the expanding edge of a (re)colonizing population, incorporates both local population growth and explicit spatial dispersal rather than solely growth across a population, as is the case for R0. The calculations for generational spreading speed are often simpler than those for annual spreading speed, and analytic or partial analytic solutions can yield insight into the processes that facilitate or slow a population's spatial spread. We analyze the spatial dynamics of green crabs, sea otters, and teasel as examples to demonstrate the flexibility of our methods and the intuitive insights that they afford.

Habitat damage, marine reserves, and the value of spatial management

The biological benefits of marine reserves have garnered favor in the conservation community, but "no-take" reserve implementation is complicated by the economic interests of fishery stakeholders. There are now a number of studies examining the conditions under which marine reserves can provide both economic and ecological benefits. A potentially important reality of fishing that these studies overlook is that fishing can damage the habitat of the target stock. Here, we construct an equilibrium bioeconomic model that incorporates this habitat damage and show that the designation of marine reserves, coupled with the implementation of a tax on fishing effort, becomes both biologically and economically favorable as habitat sensitivity increases. We also study the effects of varied degrees of spatial control on fisheries management. Together, our results provide further evidence for the potential monetary and biological value of spatial management, and the possibility of a mutually beneficial resolution to the fisherman-conservationist marine reserve designation dilemma.

Sex-biased dispersal and the speed of two-sex invasions

Population models that combine demography and dispersal are important tools for forecasting the spatial spread of biological invasions. Current models describe the dynamics of only one sex (typically females). Such models cannot account for the sex-related biases in dispersal and mating behavior that are typical of many animal species. In this article, we construct a two-sex integrodifference equation model that overcomes these limitations. We derive an explicit formula for the invasion speed from the model and use it to show that sex-biased dispersal may significantly increase or decrease the invasion speed by skewing the operational sex ratio at the invasion's low-density leading edge. Which of these possible outcomes occurs depends sensitively on complex interactions among the direction of dispersal bias, the magnitude of bias, and the relative contributions of females and males to local population growth.

A model for energetics and bioaccumulation in marine mammals with applications to the right whale

We present a dynamic energy budget (DEB) model for marine mammals, coupled with a pharmacokinetic model of a lipophilic persistent toxicant. Inputs to the model are energy availability and lipid-normalized toxicant concentration in the environment. The model predicts individual growth, reproduction, bioaccumulation, and transfer of energy and toxicant from mothers to their young. We estimated all model parameters for the right whale; with these parameters, reduction in energy availability increases the age at first parturition, increases intervals between reproductive events, reduces the organisms' ability to buffer seasonal fluctuations, and increases its susceptibility to temporal shifts in the seasonal peak of energy availability. Reduction in energy intake increases bioaccumulation and the amount of toxicant transferred from mother to each offspring. With high energy availability, the toxicant load of offspring decreases with birth order. Contrary to expectations, this ordering may be reversed with lower energy availability. Although demonstrated with parameters for the right whale, these relationships between energy intake and energetics and pharmacokinetics of organisms are likely to be much more general. Results specific to right whales include energy assimilation estimates for the North Atlantic and southern right whale, influences of history of energy availability on reproduction, and a relationship between ages at first parturition and calving intervals. Our model provides a platform for further analyses of both individual and population responses of marine mammals to pollution, and to changes in energy availability, including those likely to arise through climate change.

Dispersal delays, predator-prey stability, and the paradox of enrichment

It takes time for individuals to move from place to place. This travel time can be incorporated into metapopulation models via a delay in the interpatch migration term. Such a term has been shown to stabilize the positive equilibrium of the classical Lotka-Volterra predator-prey system with one species (either the predator or the prey) dispersing. We study a more realistic, Rosenzweig-MacArthur, model that includes a carrying capacity for the prey, and saturating functional response for the predator. We show that dispersal delays can stabilize the predator-prey equilibrium point despite the presence of a Type II functional response that is known to be destabilizing. We also show that dispersal delays reduce the amplitude of oscillations when the equilibrium is unstable, and therefore may help resolve the paradox of enrichment.

When can herbivores slow or reverse the spread of an invading plant? A test case from Mount St. Helens

Here we study the spatial dynamics of a coinvading consumer-resource pair. We present a theoretical treatment with extensive empirical data from a long-studied field system in which native herbivorous insects attack a population of lupine plants recolonizing a primary successional landscape created by the 1980 volcanic eruption of Mount St. Helens. Using detailed data on the life history and interaction strengths of the lupine and one of its herbivores, we develop a system of integrodifference equations to study plant-herbivore invasion dynamics. Our analyses yield several new insights into the spatial dynamics of coinvasions. In particular, we demonstrate that aspects of plant population growth and the intensity of herbivory under low-density conditions can determine whether the plant population spreads across a landscape or is prevented from doing so by the herbivore. In addition, we characterize the existence of threshold levels of spatial extent and/or temporal advantage for the plant that together define critical values of "invasion momentum," beyond which herbivores are unable to reverse a plant invasion. We conclude by discussing the implications of our findings for successional dynamics and the use of biological control agents to limit the spread of pest species.

Projecting rates of spread for invasive species

All else being equal, the faster an invading species spreads, the more dangerous its invasion. The projection of spread rate therefore ought to be a central part of the determination of invasion risk. Originally formulated in the 1970s to describe the spatial spread of advantageous alleles, integrodifference equation (IDE) models have since been co-opted by population biologists to describe the spread of populations. More recently, they have been modified to include population structure and environmental variability. We review how IDE models are formulated, how they are parameterized, and how they can be analyzed to project spread rates and the sensitivity of those rates to changes in model parameters. For illustrative purposes, we apply these models to Cytisus scoparius, a large shrub in the legume family that is considered a noxious invasive species in eastern and western North America, Chile, Australia, and New Zealand.

Habitat destruction, fragmentation, and disturbance promote invasion by habitat generalists in a multispecies metapopulation

Species invasions are extremely common and are vastly outpacing the ability of resource agencies to address each invasion, one species at a time. Management actions that target the whole landscape or ecosystem may provide more cost-effective protection against the establishment of invasive species than a species-by-species approach. To explore what ecosystem-level actions might effectively reduce invasions, we developed a multispecies, multihabitat metapopulation model. We assume that species that successfully establish themselves outside their native range tend to be habitat generalists and that a tradeoff exists between competitive ability and habitat breadth, such that habitat specialists are competitively superior to habitat generalists. In this model, habitat destruction, fragmentation, and short-term disturbances all favor invasion by habitat generalists, despite the inferior competitive abilities of generalist species. Our model results illustrate that providing relatively undisturbed habitat and preventing further habitat degradation and fragmentation can provide a highly cost-effective defense against invasive species.

Stabilizing dispersal delays in predator-prey metapopulation models

Time delays produced by dispersal are shown to stabilize Lotka-Volterra predator-prey models. The models are formulated as integrodifferential equations that describe local predator-prey dynamics and either intrapatch or interpatch dispersal. Dispersing individuals may (or may not) differ in the duration of their trips; these differences are captured via a distributed delay in the models. Our results include those of previous studies as special cases and show that the stabilizing effect continues to operate when the dispersal process is modeled more realistically.

Integrodifference equations, Allee effects, and invasions

Models that describe the spread of invading organisms often assume no Allee effect. In contrast, abundant observational data provide evidence for Allee effects. We study an invasion model based on an integrodifference equation with an Allee effect. We derive a general result for the sign of the speed of invasion. We then examine a special, linear-constant, Allee function and introduce a numerical scheme that allows us to estimate the speed of traveling wave solutions.

Transient dynamics and pattern formation: reactivity is necessary for Turing instabilities

The theory of spatial pattern formation via Turing bifurcations - wherein an equilibrium of a nonlinear system is asymptotically stable in the absence of dispersal but unstable in the presence of dispersal - plays an important role in biology, chemistry and physics. It is an asymptotic theory, concerned with the long-term behavior of perturbations. In contrast, the concept of reactivity describes the short-term transient behavior of perturbations to an asymptotically stable equilibrium. In this article we show that there is a connection between these two seemingly disparate concepts. In particular, we show that reactivity is necessary for Turing instability in multispecies systems of reaction-diffusion equations, integrodifference equations, coupled map lattices, and systems of ordinary differential equations.

Invasion speeds in fluctuating environments

Biological invasions are increasingly frequent and have dramatic ecological and economic consequences. A key to coping with invasive species is our ability to predict their rates of spread. Traditional models of biological invasions assume that the environment is temporally constant. We examine the consequences for invasion speed of periodic and stochastic fluctuations in population growth rates and in dispersal distributions.

Density-dependent vital rates and their population dynamic consequences

We explore a set of simple, nonlinear, two-stage models that allow us to compare the effects of density dependence on population dynamics among different kinds of life cycles. We characterize the behavior of these models in terms of their equilibria, bifurcations. and nonlinear dynamics, for a wide range of parameters. Our analyses lead to several generalizations about the effects of life history and density dependence on population dynamics. Among these are: (1) iteroparous life histories are more likely to be stable than semelparous life histories; (2) an increase in juvenile survivorship tends to be stabilizing; (3) density-dependent adult survival cannot control population growth when reproductive output is high: (4) density-dependent reproduction is more likely to cause chaotic dynamics than density dependence in other vital rates; and (5) changes in development rate have only small effects on bifurcation patterns.

A simple population model with qualitatively uncertain dynamics

Basins of attraction of nonlinear systems can be "riddled"; arbitrarily close to any point in a riddled basin are neighboring points which go to a different attractor. I present two chaotically forced single-species population models with riddled basins of attraction. As a result of this complex basin structure, the ultimate survival of these populations is effectively unpredictable. Riddled basins produce a level of unpredictability qualitatively greater than the familiar sensitive dependence on initial conditions within a single chaotic attractor, or the unpredictability caused by multiple attractors with fractal basin boundaries.Copyright 1997 Academic Press Limited Copyright 1997 Academic Press Limited

The subcritical collapse of predator populations in discrete-time predator-prey models

Many discrete-time predator-prey models possess three equilibria, corresponding to (1) extinction of both species, (2) extinction of the predator and survival of the prey at its carrying capacity, or (3) coexistence of both species. For a variety of such models, the equilibrium corresponding to coexistence may lose stability via a Hopf bifurcation, in which case trajectories approach an invariant circle. Alternatively, the equilibrium may undergo a subcritical flip bifurcation with a concomitant crash in the predator's population. We review a technique for distinguishing between subcritical and supercritical flip bifurcations and provide examples of predator-prey systems with a subcritical flip bifurcation.