Disentangling different types of El Niño episodes by evolving climate network analysis

Reconfiguration of Amazon's connectivity in the climate system

With the recent increase in deforestation, forest fires, and regional temperatures, the concerns around the rapid and complete collapse of the Amazon rainforest ecosystem have heightened. The thresholds of deforestation and the temperature increase required for such a catastrophic event are still uncertain. However, our analysis presented here shows that signatures of changing Amazon are already apparent in historical climate data sets. Here, we extend the methods of climate network analysis and apply them to study the temporal evolution of the connectivity between the Amazon rainforest and the global climate system. We observe that the Amazon rainforest is losing short-range connectivity and gaining more long-range connections, indicating shifts in regional-scale processes. Using embeddings inspired by manifold learning, we show that the Amazon connectivity patterns have undergone a fundamental shift in the 21st century. By investigating edge-based network metrics on similar regions to the Amazon, we see the changing properties of the Amazon are noticeable in comparison. Furthermore, we simulate diffusion and random walks on these networks and observe a faster spread of perturbations from the Amazon in recent decades. Our methodology innovations can act as a template for examining the spatiotemporal patterns of regional climate change and its impact on global climate using the toolbox of climate network analysis.

Network percolation provides early warnings of abrupt changes in coupled oscillatory systems: An explanatory analysis

Functional networks are powerful tools to study statistical interdependency structures in spatially extended or multivariable systems. They have been used to get insights into the dynamics of complex systems in various areas of science. In particular, percolation properties of correlation networks have been employed to identify early warning signals of critical transitions. In this work, we further investigate the corresponding potential of percolation measures for the anticipation of different types of sudden shifts in the state of coupled irregularly oscillating systems. As a paradigmatic model system, we study the dynamics of a ring of diffusively coupled noisy FitzHugh-Nagumo oscillators and show that, when the oscillators are nearly completely synchronized, the percolation-based precursors successfully provide very early warnings of the rapid switches between the two states of the system. We clarify the mechanisms behind the percolation transition by separating global trends given by the mean-field behavior from the synchronization of individual stochastic fluctuations. We then apply the same methodology to real-world data of sea surface temperature anomalies during different phases of the El Niño-Southern Oscillation. This leads to a better understanding of the factors that make percolation precursors effective as early warning indicators of incipient El Niño and La Niña events.

Early warning of the Indian Ocean Dipole using climate network analysis

The Indian Ocean Dipole (IOD), an air–sea coupled phenomenon over the tropical Indian Ocean, has substantial impacts on the climate, ecosystems, and society. Due to the winter predictability barrier, however, a reliable prediction of the IOD has been limited to 3 or 4 mo in advance. Our work approaches this problem from a new data-driven perspective: the climate network analysis. Using this network-based method, an efficient early warning signal for the IOD event was revealed in boreal winter. Our approach can correctly predict the IOD events one calendar year in advance (from December of the previous year) with a hit rate of higher than 70%, which strongly outperforms current dynamic models.

In recent years, the Indian Ocean Dipole (IOD) has received much attention in light of its substantial impacts on both the climate system and humanity. Due to its complexity, however, a reliable prediction of the IOD is still a great challenge. In this study, climate network analysis was employed to investigate whether there are early warning signals prior to the start of IOD events. An enhanced seesaw tendency in sea surface temperature (SST) among a large number of grid points between the dipole regions in the tropical Indian Ocean was revealed in boreal winter, which can be used to forewarn the potential occurrence of the IOD in the coming year. We combined this insight with the indicator of the December equatorial zonal wind in the tropical Indian Ocean to propose a network-based predictor that clearly outperforms the current dynamic models. Of the 15 IOD events over the past 37 y (1984 to 2020), 11 events were correctly predicted from December of the previous year, i.e., a hit rate of higher than 70%, and the false alarm rate was around 35%. This network-based approach suggests a perspective for better understanding and predicting the IOD.

A Novel Information Theoretical Criterion for Climate Network Construction

