Source author record

Sebastian J. Schreiber

Sebastian J. Schreiber appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Permanence via invasion graphs: Incorporating community assembly into Modern Coexistence Theory

To understand the mechanisms underlying species coexistence, ecologists often study invasion growth rates of theoretical and data-driven models. These growth rates correspond to average per-capita growth rates of one species with respect to an ergodic measure supporting other species. In the ecological literature, coexistence often is equated with the invasion growth rates being positive. Intuitively, positive invasion growth rates ensure that species recover from being rare. To provide a mathematically rigorous framework for this approach, we prove theorems that answer two questions: (i) When do the signs of the invasion growth rates determine coexistence? (ii) When signs are sufficient, which invasion growth rates need to be positive? We focus on deterministic models and equate coexistence with permanence, i.e., a global attractor bounded away from extinction. For models satisfying certain technical assumptions, we introduce invasion graphs where vertices correspond to proper subsets of species (communities) supporting an ergodic measure and directed edges correspond to potential transitions between communities due to invasions by missing species. These directed edges are determined by the signs of invasion growth rates. When the invasion graph is acyclic (i.e. there is no sequence of invasions starting and ending at the same community), we show that permanence is determined by the signs of the invasion growth rates. In this case, permanence is characterized by the invasibility of all -i communities, i.e., communities without species i where all other missing species having negative invasion growth rates. We show that dissipative Lotka-Volterra models satisfy our technical assumptions and computing their invasion graphs reduces to solving systems of linear equations. We provide additional applications of the results and discuss open problems.

preprint2020arXiv

Extinction and quasi-stationarity for discrete-time, endemic SIS and SIR models

Stochastic discrete-time SIS and SIR models of endemic diseases are introduced and analyzed. For the deterministic, mean-field model, the basic reproductive number $R_0$ determines their global dynamics. If $R_0\le 1$, then the frequency of infected individuals asymptotically converges to zero. If $R_0>1$, then the infectious class uniformly persists for all time; conditions for a globally stable, endemic equilibrium are given. In contrast, the infection goes extinct in finite time with probability one in the stochastic models for all $R_0$ values. To understand the length of the transient prior to extinction as well as the behavior of the transients, the quasi-stationary distributions and the associated mean time to extinction are analyzed using large deviation methods. When $R_0>1$, these mean times to extinction are shown to increase exponentially with the population size $N$. Moreover, as $N$ approaches $\infty$, the quasi-stationary distributions are supported by a compact set bounded away from extinction; sufficient conditions for convergence to a Dirac measure at the endemic equilibrium of the deterministic model are also given. In contrast, when $R_0<1$, the mean times to extinction are bounded above $1/(1-α)$ where $α<1$ is the geometric rate of decrease of the infection when rare; as $N$ approaches $\infty$, the quasi-stationary distributions converge to a Dirac measure at the disease-free equilibrium for the deterministic model. For several special cases, explicit formulas for approximating the quasi-stationary distribution and the associated mean extinction are given. These formulas illustrate how for arbitrarily small $R_0$ values, the mean time to extinction can be arbitrarily large, and how for arbitrarily large $R_0$ values, the mean time to extinction can be arbitrarily large.

preprint2020arXiv

Multiple attractors and long transients in spatially structured populations with an Allee effect

We present a discrete-time model of a spatially structured population and explore the effects of coupling when the local dynamics contain a strong Allee effect and overcompensation. While an isolated population can exhibit only bistability and essential extinction, a spatially structured population can exhibit numerous coexisting attractors. We identify mechanisms and parameter ranges that can protect the spatially structured population from essential extinction, whereas it is inevitable in the local system. In the case of weak coupling, a state where one subpopulation density lies above and the other one below the Allee threshold can prevent essential extinction. Strong coupling, on the other hand, enables both populations to persist above the Allee threshold when dynamics are (approximately) out-of-phase. In both cases, attractors have fractal basin boundaries. Outside of these parameter ranges, dispersal was not found to prevent essential extinction. We also demonstrate how spatial structure can lead to long transients of persistence before the population goes extinct.

preprint2014arXiv

Protected polymorphisms and evolutionary stability of patch-selection strategies in stochastic environments

