Search

Spain is the largest producer of avocado and mango fruits in Europe. The majority of production is concentrated in the Axarqu\'ia region in the south, where subtropical fruit plantations and associated water demands have steadily increased over the last two decades. Between 2019-2024, the region underwent an extreme water crisis. Reservoir reserves became nearly depleted and groundwater levels dropped to sea level in several locations, where seawater intrusion is likely, causing large socioeconomic impacts including short-term harvest losses and a long-term loss in economic centrality. We examine the causal pathway that led to this crisis using a mixed-methods approach, combining data from key informant interviews, an exhaustive review of legal documents, and quantitative analysis of time series and spatially explicit data. In particular, we analyze dam water use for irrigation and urban use, meteorological data, reservoir and groundwater levels, and irrigation land cover maps. Our findings show that an unusual meteorological drought was the immediate cause for the decline in reservoir and groundwater reserves (hydrological drought), but the underlying cause was a chronic and structural long-term imbalance between water demand and resources resulting from several structural governance shortcomings: large uncertainties in water resource availability and use hampering effective planning, lack of enforcement of individual water quotas, and the absence of regulatory mechanisms to flexibly impose resource use restrictions at both micro and macro levels based on the overall resources of the system. We propose concrete policy interventions aimed at sustainably enhancing the resilience of the system that can be useful to efficiently manage water shortages in other regions with similar problems.

Risk-driven behavior provides a feedback mechanism through which individuals both shape and are collectively affected by an epidemic. We introduce a general and flexible compartmental model to study the effect of heterogeneity in the population with regards to risk tolerance. The interplay between behavior and epidemiology leads to a rich set of possible epidemic dynamics. Depending on the behavioral composition of the population, we find that increasing heterogeneity in risk tolerance can either increase or decrease the epidemic size. We find that multiple waves of infection can arise due to the interplay between transmission and behavior, even without the replenishment of susceptibles. We find that increasing protective mechanisms such as the effectiveness of interventions, the number of risk-averse people in the population, and the duration of intervention usage reduces the epidemic overshoot. When the protection is pushed past a critical threshold, the epidemic dynamics enter an underdamped regime where the epidemic size exactly equals the herd immunity threshold. Lastly, we can find regimes where epidemic size does not monotonically decrease with a population that becomes increasingly risk-averse., Comment: 12 pages, 7 figures (main text)

Based on the success of a well-known method for solving higher order linear differential equations, a study of two of the most important mathematical features of that method, viz. the null spaces and commutativity of the product of unbounded linear operators on a Hilbert space is carried out. A principle is proved describing solutions for the product of such operators in terms of the solutions for each of the factors when the null spaces of those factors satisfy a certain geometric relation to one another. Another geometric principle equating commutativity of a closed densely defined operator and a projection to stability of the range of the projection under the closed operator is proved.

We present and analyze an actively controlled Susceptible-Infected-Susceptible (actSIS) model of interconnected populations to study how risk aversion strategies, such as social distancing, affect network epidemics. A population using a risk aversion strategy reduces its contact rate with other populations when it perceives an increase in infection risk. The network actSIS model relies on two distinct networks. One is a physical contact network that defines which populations come into contact with which other populations and thus how infection spreads. The other is a communication network, such as an online social network, that defines which populations observe the infection level of which other populations and thus how information spreads. We prove that the model, with these two networks and populations using risk aversion strategies, exhibits a transcritical bifurcation in which an endemic equilibrium emerges. For regular graphs, we prove that the endemic infection level is uniform across populations and reduced by the risk aversion strategy, relative to the network SIS endemic level. We show that when communication is sufficiently sparse, this initially stable equilibrium loses stability in a secondary bifurcation. Simulations show that a new stable solution emerges with nonuniform infection levels.

