arxiv: v2 [physics.soc-ph] 22 Apr 2018

Transcription

1 Resilience of networks to environmentl stress: From regulr to rndom networks Young-Ho Eom Deprtment of Mthemtics nd Sttistics, University of Strthclyde, 26 Richmond Street, Glsgow G1 1H, United Kingdom nd Deprtmento de Mtemátics, Universidd Crlos III de Mdrid, Legnés, Spin riv: v2 [physics.soc-ph] 22 Apr 218 Despite the huge interest in network resilience to stress, most of the studies hve concentrted on internl stress dmging network structure (e.g., node removls). Here we study how networks respond to environmentl stress deteriorting their externl conditions. We show tht, when regulr networks grdully disintegrte s environmentl stress increses, disordered networks cn suddenly collpse t criticl stress with hysteresis nd vulnerbility to perturbtions. We demonstrte tht this difference results from trde-off between node resilience nd network resilience to environmentl stress. The nodes in the disordered networks cn suppress their collpses due to the smll-world topology of the networks but eventully collpse ll together in return. Our findings indicte tht some rel networks cn be highly resilient ginst environmentl stress to threshold yet extremely vulnerble to the stress bove the threshold becuse of their smll-world topology. PACS numbers: , Cc, Hc I. INTRODUCTION Networks usully mintin their essentil functions owing to their resilience lthough they re under constnt stress such s dmges, filures, nd environmentl chnges [1 3]. But networks cn lose their resilience nd collpse when stress exceeds criticl level (i.e., threshold). How networks lose their resilience relies on their topology. Homogeneous networks re frgmented when the frction of rndomly removed nodes or links reches brekdown thresholds wheres scle-free networks re robust to such removls [3 5]. The disintegrting pttern of networks t thresholds lso depends on their topology. When nodes fil, brupt trnsitions to disintegrted stte cn rise in interdependent networks but smooth trnsitions occur in isolted networks [6 1]. Identifying the reltion between the response of networks to stress nd their topology is vitl tsk for their stble function. Stress cn cuse normlly two types of effects on networks. Accordingly one cn define two types of stress. Internl stress dmges network structure nd cuses chnges in the weighted djcency mtrix of networks (e.g., loss of nodes, links, or coupling strength), wheres environmentl stress deteriortes externl conditions of networks (e.g., chnges in externl prmeters determining the stte of nodes). Both types of stress cn hve profound impct on networks. In network of ptchy ecosystems, for instnce, not only loss of ptches but lso climte chnges nd diminishing resources cn be serious threts. In the Internet or power grids, growing demnd in the systems cn cuse the systemic collpse s the locl filures cn. However, the role of network topology in network resilience to environmentl stress hs been neglected, in contrst to the gret ttention to internl stress [2 14], despite existing works on net- Electronic ddress: youngho.eom@strth.c.uk works under environmentl chnges [15 18]. For exmple, recent studies on mutulistic networks of intercting species showed how such networks behve when interction strength [16, 17] or growth rte [18] re reduced by environmentl chnges. The role of their nested network topology ws lso highlighted [18]. But these works hve focused on mutulistic networks, specific type of ecologicl networks. How other, possibly more generl, networked systems respond to environmentl stress remins n open question. To ddress this, we consider simple network model where ll nodes re subject to common externl condition. Ech node, when isolted, cn chnge bruptly its stte (e.g., from ctive to collpsed) vi sddle-node bifurction if the condition crosses its own threshold. But the coupling mong the nodes cn ffect their collective response to the condition. This model is motivted by the models for ecosystems exposed to environmentl stress. Ecosystems cn undergo sudden stte chnges when environmentl stress exceeds ctstrophic thresholds, often clled tipping points [19 21]. Recent studies show tht tipping points rise in vrious systems rnging from ecosystems nd popultions [22 29] to climte systems [3] nd finncil mrkets [31]. To nticipte whether these lrge-scle systems respond to environmentl stress bruptly or grdully, mny studies hve considered theoreticl models composed of connected components which shift their sttes vi bifurction in response to chnging externl conditions [32 39]. In prticulr, some of these studies reported tht brupt responses to environmentl chnges cn be smoothened if there exist strong heterogeneity in the response of the components, limited coupling mong the components, or sufficient noise [38, 39] (But note tht recent study [4] climed tht there is bis in the numeric simultions of the Ref. [39].) These models hve ssumed tht their connection topology mong the components is perfectly regulr, wheres most rel networks hve disorder in their topology.

