Non-autonomous Dynamics in Complex Systems: Theory and Applications to Critical Transitions

These lectures are designed to provide a brief and limited introduction to the theoretical treatment of dynamical systems induced by nonautonomous differential equations and their critical transitions in the form of nonautonomous bifurcations. We will present the mathematical formalisms of processes and skew-product flows and introduce the notions of pullback attraction, exponential dichotomy and hyperbolic solutions and understand their meaning via several examples. Then we shall focus on the nonautonomous fold bifurcation as a mechanism to generate a critical transition in certain classes of scalar nonautonomous ODEs. During the "hands-on" session we shall review all the material by analysing the dynamical behaviour of a simple zero-dimensional energy balance model for a planet's climate.

These lectures are designed to provide a brief and limited introduction to the theoretical treatment of dynamical systems induced by nonautonomous differential equations and their critical transitions in the form of nonautonomous bifurcations. We will present the mathematical formalisms of processes and skew-product flows and introduce the notions of pullback attraction, exponential dichotomy and hyperbolic solutions and understand their meaning via several examples. Then we shall focus on the nonautonomous fold bifurcation as a mechanism to generate a critical transition in certain classes of scalar nonautonomous ODEs. During the "hands-on" session we shall review all the material by analysing the dynamical behaviour of a simple zero-dimensional energy balance model for a planet's climate.

In this lecture, we will discuss ways to detect critical transitions in complex systems based on analyses of observations of their non-autonomous dynamics. Topics include segmentation techniques for investigating non-stationary dynamics, univariate and multivariate time-series-analysis techniques as well as surrogate-based concepts to evaluate reliability of detections. Special attention is paid to sources of errors and potential remedies.

The study of dynamical systems concerns the long-term behaviour of a set of variables (the state of the system) based only on the knowledge of an initial condition and a set of rules (the evolution equations). Most studied examples are differential equations that model large classes of applications such as planetary motion, chemical reactions, population dynamics, consumer behaviour or stock markets. In contrast to an autonomous dynamical system, the time evolution of a nonautonomous dynamical system is influenced by an independent process (that reflects changes in the rules governing the system). The importance of nonautonomous dynamical systems is illustrated by the fact that a significant number of real-world applications are governed by time-dependent inputs, and the traditional mathematical theory fails to address this. There has been rapid progress concerning the development of a theory for nonautonomous dynamical systems in the last thirty years. In this talk, an introduction to some elementary aspects of the theory is given, and similarities and differences to the autonomous case will be highlighted. The theoretical results are illustrated by an application to the terrestrial carbon cycle, which is modeled by a nonautonomous compartmental system.

Heteroclinic dynamics is a suitable framework for describing transient and reproducible dynamics such as cognitive processes in the brain. We demonstrate how heteroclinic units can act as pacemakers to entrain larger sets of units from a resting state to hierarchical heteroclinic motion such as fast oscillations modulated by slow oscillations. The entrainment range depends on the type of coupling, the spatial location of the pacemaker and the individual bifurcation parameters of the pacemaker and the driven units. Noise as well as a small back-coupling to the pacemaker considerably facilitate synchronization. Units can be synchronously entrained to different temporal patterns, depending on the selected path in the hierarchical heteroclinic network. Such patterns are believed to code information in brain dynamics. Depending on the number and the location of pacemakers on a two-dimensional grid, synchronization can be maintained in the presence of a large number of resting state units and mediated via target waves when the pacemakers are concentrated to a small area of such a grid. In view of brain dynamics, our results indicate a possibly ample repertoire for coding information in temporal patterns, produced by sets of synchronized units entrained by pacemakers, without finetuning of the parameters, and distributed in space.

Most complex and networked systems are out of equilibrium and externally driven such that their collective dynamics may change drastically. They may exhibit tipping, here taken to be a qualitative change induced by sufficiently strong driving, that often disrupts the systems' intended or desired functionality. Yet, standard linear response theory generically characterizes systems responses to such fluctuations only for small driving amplitudes and cannot capture response properties that are either due to strong driving signals or intrinsically nonlinear. Here we present an alternative perturbative technique that exploits a self-consistency criterion, captures the full nonlinear response dynamics and enables the approximate prediction of tipping points and the response dynamics before tipping.

