Source author record

Davide Faranda

Davide Faranda appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
18topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

23 published item(s)

preprint2022arXiv

Stochastic Chaos and Predictability in Laboratory Earthquakes

Laboratory earthquakes exhibit characteristics of a low dimensional random attractor with a dimension similar to that of natural slow earthquakes. A model of stochastic differential equations based on rate and state-dependent friction explains the laboratory observations. We study the transition from stable sliding to stickslip events and find that aperiodic behavior can be explained by small perturbations in the stress state. Friction's nonlinear nature amplifies small scale perturbations, reducing the predictability of the otherwise periodic macroscopic dynamics.

preprint2020arXiv

Detecting regime transitions of the nocturnal and Polar near-surface temperature inversion

Many natural systems undergo critical transitions, i.e. sudden shifts from one dynamical regime to another. In the climate system, the atmospheric boundary layer can experience sudden transitions between fully turbulent states and quiescent, quasi-laminar states. Such rapid transitions are observed in Polar regions or at night when the atmospheric boundary layer is stably stratified, and they have important consequences in the strength of mixing with the higher levels of the atmosphere. To analyze the stable boundary layer, many approaches rely on the identification of regimes that are commonly denoted as weakly and very stable regimes. Detecting transitions between the regimes is crucial for modeling purposes. In this work a combination of methods from dynamical systems and statistical modeling is applied to study these regime transitions and to develop an early-warning signal that can be applied to non-stationary field data. The presented metric aims at detecting nearing transitions by statistically quantifying the deviation from the dynamics expected when the system is close to a stable equilibrium. An idealized stochastic model of near-surface inversions is used to evaluate the potential of the metric as an indicator of regime transitions. In this stochastic system, small-scale perturbations can be amplified due to the nonlinearity, resulting in transitions between two possible equilibria of the temperature inversion. The simulations show such noise-induced regime transitions, successfully identified by the indicator. The indicator is further applied to time series data from nocturnal and Polar meteorological measurements.

preprint2020arXiv

On the uncertainty of real-time predictions of epidemic growths: a COVID-19 case study for China and Italy

While COVID-19 is rapidly propagating around the globe, the need for providing real-time forecasts of the epidemics pushes fits of dynamical and statistical models to available data beyond their capabilities. Here we focus on statistical predictions of COVID-19 infections performed by fitting asymptotic distributions to actual data. By taking as a case-study the epidemic evolution of total COVID-19 infections in Chinese provinces and Italian regions, we find that predictions are characterized by large uncertainties at the early stages of the epidemic growth. Those uncertainties significantly reduce after the epidemics peak is reached. Differences in the uncertainty of the forecasts at a regional level can be used to highlight the delay in the spread of the virus. Our results warn that long term extrapolation of epidemics counts must be handled with extreme care as they crucially depend not only on the quality of data, but also on the stage of the epidemics, due to the intrinsically non-linear nature of the underlying dynamics. These results suggest that real-time epidemiological projections should include wide uncertainty ranges and urge for the needs of compiling high-quality datasets of infections counts, including asymptomatic patients.

preprint2016arXiv

Extremes and Recurrence in Dynamical Systems

This book provides a comprehensive introduction for the study of extreme events in the context of dynamical systems. The introduction provides a broad overview of the interdisciplinary research area of extreme events, underlining its relevance for mathematics, natural sciences, engineering, and social sciences. After exploring the basics of the classical theory of extreme events, the book presents a careful examination of how a dynamical system can serve as a generator of stochastic processes, and explores in detail the relationship between the hitting and return time statistics of a dynamical system and the possibility of constructing extreme value laws for given observables. Explicit derivation of extreme value laws are then provided for selected dynamical systems. The book then discusses how extreme events can be used as probes for inferring fundamental dynamical and geometrical properties of a dynamical system and for providing a novel point of view in problems of physical and geophysical relevance. A final summary of the main results is then presented along with a discussion of open research questions. Finally, an appendix with software in Matlab programming language allows the readers to develop further understanding of the presented concepts.

preprint2015arXiv

A non equilibrium Ising model of turbulence

We introduce a model of interacting lattices at different resolutions driven by the two-dimensional Ising dynamics with a nearest-neighbor interaction. We study this model both with tools borrowed from equilibrium statistical mechanics as well as non-equilibrium thermodynamics. Our findings show that this model keeps the signature of the equilibrium phase transition. Moreover the critical temperature of the equilibrium models correspond to the state maximizing the entropy and delimits two out-of-equilibrium regimes, one satisfying the Onsager relations for systems close to equilibrium and one resembling convective turbulent states. Since the model preserves the entropy and energy fluxes in the scale space, it seems a good candidate for parametric studies of out-of-equilibrium turbulent systems.