2 2 In this work we study the effects of topologicl disorder on the response of networks to environmentl stress. Specificlly, we investigte whether topologicl disorder cn lter the nture of response s the heterogeneity or the coupling cn. We show tht s the externl condition deteriortes, regulr networks grdully crumble, wheres disordered networks suddenly collpse into distinct sttes with strong hysteresis nd vulnerbility to perturbtions. This mens tht tipping points rise in the disordered networks. We demonstrte tht the nodes in the disordered networks cn keep their functions owing to the networkwide support cused by the smll-world topology of the networks but ultimtely they collpse ll together in return. This indictes tht there is trdeoff, determined by network topology, between node resilience nd network resilience to environmentl stress nd rel networks my suddenly collpse under environmentl stress owing to their smll-world topology. II. THE NETWORK MODEL A. Nodes: Miniml systems We tke cusp ctstrophe system for the dynmics of the nodes in our model. It is miniml mthemticl model for systems chnging bruptly their sttes vi sddle-node bifurction [41, 42]. The stte of this miniml system is evolved by the following differentil eqution dx dt = + bx x3, (1) where x is the stte of the system (e.g., ctivity). We consider the system s ctive (collpsed) for x > (x < ). The prmeter nd b denote the environmentl stress on the system nd the growth rte of the system, respectively. The cubic term x 3 is to prevent the stte from x () b< x - T2 (b) b> (ii) (i) T FIG. 1: (color online). Bifurction digrms of the miniml system whose dynmics is described by Eq. (1) with () b = nd (b) b =. The red lines represent the stble sttes nd the dotted gry line represents the unstble sttes. The blue filled circles (i.e., T 1 nd T 2) in (b) denote the tipping points where sddle-node bifurction occurs. Two wys of ctstrophic shifts between distinct stble sttes re illustrted in (b): (i) the prmeter is pssed the thresholds T 1 (the horizontl rrow) nd (ii) the system s stte is crossed the unstble sttes by perturbtions (the blck verticl rrow). diverging. This miniml system cn be considered s men-field model for systems with tipping points becuse sddle-node bifurction cn correspond to the simplest form of ctstrophic shifts [43 45]. Figure 1 represents the bifurction digrms of the system. When b <, the equilibrium stte curve is grdul (Fig. 1()). When b >, the equilibrium stte curve is folded t the tipping points (i.e., T 1 nd T 2 ) where one of the stble sttes disppers vi sddle-node bifurction (Fig. 1(b)). If the system is close to tipping point, then smll chnges in the stress (the horizontl rrow in Fig. 1(b)) or in the stte (the blck verticl rrow in Fig. 1(b)) cn rdiclly move the system to nother stble stte (e.g., from ctive to collpsed). Once the system collpsed t T 1, we should restore the stress to T 2 < T 1 for the recovery (i.e., hysteresis). B. Networks of the miniml systems The stte of node i in network of N miniml systems is evolved by dt = + b ix i x i 3 + D k i A ij (x j x i ), (2) where x i is the stte of node i (i = 1,..., N). The prmeter is the environmentl stress nd b i is the growth rte of node i rndomly drwn from uniform distribution B. D is the coupling strength, k i is the number of connections (i.e., the degree) of node i, nd A ij is the element of the djcency mtrix of the network which is unity when nodes i nd j re connected nd zero otherwise. The coupling term denotes the dispersion of ctivity mong the connected nodes [35, 38]. The network s stte is defined s = 1 N i x i. The network is ctive (collpsed) if > ( < ) s the nodes. We consider networks with the sme verge degree k = (1/N) i k i nd the sme size N = 1 4 generted by three network models: the Wtts-Strogtz (WS) model [46], the Erdős- Réyni () model [47], nd the Brbási-Albert () model [48]. The WS model llows us to construct networks rnging from regulr lttice (P = ) to disordered lttice (P = 1) by tuning link rewiring probbility P (i.e., the topologicl disorder in the networks). We consider two-dimensionl lttices with periodic boundries for the WS model. The model nd the model generte rndom networks with Poisson degree distributions nd with power-lw degree distributions, respectively. III. RESULTS Here we im t checking whether the grdul response of regulr networks resulting from their strong heterogeneity or limited coupling cn chnge to n brupt response when the disorder in network topology increses. Accordingly we choose the prmeter vlues D =,

3 3 - P= P=5 P=.1 - () P= P= - (b) (c) P= (d) (e) (f) FIG. 2: (color online). Stte digrms () of the networks of the miniml systems. () The Wtts-Strogtz (WS) networks with P.1 respond grdully to the incresing stress. (b) The WS networks with P.25, the Erdős-Réyni (), nd the Brbási-Albert () respond bruptly to the incresing stress. Incresing nd then decresing the stress, we observed hysteresis in (c) the WS network with P =, (d) the WS network with P =.25, (e) the network, nd (f) the network. Ech line with symbols is verged over 5 reliztions. The dotted lines in (c)-(f) show tht the hysteresis cn be mitigted in the WS network with P = but not in the WS network with P =.25, the network, nd the network. Ech dotted line shows single reliztion obtined by incresing the stress from =.8 nd then decresing the stress when it reches =.41,.43,.45,.47,.49, 1, 3, or 5. B = [, 3], nd k = 4 in our numeric simultions. Such prmeter vlues generte networks with sufficient heterogeneity (i.e., brod B) nd limited coupling (i.e., low D nd k ), so tht one cn observe grdul response to environmentl stress in regulr networks (i.e., the WS network with P 1). But note tht, lthough networks re perfectly regulr, they cn respond bruptly to environmentl stress if we consider the higher vlue of D, the nrrower B, or the higher k s predicted in the previous studies [38, 39] (see Appendix A for the observed brupt responses of the regulr networks in such cses). Here we did not consider noise to focus on the role of network topology in the response to environmentl stress. Thus the pplicbility of our results is limited to wekly noisy situtions becuse noise cn ffect the nture of the response [39]. We focus on synthetic networks becuse we cn tune the topologicl disorder. All results re obtined by fourth-order Runge-Kutt method with time step of 2. To identify how these networks respond to the environmentl stress, we report the stte digrm () of the - - networks. We consider deteriorting scenrio where the sttionry vlue of () is obtined by incresing the stress with n increment δ = 1 from =.8 to =.8. At =.8, initilly ll nodes hve the sme ctive stte = to obtin the stedy stte. Subsequently the stedy stte of the previous step is ssigned s the initil condition of the next step. As depicted in Figs. 2() nd 2(b), ll the networks under the deteriorting scenrio disintegrte with certin threshold effects. The shpe of their stte curves, however, vries depending on their topology. The regulr networks (i.e., the WS networks with P 1) respond grdully to the incresing stress s expected in the previous works for regulr lttices [38, 39]. The stte curves, however, become brupt s P increses (Fig. 2()). The stte curves of the disordered networks (i.e., the WS network with P.25, the network, nd the network) re order of mgnitude more brupt thn the curve of the WS network with P = (Fig. 2(b)). We lso consider n meliorting scenrio where the sttionry vlue of () is obtined by decresing the stress with decrement δ = 1 from =.8 to =.8 fter the deteriorting scenrio. Hysteresis is observed in ll the networks (Figs. 2(c)-2(f)). But we cn mitigte the hysteresis in the regulr networks if we strt to decrese the stress before it reches =.8 (The dotted lines in Fig. 2(c)). No such mitigtion is observed in the disordered networks (The dotted lines in Figs. 2(d)-2(f)). This observtion suggests tht lthough the stress exceeds the thresholds, mny nodes re still ctive in the regulr networks wheres most nodes immeditely collpse in the disordered networks. The bove results imply tht the thresholds in the disordered networks re tipping points wheres the ones in the regulr networks re not. To clrify the nture of the thresholds, we estimte the sttionry stte distribution of ech network which mnifests the network s ttrctor. We generte lrge number of sttionry sttes of the networks with the stress nd the initil stte rndomly drwn from uniform distributions [.4,.6] nd [ 1.6, 1.6], respectively. The stte is ssigned to ll the nodes initilly. The probbility distribution of the generted sttionry sttes in ech network indictes tht the ttrctors of the regulr networks (Figs. 3()- 3(c)) re fundmentlly different from those of the disordered networks (Figs. 3(d)-3(f)). In the regulr networks, the upper brnch ttrctors continue to exist fter the thresholds. But the density of sttionry sttes hving the vlues from to.2 becomes lower s P increses in the WS networks (Figs. 3(b) nd 3(c)). In the disordered networks, the upper ttrctors suddenly dispper t the thresholds (Figs. 3(d)-(f)). This shows tht the brupt stte chnges in the disordered networks re lso trnsitions between distinct ttrctors. Another chrcteristic of tipping points is vulnerbility to perturbtions s depicted in Fig. 1(b). Smll stte chnges cn trigger to move systems to contrsting stte ner ctstrophic thresholds. Cn we expect the