The analysis of the warming trend of surface temperature anomalies shows an abrupt acceleration of warming in the past 15 years in central Europe (ERA5 2mT). Employing change point analysis and quantifying the significance of trends in the presence of long range temporal correlations, we try to answer the question of whether this change of warming rate is related to overcoming tipping points in the release of greenhouses gases by correlating these results to greenhouse gas concentrations and emissions.

In this talk, we show that the celebrated Embedding Theorem of Takens is a particular case of a much more general statement according to which, randomly generated linear state-space representations of generic observations of an invertible dynamical system carry in their wake an embedding of the phase space dynamics into the chosen Euclidean state space. This embedding coincides with a natural generalized synchronization that arises in this setup and which yields a topological conjugacy between the state-space dynamics driven by the generic observations of the dynamical system and the dynamical system itself. This result provides additional tools for the representation, learning, and analysis of chaotic attractors and sheds additional light on the reservoir computing phenomenon that appears in the context of recurrent neural networks. This work is a collaboration with Lyudmila Grigoryeva (St. Gallen, Switzerland) and Allen Hart (Exeter, UK).

In the age of climate warming, comprehension of ecosystems’ future is one of the pressing challenges to humanity. Whilst most studies on climate warming focus on the ‘magnitude of change’ of the Earth’s temperature, the ‘rate’ at which it is increasing cannot be ruled out. Rapid warming has already caused sudden ecosystem transitions at numerous biodiversity hot spots; a mechanistic understanding of such transitions is crucial. Here, we study a slow–fast consumer–resource ecosystem interacting in rapid warming scenarios. We find that while a gradual change in mean temperature may accord population persistence, a critical warming rate can drive the resource’s sudden collapse, termed a warming-induced abrupt transition. This further triggers the bottom-up effect, resulting in the extinction of the consumer. The difference between the optimum temperature of the resource’s growth rate and the habitat temperature is crucial in deciding the critical rate of warming. Consequently, species inhabiting extreme temperature regions are more susceptible to warming-induced collapse than those within intermediate temperature ranges. We find that stochastic fluctuations in the system can advance warming-induced transitions, and the efficacy of generic early warning signals to anticipate sudden transitions is challenged.

Various physical phenomena can be modelled through the use of fast-slow systems in the form of PDEs perturbed by noise. A common technique for a proper study of their properties is the focus on the fast variable, hence fixing the slow one and considering it as a parameter whose value can affect the qualitative behaviour of the solution of the equations. Of particular importance are the bifurcation thresholds of the deterministic system and the early-warning signs that can be observed as the parameter approaches them. I will compare different types of early warning-signs for SPDEs (stochastic partial differential equations) and discuss how they are affected by the choice of the drift component. The nature of the noise is a generalized Q-Wiener process in order to replicate the heterogeneity in space often displayed in applications.

Kernel methods are a family of machine learning algorithms that have been applied successfully to a wide range of fields. Kernel methods have primarily been applied to Euclidean data. However in many applications - such as climate science, geometric mechanics and computer vision - the data naturally lies on a non-Euclidean manifold. We investigate the generalisation of kernel based algorithms to non-Euclidean spaces. Kernel methods often rely on embedding the data space into a reproducing kernel Hilbert space (RKHS), which in turn relies on a positive definite (PD) kernel. We find that the geodesic Gaussian kernel, that generalises the Euclidean Gaussian kernel to Riemannian manifolds, is usually not (PD), preventing its use for learning with RKHS. We provide alternative ways of finding PD kernels on Riemannian symmetric spaces. We also show that, on many Riemannian symmetric spaces, the geodesic Gaussian kernel satisfies the weaker condition of admitting a positive decomposition. This gives an embedding of the data space into a reproducing kernel Krein space (RKKS) instead, providing an alternative route for learning with this kernel.

