5ee57354c91b85a7b0acb3bdab3ad7c2.ppt

- Количество слайдов: 74

Climate Networks A short journey of complex systems science in climate studies Jürgen Kurths Potsdam-Institut für Klimafolgenforschung Institut für Physik, Humboldt-Universität zu Berlin RAS Inst. Applied Physics, Nizhny Novgorod juergen. [email protected] de

Contents • Introduction • Fractals & Deterministic chaos – utopistic expectations and applications to climate (70 ies, 80 ies) • Structure formation in coupled climate subsystems (90 ies) • Complex networks and climate (1998 - ? ) • Conclusions & Open Problems (2016 - ? )

Nonlinear Sciences Start in 1665 by Christiaan Huygens: Discovery of phase synchronization, called sympathy

Pendulum Clocks • Christiaan Huygens: Pendelum clocks hanging at the same wooden beam (half-timber house) It is quite worth noting that when we suspended two clocks so constructed from two hooks imbedded in the same wooden beam, the motions of each pendulum in opposite swings were so much in agreement that they never receded the last bit from each other. . . Further, if this agreement was disturbed by some interference, it reastablished itself in a short time. . . after a careful examination I finally found that the cause of this is due to the motion of the beam, even though this is hardly perceptible (Huygens, 1673)

Basic feature of a complex system ability for self-organization non-trivial (some times unexpected) structure formation (here synchronization due to almost hidden coupling)