4 4 - - () P= (c) P= (e) (b) P= (d) (f) FIG. 3: (color online). Sttionry stte distributions of the networks of the miniml systems. The Wtts-Strogtz networks with () P =, (b) P = 5, (c) P =.1, nd (d) P =.25. (e) The Erdős-Réyni network. (f) The Brbási- Albert network. For ech pnel, we generte sttionry sttes with the stress nd the initil condition rndomly selected from uniform distributions [.4,.6] nd [ 1.6, 1.6], respectively. ()-(c) The upper ttrctors exist fter the thresholds but the density of the sttionry sttes hving < <.2 decreses s P increses in the regulr networks. (d)-(f) The upper ttrctors dispper t the thresholds in the disordered networks. sme vulnerbility in the disordered networks? To nswer this, we perform perturbtion experiments numericlly. We obtin sttionry sttes B of network for given stress nd the initil condition rndomly chosen from uniform distribution [, 1.6] nd ssigned eqully to ll the nodes. If B > (i.e., the network is ctive), then we give smll shock S to the network nd obtin the sttionry stte A fter the perturbtion. We consider two types of perturbtions. For the globl perturbtion, we give the sme shock S G to ll the nodes wheres we give the shock S L to rndomly chosen 1% of the nodes for the locl perturbtion. Figure 4 represents the distributions of stte chnges = A B in the networks fter the perturbtions for the given stress nd the perturbtions. Regrdless of the perturbtion types, we found consistent pttern. The regulr networks re resilient ginst the perturbtions in ll the considered cses. The disordered networks become more likely to shift to different sttes fter the perturbtions s the stress pproches their thresholds. Our findings showed tht the regulr nd the disordered networks show qulittively different responses to environmentl stress. Where does the difference come from? The coupling term of Eq. (2) cn bring node the P( ) P( ) P( ) P( ) 1 () =.44, Globl (b) =.45, Globl (c) =.46, Globl 1-1 P= (d) =.44, Locl (e) =.45, Locl (f) =.46, Locl (g) =.44, Globl (h) =.45, Globl (i) =.46, Globl (j) =.44, Locl (k) =.45, Locl (l) =.46, Locl FIG. 4: (color online). Distributions of stte chnges in the networks of the miniml systems fter perturbtions. ()- (c) Globl perturbtions (S G =.1) were given to ll the nodes in the Wtts-Strogtz (WS) networks with P = nd P =.25. (d)-(f) Locl perturbtions (S L = ) were given to rndomly selected 1% of the nodes in the WS networks with P = nd P =.25. (g)-(i) Globl perturbtions (S G =.1) were given to ll the nodes in the Erdős-Réyni () network nd the Brbási-Albert () network. (j)-(l) Locl perturbtions (S L = ) were given to rndomly selected 1% of the nodes in the nd the networks. subsidy of ctivity from its neighboring nodes. This subsidy cn suppress the collpse of individul nodes [45]. In the regulr networks, the nodes cn get only limited subsidy since the number of their neighboring nodes is few. However, the disordered networks hve smll-world topology [46] where the nodes cn get network-wide subsidy becuse the verge distnce between the nodes is substntilly short. Hence most nodes cn suppress their collpses but eventully ll the nodes collpse together. To vlidte such trde-off between node nd network resilience, we check the rtio of ctive nodes r + to the entire popultion s function of the stress under the deteriorting scenrio. As depicted in Fig. 5(), the nodes in the regulr networks strted to crumble grdully t the lower thresholds wheres those in the disordered networks collpsed simultneously t the higher thresholds. To check the effect of reduced distnce, we divide WS network into four equivlent sublttices before rndom link rewiring nd rewire only the links between the nodes in one of the sublttices. This loclly rewired (LR) network hs the sme number of rewired links s the originl WS network with P =.25 hs, yet the verge distnce of the former is longer thn tht of the ltter since the loclly rewired links cnnot provide long-rnge short-cuts. The nodes in the LR network grdully collpse (Fig. 5()). We generte sttionry sttes of the LR network to obtin its ttrctor with the stress nd the initil stte rndomly drwn from uniform distributions [.4,.6] nd [ 1.6, 1.6], respectively. The upper ttrctors of the LR network hs smll gp indicting the collpse of the rewired sub-lttice (Fig. 5(b)). Note