We consider a population living in a patchy environment that varies stochastically in space and time. The population is composed of two morphs (that is, individuals of the same species with different genotypes). In terms of survival and reproductive success, the associated phenotypes differ only in their habitat selection strategies. We compute invasion rates corresponding to the rates at which the abundance of an initially rare morph increases in the presence of the other morph established at equilibrium. If both morphs have positive invasion rates when rare, then there is an equilibrium distribution such that the two morphs coexist; that is, there is a protected polymorphism for habitat selection. Alternatively, if one morph has a negative invasion rate when rare, then it is asymptotically displaced by the other morph under all initial conditions where both morphs are present. We refine the characterization of an evolutionary stable strategy for habitat selection from [Schreiber, 2012] in a mathematically rigorous manner. We provide a necessary and sufficient condition for the existence of an ESS that uses all patches and determine when using a single patch is an ESS. We also provide an explicit formula for the ESS when there are two habitat types. We show that adding environmental stochasticity results in an ESS that, when compared to the ESS for the corresponding model without stochasticity, spends less time in patches with larger carrying capacities and possibly makes use of sink patches, thereby practicing a spatial form of bet hedging.

preprint2014arXiv

Quasi-stationary distributions for randomly perturbed dynamical systems

We analyze quasi-stationary distributions $\{μ^{\varepsilon}\}_{\varepsilon>0}$ of a family of Markov chains $\{X^{\varepsilon}\}_{\varepsilon>0}$ that are random perturbations of a bounded, continuous map $F:M\to M$, where $M$ is a closed subset of $\mathbb{R}^k$. Consistent with many models in biology, these Markov chains have a closed absorbing set $M_0\subset M$ such that $F(M_0)=M_0$ and $F(M\setminus M_0)=M\setminus M_0$. Under some large deviations assumptions on the random perturbations, we show that, if there exists a positive attractor for $F$ (i.e., an attractor for $F$ in $M\setminus M_0$), then the weak* limit points of $μ_{\varepsilon}$ are supported by the positive attractors of $F$. To illustrate the broad applicability of these results, we apply them to nonlinear branching process models of metapopulations, competing species, host-parasitoid interactions and evolutionary games.

preprint2013arXiv

Persistence in fluctuating environments for interacting structured populations

Individuals within any species exhibit differences in size, developmental state, or spatial location. These differences coupled with environmental fluctuations in demographic rates can have subtle effects on population persistence and species coexistence. To understand these effects, we provide a general theory for coexistence of structured, interacting species living in a stochastic environment. The theory is applicable to nonlinear, multi species matrix models with stochastically varying parameters. The theory relies on long-term growth rates of species corresponding to the dominant Lyapunov exponents of random matrix products. Our coexistence criterion requires that a convex combination of these long-term growth rates is positive with probability one whenever one or more species are at low density. When this condition holds, the community is stochastically persistent: the fraction of time that a species density goes below $δ>0$ approaches zero as $δ$ approaches zero. Applications to predator-prey interactions in an autocorrelated environment, a stochastic LPA model, and spatial lottery models are provided. These applications demonstrate that positive autocorrelations in temporal fluctuations can disrupt predator-prey coexistence, fluctuations in log-fecundity can facilitate persistence in structured populations, and long-lived, relatively sedentary competing populations are likely to coexist in spatially and temporally heterogenous environments.

preprint2012arXiv

Spatial heterogeneity promotes coexistence of rock-paper-scissor metacommunities

The rock-paper-scissor game -- which is characterized by three strategies R,P,S, satisfying the non-transitive relations S excludes P, P excludes R, and R excludes S -- serves as a simple prototype for studying more complex non-transitive systems. For well-mixed systems where interactions result in fitness reductions of the losers exceeding fitness gains of the winners, classical theory predicts that two strategies go extinct. The effects of spatial heterogeneity and dispersal rates on this outcome are analyzed using a general framework for evolutionary games in patchy landscapes. The analysis reveals that coexistence is determined by the rates at which dominant strategies invade a landscape occupied by the subordinate strategy (e.g. rock invades a landscape occupied by scissors) and the rates at which subordinate strategies get excluded in a landscape occupied by the dominant strategy (e.g. scissor gets excluded in a landscape occupied by rock). These invasion and exclusion rates correspond to eigenvalues of the linearized dynamics near single strategy equilibria. Coexistence occurs when the product of the invasion rates exceeds the product of the exclusion rates. Provided there is sufficient spatial variation in payoffs, the analysis identifies a critical dispersal rate $d^*$ required for regional persistence. For dispersal rates below $d^*$, the product of the invasion rates exceed the product of the exclusion rates and the rock-paper-scissor metacommunities persist regionally despite being extinction prone locally. For dispersal rates above $d^*$, the product of the exclusion rates exceed the product of the invasion rates and the strategies are extinction prone. These results highlight the delicate interplay between spatial heterogeneity and dispersal in mediating long-term outcomes for evolutionary games.