In the age of information abundance, attention is a coveted resource. Social media platforms vigorously compete for users' engagement, influencing the evolution of their opinions on a variety of topics. With recommendation algorithms often accused of creating "filter bubbles", where like-minded individuals interact predominantly with one another, it's crucial to understand the consequences of this unregulated attention market. To address this, we present a model of opinion dynamics on a multiplex network. Each layer of the network represents a distinct social media platform, each with its unique characteristics. Users, as nodes in this network, share their opinions across platforms and decide how much time to allocate in each platform depending on its perceived quality. Our model reveals two key findings. i) When examining two platforms - one with a neutral recommendation algorithm and another with a homophily-based algorithm - we uncover that even if users spend the majority of their time on the neutral platform, opinion polarization can persist. ii) By allowing users to dynamically allocate their social energy across platforms in accordance to their homophilic preferences, a further segregation of individuals emerges. While network fragmentation is usually associated with "echo chambers", the emergent multi-platform segregation leads to an increase in users' satisfaction without the undesired increase in polarization. These results underscore the significance of acknowledging how individuals gather information from a multitude of sources. Furthermore, they emphasize that policy interventions on a single social media platform may yield limited impact., Comment: 13 pages, 7 figures

Complex adaptive systems (CASs), from ecosystems to economies, are open systems and inherently dependent on external conditions. While a system can transition from one state to another based on the magnitude of change in external conditions, the rate of change -- irrespective of magnitude -- may also lead to system state changes due to a phenomenon known as a rate-induced transition (RIT). This study presents a novel framework that captures RITs in CASs through a local model and a network extension where each node contributes to the structural adaptability of others. Our findings reveal how RITs occur at a critical environmental change rate, with lower-degree nodes tipping first due to fewer connections and reduced adaptive capacity. High-degree nodes tip later as their adaptability sources (lower-degree nodes) collapse. This pattern persists across various network structures. Our study calls for an extended perspective when managing CASs, emphasizing the need to focus not only on thresholds of external conditions but also the rate at which those conditions change, particularly in the context of the collapse of surrounding systems that contribute to the focal system's resilience. Our analytical method opens a path to designing management policies that mitigate RIT impacts and enhance resilience in ecological, social, and socioecological systems. These policies could include controlling environmental change rates, fostering system adaptability, implementing adaptive management strategies, and building capacity and knowledge exchange. Our study contributes to the understanding of RIT dynamics and informs effective management strategies for complex adaptive systems in the face of rapid environmental change., Comment: 25 pages, 4 figures, 1 box, supplementary information

We analyze a spatially extended version of a well-known model of forest-savanna dynamics, which presents as a system of nonlinear partial integro-differential equations, and study necessary conditions for pattern-forming bifurcations. Analytically, we show that homogeneous solutions dominate the dynamics of the standard forest-savanna model, regardless of the length scales of the various spatial processes considered. However, several different pattern-forming scenarios are possible upon including spatial resource limitation, such as competition for water, soil nutrients, or herbivory effects. Using numerical simulations and continuation, we study the nature of the resulting patterns as a function of system parameters and length scales, uncovering subcritical pattern-forming bifurcations and observing significant regions of multistability for realistic parameter regimes. Finally, we discuss our results in the context of extant savanna-forest modeling efforts and highlight ongoing challenges in building a unifying mathematical model for savannas across different rainfall levels., Comment: 27 pages, 7 figures

Spatial systems with heterogeneities are ubiquitous in nature, from precipitation, temperature and soil gradients controlling vegetation growth to morphogen gradients controlling gene expression in embryos. Such systems, generally described by nonlinear dynamical systems, often display complex parameter dependence and exhibit bifurcations. The dynamics of heterogeneous spatially extended systems passing through bifurcations are still relatively poorly understood, yet recent theoretical studies and experimental data highlight the resulting complex behaviors and their relevance to real-world applications. We explore the consequences of spatial heterogeneities passing through bifurcations via two examples strongly motivated by applications. These model systems illustrate that studying heterogeneity-induced behaviors in spatial systems is crucial for a better understanding of ecological transitions and functional organization in brain development., Comment: 28 pages

We derive an exact upper bound on the epidemic overshoot for the Kermack-McKendrick SIR model. This maximal overshoot value of 0.2984... occurs at $R_0^*$ = 2.151... . In considering the utility of the notion of overshoot, a rudimentary analysis of data from the first wave of the COVID-19 pandemic in Manaus, Brazil highlights the public health hazard posed by overshoot for epidemics with $R_0$ near 2. Using the general analysis framework presented within, we then consider more complex SIR models that incorporate vaccination., Comment: Main text: 7 pages + 2 figures. Supplement: 7 pages + 6 figures + code