5 5 () r + P= LR (b) FIG. 5: (color online). Suppressed collpse of individul nodes in the networks of the miniml systems. () The rtio of ctive nodes r + s function of the stress shows the synchronized collpses of the nodes in the Wtts-Strogtz (WS) network with P =.25, the Erdős-Réyni network, nd the Brbási-Albert network but the grdul collpses of the nodes in the WS network with P = nd the loclly rewired (LR) networks. To generte the LR network, we divide WS network into four equivlent sub-lttice before link rewiring nd rewire only the links between the nodes in one of the sub-lttice. Thus the sme number of links is rewired in the LR network s in the WS network with P =.25. Ech line is verged over 5 reliztions. (b) The sttionry stte distribution of the LR network shows smll gp indicting the collpse of the rewired sub-lttice in the upper ttrctor. tht suppressing the onset of trnsitions lso leds to other brupt trnsitions in networks such s explosive percoltion [49] nd explosive synchroniztion [5]. Although we provide results only for the cse of D = nd B = [, 3], nd we consider the cusp ctstrophe system s the dynmicl model of the nodes here, our results re more generl. To support the generlity, we consider three dditionl cses. First we consider the sme popultions on the model networks but with the higher coupling strength D = 1.2 (See Appendix B). Second we consider the mixed miniml systems whose growth rte b is rndomly drwn from uniform distribution B = [, ] (see Appendix C) rther thn B = [, 3] of the bove nlysis. As shown in Fig. 1, the miniml system responds bruptly when b > but grdully when b <. Hence the hlf of the nodes responds grdully in these mixed popultions when they re isolted. Third we consider the grzing system [51, 52], prdigmtic model for logisticlly growing vegettion under grzing stress (see Appendix D). The stte of node i in network of N grzing systems is evolved by dt = x i(1 x i ) x2 i b i x 2 i D k i A ij (x j x i ), (3) where x i ( ) is the biomss density of node i, is the grzing stress, b i is the crrying cpcity of node i rndomly drwn from uniform distribution B = [5., 1], D is the coupling strength (D =.3). Hence the grzing systems with ecologicl implictions re more relistic thn the miniml systems. We consider these three dditionl cses on networks with N = 1 4 nd k = 4. As depicted in Appendices B, C, nd D, we obtin qulittively the sme results demonstrting tht, regrdless of the some chnges in the coupling strength, the heterogeneity, or the dynmicl model of the nodes, the topologicl disorder of the networks cn mke their response to environmentl stress more brupt nd led to tipping points by inducing smll-world topology in the networks. Lst one my wonder the impct of the lttice dimension on the response of regulr networks s we considered only the cse of two-dimensionl lttice. We cn infer the impct of lttice dimension from our findings. The verge distnce, d r, mong the nodes in regulr lttices is given by d r N 1/α, where N is the number of the nodes in the lttice nd α is the lttice dimension [1]. Hence, s α increses, the verge distnce decreses. Therefore, it is expected tht high-dimensionl regulr networks cn respond more bruptly to environmentl stress thn low-dimensionl regulr networks. IV. DISCUSSION AND CONCLUSIONS In conclusion, using simple model, we hve investigted the role of network topology in the response of networks to environmentl stress, which hs profound impcts on network resilience. We showed tht network topology cn determine how networks under environmentl stress collpse t thresholds by regulting the trdeoff between node resilience nd network resilience to such stress. In the disordered networks, when the environments deteriorte, the nodes cn persist in functioning owing to the networkwide subsidy cused by the smllworld topology of the networks but eventully the whole networks bruptly collpse t tipping points in return. In the regulr networks, some nodes fil in response to smll environmentl stress without the networkwide support becuse the networks re not smll-world. Hence the regulr networks grdully crumble s the stress increses. These results revel the importnt role of smll-world topology in network resilience to environmentl stress. Note tht network resilience to internl stress is minly influenced by degree distributions or interdependency between networks [3 14]. Such difference between the response to environmentl stress nd the response to internl stress indicte tht network my respond grdully to one type of stress but it cn respond bruptly to the other stress. For exmple, network is known to be robust to the rndom filures of its nodes owing to its brod degree distribution wheres it cn collpse bruptly s result of smll chnges in environmentl stress becuse of its smll-world topology. Therefore it is necessry to tke into ccount the impct of both internl nd environmentl stress for deeper understnding of network resilience. Our findings indicte tht some rel networks cn be highly resilient ginst environmentl stress to threshold yet extremely vulnerble to the stress bove the threshold becuse they hve smll-world topology. In such cses, the resilience to certin level of the stress cnnot gurntee the resilience to the slightly lrger stress. Insted, network cn collpse bruptly t the

6 6 () D=.8 (b) D= (c) D= (d) D=1.4 (e) D=1.6 (f) D= () B=[.1].2.4 (d) B=[,4] (b) B=[,2] (e) B=[,5] (c) B=[,3] FIG. 6: (color online). Stte digrms () of the regulr networks of the miniml systems. All networks re generted from the two-dimensionl Wtts-Strogtz model with P =. As the coupling strength D increses, the response of network sttes to environmentl stress becomes brupt. The prmeter vlues re N = 1 4, B = [, 3], nd k = 4. () D =.8. (b) D =. (c) D = 1.2. (d) D = 1.4. (e) D = 1.6. (f) D = 1.8. Ech curve shows single reliztion. Ten reliztions re represented for ech pnel. FIG. 7: (color online). Stte digrms () of the regulr networks of the miniml systems. All networks re generted from the two-dimensionl Wtts-Strogtz model with P =. As the width of the intervl B decreses, the response of network sttes to environmentl stress becomes brupt. The prmeter vlues re N = 1 4, D =, nd k = 4. () B = [, 1]. (b) B = [, 2]. (c) B = [, 3]. (d) B = [, 4]. (e) B = [, 5]. Ech curve shows single reliztion. Ten reliztions re represented for ech pnel. slightly lrger stress. Erly-wrning indictors for such vulnerbility need to be developed becuse it is difficult to reverse once the networks collpsed. Criticl slowing down (i.e., slower recovery from perturbtions in vicinity of bifurction points) nd its relted indictors hve been suggested s erly-wrning indictors for ctstrophic shifts in some ecosystems, microorgnism popultions, nd mutulistic networks [17, 2, 21, 26 28]. This prompts us to check whether similr phenomenon cn rise in our model to foresee the sudden collpse of networks s future work. Lst our results suggest tht whether growing environmentl stress is crucil thret or not cn depend on other topologicl properties of networks. In fct, rel networks hve not only smll-world topology but lso other common properties such s degree-correltions, community structure, nd brod degree distributions, to nme few [1]. We wish to study how such chrcteristics determine the response of rel networks to environmentl stress to design or reshpe network topology for more resilient networks s future works. Appendix A: The response of regulr networks with different prmeters We consider networks of the sme miniml systems defined in the min text. Hence the stte of node i in network of N miniml systems is evolved by dt = + b ix i x i 3 + D k i A ij (x j x i ), (A1) where x i is the stte of node i, is the environmentl stress, b i is the growth rte of node i drwn from uniform distribution B, D is the overll coupling strength mong the nodes, k i is the degree of node i, nd A ij is the element of the djcency mtrix of the network. The net- work s stte is defined s = 1 N i x i. In this section we focus on the response of regulr networks to environmentl stress. Regulr networks cn respond bruptly to environmentl stress when the heterogeneity in the response of the individul components is wek (i.e., the B is nrrow) or the connectivity between the components re high (i.e., the coupling strength D or the verge degree k re high) [38, 39]. We consider deteriorting scenrio where the sttionry vlue of () is obtined by incresing the stress with n increment δ = 1 from =.8 to =.8. We observed tht the response of the networks generted by the Wtts-Strogtz model with probbility P = gets more brupt when the coupling strength D increses (Fig. 6), the heterogeneity of the growth rtes decreses (Fig. 7), or the verge degree - - () <k>=4.3.4 (d) <k>= (b) <k>=1.3.4 (e) <k>= (c) <k>=2.3.4 (f) <k>= FIG. 8: (color online). Stte digrms () of the regulr networks of the miniml systems. All networks re generted from the one-dimensionl Wtts-Strogtz model with P =. Note tht, unlike the other cses, we consider onedimensionl Wtts-Strogtz model to chnge the verge degree in the networks. As the verge degree increses, the response of network sttes to environmentl stress becomes brupt. The prmeter vlues re N = 1 4, B = [, 3], nd D =. () k = 4. (b) k = 1. (c) k = 2. (d) k = 4. (e) k = 6. (f) k = 8. Ech curve shows single reliztion. Ten reliztions re represented for ech pnel.

