Catalog footprint

What is connected

133works
34topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

133 published item(s)

preprint2026arXiv

Long-term prediction of ENSO with physics-guided Deep Echo State Networks

The El Niño-Southern Oscillation (ENSO) is a dominant mode of interannual climate variability, yet the mechanisms limiting its long-lead predictability remain unclear. Here we develop a physics-guided Deep Echo State Network (DESN) that operates on physically interpretable climate modes selected from the extended recharge oscillator (XRO) framework. DESN achieves skillful Niño3.4 predictions up to 16-20 months ahead with minimal computational cost. Mechanistic experiments show that extended predictability arises from nonlinear coupling between warm water volume and inter-basin climate modes. Error-growth analysis further indicates a finite ENSO predictability horizon of approximately 30 months. These results demonstrate that physics-guided reservoir computing provides an efficient and interpretable framework for diagnosing and predicting ENSO at long lead times.

preprint2022arXiv

Delayed Impact of Interdisciplinary Research

Interdisciplinary research increasingly fuels innovation, and is considered to be a key to tomorrow breakthrough. Yet little is known about whether interdisciplinary research manifests delayed impact. Here, we use the time to reach the citation peak to quantify the highest impact time and citation dynamics, and examine its relationship with interdisciplinarity. Using large scale publication datasets, our results suggest that interdisciplinary papers show significant delayed impact both microscopically per paper and macroscopically collectively, as it takes longer time for interdisciplinary papers to reach their citation peak. Furthermore, we study the underlying forces of such delayed impact, finding that the effect goes beyond the Matthew effect (i.e., the rich-get-richer effect). Finally, we find that team size and content conventionality only partly account for this effect. Overall, our results suggest that governments, research administrators, funding agencies should be aware of this general feature of interdisciplinary science, which may have broad policy implications.

preprint2022arXiv

Epidemics on evolving networks with varying degrees

Epidemics on complex networks is a widely investigated topic in the last few years, mainly due to the last pandemic events. Usually, real contact networks are dynamic, hence much effort has been invested in studying epidemics on evolving networks. Here we propose and study a model for evolving networks based on varying degrees, where at each time step a node might get, with probability $r$, a new degree and new neighbors according to a given degree distribution, instead of its former neighbors. We find analytically, using the generating functions framework, the epidemic threshold and the probability for a macroscopic spread of disease depending on the rewiring rate $r$. Our analytical results are supported by numerical simulations. We find surprisingly that the impact of the rewiring rate $r$ has qualitative different trends for networks having different degree distributions. That is, in some structures, such as random regular networks the dynamics enhances the epidemic spreading while in others such as scale free the dynamics reduces the spreading. In addition, for scale-free networks, we reveal that fast dynamics of the network, $r=1$, changes the epidemic threshold to nonzero rather than zero found for $r<1$, which is similar to the known case of $r=0$, i.e., a static network. Finally, we find the epidemic threshold also for a general distribution of the recovery time.

preprint2022arXiv

Fractal fluctuations at mixed-order transitions in interdependent networks

We study the geometrical features of the order parameter's fluctuations near the critical point of mixed-order phase transitions in randomly interdependent spatial networks. In contrast to continuous transitions, where the structure of the order parameter at criticality is fractal, in mixed-order transitions the structure of the order parameter is known to be compact. Remarkably, we find that although being compact, the fluctuations of the order parameter close to mixed-order transitions are fractal up to a well-defined correlation length $ξ'$, which diverges when approaching the critical threshold. We characterize the self-similar nature of these critical fluctuations through their fractal dimension, $d_f'=3d/4$, and correlation length exponent, $ν'=2/d$, where $d$ is the dimension of the system. By means of percolation and magnetization, we demonstrate that $d_f'$ and $ν'$ are independent on the symmetry of the underlying process for any $d$ of the underlying networks.

preprint2022arXiv

Impactful scientists have higher tendency to involve collaborators in new topics

In scientific research, collaboration is one of the most effective ways to take advantage of new ideas, skills, resources, and for performing interdisciplinary research. Although collaboration networks have been intensively studied, the question of how individual scientists choose collaborators to study a new research topic remains almost unexplored. Here, we investigate the statistics and mechanisms of collaborations of individual scientists along their careers, revealing that, in general, collaborators are involved in significantly fewer topics than expected from controlled surrogate. In particular, we find that highly productive scientists tend to have higher fraction of single-topic collaborators, while highly cited, i.e., impactful, scientists have higher fraction of multi-topic collaborators. We also suggest a plausible mechanism for this distinction. Moreover, we investigate the cases where scientists involve existing collaborators into a new topic. We find that compared to productive scientists, impactful scientists show strong preference of collaboration with high impact scientists on a new topic. Finally, we validate our findings by investigating active scientists in different years and across different disciplines.

preprint2022arXiv

Percolation on spatial anisotropic networks

Many realistic systems such as infrastructures are characterized by spatial structure and anisotropic alignment. Here we propose and study a model for dealing with such characteristics by introducing a parameter that controls the strength of the anisotropy in the spatial network. This parameter is added to an existing isotropic model used to describe networks under spatial constraints, thus generalizing the spatial model to take into account both spatial and anisotropic features. We study the resilience of such networks by using a percolation process and find that anisotropy has a negative impact on a network's robustness. In addition, our results suggest that the anisotropy in this model does not affect the critical exponent of the correlation length, $ν$, which remains the same as the known $ν$ in 2D isotropic lattices.

preprint2022arXiv

Reviving a failed network through microscopic interventions

From mass extinction to cell death, complex networked systems often exhibit abrupt dynamic transitions between desirable and undesirable states. Such transitions are often caused by topological perturbations, such as node or link removal, or decreasing link strengths. The problem is that reversing the topological damage, namely retrieving the lost nodes or links, or reinforcing the weakened interactions, does not guarantee the spontaneous recovery to the desired functional state. Indeed, many of the relevant systems exhibit a hysteresis phenomenon, remaining in the dysfunctional state, despite reconstructing their damaged topology. To address this challenge, we develop a two-step recovery scheme: first - topological reconstruction to the point where the system can be revived, then dynamic interventions, to reignite the system's lost functionality. Applying this method to a range of nonlinear network dynamics, we identify the recoverable phase of a complex system, a state in which the system can be reignited by microscopic interventions, for instance, controlling just a single node. Mapping the boundaries of this dynamical phase, we obtain guidelines for our two-step recovery.

preprint2022arXiv

Sustaining a network by controlling a fraction of nodes

Multi-stability is a widely observed phenomenon in real complex networked systems, such as technological infrastructures, ecological systems, gene regulation, transportation and more. When a system functions normally but there exists also a potential state with abnormal low activity, although the system is at equilibrium it might make a transition into the low activity undesired state due to external disturbances and perturbations. Thus, such a system can be regarded as unsustainable, due to the danger of falling into the potential inactive state. Here we explore, analytically and by simulations, how supporting the activity of a fraction $ρ$ of nodes can turn an unsustainable system to be sustainable by eliminating the inactive potential stable state. We thus unveil a new sustainability phase diagram in the presence of a fraction of controlled nodes $ρ$. This phase diagram could provide guidelines to sustain a network by external intervention and/or by strengthening the connectivity of the network.

preprint2022arXiv

Teleconnections among Tipping Elements in the Earth System

Tipping elements of the Earth system may shift abruptly and irreversibly from one state to another at tipping points, resulting in a growing threat to our society. Yet, it is not fully clear how to assess and quantify the influence of a tipping element and how to explore the teleconnections between different tipping elements. To fill this knowledge gap, we propose a climate network approach to quantitatively analyze the global impacts of a prominent tipping element, the Amazon Rainforest Area (ARA). We find that regions, such as, the Tibetan Plateau (TP) and West Antarctic ice sheet, are characterized by higher network weighted links and exhibit strong correlations with the ARA. We then identify a teleconnection propagation path between the ARA and the TP. This path is robust under climate change as simulated by various climate models of CMIP5 and CMIP6. In addition, we detect early warning signals for critical transition in the snow cover extent on the Tibetan Plateau by applying critical slowing down indicators, lag-1 autocorrelation and detrended fluctuation analysis. We find that the snow cover of the TP has been losing stability since 2008, revealing that the TP is operating like a tipping element and approaching a potential tipping point. We further uncover that various climate extremes between the ARA and the TP are significantly synchronized under climate change. Our framework provides new insights into how tipping elements are linked to each other and into the potential predictability of cascading tipping dynamics.

preprint2022arXiv

Topological-temporal properties of evolving networks

Many real-world complex systems including human interactions can be represented by temporal (or evolving) networks, where links activate or deactivate over time. Characterizing temporal networks is crucial to compare such systems and to study the dynamical processes unfolding on them. A systematic method to characterize simultaneously the temporal and topological relations of active links (also called contacts or events), in order to compare different real-world networks and to detect their common patterns or differences is still missing. In this paper, we propose a method to characterize to what extent contacts that happen close in time occur also close in topology. Specifically, we study the interrelation between temporal and topological properties of contacts from three perspectives: (1) the autocorrelation of the time series recording the total number of contacts happened at each time step in a network; (2) the interplay between the topological distance and interevent time of two contacts; (3) the temporal correlation of contacts within local neighborhoods beyond a node pair. By applying our method on 13 real-world temporal networks, we found that temporal-topological correlation of contacts is more evident in virtual contact networks than in physical contact ones. This could be due to the lower cost and easier access of online communications than physical interactions, allowing and possibly facilitating social contagion, i.e., interactions of one individual may influence the activity of its neighbors. We also identify different patterns between virtual and physical networks and among physical contact networks at, e.g., school and workplace, in the formation of correlation in local neighborhoods. Detected patterns and differences may further inspire the development of more realistic temporal network models, that could reproduce jointly temporal and topological properties of contacts.

preprint2021arXiv

Optimization of robustness based on reinforced nodes in a modular network

Many systems such as critical infrastructure exhibit a modular structure with many links within the modules and few links between them. One approach to increase the robustness of these systems is to reinforce a fraction of the nodes in each module, so that the reinforced nodes provide additional needed sources for themselves as well as for their nearby neighborhood. Since reinforcing a node can be an expensive task, the efficiency of the decentralization process by reinforced nodes is vital. In our study we analyze a new model which combines both above mentioned features of real complex systems - modularity and reinforced nodes. Using tools from percolation theory, we derived an analytical solution for any partition of reinforced nodes; between nodes which have links that connect them to other modules ("inter-nodes") and nodes which have connections only within their modules ("intra-nodes"). Among our results, we find that near the critical percolation point ($p\approx p_c$) the robustness is greatly affected by the distribution. In particular, we find a partition of reinforced nodes which yields an optimal robustness and we show that the optimal partition remains constant for high average degrees.

preprint2020arXiv

Cascading Failures in Complex Networks

Cascading failure is a potentially devastating process that spreads on real-world complex networks and can impact the integrity of wide-ranging infrastructures, natural systems, and societal cohesiveness. One of the essential features that create complex network vulnerability to failure propagation is the dependency among their components, exposing entire systems to significant risks from destabilizing hazards such as human attacks, natural disasters or internal breakdowns. Developing realistic models for cascading failures as well as strategies to halt and mitigate the failure propagation can point to new approaches to restoring and strengthening real-world networks. In this review, we summarize recent progress on models developed based on physics and complex network science to understand the mechanisms, dynamics and overall impact of cascading failures. We present models for cascading failures in single networks and interdependent networks and explain how different dynamic propagation mechanisms can lead to an abrupt collapse and a rich dynamic behavior. Finally, we close the review with novel emerging strategies for containing cascades of failures and discuss open questions that remain to be addressed.

preprint2020arXiv

Discrimination universally determines reconstruction of multiplex networks

Network reconstruction is fundamental to understanding the dynamical behaviors of the networked systems. Many systems, modeled by multiplex networks with various types of interactions, display an entirely different dynamical behavior compared to the corresponding aggregated network. In many cases, unfortunately, only the aggregated topology and partial observations of the network layers are available, raising an urgent demand for reconstructing multiplex networks. We fill this gap by developing a mathematical and computational tool based on the Expectation-Maximization framework to reconstruct multiplex layer structures. The reconstruction accuracy depends on the various factors, such as partial observation and network characteristics, limiting our ability to predict and allocate observations. Surprisingly, by using a mean-field approximation, we discovered that a discrimination indicator that integrates all these factors universally determines the accuracy of reconstruction. This discovery enables us to design the optimal strategies to allocate the fixed budget for deriving the partial observations, promoting the optimal reconstruction of multiplex networks. To further evaluate the performance of our method, we predict beside structure also dynamical behaviors on the multiplex networks, including percolation, random walk, and spreading processes. Finally, applying our method on empirical multiplex networks drawn from biological, transportation, and social domains, corroborate the theoretical analysis.

preprint2020arXiv

Distance Distribution in Extreme Modular Networks

Modularity is a key organizing principle in real-world large-scale complex networks. Many real-world networks exhibit modular structures such as transportation infrastructures, communication networks and social media. Having the knowledge of the shortest paths length distribution (DSPL) between random pairs of nodes in such networks is important for understanding many processes, including diffusion or flow. Here, we provide analytical methods which are in good agreement with simulations on large scale networks with an extreme modular structure. By extreme modular, we mean that two modules or communities may be connected by maximum one link. As a result of the modular structure of the network, we obtain a distribution showing many peaks that represent the number of modules a typical shortest path is passing through. We present theory and results for the case where inter-links are weighted, as well as cases in which the inter-links are spread randomly across nodes in the community or limited to a specific set of nodes.

preprint2020arXiv

Efficient network immunization under limited knowledge

Targeted immunization or attacks of large-scale networks has attracted significant attention by the scientific community. However, in real-world scenarios, knowledge and observations of the network may be limited thereby precluding a full assessment of the optimal nodes to immunize (or remove) in order to avoid epidemic spreading such as that of current COVID-19 epidemic. Here, we study a novel immunization strategy where only $n$ nodes are observed at a time and the most central between these $n$ nodes is immunized (or attacked). This process is continued repeatedly until $1-p$ fraction of nodes are immunized (or attacked). We develop an analytical framework for this approach and determine the critical percolation threshold $p_c$ and the size of the giant component $P_{\infty}$ for networks with arbitrary degree distributions $P(k)$. In the limit of $n\to\infty$ we recover prior work on targeted attack, whereas for $n=1$ we recover the known case of random failure. Between these two extremes, we observe that as $n$ increases, $p_c$ increases quickly towards its optimal value under targeted immunization (attack) with complete information. In particular, we find a new scaling relationship between $|p_c(\infty)-p_c(n)|$ and $n$ as $|p_c(\infty)-p_c(n)|\sim n^{-1}\exp(-αn)$. For Scale-free (SF) networks, where $P(k)\sim k^{-γ}, 2<γ<3$, we find that $p_c$ has a transition from zero to non-zero when $n$ increases from $n=1$ to order of $\log N$ ($N$ is the size of network). Thus, for SF networks, knowledge of order of $\log N$ nodes and immunizing them can reduce dramatically an epidemics.

preprint2020arXiv

Epidemic spreading and control strategies in spatial modular network

Epidemic spread on networks is one of the most studied dynamics in network science and has important implications in real epidemic scenarios. Nonetheless, the dynamics of real epidemics and how it is affected by the underline structure of the infection channels are still not fully understood. Here we apply the SIR model and study analytically and numerically the epidemic spread on a recently developed spatial modular model imitating the structure of cities in a country. The model assumes that inside a city the infection channels connect many different locations, while the infection channels between cities are less and usually directly connect only a few nearest neighbor cities in a two-dimensional plane. We find that the model experience two epidemic transitions. The first lower threshold represents a local epidemic spread within a city but not to the entire country and the second higher threshold represents a global epidemic in the entire country. Based on our analytical solution we proposed several control strategies and how to optimize them. We also show that while control strategies can successfully control the disease, early actions are essentials to prevent the disease global spread.