A key scientific challenge during the outbreak of novel infectious diseases is to predict how the course of the epidemic changes under different countermeasures that limit interaction in the population. Most epidemiological models do not consider the role of mutations and heterogeneity in the type of contact events. However, pathogens have the capacity to mutate in response to changing environments, especially caused by the increase in population immunity to existing strains and the emergence of new pathogen strains poses a continued threat to public health. Further, in light of differing transmission risks in different congregate settings (e.g., schools and offices), different mitigation strategies may need to be adopted to control the spread of infection. We analyze a multi-layer multi-strain model by simultaneously accounting for i) pathways for mutations in the pathogen leading to the emergence of new pathogen strains, and ii) differing transmission risks in different congregate settings, modeled as network-layers. Assuming complete cross-immunity among strains, namely, recovery from any infection prevents infection with any other (an assumption that will need to be relaxed to deal with COVID-19 or influenza), we derive the key epidemiological parameters for the proposed multi-layer multi-strain framework. We demonstrate that reductions to existing network-based models that discount heterogeneity in either the strain or the network layers can lead to incorrect predictions for the course of the outbreak. In addition, our results highlight that the impact of imposing/lifting mitigation measures concerning different contact network layers (e.g., school closures or work-from-home policies) should be evaluated in connection with their effect on the likelihood of the emergence of new pathogen strains.

Biological and social systems are structured at multiple scales, and the incentives of individuals who interact in a group may diverge from the collective incentive of the group as a whole. Mechanisms to resolve this tension are responsible for profound transitions in evolutionary history, including the origin of cellular life, multi-cellular life, and even societies. Here we synthesize a growing literature that extends evolutionary game theory to describe multilevel evolutionary dynamics, using nested birth-death processes and partial differential equations to model natural selection acting on competition within and among groups of individuals. We apply this theory to analyze how mechanisms known to promote cooperation within a single group -- including assortment, reciprocity, and population structure -- alter evolutionary outcomes in the presence of competition among groups. We find that population structures most conducive to cooperation in multi-scale systems may differ from those most conducive within a single group. Likewise, for competitive interactions with a continuous range of strategies we find that among-group selection may fail to produce socially optimal outcomes, but it can nonetheless produce second-best solutions that balance individual incentives to defect with the collective incentives for cooperation. We conclude by describing the broad applicability of multi-scale evolutionary models to problems ranging from the production of diffusible metabolites in microbes to the management of common-pool resources in human societies., Comment: 48 pages, 8 figures

Anthropogenic emissions of CO2 must soon approach net-zero to stabilize the global mean temperature. Although several international agreements have advocated for coordinated climate actions, their implementation has remained below expectations. One of the main challenges of international cooperation is the different degrees of socio-political acceptance of decarbonization. Here we interrogate a minimalistic model of the coupled natural-human system representing the impact of such socio-political acceptance on clean energy investments and the path to net-zero emissions. We show that incentives and carbon pricing are essential tools to achieve net-zero before critical CO2 concentrations are reached, and deep international coordination is necessary for a rapid and effective transition. Although a perfect coordination scenario remains unlikely, as investments in clean energy are ultimately limited by myopic economic strategies and a policy system that promotes free-riding, more realistic decentralized cooperation with partial efforts from each actor could still lead to significant emissions cuts.

Levels of sociality in nature vary widely. Some species are solitary; others live in family groups; some form complex multi-family societies. Increased levels of social interaction can allow for the spread of useful innovations and beneficial information, but can also facilitate the spread of harmful contagions, such as infectious diseases. It is natural to assume that these contagion processes shape the evolution of complex social systems, but an explicit account of the dynamics of sociality under selection pressure imposed by contagion remains elusive. We consider a model for the evolution of sociality strategies in the presence of both a beneficial and costly contagion. We study the dynamics of this model at three timescales: using a susceptible-infectious-susceptible (SIS) model to describe contagion spread for given sociality strategies, a replicator equation to study the changing fractions of two different levels of sociality, and an adaptive dynamics approach to study the long-time evolution of the population level of sociality. For a wide range of assumptions about the benefits and costs of infection, we identify a social dilemma: the evolutionarily-stable sociality strategy (ESS) is distinct from the collective optimum -- the level of sociality that would be best for all individuals. In particular, the ESS level of social interaction is greater (respectively less) than the social optimum when the good contagion spreads more (respectively less) readily than the bad contagion. Our results shed light on how contagion shapes the evolution of social interaction, but reveals that evolution may not necessarily lead populations to social structures that are good for any or all., Comment: 42 pages, 9 figures, equal contribution from DBC and DHM, v2: Revised sections on general utility and linear utility functions and added new section in appendix on Constant Elasticity of Substitution (CES) utility