If we knew exactly the laws of nature and the situation of the universe at the initial moment, we could predict exactly the situation of that same universe at the succeeding moment. but even if it were the case that the natural laws had no longer any secret for us, we could still only know the initial situation approximately. If that enabled us to predict the succeeding situation with the same approximation, that is all we require, and we should say that the phenomenon had been predicted, that it is governed by laws. But it is not always so; it may happen that small differences in the initial conditions produce very great ones in the final phenomena. A small error in the former will produce an enormous error in the latter. Prediction becomes impossible, and we have the fortuitous phenomenon. Concept of Weak Causality (Poincare`, 1903)

Achievements of Nonlinear Science/ Complex Systems Science Potentials & Limits

Basic Approaches 1) Fractals

Mathematical Monsters: Koch Curve (Snowflake) - 1904 5. Iteration

Problem in traditional Understanding • Length n-th Iteration: U(n) = (4 / 3)ⁿ Limit: U(n) → ∞ für n → ∞ infinitely long • But: Area finite (bounded) • Fractal – Fractal Dimension (Felix Hausdorff, 1919)

Fractals (B. B. Mandelbrot) (70 ies, 80 ies) New kind of connection between Complexity and Infinity New Paradigm: Self-Similarity

Coastline South Norway D = 1. x

2) Deterministic Chaos (Yorke, 1975) Deterministic Chaos (Schuster, 1984) Originally: Henri Poincare´ (1854 – 1912)

Simple map generates huge variety of dynamics (in dependence on parameter r) xn+1 = r xn ( 1 – xn )

CHAOS: Sensitivity to initial conditions r = 3. 8, x(0) = 0. 4, x(0) + Δ = 0. 41 strong growth of initial error Exponential growth – Lyapunov Exponent λ > 0 Deterministic Chaos r = 4 λ = 1 (fully developed chaos)

Periodic windows, inner crisis, intermittency

Intermittency alternation of laminar and bursting regimes (chaos „partly visible“ in periodic phase)

Techniques of nonlinear time series analysis developed • Estimates of fractal dimensions (Grassberger/Procaccia, 1983) and • max Lyapunov exponents (Wolf, 1984) • Successful applications to lab experiments (laser, fluid dynamics) (utopistic) expectation (hype): dynamics of many real-world systems is low-dimensional chaos (2 numbers!!!)

„Hunting“ for a low-dimensional climatic attractor Climate • D = 3. 1 Nicolis… Nature 1984 • D= 4. 4 – 4. 8 Fraedrich 1986 Weather • D = 5 - 6 • 7. 3 • 7 – 8 Essex, Nature 1987 Tsonis/Elsner, Nature 1988 Keppene 1989

All these „results“ are artifacts of the method (intrinsic correlations) No statistical evidence for a lowdimensional climatic attractor Climate cannot be described by a few ordinary differential equations Grassberger, … late 80 ies, early 90 ies

Depression time (short) after too strong promises

New Phenomena and Techniques (golden 90 ies and later) • Synchronization of chaotic systems (complete, generalized, phase; cluster) • Amplitude/ Oscillation Death, Chimera • Control of Chaos • Communication with chaos (carrier; secure, efficient) • Noise-induced order (stochastic resonance, coherence resonance, synchronization…)

Applications of complex synchronization • lasers, circuits, engineering, chemical oscillators • Physiology: mostly two coupled: heart and breathing rhythms • Neuroscience: pairwise brain regions… • Climate: ENSO vs. Indian Monsoon

El Niño vs. Indian Monsoon

Wavelet Analysis Thick black curve – statistically significant Thin black curve – not-correct statistical significance region GRL 2005, 2011 (with D. Smirnov and I. Mokhov)

Next big step (1998 - ) Networks with complex topology approaching real-world much better

Erdös-Renyi Graphs – Random Networks Distribution of degrees - Poisson

Small-world Networks k nearest neighbour connections Nearest neighbour and a few long-range connections Regular Complex Topology

Evolving Networks

Network of Networks Interconnected Networks Interdependent Networks Multiplex Networks Multilayer Networks…

Basic Studies

Fields of applications • Spreading of diseases (humans, computers…) • Spreading of information - Opinion formation (social networks) • Structure formation in socio-economic networks • Cascading failure (transport systems) • Functioning of the brain • Stability of power grids

System Earth

Increasing CO 2 -Concentration 2014: 400 ppm IPCC 2007

Is there a break of global warming? Source NASA GISS

No weakening of global warming

Complex Network Approach to Climate

Network Reconstruction from a continuous dynamic system (structure vs. functionality) New (inverse) problems arise! Is there a backbone underlying the climate system?

Basic Idea: Use of rich instrumentarium of complex network (graph) theory for system Earth and sustainability Hope: Deepened understanding of system Earth (with other techniques NOT possible)

Climate Networks Observation sites Earth system Climate network Network analysis Time series

Infer long-range connections – Teleconnections

Artifacts and Interpretation of (Climate) Network Approach

Reconstructing causality from data ? Artefacts due to - Common drivers - Indirect links Achievements 1. Causal algorithm to efficiently detect linear and nonlinear links (Phys. Rev. Lett. 2012) 2. Quantifying causal strength with Momentary Information Transfer (Phys. Rev. E 2012) 3. Reconstructing Walker Circulation from data (J. Climate 2014, )

Reconstructing causality from data Classic techniques Correlation/regression Advanced method conditional independencies

Identifying causal gateways and mediators in complex spatio-temporal systems • Step 1: Dimension reduction via VARIMAX (principal components, rotation, significance) • Step 2: Causal reconstruction: identify causalities based on conditional dependencies (different time lags) • Step 3: Causal interaction quantification: identify strongest paths • Step 4: Hypothesis testing of causal mechanisms Nature Commun, 6, 8502 (2015)

Each node corresponds to a lagged subprocess and owing to stationarity links for t imply links for all t-1, t-2, . . Process (black node) has two parents (gray nodes, connected via incoming links from the past) and one neighbor (hatched node, connected with an undirected contemporaneous link). In the case of a linear model, each multivariate regression coefficient a, b, and c can be attributed to a link. Parents

Atmospheric data • • • Reanalysis data – NCEP/NCAR (Boulder) surface pressure 1948 – 2012 Spatial resolution: 2. 5º → 10, 512 grid points Weekly data: each node time series of 3, 339 points • Big data

60 strongest VARIMAX components refer to main climatic patterns • ENSO: “ 0” – western uplift, “ 1” – eastern downdraft limbs • Monsoon: “ 33” Arabian Sea high-surface-pressure sector, “ 26” tropical Atlantic West African Monsoon system

Identification of causal pathways Effects of sea level pressure anomalies in ENSO region to pressure variability in the Arabic Sea via the Indonesian Archipelago

Extreme Events Strong Rainfall during Monsoon Challenge: Predictability

Motivation: the predictability of the Indian monsoon remains a problem of vital importance Objectives: to reveal spatial structures in network of extreme events over the Indian subcontinent and their seasonal evolution during the year.

New Technique: Event Synchronization

2. Event synchronization 1. Network approach Nodes: geographical locations Links: synchronization of extreme rainfall events between nodes Step 1. Apply a threashold to time series of each grid point to obtain event series Rainfall amount (mm/day) METHOD Time (days) Step 2. Event synchronization – use time lags to compare individual events between two grid points Step 3. Construct the network by creating links between points with the highest synchronization values Network measures degree Average link lengths betweenness Network Quiroga et. al. 2002 59 Malik et. al. 2011 Boers et. al. 2013

Extreme Rainfall Events of the South American Monsoon System • • TRMM 3 B 42 V 7 daily satellite data Measured: Jan 1, 1998 – March 31, 2012 Spatial resolution: 0. 25 x 0. 25 Spatial coverage: Method: event synchronization Extreme event: > 99 % percentile Dec-Feb (DJF) – summer monsoon months

Non-symmetric Adjacency Matrix (in – out) > 0 – sink: extreme events here preceded by those at another location < 0 – source: extreme events follow at another location

• > 60 % (90 % during El Nino conditions) of extreme rainfall events in Eastern Central Andes (ECA) are preceded by those in Southeastern South America (SESA) • Low pressure anomaly from Rossbywave activity propagates northwards (cold front) and low-level wind channel from Amazon Nature Communication (2014), GRL (2014), J. Clim. (2014), Clim. Dyn. (2015)

Analysis of Simulation (Model) Data Case Study: Atlantic Meridional Overturning Circulation (MOC) With Henk Dijkstra & group GRL 40, 4386 (2013)

How stable is the conveyor belt? Changing parameter: anomalous fresh water flux (from North Atlantic)

2 D meridional depth model of the Atlantic MOC (Den Toom, Dijkstra, Wubs, 2011) 32 x 16 (x, z) spatial grid Network reconstruction

Our papers on climate networks • Europhys. Lett. 87, 48007 (2009) • Phys. Rev. E 81, 015101 R (2010) • Climate Dynamics 39, 971 (2012) • PNAS 108, 20422 (2011) • Phys. Rev. Lett. 106, 258701 (2012) • Europhys. Lett. 97, 40009 (2012) • Climate Past 8, 1765 (2012) • Geophys. Res. Lett. 40, 2714 (2013) • Climate Dynamics 41, 3 (2013) • J. Climate 27, 720 (2014) • Nature Scientific Reports 4, 4119 (2014) • Climate Dynamics (2014) • Geophys. Res. Lett. 41, 7397 (2014) • Nature Commun. 5, 5199 (2014) • Climate Dynamics 44, 1567 (2015) • J. Climate 28, 1031 (2015) • Climate Past 11, 709 (2015) • Climate Dynamics (online 2015) • Nature Commun. 6, 8502 (2015)

Codes available • Unified functional network and nonlinear time series analysis for complex systems science: The pyunicorn software package, CHAOS 25, 113101 (2015) https: //github. com/pik-copan/pyunicorn • Causal network identification: Python software script by J. Runge http: //tocsy. pik-potsdam. de/tigramite. php

Summary • Complex systems science promising approach in climate science, mitigation and adaptation • But in its infancy – many open problems large room for new research - new basic phenomena - new techniques - challenging application