preprint2020arXiv

Geometric characterization of SARS-CoV-2 pandemic events

While the SARS-CoV-2 keeps spreading world-wide, comparing its evolution across different nations is a timely challenge of both theoretical and practical importance. The large variety of dissimilar and country-dependent epidemiological factors, in fact, makes extremely difficult to understand their influence on the epidemic trends within a unique and coherent framework. We present a geometric framework to characterize, in an integrated and low-dimensional fashion, the epidemic plume-like trajectories traced by the infection rate, $I$, and the fatality rate, $D$, in the $(I,D)$ plane. Our analysis enables the definition of an epidemiometric system based on three geometric observables rating the SARS-CoV-2 pandemic events via scales analogous to those for the magnitude and the intensity of seismic events. Being exquisitely geometric, our framework can be applied to classify other epidemic data and secondary waves, raising the possibility of designing epidemic alerts or early warning systems to enhance public and governmental responses to a rapidly emerging outbreak.

preprint2020arXiv

Globalization emergence in the European Patent Office (EPO) patent network

We studied the evolution of the European Patent Office (EPO) patents applicants' collaborations network, within a 35 years span of data (1978-2013). Focusing on the Giant Component (GC) formation process over many time-windows, distributed throughout the data timeline, we found that the features governing this phenomenon are indicative of emerging globalization in the applicants' collaborations. The timeline appears to be divided into three regimes, corresponding to three states of the network's evolution. In the early years state, the GC takes long to form and the instant of its creation is easily pinpointed, while it features geographically segregated groups of applicants with technologically similar activities. In contrast, in the late years state, the GC forms quickly, the exact point of its creation is harder to spot, the applicants' activities are more disparate technologically, while their inter-regional collaborations are significantly increased. The middle years are an intermediate state between the two extreme of early and late years. Moreover, we concluded that the critical patents, which induce the GC's formation, are typically introduced by large-sized applicants and also that top patent-producing applicants are likely to submit critical patents, albeit at a lower rate than their overall patent submission. Lastly, we uncovered the crucial role that Japan plays in the network's coherence, through its prominent participation in the GC and the critical patents.

preprint2020arXiv

Improved Earthquake Forecasting Model Based on Long-term Memory in Earthquake

A prominent feature of earthquakes is their empirical laws including memory (clustering) in time and space. Several earthquake forecasting models, like the EpidemicType Aftershock Sequence (ETAS) model, were developed based on earthquake empirical laws. Yet, a recent study showed that the ETAS model fails in reproducing significant long-term memory characteristics found in real earthquake catalogs. Here we modify and generalize the ETAS model to include short- and long-term triggering mechanisms, to account for the short- and long-time memory (exponents) recently discovered in the data. Our generalized ETAS model reproduces accurately the short- and long-term/distance memory observed in the Italian and South California earthquake catalogs. The revised ETAS model is also found to significantly improve earthquake forecasting.

preprint2020arXiv

Interdependent transport via percolation backbones in spatial networks

The functionality of nodes in a network is often described by the structural feature of belonging to the giant component. However, when dealing with problems like transport, a more appropriate functionality criterion is for a node to belong to the network's backbone, where the flow of information and of other physical quantities (such as current) occurs. Here we study percolation in a model of interdependent resistor networks and show the effect of spatiality on their coupled functioning. We do this on a realistic model of spatial networks, featuring a Poisson distribution of link-lengths. We find that interdependent resistor networks are significantly more vulnerable than their percolation-based counterparts, featuring first-order phase transitions at link-lengths where the mutual giant component still emerges continuously. We explain this apparent contradiction by tracing the origin of the increased vulnerability of interdependent transport to the crucial role played by the dandling ends. Moreover, we interpret these differences by considering an heterogeneous $k$-core percolation process which enables to define a one-parameter family of functionality criteria whose constraints become more and more stringent. Our results highlight the importance that different definitions of nodes functionality have on the collective properties of coupled processes, and provide better understanding of the problem of interdependent transport in many real-world networks.

preprint2020arXiv

Network Geometry

Real networks are finite metric spaces. Yet the geometry induced by shortest path distances in a network is definitely not its only geometry. Other forms of network geometry are the geometry of latent spaces underlying many networks, and the effective geometry induced by dynamical processes in networks. These three approaches to network geometry are all intimately related, and all three of them have been found to be exceptionally efficient in discovering fractality, scale-invariance, self-similarity, and other forms of fundamental symmetries in networks. Network geometry is also of great utility in a variety of practical applications, ranging from the understanding how the brain works, to routing in the Internet. Here, we review the most important theoretical and practical developments dealing with these approaches to network geometry in the last two decades, and offer perspectives on future research directions and challenges in this novel frontier in the study of complexity.

preprint2020arXiv

Percolation framework reveals limits of privacy in Conspiracy, Dark Web, and Blockchain networks

We consider the privacy of interactions between individuals in a network. For many networks, while nodes are anonymous to outside observers, the existence of a link between individuals implies the possibility of one node revealing identifying information about its neighbor. Moreover, while the identities of the accounts are likely hidden to an observer, the network of interaction between two anonymous accounts is often available. For example, in blockchain cryptocurrencies, transactions between two anonymous accounts are published openly. Here we consider what happens if one (or more) parties in such a network are deanonymized by an outside identity. These compromised individuals could leak information about others with whom they interacted, which could then cascade to more and more nodes' information being revealed. We use a percolation framework to analyze the scenario outlined above and show for different likelihoods of individuals possessing information on their counter-parties, the fraction of accounts that can be identified and the idealized minimum number of steps from a deanonymized node to an anonymous node (a measure of the effort required to deanonymize that individual). We further develop a greedy algorithm to estimate the \emph{actual} number of steps that will be needed to identify a particular node based on the noisy information available to the attacker. We apply our framework to three real-world networks: (1) a blockchain transaction network, (2) a network of interactions on the dark web, and (3) a political conspiracy network. We find that in all three networks, beginning from one compromised individual, it is possible to deanonymize a significant fraction of the network ($>50$%) within less than 5 steps. Overall these results provide guidelines for investigators seeking to identify actors in anonymous networks, as well as for users seeking to maintain their privacy.

preprint2020arXiv

Predictability of real temporal networks

Links in most real networks often change over time. Such temporality of links encodes the ordering and causality of interactions between nodes and has a profound effect on network dynamics and function. Empirical evidences have shown that the temporal nature of links in many real-world networks is not random. Nonetheless, it is challenging to predict temporal link patterns while considering the entanglement between topological and temporal link patterns. Here we propose an entropy-rate based framework, based on combined topological-temporal regularities, for quantifying the predictability of any temporal network. We apply our framework on various model networks, demonstrating that it indeed captures the intrinsic topological-temporal regularities whereas previous methods considered only temporal aspects. We also apply our framework on 18 real networks of different types and determine their predictability. Interestingly, we find that for most real temporal networks, despite the greater complexity of predictability brought by the increase in dimension the combined topological-temporal predictability is higher than the temporal predictability. Our results demonstrate the necessity of incorporating both temporal and topological aspects of networks in order to improve predictions of dynamical processes.

preprint2020arXiv

Statistical physics approaches to the complex Earth system

Global climate change, extreme climate events, earthquakes and their accompanying natural disasters pose significant risks to humanity. Yet due to the nonlinear feedbacks, strategic interactions and complex structure of the Earth system, the understanding and in particular the predicting of such disruptive events represent formidable challenges for both scientific and policy communities. During the past years, the emergence and evolution of Earth system science has attracted much attention and produced new concepts and frameworks. Especially, novel statistical physics and complex networks-based techniques have been developed and implemented to substantially advance our knowledge for a better understanding of the Earth system, including climate extreme events, earthquakes and Earth geometric relief features, leading to substantially improved predictive performances. We present here a comprehensive review on the recent scientific progress in the development and application of how combined statistical physics and complex systems science approaches such as, critical phenomena, network theory, percolation, tipping points analysis, as well as entropy can be applied to complex Earth systems (climate, earthquakes, etc.). Notably, these integrating tools and approaches provide new insights and perspectives for understanding the dynamics of the Earth systems. The overall aim of this review is to offer readers the knowledge on how statistical physics approaches can be useful in the field of Earth system science.

preprint2020arXiv

The critical role of fresh teams in creating original and multi-disciplinary research

Teamwork is one of the most prominent features in modern science. It is now well-understood that the team size is an important factor that affects team creativity. However, the crucial question of how the character of research studies is influenced by the freshness of the team remains unclear. In this paper, we quantify the team freshness according to the absent of prior collaboration among team members. Our results suggest that fresher teams tend to produce works of higher originality and more multi-disciplinary impact. These effects are even magnified in larger teams. Furthermore, we find that freshness defined by new team members in a paper is a more effective indicator of research originality and multi-disciplinarity compared to freshness defined by new collaboration relations among team members. Finally, we show that career freshness of members also plays an important role in increasing the originality and multi-disciplinarity of produced papers.

preprint2020arXiv

Two transitions in spatial modular networks

Understanding the resilience of infrastructures such as transportation network has significant importance for our daily life. Recently, a homogeneous spatial network model was developed for studying spatial embedded networks with characteristic link length such as power-grids and the brain. However, although many real-world networks are spatially embedded and their links have characteristics length such as pipelines, power lines or ground transportation lines they are not homogeneous but rather heterogeneous. For example, density of links within cities are significantly higher than between cities. Here we present and study numerically and analytically a similar realistic heterogeneous spatial modular model using percolation process to better understand the effect of heterogeneity on such networks. The model assumes that inside a city there are many lines connecting different locations, while long lines between the cities are sparse and usually directly connecting only a few nearest neighbours cities in a two dimensional plane. We find that this model experiences two distinct continues transitions, one when the cities disconnect from each other and the second when each city breaks apart. Although the critical threshold for site percolation in 2D grid remains an open question we analytically find the critical threshold for site percolation in this model. In addition, while the homogeneous model experience a single transition having a unique phenomenon called \textit{critical stretching} where a geometric crossover from random to spatial structure in different scales found to stretch non-linearly with the characteristic length at criticality. Here we show that the heterogeneous model does not experience such a phenomenon indicating that critical stretching strongly depends on the network structure.

preprint2020arXiv

Very early warning signal for El Niño in 2020 with a 4 in 5 likelihood

The El Niño Southern Oscillation (ENSO) is the most important driver of climate variability and can trigger extreme weather events and disasters in various parts of the globe. Recently we have developed a network approach, which allows forecasting an El Niño event about 1 year ahead. Here we communicate that since 2012 this network approach, which does not involve any fit parameter, correctly predicted the absence of El Niño events in 2012, 2013 and 2017 as well as the onset of the large El Niño event that started in 2014 and ended in 2016. Our model also correctly forecasted the onset of the last El Niño event in 2018. In September 2019, the model indicated the return of El Niño in 2020 with an 80% probability.

preprint2019arXiv

Interconnections between networks act like an external field in first-order percolation transitions

Many interdependent, real-world infrastructures involve interconnections between different communities or cities. Here we study if and how the effects of such interconnections can be described as an external field for interdependent networks experiencing first-order percolation transitions. We find that the critical exponents $γ$ and $δ$, related to the external field can also be defined for first-order transitions but that they have different values than those found for second-order transitions. Surprisingly, we find that both sets of different exponents can be found even within a single model of interdependent networks, depending on the dependency coupling strength. Specifically, the exponent $γ$ in the first-order regime (high coupling) does not obey the fluctuation dissipation theorem, whereas in the continuous regime (for low coupling) it does. Nevertheless, in both cases they satisfy Widom's identity, $δ- 1 = γ/ β$ which further supports the validity of their definitions. Our results provide physical intuition into the nature of the phase transition in interdependent networks and explain the underlying reasons for two distinct sets of exponents.

preprint2019arXiv

Scaling Laws in Earthquake Memory for Interevent Times and Distances

Over the past decades much effort has been devoted towards understanding and forecasting natural hazards. However, earthquake forecasting skill is still very limited and remains a great scientific challenge. The limited earthquake predictability is partly due to the erratic nature of earthquakes and partly to the lack of understanding the underlying mechanisms of earthquakes. To improve our understanding and potential forecasting, here we study the spatial and temporal long-term memory of interevent earthquakes above a certain magnitude using lagged conditional probabilities. We find, in real data, that the lagged conditional probabilities show long-term memory for both the interevent times and interevent distances and that the memory functions obey scaling and decay slowly with time, while, at a characteristic time, the decay crossesover to a fast decay. We also show that the ETAS model, which is often used to forecast earthquake events, yields scaling functions of the temporal and spatial interevent intervals which are not consistent with those of real data.

preprint2019arXiv

Significant Impact of Rossby Waves on Air Pollution Detected by Network Analysis

Air pollution is associated with human diseases and has been found to be related to premature mortality. In response, environmental policies have been adopted in many countries, to decrease anthropogenic air pollution for the improvement of long-term air quality, since most air pollutant sources are anthropogenic. However, air pollution fluctuations have been found to strongly depend on the weather dynamics. This raises a fundamental question: What are the significant atmospheric processes that affect the local daily variability of air pollution? For this purpose, we develop here a multi-layered network analysis to detect the interlinks between the geopotential height of the upper air (~5 km) and surface air pollution in both China and the USA. We find that Rossby waves significantly affect air pollution fluctuations through the development of cyclone and anticyclone systems, and further affect the local stability of the air and the winds. The significant impacts of Rossby waves on air pollution are found to underlie most of the daily fluctuations in air pollution. Thus, the impact of Rossby waves on human life is greater than previously assumed. The rapid warming of the Arctic could slow down Rossby waves, thus increasing human health risks. Our method can help to determine the risk assessment of such extreme events and can improve potential predictability.

preprint2016arXiv

$k$-core percolation on complex networks: Comparing random, localized and targeted attacks

The type of malicious attack inflicting on networks greatly influences their stability under ordinary percolation in which a node fails when it becomes disconnected from the giant component. Here we study its generalization, $k$-core percolation, in which a node fails when it loses connection to a threshold $k$ number of neighbors. We study and compare analytically and by numerical simulations of $k$-core percolation the stability of networks under random attacks (RA), localized attacks (LA) and targeted attacks (TA), respectively. By mapping a network under LA or TA into an equivalent network under RA, we find that in both single and interdependent networks, TA exerts the greatest damage to the core structure of a network. We also find that for Erdős-Rényi (ER) networks, LA and RA exert equal damage to the core structure whereas for scale-free (SF) networks, LA exerts much more damage than RA does to the core structure.

preprint2016arXiv

Climate network suggests enhanced El Niño global impacts in localized areas

We construct directed and weighted climate networks based on near surface air temperature to investigate the global impacts of El Nino and La Nina. We find that regions which are characterized by higher positive or negative network in weighted links, are exhibiting stronger correlations with the El Nino basin and are warmer or cooler during El Nino or La Nina periods. These stronger in-weighted activities are found to be concentrated in localized areas, as compared to non-El Nino periods, whereas a large fraction of the globe is not influenced by the events. The regions of localized activity vary from one El Nino (La Nina) event to another; still some El Nino (La Nina) events are more similar to each other. We quantify this similarity using network community structure. The results and methodology reported here may be used to improve the understanding and prediction of El Nino or La Nina events and also may be applied in the investigation of other climate variables.

preprint2016arXiv

Collective behavior in the spatial spreading of obesity