This paper presents a novel methodology for Climate Network (CN) construction based on the Kullback-Leibler divergence (KLD) among Membership Probability (MP) distributions, obtained from the Second Order Data-Coupled Clustering (SODCC) algorithm. The proposed method is able to obtain CNs with emergent behaviour adapted to the variables being analyzed, and with a low number of spurious or missing links. We evaluate the proposed method in a problem of CN construction to assess differences in wind speed prediction at different wind farms in Spain. The considered problem presents strong local and mesoscale relationships, but low synoptic scale relationships, which have a direct influence in the CN obtained. We carry out a comparison of the proposed approach with a classical correlation-based CN construction method. We show that the proposed approach based on the SODCC algorithm and the KLD constructs CNs with an emergent behaviour according to underlying wind speed prediction data physics, unlike the correlation-based method that produces spurious and missing links. Furthermore, it is shown that the climate network construction method facilitates the evaluation of symmetry properties in the resulting complex networks.

Event synchrony measures for functional climate network analysis: A case study on South American rainfall dynamics

Understanding spatiotemporal patterns of climate extremes has gained considerable relevance in the context of ongoing climate change. With enhanced computational capacity, data driven methods such as functional climate networks have been proposed and have already contributed to significant advances in understanding and predicting extreme events, as well as identifying interrelations between the occurrences of various climatic phenomena. While the (in its basic setting) parameter free event synchronization (ES) method has been widely applied to construct functional climate networks from extreme event series, its original definition has been realized to exhibit problems in handling events occurring at subsequent time steps, which need to be accounted for. Along with the study of this conceptual limitation of the original ES approach, event coincidence analysis (ECA) has been suggested as an alternative approach that incorporates an additional parameter for selecting certain time scales of event synchrony. In this work, we compare selected features of functional climate network representations of South American heavy precipitation events obtained using ES and ECA without and with the correction for temporal event clustering. We find that both measures exhibit different types of biases, which have profound impacts on the resulting network structures. By combining the complementary information captured by ES and ECA, we revisit the spatiotemporal organization of extreme events during the South American Monsoon season. While the corrected version of ES captures multiple time scales of heavy rainfall cascades at once, ECA allows disentangling those scales and thereby tracing the spatiotemporal propagation more explicitly.

Areawise significance tests for windowed recurrence network analysis

Many time-series analysis techniques use sliding window approaches or are repeatedly applied over a continuous range of parameters. When combined with a significance test, intrinsic correlations among the pointwise analysis results can make falsely positive significant points appear as continuous patches rather than as isolated points. To account for this effect, we present an areawise significance test that identifies such false-positive patches. For this purpose, we numerically estimate the decorrelation length of the statistic of interest by calculating correlation functions between the analysis results and require an areawise significant point to belong to a patch of pointwise significant points that is larger than this decorrelation length. We apply our areawise test to results from windowed traditional and scale-specific recurrence network analysis in order to identify dynamical anomalies in time series of a non-stationary Rössler system and tree ring width index values from Eastern Canada. Especially, in the palaeoclimate context, the areawise testing approach markedly reduces the number of points that are identified as significant and therefore highlights only the most relevant features in the data. This provides a crucial step towards further establishing recurrence networks as a tool for palaeoclimate data analysis.

Disentangling the multi-scale effects of sea-surface temperatures on global precipitation: A coupled networks approach

The oceans and atmosphere interact via a multiplicity of feedback mechanisms, shaping to a large extent the global climate and its variability. To deepen our knowledge of the global climate system, characterizing and investigating this interdependence is an important task of contemporary research. However, our present understanding of the underlying large-scale processes is greatly limited due to the manifold interactions between essential climatic variables at different temporal scales. To address this problem, we here propose to extend the application of complex network techniques to capture the interdependence between global fields of sea-surface temperature (SST) and precipitation (P) at multiple temporal scales. For this purpose, we combine time-scale decomposition by means of a discrete wavelet transform with the concept of coupled climate network analysis. Our results demonstrate the potential of the proposed approach to unravel the scale-specific interdependences between atmosphere and ocean and, thus, shed light on the emerging multiscale processes inherent to the climate system, which traditionally remain undiscovered when investigating the system only at the native resolution of existing climate data sets. Moreover, we show how the relevant spatial interdependence structures between SST and P evolve across time-scales. Most notably, the strongest mutual correlations between SST and P at annual scale (8-16 months) concentrate mainly over the Pacific Ocean, while the corresponding spatial patterns progressively disappear when moving toward longer time-scales.