Zombie fires in peatlands, that disappear from the surface, smoulder underground during the winter, and `come back to life' the following spring, release Gigatonnes of carbon into the atmosphere. These fires are generally believed to be caused by surface wildfires. However, we found that rate-induced tipping to a subsurface hot metastable state in bioactive peat soils, due to too-fast rates of atmospheric warming, can be a main cause of Zombie fires. Our hypothesis is based on a soil-carbon ODE model subjected to realistic global warming scenarios and summer heatwaves. Mathematically speaking, R-tipping to the hot metastable state is a nonautonomous instability, due to crossing an elusive quasithreshold, in a multiple timescale dynamical system. To explain this instability, we provide a mathematical framework that combines a special compactification technique with concepts from geometric singular perturbation theory. This framework allows us to reduce an R-tipping problem due to crossing a quasithreshold to an `exotic' heteroclinic orbit problem in a singular limit. We identify generic cases of such R-tipping via unfolding of a codimension-two heteroclinic folded saddle-node type-I singularity, which in turn reveals new types of excitability quasithresholds.

In systems of stochastic dynamical systems, noise-induced transitions occur when the presence of noise prompts the state of the system to transition from one (meta)stable state to another, potentially as a rare event. Understanding a system’s behavior through this transition is a significant area of current interest. In this talk, we will describe recent analytical and computational tools to identify Most Probable Escape Paths (MPEPs) in Stochastic Differential Equations (SDEs) over periodic boundaries. The study of noise-induced tipping has been dominated by Freidlin-Wentzell (FW) theory. One limitation of the FW theory is that it necessitates vanishingly small noise. However, for several applications of interest, particularly in environmental, biological, or social contexts, intermediate noise is of more relevance. In focusing on the intermediate noise regime, the transient behavior of the underlying deterministic system will play a key role. We will thus use a dynamical system approach to identify MPEPs in SDEs. The Maslov index will help us distinguish which critical points of the FW functional are minimizers and help explain the effect of the interaction of noise and transient dynamics. The Onsager-Machlup (OM) functional, which is treated as a perturbation of the FW functional, will provide a selection mechanism to pick out a specific MPEP. Subsequently, our computations will be compared with Monte Carlo simulations to verify theoretical predictions in the Inverted Van der Pol system.

The Atlantic Meridional Overturning Circulation is responsible for the comparatively temperate climate found in Western Europe, and its previous collapse thought to have triggered glacial periods seen in the paleo data. This is a system that has multiple stable states- referred to as ‘on’ when the circulation is strong as in the current climate, and ‘off’ when it is much weaker. The AMOC has tipping points between these states. Tipping points occur when a rapid shift in dynamics happens in response to a relatively small change in a parameter. Making future projections of AMOC response to the climate is essential for avoiding any anthropogenic caused tipping, but it is computationally expensive to calculate the full hysteresis for different scenarios. This work looks at a conceptual five box model of the AMOC [1] which is easy to understand and cheap to implement. Previous work has considered bifurcation and rate-dependent tipping [2]. This current work looks to estimate a realistic amount of noise from various GCM data sets and apply this to the model. We compare the covariance of the salinity data for a variety of CMIP6 models, and we compare the amount of noise covariance found in each data set, and how this can be input back into the box model. Ongoing work looks to find first passage times for the box models with the estimated noise and a small amount of freshwater forcing. (1) Wood, R et.al. (2019), Climate Dynamics, 53(11), 6815-6834 (2) Alkhayuon, H et.al. (2019), Proc. R. Soc. A, 475(2225)

We consider discrete-time random dynamical systems with uniform spherical noise and introduce a geometric approach to study the compound behaviour of such systems. The key insight of this novel approach is that the boundary of attractors can be represented as invariant sets of a deterministic finite-dimensional mapping, the so-called boundary mapping, that acts on the unit tangent bundle of the phase space. We address questions regarding the persistence of attractors and the nature of bifurcations in this context.