preprint2012arXiv

Stochastic population growth in spatially heterogeneous environments

Classical ecological theory predicts that environmental stochasticity increases extinction risk by reducing the average per-capita growth rate of populations. To understand the interactive effects of environmental stochasticity, spatial heterogeneity, and dispersal on population growth, we study the following model for population abundances in $n$ patches: the conditional law of $X_{t+dt}$ given $X_t=x$ is such that when $dt$ is small the conditional mean of $X_{t+dt}^i-X_t^i$ is approximately $[x^iμ_i+\sum_j(x^j D_{ji}-x^i D_{ij})]dt$, where $X_t^i$ and $μ_i$ are the abundance and per capita growth rate in the $i$-th patch respectivly, and $D_{ij}$ is the dispersal rate from the $i$-th to the $j$-th patch, and the conditional covariance of $X_{t+dt}^i-X_t^i$ and $X_{t+dt}^j-X_t^j$ is approximately $x^i x^j σ_{ij}dt$. We show for such a spatially extended population that if $S_t=(X_t^1+...+X_t^n)$ is the total population abundance, then $Y_t=X_t/S_t$, the vector of patch proportions, converges in law to a random vector $Y_\infty$ as $t\to\infty$, and the stochastic growth rate $\lim_{t\to\infty}t^{-1}\log S_t$ equals the space-time average per-capita growth rate $\sum_iμ_i\E[Y_\infty^i]$ experienced by the population minus half of the space-time average temporal variation $\E[\sum_{i,j}σ_{ij}Y_\infty^i Y_\infty^j]$ experienced by the population. We derive analytic results for the law of $Y_\infty$, find which choice of the dispersal mechanism $D$ produces an optimal stochastic growth rate for a freely dispersing population, and investigate the effect on the stochastic growth rate of constraints on dispersal rates. Our results provide fundamental insights into "ideal free" movement in the face of uncertainty, the persistence of coupled sink populations, the evolution of dispersal rates, and the single large or several small (SLOSS) debate in conservation biology.

preprint2012arXiv

Temporally variable dispersal and demography can accelerate the spread of invading species

We analyze how temporal variability in local demography and dispersal combine to affect the rate of spread of an invading species. Our model combines state-structured local demography (specified by an integral or matrix projection model) with general dispersal distributions that may depend on the state of the individual or its parent, and it allows very general patterns of stationary temporal variation in both local demography and in the frequency and distribution of dispersal distances. We show that expressions for the asymptotic spread rate and its sensitivity to parameters, that have been derived previously for less general models, continue to hold. Using these results, we show that random temporal variability in dispersal can accelerate population spread. Demographic variability can further accelerate spread if it is positively correlated with dispersal variability, for example if high-fecundity years are also years in which juveniles tend to settle further away from their parents. A simple model for the growth and spread of patches of an invasive plant (perennial pepperweed, Lepidium latifolium) illustrates these effects and shows that they can have substantial impacts on the predicted speed of an invasion wave. Temporal variability in dispersal has gotten very little attention in both the theoretical and empirical literatures on invasive species spread. Our results suggest that this needs to change.

preprint2011arXiv

Persistence for stochastic difference equations: A mini-review

Understanding under what conditions populations, whether they be plants, animals, or viral particles, persist is an issue of theoretical and practical importance in population biology. Both biotic interactions and environmental fluctuations are key factors that can facilitate or disrupt persistence. One approach to examining the interplay between these deterministic and stochastic forces is the construction and analysis of stochastic difference equations $X_{t+1}=F(X_t,ξ_{t+1})$ where $X_t \in \R^k$ represents the state of the populations and $ξ_1,ξ_2,...$ is a sequence of random variables representing environmental stochasticity. In the analysis of these stochastic models, many theoretical population biologists are interested in whether the models are bounded and persistent. Here, boundedness asserts that asymptotically $X_t$ tends to remain in compact sets. In contrast, persistence requires that $X_t$ tends to be "repelled" by some "extinction set" $S_0\subset \R^k$. Here, results on both of these proprieties are reviewed for single species, multiple species, and structured population models. The results are illustrated with applications to stochastic versions of the Hassell and Ricker single species models, Ricker, Beverton-Holt, lottery models of competition, and lottery models of rock-paper-scissor games. A variety of conjectures and suggestions for future research are presented.