Critical transitions, or large changes in the state of a system after a small change in the system's external conditions or parameters, commonly occur in a wide variety of disciplines, from the biological and social sciences to physics. Statistical physics first confronted the problem of emergent phenomena such as critical transitions in the 1800s and 1900s, culminating in the theory of phase transitions. However, although phase transitions show a strong resemblance to critical transitions, the theoretical connections between the two sets of phenomena are tenuous at best, and it would be advantageous to make them more concrete in order to take advantage of the theoretical methods developed by physicists to study phase transitions. Here we attempt to explicitly connect the theory of critical transitions to phase transitions in physics. We initially find something paradoxical, that many critical transitions closely resemble first-order phase transitions, but that many of the early warning indicators developed to anticipate critical transitions, such as critical slowing down or increasing spatial correlations, occur instead in second-order phase transitions. We attempt to reconcile these disparities by making the connection with other phenomena associated with first-order phase transitions, such as spinodal instabilities and metastable states., Comment: 16 pages, 5 figures, To appear as book chapter in Global Systemic Risk

Consider a cooperation game on a spatial network of habitat patches, where players can relocate between patches if they judge the local conditions to be unfavorable. In time, the relocation events may lead to a homogeneous state where all patches harbor the same relative densities of cooperators and defectors or they may lead to self-organized patterns, where some patches become safe havens that maintain an elevated cooperator density. Here we analyze the transition between these states mathematically. We show that safe havens form once a certain threshold in connectivity is crossed. This threshold can be analytically linked to the structure of the patch network and specifically to certain network motifs. Surprisingly, a forgiving defector avoidance strategy may be most favorable for cooperators. Our results demonstrate that the analysis of cooperation games in ecological metacommunity models is mathematically tractable and has the potential to link topics such as macroecological patterns, behavioral evolution, and network topology., Comment: 14 pages, 4 figures

Masks have remained an important mitigation strategy in the fight against COVID-19 due to their ability to prevent the transmission of respiratory droplets between individuals. In this work, we provide a comprehensive quantitative analysis of the impact of mask-wearing. To this end, we propose a novel agent-based model of viral spread on networks where agents may either wear no mask or wear one of several types of masks with different properties (e.g., cloth or surgical). We derive analytical expressions for three key epidemiological quantities: the probability of emergence, the epidemic threshold, and the expected epidemic size. In particular, we show how the aforementioned quantities depend on the structure of the contact network, viral transmission dynamics, and the distribution of the different types of masks within the population. Through extensive simulations, we then investigate the impact of different allocations of masks within the population and trade-offs between the outward efficiency and inward efficiency of the masks. Interestingly, we find that masks with high outward efficiency and low inward efficiency are most useful for controlling the spread in the early stages of an epidemic, while masks with high inward efficiency but low outward efficiency are most useful in reducing the size of an already large spread. Lastly, we study whether degree-based mask allocation is more effective in reducing the probability of epidemic as well as epidemic size compared to random allocation. The result echoes the previous findings that mitigation strategies should differ based on the stage of the spreading process, focusing on source control before the epidemic emerges and on self-protection after the emergence., Comment: Accepted at Physical Review E

The evolution of complex cellular life involved two major transitions: the encapsulation of self-replicating genetic entities into cellular units and the aggregation of individual genes into a collectively replicating genome. In this paper, we formulate a minimal model of the evolution of proto-chromosomes within protocells. We model a simple protocell composed of two types of genes: a "fast gene" with an advantage for gene-level self-replication and a "slow gene" that replicates more slowly at the gene level, but which confers an advantage for protocell-level reproduction. Protocell-level replication capacity depends on cellular composition of fast and slow genes. We use a partial differential equation to describe how the composition of genes within protocells evolves over time under within-cell and between-cell competition. We find that the gene-level advantage of fast replicators casts a long shadow on the multilevel dynamics of protocell evolution: no level of between-protocell competition can produce coexistence of the fast and slow replicators when the two genes are equally needed for protocell-level reproduction. By introducing a "dimer replicator" consisting of a linked pair of the slow and fast genes, we show analytically that coexistence between the two genes can be promoted in pairwise multilevel competition between fast and dimer replicators, and provide numerical evidence for coexistence in trimorphic competition between fast, slow, and dimer replicators. Our results suggest that dimerization, or the formation of a simple chromosome-like dimer replicator, can help to overcome the shadow of lower-level selection and work in concert with deterministic multilevel selection to allow for the coexistence of two genes that are complementary at the protocell-level but compete at the level of individual gene-level replication., Comment: 75 pages, 22 figures