preprint2015arXiv

Global vs local energy dissipation: the energy cycle of the turbulent von Kármán flow

In this paper, we investigate the relations between global and local energy transfers in a turbulent von Kármán flow. The goal is to understand how and where energy is dissipated in such a flow and to reconstruct the energy cycle in an experimental device where local as well as global quantities can be measured. We use PIV measurements and we model the Reynolds stress tensor to take subgrid scales into account. This procedure involves a free parameter that is calibrated using angular momentum balance. We then estimate the local and global mean injected and dissipated power for several types of impellers, for various Reynolds numbers and for various flow topologies. These PIV estimates are then compared with direct injected power estimates provided by torque measurements at the impellers. The agreement between PIV estimates and direct measurements depends on the flow topology. In symmetric situations, we are able to capture up to 90% of the actual global energy dissipation rate. However, our results become increasingly inaccurate as the shear layer responsible for most of the dissipation approaches one of the impellers, and cannot be resolved by our PIV set-up. Finally, we show that a very good agreement between PIV estimates and direct measurements is obtained using a new method based on the work of Duchon and Robert which generalizes the Kármán-Howarth equation to nonisotropic, nonhomogeneous flows. This method provides parameter-free estimates of the energy dissipation rate as long as the smallest resolved scale lies in the inertial range. These results are used to evidence a well-defined stationary energy cycle within the flow in which most of the energy is injected at the top and bottom impellers, and dissipated within the shear layer. The influence of the mean flow geometry and the Reynolds number on this energy cycle is studied for a wide range of parameters.

preprint2015arXiv

Maximum Kolmogorov-Sinai entropy vs minimum mixing time in Markov chains

Many modern techniques employed in physics, such a computation of path integrals, rely on random walks on graphs that can be represented as Markov chains. Traditionally, estimates of running times of such sampling algorithms are computed using the number of steps in the chain needed to reach the stationary distribution. This quantity is generally defined as mixing time and is often difficult to compute. In this paper, we suggest an alternative estimate based on the Kolmogorov-Sinai entropy, by establishing a link between the maximization of KSE and the minimization of the mixing time. Since KSE are easier to compute in general than mixing time, this link provides a new faster method to approximate the minimum mixing time that could be interesting in computer sciences and statistical physics. Beyond this, our finding will also be of interest to the out-of-equilibrium community, by providing a new rational to select stationary states in out-of-equilibrium physics: it seems reasonable that in a physical system with two simultaneous equiprobable possible dynamics, the final stationary state will be closer to the stationary state corresponding to the fastest dynamics (smallest mixing time).Through the empirical link found in this letter, this state will correspond to a state of maximal Kolmogorov-Sinai entropy. If this is true, this would provide a more satisfying rule for selecting stationary states in complex systems such as climate than the maximization of the entropy production.

preprint2015arXiv

Statistical optimization for passive scalar transport: maximum entropy production vs maximum Kolmogorov-Sinay entropy

We derive rigorous results on the link between the principle of maximum entropy production and the principle of maximum Kolmogorov-Sinai entropy using a Markov model of the passive scalar diffusion called the Zero Range Process. We show analytically that both the entropy production and the Kolmogorov-Sinai entropy seen as functions of f admit a unique maximum denoted fmaxEP and fmaxKS. The behavior of these two maxima is explored as a function of the system disequilibrium and the system resolution N. The main result of this article is that fmaxEP and fmaxKS have the same Taylor expansion at _rst order in the deviation of equilibrium. We find that fmaxEP hardly depends on N whereas fmaxKS depends strongly on N. In particular, for a fixed difference of potential between the reservoirs, fmaxEP (N) tends towards a non-zero value, while fmaxKS (N) tends to 0 when N goes to infinity. For values of N typical of that adopted by Paltridge and climatologists we show that fmaxEP and fmaxKS coincide even far from equilibrium. Finally, we show that one can find an optimal resolution N_ such that fmaxEP and fmaxKS coincide, at least up to a second order parameter proportional to the non-equilibrium uxes imposed to the boundaries.

preprint2014arXiv

A recurrence-based technique for detecting genuine extremes in instrumental temperature records

In this paper, we analyze several instrumental records of temperatures at different locations by using new techniques originally developed for the analysis of extreme values of dynamical systems. We show that they have the same recurrence time statistics as a chaotic dynamical system perturbed with dynamical noise and by instrument errors. The technique provides a criterion to discriminate whether the recurrence of a certain temperature belongs to the normal variability or can be considered as a genuine extreme event with respect to a specific timescale fixed as parameter. The method gives a self-consistent estimation of the convergence of the statistics of recurrences toward the theoretical extreme value laws.

