Search

The Amazon rainforest is considered one of the Earth's tipping elements and may lose stability under ongoing climate change. Recently a decrease in tropical rainforest resilience has been identified globally from remotely sensed vegetation data. However, the underlying theory assumes a Gaussian distribution of forest disturbances, which is different from most observed forest stressors such as fires, deforestation, or windthrow. Those stressors often occur in power-law-like distributions and can be approximated by $\alpha$-stable L\'evy noise. Here, we show that classical critical slowing down indicators to measure changes in forest resilience are robust under such power-law disturbances. To assess the robustness of critical slowing down indicators, we simulate pulse-like perturbations in an adapted and conceptual model of a tropical rainforest. We find few missed early warnings and few false alarms are achievable simultaneously if the following steps are carried out carefully: First, the model must be known to resolve the timescales of the perturbation. Second, perturbations need to be filtered according to their absolute temporal autocorrelation. Third, critical slowing down has to be assessed using the non-parametric Kendall-$\tau$ slope. These prerequisites allow for an increase in the sensitivity of early warning signals. Hence, our findings imply improved reliability of the interpretation of empirically estimated rainforest resilience through critical slowing down indicators.

Based on suggested interactions of potential tipping elements in the Earth's climate and in ecological systems, tipping cascades as possible dynamics are increasingly discussed and studied as their activation would impose a considerable risk for human societies and biosphere integrity. However, there are ambiguities in the description of tipping cascades within the literature so far. Here we illustrate how different patterns of multiple tipping dynamics emerge from a very simple coupling of two previously studied idealized tipping elements. In particular, we distinguish between a two phase cascade, a domino cascade and a joint cascade. While a mitigation of an unfolding two phase cascade may be possible and common early warning indicators are sensitive to upcoming critical transitions to a certain degree, the domino cascade may hardly be stopped once initiated and critical slowing down--based indicators fail to indicate tipping of the following element. These different potentials for intervention and anticipation across the distinct patterns of multiple tipping dynamics should be seen as a call to be more precise in future analyses on cascading dynamics arising from tipping element interactions in the Earth system.

Tipping elements in the Earth System receive increased scientific attention over the recent years due to their nonlinear behavior and the risks of abrupt state changes. While being stable over a large range of parameters, a tipping element undergoes a drastic shift in its state upon an additional small parameter change when close to its tipping point. Recently, the focus of research broadened towards emergent behavior in networks of tipping elements, like global tipping cascades triggered by local perturbations. Here, we analyze the response to the perturbation of a single node in a system that initially resides in an unstable equilibrium. The evolution is described in terms of coupled nonlinear equations for the cumulants of the distribution of the elements. We show that drift terms acting on individual elements and offsets in the coupling strength are sub-dominant in the limit of large networks, and we derive an analytical prediction for the evolution of the expectation (i.e., the first cumulant). It behaves like a single aggregated tipping element characterized by a dimensionless parameter that accounts for the network size, its overall connectivity, and the average coupling strength. The resulting predictions are in excellent agreement with numerical data for Erd\"os-R\'enyi, Barab\'asi-Albert and Watts-Strogatz networks of different size and with different coupling parameters.

Tipping elements occur in various systems such as in socio-economics, ecology and the climate system. In many cases, the individual tipping elements are not independent from each other, but they interact across scales in time and space. To model systems of interacting tipping elements, we here introduce the PyCascades open source software package for studying interacting tipping elements (doi: 10.5281/zenodo.4153102). PyCascades is an object-oriented and easily extendable package written in the programming language Python. It allows for investigating under which conditions potentially dangerous cascades can emerge between interacting dynamical systems, with a focus on tipping elements. With PyCascades it is possible to use different types of tipping elements such as double-fold and Hopf types and interactions between them. PyCascades can be applied to arbitrary complex network structures and has recently been extended to stochastic dynamical systems. This paper provides an overview of the functionality of PyCascades by introducing the basic concepts and the methodology behind it. In the end, three examples are discussed, showing three different applications of the software package. First, the moisture recycling network of the Amazon rainforest is investigated. Second, a model of interacting Earth system tipping elements is discussed. And third, the PyCascades modelling framework is applied to a global trade network.

In ecology, climate and other fields, (sub)systems have been identified that can transition into a qualitatively different state when a critical threshold or tipping point in a driving process is crossed. An understanding of those tipping elements is of great interest given the increasing influence of humans on the biophysical Earth system. Complex interactions exist between tipping elements, e.g. physical mechanisms connect subsystems of the climate system. Based on earlier work on such coupled nonlinear systems, we systematically assessed the qualitative long-term behaviour of interacting tipping elements. We developed an understanding of the consequences of interactions on the tipping behaviour allowing for tipping cascades to emerge under certain conditions. The (narrative) application of these qualitative results to real-world examples of interacting tipping elements indicates that tipping cascades with profound consequences may occur: the interacting Greenland ice sheet and thermohaline ocean circulation might tip before the tipping points of the isolated subsystems are crossed. The eutrophication of the first lake in a lake chain might propagate through the following lakes without a crossing of their individual critical nutrient input levels. The possibility of emerging cascading tipping dynamics calls for the development of a unified theory of interacting tipping elements and the quantitative analysis of interacting real-world tipping elements., Comment: 22 pages, 8 figures