7 7 - P= P=5 P=.1 - () P= P= - (b) (c) P= (d) (e) (f) FIG. 9: (color online). Stte digrms () of the networks of the miniml systems with D = 1.2. () The Wtts-Strogtz (WS) network with P = responds grdully to the incresing stress. But, s P increses, the response becomes more brupt. (b) The WS networks with P.25, the Erdős- Réyni (), nd the Brbási-Albert () respond bruptly to the incresing stress. Incresing nd then decresing the stress, we observed hysteresis ptterns in (c) the WS network with P =, (d) the WS network with P =.25, (e) the network, nd (f) the network. Ech line with symbols is verged over 5 reliztions. The dotted lines in (c)-(f) show tht the hysteresis cn be mitigted in the WS network with P = but not in the WS network with P =.25, the network, nd the network. Ech dotted line shows single reliztion obtined by incresing nd then decresing the stress when it reches =.41,.43,.45,.47,.49, 1, 3, or 5. increses (Fig. 8). Note tht we consider one-dimensionl Wtts-Strogtz model in Fig. 8 unlike in Figs. 6 nd 7 since one-dimensionl cse is esy to control the verge degree thn two-dimensionl cse in the min text. Appendix B: Networks of the miniml systems with different coupling strength We consider networks of the sme miniml systems defined in the min text. Hence the stte of node i in network of N miniml systems is evolved by dt = + b ix i x i 3 + D k i A ij (x j x i ), (B1) where x i is the stte of node i, is the environmentl stress, b i is the growth rte of node i drwn from uniform distribution B, D is the overll coupling strength mong the nodes, k i is the degree of node i, nd A ij is the element of the djcency mtrix of the network. The network s stte is defined s = 1 N i x i. We consider the networks with N = 1 4, B = [, 3], nd k = 4 sme s the ones in the min text. The only difference is tht the coupling strength is chnged to D = 1.2 wheres we took D = in the min text. In Fig. 9 we consider deteriorting (meliorting) scenrio where the sttionry vlue of () is obtined by incresing the stress with n increment (decrement) δ = 1 from =.8 ( =.8) to =.8 ( =.8). The informtion on the other numeric experiments is given in the cption of - - () P= (d) (b) P= (e) (c) P= (f) FIG. 1: (color online). Sttionry stte distributions of the networks of the miniml systems with D = 1.2. The Wtts- Strogtz networks with () P =, (b) P = 5, (c) P =.1, nd (d) P =.25. (e) The Erdős-Réyni network. (f) The Brbási-Albert network. For ech pnel, we generte sttionry sttes with the stress nd the initil condition rndomly selected from uniform distributions [.4,.6] nd [ 1.6, 1.6], respectively. ()-(b) The upper ttrctors exist fter the thresholds but the density of the sttionry sttes hving < <.2 decreses s P increses in the regulr networks. (c)-(f) The upper ttrctors dispper t the thresholds in the disordered networks. ech figure. Figure 1 represents the sttionry stte distributions of the networks of the miniml systems with D = 1.2. Figure 11 depicts the distributions of stte P( ) P( ) P( ) P( ) 1 () =.45, Globl (b) =.46, Globl (c) =.47, Globl 1-1 P= (d) =.45, Locl (e) =.46, Locl (f) =.47, Locl (g) =.47, Globl (h) =.48, Globl (i) =.49, Globl (j) =.47, Locl (k) =.48, Locl (l) =.49, Locl FIG. 11: (color online). Distributions of stte chnges = A B in the networks of the miniml systems with D = 1.2 fter perturbtions. We obtin sttionry sttes B of network for given stress nd the initil condition rndomly chosen from uniform distribution [, 1.6] nd ssigned eqully to ll the nodes. If B > (i.e., the network is ctive), we give smll shock S to the network nd obtin the sttionry stte A fter the perturbtion. ()-(c) Globl perturbtions (S G =.1) were given to ll the nodes in the Wtts-Strogtz (WS) networks with P = nd P =.25. (d)- (f) Locl perturbtions (S L = ) were given to rndomly selected 1% of the nodes in the WS networks with P = nd P =.25. (g)-(i) Globl perturbtions (S G =.1) were given to ll the nodes in the Erdős-Réyni () network nd the Brbási-Albert () network. (j)-(l) Locl perturbtions (S L = ) were given to rndomly selected 1% of the nodes in the nd the networks