Edge directionality properties in complex spherical networks

Spatially embedded networks have attracted increasing attention in the past decade. In this context, network characteristics have been introduced which explicitly take spatial information into account. Among others, edge directionality properties have recently gained particular interest. In this work, we investigate the applicability of mean edge direction, anisotropy, and local mean angle as geometric characteristics in complex spherical networks. By studying these measures, both analytically and numerically, we demonstrate the existence of a systematic bias in spatial networks where individual nodes represent different shares on a spherical surface, and we describe a strategy for correcting for this effect. Moreover, we illustrate the application of the mentioned edge directionality properties to different examples of real-world spatial networks in spherical geometry (with or without the geometric correction depending on each specific case), including functional climate networks, transportation, and trade networks. In climate networks, our approach highlights relevant patterns, such as large-scale circulation cells, the El Niño-Southern Oscillation, and the Atlantic Niño. In an air transportation network, we are able to characterize distinct air transportation zones, while we confirm the important role of the European Union for the global economy by identifying convergent edge directionality patterns in the world trade network.

Evaluation of ENSO simulations in CMIP5 models: A new perspective based on percolation phase transition in complex networks

In this study, the performance of CMIP5 models in simulating the El Niño-Southern Oscillation (ENSO) is evaluated by using a new metric based on percolation theory. The surface air temperatures (SATs) over the tropical Pacific Ocean are constructed as a SAT network, and the nodes within the network are linked if they are highly connected (e.g., high correlations). It has been confirmed from reanalysis datasets that the SAT network undergoes an abrupt percolation phase transition when the influences of the sea surface temperature anomalies (SSTAs) below are strong enough. However, from simulations of the CMIP5 models, most models are found incapable of capturing the observed phase transition at a proper critical point P_c. For the 15 considered models, four even miss the phase transition, indicating that the simulated SAT network is too stable to be significantly changed by the SSTA below. Only four models can be considered cautiously with some skills in simulating the observed phase transition of the SAT network. By comparing the simulated SSTA patterns with the node vulnerabilities, which is the chance of each node being isolated during a ENSO event, we find that the improperly simulated sea-air interactions are responsible for the missing of the observed percolation phase transition. Accordingly, a careful study of the sea-air couplers, as well as the atmospheric components of the CMIP5 models is suggested. Since the percolation phase transition of the SAT network is a useful phenomenon to indicate whether the ENSO impacts can be transferred remotely, it deserves more attention for future model development.

Constructing regional climate networks in the Amazonia during recent drought events

Climate networks are powerful approaches to disclose tele-connections in climate systems and to predict severe climate events. Here we construct regional climate networks from precipitation data in the Amazonian region and focus on network properties under the recent drought events in 2005 and 2010. Both the networks of the entire Amazon region and the extreme networks resulted from locations severely affected by drought events suggest that network characteristics show slight difference between the two drought events. Based on network degrees of extreme drought events and that without drought conditions, we identify regions of interest that are correlated to longer expected drought period length. Moreover, we show that the spatial correlation length to the regions of interest decayed much faster in 2010 than in 2005, which is because of the dual roles played by both the Pacific and Atlantic oceans. The results suggest that hub nodes in the regional climate network of Amazonia have fewer long-range connections when more severe drought conditions appeared in 2010 than that in 2005.

ENSO teleconnections in the southern hemisphere: A climate network view

Using functional network analysis, we study the seasonality of atmospheric connectivity and its interannual variability depending on the different phases of the El Niño-Southern Oscillation (ENSO) phenomenon. We find a strong variability of the connectivity on seasonal and interannual time scales both in the tropical and extratropical regions. In particular, there are significant changes in the southern hemisphere extratropical atmospheric connectivity during austral spring within the different stages of ENSO: We find that the connectivity patterns due to stationary Rossby waves differ during El Niño and La Niña, showing a very clear wave train originating close to Australia in the former case, as opposed to La Niña that seems to generate a wave train from the central Pacific. An attempt to understand these differences in terms of changes in the frequency of intraseasonal weather regimes cannot fully explain the differences in connectivity, even though we found the prevalence of different intraseasonal regimes in each phase of ENSO. We conclude that the differential response to extreme phases of ENSO during austral springtime is related to the forcing of waves of different tropical origins.