preprint2010arXiv

Evolution of unconditional dispersal in periodic environments

Organisms modulate their fitness in heterogeneous environments by dispersing. Prior work shows that there is selection against "unconditional" dispersal in spatially heterogeneous environments. "Unconditional" means individuals disperse at a rate independent of their location. We prove that if within-patch fitness varies spatially and between two values temporally, then there is selection for unconditional dispersal: any evolutionarily stable strategy (ESS) or evolutionarily stable coalition (ESC) includes a dispersive phenotype. Moreover, at this ESS or ESC, there is at least one sink patch (i.e. geometric mean of fitness less than one) and no sources patches (i.e. geometric mean of fitness greater than one). These results coupled with simulations suggest that spatial-temporal heterogeneity due to abiotic forcing result in either an ESS with a dispersive phenotype or an ESC with sedentary and dispersive phenotypes. In contrast, spatial-temporal heterogeneity due to biotic interactions can select for higher dispersal rates that ultimately spatially synchronize population dynamics.

preprint2010arXiv

Invasion speeds for structured populations in fluctuating environments

We live in a time where climate models predict future increases in environmental variability and biological invasions are becoming increasingly frequent. A key to developing effective responses to biological invasions in increasingly variable environments will be estimates of their rates of spatial spread and the associated uncertainty of these estimates. Using stochastic, stage-structured, integro-difference equation models, we show analytically that invasion speeds are asymptotically normally distributed with a variance that decreases in time. We apply our methods to a simple juvenile-adult model with stochastic variation in reproduction and an illustrative example with published data for the perennial herb, \emph{Calathea ovandensis}. These examples buttressed by additional analysis reveal that increased variability in vital rates simultaneously slow down invasions yet generate greater uncertainty about rates of spatial spread. Moreover, while temporal autocorrelations in vital rates inflate variability in invasion speeds, the effect of these autocorrelations on the average invasion speed can be positive or negative depending on life history traits and how well vital rates ``remember'' the past.

preprint2010arXiv

Persistence in fluctuating environments

Understanding under what conditions interacting populations, whether they be plants, animals, or viral particles, coexist is a question of theoretical and practical importance in population biology. Both biotic interactions and environmental fluctuations are key factors that can facilitate or disrupt coexistence. To better understand this interplay between these deterministic and stochastic forces, we develop a mathematical theory extending the nonlinear theory of permanence for deterministic systems to stochastic difference and differential equations. Our condition for coexistence requires that there is a fixed set of weights associated with the interacting populations and this weighted combination of populations' invasion rates is positive for any (ergodic) stationary distribution associated with a subcollection of populations. Here, an invasion rate corresponds to an average per-capita growth rate along a stationary distribution. When this condition holds and there is sufficient noise in the system, we show that the populations approach a unique positive stationary distribution. Moreover, we show that our coexistence criterion is robust to small perturbations of the model functions. Using this theory, we illustrate that (i) environmental noise enhances or inhibits coexistence in communities with rock-paper-scissor dynamics depending on correlations between interspecific demographic rates, (ii) stochastic variation in mortality rates has no effect on the coexistence criteria for discrete-time Lotka-Volterra communities, and (iii) random forcing can promote genetic diversity in the presence of exploitative interactions.

preprint2010arXiv

Robust permanence for interacting structured populations

The dynamics of interacting structured populations can be modeled by $\frac{dx_i}{dt}= A_i (x)x_i$ where $x_i\in \R^{n_i}$, $x=(x_1,\dots,x_k)$, and $A_i(x)$ are matrices with non-negative off-diagonal entries. These models are permanent if there exists a positive global attractor and are robustly permanent if they remain permanent following perturbations of $A_i(x)$. Necessary and sufficient conditions for robust permanence are derived using dominant Lyapunov exponents $λ_i(μ)$ of the $A_i(x)$ with respect to invariant measures $μ$. The necessary condition requires $\max_i λ_i(μ)>0$ for all ergodic measures with support in the boundary of the non-negative cone. The sufficient condition requires that the boundary admits a Morse decomposition such that $\max_i λ_i(μ)>0$ for all invariant measures $μ$ supported by a component of the Morse decomposition. When the Morse components are Axiom A, uniquely ergodic, or support all but one population, the necessary and sufficient conditions are equivalent. Applications to spatial ecology, epidemiology, and gene networks are given.