Non-communicable diseases like diabetes, obesity and certain forms of cancer have been increasing in many countries at alarming levels. A difficulty in the conception of policies to reverse these trends is the identification of the drivers behind the global epidemics. Here, we implement a spatial spreading analysis to investigate whether diabetes, obesity and cancer show spatial correlations revealing the effect of collective and global factors acting above individual choices. We adapt a theoretical framework for critical physical systems displaying collective behavior to decipher the laws of spatial spreading of diseases. We find a regularity in the spatial fluctuations of their prevalence revealed by a pattern of scale-free long-range correlations. The fluctuations are anomalous, deviating in a fundamental way from the weaker correlations found in the underlying population distribution. This collective behavior indicates that the spreading dynamics of obesity, diabetes and some forms of cancer like lung cancer are analogous to a critical point of fluctuations, just as a physical system in a second-order phase transition. According to this notion, individual interactions and habits may have negligible influence in shaping the global patterns of spreading. Thus, obesity turns out to be a global problem where local details are of little importance. Interestingly, we find the same critical fluctuations in obesity and diabetes, and in the activities of economic sectors associated with food production such as supermarkets, food and beverage stores--- which cluster in a different universality class than other generic sectors of the economy. These results motivate future interventions to investigate the causality of this relation providing guidance for the implementation of preventive health policies.

preprint2016arXiv

Gravitational scaling in Beijing Subway Network

Recently, with the availability of various traffic datasets, human mobility has been studied in different contexts. Researchers attempt to understand the collective behaviors of human movement with respect to the spatio-temporal distribution in traffic dynamics, from which a gravitational scaling law characterizing the relation between the traffic flow, population and distance has been found. However, most studies focus on the integrated properties of gravitational scaling, neglecting its dynamical evolution during different hours of a day. Investigating the hourly traffic flow data of Beijing subway network, based on the hop-count distance of passengers, we find that the scaling exponent of the gravitational law is smaller in Beijing subway system compared to that reported in Seoul subway system. This means that traffic demand in Beijing is much stronger and less sensitive to the travel distance. Furthermore, we analyzed the temporal evolution of the scaling exponents in weekdays and weekends. Our findings may help to understand and improve the traffic congestion control in different subway systems.

preprint2016arXiv

Stochastic model of financial markets reproducing scaling and memory in volatility return intervals

We investigate the volatility return intervals in the NYSE and FOREX markets. We explain previous empirical findings using a model based on the interacting agent hypothesis instead of the widely-used efficient market hypothesis. We derive macroscopic equations based on the microscopic herding interactions of agents and find that they are able to reproduce various stylized facts of different markets and different assets with the same set of model parameters. We show that the power-law properties and the scaling of return intervals and other financial variables have a similar origin and could be a result of a general class of non-linear stochastic differential equations derived from a master equation of an agent system that is coupled by herding interactions. Specifically, we find that this approach enables us to recover the volatility return interval statistics as well as volatility probability and spectral densities for the NYSE and FOREX markets, for different assets, and for different time-scales. We find also that the historical S\&P500 monthly series exhibits the same volatility return interval properties recovered by our proposed model. Our statistical results suggest that human herding is so strong that it persists even when other evolving fluctuations perturbate the financial system.

preprint2015arXiv

Debunking in a World of Tribes

Recently a simple military exercise on the Internet was perceived as the beginning of a new civil war in the US. Social media aggregate people around common interests eliciting a collective framing of narratives and worldviews. However, the wide availability of user-provided content and the direct path between producers and consumers of information often foster confusion about causations, encouraging mistrust, rumors, and even conspiracy thinking. In order to contrast such a trend attempts to \textit{debunk} are often undertaken. Here, we examine the effectiveness of debunking through a quantitative analysis of 54 million users over a time span of five years (Jan 2010, Dec 2014). In particular, we compare how users interact with proven (scientific) and unsubstantiated (conspiracy-like) information on Facebook in the US. Our findings confirm the existence of echo chambers where users interact primarily with either conspiracy-like or scientific pages. Both groups interact similarly with the information within their echo chamber. We examine 47,780 debunking posts and find that attempts at debunking are largely ineffective. For one, only a small fraction of usual consumers of unsubstantiated information interact with the posts. Furthermore, we show that those few are often the most committed conspiracy users and rather than internalizing debunking information, they often react to it negatively. Indeed, after interacting with debunking posts, users retain, or even increase, their engagement within the conspiracy echo chamber.

preprint2015arXiv

Long-Range Correlations and Memory in the Dynamics of Internet Interdomain Routing

Data transfer is one of the main functions of the Internet. The Internet consists of a large number of interconnected subnetworks or domains, known as Autonomous Systems. Due to privacy and other reasons the information about what route to use to reach devices within other Autonomous Systems is not readily available to any given Autonomous System. The Border Gateway Protocol is responsible for discovering and distributing this reachability information to all Autonomous Systems. Since the topology of the Internet is highly dynamic, all Autonomous Systems constantly exchange and update this reachability information in small chunks, known as routing control packets or Border Gateway Protocol updates. Motivated by scalability and predictability issues with the dynamics of these updates in the quickly growing Internet, we conduct a systematic time series analysis of Border Gateway Protocol update rates. We find that Border Gateway Protocol update time series are extremely volatile, exhibit long-term correlations and memory effects, similar to seismic time series, or temperature and stock market price fluctuations. The presented statistical characterization of Border Gateway Protocol update dynamics could serve as a ground truth for validation of existing and developing better models of Internet interdomain routing.

preprint2015arXiv

Multiple Tipping Points and Optimal Repairing in Interacting Networks

Systems that comprise many interacting dynamical networks, such as the human body with its biological networks or the global economic network consisting of regional clusters, often exhibit complicated collective dynamics. To understand the collective behavior of such systems, we investigate a model of interacting networks exhibiting the fundamental processes of failure, damage spread, and recovery. We find a very rich phase diagram that becomes exponentially more complex as the number of networks is increased. In the simplest example of $n=2$ interacting networks we find two critical points, 4 triple points, 10 allowed transitions, and two "forbidden" transitions, as well as complex hysteresis loops. Remarkably, we find that triple points play the dominant role in constructing the optimal repairing strategy in damaged interacting systems. To support our model, we analyze an example of real interacting financial networks and find evidence of rapid dynamical transitions between well-defined states, in agreement with the predictions of our model.

preprint2015arXiv

Oceanic El-Niño wave dynamics and climate networks

Oceanic Kelvin and Rossby waves play an important role in tropical climate and \en dynamics. Here we develop and apply a climate network approach to quantify the characteristics of \en related oceanic waves, based on sea surface height satellite data. We associate the majority of dominant long distance ($\geq 500$ km) links of the network with (i) equatorial Kelvin waves, (ii) off-equatorial Rossby waves, and (iii) tropical instability waves. Notably, we find that the location of the hubs of out-going ($\sim 180^{\circ}\rm{E}$) and in-coming ($\sim 140^{\circ}\rm{W}$) links of the climate network coincide with the locations of the Kelvin wave initiation and dissipation, respectively. We also find that this dissipation at $\sim 140^{\circ}\rm{W}$ is much weaker during \en times. Moreover, the hubs of the off-equatorial network coincide with the locations of westerly wind burst activity and high wind vorticity, two mechanisms that were associated with Rossby waves activity. The quantitative methodology and measures developed here can improve the understanding of \en dynamics and possibly its prediction.

preprint2015arXiv

Percolation properties in a traffic model

As a dynamical complex system, traffic is characterized by a transition from free flow to congestions, which is mostly studied in highways. However, despite its importance in developing congestion mitigation strategies, the understanding of this common traffic phenomenon in a city-scale is still missing. An open question is how the traffic in the network collapses from a global efficient traffic to isolated local flows in small clusters, i.e. the question of traffic percolation. Here we study the traffic percolation properties on a lattice by simulation of an agent-based model for traffic. A critical traffic volume in this model distinguishes the free-state from congested state of traffic. Our results show that the threshold of traffic percolation decreases with increasing traffic volume and reaches a minimum value at the critical traffic volume. We show that this minimal threshold is the result of longest spatial correlation between traffic flows at the critical traffic volume. These findings may help to develop congestion mitigation strategies in a network view.

preprint2015arXiv

Resilience of Networks Formed of Interdependent Modular Networks

Many infrastructure networks have a modular structure and are also interdependent. While significant research has explored the resilience of interdependent networks, there has been no analysis of the effects of modularity. Here we develop a theoretical framework for attacks on interdependent modular networks and support our results by simulations. We focus on the case where each network has the same number of communities and the dependency links are restricted to be between pairs of communities of different networks. This is very realistic for infrastructure across cities. Each city has its own infrastructures and different infrastructures are dependent within the city. However, each infrastructure is connected within and between cities. For example, a power grid will connect many cities as will a communication network, yet a power station and communication tower that are interdependent will likely be in the same city. It has been shown that single networks are very susceptible to the failure of the interconnected nodes (between communities) Shai et al. and that attacks on these nodes are more crippling than attacks based on betweenness da Cunha et al. In our example of cities these nodes have long range links which are more likely to fail. For both treelike and looplike interdependent modular networks we find distinct regimes depending on the number of modules, $m$. (i) In the case where there are fewer modules with strong intraconnections, the system first separates into modules in an abrupt first-order transition and then each module undergoes a second percolation transition. (ii) When there are more modules with many interconnections between them, the system undergoes a single transition. Overall, we find that modular structure can influence the type of transitions observed in interdependent networks and should be considered in attempts to make interdependent networks more resilient.

preprint2015arXiv

Spatio-temporal propagation of cascading overload failures

Different from the direct contact in epidemics spread, overload failures propagate through hidden functional dependencies. Many studies focused on the critical conditions and catastrophic consequences of cascading failures. However, to understand the network vulnerability and mitigate the cascading overload failures, the knowledge of how the failures propagate in time and space is essential but still missing. Here we study the spatio-temporal propagation behavior of cascading overload failures analytically and numerically. The cascading overload failures are found to spread radially from the center of the initial failure with an approximately constant velocity. The propagation velocity decreases with increasing tolerance, and can be well predicted by our theoretical framework with one single correction for all the tolerance values. This propagation velocity is found similar in various model networks and real network structures. Our findings may help to predict and mitigate the dynamics of cascading overload failures in realistic systems.

preprint2015arXiv

The influence of the broadness of the degree distribution on network's robustness: comparing localized attack and random attack

The stability of networks is greatly influenced by their degree distributions and in particular by their broadness. Networks with broader degree distributions are usually more robust to random failures but less robust to localized attacks. To better understand the effect of the broadness of the degree distribution we study here two models where the broadness is controlled and compare their robustness against localized attacks (LA) and random attacks (RA). We study analytically and by numerical simulations the cases where the degrees in the networks follow a Bi-Poisson distribution $P(k)=αe^{-λ_1}\frac{λ_1^k}{k!}+(1-α) e^{-λ_2}\frac{λ_2^k}{k!},α\in[0,1]$, and a Gaussian distribution $P(k)=A \cdot exp{(-\frac{(k-μ)^2}{2σ^2})}$ with a normalization constant $A$ where $k\geq 0$. In the Bi-Poisson distribution the broadness is controlled by the values of $α$, $λ_1$ and $λ_2$, while in the Gaussian distribution it is controlled by the standard deviation, $σ$. We find that only for $α=0$ or $α=1$, namely degrees obeying a pure Poisson distribution, LA and RA are the same but for all other cases networks are more vulnerable under LA compared to RA. For Gaussian distribution, with an average degree $μ$ fixed, we find that when $σ^2$ is smaller than $μ$ the network is more vulnerable against random attack. However, when $σ^2$ is larger than $μ$ the network becomes more vulnerable against localized attack. Similar qualitative results are also shown for interdependent networks.

preprint2015arXiv

Two distinct transitions in spatially embedded multiplex networks

Multilayer infrastructure is often interdependent, with nodes in one layer depending on nearby nodes in another layer to function. The links in each layer are often of limited length, due to the construction cost of longer links. Here, we model such systems as a multiplex network composed of two or more layers, each with links of characteristic geographic length, embedded in 2-dimensional space. This is equivalent to a system of interdependent spatially embedded networks in two dimensions in which the connectivity links are constrained in length but varied while the length of the dependency links is always zero. We find two distinct percolation transition behaviors depending on the characteristic length, $ζ$, of the links. When $ζ$ is longer than a certain critical value, $ζ_c$, abrupt, first-order transitions take place, while for $ζ<ζ_c$ the transition is continuous. We show that, though in single-layer networks increasing $ζ$ decreases the percolation threshold $p_c$, in multiplex networks it has the opposite effect: increasing $p_c$ to a maximum at $ζ=ζ_c$. By providing a more realistic topological model for spatially embedded interdependent and multiplex networks and highlighting its similarities to lattice-based models, we provide a new direction for more detailed future studies.

preprint2014arXiv

Classical mechanics of economic networks

Financial networks are dynamic. To assess their systemic importance to the world-wide economic network and avert losses we need models that take the time variations of the links and nodes into account. Using the methodology of classical mechanics and Laplacian determinism we develop a model that can predict the response of the financial network to a shock. We also propose a way of measuring the systemic importance of the banks, which we call BankRank. Using European Bank Authority 2011 stress test exposure data, we apply our model to the bipartite network connecting the largest institutional debt holders of the troubled European countries (Greece, Italy, Portugal, Spain, and Ireland). From simulating our model we can determine whether a network is in a "stable" state in which shocks do not cause major losses, or a "unstable" state in which devastating damages occur. Fitting the parameters of the model, which play the role of physical coupling constants, to Eurozone crisis data shows that before the Eurozone crisis the system was mostly in a "stable" regime, and that during the crisis it transitioned into an "unstable" regime. The numerical solutions produced by our model match closely the actual time-line of events of the crisis. We also find that, while the largest holders are usually more important, in the unstable regime smaller holders also exhibit systemic importance. Our model also proves useful for determining the vulnerability of banks and assets to shocks. This suggests that our model may be a useful tool for simulating the response dynamics of shared portfolio networks.

preprint2014arXiv

Competing for Attention in Social Media under Information Overload Conditions

Although the many forms of modern social media have become major channels for the dissemination of information, they are becoming overloaded because of the rapidly-expanding number of information feeds. We analyze the expanding user-generated content in Sina Weibo, the largest micro-blog site in China, and find evidence that popular messages often follow a mechanism that differs from that found in the spread of disease, in contrast to common believe. In this mechanism, an individual with more friends needs more repeated exposures to spread further the information. Moreover, our data suggest that in contrast to epidemics, for certain messages the chance of an individual to share the message is proportional to the fraction of its neighbours who shared it with him/her. Thus the greater the number of friends an individual has the greater the number of repeated contacts needed to spread the message, which is a result of competition for attention. We model this process using a fractional susceptible infected recovered (FSIR) model, where the infection probability of a node is proportional to its fraction of infected neighbors. Our findings have dramatic implications for information contagion. For example, using the FSIR model we find that real-world social networks have a finite epidemic threshold. This is in contrast to the zero threshold that conventional wisdom derives from disease epidemic models. This means that when individuals are overloaded with excess information feeds, the information either reaches out the population if it is above the critical epidemic threshold, or it would never be well received, leading to only a handful of information contents that can be widely spread throughout the population.

preprint2014arXiv

Conditions for viral influence spreading through multiplex correlated social networks