Percolation Phase Transition of Surface Air Temperature Networks: A new test bed for El Niño/La Niña simulations

In this work, we studied the air-sea interaction over the tropical central eastern Pacific from a new perspective, climate network. The surface air temperatures over the tropical Pacific were constructed as a network, and the nodes within this network were linked if they have a similar temporal varying pattern. Using three different reanalysis datasets, we verified the percolation phase transition. That is, when the influences of El Niño/La Niña are strong enough to isolate more than 48% of the nodes, the network may abruptly be divided into many small pieces, indicating a change of the network state. This phenomenon was reproduced successfully by a coupled general circulation model, Flexible Global Ocean-Atmosphere-Land System Model Spectral Version 2, but another model, Flexible Global Ocean-Atmosphere-Land System Model Grid-point Version 2, failed. As both models have the same oceanic component, but are with different atmospheric components, the improperly used atmospheric component should be responsible for the missing of the percolation phase transition. Considering that this new phenomenon is only recently noticed, current state-of-the-art models may ignore this process and induce unrealistic simulations. Accordingly, percolation phase transition is proposed as a new test bed, which deserves more attention in the future.

Network analysis reveals strongly localized impacts of El Niño

Climatic conditions influence the culture and economy of societies and the performance of economies. Specifically, El Niño as an extreme climate event is known to have notable effects on health, agriculture, industry, and conflict. Here, we construct directed and weighted climate networks based on near-surface air temperature to investigate the global impacts of El Niño and La Niña. We find that regions that are characterized by higher positive/negative network "in"-weighted links are exhibiting stronger correlations with the El Niño basin and are warmer/cooler during El Niño/La Niña periods. In contrast to non-El Niño periods, these stronger in-weighted activities are found to be concentrated in very localized areas, whereas a large fraction of the globe is not influenced by the events. The regions of localized activity vary from one El Niño (La Niña) event to another; still, some El Niño (La Niña) 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 Niño/La Niña events and also may be applied in the investigation of other climate variables.

The theory of parallel climate realizations as a new framework for teleconnection analysis

Teleconnections are striking features of the Earth climate system which appear as statistically correlated climate-related patterns between remote geographical regions of the globe. In a changing climate, however, the strength of teleconnections might change, and an appropriate characterization of these correlations and their change (more appropriate than detrending the time series) is lacking in the literature. Here we present a novel approach, based on the theory of snapshot attractors, corresponding in our context to studying parallel climate realizations. Imagining an ensemble of parallel Earth systems, instead of the single one observed (i.e., the real Earth), the ensemble, after some time, characterizes the appropriate probabilities of all options permitted by the climate dynamics, reflecting the internal variability of the climate. We claim that the relevant quantities for characterizing teleconnections in a changing climate are correlation coefficients taken over the temporally evolving ensemble in any time instant. As a particular example, we consider the teleconnections of the North Atlantic Oscillation (NAO). In a numerical climate model, we demonstrate that this approach provides the only statistically correct characterization, in contrast to commonly used temporal correlations evaluated along single detrended time series. The teleconnections of the NAO are found to survive the climate change, but their strength might be time-dependent.

Multiple-node basin stability in complex dynamical networks

Dynamical entities interacting with each other on complex networks often exhibit multistability. The stability of a desired steady regime (e.g., a synchronized state) to large perturbations is critical in the operation of many real-world networked dynamical systems such as ecosystems, power grids, the human brain, etc. This necessitates the development of appropriate quantifiers of stability of multiple stable states of such systems. Motivated by the concept of basin stability (BS) [P. J. Menck et al., Nat. Phys. 9, 89 (2013)1745-247310.1038/nphys2516], we propose here the general framework of multiple-node basin stability for gauging the global stability and robustness of networked dynamical systems in response to nonlocal perturbations simultaneously affecting multiple nodes of a system. The framework of multiple-node BS provides an estimate of the critical number of nodes that, when simultaneously perturbed, significantly reduce the capacity of the system to return to the desired stable state. Further, this methodology can be applied to estimate the minimum number of nodes of the network to be controlled or safeguarded from external perturbations to ensure proper operation of the system. Multiple-node BS can also be utilized for probing the influence of spatially localized perturbations or targeted attacks to specific parts of a network. We demonstrate the potential of multiple-node BS in assessing the stability of the synchronized state in a deterministic scale-free network of Rössler oscillators and a conceptual model of the power grid of the United Kingdom with second-order Kuramoto-type nodal dynamics.