The Atlantic Meridional Overturning Circulation (AMOC) was identified as a potential tipping element in the earth system. In a hierarchy of models of increasing complexity, ranging from simple box models to intermediate-complexity models to state-of-the-art climate models, the AMOC exhibits bistability with an 'on-state' characterised by a functioning overturning circulation and an 'off-state' in which the overturning circulation shuts down. The proximity of the present AMOC to a critical transition threshold is under debate. The impact of an AMOC collapse on the climate in many parts of the world would be enormous. In this study, a data-driven methodology for identifying, anticipating and predicting critical transitions in high-dimensional model or observational data sets is introduced, based on explicit non-stationary low-order modelling of the tipping dynamics. Unlike the more traditional early-warning signs such as rising autocorrelation (critical slowing down) or rising variance, this allows for dynamical understanding of the underlying tipping mechanism and genuine prediction of the future system state by extrapolation. A set of spatial modes carrying the tipping dynamics are identified and a stochastic model of appropriate complexity is estimated in the subspace spanned by these modes. Analysis of the reconstructed dynamics allows to determine the proximity to a bifurcation point and the type of the impending bifurcation. Different competing tipping mechanisms can be compared and assessed using likelihood inference and information criteria. The technique is inherently probabilistic with the randomness stemming from two sources: the statistical uncertainty in the model estimation and the inevitably stochastic nature of the low-order tipping model. As such, the method allows to quantify the likelihood or risk of a critical transition at some point in the future having observed a certain amount of data up to present. The methodology is here applied to a data set from a simulation of AMOC collapse with a complex climate model, actually a freshwater hosing experiment with the FAMOUS GCM. The AMOC on-state is found to lose stability via a subcritical Hopf bifurcation; however, the transition to the off-state occurs far ahead of the bifurcation point. The early collapse can be explained by a combination of rate-induced and noise-induced tipping. The finding is relevant in terms of safe operating zones as the risk of AMOC collapse in such a scenario would be higher than anticipated with conventional early-warning signs.

In this study, we propose a simplified model of a socio-environmental system that accounts for population, resources, and wealth, with a quadratic population contribution in the resource extraction term. Given its structure, an analytical treatment of attractors and bifurcations is possible. In particular, a Hopf bifurcation from a stable fixed point to a limit cycle emerges above a critical value of the extraction rate parameter. The stable fixed-point attractor can be interpreted as a sustainable regime, and a large-amplitude limit cycle as an unsustainable regime. The model is generalized to multiple interacting systems, with chaotic dynamics emerging for small non-uniformities in the interaction matrix. In contrast to systems where a specific parameter choice or high dimensionality is necessary for chaos to emerge, chaotic dynamics here appears as a generic feature of the system. In addition, we show that diffusion can stabilize networks of sustainable and unsustainable societies, and thus, interconnection could be a way of increasing resilience in global networked systems. Overall, the multi-systems model provides a timescale of predictability (300-1000 years) for societal dynamics comparable to results from other studies, while indicating that the emergent dynamics of networks of interacting societies over longer time spans is likely chaotic and hence unpredictable. In addition, I will give an introduction to societal collapse and contextualize the above work in my broader research area.

Climate change and the development of drier climates threaten ecosystems’ health and the services they provide to humans. Understanding the response of ecosystems to drier climates may provide clues on how to improve their resilience. Two robust mechanisms that enhance ecosystem resilience are phenotypic changes in individual plants and partial mortality of plant populations to form vegetation patterns. In nature, these mechanisms are likely to act in concert, but their interplay has escaped consideration. We demonstrate the need for a theory that integrates these plant-level and population-level mechanisms by addressing the fascinating fairy-circle phenomenon in Namibia. We show that such an integration resolves two outstanding puzzles in the current theory: observations of multi-scale patterns and the absence of theoretically predicted large-scale stripe and spot patterns along the rainfall gradient. Importantly, we find that multi-level responses to stress unveil a wide variety of more effective stress relaxation pathways, compared to single-level responses, implying a previously underestimated resilience of dryland ecosystems.