We explore the commonalities between methods for assuring the security of computer systems (cybersecurity) and the mechanisms that have evolved through natural selection to protect vertebrates against pathogens, and how insights derived from studying the evolution of natural defenses can inform the design of more effective cybersecurity systems. More generally, security challenges are crucial for the maintenance of a wide range of complex adaptive systems, including financial systems, and again lessons learned from the study of the evolution of natural defenses can provide guidance for the protection of such systems., Comment: 20 pages

We develop a landscape-flux framework to investigate observed frequency distributions of vegetation and the stability of these ecological systems under fluctuations. The frequency distributions can characterize the population-potential landscape related to the stability of ecological states. We illustrate the practical utility of this approach by analyzing a forest-savanna model. Savanna, and Forest states coexist under certain conditions, consistent with past theoretical work and empirical observations. However, a new Grassland state, unseen in the corresponding deterministic model, emerges as an alternative quasi-stable state under fluctuations, providing a novel theoretical basis for the appearance of widespread grasslands in some empirical analyses. The ecological dynamics are determined by both the population-potential landscape gradient and the steady-state probability flux. The flux quantifies the net input/output to the ecological system and therefore the degree of nonequilibriumness. Landscape and flux together determine the transitions between stable states characterized by dominant paths and switching rates. The intrinsic potential landscape admits a Lyapunov function, which provides a quantitative measure of global stability. We find that the average flux, entropy production rate, and free energy have significant changes near bifurcations under both finite and zero fluctuation. These may provide both dynamical and thermodynamic origins of the bifurcations. We identified the variances in observed frequency time traces, fluctuations and time irreversibility as kinematic measures for bifurcations. This new framework opens the way to characterize ecological systems globally, to uncover how they change among states, and to quantify the emergence of new quasi-stable states under stochastic fluctuations.

The R package "sensobol" provides several functions to conduct variance-based uncertainty and sensitivity analysis, from the estimation of sensitivity indices to the visual representation of the results. It implements several state-of-the-art first and total-order estimators and allows the computation of up to third-order effects, as well as of the approximation error, in a swift and user-friendly way. Its flexibility makes it also appropriate for models with either a scalar or a multivariate output. We illustrate its functionality by conducting a variance-based sensitivity analysis of three classic models: the Sobol' (1998) G function, the logistic population growth model of Verhulst (1845), and the spruce budworm and forest model of Ludwig, Jones and Holling (1976)., Comment: The paper has been accepted for publication in the Journal of Statistical Software (expected for 2022)

A central need in the field of astrobiology is generalized perspectives on life that make it possible to differentiate abiotic and biotic chemical systems. A key component of many past and future astrobiological measurements is the elemental ratio of various samples. Classic work on Earth's oceans has shown that life displays a striking regularity in the ratio of elements as originally characterized by Redfield. The body of work since the original observations has connected this ratio with basic ecological dynamics and cell physiology, while also documenting the range of elemental ratios found in a variety of environments. Several key questions remain in considering how to best apply this knowledge to astrobiological contexts: How can the observed variation of the elemental ratios be more formally systematized using basic biological physiology and ecological or environmental dynamics? How can these elemental ratios be generalized beyond the life that we have observed on our own planet? Here we expand recently developed generalized physiological models to create a simple framework for predicting the variation of elemental ratios found in various environments. We then discuss further generalizing the physiology for astrobiological applications. Much of our theoretical treatment is designed for in situ measurements applicable to future planetary missions. We imagine scenarios where three measurements can be made - particle/cell sizes, particle/cell stoichiometry, and fluid or environmental stoichiometry - and develop our theory in connection with these often deployed measurements., Comment: 18 pages, 5 figures