A fundamental problem in network science is to predict how certain individuals are able to initiate new networks to spring up "new ideas". Frequently, these changes in trends are triggered by a few innovators who rapidly impose their ideas through "viral" influence spreading producing cascades of followers fragmenting an old network to create a new one. Typical examples include the raise of scientific ideas or abrupt changes in social media, like the raise of Facebook.com to the detriment of Myspace.com. How this process arises in practice has not been conclusively demonstrated. Here, we show that a condition for sustaining a viral spreading process is the existence of a multiplex correlated graph with hidden "influence links". Analytical solutions predict percolation phase transitions, either abrupt or continuous, where networks are disintegrated through viral cascades of followers as in empirical data. Our modeling predicts the strict conditions to sustain a large viral spreading via a scaling form of the local correlation function between multilayers, which we also confirm empirically. Ultimately, the theory predicts the conditions for viral cascading in a large class of multiplex networks ranging from social to financial systems and markets.

preprint2014arXiv

Efficiency of message transmission using biased random walks in complex networks in the presence of traps

We study the problem of a particle/message that travels as a biased random walk towards a target node in a network in the presence of traps. The bias is represented as the probability $p$ of the particle to travel along the shortest path to the target node. The efficiency of the transmission process is expressed through the fraction $f_g$ of particles that succeed to reach the target without being trapped. By relating $f_g$ with the number $S$ of nodes visited before reaching the target, we firstly show that, for the unbiased random walk, $f_g$ is inversely proportional to both the concentration $c$ of traps and the size $N$ of the network. For the case of biased walks, a simple approximation of $S$ provides an analytical solution that describes well the behavior of $f_g$, especially for $p>0.5$. Also, it is shown that for a given value of the bias $p$, when the concentration of traps is less than a threshold value equal to the inverse of the Mean First Passage Time (MFPT) between two randomly chosen nodes of the network, the efficiency of transmission is unaffected by the presence of traps and almost all the particles arrive at the target. As a consequence, for a given concentration of traps, we can estimate the minimum bias that is needed to have unaffected transmission, especially in the case of Random Regular (RR), Erdős-Rényi (ER) and Scale-Free (SF) networks, where an exact expression (RR and ER) or an upper bound (SF) of the MFPT is known analytically. We also study analytically and numerically, the fraction $f_g$ of particles that reach the target on SF networks, where a single trap is placed on the highest degree node. For the unbiased random walk, we find that $f_g \sim N^{-1/(γ-1)}$, where $γ$ is the power law exponent of the SF network.

preprint2014arXiv

How does public opinion become extreme?

We investigate the emergence of extreme opinion trends in society by employing statistical physics modeling and analysis on polls that inquire about a wide range of issues such as religion, economics, politics, abortion, extramarital sex, books, movies, and electoral vote. The surveys lay out a clear indicator of the rise of extreme views. The precursor is a nonlinear relation between the fraction of individuals holding a certain extreme view and the fraction of individuals that includes also moderates, e.g., in politics, those who are "very conservative" versus "moderate to very conservative" ones. We propose an activation model of opinion dynamics with interaction rules based on the existence of individual "stubbornness" that mimics empirical observations. According to our modeling, the onset of nonlinearity can be associated to an abrupt bootstrap-percolation transition with cascades of extreme views through society. Therefore, it represents an early-warning signal to forecast the transition from moderate to extreme views. Moreover, by means of a phase diagram we can classify societies according to the percolative regime they belong to, in terms of critical fractions of extremists and people's ties.

preprint2014arXiv

Influence of autocorrelation on the topology of the climate network

Different definitions of links in climate networks may lead to considerably different network topologies. We construct a network from climate records of surface level atmospheric temperature in different geographical sites around the globe using two commonly used definitions of links. Utilizing detrended fluctuation analysis, shuffled surrogates and separation analysis of maritime and continental records, we find that one of the major influences on the structure of climate networks is due to the auto-correlation in the records, that may introduce spurious links. This may explain why different methods could lead to different climate network topologies.

preprint2014arXiv

Non-consensus opinion model on directed networks

Dynamic social opinion models have been widely studied on undirected networks, and most of them are based on spin interaction models that produce a consensus. In reality, however, many networks such as Twitter and the World Wide Web are directed and are composed of both unidirectional and bidirectional links. Moreover, from choosing a coffee brand to deciding who to vote for in an election, two or more competing opinions often coexist. In response to this ubiquity of directed networks and the coexistence of two or more opinions in decision-making situations, we study a non-consensus opinion model introduced by Shao et al. \cite{shao2009dynamic} on directed networks. We define directionality $ξ$ as the percentage of unidirectional links in a network, and we use the linear correlation coefficient $ρ$ between the indegree and outdegree of a node to quantify the relation between the indegree and outdegree. We introduce two degree-preserving rewiring approaches which allow us to construct directed networks that can have a broad range of possible combinations of directionality $ξ$ and linear correlation coefficient $ρ$ and to study how $ξ$ and $ρ$ impact opinion competitions. We find that, as the directionality $ξ$ or the indegree and outdegree correlation $ρ$ increases, the majority opinion becomes more dominant and the minority opinion's ability to survive is lowered.

preprint2014arXiv

Partial correlation analysis: Applications for financial markets

The presence of significant cross-correlations between the synchronous time evolution of a pair of equity returns is a well-known empirical fact. The Pearson correlation is commonly used to indicate the level of similarity in the price changes for a given pair of stocks, but it does not measure whether other stocks influence the relationship between them. To explore the influence of a third stock on the relationship between two stocks, we use a partial correlation measurement to determine the underlying relationships between financial assets. Building on previous work, we present a statistically robust approach to extract the underlying relationships between stocks from four different financial markets: the United States, the United Kingdom, Japan, and India. This methodology provides new insights into financial market dynamics and uncovers implicit influences in play between stocks. To demonstrate the capabilities of this methodology, we (i) quantify the influence of different companies and, by studying market similarity across time, present new insights into market structure and market stability, and (ii) we present a practical application, which provides information on the how a company is influenced by different economic sectors, and how the sectors interact with each other. These examples demonstrate the effectiveness of this methodology in uncovering information valuable for a range of individuals, including not only investors and traders but also regulators and policy makers.

preprint2014arXiv

Percolation of localized attack on complex networks

The robustness of complex networks against node failure and malicious attack has been of interest for decades, while most of the research has focused on random attack or hub-targeted attack. In many real-world scenarios, however, attacks are neither random nor hub-targeted, but localized, where a group of neighboring nodes in a network are attacked and fail. In this paper we develop a percolation framework to analytically and numerically study the robustness of complex networks against such localized attack. In particular, we investigate this robustness in Erdős-Rényi networks, random-regular networks, and scale-free networks. Our results provide insight into how to better protect networks, enhance cybersecurity, and facilitate the design of more robust infrastructures.

preprint2014arXiv

Ranking the Economic Importance of Countries and Industries

In the current era of worldwide stock market interdependencies, the global financial village has become increasingly vulnerable to systemic collapse. The recent global financial crisis has highlighted the necessity of understanding and quantifying interdependencies among the world's economies, developing new effective approaches to risk evaluation, and providing mitigating solutions. We present a methodological framework for quantifying interdependencies in the global market and for evaluating risk levels in the world-wide financial network. The resulting information will enable policy and decision makers to better measure, understand, and maintain financial stability. We use the methodology to rank the economic importance of each industry and country according to the global damage that would result from their failure. Our quantitative results shed new light on China's increasing economic dominance over other economies, including that of the USA, to the global economy.

preprint2014arXiv

Resilience of modular complex networks

Complex networks often have a modular structure, where a number of tightly- connected groups of nodes (modules) have relatively few interconnections. Modularity had been shown to have an important effect on the evolution and stability of biological networks, on the scalability and efficiency of large-scale infrastructure, and the development of economic and social systems. An analytical framework for understanding modularity and its effects on network vulnerability is still missing. Through recent advances in the understanding of multilayer networks, however, it is now possible to develop a theoretical framework to systematically study this critical issue. Here we study, analytically and numerically, the resilience of modular networks under attacks on interconnected nodes, which exhibit high betweenness values and are often more exposed to failure. Our model provides new understandings into the feedback between structure and function in real world systems, and consequently has important implications as diverse as developing efficient immunization strategies, designing robust large-scale infrastructure, and understanding brain function.

preprint2014arXiv

Robustness of a Network Formed of Spatially Embedded Networks

We present analytic and numeric results for percolation in a network formed of interdependent spatially embedded networks. We show results for a treelike and a random regular network of networks each with $(i)$ unconstrained interdependent links and $(ii)$ interdependent links restricted to a maximum length, $r$. Analytic results are given for each network of networks with unconstrained dependency links and compared with simulations. For the case of two spatially embedded networks it was found that only for $r>r_c\approx8$ does the system undergo a first order phase transition. We find that for treelike networks of networks $r_c$ significantly decreases as $n$ increases and rapidly reaches its limiting value, $r=1$. For cases where the dependencies form loops, such as in random regular networks, we show analytically and confirm through simulations, that there is a certain fraction of dependent nodes, $q_{max}$, above which the entire network structure collapses even if a single node is removed. This $q_{max}$ decreases quickly with $m$, the degree of the random regular network of networks. Our results show the extreme sensitivity of coupled spatial networks and emphasize the susceptibility of these networks to sudden collapse. The theory derived here can be used to find the robustness of any network of networks where the profile of percolation of a single network is known.

preprint2014arXiv

Robustness of networks with topologies of dependency links

The robustness of complex networks with dependencies has been studied in recent years. However, previous studies focused on the robustness of networks composed of dependency links without network topology. In this study, we will analyze the percolation properties of a realistic network model where dependency links follow certain network topology. We perform the theoretical analysis and numerical simulations to show the critical effects of topology of dependency links on robustness of complex networks. For Erdös-Rényi (ER) connectivity network, we find that the system with dependency of RR topology is more vulnerable than system with dependency of ER topology. And RR-RR (i.e. random-regular (RR) network with dependency of RR topology) disintegrates in an abrupt transition. In particular, we find that the system of RR-ER shows different types of phase transitions. For system of different combinations, the type of percolation depends on the interaction between connectivity topology and dependency topology.

preprint2014arXiv

Simultaneous first and second order percolation transitions in interdependent networks

In a system of interdependent networks, an initial failure of nodes invokes a cascade of iterative failures that may lead to a total collapse of the whole system in a form of an abrupt first order transition. When the fraction of initial failed nodes $1-p$ reaches criticality, $p=p_c$, the abrupt collapse occurs by spontaneous cascading failures. At this stage, the giant component decreases slowly in a plateau form and the number of iterations in the cascade, $τ$, diverges. The origin of this plateau and its increasing with the size of the system remained unclear. Here we find that simultaneously with the abrupt first order transition a spontaneous second order percolation occurs during the cascade of iterative failures. This sheds light on the origin of the plateau and on how its length scales with the size of the system. Understanding the critical nature of the dynamical process of cascading failures may be useful for designing strategies for preventing and mitigating catastrophic collapses.

preprint2013arXiv

Anomalous biased diffusion in networks

We study diffusion with a bias towards a target node in networks. This problem is relevant to efficient routing strategies in emerging communication networks like optical networks. Bias is represented by a probability $p$ of the packet/particle to travel at every hop towards a site which is along the shortest path to the target node. We investigate the scaling of the mean first passage time (MFPT) with the size of the network. We find by using theoretical analysis and computer simulations that for Random Regular (RR) and Erdős-Rényi (ER) networks, there exists a threshold probability, $p_{th}$, such that for $p<p_{th}$ the MFPT scales anomalously as $N^α$, where $N$ is the number of nodes, and $α$ depends on $p$. For $p>p_{th}$ the MFPT scales logarithmically with $N$. The threshold value $p_{th}$ of the bias parameter for which the regime transition occurs is found to depend only on the mean degree of the nodes. An exact solution for every value of $p$ is given for the scaling of the MFPT in RR networks. The regime transition is also observed for the second moment of the probability distribution function, the standard deviation.

preprint2013arXiv

Cascading Failures in Bi-partite Graphs: Model for Systemic Risk Propagation

As economic entities become increasingly interconnected, a shock in a financial network can provoke significant cascading failures throughout the system. To study the systemic risk of financial systems, we create a bi-partite banking network model composed of banks and bank assets and propose a cascading failure model to describe the risk propagation process during crises. We empirically test the model with 2007 US commercial banks balance sheet data and compare the model prediction of the failed banks with the real failed banks after 2007. We find that our model efficiently identifies a significant portion of the actual failed banks reported by Federal Deposit Insurance Corporation. The results suggest that this model could be useful for systemic risk stress testing for financial systems. The model also identifies that commercial rather than residential real estate assets are major culprits for the failure of over 350 US commercial banks during 2008-2011.

preprint2013arXiv

Effect of the Interconnected Network Structure on the Epidemic Threshold

Most real-world networks are not isolated. In order to function fully, they are interconnected with other networks, and this interconnection influences their dynamic processes. For example, when the spread of a disease involves two species, the dynamics of the spread within each species (the contact network) differs from that of the spread between the two species (the interconnected network). We model two generic interconnected networks using two adjacency matrices, A and B, in which A is a 2N*2N matrix that depicts the connectivity within each of two networks of size N, and B a 2N*2N matrix that depicts the interconnections between the two. Using an N-intertwined mean-field approximation, we determine that a critical susceptable-infected-susceptable (SIS) epidemic threshold in two interconnected networks is 1/λ1(A+αB), where the infection rate is βwithin each of the two individual networks and αβin the interconnected links between the two networks and λ1(A+αB) is the largest eigenvalue of the matrix A+αB. In order to determine how the epidemic threshold is dependent upon the structure of interconnected networks, we analytically derive λ1(A+αB) using perturbation approximation for small and large α, the lower and upper bound for any αas a function of the adjacency matrix of the two individual networks, and the interconnections between the two and their largest eigenvalues/eigenvectors. We verify these approximation and boundary values for λ1(A+αB) using numerical simulations, and determine how component network features affect λ1(A+αB).

preprint2013arXiv

Efficient learning strategy of Chinese characters based on network approach

Based on network analysis of hierarchical structural relations among Chinese characters, we develop an efficient learning strategy of Chinese characters. We regard a more efficient learning method if one learns the same number of useful Chinese characters in less effort or time. We construct a node-weighted network of Chinese characters, where character usage frequencies are used as node weights. Using this hierarchical node-weighted network, we propose a new learning method, the distributed node weight (DNW) strategy, which is based on a new measure of nodes' importance that takes into account both the weight of the nodes and the hierarchical structure of the network. Chinese character learning strategies, particularly their learning order, are analyzed as dynamical processes over the network. We compare the efficiency of three theoretical learning methods and two commonly used methods from mainstream Chinese textbooks, one for Chinese elementary school students and the other for students learning Chinese as a second language. We find that the DNW method significantly outperforms the others, implying that the efficiency of current learning methods of major textbooks can be greatly improved.

preprint2013arXiv

Global Climate network evolves with North Atlantic Oscillation phases: Coupling to Southern Pacific Ocean

We construct a network from climate records of atmospheric temperature at surface level, at different geographical sites in the globe, using reanalysis data from years 1948-2010. We find that the network correlates with the North Atlantic Oscillation (NAO), both locally in the north Atlantic, and through coupling to the southern Pacific Ocean. The existence of tele-connection links between those areas and their stability over time allows us to suggest a possible physical explanation for this phenomenon.

preprint2013arXiv

Improved El Niño-Forecasting by Cooperativity Detection

Although anomalous episodical warming of the eastern equatorial Pacific, dubbed El Niño by Peruvian fishermen, has major (and occasionally devastating) impacts around the globe, robust forecasting is still limited to about six months ahead. A significant extension of the pre-warming time would be instrumental for avoiding some of the worst damages such as harvest failures in developing countries. Here we introduce a novel avenue towards El Niño-prediction based on network methods inspecting emerging teleconnections. Our approach starts from the evidence that a large-scale cooperative mode - linking the El Niño-basin (equatorial Pacific corridor) and the rest of the ocean - builds up in the calendar year before the warming event. On this basis, we can develop an efficient 12 months-forecasting scheme, i.e., achieve some doubling of the early-warning period. Our method is based on high-quality observational data as available since 1950 and yields hit rates above 0.5, while false-alarm rates are below 0.1.