Smooth information flow in temperature climate network reflects mass transport

A directed climate network is constructed by Granger causality analysis of air temperature time series from a regular grid covering the whole Earth. Using winner-takes-all network thresholding approach, a structure of a smooth information flow is revealed, hidden to previous studies. The relevance of this observation is confirmed by comparison with the air mass transfer defined by the wind field. Their close relation illustrates that although the information transferred due to the causal influence is not a physical quantity, the information transfer is tied to the transfer of mass and energy.

Spatio-temporal organization of dynamics in a two-dimensional periodically driven vortex flow: A Lagrangian flow network perspective

We study the Lagrangian dynamics of passive tracers in a simple model of a driven two-dimensional vortex resembling real-world geophysical flow patterns. Using a discrete approximation of the system's transfer operator, we construct a directed network that describes the exchange of mass between distinct regions of the flow domain. By studying different measures characterizing flow network connectivity at different time-scales, we are able to identify the location of dynamically invariant structures and regions of maximum dispersion. Specifically, our approach allows us to delimit co-existing flow regimes with different dynamics. To validate our findings, we compare several network characteristics to the well-established finite-time Lyapunov exponents and apply a receiver operating characteristic analysis to identify network measures that are particularly useful for unveiling the skeleton of Lagrangian chaos.

A perturbation-theoretic approach to Lagrangian flow networks

Complex network approaches have been successfully applied for studying transport processes in complex systems ranging from road, railway, or airline infrastructures over industrial manufacturing to fluid dynamics. Here, we utilize a generic framework for describing the dynamics of geophysical flows such as ocean currents or atmospheric wind fields in terms of Lagrangian flow networks. In this approach, information on the passive advection of particles is transformed into a Markov chain based on transition probabilities of particles between the volume elements of a given partition of space for a fixed time step. We employ perturbation-theoretic methods to investigate the effects of modifications of transport processes in the underlying flow for three different problem classes: efficient absorption (corresponding to particle trapping or leaking), constant input of particles (with additional source terms modeling, e.g., localized contamination), and shifts of the steady state under probability mass conservation (as arising if the background flow is perturbed itself). Our results demonstrate that in all three cases, changes to the steady state solution can be analytically expressed in terms of the eigensystem of the unperturbed flow and the perturbation itself. These results are potentially relevant for developing more efficient strategies for coping with contaminations of fluid or gaseous media such as ocean and atmosphere by oil spills, radioactive substances, non-reactive chemicals, or volcanic aerosols.

Percolation framework to describe El Niño conditions

Complex networks have been used intensively to investigate the flow and dynamics of many natural systems including the climate system. Here, we develop a percolation based measure, the order parameter, to study and quantify climate networks. We find that abrupt transitions of the order parameter usually occur ∼1 year before El Niño events, suggesting that they can be used as early warning precursors of El Niño. Using this method, we analyze several reanalysis datasets and show the potential for good forecasting of El Niño. The percolation based order parameter exhibits discontinuous features, indicating a possible relation to the first order phase transition mechanism.

Edge anisotropy and the geometric perspective on flow networks

Spatial networks have recently attracted great interest in various fields of research. While the traditional network-theoretic viewpoint is commonly restricted to their topological characteristics (often disregarding the existing spatial constraints), this work takes a geometric perspective, which considers vertices and edges as objects in a metric space and quantifies the corresponding spatial distribution and alignment. For this purpose, we introduce the concept of edge anisotropy and define a class of measures characterizing the spatial directedness of connections. Specifically, we demonstrate that the local anisotropy of edges incident to a given vertex provides useful information about the local geometry of geophysical flows based on networks constructed from spatio-temporal data, which is complementary to topological characteristics of the same flow networks. Taken both structural and geometric viewpoints together can thus assist the identification of underlying flow structures from observations of scalar variables.