Recent attempts at cooperating on climate change mitigation highlight the limited efficacy of large-scale agreements, when commitment to mitigation is costly and initially rare. Bottom-up approaches using region-specific mitigation agreements promise greater success, at the cost of slowing global adoption. Here, we show that a well-timed switch from regional to global negotiations dramatically accelerates climate mitigation compared to using only local, only global, or both agreement types simultaneously. This highlights the scale-specific roles of mitigation incentives: local incentives capitalize on regional differences (e.g., where recent disasters incentivize mitigation) by committing early-adopting regions, after which global agreements draw in late-adopting regions. We conclude that global agreements are key to overcoming the expenses of mitigation and economic rivalry among regions but should be attempted once regional agreements are common. Gradually up-scaling efforts could likewise accelerate mitigation at smaller scales, for instance when costly ecosystem restoration initially faces limited public and legislative support.

An actively controlled Susceptible-Infected-Susceptible (actSIS) contagion model is presented for studying epidemic dynamics with continuous-time feedback control of infection rates. Our work is inspired by the observation that epidemics can be controlled through decentralized disease-control strategies such as quarantining, sheltering in place, social distancing, etc., where individuals actively modify their contact rates with others in response to observations of infection levels in the population. Accounting for a time lag in observations and categorizing individuals into distinct sub-populations based on their risk profiles, we show that the actSIS model manifests qualitatively different features as compared with the SIS model. In a homogeneous population of risk-averters, the endemic equilibrium is always reduced, although the transient infection level can exhibit overshoot or undershoot. In a homogeneous population of risk-tolerating individuals, the system exhibits bistability, which can also lead to reduced infection. For a heterogeneous population comprised of risk-tolerators and risk-averters, we prove conditions on model parameters for the existence of a Hopf bifurcation and sustained oscillations in the infected population.

Strong social distancing restrictions have been crucial to controlling the COVID-19 outbreak thus far, and the next question is when and how to relax these restrictions. A sequential timing of relaxing restrictions across groups is explored in order to identify policies that simultaneously reduce health risks and economic stagnation relative to current policies. The goal will be to mitigate health risks, particularly among the most fragile sub-populations, while also managing the deleterious effect of restrictions on economic activity. The results of this paper show that a properly constructed sequential release of age-defined subgroups from strict social distancing protocols can lead to lower overall fatality rates than the simultaneous release of all individuals after a lockdown. The optimal release policy, in terms of minimizing overall death rate, must be sequential in nature, and it is important to properly time each step of the staggered release. This model allows for testing of various timing choices for staggered release policies, which can provide insights that may be helpful in the design, testing, and planning of disease management policies for the ongoing COVID-19 pandemic and future outbreaks., Comment: 22 pages (including Appendix), 10 figures

In the absence of drugs and vaccines, policymakers use non-pharmaceutical interventions such as social distancing to decrease rates of disease-causing contact, with the aim of reducing or delaying the epidemic peak. These measures carry social and economic costs, so societies may be unable to maintain them for more than a short period of time. Intervention policy design often relies on numerical simulations of epidemic models, but comparing policies and assessing their robustness demands clear principles that apply across strategies. Here we derive the theoretically optimal strategy for using a time-limited intervention to reduce the peak prevalence of a novel disease in the classic Susceptible-Infectious-Recovered epidemic model. We show that broad classes of easier-to-implement strategies can perform nearly as well as the theoretically optimal strategy. But neither the optimal strategy nor any of these near-optimal strategies is robust to implementation error: small errors in timing the intervention produce large increases in peak prevalence. Our results reveal fundamental principles of non-pharmaceutical disease control and expose their potential fragility. For robust control, an intervention must be strong, early, and ideally sustained., Comment: 23 pages plus appendix, 3 figures, 1 appendix figure

Deterministic models of vegetation often summarize, at a macroscopic scale, a multitude of intrinsically random events occurring at a microscopic scale. We bridge the gap between these scales by demonstrating convergence to a mean-field limit for a general class of stochastic models representing each individual ecological event in the limit of large system size. The proof relies on classical stochastic coupling techniques that we generalize to cover spatially extended interactions. The mean-field limit is a spatially extended non-Markovian process characterized by nonlocal integro-differential equations describing the evolution of the probability for a patch of land to be in a given state (the generalized Kolmogorov equations of the process, GKEs). We thus provide an accessible general framework for spatially extending many classical finite-state models from ecology and population dynamics. We demonstrate the practical effectiveness of our approach through a detailed comparison of our limiting spatial model and the finite-size version of a specific savanna-forest model, the so-called Staver-Levin model. There is remarkable dynamic consistency between the GKEs and the finite-size system, in spite of almost sure forest extinction in the finite-size system. To resolve this apparent paradox, we show that the extinction rate drops sharply when nontrivial equilibria emerge in the GKEs, and that the finite-size system's quasi-stationary distribution (stationary distribution conditional on non-extinction) closely matches the bifurcation diagram of the GKEs. Furthermore, the limit process can support periodic oscillations of the probability distribution, thus providing an elementary example of a jump process that does not converge to a stationary distribution. In spatially extended settings, environmental heterogeneity can lead to waves of invasion and front-pinning phenomena., Comment: 38 pages