preprint2013arXiv

Origins of power-law degree distribution in the heterogeneity of human activity in social networks

The probability distribution of number of ties of an individual in a social network follows a scale-free power-law. However, how this distribution arises has not been conclusively demonstrated in direct analyses of people's actions in social networks. Here, we perform a causal inference analysis and find an underlying cause for this phenomenon. Our analysis indicates that heavy-tailed degree distribution is causally determined by similarly skewed distribution of human activity. Specifically, the degree of an individual is entirely random - following a "maximum entropy attachment" model - except for its mean value which depends deterministically on the volume of the users' activity. This relation cannot be explained by interactive models, like preferential attachment, since the observed actions are not likely to be caused by interactions with other people.

preprint2013arXiv

Percolation of a general network of networks

Percolation theory is an approach to study vulnerability of a system. We develop analytical framework and analyze percolation properties of a network composed of interdependent networks (NetONet). Typically, percolation of a single network shows that the damage in the network due to a failure is a continuous function of the fraction of failed nodes. In sharp contrast, in NetONet, due to the cascading failures, the percolation transition may be discontinuous and even a single node failure may lead to abrupt collapse of the system. We demonstrate our general framework for a NetONet composed of $n$ classic Erdős-Rényi (ER) networks, where each network depends on the same number $m$ of other networks, i.e., a random regular network of interdependent ER networks. In contrast to a \emph{treelike} NetONet in which the size of the largest connected cluster (mutual component) depends on $n$, the loops in the RR NetONet cause the largest connected cluster to depend only on $m$. We also analyzed the extremely vulnerable feedback condition of coupling. In the case of ER networks, the NetONet only exhibits two phases, a second order phase transition and collapse, and there is no first phase transition regime unlike the no feedback condition. In the case of NetONet composed of RR networks, there exists a first order phase transition when $q$ is large and second order phase transition when $q$ is small. Our results can help in designing robust interdependent systems.

preprint2013arXiv

Percolation of Interdependent Networks with Inter-similarity

Real data show that interdependent networks usually involve inter-similarity. Intersimilarity means that a pair of interdependent nodes have neighbors in both networks that are also interdependent (Parshani et al \cite{PAR10B}). For example, the coupled world wide port network and the global airport network are intersimilar since many pairs of linked nodes (neighboring cities), by direct flights and direct shipping lines exist in both networks. Nodes in both networks in the same city are regarded as interdependent. If two neighboring nodes in one network depend on neighboring nodes in the another we call these links common links. The fraction of common links in the system is a measure of intersimilarity. Previous simulation results suggest that intersimilarity has considerable effect on reducing the cascading failures, however, a theoretical understanding on this effect on the cascading process is currently missing. Here, we map the cascading process with inter-similarity to a percolation of networks composed of components of common links and non common links. This transforms the percolation of inter-similar system to a regular percolation on a series of subnetworks, which can be solved analytically. We apply our analysis to the case where the network of common links is an Erdős-Rényi (ER) network with the average degree $K$, and the two networks of non-common links are also ER networks. We show for a fully coupled pair of ER networks, that for any $K\geq0$, although the cascade is reduced with increasing $K$, the phase transition is still discontinuous. Our analysis can be generalized to any kind of interdependent random networks system.

preprint2013arXiv

Percolation of Partially Interdependent Scale-free Networks

We study the percolation behavior of two interdependent scale-free (SF) networks under random failure of 1-$p$ fraction of nodes. Our results are based on numerical solutions of analytical expressions and simulations. We find that as the coupling strength between the two networks $q$ reduces from 1 (fully coupled) to 0 (no coupling), there exist two critical coupling strengths $q_1$ and $q_2$, which separate three different regions with different behavior of the giant component as a function of $p$. (i) For $q \geq q_1$, an abrupt collapse transition occurs at $p=p_c$. (ii) For $q_2<q<q_1$, the giant component has a hybrid transition combined of both, abrupt decrease at a certain $p=p^{jump}_c$ followed by a smooth decrease to zero for $p < p^{jump}_c$ as $p$ decreases to zero. (iii) For $q \leq q_2$, the giant component has a continuous second-order transition (at $p=p_c$). We find that $(a)$ for $λ\leq 3$, $q_1 \equiv 1$; and for $λ> 3$, $q_1$ decreases with increasing $λ$. $(b)$ In the hybrid transition, at the $q_2 < q < q_1$ region, the mutual giant component $P_{\infty}$ jumps discontinuously at $p=p^{jump}_c$ to a very small but non-zero value, and when reducing $p$, $P_{\infty}$ continuously approaches to 0 at $p_c = 0$ for $λ< 3$ and at $p_c > 0$ for $λ> 3$. Thus, the known theoretical $p_c=0$ for a single network with $λ\leqslant 3$ is expected to be valid also for strictly partial interdependent networks.

preprint2013arXiv

Random walk with priorities in communication-like networks

We study a model for a random walk of two classes of particles (A and B). Where both species are present in the same site, the motion of A's takes precedence over that of B's. The model was originally proposed and analyzed in Maragakis et al., Phys. Rev. E 77, 020103 (2008); here we provide additional results. We solve analytically the diffusion coefficients of the two species in lattices for a number of protocols. In networks, we find that the probability of a B particle to be free decreases exponentially with the node degree. In scale-free networks, this leads to localization of the B's at the hubs and arrest of their motion. To remedy this, we investigate several strategies to avoid trapping of the B's: moving an A instead of the hindered B; allowing a trapped B to hop with a small probability; biased walk towards non-hub nodes; and limiting the capacity of nodes. We obtain analytic results for lattices and networks, and discuss the advantages and shortcomings of the possible strategies.

preprint2013arXiv

Robustness of partially interdependent network formed of clustered networks

Clustering, or transitivity has been observed in real networks and its effects on their structure and function has been discussed extensively. The focus of these studies has been on clustering of single networks while the effect of clustering on the robustness of coupled networks received very little attention. Only the case of a pair of fully coupled networks with clustering has been studied recently. Here we generalize the study of clustering of a fully coupled pair of networks to the study of partially interdependent network of networks with clustering within the network components. We show both analytically and numerically, how clustering within the networks, affects the percolation properties of interdependent networks, including percolation threshold, size of giant component and critical coupling point where first order phase transition changes to second order phase transition as the coupling between the networks reduces. We study two types of clustering: one type proposed by Newman where the average degree is kept constant while changing the clustering and the other proposed by Hackett $et$ $al.$ where the degree distribution is kept constant. The first type of clustering is treated both analytically and numerically while the second one is treated only numerically.

preprint2013arXiv

Spatially localized attacks on interdependent networks: the existence of a finite critical attack size

Many real world complex systems such as infrastructure, communication and transportation networks are embedded in space, where entities of one system may depend on entities of other systems. These systems are subject to geographically localized failures due to malicious attacks or natural disasters. Here we study the resilience of a system composed of two interdependent spatially embedded networks to localized geographical attacks. We find that if an attack is larger than a finite (zero fraction of the system) critical size, it will spread through the entire system and lead to its complete collapse. If the attack is below the critical size, it will remain localized. In contrast, under random attack a finite fraction of the system needs to be removed to initiate system collapse. We present both numerical simulations and a theoretical approach to analyze and predict the effect of local attacks and the critical attack size. Our results demonstrate the high risk of local attacks on interdependent spatially embedded infrastructures and can be useful for designing more resilient systems.

preprint2013arXiv

The dominant imprint of Rossby waves in the climate network

The connectivity pattern of networks, which are based on a correlation between ground level temperature time series, shows a dominant dense stripe of links in the southern ocean. We show that statistical categorization of these links yields a clear association with the pattern of an atmospheric Rossby wave, one of the major mechanisms associated with the weather system and with planetary scale energy transport. It is shown that alternating densities of negative and positive links (correlations) are arranged in half Rossby wave distances around 3,500 km, 7,000 km and 10,000 km and are aligned with the expected direction of energy flow, distribution of time delays and the seasonality of these waves. It is also shown that long distance links (i.e., of distances larger than 2,000 km) that are associated with Rossby waves are the most dominant in the climate network. Climate networks may thus be used as an efficient new way to detect and analyze Rossby waves, based on reliable and available ground level measurements, in addition to the frequently used 300 hPa reanalysis meridional wind data.

preprint2013arXiv

Towards designing robust coupled networks

Natural and technological interdependent systems have been shown to be highly vulnerable due to cascading failures and an abrupt collapse of global connectivity under initial failure. Mitigating the risk by partial disconnection endangers their functionality. Here we propose a systematic strategy of selecting a minimum number of autonomous nodes that guarantee a smooth transition in robustness. Our method which is based on betweenness is tested on various examples including the famous 2003 electrical blackout of Italy. We show that, with this strategy, the necessary number of autonomous nodes can be reduced by a factor of five compared to a random choice. We also find that the transition to abrupt collapse follows tricritical scaling characterized by a set of exponents which is independent on the protection strategy.

preprint2012arXiv

A k-shell decomposition method for weighted networks

We present a generalized method for calculating the k-shell structure of weighted networks. The method takes into account both the weight and the degree of a network, in such a way that in the absence of weights we resume the shell structure obtained by the classic k-shell decomposition. In the presence of weights, we show that the method is able to partition the network in a more refined way, without the need of any arbitrary threshold on the weight values. Furthermore, by simulating spreading processes using the susceptible-infectious-recovered model in four different weighted real-world networks, we show that the weighted k-shell decomposition method ranks the nodes more accurately, by placing nodes with higher spreading potential into shells closer to the core. In addition, we demonstrate our new method on a real economic network and show that the core calculated using the weighted k-shell method is more meaningful from an economic perspective when compared with the unweighted one.

preprint2012arXiv

Cascading Failures in Interdependent Lattice Networks: The Critical Role of the Length of Dependency Links

We study the cascading failures in a system composed of two interdependent square lattice networks A and B placed on the same Cartesian plane, where each node in network A depends on a node in network B randomly chosen within a certain distance $r$ from the corresponding node in network A and vice versa. Our results suggest that percolation for small $r$ below $r_{\rm max}\approx 8$ (lattice units) is a second-order transition, and for larger $r$ is a first-order transition. For $r<r_{\rm max}$, the critical threshold increases linearly with $r$ from 0.593 at $r=0$ and reaches a maximum, 0.738 for $r=r_{\rm max}$ and then gradually decreases to 0.683 for $r=\infty$. Our analytical considerations are in good agreement with simulations. Our study suggests that interdependent infrastructures embedded in Euclidean space become most vulnerable when the distance between interdependent nodes is in the intermediate range, which is much smaller than the size of the system.

preprint2012arXiv

Communication activity in a social network: relation between long-term correlations and inter-event clustering

The timing patterns of human communication in social networks is not random. On the contrary, communication is dominated by emergent statistical laws such as non-trivial correlations and clustering. Recently, we found long-term correlations in the user's activity in social communities. Here, we extend this work to study collective behavior of the whole community. The goal is to understand the origin of clustering and long-term persistence. At the individual level, we find that the correlations in activity are a byproduct of the clustering expressed in the power-law distribution of inter-event times of single users. On the contrary, the activity of the whole community presents long-term correlations that are a true emergent property of the system, i.e. they are not related to the distribution of inter-event times. This result suggests the existence of collective behavior, possible arising from nontrivial communication patterns through the embedding social network.

preprint2012arXiv

Complex networks embedded in space: Dimension and scaling relations between mass, topological distance and Euclidean distance

Many real networks are embedded in space, where in some of them the links length decay as a power law distribution with distance. Indications that such systems can be characterized by the concept of dimension were found recently. Here, we present further support for this claim, based on extensive numerical simulations for model networks embedded on lattices of dimensions $d_e=1$ and $d_e=2$. We evaluate the dimension $d$ from the power law scaling of (a) the mass of the network with the Euclidean radius $r$ and (b) the probability of return to the origin with the distance $r$ travelled by the random walker. Both approaches yield the same dimension. For networks with $δ< d_e$, $d$ is infinity, while for $δ> 2d_e$, $d$ obtains the value of the embedding dimension $d_e$. In the intermediate regime of interest $d_e \leq δ< 2 d_e$, our numerical results suggest that $d$ decreases continously from $d = \infty$ to $d_e$, with $d - d_e \sim (δ- d_e)^{-1}$ for $δ$ close to $d_e$. Finally, we discuss the scaling of the mass $M$ and the Euclidean distance $r$ with the topological distance $\ell$. Our results suggest that in the intermediate regime $d_e \leq δ< 2 d_e$, $M(\ell)$ and $r(\ell)$ do not increase with $\ell$ as a power law but with a stretched exponential, $M(\ell) \sim \exp [A \ell^{δ' (2 - δ')}]$ and $r(\ell) \sim \exp [B \ell^{δ' (2 - δ')}]$, where $δ' = δ/d_e$. The parameters $A$ and $B$ are related to $d$ by $d = A/B$, such that $M(\ell) \sim r(\ell)^d$. For $δ< d_e$, $M$ increases exponentially with $\ell$, as known for $δ=0$, while $r$ is constant and independent of $\ell$. For $δ\geq 2d_e$, we find power law scaling, $M(\ell) \sim \ell^{d_\ell}$ and $r(\ell) \sim \ell^{1/d_{min}}$, with $d_\ell \cdot d_{min} = d$.

preprint2012arXiv

Earthquake networks based on similar activity patterns

Earthquakes are a complex spatiotemporal phenomenon, the underlying mechanism for which is still not fully understood despite decades of research and analysis. We propose and develop a network approach to earthquake events. In this network, a node represents a spatial location while a link between two nodes represents similar activity patterns in the two different locations. The strength of a link is proportional to the strength of the cross-correlation in activities of two nodes joined by the link. We apply our network approach to a Japanese earthquake catalog spanning the 14-year period 1985-1998. We find strong links representing large correlations between patterns in locations separated by more than 1000 km, corroborating prior observations that earthquake interactions have no characteristic length scale. We find network characteristics not attributable to chance alone, including a large number of network links, high node assortativity, and strong stability over time.

preprint2012arXiv

Languages cool as they expand: Allometric scaling and the decreasing need for new words

We analyze the occurrence frequencies of over 15 million words recorded in millions of books published during the past two centuries in seven different languages. For all languages and chronological subsets of the data we confirm that two scaling regimes characterize the word frequency distributions, with only the more common words obeying the classic Zipf law. Using corpora of unprecedented size, we test the allometric scaling relation between the corpus size and the vocabulary size of growing languages to demonstrate a decreasing marginal need for new words, a feature that is likely related to the underlying correlations between words. We calculate the annual growth fluctuations of word use which has a decreasing trend as the corpus size increases, indicating a slowdown in linguistic evolution following language expansion. This "cooling pattern" forms the basis of a third statistical regularity, which unlike the Zipf and the Heaps law, is dynamical in nature.

preprint2012arXiv

Network Physiology reveals relations between network topology and physiological function

The human organism is an integrated network where complex physiologic systems, each with its own regulatory mechanisms, continuously interact, and where failure of one system can trigger a breakdown of the entire network. Identifying and quantifying dynamical networks of diverse systems with different types of interactions is a challenge. Here, we develop a framework to probe interactions among diverse systems, and we identify a physiologic network. We find that each physiologic state is characterized by a specific network structure, demonstrating a robust interplay between network topology and function. Across physiologic states the network undergoes topological transitions associated with fast reorganization of physiologic interactions on time scales of a few minutes, indicating high network flexibility in response to perturbations. The proposed system-wide integrative approach may facilitate the development of a new field, Network Physiology.