preprint2014arXiv

Extreme Value laws for dynamical systems under observational noise

In this paper we prove the existence of Extreme Value Laws for dynamical systems perturbed by instrument-like-error, also called observational noise. An orbit perturbed with observational noise mimics the behavior of an instrumentally recorded time series. Instrument characteristics - defined as precision and accuracy - act both by truncating and randomly displacing the real value of a measured observable. Here we analyze both these effects from a theoretical and numerical point of view. First we show that classical extreme value laws can be found for orbits of dynamical systems perturbed with observational noise. Then we present numerical experiments to support the theoretical findings and give an indication of the order of magnitude of the instrumental perturbations which cause relevant deviations from the extreme value laws observed in deterministic dynamical systems. Finally, we show that the observational noise preserves the structure of the deterministic attractor. This goes against the common assumption that random transformations cause the orbits asymptotically fill the ambient space with a loss of information about any fractal structures present on the attractor.

preprint2014arXiv

Mixing properties in the advection of passive tracers via recurrences and extreme value theory

In this paper we characterize the mixing properties in the advection of passive tracers by exploiting the extreme value theory for dynamical systems. With respect to classical techniques directly related to the Poincaré recurrences analysis, our method provides reliable estimations of the characteristic mixing times and distinguishes between barriers and unstable fixed points. The method is based on a check of convergence for extreme value laws on finite datasets. We define the mixing times in terms of the shortest time intervals such that extremes converge to the asymptotic (known) parameters of the Generalized Extreme Value distribution. Our technique is suitable for applications in the analysis of other systems where mixing time scales need to be determined and limited datasets are available.

preprint2014arXiv

On using Extreme Values to detect global stability thresholds in multi-stable systems: The case of transitional plane Couette flow

Extreme Value Theory (EVT) is exploited to determine the global stability threshold $R_g$ of plane Couette flow --the flow of a viscous fluid in the space between two parallel plates-- whose laminar or turbulent behavior depends on the Reynolds number R. Even if the existence of a global stability threshold has been detected in simulations and experiments, its numerical value has not been unequivocally defined. $R_g$ is the value such that for $R > R_g$ turbulence is sustained, whereas for $R < R_g$ it is transient and eventually decays. We address the problem of determining $R_g$ by using the extremes - maxima and minima - of the perturbation energy fluctuations. When $R_g$, both the positive and negative extremes are bounded. As the critical Reynolds number is approached from above, the probability of observing a very low minimum increases causing asymmetries in the distributions of maxima and minima. On the other hand, the maxima distribution is unaffected as the fluctuations towards higher values of the perturbation energy remain bounded. This tipping point can be detected by fitting the data to the Generalized Extreme Value (GEV) distribution and by identifying $R_g$ as the value of $R$ such that the shape parameter of the GEV for the minima changes sign from negative to positive. The results are supported by the analysis of theoretical models which feature a bistable behavior.

preprint2014arXiv

Probing turbulence intermittency via Auto-Regressive Moving-Average models

We suggest a new approach to probing intermittency corrections to the Kolmogorov law in turbulent flows based on the Auto-Regressive Moving-Average modeling of turbulent time series. We introduce a new index $Υ$ that measures the distance from a Kolmogorov-Obukhov model in the Auto-Regressive Moving-Average models space. Applying our analysis to Particle Image Velocimetry and Laser Doppler Velocimetry measurements in a von Kármán swirling flow, we show that $Υ$ is proportional to the traditional intermittency correction computed from the structure function. Therefore it provides the same information, using much shorter time series. We conclude that $Υ$ is a suitable index to reconstruct the spatial intermittency of the dissipation in both numerical and experimental turbulent fields.

preprint2014arXiv

Robust estimate of dynamo thresholds in the von Kármán sodium experiment using the Extreme Value Theory

We apply a new threshold detection method based on the extreme value theory to the von Kármán sodium (VKS) experiment data. The VKS experiment is a successful attempt to get a dynamo magnetic field in a laboratory liquid-metal experiment. We first show that the dynamo threshold is associated to a change of the probability density function of the extreme values of the magnetic field. This method does not require the measurement of response functions from applied external perturbations, and thus provides a simple threshold estimate. We apply our method to different configurations in the VKS experiment showing that it yields a robust indication of the dynamo threshold as well as evidence of hysteretic behaviors. Moreover, for the experimental configurations in which a dynamo transition is not observed, the method provides a way to extrapolate an interval of possible threshold values.

preprint2014arXiv