While the benefits of common and public goods are shared, they tend to be scarce when contributions are provided voluntarily. Failure to cooperate in the provision or preservation of these goods is fundamental to sustainability challenges, ranging from local fisheries to global climate change. In the real world, such cooperative dilemmas occur in multiple interactions with complex strategic interests and frequently without full information. We argue that voluntary cooperation enabled across multiple coalitions (akin to polycentricity) not only facilitates greater generation of non-excludable public goods, but may also allow evolution toward a more cooperative, stable, and inclusive approach to governance. Contrary to any previous study, we show that these merits of multi-coalition governance are far more general than the singular examples occurring in the literature, and are robust under diverse conditions of excludability, congestability of the non-excludable public good, and arbitrary shapes of the return-to-contribution function. We first confirm the intuition that a single coalition without enforcement and with players pursuing their self-interest without knowledge of returns to contribution is prone to cooperative failure. Next, we demonstrate that the same pessimistic model but with a multi-coalition structure of governance experiences relatively higher cooperation by enabling recognition of marginal gains of cooperation in the game at stake. In the absence of enforcement, public-goods regimes that evolve through a proliferation of voluntary cooperative forums can maintain and increase cooperation more successfully than singular, inclusive regimes., Comment: 16 pages (2 figures) plus supplementary material (1 figure)

Ephemeral aggregations of bacteria are ubiquitous in the environment, where they serve as hotbeds of metabolic activity, nutrient cycling, and horizontal gene transfer. In many cases, these regions of high bacterial concentration are thought to form when motile cells use chemotaxis to navigate to chemical hotspots. However, what governs the dynamics of bacterial aggregations is unclear. Here, we use a novel experimental platform to create realistic sub-millimeter scale nutrient pulses with controlled nutrient concentrations. By combining experiments, mathematical theory and agent-based simulations, we show that individual \textit{Vibrio ordalii} bacteria begin chemotaxis toward hotspots of dissolved organic matter (DOM) when the magnitude of the chemical gradient rises sufficiently far above the sensory noise that is generated by stochastic encounters with chemoattractant molecules. Each DOM hotspot is surrounded by a dynamic ring of chemotaxing cells, which congregate in regions of high DOM concentration before dispersing as DOM diffuses and gradients become too noisy for cells to respond to. We demonstrate that \textit{V. ordalii} operates close to the theoretical limits on chemotactic precision. Numerical simulations of chemotactic bacteria, in which molecule counting noise is explicitly taken into account, point at a tradeoff between nutrient acquisition and the cost of chemotactic precision. More generally, our results illustrate how limits on sensory precision can be used to understand the location, spatial extent, and lifespan of bacterial behavioral responses in ecologically relevant environments., Comment: 6 pages, 5 figures. PNAS first published May 16, 2019 https://doi.org/10.1073/pnas.1816621116

The rate of adoption of new information depends on reinforcement from multiple sources in a way that often cannot be described by simple contagion processes. In such cases, contagion is said to be complex. Complex contagion happens in the diffusion of human behaviours, innovations, and knowledge. Based on that evidence, we propose a model that considers multiple, potentially asymmetric, and competing contagion processes and analyse its respective population-wide dynamics, bringing together ideas from complex contagion, opinion dynamics, evolutionary game theory, and language competition by shifting the focus from individuals to the properties of the diffusing processes. We show that our model spans a dynamical space in which the population exhibits patterns of consensus, dominance, and, importantly, different types of polarisation, a more diverse dynamical environment that contrasts with single simple contagion processes. We show how these patterns emerge and how different population structures modify them through a natural development of spatial correlations: structured interactions increase the range of the dominance regime by reducing that of dynamic polarisation, tight modular structures can generate structural polarisation, depending on the interplay between fundamental properties of the processes and the modularity of the interaction network., Comment: 18 Pages, 5 Figures