Magnetic storms constitute the most remarkable large-scale phenomena of nonlinear magnetospheric dynamics. Studying the dynamical organization of macroscopic variability in terms of geomagnetic activity index data by means of complexity measures provides a promising approach for identifying the underlying processes and associated time-scales. Here, we apply a suite of characteristics from recurrence quantification analysis (RQA) and recurrence network analysis (RNA) in order to unveil some key nonlinear features of the hourly Disturbance storm-time (Dst) index during periods with magnetic storms and such of normal variability. Our results demonstrate that recurrence-based measures can serve as excellent tracers for changes in the dynamical complexity along non-stationary records of geomagnetic activity. In particular, trapping time (characterizing the typical length of "laminar phases" in the observed dynamics) and recurrence network transitivity (associated with the number of the system's effective dynamical degrees of freedom) allow for a very good discrimination between magnetic storm and quiescence phases. In general, some RQA and RNA characteristics distinguish between storm and non-storm times equally well or even better than other previously considered nonlinear characteristics like Hurst exponent or symbolic dynamics based entropy concepts. Our results point to future potentials of recurrence characteristics for unveiling temporal changes in the dynamical complexity of the magnetosphere.

The effects of El Ni\~no's two distinct flavors, East Pacific (EP) and Central Pacific (CP)/Modoki El Ni\~no, on global climate variability have been studied intensively in recent years. Most of these studies have made use of linear multivariate statistics or composite analysis. Especially the former assumes the same type of linear statistical dependency to apply across different phases of the El Ni\~no--Southern Oscillation (ENSO), which appears not necessarily a justified assumption. Here, we statistically evaluate the likelihood of co-occurrences between very high or very low seasonal precipitation sums over vast parts of the global land surface and the presence of the respective EP and CP types of both, El Ni\~no and La Ni\~na. By employing event coincidence analysis, we uncover differential imprints of both flavors on very low and very high seasonal precipitation patterns over distinct regions across the globe, which may severely affect, among others, agricultural and biomass production or public health. We particularly find that EP periods exhibit statistically significant event coincidence rates with hydrometeorological anomalies at larger spatial scales, whereas sparser patterns emerge along with CP periods. Our statistical analysis confirms previously reported interrelations for EP periods and uncovers additional distinct regional patterns of very high/low seasonal precipitation, such as increased rainfall over Central Asia alongside CP periods that have to our knowledge not been reported so far. Our results demonstrate that a thorough distinction of El Ni\~no and La Ni\~na into their two respective flavors could be crucial for understanding the emergence of strong regional hydrometeorological anomalies and anticipating their associated ecological and socioeconomic impacts.

El Ni\~no exhibits distinct Eastern Pacific (EP) and Central Pacific (CP) types which are commonly, but not always consistently, distinguished from each other by different signatures in equatorial climate variability. Here, we propose an index based on evolving climate networks to objectively discriminate between both flavors by utilizing a scalar-valued evolving climate network measure that quantifies spatial localization and dispersion in El Ni\~no's associated teleconnections. Our index displays a sharp peak (high localization) during EP events, whereas during CP events (larger dispersion) it remains close to the baseline observed during normal periods. In contrast to previous classification schemes, our approach specifically account for El Ni\~no's global impacts. We confirm recent El Ni\~no classifications for the years 1951 to 2014 and assign types to those cases were former works yielded ambiguous results. Ultimately, we study La Ni\~na episodes and demonstrate that our index provides a similar discrimination into two types.

We introduce the \texttt{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. \texttt{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, \texttt{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 (RQA), 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., Comment: 28 pages, 17 figures

In recent years extensive studies on the Earth's climate system have been carried out by means of advanced complex network statistics. The great majority of these studies, however, have been focusing on investigating correlation structures within single climatic fields directly on or parallel to the Earth's surface. Here, we develop a novel approach of node weighted coupled network measures to study correlations between ocean and atmosphere in the Northern Hemisphere extratropics and construct 18 coupled climate networks, each consisting of two subnetworks. In all cases, one subnetwork represents monthly sea-surface temperature (SST) anomalies, while the other is based on the monthly geopotential height (HGT) of isobaric surfaces at different pressure levels covering the troposphere as well as the lower stratosphere. The weighted cross-degree density proves to be consistent with the leading coupled pattern obtained from maximum covariance analysis. Network measures of higher order allow for a further analysis of the correlation structure between the two fields and consistently indicate that in the Northern Hemisphere extratropics the ocean is correlated with the atmosphere in a hierarchical fashion such that large areas of the ocean surface correlate with multiple statistically dissimilar regions in the atmosphere. Ultimately we show that, this observed hierarchy is linked to large-scale atmospheric variability patterns, such as the Pacific North American pattern, forcing the ocean on monthly time scales.

Complex network theory provides a powerful toolbox for studying the structure of statistical interrelationships between multiple time series in various scientific disciplines. In this work, we apply the recently proposed climate network approach for characterizing the evolving correlation structure of the Earth's climate system based on reanalysis data of surface air temperatures. We provide a detailed study on the temporal variability of several global climate network characteristics. Based on a simple conceptual view on red climate networks (i.e., networks with a comparably low number of edges), we give a thorough interpretation of our evolving climate network characteristics, which allows a functional discrimination between recently recognized different types of El Ni\~no episodes. Our analysis provides deep insights into the Earth's climate system, particularly its global response to strong volcanic eruptions and large-scale impacts of different phases of the El Ni\~no Southern Oscillation (ENSO)., Comment: 20 pages, 12 figures

Eigen techniques such as empirical orthogonal function (EOF) or coupled pattern (CP) / maximum covariance analysis have been frequently used for detecting patterns in multivariate climatological data sets. Recently, statistical methods originating from the theory of complex networks have been employed for the very same purpose of spatio-temporal analysis. This climate network (CN) analysis is usually based on the same set of similarity matrices as is used in classical EOF or CP analysis, e.g., the correlation matrix of a single climatological field or the cross-correlation matrix between two distinct climatological fields. In this study, formal relationships as well as conceptual differences between both eigen and network approaches are derived and illustrated using exemplary global precipitation, evaporation and surface air temperature data sets. These results allow to pinpoint that CN analysis can complement classical eigen techniques and provides additional information on the higher-order structure of statistical interrelationships in climatological data. Hence, CNs are a valuable supplement to the statistical toolbox of the climatologist, particularly for making sense out of very large data sets such as those generated by satellite observations and climate model intercomparison exercises., Comment: 18 pages, 11 figures

Network theory provides various tools for investigating the structural or functional topology of many complex systems found in nature, technology and society. Nevertheless, it has recently been realised that a considerable number of systems of interest should be treated, more appropriately, as interacting networks or networks of networks. Here we introduce a novel graph-theoretical framework for studying the interaction structure between subnetworks embedded within a complex network of networks. This framework allows us to quantify the structural role of single vertices or whole subnetworks with respect to the interaction of a pair of subnetworks on local, mesoscopic and global topological scales. Climate networks have recently been shown to be a powerful tool for the analysis of climatological data. Applying the general framework for studying interacting networks, we introduce coupled climate subnetworks to represent and investigate the topology of statistical relationships between the fields of distinct climatological variables. Using coupled climate subnetworks to investigate the terrestrial atmosphere's three-dimensional geopotential height field uncovers known as well as interesting novel features of the atmosphere's vertical stratification and general circulation. Specifically, the new measure "cross-betweenness" identifies regions which are particularly important for mediating vertical wind field interactions. The promising results obtained by following the coupled climate subnetwork approach present a first step towards an improved understanding of the Earth system and its complex interacting components from a network perspective.

We propose a method to reconstruct and analyze a complex network from data generated by a spatio-temporal dynamical system, relying on the nonlinear mutual information of time series analysis and betweenness centrality of complex network theory. We show, that this approach reveals a rich internal structure in complex climate networks constructed from reanalysis and model surface air temperature data. Our novel method uncovers peculiar wave-like structures of high energy flow, that we relate to global surface ocean currents. This points to a major role of the oceanic surface circulation in coupling and stabilizing the global temperature field in the long term mean (140 years for the model run and 60 years for reanalysis data). We find that these results cannot be obtained using classical linear methods of multivariate data analysis, and have ensured their robustness by intensive significance testing., Comment: 6 pages, 5 figures

Complex network theory provides a powerful framework to statistically investigate the topology of local and non-local statistical interrelationships, i.e. teleconnections, in the climate system. Climate networks constructed from the same global climatological data set using the linear Pearson correlation coefficient or the nonlinear mutual information as a measure of dynamical similarity between regions, are compared systematically on local, mesoscopic and global topological scales. A high degree of similarity is observed on the local and mesoscopic topological scales for surface air temperature fields taken from AOGCM and reanalysis data sets. We find larger differences on the global scale, particularly in the betweenness centrality field. The global scale view on climate networks obtained using mutual information offers promising new perspectives for detecting network structures based on nonlinear physical processes in the climate system., Comment: 24 pages, 10 figures