We consider systems, which are governed by typical scalar differential equations with a stabilizing linear and a delayed nonlinear term. We investigate the influence of periodically modulating the delay time. We show that such an explicit time-dependence can fundamentally change the system properties. This is demonstrated for nonlinearities which allow for a chaotic, diffusive motion of the state variable. In this case a change of the mean delay or the modulation amplitude can lead to changes of the diffusion constant by orders of magnitude [1,2], or, for certain nonlinearities, can even turn the system from being normally diffusive to one showing anomalous diffusion and weak ergodicity breaking. [1] T. Albers, D. Müller-Bender, L. Hille, and G. Radons, Chaotic diffusion in delay systems: Giant enhancement by time lag modulation, Phys. Rev. Lett. 128, 074101, Editors' Suggestion (2022) [2] T. Albers, D. Müller-Bender, and G. Radons, Anti-persistent random walks in time-delayed systems, Phys. Rev. E 105, 064212 (2022)

We discuss the suppression of a collective rhythm produced by a population of interactive oscillatory units by pulsatile stimulation. The problem is motivated by the medical technique known as deep brain stimulation. To imitate beta-band brain activity observed in Parkinsonian patients, we introduce time-dependent connectivity and consider a feedback-based approach for controlling this non-autonomous system. We pay attention to related technical problems like a causal real-time estimation of the phase and amplitude of a signal mixed with strong observational noise. We demonstrate that an adaptive control scheme can successfully suppress an undesired rhythm.

Parametrized state space models in the form of recurrent neural networks are often used in machine learning to learn from data streams exhibiting temporal dependencies. I will talk about a framework for rigorous analysis of such state space models known as Echo State Networks (ESN). In particular, I will show that there is a deep connection between the state space and the notion of feature space in kernel machines (a well known class of machine learning methods). This viewpoint will lead to several rather surprising results linking the structure of the dynamical coupling in ESN with properties of the associated temporal feature space. I will also present statements regarding universality of ESNs as capable of approximating a large class of fading memory filters with a minimal structural complexity.

The spectral radius R of the synaptic weight matrix determines whether the activity of a recurrent network dies out, or explodes, modulo non-linear effects, with an absorbing phase transition occurring at R=1. With the spectral radius one denotes the norm of the largest eigenvalue of the synaptic weight matrix, which is a global property. The question is then, whether there exists biologically plausible adaption rules for the online regulation of the spectral radius, where online means that adaption is functional in the non-autonomous regime when the network is actively performing tasks. We show that the spectral radius can be set using a rule which is based on the regulation of the local flow of activity. Neurons compare the variance of their own activity to the variance of the input signals received. Regulating the respective ratio by adapting their internal gain is shown to be equivalent to the regulation of the spectral radius whenever inter-neuron correlations are weak. Importantly, this principle works also in the presence of external inputs, viz in the non-autonomous regime. In addition we present an explicit analytic expression for absorbing phase transitions in the presence of an external inputs which becomes exact in the absence of spatial correlations.

In the era of the data science, model-free techniques are developed overwhelmingly. When the experimentally-collected data are generated by dynamical systems, the missions of reconstruction, prediction, and modulation only based on these data are highly anticipated to be achieved for these systems. Here, we introduce several directions of progresses made by our research group in developing the model-free techniques using machine learning techniques and dynamical systems theory. We use representative systems of physical or/and biological significance to demonstrate the developed techniques. We hope that all the methods can shed a light on deciphering and controlling the hidden dynamics that dominate the evolutions of any systems in real-world.

The human brain is an open, dissipative, and adaptive nonstationary dynamical system composed of a large number of interacting subsystems. Its complicated spatial–temporal dynamics is still poorly understood. Epilepsy is a malfunction of the brain that affects about 50 million people worldwide. Epileptic seizures are the cardinal symptom of the disease and often appear as a transformation of otherwise normal brain dynamics. The exact mechanisms underlying seizure generation are still as uncertain as are mechanisms underlying seizure spreading and termination. Identifying early warning signals of seizures from brain dynamics could drastically improve therapeutic possibilities and thus, the quality of life of people with epilepsy. Methods from nonlinear dynamics, statistical physics, synchronization and network theory are capable of identifying seizure precursors from EEG recordings in a large number of people with epilepsy and with high sensitivity and specificity. In this talk, I will provide a brief overview of the current status of the field, will highlight shortcomings of recent approaches as well as unsovled issues, and will discuss possible research directions that may help to find better solutions.