Statistical early-warning indicators based on Auto-Regressive Moving-Average processes

We address the problem of defining early warning indicators of critical transition. To this purpose, we fit the relevant time series through a class of linear models, known as Auto-Regressive Moving-Average (ARMA(p,q)) models. We define two indicators representing the total order and the total persistence of the process, linked, respectively, to the shape and to the characteristic decay time of the autocorrelation function of the process. We successfully test the method to detect transitions in a Langevin model and a 2D Ising model with nearest-neighbour interaction. We then apply the method to complex systems, namely for dynamo thresholds and financial crisis detection.

preprint2013arXiv

Extreme value statistics for dynamical systems with noise

We study the distribution of maxima (Extreme Value Statistics) for sequences of observables computed along orbits generated by random transformations. The underlying, deterministic, dynamical system can be regular or chaotic. In the former case, we will show that by perturbing rational or irrational rotations with additive noise, an extreme value law appears, regardless of the intensity of the noise, while unperturbed rotations do not admit such limiting distributions. In the case of deterministic chaotic dynamics, we will consider observables specially designed to study the recurrence properties in the neighbourhood of periodic points. Hence, the exponential limiting law for the distribution of maxima is modified by the presence of the extremal index, a positive parameter not larger than one, whose inverse gives the average size of the clusters of extreme events. The theory predicts that such a parameter is unitary when the system is perturbed randomly. We perform sophisticated numerical tests to assess how strong is the impact of noise level, when finite time series are considered. We find agreement with the asymptotic theoretical results but also non-trivial behaviour in the finite range. In particular our results suggest that in many applications where finite datasets can be produced or analysed one must be careful in assuming that the smoothing nature of noise prevails over the underlying deterministic dynamics.

preprint2013arXiv

Towards a General Theory of Extremes for Observables of Chaotic Dynamical Systems

In this paper we provide a connection between the geometrical properties of a chaotic dynamical system and the distribution of extreme values. We show that the extremes of so-called physical observables are distributed according to the classical generalised Pareto distribution and derive explicit expressions for the scaling and the shape parameter. In particular, we derive that the shape parameter does not depend on the chosen observables, but only on the partial dimensions of the invariant measure on the stable, unstable, and neutral manifolds. The shape parameter is negative and is close to zero when high-dimensional systems are considered. This result agrees with what was derived recently using the generalized extreme value approach. Combining the results obtained using such physical observables and the properties of the extremes of distance observables, it is possible to derive estimates of the partial dimensions of the attractor along the stable and the unstable directions of the flow. Moreover, by writing the shape parameter in terms of moments of the extremes of the considered observable and by using linear response theory, we relate the sensitivity to perturbations of the shape parameter to the sensitivity of the moments, of the partial dimensions, and of the Kaplan-Yorke dimension of the attractor. Preliminary numerical investigations provide encouraging results on the applicability of the theory presented here. The results presented here do not apply for all combinations of Axiom A systems and observables, but the breakdown seems to be related to very special geometrical configurations.

preprint2012arXiv

Analysis of Round Off Errors with Reversibility Test as a Dynamical Indicator

We compare the divergence of orbits and the reversibility error for discrete time dynamical systems. These two quantities are used to explore the behavior of the global error induced by round off in the computation of orbits. The similarity of results found for any system we have analysed suggests the use of the reversibility error, whose computation is straightforward since it does not require the knowledge of the exact orbit, as a dynamical indicator. The statistics of fluctuations induced by round off for an ensemble of initial conditions has been compared with the results obtained in the case of random perturbations. Significant differences are observed in the case of regular orbits due to the correlations of round off error, whereas the results obtained for the chaotic case are nearly the same. Both the reversibility error and the orbit divergence computed for the same number of iterations on the whole phase space provide an insight on the local dynamical properties with a detail comparable with other dynamical indicators based on variational methods such as the finite time maximum Lyapunov characteristic exponent, the mean exponential growth factor of nearby orbits and the smaller alignment index. For 2D symplectic maps the differentiation between regular and chaotic regions is well full-filled. For 4D symplectic maps the structure of the resonance web as well as the nearby weakly chaotic regions are accurately described.

preprint2011arXiv

Bistable systems with Stochastic Noise: Virtues and Limits of effective Langevin equations for the Thermohaline Circulation strength

