Abstract
The human brain1,2, power grids3, arrays of coupled lasers4 and the Amazon rainforest5,6 are all characterized by multistability7. The likelihood that these systems will remain in the most desirable of their many stable states depends on their stability against significant perturbations, particularly in a state space populated by undesirable states. Here we claim that the traditional linearization-based approach to stability is too local to adequately assess how stable a state is. Instead, we quantify it in terms of basin stability, a new measure related to the volume of the basin of attraction. Basin stability is non-local, nonlinear and easily applicable, even to high-dimensional systems. It provides a long-sought-after explanation for the surprisingly regular topologies8,9,10 of neural networks and power grids, which have eluded theoretical description based solely on linear stability11,12,13. We anticipate that basin stability will provide a powerful tool for complex systems studies, including the assessment of multistable climatic tipping elements14.
Similar content being viewed by others
Main
Complex systems science relies heavily on linear stability analysis, in which state of a dynamic system (more correctly, its dynamic regime) is assessed basically by inspecting the dominant curvature of the potential energy function in the state’s surroundings (as expressed by Lyapunov exponents). The absolute value of the curvature measures the speed of convergence or divergence after a small perturbation, and its sign qualifies the state as stable or unstable. Such linearization-based considerations are inherently local; therefore, they are not sufficient to quantify how stable a state is against non-small perturbations. Quantification of stability in this sense requires a global concept: the basin of attraction of a state is the set of initial points in state space from which the system converges to this state. Complete knowledge of the basin would allow us to fully assess the state’s stability: one could classify perturbations into the permissible and the impermissible. See Fig. 1.
However, basins are intricate entities15 and especially hard to explore in high dimensions. Here we therefore focus on a single but fundamental property: the basin’s volume. The authors of ref. 16 interpret the volume of a state’s basin of attraction as a measure of the likelihood of arrival at this state, that is, as a measure of the state’s relevance. Almost equivalently, we understand the volume of the basin as an expression of the likelihood of return to the state after any random—possibly non-small—perturbation. This yields a second interpretation: the basin’s volume quantifies how stable a state is. To the best of our knowledge, this interpretation has not yet been employed in complex systems science. We refer to the quantification of stability based on the basin’s volume as basin stability .
For climatic tipping elements14 it would be particularly useful to know how stable the desirable (that is, present) state is against perturbations. One such tipping element is the Amazon rainforest, which presumably possesses two stable states: the present fertile forest state and a barren savanna state5,6. A transition would emit huge amounts of carbon dioxide captured in the rich vegetation. Amazonian bistability arises from a positive feedback: deep-rooting trees take up water stored in the soil and transpire it to the atmosphere. Thereby, forest cover in an area increases overall precipitation and improves its own growing conditions. Consequently, a rather arid area (that is, an area with weak precipitation inflow) may still be supportive of forest growth if its forest cover exceeds a certain critical threshold; were forest cover pushed below this threshold, the area would lose all of its trees.
This is summarized in a conceptual model (see Supplementary Information),
Here, C is the relative forest cover that grows with the saturating rate r if C>Ccrit and dies with rate x (assuming r>x>0). Ccrit is the critical forest cover threshold. This model has two equilibria, the forest state CF = 1−x/r and the savanna state CS = 0. The equilibrium CF (respectively, CS) exists and is stable if CF>Ccrit (respectively, Ccrit>0). Assuming that Ccrit increases linearly with aridity A, we obtain the state diagram shown in Fig. 2.
Global warming may drive up aridity, pushing it eventually beyond the bifurcation point Acrit where the forest state vanishes. As aridity increases, the volume of the CF state’s basin of attraction shrinks, indicating that the forest state becomes less stable against perturbations. Indeed, owing to this reduced basin stability, a large perturbation such as strong deforestation might push the system to the savanna state long before aridity reaches Acrit. Crucially, none of this is reflected by linear stability, which judges the forest state CF only by the stability coefficient F′(CF) = x−r that remains constant as aridity goes up.
This implies, first, that there is no critical slowing down. Indeed, linear early-warning signals17 can be absent in systems with strong nonlinearities such as this model. Second, because of their highly local perspective, linear stability and hence the small-perturbation convergence rate do not indicate how stable CF is against non-small perturbations. They are unreliable proxies, in contradiction to previous observations17,18.
Note that both effects may also emerge in such simple mechanical systems as a damped driven pendulum (see Supplementary Information). Clearly, global stability concepts are needed.
Such global stability concepts could incorporate potential energy functions: deeper valleys would correspond to more stable states (Fig. 1). However, energy functions may not be available for many relevant (dissipative) systems. Furthermore, the estimation of energy levels on the basin boundary is numerically costly in high dimensions. An alternative approach was suggested, globally quantifying stability (or resilience, in the original nomenclature) in terms of the width of the basin of attraction in a particular direction19,20. This concept has inspired a host of studies on complex socioecological systems. Yet the measure of stability it implies “is not so easily quantified, even in models”18. One reason is that, in systems with many state variables, it is impossible to identify the single most relevant direction along which the width of the basin should be gauged. Our basin stability offers two important improvements: first, it follows a volume-based probabilistic approach that is compatible with the natural uncertainty about the strength and direction of perturbations; second, it provides a measure of stability that is clearly defined and easily quantified even in high-dimensional systems (see Methods).
Basin stability’s applicability to high-dimensional systems allows us to tackle a puzzle that has long haunted complex networks science. Researchers in this field strive to understand how a network’s topology serves its function and robustness21,22. Special effort has been put into multistable dynamic networks in which a synchronous state competes with alternative non-synchronous or partly synchronous states11,12,13,23. In a power grid, for instance, all components have to be operated at the same synchronous frequency to achieve steady power flows and to avoid damaging resonance effects3. In the brain, both neural communication24 and memory processes25,26 rely vitally on the synchronous firing of neurons. This means that, although synchronization is also associated with pathological states such as Parkinsonian tremor27, the functional ability to support synchrony is as pivotal for the brain as it is for power grids. The problem with both kinds of network is that their real-world topologies look completely different from what the theory predicts.
The theory is as follows. If the synchronous state of a dynamic network is to be maintained, it must be stable against perturbations. A groundbreaking study11 based on linearization revealed that, for a network of identical oscillators, the stability of the synchronous state can directly be inferred from the Laplacian, a matrix that reflects the coupling topology (see Supplementary Information). Indeed, for many types of oscillator, the synchronous state is stable if the ratio of the Laplacian’s maximum and minimum non-zero eigenvalues, R = λmax/λmin, is smaller than an oscillator-specific stability threshold, β = α2/α1, provided the coupling strength is chosen from the stability interval, Is = (α1/λmin,α2/λmax). The ratio R is known as the synchronizability of a network. Networks with smaller R are considered more synchronizable12. To determine what particularly synchronizable networks look like, researchers13 employed the Watts–Strogatz graph generation model8 and found that, as the model is tuned from regular lattices (model parameter p = 0) to random graphs (p = 1), synchronizability shows a strong, monotonical improvement (Fig. 3a).
Consequently, according to linear-stability-based synchronizability, real-world networks whose function relies on synchronization should ideally look like random graphs. However, neural networks and power grids exhibit small-world topologies that, from the Watts–Strogatz model’s perspective, are far more regular than random graphs8,9,10. Indeed, when building well-functioning synchronizable networks, nature and civilization seem to shun the predicted randomness. This discrepancy between theory and observation has left networks research with a long-standing puzzle.
Attempting to complement the theory, we applied basin stability to ensembles of Watts–Strogatz networks consisting of paradigmatic Rössler oscillators, in which the dynamics at node i obey
with coupling constant K, Laplacian matrix L, a = b = 0.2, and c = 7.0. Every such network has a synchronous state in which all nodes follow the same trajectory. A network’s synchronous state is stable if its synchronizability R<α2/α1 = 37.85 and K∈Is = (α1/λmin,α2/λmax), where α1 = 0.1232 and α2 = 4.663. However, the level of R does not quantify how stable the synchronous state is against perturbations. To address this yet unasked question, for each network we estimated the synchronous state’s basin stability for several K∈Is and computed the mean (see Methods). Finally, we averaged over the ensemble to obtain the expected basin stability . We found that, in sharp contrast to synchronizability, expected basin stability declines exponentially fast as networks become more random (Fig. 3a,b; for a qualitative explanation, see Supplementary Information). Therefore, the synchronous state is much more stable in networks that are more regular.
This adds a crucial piece to the puzzle and, we conjecture, makes its solution emerge (Fig. 4): in synchronizing networks, the functional need for the synchronous state to be as stable as possible promotes topological regularity. Thus, during network evolution, the optimization of synchronizability and the simultaneous optimization of basin stability have acted as two opposing forces. Their contest ended in a topological tradeoff: small-worldness.
Here, we have introduced basin stability, a new universal concept of stability. We see many important applications, notably cell regulatory networks, whose carcinogenic gene expression profiles have been related to cancer attractors in a high-dimensional multistable state space28. Like linear stability, basin stability is a property of a deterministic system and contains no information on the external perturbations that may affect it. Hence, it may often be promising to conceptually combine basin stability with a non-uniform distribution of perturbations. Basin stability could also be applied to stochastic systems by incorporating a suitable probabilistic notion of basins.
Methods
Estimation of basin stability in Watts–Strogatz networks of Rössler oscillators.
A Watts–Strogatz network is constructed as follows8: starting from a one-dimensional ring of N nodes in which every node is connected to its k nearest neighbours, each edge is rewired independently with probability p by re-choosing one of its endpoints randomly. A Watts–Strogatz network generated with p = 0 (p = 1) has a regular (random) topology. An intermediate value of p yields a small-world topology.
In a given network, the dynamics at node i (i = 1,…,N) are governed by the coupled Rössler equations (1)–(3). For K∈Is: = (0.1232/λmin,4.663/λmax), we want to estimate the volume of the synchronous state’s basin of attraction . High dimensionality poses challenges. If were a bounded convex set, its volume could be estimated in time steps with today’s best algorithm29, where n = 3N is the dimension of state space. We have N≥100, so this would be numerically very expensive. In any case, is not convex in Rössler networks (see Supplementary Information). Thus, we retreat to something feasible: we estimate the volume of in a relative sense, measuring basin stability as , where is a subset of state space that has finite volume.
More specifically, we integrate the system equations for T initial conditions drawn uniformly at random from . We count the number M of initial conditions that arrive at the synchronous state (the other possible attractor being infinity) and estimate as M/T. Observing that this is a repeated Bernoulli experiment, we infer that the estimate carries a standard error of
At T = 500, for example, e<0.023 in absolute terms.
In relative terms, for . should be chosen such that values of typically surpass this level. Our choice underlying the results presented above is with q = [−15,15]×[−15,15]×[−5,35]. Note that the Rössler attractor is included in q. We also studied other choices of —for instance ([−8,8]×[−8,8]×[−8,8])N—yet observed no qualitative difference in the outcomes. Therefore, we suppress the subscript when stating in Fig. 3.
For each network in the ensemble, we estimate basin stability for ten different equally spaced values K∈Is and average to obtain mean basin stability .
The results are not qualitatively different for networks produced by a two-dimensional Watts–Strogatz model and another model that varies the link length distribution30. Details on this are provided in the Supplementary Information.
Topological comparison of small-world networks of different sizes.
The axes of Fig. 4 were chosen so that different real-world networks can be compared, without the distorting effects of network size, with respect to average shortest path L and clustering coefficient C. L and C have been widely used to characterize small-worldness8. We plot ξL against ξC, where ξX = 1−log(X/XR)/log(XL/XR) with X = L or C. LR,CR (or LL,CL) are the values of the respective quantities in random networks (or regular lattices) of the same size. ξX counts how many orders of magnitude X is away from XL in relation to the count of orders of magnitude between XL and XR. This way (ξL,ξC) = (0,0) labels a regular network whereas (ξL,ξC) = (1,1) labels a random network. Small-world networks have8 L≈LR and C≫CR and therefore reside in the top-left quadrant. See Supplementary Information.
References
Babloyantz, A. & Destexhe, A. Low-dimensional chaos in an instance of epilepsy. Proc. Natl. Acad. Sci. USA 83, 3513–3517 (1986).
Lytton, W. W. Computer modelling of epilepsy. Nature Rev. Neurosci. 9, 626–637 (2008).
Machowski, J., Bialek, J. W. & Bumby, J. R. Power System Dynamics: Stability and Control (Wiley, 2008).
Erzgräber, H. et al. Mutually delay-coupled semiconductor lasers: Mode bifurcation scenarios. Opt. Commun. 255, 286–296 (2005).
Da Silveira Lobo Sternberg, L. Savanna-forest hysteresis in the tropics. Glob. Ecol. Biogeogr. 10, 369–378 (2001).
Hirota, M., Holmgren, M., Van Nes, E. H. & Scheffer, M. Global resilience of tropical forest and savanna to critical transitions. Science 334, 232–235 (2011).
May, R. M. Thresholds and breakpoints in ecosystems with a multiplicity of stable states. Nature 269, 471–477 (1977).
Watts, D. J. & Strogatz, S. H. Collective dynamics of ‘small-world’ networks. Nature 393, 440–442 (1998).
Sporns, O. & Zwi, J. The small world of the cerebral cortex. Neuroinformatics 2, 145–162 (2004).
Wang, Z., Scaglione, A. & Thomas, R. J. Generating statistically correct random topologies for testing smart grid communication and control networks. IEEE Trans. Smart Grid 1, 28–39 (2010).
Pecora, L. M. & Carroll, T. L. Master stability functions for synchronized coupled systems. Phys. Rev. Lett. 80, 2109–2112 (1998).
Barahona, M. & Pecora, L. M. Synchronization in small-world systems. Phys. Rev. Lett. 89, 054101 (2002).
Hong, H., Kim, B. J., Choi, M. Y. & Park, H. Factors that predict better synchronizability on complex networks. Phys. Rev. E 69, 067105 (2004).
Lenton, T. M. et al. Tipping elements in the earth’s climate system. Proc. Natl. Acad. Sci. USA 105, 1786–1793 (2008).
Nusse, H. E. & Yorke, J. A. Basins of attraction. Science 271, 1376–1380 (1996).
Wiley, D. A., Strogatz, S. H. & Girvan, M. The size of the sync basin. Chaos 16, 015103 (2006).
Scheffer, M. et al. Early-warning signals for critical transitions. Nature 461, 53–59 (2009).
Van Nes, E. H. & Scheffer, M. Slow recovery from perturbations as a generic indicator of a nearby catastrophic shift. Am. Nature 169, 738–747 (2007).
Holling, C. S. Resilience and stability of ecological systems. Annu. Rev. Ecol. Syst. 4, 1–23 (1973).
Scheffer, M. Critical Transitions in Nature and Society (Princeton Univ. Press, 2009).
Boccaletti, S., Latora, V., Moreno, Y., Chavez, M. & Hwang, D-U. Complex networks: Structure and dynamics. Phys. Rep. 424, 175–308 (2006).
Buldyrev, S. V., Parshani, R., Paul, G., Stanley, H. E. & Havlin, S. Catastrophic cascade of failures in interdependent networks. Nature 464, 1025–1028 (2010).
Arenas, A., Díaz-Guilera, A., Kurths, J., Moreno, Y. & Zhou, C. Synchronization in complex networks. Phys. Rep. 469, 93–153 (2008).
Fries, P. A mechanism for cognitive dynamics: neuronal communication through neuronal coherence. Trends Cogn. Sci. 9, 474–480 (2005).
Fell, J. et al. Human memory formation is accompanied by rhinal-hippocampal coupling and decoupling. Nature Neurosci. 4, 1259–1264 (2001).
Fell, J. & Axmacher, N. The role of phase synchronization in memory processes. Nature Rev. Neurosci. 12, 105–118 (2011).
Hammond, C., Bergman, H. & Brown, P. Pathological synchronization in Parkinson’s disease: Networks, models and treatments. Trends Neurosci. 30, 357–364 (2007).
Huang, S. & Ingber, D.E. A non-genetic basis for cancer progression and metastasis: Self-organizing attractors in cell regulatory networks. Breast Disease 26, 27–54 (2007).
Lovász, L. & Vempala, S. Simulated annealing in convex bodies and an volume algorithm. J. Comput. Syst. Sci. 72, 392–417 (2006).
Li, G. et al. Towards design principles for optimal transport networks. Phys. Rev. Lett. 104, 018701 (2010).
Acknowledgements
The authors acknowledge financial support from IRTG 1740 (Deutsche Forschungsgemeinschaft), the SUMO-project (European Union), the ECONS-project (Leibniz Association) and Konrad-Adenauer-Stiftung. They thank N. Fujiwara and A. Rammig for inspiring discussions.
Author information
Authors and Affiliations
Contributions
P.J.M. conceived the study, performed the numerics, and prepared the manuscript. All authors discussed the results, drew conclusions and edited the manuscript. J.K. supervised the study.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Information (PDF 1035 kb)
Rights and permissions
About this article
Cite this article
Menck, P., Heitzig, J., Marwan, N. et al. How basin stability complements the linear-stability paradigm. Nature Phys 9, 89–92 (2013). https://doi.org/10.1038/nphys2516
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys2516
This article is cited by
-
Fractal basins as a mechanism for the nimble brain
Scientific Reports (2023)
-
Structure of parameter space of a three-species food chain model with immigration and emigration
Nonlinear Dynamics (2023)
-
Complex dynamics of a conceptual airfoil structure with consideration of extreme flight conditions
Nonlinear Dynamics (2023)
-
Vulnerable Waters are Essential to Watershed Resilience
Ecosystems (2023)
-
Alternative stable ecological states observed after a biological invasion
Scientific Reports (2022)