8 8 () r + P= LR (b) P= P=5 P=.1 - () P= P= - (b) (c) P= FIG. 12: (color online). Suppressed collpse of individul nodes in the networks of the miniml systems with D = 1.2. () The rtio of ctive nodes r + s function of the stress shows the synchronized collpses of the nodes in the Wtts- Strogtz (WS) network with P =.25, the Erdős-Réyni network, nd the Brbási-Albert network but the grdul collpses of the nodes in the WS network with P = nd the loclly rewired (LR) networks. To generte the LR network, we divide WS network into four equivlent sub-lttice before link rewiring nd rewire only the links between the nodes in one of the sub-lttice. Thus the sme number of links is rewired in the LR network s in the WS network with P =.25. Ech line is verged over 5 reliztions. (b) The sttionry stte distribution of the LR network shows smll gp indicting the collpse of the rewired region in the upper ttrctor. We generte sttionry sttes with the stress nd the initil condition rndomly selected from uniform distributions [.4,.6] nd [ 1.6, 1.6], respectively (d) (e) (f) FIG. 13: (color online). Stte digrms () of the networks of the mixed miniml systems with B = [, ]. () The Wtts-Strogtz (WS) networks with P.1 respond grdully to the incresing stress. (b) The WS networks with P.25, the Erdős-Réyni (), nd the Brbási-Albert () respond bruptly to the incresing stress. Incresing nd then decresing the stress, we observed hysteresis ptterns in (c) the WS network with P =, (d) the WS network with P =.25, (e) the network, nd (f) the network. Ech line with symbols is verged over 5 reliztions. The dotted lines in (c)-(f) show tht the hysteresis cn be mitigted in the WS network with P = but not in the WS network with P =.25, the network, nd the network. Ech dotted line shows single reliztion obtined by incresing nd then decresing the stress when it reches = 4, 5, 6, 7, 8, 9,.1, or.11. chnges = A B in the networks of the miniml systems with D = 1.2 fter perturbtions. Figure 12 shows the suppressed collpse of individul nodes in the networks of the miniml systems with D = 1.2. Appendix C: Networks of the mixed miniml systems We consider networks of the sme miniml systems defined in the min text. Hence the stte of node i in network of N miniml systems is evolved by dt = + b ix i x i 3 + D k i A ij (x j x i ), (C1) where x i is the stte of node i, is the environmentl stress, b i is the growth rte of node i drwn from uniform distribution B, D is the overll coupling strength mong the nodes, k i is the degree of node i, nd A ij is the element of the djcency mtrix of the network. The network s stte is defined s = 1 N i x i. We consider the networks with N = 1 4, D =, nd k = 4 sme s the ones in the min text. The only difference is the uniform distribution B for the growth rte of the nodes is chnged to B = [, ] wheres B = [, 3] in the min text. This mens tht ll nodes respond bruptly when isolted in the min text. On the other hnd, here, the hlf of the nodes responds grdully nd the other hlf responds bruptly when the nodes re isolted. In Fig. 13 we consider deteriorting (meliorting) scenrio where the sttionry vlue of () is obtined by incresing the stress with n increment (decrement) δ = 5 from =.2 ( =.2) to =.2 ( =.2). The informtion on the other numeric experiments is given in the cption of ech figure. Figure 14 represents the sttionry stte distributions of the networks of the mixed miniml systems with B = [, ]. Figure 15 depicts the distributions of stte chnges = A B in the networks of the mixed miniml systems with B = [, ] () P= (d) (b) P= (e) (c) P= (f) FIG. 14: (color online). Sttionry stte distributions of the networks of the mixed miniml systems with B = [, ]. The Wtts-Strogtz networks with () P =, (b) P = 5, (c) P =.1, nd (d) P =.25. (e) The Erdős-Réyni network. (f) The Brbási-Albert network. For ech pnel, we generte sttionry sttes with the stress nd the initil condition rndomly selected from uniform distributions [.2,.2] nd [ 1.6, 1.6], respectively. ()-(c) The upper ttrctors exist fter the thresholds but the density of the sttionry sttes hving.3 < <.3 decreses s P increses in the regulr networks. (d)-(f) The upper ttrctors dispper t the thresholds in the disordered networks

9 9 P( ) P( ) P( ) P( ) 1 () =65, Globl (b) =7, Globl (c) =75, Globl 1-1 P= (d) =65, Locl (e) =7, Locl (f) =75, Locl (g) =65, Globl (h) =7, Globl (i) =75, Globl (j) =65, Locl (k) =7, Locl (l) =75, Locl FIG. 15: (color online). Distributions of stte chnges = A B in the networks of the mixed miniml systems with B = [, ] fter perturbtions. We obtin sttionry sttes B of network for given stress nd the initil condition rndomly chosen from uniform distribution [, 1.6] nd ssigned eqully to ll the nodes. If B > (i.e., the network is ctive), then we give smll shock S to the network nd obtin the sttionry stte A fter the perturbtion. ()-(c) Globl perturbtions (S G =.1) were given to ll the nodes in the Wtts-Strogtz (WS) networks with P = nd P =.25. (d)-(f) Locl perturbtions (S L = ) were given to rndomly selected 1% of the nodes in the WS networks with P = nd P =.25. (g)-(i) Globl perturbtions (S G =.1) were given to ll the nodes in the Erdős-Réyni () network nd the Brbási-Albert () network. (j)- (l) Locl perturbtions (S L = ) were given to rndomly selected 1% of the nodes in the nd the networks. fter perturbtions. Figure 16 shows the suppressed collpse of individul nodes in the networks of the mixed miniml systems with B = [, ]. Appendix D: Networks of the grzing systems We consider the grzing system [51, 52], prdigmtic model for logisticlly growing vegettion under grzing stress, for the behvior of the nodes. The stte of node i in network of N grzing systems is evolved by dt = x i(1 x i ) x2 i b i x 2 i D k i A ij (x j x i ), (D1) where x i ( ) is the biomss density of node i, is the mximum consumption rte (i.e., grzing stress), b i is the crrying cpcity of node i rndomly drwn from uniform distribution B = [5., 1], nd D is the overll coupling strength (D =.3). The first term represents the logistic growth of the vegettion, the second term represents the grzing stress, nd the third term represents the dispersion of the biomss. The network s stte is defined s = 1 N i x i. We consider the networks () r + P= LR (b) FIG. 16: (color online). Suppressed collpse of individul nodes in the networks of the mixed miniml systems with B = [, ]. () The rtio of ctive nodes r + s function of the stress shows the synchronized collpses of the nodes in the Wtts-Strogtz (WS) network with P =.25, the Erdős-Réyni network, nd the Brbási-Albert network but the grdul collpses of the nodes in the WS network with P = nd the loclly rewired (LR) networks. To generte the LR network, we divide WS network into four equivlent sublttice before link rewiring nd rewire only the links between the nodes in one of the sublttice. Thus the sme number of links is rewired in the LR network s in the WS network with P =.25. Ech line is verged over 5 reliztions. (b) The sttionry stte distribution of the LR network shows smll gp indicting the collpse of the rewired region in the upper ttrctor. We generte sttionry sttes with the stress nd the initil condition rndomly selected from uniform distributions [5,.15] nd [ 1.6, 1.6], respectively. with N = 1 4 nd k = 4. The bifurction digrm of the grzing systems when they re isolted is depicted in Fig. 17. We consider the grzing system i s ctive (collpsed) for x i > (x i < ). In Fig. 18 we consider deteriorting (meliorting) scenrio where the sttionry vlue of () is obtined by incresing the stress with n increment (decrement) δ = 1 from = 1.6 ( = 2.3) to = 2.3 ( = 1.6). The informtion on the other numeric experiments is given in the cption of ech figure. Figure 19 represents the sttionry stte distributions of the networks of the grzing systems. Figure 2 depicts the distributions of stte chnges = A B in the networks of the grzing systems fter perturbtions. Figure 21 shows the suppressed collpse of individul nodes in the networks of the grzing systems. x b=5. b=7.5 b= FIG. 17: (color online). Bifurction digrm of grzing system whose dynmics is described by dx/dt = x(1 x/b) x 2 /(x 2 + 1). The red lines represent the stble sttes nd the dotted gry line represents the unstble sttes.