Introduction to Focus Issue: Complex network perspectives on flow systems

During the last few years, complex network approaches have demonstrated their great potentials as versatile tools for exploring the structural as well as dynamical properties of dynamical systems from a variety of different fields. Among others, recent successful examples include (i) functional (correlation) network approaches to infer hidden statistical interrelationships between macroscopic regions of the human brain or the Earth's climate system, (ii) Lagrangian flow networks allowing to trace dynamically relevant fluid-flow structures in atmosphere, ocean or, more general, the phase space of complex systems, and (iii) time series networks unveiling fundamental organization principles of dynamical systems. In this spirit, complex network approaches have proven useful for data-driven learning of dynamical processes (like those acting within and between sub-components of the Earth's climate system) that are hidden to other analysis techniques. This Focus Issue presents a collection of contributions addressing the description of flows and associated transport processes from the network point of view and its relationship to other approaches which deal with fluid transport and mixing and/or use complex network techniques.

New Insights into Signed Path Coefficient Granger Causality Analysis

Granger causality analysis, as a time series analysis technique derived from econometrics, has been applied in an ever-increasing number of publications in the field of neuroscience, including fMRI, EEG/MEG, and fNIRS. The present study mainly focuses on the validity of “signed path coefficient Granger causality,” a Granger-causality-derived analysis method that has been adopted by many fMRI researches in the last few years. This method generally estimates the causality effect among the time series by an order-1 autoregression, and defines a positive or negative coefficient as an “excitatory” or “inhibitory” influence. In the current work we conducted a series of computations from resting-state fMRI data and simulation experiments to illustrate the signed path coefficient method was flawed and untenable, due to the fact that the autoregressive coefficients were not always consistent with the real causal relationships and this would inevitablely lead to erroneous conclusions. Overall our findings suggested that the applicability of this kind of causality analysis was rather limited, hence researchers should be more cautious in applying the signed path coefficient Granger causality to fMRI data to avoid misinterpretation.

Percolation Phase Transition of Surface Air Temperature Networks under Attacks of El Niño/La Niña

In this study, sea surface air temperature over the Pacific is constructed as a network, and the influences of sea surface temperature anomaly in the tropical central eastern Pacific (El Niño/La Niña) are regarded as a kind of natural attack on the network. The results show that El Niño/La Niña leads an abrupt percolation phase transition on the climate networks from stable to unstable or metastable phase state, corresponding to the fact that the climate condition changes from normal to abnormal significantly during El Niño/La Niña. By simulating three different forms of attacks on an idealized network, including Most connected Attack (MA), Localized Attack (LA) and Random Attack (RA), we found that both MA and LA lead to stepwise phase transitions, while RA leads to a second-order phase transition. It is found that most attacks due to El Niño/La Niña are close to the combination of MA and LA, and a percolation critical threshold Pc can be estimated to determine whether the percolation phase transition happens. Therefore, the findings in this study may renew our understandings of the influence of El Niño/La Niña on climate, and further help us in better predicting the subsequent events triggered by El Niño/La Niña.

Spatial network surrogates for disentangling complex system structure from spatial embedding of nodes

Networks with nodes embedded in a metric space have gained increasing interest in recent years. The effects of spatial embedding on the networks' structural characteristics, however, are rarely taken into account when studying their macroscopic properties. Here, we propose a hierarchy of null models to generate random surrogates from a given spatially embedded network that can preserve certain global and local statistics associated with the nodes' embedding in a metric space. Comparing the original network's and the resulting surrogates' global characteristics allows one to quantify to what extent these characteristics are already predetermined by the spatial embedding of the nodes and links. We apply our framework to various real-world spatial networks and show that the proposed models capture macroscopic properties of the networks under study much better than standard random network models that do not account for the nodes' spatial embedding. Depending on the actual performance of the proposed null models, the networks are categorized into different classes. Since many real-world complex networks are in fact spatial networks, the proposed approach is relevant for disentangling the underlying complex system structure from spatial embedding of nodes in many fields, ranging from social systems over infrastructure and neurophysiology to climatology.

Unified functional network and nonlinear time series analysis for complex systems science: The pyunicorn package