Motion of stars is usually considered in time independent galactic potential. However, for various circumstances (e.g. encounters between the individual stars; stars of galaxies and globular clusters lose substantial quantities of mass) one have to take into account potential variations that are slow (or adiabatic) compared to a typical orbital frequency. In this case the action variables of stars are constant during such adiabatic changes of the potential. For this reason actions are often called adiabatic invariants. The nature of the problem posed by a time-varying potential depends on the speed with which the potential evolves. In slowly evolving potentials, angle-action variables enable us to predict how a stellar system will respond to changes in the gravitational field that confines it. We seek in our research a novel approach to describe the dynamics beyond the adiabatic limit, i.e. when the change of potential is comparable to the characteristic time scale of the problem. To our knowledge this is the first attempt to use contemporary results of nonautonomous Hamiltonian dynamics in celestial mechanics. We hope that this idea will enable us to explore the peculiar nature of stellar dynamics subjected to external forces.

The Atlantic Meridional Overturning Circulation (AMOC) is believed to possess two stable flow regimes, thus qualifying as a proposed tipping element in the Earth system. Abrupt climate changes in paleo records are often associated with shifts between different AMOC flow patterns. To better understand the mechanisms and probability of such transitions, it is important to understand their transition pathway in the state space of climate models. Specifically, in weakly perturbed dynamical systems, it is expected that sample transitions will cross from one to another basin of attraction via special regions on the basin boundary, termed Melancholia states. Here we implement an edge-tracking algorithm in PlaSim-LSG, a global climate model of intermediate complexity, to identify a Melancholia state between the vigorous and collapsed states of the AMOC. We discuss the properties of such a state, its role for critical transitions of the AMOC, and its dependence on time-dependent forcing. Knowing what lies in between the competing metastable states could help elucidate the risk and most likely scenario of potential future AMOC tipping.

Syamal Kumar Dana Centre for Mathematical Biology and Ecology, Department of Mathematics, Jadavpur University, Kolkata 700032, India Abstract: We report tipping in an insect outbreak model near two saddle-node bifurcation points of the model system. The model defines outbreak of budworm population from a desirable low population state against the carrying capacity. If the carrying capacity is varied at a linear rate, the system does not show sharp transitions immediately at the bifurcation points but tips to the alternate states after a lapse of time. This delay in tipping decreases with faster rate of change of the carrying capacity. We check the impact of a shock (external effect) modelled by a triangular shape impulse that slowly varies the carrying capacity at a constant rate, however, withdrawn at an identical or different constant rate. The tipping occurs in such cases but shows a dependence on the falling and rising rates of the impulse that influences the carrying capacity. We identify the rate of rising and falling in a phase diagram that provides information when a tipping may occur. An interplay of the dissipation rate of the system and the rate parameter decides the occurrence of tipping. In case of insufficient strength of a single impulse that may fail to induce any tipping to the desirable low population state from a large population or an outbreak state, a second impulse is applied, but weaker in strength compared to the first one and it is able to induce a tipping that prevents an insect outbreak. Therefore, in a situation of multiple shocks or impulses applied on the carrying capacity, tipping may occur as a consequence of past environmental changes. We explain the scenarios with numerical experiments and using the dynamical change in the potential function of the system.

Abrupt and irreversible winter Arctic sea ice loss may occur under anthropogenic warming due to the disappearance of a sea ice equilibrium at a threshold value of CO$_2$, commonly referred to as a tipping point. Previous work has been unable to conclusively identify whether a tipping point in winter Arctic sea ice exists because fully coupled climate models are too computationally expensive to run to equilibrium for many CO2 values. Here, we explore the non-autonomous behavior of sea ice under realistic rates of CO2 increase to demonstrate for the first time how a few time-dependent CO2 experiments can be used to predict the existence and timing of sea ice tipping points without running the model to steady state. This work highlights the inefficacy of using a single experiment with slow-changing CO2 to discover changes in the sea-ice steady state and provides a novel alternate method that can be developed for the identification of tipping points in realistic climate models.