10 P= P=5 P=.1 () P= P= (b) (c) P= (d) (e) (f) FIG. 18: (color online). Stte digrms () of the networks of the grzing systems. () The Wtts-Strogtz (WS) networks with P.1 respond grdully to the incresing stress. (b) The WS networks with P.25, the Erdős- Réyni (), nd the Brbási-Albert () respond bruptly to the incresing stress. Incresing nd then decresing the stress, we observed hysteresis ptterns in (c) the WS network with P =, (d) the WS network with P =.25, (e) the network, nd (f) the network. Ech line with symbols is verged over 5 reliztions. The dotted lines in (c)-(f) show tht the hysteresis cn be substntilly mitigted in the WS network with P =. But only wek mitigtion is observed in the WS network with P =.25, the network, nd the network becuse some nodes with high growth rtes (e.g., nodes with b > 9) re still ctive fter the thresholds nd they cn provide strong subsidy to other nodes s shown in Fig. 17. Ech dotted line shows single reliztion obtined by incresing nd then decresing the stress when it reches = 1.96, 1.97, 1.98, 1.99,, 1, or 2. () P= (d) (b) P= (e) (c) P= (f) FIG. 19: (color online). Sttionry stte distributions of the networks of the grzing systems. The Wtts-Strogtz networks with () P =, (b) P = 5, (c) P =.1, nd (d) P =.25. (e) The Erdős-Réyni network. (f) The Brbási- Albert network. For ech pnel, we generte sttionry sttes with the stress nd the initil condition rndomly selected from uniform distributions [1.9, ] nd [.1, 5.1], respectively. ()-(c) The upper ttrctors exist fter the thresholds but the density of the sttionry sttes hving < < decreses s P increses in the regulr networks. (d)-(f) The upper ttrctors dispper t the thresholds in the disordered networks. Acknowledgments We thnk Hng-Hyun Jo, Seung Ki Bek, nd Seung- Woo Son for creful reding of the preliminry P( ) P( ) P( ) P( ) 1 () =1.97, Globl (b) =1.98, Globl (c) =1.99, Globl 1-1 P= (d) =1.97, Locl (e) =1.98, Locl (f) =1.99, Locl (g) =1.97, Globl (h) =1.98, Globl (i) =1.99, Globl (j) =1.97, Locl (k) =1.98, Locl (l) =1.99, Locl FIG. 2: (color online). Distributions of stte chnges = A B in the networks of the grzing systems fter perturbtions. We obtin sttionry sttes B of network for given stress nd the initil condition rndomly chosen from uniform distribution [4., 6.] nd ssigned eqully to ll the nodes. If B > (i.e., the network is ctive), then we give smll shock S to the network nd obtin the sttionry stte A fter the perturbtion. ()-(c) Globl perturbtions (S G =.2) were given to ll the nodes in the Wtts-Strogtz (WS) networks with P = nd P =.25. (d)-(f) Locl perturbtions were given to rndomly selected 1% of the nodes in the WS networks with P = nd P =.25. (g)-(i) Globl perturbtions (S G =.2) were given to ll the nodes in the Erdős-Réyni () network nd the Brbási-Albert () network. (j)-(l) Locl perturbtions were given to rndomly selected 1% of the nodes in the nd the networks. For the locl perturbtions, we rndomly selected 1% of the nodes nd set their stte to. () r + P= LR (b) FIG. 21: (color online). Suppressed collpse of individul nodes in the networks of the grzing systems. Note tht we consider the grzing system i s ctive when x i >. () The rtio of ctive nodes r + s function of the stress shows the synchronized collpses of the nodes in the Wtts-Strogtz (WS) network with P =.25, the Erdős-Réyni network, nd the Brbási-Albert network but the grdul collpses of the nodes in the WS network with P = nd the loclly rewired (LR) networks. To generte the LR network, we divide WS network into four equivlent sublttice before link rewiring nd rewire only the links between the nodes in one of the sublttice. Thus the sme number of links is rewired in the LR network s in the WS network with P =.25. Ech line is verged over 5 reliztions. (b) The sttionry stte distribution of the LR network shows smll gp indicting the collpse of the rewired region in the upper ttrctor. We generte sttionry sttes with the stress nd the initil condition rndomly selected from uniform distributions [1.9, ] nd [.1, 5.1], respectively.