preprint2012arXiv

Non-consensus opinion models on complex networks

We focus on non-consensus opinion models in which above a certain threshold two opinions coexist in a stable relationship. We revisit and extend the non-consensus opinion (NCO) model introduced by Shao. We generalize the NCO model by adding a weight factor W to individual's own opinion when determining its future opinion (NCOW model). We find that as W increases the minority opinion holders tend to form stable clusters with a smaller initial minority fraction compared to the NCO model. We also revisit another non-consensus opinion, the inflexible contrarian opinion (ICO) model, which introduces inflexible contrarians to model a competition between two opinions in the steady state. In the ICO model, the inflexible contrarians effectively decrease the size of the largest cluster of the rival opinion. All of the above models have previously been explored in terms of a single network. However opinions propagate not only within single networks but also between networks, we study here the opinion dynamics in coupled networks. We apply the NCO rule on each individual network and the global majority rule on interdependent pairs. We find that the interdependent links effectively force the system from a second order phase transition, which is characteristic of the NCO model on a single network, to a hybrid phase transition, i.e., a mix of second-order and abrupt jump-like transitions that ultimately becomes, as we increase the percentage of interdependent agents, a pure abrupt transition. We conclude that for the NCO model on coupled networks, interactions through interdependent links could push the non-consensus opinion type model to a consensus opinion type model, which mimics the reality that increased mass communication causes people to hold opinions that are increasingly similar.

preprint2012arXiv

Robustness of onion-like correlated networks against targeted attacks

Recently, it was found by Schneider et al. [Proc. Natl. Acad. Sci. USA, 108, 3838 (2011)], using simulations, that scale-free networks with "onion structure" are very robust against targeted high degree attacks. The onion structure is a network where nodes with almost the same degree are connected. Motivated by this work, we propose and analyze, based on analytical considerations, an onion-like candidate for a nearly optimal structure against simultaneous random and targeted high degree node attacks. The nearly optimal structure can be viewed as a hierarchically interconnected random regular graphs, the degrees and populations of which are specified by the degree distribution. This network structure exhibits an extremely assortative degree-degree correlation and has a close relationship to the "onion structure." After deriving a set of exact expressions that enable us to calculate the critical percolation threshold and the giant component of a correlated network for an arbitrary type of node removal, we apply the theory to the cases of random scale-free networks that are highly vulnerable against targeted high degree node removal. Our results show that this vulnerability can be significantly reduced by implementing this onion-like type of degree-degree correlation without much undermining the almost complete robustness against random node removal. We also investigate in detail the robustness enhancement due to assortative degree-degree correlation by introducing a joint degree-degree probability matrix that interpolates between an uncorrelated network structure and the onion-like structure proposed here by tuning a single control parameter. The optimal values of the control parameter that maximize the robustness against simultaneous random and targeted attacks are also determined. Our analytical calculations are supported by numerical simulations.

preprint2012arXiv

Statistical Laws Governing Fluctuations in Word Use from Word Birth to Word Death

We analyze the dynamic properties of 10^7 words recorded in English, Spanish and Hebrew over the period 1800--2008 in order to gain insight into the coevolution of language and culture. We report language independent patterns useful as benchmarks for theoretical models of language evolution. A significantly decreasing (increasing) trend in the birth (death) rate of words indicates a recent shift in the selection laws governing word use. For new words, we observe a peak in the growth-rate fluctuations around 40 years after introduction, consistent with the typical entry time into standard dictionaries and the human generational timescale. Pronounced changes in the dynamics of language during periods of war shows that word correlations, occurring across time and between words, are largely influenced by coevolutionary social, technological, and political factors. We quantify cultural memory by analyzing the long-term correlations in the use of individual words using detrended fluctuation analysis.

preprint2012arXiv

The extreme vulnerability of interdependent spatially embedded networks

Recent studies show that in interdependent networks a very small failure in one network may lead to catastrophic consequences. Above a critical fraction of interdependent nodes, even a single node failure can invoke cascading failures that may abruptly fragment the system, while below this "critical dependency" (CD) a failure of few nodes leads only to small damage to the system. So far, the research has been focused on interdependent random networks without space limitations. However, many real systems, such as power grids and the Internet, are not random but are spatially embedded. Here we analytically and numerically analyze the stability of systems consisting of interdependent spatially embedded networks modeled as lattice networks. Surprisingly, we find that in lattice systems, in contrast to non-embedded systems, there is no CD and \textit{any} small fraction of interdependent nodes leads to an abrupt collapse. We show that this extreme vulnerability of very weakly coupled lattices is a consequence of the critical exponent describing the percolation transition of a single lattice. Our results are important for understanding the vulnerabilities and for designing robust interdependent spatial embedded networks.

preprint2012arXiv

The robustness of interdependent clustered networks

It was recently found that cascading failures can cause the abrupt breakdown of a system of interdependent networks. Using the percolation method developed for single clustered networks by Newman [Phys. Rev. Lett. {\bf 103}, 058701 (2009)], we develop an analytical method for studying how clustering within the networks of a system of interdependent networks affects the system's robustness. We find that clustering significantly increases the vulnerability of the system, which is represented by the increased value of the percolation threshold $p_c$ in interdependent networks.

preprint2011arXiv

A dynamical phase transition in ferromagnetic granular materials

We study, using simulations the dynamical properties of complex ferromagnetic granular materials. The system of grains is modeled by a disordered two-dimensional lattice in which the grains are embedded, while the magnitude and direction of the easy axis are random. Using the monte-carlo method we track the dynamics of the magnetic moments of the grains. We observe a transition of the system from a macroscopic blocked (ferromagnetic) phase at low temperature in which the grain's magnetic moment do not flip to the other direction to an unblocked (superparamagnetic) phase at high temperature in which the magnetic moment is free to rotate. Our results suggest that this transition exhibits the characteristics of a second order phase transition such as the appearance of a giant cluster of unblocked grains which is fractal at the critical temperature, a peak in the size of the second largest cluster at the same temperature and a power law distribution of cluster sizes near the criticality.

preprint2011arXiv

Communication activity in social networks: growth and correlations

We investigate the timing of messages sent in two online communities with respect to growth fluctuations and long-term correlations. We find that the timing of sending and receiving messages comprises pronounced long-term persistence. Considering the activity of the community members as growing entities, i.e. the cumulative number of messages sent (or received) by the individuals, we identify non-trivial scaling in the growth fluctuations which we relate to the long-term correlations. We find a connection between the scaling exponents of the growth and the long-term correlations which is supported by numerical simulations based on peaks over threshold. In addition, we find that the activity on directed links between pairs of members exhibits long-term correlations, indicating that communication activity with the most liked partners may be responsible for the long-term persistence in the timing of messages. Finally, we show that the number of messages, $M$, and the number of communication partners, $K$, of the individual members are correlated following a power-law, $K\sim M^λ$, with exponent $λ\approx 3/4$.

preprint2011arXiv

Financial factor influence on scaling and memory of trading volume in stock market

We study the daily trading volume volatility of 17,197 stocks in the U.S. stock markets during the period 1989--2008 and analyze the time return intervals $τ$ between volume volatilities above a given threshold q. For different thresholds q, the probability density function P_q(τ) scales with mean interval <τ> as P_q(τ)=<τ>^{-1}f(τ/<τ>) and the tails of the scaling function can be well approximated by a power-law f(x)~x^{-γ}. We also study the relation between the form of the distribution function P_q(τ) and several financial factors: stock lifetime, market capitalization, volume, and trading value. We find a systematic tendency of P_q(τ) associated with these factors, suggesting a multi-scaling feature in the volume return intervals. We analyze the conditional probability P_q(τ|τ_0) for $τ$ following a certain interval τ_0, and find that P_q(τ|τ_0) depends on τ_0 such that immediately following a short/long return interval a second short/long return interval tends to occur. We also find indications that there is a long-term correlation in the daily volume volatility. We compare our results to those found earlier for price volatility.

preprint2011arXiv

How people interact in evolving online affiliation networks

The study of human interactions is of central importance for understanding the behavior of individuals, groups and societies. Here, we observe the formation and evolution of networks by monitoring the addition of all new links and we analyze quantitatively the tendencies used to create ties in these evolving online affiliation networks. We first show that an accurate estimation of these probabilistic tendencies can only be achieved by following the time evolution of the network. For example, actions that are attributed to the usual friend of a friend mechanism through a static snapshot of the network are overestimated by a factor of two. A detailed analysis of the dynamic network evolution shows that half of those triangles were generated through other mechanisms, in spite of the characteristic static pattern. We start by characterizing every single link when the tie was established in the network. This allows us to describe the probabilistic tendencies of tie formation and extract sociological conclusions as follows. The tendencies to add new links differ significantly from what we would expect if they were not affected by the individuals' structural position in the network, i.e., from random link formation. We also find significant differences in behavioral traits among individuals according to their degree of activity, gender, age, popularity and other attributes. For instance, in the particular datasets analyzed here, we find that women reciprocate connections three times as much as men and this difference increases with age. Men tend to connect with the most popular people more often than women across all ages. On the other hand, triangular ties tendencies are similar and independent of gender. Our findings can be useful to build models of realistic social network structures and discover the underlying laws that govern establishment of ties in evolving social networks.

preprint2011arXiv

Identification of influential spreaders in complex networks

Networks portray a multitude of interactions through which people meet, ideas are spread, and infectious diseases propagate within a society. Identifying the most efficient "spreaders" in a network is an important step to optimize the use of available resources and ensure the more efficient spread of information. Here we show that, in contrast to common belief, the most influential spreaders in a social network do not correspond to the best connected people or to the most central people (high betweenness centrality). Instead, we find: (i) The most efficient spreaders are those located within the core of the network as identified by the k-shell decomposition analysis. (ii) When multiple spreaders are considered simultaneously, the distance between them becomes the crucial parameter that determines the extend of the spreading. Furthermore, we find that-- in the case of infections that do not confer immunity on recovered individuals-- the infection persists in the high k-shell layers of the network under conditions where hubs may not be able to preserve the infection. Our analysis provides a plausible route for an optimal design of efficient dissemination strategies.

preprint2011arXiv

Mitigation of Malicious Attacks on Networks

Terrorist attacks on transportation networks have traumatized modern societies. With a single blast, it has become possible to paralyze airline traffic, electric power supply, ground transportation or Internet communication. How and at which cost can one restructure the network such that it will become more robust against a malicious attack? We introduce a unique measure for robustness and use it to devise a method to mitigate economically and efficiently this risk. We demonstrate its efficiency on the European electricity system and on the Internet as well as on complex networks models. We show that with small changes in the network structure (low cost) the robustness of diverse networks can be improved dramatically while their functionality remains unchanged. Our results are useful not only for improving significantly with low cost the robustness of existing infrastructures but also for designing economically robust network systems.

preprint2011arXiv

Percolation in Interdependent and Interconnected Networks: Abrupt Change from Second to First Order Transition