We introduce the pyunicorn (Pythonic unified complex network and recurrence analysis toolbox) open source software package for applying and combining modern methods of data analysis and modeling from complex network theory and nonlinear time series analysis. pyunicorn is a fully object-oriented and easily parallelizable package written in the language Python. It allows for the construction of functional networks such as climate networks in climatology or functional brain networks in neuroscience representing the structure of statistical interrelationships in large data sets of time series and, subsequently, investigating this structure using advanced methods of complex network theory such as measures and models for spatial networks, networks of interacting networks, node-weighted statistics, or network surrogates. Additionally, pyunicorn provides insights into the nonlinear dynamics of complex systems as recorded in uni- and multivariate time series from a non-traditional perspective by means of recurrence quantification analysis, recurrence networks, visibility graphs, and construction of surrogate time series. The range of possible applications of the library is outlined, drawing on several examples mainly from the field of climatology.

A Network of Networks Perspective on Global Trade

Mutually intertwined supply chains in contemporary economy result in a complex network of trade relationships with a highly non-trivial topology that varies with time. In order to understand the complex interrelationships among different countries and economic sectors, as well as their dynamics, a holistic view on the underlying structural properties of this network is necessary. This study employs multi-regional input-output data to decompose 186 national economies into 26 industry sectors and utilizes the approach of interdependent networks to analyze the substructure of the resulting international trade network for the years 1990-2011. The partition of the network into national economies is observed to be compatible with the notion of communities in the sense of complex network theory. By studying internal versus cross-subgraph contributions to established complex network metrics, new insights into the architecture of global trade are obtained, which allow to identify key elements of global economy. Specifically, financial services and business activities dominate domestic trade whereas electrical and machinery industries dominate foreign trade. In order to further specify each national sector's role individually, (cross-)clustering coefficients and cross-betweenness are obtained for different pairs of subgraphs. The corresponding analysis reveals that specific industrial sectors tend to favor distinct directionality patterns and that the cross-clustering coefficient for geographically close country pairs is remarkably high, indicating that spatial factors are still of paramount importance for the organization of trade patterns in modern economy. Regarding the evolution of the trade network's substructure, globalization is well-expressed by trends of several structural characteristics (e.g., link density and node strength) in the interacting network framework. Extreme events, such as the financial crisis 2008/2009, are manifested as anomalies superimposed to these trends. The marked reorganization of trade patterns, associated with this economic crisis in comparison to "normal" annual fluctuations in the network structure is traced and quantified by a new widely applicable generalization of the Hamming distance to weighted networks.

Local difference measures between complex networks for dynamical system model evaluation

A faithful modeling of real-world dynamical systems necessitates model evaluation. A recent promising methodological approach to this problem has been based on complex networks, which in turn have proven useful for the characterization of dynamical systems. In this context, we introduce three local network difference measures and demonstrate their capabilities in the field of climate modeling, where these measures facilitate a spatially explicit model evaluation.Building on a recent study by Feldhoff et al. [8] we comparatively analyze statistical and dynamical regional climate simulations of the South American monsoon system [corrected]. types of climate networks representing different aspects of rainfall dynamics are constructed from the modeled precipitation space-time series. Specifically, we define simple graphs based on positive as well as negative rank correlations between rainfall anomaly time series at different locations, and such based on spatial synchronizations of extreme rain events. An evaluation against respective networks built from daily satellite data provided by the Tropical Rainfall Measuring Mission 3B42 V7 reveals far greater differences in model performance between network types for a fixed but arbitrary climate model than between climate models for a fixed but arbitrary network type. We identify two sources of uncertainty in this respect. Firstly, climate variability limits fidelity, particularly in the case of the extreme event network; and secondly, larger geographical link lengths render link misplacements more likely, most notably in the case of the anticorrelation network; both contributions are quantified using suitable ensembles of surrogate networks. Our model evaluation approach is applicable to any multidimensional dynamical system and especially our simple graph difference measures are highly versatile as the graphs to be compared may be constructed in whatever way required. Generalizations to directed as well as edge- and node-weighted graphs are discussed.

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 autocorrelation in the records, which may introduce spurious links. This may explain why different methods could lead to different climate network topologies.