The understanding of the statistical properties and of the dynamics of multistable systems is gaining more and more importance in a vast variety of scientific fields. This is especially relevant for the investigation of the tipping points of complex systems. Sometimes, in order to understand the time series of given observables exhibiting bimodal distributions, simple one-dimensional Langevin models are fitted to reproduce the observed statistical properties, and used to investing-ate the projected dynamics of the observable. This is of great relevance for studying potential catastrophic changes in the properties of the underlying system or resonant behaviours like those related to stochastic resonance-like mechanisms. In this paper, we propose a framework for encasing this kind of studies and show, using simple box models of the oceanic circulation and choosing as observable the strength of the thermohaline circulation. We study the statistical properties of the transitions between the two modes of operation of the thermohaline circulation under symmetric boundary forcing and test their agreement with simplified one-dimensional phenomenological theories. We extend our analysis to include stochastic resonance-like amplification processes. We conclude that fitted one-dimensional Langevin models, when closely scrutinised, may result to be more ad-hoc than they seem, lacking robustness and/or well-posedness. They should be treated with care, more as an empiric descriptive tool than as methodology with predictive power.

preprint2011arXiv

Extreme Value distribution for singular measures

In this paper we perform an analytical and numerical study of Extreme Value distributions in discrete dynamical systems that have a singular measure. Using the block maxima approach described in Faranda et al. [2011] we show that, numerically, the Extreme Value distribution for these maps can be associated to the Generalised Extreme Value family where the parameters scale with the information dimension. The numerical analysis are performed on a few low dimensional maps. For the middle third Cantor set and the Sierpinskij triangle obtained using Iterated Function Systems, experimental parameters show a very good agreement with the theoretical values. For strange attractors like Lozi and Hènon maps a slower convergence to the Generalised Extreme Value distribution is observed. Even in presence of large statistics the observed convergence is slower if compared with the maps which have an absolute continuous invariant measure. Nevertheless and within the uncertainty computed range, the results are in good agreement with the theoretical estimates.

preprint2011arXiv

Numerical convergence of the block-maxima approach to the Generalized Extreme Value distribution

In this paper we perform an analytical and numerical study of Extreme Value distributions in discrete dynamical systems. In this setting, recent works have shown how to get a statistics of extremes in agreement with the classical Extreme Value Theory. We pursue these investigations by giving analytical expressions of Extreme Value distribution parameters for maps that have an absolutely continuous invariant measure. We compare these analytical results with numerical experiments in which we study the convergence to limiting distributions using the so called block-maxima approach, pointing out in which cases we obtain robust estimation of parameters. In regular maps for which mixing properties do not hold, we show that the fitting procedure to the classical Extreme Value Distribution fails, as expected. However, we obtain an empirical distribution that can be explained starting from a different observable function for which Nicolis et al. [2006] have found analytical results.

preprint2011arXiv

Relevance of sampling schemes in light of Ruelle's linear response theory

We reconsider the theory of the linear response of non-equilibrium steady states to perturbations. We first show that by using a general functional decomposition for space-time dependent forcings, we can define elementary susceptibilities that allow to construct the response of the system to general perturbations. Starting from the definition of SRB measure, we then study the consequence of taking different sampling schemes for analysing the response of the system. We show that only a specific choice of the time horizon for evaluating the response of the system to a general time-dependent perturbation allows to obtain the formula first presented by Ruelle. We also discuss the special case of periodic perturbations, showing that when they are taken into consideration the sampling can be fine-tuned to make the definition of the correct time horizon immaterial. Finally, we discuss the implications of our results in terms of strategies for analyzing the outputs of numerical experiments by providing a critical review of a formula proposed by Reick.

preprint2011arXiv

Universal behavior of extreme value statistics for selected observables of dynamical systems

The main results of the extreme value theory developed for the investigation of the observables of dynamical systems rely, up to now, on the Gnedenko approach. In this framework, extremes are basically identified with the block maxima of the time series of the chosen observable, in the limit of infinitely long blocks. It has been proved that, assuming suitable mixing conditions for the underlying dynamical systems, the extremes of a specific class of observables are distributed according to the so called Generalized Extreme Value (GEV) distribution. Direct calculations show that in the case of quasi-periodic dynamics the block maxima are not distributed according to the GEV distribution. In this paper we show that, in order to obtain a universal behaviour of the extremes, the requirement of a mixing dynamics can be relaxed if the Pareto approach is used, based upon considering the exceedances over a given threshold. Requiring that the invariant measure locally scales with a well defined exponent - the local dimension -, we show that the limiting distribution for the exceedances of the observables previously studied with the Gnedenko approach is a Generalized Pareto distribution where the parameters depends only on the local dimensions and the value of the threshold. This result allows to extend the extreme value theory for dynamical systems to the case of regular motions. We also provide connections with the results obtained with the Gnedenko approach. In order to provide further support to our findings, we present the results of numerical experiments carried out considering the well-known Chirikov standard map.