11 11 mnuscript nd for helpful comments. Y.-H.E. cknowledges support from the University of Strthclyde, funding from Ministerio de Economí y Competividd (Spin) through project FIS C3-3-P, nd support from the Universidd Crlos III de Mdrid, the Europen Union s Seventh Frmework Progrmme for reserch, technologicl development nd demonstrtion under grnt greement (no. 6371), el Ministerio de Economí y Competividd (COFUND ) nd Bnco Sntnder. [1] A.-L. Brbási, Network science (Cmbridge University Press, Cmbridge, UK, 216). [2] J. Go, B. Brzel, nd A.-L. Brbási, Nture 53, 37 (216). [3] R. Albert, H. Jeong, nd A.-L. Brbási, Nture 46, 378 (2). [4] R. Cohen, K. Erez, D. ben Avrhm, nd S. Hvlin, Phys. Rev. Lett. 85, 4626 (2). [5] D. S. Cllwy, M. E. J. Newmn, S. H. Strogtz, nd D. J. Wtts, Phys. Rev. Lett. 85, 5468 (2). [6] S. V. Buldyrev, R. Prshni, G. Pul, H. E. Stnley, nd S. Hvlin, Nture 464, 125 (21). [7] R. Prshni, S. V. Buldyrev, nd S. Hvlin, Phys. Rev. Lett. 15, 4871 (21). [8] J. Go, S. V. Buldyrev, S. Hvlin, nd H. E. Stnley, Phys. Rev. Lett. 17, (211). [9] J. Go, S. V. Buldyrev, H. E. Stnley, nd S. Hvlin, Nt. Phys. 8, 4 (211). [1] A. Bshn, Y. Berezin, S. V. Buldyrev, nd S. Hvlin, Nt. Phys. 9, 667 (213). [11] R. Cohen, K. Erez, D. ben Avrhm, nd S. Hvlin, Phys. Rev. Lett. 86, 3682 (21). [12] A. Vázquez nd Y. Moreno, Phys. Rev. E 67, 1511(R) (23). [13] A.. C. N. Vlente, A. Srkr, nd H. A. Stone, Phys. Rev. Lett. 92, (24). [14] L. K. Gllos, R. Cohen, P. Argyrkis, A. Bunde, nd S. Hvlin, Phys. Rev. Lett. 94, (25). [15] S. Lozno nd A. Arens, J. Artif. Soc. Soc. Simul. 1, 8 (27). [16] S. Svedr, R. P. Rohr, V. Dkos, nd J. Bscompte, Nt. Commun. 4, 235 (213). [17] V. Dkos nd J. Bscompte, Proc. Ntl. Acd. Sci. USA 111, (214). [18] J. J. Lever, E. H. vn Nes, M. Scheffer, nd J. Bscompte, Ecol. Lett. 17, 35 (214). [19] A. D. Brnosky, E. A. Hdly, J. Bscompte, E. L. Berlow, J. H. Brown, M. Fortelius, W. M. Getz, J. Hrte, A. Hstings, P. A. Mrquet, et l., Nture 486, 52 (212). [2] M. Scheffer, S. Crpenter, T. M. Lenton, J. Bscompte, W. Brock, V. Dkos, J. vn de Koppel, I. A. vn de Leemput, S. A. Levin, E. H. vn Nes, et l., Science 338, 344 (212). [21] M. Scheffer, J. Bscompte, W. A. Brock, V. Brovkin, S. Crpenter, V. Dkos, H. Held, E. H. vn Nes, M. Rietkerk, nd G. Sugihr, Nture 461, 53 (29). [22] M. Scheffer, S. Crpenter, J. A. Foley, nd B. Wlker, Nture 413, 591 (21). [23] S. R. Crpenter, Proc. Ntl. Acd. Sci. USA 12, 12 (25). [24] S. Kéfi, M. Rietkerk, C. L. Aldos, Y. Pueyo, V. P. Ppnstsis, A. ElAich, nd P. C. de Ruiter, Nture 449, 213 (27). [25] M. Scheffer, Criticl trnsitions in nture nd society (Princeton Studiess in Complexity, Princeton University Press, Princeton, NJ, USA, 29). [26] J. M. Drke nd B. D. Griffen, Nture 467, 456 (21). [27] S. R. Crpenter, J. J. Cole, M. L. Pce, R. Btt, B. A. Brock, T. Cline, J. Colosso, J. R. Hodgson, J. F. Kitchell, D. A. Seekell, et l., Science 332, 179 (211). [28] L. Di, D. Vorselen, K. S. Korolev, nd J. Gore, Science 336, 1175 (212). [29] A. J. Vert, E. J. Fssen, V. Dkos, E. H. vn Nes, M. Lürling, nd M. Scheffer, Nture 481, 357 (212). [3] T. M. Lenton, H. Held, E. Kriegler, J. W. Hll, W. Lucht, S. Rhmstorf, nd H. J. Schellnhuber, Proc. Ntl. Acd. Sci. USA 15, 1786 (28). [31] R. M. My, S. A. Levin, nd G. Sugihr, Nture 451, 893 (28). [32] M. Rietkerk, S. C. Dekker, P. C. de Ruiter, nd J. vn de Koppel, Science 35, 1926 (24). [33] V. Guttl nd C. Jyprksh, Theor. Ecol. 2, 3 (29). [34] A. Fernández nd H. Fort, J. Stt. Mech. P914 (29). [35] V. Dkos, E. H. vn Nes, R. Donngelo, H. Fort, nd M. Scheffer, Theor. Ecol. 3, 163 (21). [36] L. Di, K. S. Korolev, nd J. Gore, Nture 496, 355 (213). [37] H. Weissmnn nd N. M. Shnerb, J. Theor. Biol. 397, 128 (216). [38] E. H. vn Nes nd M. Scheffer, Ecology 86, 1797 (25). [39] P. V. Mrtín, J. A. Bonchel, S. A. Levin, nd M. A. Muñoz, Proc. Ntl. Acd. Sci. USA 112, E1828 (215). [4] H. Weissmnn, N. M. Shnerb, nd D. A. Kessler, preprint on webpge t rxiv.org/bs/ [41] E. C. Zeemn, Sci. Am. 234, 65 (1976). [42] S. H. Strogtz, Nonliner dynmics nd chos: with pplictions to physics, biology, chemistry, nd engineering (Addison Wesley, Reding, MA, USA, 1994). [43] E. H. vn Nes, B. M. S. Arni, A. Stll, B. vn der Bolt, B. M. Flores, S. Bthiny, nd M. Scheffer, Trends Eco. Evol. 31, 92 (216). [44] C. Kuehn, Phys. D 24, 12 (211). [45] C. D. Brumitt, G. Brnett, nd R. M. D Souz, J. R. Soc. Interfce 12, (215). [46] D. J. Wtts nd S. H. Strogtz, Nture 393, 44 (1998). [47] B. Bollobás, Rndom Grphs (Cmbridge University Press, Cmbridge, UK, 28). [48] A.-L. Brbási nd R. Albert, Science 286, 59 (1999). [49] R. M. D Souz nd J. Ngler, Nt. Phys. 11, 531 (215). [5]. Zhng, S. Boccletti, S. Gun, nd Z. Liu, Phys. Rev. Lett. 114, 3871 (215). [51] I. Noy-Meir, J. Ecol. 63, 459 (1975). [52] R. M. My, Nture 269, 471 (1977).