Several Earth system components have been suggested to exhibit alternative stable states, and for some paleoclimate evidence of abrupt transitions indeed indicates abrupt shifts between states in past climates. Well-known examples include the polar ice sheets, the Atlantic Meridional Overturning Circulation, or the Amazon rainforest. I will present some recent results on the stability of these systems from empirical evidence and model simulations, focussing on the application of concepts from dynamical systems theory.

An excitable low-order quasigeostrophic model [1] captures key features of intrinsic low-frequency variability of the oceans’ wind-driven circulation. This double-gyre model is used as a prototype of an excitable nonlinear dynamical system with time-dependent forcing to explore basic features of climate change in the presence of natural variability. The present studies rely on the theoretical framework of nonautonomous dynamical systems and of their pullback attractors (PBAs), namely the time-dependent invariant sets that attract all trajectories initialized in the remote past. Ensemble simulations help us explore these PBAs. In an excitable dynamical system, a relaxation oscillation connects a basic state to an unstable excited state, from which spontaneous returns to the original state occur. Such back-and-forth oscillations are self-sustained in a certain parameter range of the autonomous system, and they can be excited by a suitable external time-dependent forcing outside this range. In this presentation, by using model [1] in its excitable state, we discuss: (i) the effect of a periodic forcing on synchronization scenarios and on the onset of chaos in an otherwise non-chaotic system [2,3]; (ii) the effect of an aperiodic forcing mimicking a time dependence dominated by interdecadal and high-frequency variability [4,5]; (iii) the tipping points due to parameter drift in the presence of small periodic perturbations [6]. The latter situation can be thought of as the seasonal-to-interannual variability in the wind stress, while the monotonically increasing component stands for the effect of reduction in the midlatitude winds due to anthropogenic warming. References: [1] Pierini S., 2011. J. Phys. Oceanogr. 41, 1585-1604. [2] Pierini S., 2014. J. Phys. Oceanogr. 441, 3245-3254. [3] Pierini S., Chekroun M.D., Ghil M., 2018. Nonlin. Processes Geophys. 25, 671-692. [4] Pierini S., Ghil M., Chekroun M.D., 2016. J. Climate 29, 4185-4202. [5] Pierini S., 2020. J. Stat. Phys. 179, 1475-1495. [6] Pierini S., Ghil M., 2021. Sci. Rep. 11, 11126

Any dynamical system subjected to monotonous parameter drift can be considered to undergo its own climate change if the rate of the parameter change is not adiabatically slow. Such a situation calls for ensemble simulations, proposed in chaos theory with the appearance of snapshot attractors, as early as 1990. In systems exhibiting trends, basic methods of standard chaos theory are not applicable: the efficient tool of periodic orbit expansion cannot be used since periodic orbits do not exist. Furthermore, long-time limits are ill-defined since the system might become qualitatively different from the original one even after short times. The talk deals with the question of how to identify chaos in such systems. From the point of view of phase space structures, we argue that stable and unstable foliations are easy to generate numerically, without relying on the existence of hyperbolic periodic orbits. Chaos originates from a Smale horseshoe-like pattern of the foliations, the transverse intersections of which indicate a chaotic set changing in time. We discuss in which sense such a process is related to transient chaos and time-dependent chaotic saddles. In dissipative cases, the unstable foliation is found to be part of the chaotic snapshot attractor, while the stable one turns out to be related to the basin of attraction of the time-dependent chaos. Concerning a quantitative characterization, the so-called ensemble-averaged pairwise distance is shown to provide a generalization of the concept of Lyapunov exponents, since the slope of this function can be interpreted as an instantaneous Lyapunov exponent. This is a tool by means of which a strengthening, weaking, or even disappearance of chaos can be investigated. D. Jánosi, T. Tél: Characterizing chaos in systems subjected to parameter drift, Phys. Rev. E 105, L062202 (2022)

Over the last two decades, tipping points have become a hot topic due to the devastating consequences that they may have on natural and human systems. Tipping points are typically associated with a system bifurcation when external forcing crosses a critical level, causing an abrupt transition to an alternative, and often less desirable, state. However, the rate of change in forcing is arguably of even greater relevance in the human-dominated anthropocene, but is rarely examined as a potential sole mechanism for tipping points. Thus, I will introduce the related phenomenon of rate-induced tipping: an instability that occurs when external forcing varies across some critical rate, usually without crossing any bifurcations. First, I will explain when to expect rate-induced tipping. Then, using illustrating examples of differing complexity I will highlight universal and generic properties of rate-induced tipping in a range of natural and human systems.

