In spatially heterogeneous landscapes, some habitats may be persistent sources, providing immigrants to sustain populations in unfavorable sink habitats (where extinction is inevitable without immigration). Recent theoretical and empirical studies of source-sink systems demonstrate that temporally variable local growth rates in sinks can substantially increase average abundance of a persisting population, provided that the variation is positively autocorrelated--in effect, temporal variation inflates average abundance. Here we extend these results to a metapopulation in which all habitat patches are sinks. Using numerical studies of a population with discrete generations (buttressed by analytic results), we show that temporal variation and moderate dispersal can jointly permit indefinite persistence of the metapopulation and that positive autocorrelation both lowers the magnitude of variation required for persistence and increases the average abundance of persisting metapopulations. These effects are weakened--but not destroyed--if variation in local growth rates is spatially synchronized and dispersal is localized. We show that the inflationary effect is robust to a number of extensions of the basic model, including demographic stochasticity and density dependence. Because ecological and environmental processes contributing to temporally variable growth rates in natural populations are typically autocorrelated, these observations may have important implications for species persistence.

Robust critical systems are characterized by power laws which occur over a broad range of conditions. Their robust behaviour has been explained by local interactions. While such systems could be widespread in nature, their properties are not well understood. Here, we study three robust critical ecosystem models and a null model that lacks spatial interactions. In all these models, individuals aggregate in patches whose size distributions follow power laws which melt down under increasing external stress. We propose that this power-law decay associated with the connectivity of the system can be used to evaluate the level of stress exerted on the ecosystem. We identify several indicators along the transition to extinction. These indicators give us a relative measure of the distance to extinction, and have therefore potential application to conservation biology, especially for ecosystems with self-organization and critical transitions.

Three different lattice-based models for antagonistic ecological interactions, both nonlinear and stochastic, exhibit similar power-law scalings in the geometry of clusters. Specifically, cluster size distributions and perimeter-area curves follow power-law scalings. In the coexistence regime, these patterns are robust: their exponents, and therefore the associated Korcak exponent characterizing patchiness, depend only weakly on the parameters of the systems. These distributions, in particular the values of their exponents, are close to those reported in the literature for systems associated with self-organized criticality (SOC) such as forest-fire models; however, the typical assumptions of SOC need not apply. Our results demonstrate that power-law scalings in cluster size distributions are not restricted to systems for antagonistic interactions in which a clear separation of time-scales holds. The patterns are characteristic of processes of growth and inhibition in space, such as those in predator-prey and disturbance-recovery dynamics. Inversions of these patterns, that is, scalings with a positive slope as described for plankton distributions, would therefore require spatial forcing by environmental variability.

BackgroundWith over a hundred million annual infections and rising morbidity and mortality, Plasmodium vivax malaria remains largely a neglected disease. In particular, the dependence of this malaria species on relapses and the potential significance of the dormant stage as a therapeutic target, are poorly understood.Methodology/Principal FindingsTo quantify relapse parameters and assess the population-wide consequences of anti-relapse treatment, we formulated a transmission model for P. vivax suitable for parameter inference with a recently developed statistical method based on routine surveillance data. A low-endemic region in NW India, whose strong seasonality demarcates the transmission season, provides an opportunity to apply this modeling approach. Our model gives maximum likelihood estimates of 7.1 months for the mean latency and 31% for the relapse rate, in close agreement with regression estimates and clinical evaluation studies in the area. With a baseline of prevailing treatment practices, the model predicts that an effective anti-relapse treatment of 65% of those infected would result in elimination within a decade, and that periodic mass treatment would dramatically reduce the burden of the disease in a few years.Conclusion/SignificanceThe striking dependence of P. vivax on relapses for survival reinforces the urgency to develop more effective anti-relapse treatments to replace Primaquine (PQ), the only available drug for the last fifty years. Our methods can provide alternative and simple means to estimate latency times and relapse frequency using routine epidemiological data, and to evaluate the population-wide impact of relapse treatment in areas similar to our study area.