Robustness of two coupled networks system has been studied only for dependency coupling (S. Buldyrev et. al., Nature, 2010) and only for connectivity coupling (E. A. Leicht and R. M. D'Souza, arxiv:09070894). Here we study, using a percolation approach, a more realistic coupled networks system where both interdependent and interconnected links exist. We find a rich and unusual phase transition phenomena including hybrid transition of mixed first and second order i.e., discontinuities like a first order transition of the giant component followed by a continuous decrease to zero like a second order transition. Moreover, we find unusual discontinuous changes from second order to first order transition as a function of the dependency coupling between the two networks.

preprint2011arXiv

Percolation in networks composed of connectivity and dependency links

Networks composed from both connectivity and dependency links were found to be more vulnerable compared to classical networks with only connectivity links. Their percolation transition is usually of a first order compared to the second order transition found in classical networks. We analytically analyze the effect of different distributions of dependencies links on the robustness of networks. For a random Erd$\ddot{o}$s-R$\acute{e}$nyi (ER) network with average degree $k$ that is divided into dependency clusters of size $s$, the fraction of nodes that belong to the giant component, $P_\infty$, is given by $ P_\infty=p^{s-1} [1-\exp{(-kpP_\infty)}]^s $ where $1-p$ is the initial fraction of removed nodes. Our general result coincides with the known Erd$\ddot{o}$s-R$\acute{e}$nyi equation for random networks for $s=1$ and with the result of Parshani et al (PNAS, in press, 2011) for $s=2$. For networks with Poissonian distribution of dependency links we find that $P_\infty$ is given by $P_\infty = f_{k,p}(P_\infty) e^{(<s>-1)(pf_{k,p}(P_\infty)-1)}$ where $f_{k,p}(P_\infty) \equiv 1-\exp{(-kpP_\infty)}$ and $<s>$ is the mean value of the size of dependency clusters. For networks with Gaussian distribution of dependency links we show how the average and width of the distribution affect the robustness of the networks.

preprint2011arXiv

Stability of Climate Networks with Time

We construct and analyze climate networks based on daily satellite measurements of temperatures and geopotential heights. We show that these networks are stable during time and are similar over different altitudes. Each link in our network is stable with typical 15% variability. The entire hierarchy of links is about 80% consistent during time. We show that about half of this stability is due to the spatial 2D embedding of the network, and half is due to physical coupling mechanisms. The network stability of equatorial regions is found to be lower compared to the stability of a typical network in non-equatorial regions.

preprint2011arXiv

Strategy of Competition between Two Groups based on a Contrarian Opinion Model

We introduce a contrarian opinion (CO) model in which a fraction p of contrarians within a group holds a strong opinion opposite to the opinion held by the rest of the group. At the initial stage, stable clusters of two opinions, A and B exist. Then we introduce contrarians which hold a strong B opinion into the opinion A group. Through their interactions, the contrarians are able to decrease the size of the largest A opinion cluster, and even destroy it. We see this kind of method in operation, e.g when companies send free new products to potential customers in order to convince them to adopt the product and influence others. We study the CO model, using two different strategies, on both ER and scale-free networks. In strategy I, the contrarians are positioned at random. In strategy II, the contrarians are chosen to be the highest degrees nodes. We find that for both strategies the size of the largest A cluster decreases to zero as p increases as in a phase transition. At a critical threshold value p_c the system undergoes a second-order phase transition that belongs to the same universality class of mean field percolation. We find that even for an ER type model, where the degrees of the nodes are not so distinct, strategy II is significantly more effctive in reducing the size of the largest A opinion cluster and, at very small values of p, the largest A opinion cluster is destroyed.

preprint2011arXiv

Suppressing Epidemics with a Limited Amount of Immunization Units

The way diseases spread through schools, epidemics through countries, and viruses through the Internet is crucial in determining their risk. Although each of these threats has its own characteristics, its underlying network determines the spreading. To restrain the spreading, a widely used approach is the fragmentation of these networks through immunization, so that epidemics cannot spread. Here we develop an immunization approach based on optimizing the susceptible size, which outperforms the best known strategy based on immunizing the highest-betweenness links or nodes. We find that the network's vulnerability can be significantly reduced, demonstrating this on three different real networks: the global flight network, a school friendship network, and the internet. In all cases, we find that not only is the average infection probability significantly suppressed, but also for the most relevant case of a small and limited number of immunization units the infection probability can be reduced by up to 55%.

preprint2011arXiv

The combined effect of connectivity and dependency links on percolation of networks

Percolation theory is extensively studied in statistical physics and mathematics with applications in diverse fields. However, the research is focused on systems with only one type of links, connectivity links. We review a recently developed mathematical framework for analyzing percolation properties of realistic scenarios of networks having links of two types, connectivity and dependency links. This formalism was applied to study Erd$\ddot{o}$s-R$\acute{e}$nyi (ER) networks that include also dependency links. For an ER network with average degree $k$ that is composed of dependency clusters of size $s$, the fraction of nodes that belong to the giant component, $P_\infty$, is given by $ P_\infty=p^{s-1}[1-\exp{(-kpP_\infty)}]^s $ where $1-p$ is the initial fraction of randomly removed nodes. Here, we apply the formalism to the study of random-regular (RR) networks and find a formula for the size of the giant component in the percolation process: $P_\infty=p^{s-1}(1-r^k)^s$ where $r$ is the solution of $r=p^s(r^{k-1}-1)(1-r^k)+1$. These general results coincide, for $s=1$, with the known equations for percolation in ER and RR networks respectively without dependency links. In contrast to $s=1$, where the percolation transition is second order, for $s>1$ it is of first order. Comparing the percolation behavior of ER and RR networks we find a remarkable difference regarding their resilience. We show, analytically and numerically, that in ER networks with low connectivity degree or large dependency clusters, removal of even a finite number (zero fraction) of the network nodes will trigger a cascade of failures that fragments the whole network. This result is in contrast to RR networks where such cascades and full fragmentation can be triggered only by removal of a finite fraction of nodes in the network.

preprint2010arXiv

Cascade of failures in coupled network systems with multiple support-dependent relations

We study, both analytically and numerically, the cascade of failures in two coupled network systems A and B, where multiple support-dependent relations are randomly built between nodes of networks A and B. In our model we assume that each node in one network can function only if it has at least a single support node in the other network. If both networks A and B are Erdős-Rényi networks, A and B, with (i) sizes $N^A$ and $N^B$, (ii) average degrees $a$ and $b$, and (iii) $c^{AB}_0N^B$ support links from network A to B and $c^{BA}_0N^B$ support links from network B to A, we find that under random attack with removal of fractions $(1-R^A)N^A$ and $(1-R^B)N^B$ nodes respectively, the percolating giant components of both networks at the end of the cascading failures, $μ^A_\infty$ and $μ^B_\infty$, are given by the percolation laws $μ^A_\infty = R^A [1-\exp{({-c^{BA}_0μ^B_\infty})}] [1-\exp{({-aμ^A_\infty})}]$ and $μ^B_\infty = R^B [1-\exp{({-c^{AB}_0μ^A_\infty})}] [1-\exp{({-bμ^B_\infty})}]$. In the limit of $c^{BA}_0 \to \infty$ and $c^{AB}_0 \to \infty$, both networks become independent, and the giant components are equivalent to a random attack on a single Erdős-Rényi network. We also test our theory on two coupled scale-free networks, and find good agreement with the simulations.

preprint2010arXiv

Epidemic threshold for the SIS model on networks

We derive an analytical expression for the critical infection rate r_c of the susceptible-infectious-susceptible (SIS) disease spreading model on random networks. To obtain r_c, we first calculate the probability of reinfection, pi, defined as the probability of a node to reinfect the node that had earlier infected it. We then derive r_c from pi using percolation theory. We show that pi is governed by two effects: (i) The requirement from an infecting node to recover prior to its reinfection, which depends on the disease spreading parameters; and (ii) The competition between nodes that simultaneously try to reinfect the same ancestor, which depends on the network topology.

preprint2010arXiv

Inter-similarity between coupled networks

Recent studies have shown that a system composed from several randomly interdependent networks is extremely vulnerable to random failure. However, real interdependent networks are usually not randomly interdependent, rather a pair of dependent nodes are coupled according to some regularity which we coin inter-similarity. For example, we study a system composed from an interdependent world wide port network and a world wide airport network and show that well connected ports tend to couple with well connected airports. We introduce two quantities for measuring the level of inter-similarity between networks (i) Inter degree-degree correlation (IDDC) (ii) Inter-clustering coefficient (ICC). We then show both by simulation models and by analyzing the port-airport system that as the networks become more inter-similar the system becomes significantly more robust to random failure.

preprint2010arXiv

Interdependent networks: Reducing the coupling strength leads to a change from a first to second order percolation transition

We study a system composed from two interdependent networks A and B, where a fraction of the nodes in network A depends on the nodes of network B and a fraction of the nodes in network B depends on the nodes of network A. Due to the coupling between the networks when nodes in one network fail they cause dependent nodes in the other network to also fail. This invokes an iterative cascade of failures in both networks. When a critical fraction of nodes fail the iterative process results in a percolation phase transition that completely fragments both networks. We show both analytically and numerically that reducing the coupling between the networks leads to a change from a first order percolation phase transition to a second order percolation transition at a critical point. The scaling of the percolation order parameter near the critical point is characterized by the critical exponent beta=1.

preprint2010arXiv

Market dynamics immediately before and after financial shocks: quantifying the Omori, productivity and Bath laws

We study the cascading dynamics immediately before and immediately after 219 market shocks. We define the time of a market shock T_{c} to be the time for which the market volatility V(T_{c}) has a peak that exceeds a predetermined threshold. The cascade of high volatility "aftershocks" triggered by the "main shock" is quantitatively similar to earthquakes and solar flares, which have been described by three empirical laws --- the Omori law, the productivity law, and the Bath law. We analyze the most traded 531 stocks in U.S. markets during the two-year period 2001-2002 at the 1-minute time resolution. We find quantitative relations between (i) the "main shock" magnitude M \equiv \log V(T_{c}) occurring at the time T_{c} of each of the 219 "volatility quakes" analyzed, and (ii) the parameters quantifying the decay of volatility aftershocks as well as the volatility preshocks. We also find that stocks with larger trading activity react more strongly and more quickly to market shocks than stocks with smaller trading activity. Our findings characterize the typical volatility response conditional on M, both at the market and the individual stock scale. We argue that there is potential utility in these three statistical quantitative relations with applications in option pricing and volatility trading.

preprint2010arXiv

Maximizing Entropy Yields Spatial Scaling in Social Networks

In addition to the well known common properties such as small world and community structures, recent empirical investigations suggest a universal scaling law for the spatial structure of social networks. It is found that the probability density distribution of an individual to have a friend at distance $r$ scales as $P(r)\propto r^{-1}$. The basic principle that yields this spatial scaling property is not yet understood. Here we propose a fundamental origin for this law based on the concept of entropy. We show that this spatial scaling law can result from maximization of information entropy, which means individuals seek to maximize the diversity of their friendships. Such spatial distribution can benefit individuals significantly in optimally collecting information in a social network.

preprint2010arXiv

Quantitative law describing market dynamics before and after interest-rate change

We study the behavior of U.S. markets both before and after U.S. Federal Open Market Committee (FOMC) meetings, and show that the announcement of a U.S. Federal Reserve rate change causes a financial shock, where the dynamics after the announcement is described by an analogue of the Omori earthquake law. We quantify the rate n(t) of aftershocks following an interest rate change at time T, and find power-law decay which scales as n(t-T) (t-T)^-$Ω$, with $Ω$ positive. Surprisingly, we find that the same law describes the rate n'(|t-T|) of "pre-shocks" before the interest rate change at time T. This is the first study to quantitatively relate the size of the market response to the news which caused the shock and to uncover the presence of quantifiable preshocks. We demonstrate that the news associated with interest rate change is responsible for causing both the anticipation before the announcement and the surprise after the announcement. We estimate the magnitude of financial news using the relative difference between the U. S. Treasury Bill and the Federal Funds Effective rate. Our results are consistent with the "sign effect," in which "bad news" has a larger impact than "good news." Furthermore, we observe significant volatility aftershocks, confirming a "market underreaction" that lasts at least 1 trading day.

preprint2010arXiv

Robustness of a Network of Networks

Almost all network research has been focused on the properties of a single network that does not interact and depends on other networks. In reality, many real-world networks interact with other networks. Here we develop an analytical framework for studying interacting networks and present an exact percolation law for a network of $n$ interdependent networks. In particular, we find that for $n$ Erdős-Rényi networks each of average degree $k$, the giant component, $P_{\infty}$, is given by $P_{\infty}=p[1-\exp(-kP_{\infty})]^n$ where $1-p$ is the initial fraction of removed nodes. Our general result coincides for $n=1$ with the known Erdős-Rényi second-order phase transition for a single network. For any $n \geq 2$ cascading failures occur and the transition becomes a first-order percolation transition. The new law for $P_{\infty}$ shows that percolation theory that is extensively studied in physics and mathematics is a limiting case ($n=1$) of a more general general and different percolation law for interdependent networks.

preprint2010arXiv

Robustness of interdependent networks under targeted attack

When an initial failure of nodes occurs in interdependent networks, a cascade of failure between the networks occurs. Earlier studies focused on random initial failures. Here we study the robustness of interdependent networks under targeted attack on high or low degree nodes. We introduce a general technique and show that the {\it targeted-attack} problem in interdependent networks can be mapped to the {\it random-attack} problem in a transformed pair of interdependent networks. We find that when the highly connected nodes are protected and have lower probability to fail, in contrast to single scale free (SF) networks where the percolation threshold $p_c=0$, coupled SF networks are significantly more vulnerable with $p_c$ significantly larger than zero. The result implies that interdependent networks are difficult to defend by strategies such as protecting the high degree nodes that have been found useful to significantly improve robustness of single networks.

preprint2010arXiv

The critical effect of dependency groups on the function of networks

Current network models assume one type of links to define the relations between the network entities. However, many real networks can only be correctly described using two different types of relations. Connectivity links that enable the nodes to function cooperatively as a network and dependency links that bind the failure of one network element to the failure of other network elements. Here we present for the first time an analytical framework for studying the robustness of networks that include both connectivity and dependency links. We show that the synergy between the two types of failures leads to an iterative process of cascading failures that has a devastating effect on the network stability and completely alters the known assumptions regarding the robustness of networks. We present exact analytical results for the dramatic change in the network behavior when introducing dependency links. For a high density of dependency links the network disintegrates in a form of a first order phase transition while for a low density of dependency links the network disintegrates in a second order transition. Moreover, opposed to networks containing only connectivity links where a broader degree distribution results in a more robust network, when both types of links are present a broad degree distribution leads to higher vulnerability.

preprint2010arXiv

The Emergence of El-Niño as an Autonomous Component in the Climate Network

We construct and analyze a climate network which represents the interdependent structure of the climate in different geographical zones and find that the network responds in a unique way to El-Niño events. Analyzing the dynamics of the climate network shows that when El-Niño events begin, the El-Niño basin partially loses its influence on its surroundings. After typically three months, this influence is restored while the basin loses almost all dependence on its surroundings and becomes \textit{autonomous}. The formation of an autonomous basin is the missing link to understand the seemingly contradicting phenomena of the afore--noticed weakening of the interdependencies in the climate network during El-Niño and the known impact of the anomalies inside the El-Niño basin on the global climate system.

preprint2010arXiv

Worldwide spreading of economic crisis

We model the spreading of a crisis by constructing a global economic network and applying the Susceptible-Infected-Recovered (SIR) epidemic model with a variable probability of infection. The probability of infection depends on the strength of economic relations between the pair of countries, and the strength of the target country. It is expected that a crisis which originates in a large country, such as the USA, has the potential to spread globally, like the recent crisis. Surprisingly we show that also countries with much lower GDP, such as Belgium, are able to initiate a global crisis. Using the {\it k}-shell decomposition method to quantify the spreading power (of a node), we obtain a measure of ``centrality'' as a spreader of each country in the economic network. We thus rank the different countries according to the shell they belong to, and find the 12 most central countries. These countries are the most likely to spread a crisis globally. Of these 12 only six are large economies, while the other six are medium/small ones, a result that could not have been otherwise anticipated. Furthermore, we use our model to predict the crisis spreading potential of countries belonging to different shells according to the crisis magnitude.

preprint2009arXiv

Catastrophic cascade of failures in interdependent networks

Many systems, ranging from engineering to medical to societal, can only be properly characterized by multiple interdependent networks whose normal functioning depends on one another. Failure of a fraction of nodes in one network may lead to a failure in another network. This in turn may cause further malfunction of additional nodes in the first network and so on. Such a cascade of failures, triggered by a failure of a small faction of nodes in only one network, may lead to the complete fragmentation of all networks. We introduce a model and an analytical framework for studying interdependent networks. We obtain interesting and surprising results that should significantly effect the design of robust real-world networks. For two interdependent Erdos-Renyi (ER) networks, we find that the critical average degree below which both networks collapse is <k_c>=2.445, compared to <k_c>=1 for a single ER network. Furthermore, while for a single network a broader degree distribution of the network nodes results in higher robustness to random failure, for interdependent networks, the broader the distribution is, the more vulnerable the networks become to random failure.

preprint2009arXiv

Dynamic networks and directed percolation

We introduce a model for dynamic networks, where the links or the strengths of the links change over time. We solve the model by mapping dynamic networks to the problem of directed percolation, where the direction corresponds to the evolution of the network in time. We show that the dynamic network undergoes a percolation phase transition at a critical concentration $p_c$, which decreases with the rate $r$ at which the network links are changed. The behavior near criticality is universal and independent of $r$. We find fundamental network laws are changed. (i) For Erdős-Rényi networks we find that the size of the giant component at criticality scales with the network size $N$ for all values of $r$, rather than as $N^{2/3}$. (ii) In the presence of a broad distribution of disorder, the optimal path length between two nodes in a dynamic network scales as $N^{1/2}$, compared to $N^{1/3}$ in a static network.

preprint2009arXiv

Scaling laws of human interaction activity

Even though people in our contemporary, technological society are depending on communication, our understanding of the underlying laws of human communicational behavior continues to be poorly understood. Here we investigate the communication patterns in two social Internet communities in search of statistical laws in human interaction activity. This research reveals that human communication networks dynamically follow scaling laws that may also explain the observed trends in economic growth. Specifically, we identify a generalized version of Gibrat's law of social activity expressed as a scaling law between the fluctuations in the number of messages sent by members and their level of activity. Gibrat's law has been essential in understanding economic growth patterns, yet without an underlying general principle for its origin. We attribute this scaling law to long-term correlation patterns in human activity, which surprisingly span from days to the entire period of the available data of more than one year. Further, we provide a mathematical framework that relates the generalized version of Gibrat's law to the long-term correlated dynamics, which suggests that the same underlying mechanism could be the source of Gibrat's law in economics, ranging from large firms, research and development expenditures, gross domestic product of countries, to city population growth. These findings are also of importance for designing communication networks and for the understanding of the dynamics of social systems in which communication plays a role, such as economic markets and political systems.

preprint2009arXiv

Structure of shells in complex networks

In a network, we define shell $\ell$ as the set of nodes at distance $\ell$ with respect to a given node and define $r_\ell$ as the fraction of nodes outside shell $\ell$. In a transport process, information or disease usually diffuses from a random node and reach nodes shell after shell. Thus, understanding the shell structure is crucial for the study of the transport property of networks. For a randomly connected network with given degree distribution, we derive analytically the degree distribution and average degree of the nodes residing outside shell $\ell$ as a function of $r_\ell$. Further, we find that $r_\ell$ follows an iterative functional form $r_\ell=ϕ(r_{\ell-1})$, where $ϕ$ is expressed in terms of the generating function of the original degree distribution of the network. Our results can explain the power-law distribution of the number of nodes $B_\ell$ found in shells with $\ell$ larger than the network diameter $d$, which is the average distance between all pairs of nodes. For real world networks the theoretical prediction of $r_\ell$ deviates from the empirical $r_\ell$. We introduce a network correlation function $c(r_\ell)\equiv r_{\ell+1}/ϕ(r_\ell)$ to characterize the correlations in the network, where $r_{\ell+1}$ is the empirical value and $ϕ(r_\ell)$ is the theoretical prediction. $c(r_\ell)=1$ indicates perfect agreement between empirical results and theory. We apply $c(r_\ell)$ to several model and real world networks. We find that the networks fall into two distinct classes: (i) a class of {\it poorly-connected} networks with $c(r_\ell)>1$, which have larger average distances compared with randomly connected networks with the same degree distributions; and (ii) a class of {\it well-connected} networks with $c(r_\ell)<1$.

preprint2008arXiv

Climate Networks around the Globe are Significantly Effected by El Nino

The temperatures in different zones in the world do not show significant changes due to El-Nino except when measured in a restricted area in the Pacific Ocean. We find, in contrast, that the dynamics of a climate network based on the same temperature records in various geographical zones in the world is significantly influenced by El-Nino. During El-Nino many links of the network are broken, and the number of surviving links comprises a specific and sensitive measure for El-Nino events. While during non El-Nino periods these links which represent correlations between temperatures in different sites are more stable, fast fluctuations of the correlations observed during El-Nino periods cause the links to break.

preprint2008arXiv

Fractal Boundaries of Complex Networks

We introduce the concept of boundaries of a complex network as the set of nodes at distance larger than the mean distance from a given node in the network. We study the statistical properties of the boundaries nodes of complex networks. We find that for both Erdös-Rényi and scale-free model networks, as well as for several real networks, the boundaries have fractal properties. In particular, the number of boundaries nodes {\it B} follows a power-law probability density function which scales as $B^{-2}$. The clusters formed by the boundary nodes are fractals with a fractal dimension $d_{f} \approx 2$. We present analytical and numerical evidence supporting these results for a broad class of networks. Our findings imply potential applications for epidemic spreading.

preprint2006arXiv

Anomalous electrical and frictionless flow conductance in complex networks

We study transport properties such as electrical and frictionless flow conductance on scale-free and Erdos-Renyi networks. We consider the conductance G between two arbitrarily chosen nodes where each link has the same unit resistance. Our theoretical analysis for scale-free networks predicts a broad range of values of G, with a power-law tail distribution Φ_{SF}(G) \sim G^{g_G}, where g_G = 2λ- 1, where λis the decay exponent for the scale-free network degree distribution. We confirm our predictions by simulations of scale-free networks solving the Kirchhoff equations for the conductance between a pair of nodes. The power-law tail in Φ_{SF}(G) leads to large values of G, thereby significantly improving the transport in scale-free networks, compared to Erdos-Renyi networks where the tail of the conductivity distribution decays exponentially. Based on a simple physical 'transport backbone' picture we suggest that the conductances of scale-free and Erdos-Renyi networks can be approximated by ck_Ak_B/(k_A+k_B) for any pair of nodes A and B with degrees k_A and k_B. Thus, a single quantity c, which depends on the average degree <k> of the network, characterizes transport on both scale-free and Erdos-Renyi networks. We determine that c tends to 1 for increasing <k>, and it is larger for scale-free networks. We compare the electrical results with a model for frictionless transport, where conductance is defined as the number of link-independent paths between A and B, and find that a similar picture holds. The effects of distance on the value of conductance are considered for both models, and some differences emerge. Finally, we use a recent data set for the AS (autonomous system) level of the Internet and confirm that our results are valid in this real-world example.

preprint2006arXiv

Transport of multiple users in complex networks

We study the transport properties of model networks such as scale-free and Erdős-Rényi networks as well as a real network. We consider the conductance $G$ between two arbitrarily chosen nodes where each link has the same unit resistance. Our theoretical analysis for scale-free networks predicts a broad range of values of $G$, with a power-law tail distribution $Φ_{\rm SF}(G)\sim G^{-g_G}$, where $g_G=2λ-1$, and $λ$ is the decay exponent for the scale-free network degree distribution. We confirm our predictions by large scale simulations. The power-law tail in $Φ_{\rm SF}(G)$ leads to large values of $G$, thereby significantly improving the transport in scale-free networks, compared to Erdős-Rényi networks where the tail of the conductivity distribution decays exponentially. We develop a simple physical picture of the transport to account for the results. We study another model for transport, the \emph{max-flow} model, where conductance is defined as the number of link-independent paths between the two nodes, and find that a similar picture holds. The effects of distance on the value of conductance are considered for both models, and some differences emerge. We then extend our study to the case of multiple sources, where the transport is define between two \emph{groups} of nodes. We find a fundamental difference between the two forms of flow when considering the quality of the transport with respect to the number of sources, and find an optimal number of sources, or users, for the max-flow case. A qualitative (and partially quantitative) explanation is also given.

preprint2006arXiv

Universal behavior of optimal paths in weighted networks with general disorder

We study the statistics of the optimal path in both random and scale free networks, where weights $w$ are taken from a general distribution $P(w)$. We find that different types of disorder lead to the same universal behavior. Specifically, we find that a single parameter ($S \equiv AL^{-1/ν}$ for $d$-dimensional lattices, and $S\equiv AN^{-1/3}$ for random networks) determines the distributions of the optimal path length, including both strong and weak disorder regimes. Here $ν$ is the percolation connectivity exponent, and $A$ depends on the percolation threshold and $P(w)$. For $P(w)$ uniform, Poisson or Gaussian the crossover from weak to strong does not occur, and only weak disorder exists.

preprint2005arXiv

Current Flow in Random Resistor Networks: The Role of Percolation in Weak and Strong Disorder

We study the current flow paths between two edges in a random resistor network on a $L\times L$ square lattice. Each resistor has resistance $e^{ax}$, where $x$ is a uniformly-distributed random variable and $a$ controls the broadness of the distribution. We find (a) the scaled variable $u\equiv L/a^ν$, where $ν$ is the percolation connectedness exponent, fully determines the distribution of the current path length $\ell$ for all values of $u$. For $u\gg 1$, the behavior corresponds to the weak disorder limit and $\ell$ scales as $\ell\sim L$, while for $u\ll 1$, the behavior corresponds to the strong disorder limit with $\ell\sim L^{d_{\scriptsize opt}}$, where $d_{\scriptsize opt} = 1.22\pm0.01$ is the optimal path exponent. (b) In the weak disorder regime, there is a length scale $ξ\sim a^ν$, below which strong disorder and critical percolation characterize the current path.

preprint2005arXiv

Possible Connection between the Optimal Path and Flow in Percolation Clusters

We study the behavior of the optimal path between two sites separated by a distance $r$ on a $d$-dimensional lattice of linear size $L$ with weight assigned to each site. We focus on the strong disorder limit, i.e., when the weight of a single site dominates the sum of the weights along each path. We calculate the probability distribution $P(\ell_{\rm opt}|r,L)$ of the optimal path length $\ell_{\rm opt}$, and find for $r\ll L$ a power law decay with $\ell_{\rm opt}$, characterized by exponent $g_{\rm opt}$. We determine the scaling form of $P(\ell_{\rm opt}|r,L)$ in two- and three-dimensional lattices. To test the conjecture that the optimal paths in strong disorder and flow in percolation clusters belong to the same universality class, we study the tracer path length $\ell_{\rm tr}$ of tracers inside percolation through their probability distribution $P(\ell_{\rm tr}|r,L)$. We find that, because the optimal path is not constrained to belong to a percolation cluster, the two problems are different. However, by constraining the optimal paths to remain inside the percolation clusters in analogy to tracers in percolation, the two problems exhibit similar scaling properties.

preprint2005arXiv

Scale-Free Networks Emerging from Weighted Random Graphs

We study Erdös-Rényi random graphs with random weights associated with each link. We generate a new ``Supernode network'' by merging all nodes connected by links having weights below the percolation threshold (percolation clusters) into a single node. We show that this network is scale-free, i.e., the degree distribution is $P(k)\sim k^{-λ}$ with $λ=2.5$. Our results imply that the minimum spanning tree (MST) in random graphs is composed of percolation clusters, which are interconnected by a set of links that create a scale-free tree with $λ=2.5$. We show that optimization causes the percolation threshold to emerge spontaneously, thus creating naturally a scale-free ``supernode network''. We discuss the possibility that this phenomenon is related to the evolution of several real world scale-free networks.

preprint2004arXiv

Anomalous Transport in Complex Networks

To study transport properties of complex networks, we analyze the equivalent conductance $G$ between two arbitrarily chosen nodes of random scale-free networks with degree distribution $P(k)\sim k^{-λ}$ in which each link has the same unit resistance. We predict a broad range of values of $G$, with a power-law tail distribution $Φ_{\rm SF}(G)\sim G^{-g_G}$, where $g_G=2λ-1$, and confirm our predictions by simulations. The power-law tail in $Φ_{\rm SF}(G)$ leads to large values of $G$, thereby significantly improving the transport in scale-free networks, compared to Erdős-Rényi random graphs where the tail of the conductivity distribution decays exponentially. Based on a simple physical ``transport backbone'' picture we show that the conductances are well approximated by $ck_Ak_B/(k_A+k_B)$ for any pair of nodes $A$ and $B$ with degrees $k_A$ and $k_B$. Thus, a single parameter $c$ characterizes transport on scale-free networks.

preprint2004arXiv

Volcanic forcing improves Atmosphere-Ocean Coupled General Circulation Model scaling performance

Recent Atmosphere-Ocean Coupled General Circulation Model (AOGCM) simulations of the twentieth century climate, which account for anthropogenic and natural forcings, make it possible to study the origin of long-term temperature correlations found in the observed records. We study ensemble experiments performed with the NCAR PCM for 10 different historical scenarios, including no forcings, greenhouse gas, sulfate aerosol, ozone, solar, volcanic forcing and various combinations, such as it natural, anthropogenic and all forcings. We compare the scaling exponents characterizing the long-term correlations of the observed and simulated model data for 16 representative land stations and 16 sites in the Atlantic Ocean for these scenarios. We find that inclusion of volcanic forcing in the AOGCM considerably improves the PCM scaling behavior. The scenarios containing volcanic forcing are able to reproduce quite well the observed scaling exponents for the land with exponents around 0.65 independent of the station distance from the ocean. For the Atlantic Ocean, scenarios with the volcanic forcing slightly underestimate the observed persistence exhibiting an average exponent 0.74 instead of 0.85 for reconstructed data.

preprint2003arXiv

Multifractality of river runoff and precipitation: Comparison of fluctuation analysis and wavelet methods

We study the multifractal temporal scaling properties of river discharge and precipitation records. We compare the results for the multifractal detrended fluctuation analysis method with the results for the wavelet transform modulus maxima technique and obtain agreement within the error margins. In contrast to previous studies, we find non-universal behaviour: On long time scales, above a crossover time scale of several months, the runoff records are described by fluctuation exponents varying from river to river in a wide range. Similar variations are observed for the precipitation records which exhibit weaker, but still significant multifractality. For all runoff records the type of multifractality is consistent with a modified version of the binomial multifractal model, while several precipitation records seem to require different models.

preprint2002arXiv

Nonlinearity and Multifractality of Climate Change in the Past 420,000 Years

Evidence of past climate variations are stored in ice and indicate glacial-interglacial cycles characterized by three dominant time periods of 20kyr, 40kyr, and 100kyr. We study the scaling properties of temperature proxy records of four ice cores from Antarctica and Greenland. These series are long-range correlated in the time scales of 1-100kyr. We show that these series are nonlinear as expressed by volatility correlations and a broad multifractal spectrum. We present a stochastic model that captures the scaling and the nonlinear properties observed in the data.

preprint2000arXiv

Universality of the Directed Polymer Model

The universality of the directed polymer model and the analogous KPZ equation is supported by numerical simulations using non-Gaussian random probability distributions in two, three and four dimensions. It is shown that although in the non-Gaussian cases the \emph{finite size} estimates of the energy exponents are below the persumed universal values, these estimates \emph{increase} with the system size, and the further they are below the universal values, the higher is their rate of increase. The results are explained in terms of the efficiency of variance reduction during the optimization process.

preprint1999arXiv

Directed polymers at finite temperatures in 1+1 and 2+1 dimensions

We present systematic numerical simulations for directed polymers at finite temperatures in 1+1 and 2+1 dimensions. The transverse fluctuations and free energy fluctuations tend to the strong coupling limit at any temperature in both 1+1 and 2+1 dimensions for long time t. Two different definitions for energy fluctuations at finite temperatures, which are the ensemble energy fluctuations and the internal energy fluctuations, are investigated. Apart from zero temperature, the behavior of the energy fluctuations and the free energy fluctuations for directed polymers is shown to be different. At finite temperatures, the ensemble energy fluctuations in both 1+1 and 2+1 dimensions and internal energy fluctuations in 1+1 dimensions scale as t^{1/2} where the free energy fluctuations in 1+1 dimensions and 2+1 dimensions scale as t^{1/3} and t^{0.2} respectively. As a consequence of that the specific heat in both 1+1 and 2+1 dimensions scales as t and the entropy fluctuations in 1+1 dimensions scale as t^{1/2} at any finite temperature.

preprint1999arXiv

Distribution of Dangling Ends on the Incipient Percolation Cluster

We study numerically and by scaling arguments the probability P(M)dM that a given dangling end of the incipient percolation cluster has a mass between M and M + dM. We find by scaling arguments that P(M) decays with a power law, P(M)~M^(-(1+k)), with an exponent k=dBf/df, where df and dBf are the fractal dimensions of the cluster and its backbone, respectively. Our numerical results yield k=0.83 in d=2 and k=0.74 in d=3 in very good agreement with theory.

preprint1999arXiv

Multifractality in Human Heartbeat Dynamics

Recent evidence suggests that physiological signals under healthy conditions may have a fractal temporal structure. We investigate the possibility that time series generated by certain physiological control systems may be members of a special class of complex processes, termed multifractal, which require a large number of exponents to characterize their scaling properties. We report on evidence for multifractality in a biological dynamical system --- the healthy human heartbeat. Further, we show that the multifractal character and nonlinear properties of the healthy heart rate are encoded in the Fourier phases. We uncover a loss of multifractality for a life-threatening condition, congestive heart failure.

preprint1999arXiv

Sleep-Wake Differences in Scaling Behavior of the Human Heartbeat: Analysis of Terrestrial and Long-Term Space Flight Data

We compare scaling properties of the cardiac dynamics during sleep and wake periods for healthy individuals, cosmonauts during orbital flight, and subjects with severe heart disease. For all three groups, we find a greater degree of anticorrelation in the heartbeat fluctuations during sleep compared to wake periods. The sleep-wake difference in the scaling exponents for the three groups is comparable to the difference between healthy and diseased individuals. The observed scaling differences are not accounted for simply by different levels of activity, but appear related to intrinsic changes in the neuroautonomic control of the heartbeat.

preprint1998arXiv

Spontaneous Stratification in Granular Mixtures

Granular materials size segregate when exposed to external periodic perturbations such as vibrations. Moreover, mixtures of grains of different sizes spontaneously segregate in the absence of external perturbations: when a mixture is simply poured onto a pile, the large grains are more likely to be found near the base, while the small grains are more likely to be near the top. Here, we report a spontaneous phenomenon arising when we pour a mixture between two vertical plates: the mixture spontaneously stratifies into alternating layers of small and large grains whenever the large grains are rougher than the small grains. In contrast, we find only spontaneous segregation when the large grains are more rounded than the small grains. The stratification is related to the occurrence of avalanches; during each avalanche the grains comprising the avalanche spontaneously stratify into a pair of layers through a "kink" mechanism, with the small grains forming a sublayer underneath the layer of large grains.