The effect of climate warming on species' physiological parameters, including growth rate, mortality rate, and handling time, is well established from empirical data. However, with an alarming rise in global temperature more than ever, predicting the interactive influence of these changes on mutualistic communities remains uncertain. Using 139 real plant–pollinator networks sampled across the globe and a modeling approach, we study the impact of species’ individual thermal responses on mutualistic communities. We show that at low mutualistic strength plant–pollinator networks are at potential risk of rapid transitions at higher temperatures. Evidently, generalist species play a critical role in guiding tipping points in mutualistic networks. Further, we derive stability criteria for the networks in a range of temperatures using a two-dimensional reduced model. We identify network structures that can ascertain the delay of a community collapse. Until the end of this century, on account of increasing climate warming many real mutualistic networks are likely to be under the threat of sudden collapse, and we frame strategies to mitigate this. Together, our results indicate that knowing individual species' thermal responses and network structure can improve predictions for communities facing rapid transitions.

In certain drylands around the globe vegetation is distributed in bands that are regularly spaced and oriented transverse to a gentle slope. The limited water resource ends up concentrated in the bands via feedback mechanisms that act on the fast timescales of the rare storms. In our modeling framework, we treat the storms as kicks that add water to the soil in a heterogeneous fashion that depends on the biomass distribution. These kicks then determine the initialization of a flow map that captures the slow dynamics of the vegetation distribution during the long dry periods between the storms. This modeling framework is designed so we can explore possible impacts of storm variability. In particular we introduce randomness to the kick strength and the subsequent flow duration. We are interested in how these random aspects of storm strength and timing impact the sustainability of these fragile ecosystems under climate change.

We investigate the phenomenon of extreme transient dynamics in four pendula coupled via the movable support. Our findings uncover a traveling phase state, within which the oscillators synchronize with phases evolving. The scenarios of possible traveling are discussed and the properties of patterns are determined. We report on the problem of transient dynamics, which can be observed on the road to final network synchronization, exhibiting that the accompanied transient windows can become extremely long. The behaviours are carefully studied in relation with initial conditions and parameters, showing how the system characteristics influence the transients. The results of our research underline the problem of slow dynamics evolution and its relation with mechanical setups.

I will present a new concept "dynamics-based data science" in biology and medicine for studying dynamical processes and disease progressions, including dynamic network biomarkers (DNB) for early-warning signals of critical transitions, spatial-temporal information (STI) transformation for short-term time-series prediction, and partial cross-mapping (PCM) for causal inference among variables. These methods are all data-driven or model-free approaches but based on the theoretical frameworks of nonlinear dynamics. We show the principles and advantages of dynamics-based data-driven approaches as explicable, quantifiable, and generalizable. In particular, dynamics-based data science approaches exploit the essential features of dynamical systems in terms of data, e.g. critical collective fluctuation near a tipping point, low-dimensionality of a center manifold or an attractor, and phase-space reconstruction from a single variable by delay embedding theorem, and thus are able to provide different or additional information to the traditional approaches, i.e. statistics-based data science approaches. The dynamical data science approaches will further play an important role in the systematical research of various fields in biology and medicine as well as AI.

Nonlinear dynamics teaches us that small stressors can lead to much larger catastrophes and even the collapse of entire complex systems. Modern global human civilisation is precisely such a system, and the stressors --- such as global climate change --- are not small. Catastrophic scenarios such as local or global societal collapse are increasingly being studied, yet essential ideas from nonlinear dynamics have not yet been incorporated. One of these ideas is rate-induced tipping. Systems break down when external changes overcome the capacity for adaptation; in other words, when things change too quickly. Using evidence from nature, dynamical systems, and agent-based models, I argue that rate-induced collapse is ubiquitous across complex systems. If true, rate-induced tipping is of fundamental relevance to some of the most consequential problems facing humanity this century. There is much exciting work to be done.