BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Valerio Lucarini

Valerio Lucarini contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechphysics.ao-phnlin.CDphysics.data-anphysics.geo-phphysics.flu-dynphysics.comp-phastro-ph.EPphysics.class-phphysics.hist-phphysics.soc-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2022arXiv

Decomposing the Dynamics of the Lorenz 1963 model using Unstable Periodic Orbits: Averages, Transitions, and Quasi-Invariant Sets

Unstable periodic orbits (UPOs) are a valuable tool for studying chaotic dynamical systems, as they allow one to distill their dynamical structure. We consider here the Lorenz 1963 model with the classic parameters' value. We investigate how a chaotic trajectory can be approximated using a complete set of UPOs up to symbolic dynamics' period 14. At each instant, we rank the UPOs according to their proximity to the position of the orbit in the phase space. We study this process from two different perspectives. First, we find that longer period UPOs overwhelmingly provide the best local approximation to the trajectory. Second, we construct a finite-state Markov chain by studying the scattering of the orbit between the neighbourhood of the various UPOs. Each UPO and its neighbourhood are taken as a possible state of the system. Through the analysis of the subdominant eigenvectors of the corresponding stochastic matrix we provide a different interpretation of the mixing processes occurring in the system by taking advantage of the concept of quasi-invariant sets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lévy-noise versus Gaussian-noise-induced Transitions in the Ghil-Sellers Energy Balance Model

We study the impact of applying stochastic forcing to the Ghil-Sellers energy balance climate model in the form of a fluctuating solar irradiance. Through numerical simulations, we explore the noise-induced transitions between the competing warm and snowball climate states. We consider multiplicative stochastic forcing driven by Gaussian and $α$-stable Lévy - $α\in(0,2)$ - noise laws, examine the statistics of transition times, and estimate most probable transition paths. While the Gaussian noise case has been carefully studied in a plethora of investigations on metastable systems, much less is known about the Lévy case, especially in the case of high- and infinite-dimensional systems. In the weak noise limit, the expected residence time in each metastable state scales in a fundamentally different way in the Gaussian vs. Lévy noise case with respect to the intensity of the noise. In the former case, the classical Kramers-like exponential law is recovered. In the latter case, power laws are found, with the exponent equal to $-α$, in apparent agreement with rigorous results obtained for additive noise in a related - yet different - reaction-diffusion equation as well as in simpler models. This can be better understood by treating the Lévy noise as a compound Poisson process. The transition paths are studied in a projection of the state space and remarkable differences are observed between the two different types of noise. The snowball-to-warm and the warm-to-snowball most probable transition path cross at the single unstable edge state on the basin boundary. In the case of Lévy noise, the most probable transition paths in the two directions are wholly separated, as transitions apparently take place via the closest basin boundary region to the outgoing attractor. This property can be better elucidated by considering singular perturbations to the solar irradiance.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Non-equilibrium Ensembles for the three-dimensional Navier-Stokes Equations

At the molecular level fluid motions are, by first principles, described by time reversible laws. On the other hand, the coarse grained macroscopic evolution is suitably described by the Navier-Stokes equations, which are inherently irreversible, due to the dissipation term. Here, a reversible version of three-dimensional Navier-Stokes is studied, by introducing a fluctuating viscosity constructed in such a way that enstrophy is conserved, along the lines of the paradigm of microcanonical versus canonical treatment in equilibrium statistical mechanics. Through systematic simulations we attack two important questions: (a) What are the conditions that must be satisfied in order to have a statistical equivalence between the two non-equilibrium ensembles? (b) What is the empirical distribution of the fluctuating viscosity observed by changing the Reynolds number and the number of modes used in the discretization of the evolution equation? The latter point is important also to establish regularity conditions for the reversible equations. We find that the probability to observe negative values of the fluctuating viscosity becomes very quickly extremely small when increasing the effective Reynolds number of the flow in the fully resolved hydro dynamical regime, at difference from what was observed previously.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Robustness of competing climatic states

The climate is a non-equilibrium system undergoing the continuous action of forcing and dissipation. Under the effect of a spatially inhomogeneous absorption of solar energy, all the climate components dynamically respond until an approximate steady state (or attractor) is reached. However, multiple steady states can co-exist for a given forcing and with the same boundary conditions. Here, we apply the Thermodynamic Diagnostic Tool (TheDiaTo) to investigate the statistical properties of five co-existing climates, ranging from a snowball to an ice-free aquaplanet, obtained in MITgcm coupled simulations. The aim is to explore the multistability of the climate model setup by highlighting differences in competing steady states and their characteristic signatures regarding the meridional transport of heat and water mass, the Lorenz energy cycle and the material entropy production. We also investigate how such attractors change when the model configuration is varied. We consider, in particular, the effect of changing the representation of the cloud albedo, and of implementing an improved closure of the energy budget. We find that, even if the dynamics remains on the same attractor, state variables are modified. The set of metrics in TheDiaTo quantifies such modifications and represents a valuable tool for model evaluation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Typicality of the 2021 Western North America Summer Heatwave

Elucidating the statistical properties of extreme meteo-climatic events and capturing the physical processes responsible for their occurrence are key steps for improving our understanding of climate variability and climate change and for better evaluating the associated hazards. It has recently become apparent that large deviation theory is very useful for investigating persistent extreme events, and specifically, for flexibly estimating long return periods and for introducing a notion of dynamical typicality. Using a methodological framework based on large deviation theory and taking advantage of long simulations by a state-of-the-art Earth System Model, we investigate the 2021 North America Heatwave. Indeed, our analysis shows that the 2021 event can be seen as an unlikely but possible manifestation of climate variability, whilst its probability of occurrence is greatly amplified by the ongoing climate change. We also clarify the properties of spatial coherence of the 2021 heatwave and elucidate the role played by the Rocky Mountains in modulating hot, dry, and persistent extreme events in the Western Pacific region of North America.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Eddy saturation in a reduced two-level model of the atmosphere

Eddy saturation describes the nonlinear mechanism in geophysical flows whereby, when average conditions are considered, direct forcing of the zonal flow increases the eddy kinetic energy, while the energy associated with the zonal flow does not increase. Here we present a minimal baroclinic model that exhibits complete eddy saturation. Starting from Phillips' classical quasi-geostrophic two-level model on the beta channel of the mid-latitudes, we derive a reduced order model comprising of six ordinary differential equations including parameterised eddies. This model features two physically realisable steady state solutions, one a purely zonal flow and one where, additionally, finite eddy motions are present. As the baroclinic forcing in the form of diabatic heating is increased, the zonal solution loses stability and the eddy solution becomes attracting. After this bifurcation, the zonal components of the solution are independent of the baroclinic forcing, and the excess of heat in the low latitudes is efficiently transported northwards by finite eddies, in the spirit of baroclinic adjustment.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Inferring the instability of a dynamical system from the skill of data assimilation exercises

Data assimilation (DA) aims at optimally merging observational data and model outputs to create a coherent statistical and dynamical picture of the system under investigation. Indeed, DA aims at minimizing the effect of observational and model error, and at distilling the correct ingredients of its dynamics. DA is of critical importance for the analysis of systems featuring sensitive dependence on the initial conditions, as chaos wins over any finitely accurate knowledge of the state of the system, even in absence of model error. Clearly, the skill of DA is guided by the properties of dynamical system under investigation, as merging optimally observational data and model outputs is harder when strong instabilities are present. In this paper we reverse the usual angle on the problem and show that it is indeed possible to use the skill of DA to infer some basic properties of the tangent space of the system, which may be hard to compute in very high-dimensional systems. Here, we focus our attention on the first Lyapunov exponent and the Kolmogorov-Sinai entropy, and perform numerical experiments on the Vissio-Lucarini 2020 model, a recently proposed generalisation of the Lorenz 1996 model that is able to describe in a simple yet meaningful way the interplay between dynamical and thermodynamical variables.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Global Stability Properties of the Climate: Melancholia States, Invariant Measures, and Phase Transitions

For a wide range of values of the incoming solar radiation, the Earth features at least two attracting states, which correspond to competing climates. The warm climate is analogous to the present one; the snowball climate features global glaciation and conditions that can hardly support life forms. Paleoclimatic evidences suggest that in the past our planet flipped between these two states. The main physical mechanism responsible for such an instability is the ice-albedo feedback. In a previous work, we defined the Melancholia states that sit between the two climates. Such states are embedded in the boundaries between the two basins of attraction and feature extensive glaciation down to relatively low latitudes. Here, we explore the global stability properties of the system by introducing random perturbations as modulations to the intensity of the incoming solar radiation. We observe noise-induced transitions between the competing basins of attraction. In the weak noise limit, large deviation laws define the invariant measure, the statistics of escape times, and typical escape paths called instantons. By constructing the instantons empirically, we show that the Melancholia states are the gateways for the noise-induced transitions. In the region of multistability, in the zero-noise limit, the measure is supported only on one of the competing attractors. For low (high) values of the solar irradiance, the limit measure is the snowball (warm) climate. The changeover between the two regimes corresponds to a first-order phase transition in the system. The framework we propose seems of general relevance for the study of complex multistable systems. Finally, we put forward a new method for constructing Melancholia states from direct numerical simulations, black which provides a possible alternative with respect to the edge-tracking algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Introduction to the Special Issue on the Statistical Mechanics of Climate

We introduce the special issue on the Statistical Mechanics of Climate published on the Journal of Statistical Physics by presenting an informal discussion of some theoretical aspects of climate dynamics that make it a topic of great interest for mathematicians and theoretical physicists. In particular, we briefly discuss its nonequilibrium and multiscale properties, the relationship between natural climate variability and climate change, the different regimes of climate response to perturbations, and critical transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Physics of Climate Variability and Climate Change

The climate system is a forced, dissipative, nonlinear, complex and heterogeneous system that is out of thermodynamic equilibrium. The system exhibits natural variability on many scales of motion, in time as well as space, and it is subject to various external forcings, natural as well as anthropogenic. This paper reviews the observational evidence on climate phenomena and the governing equations of planetary-scale flow, as well as presenting the key concept of a hierarchy of models as used in the climate sciences. Recent advances in the application of dynamical systems theory, on the one hand, and of nonequilibrium statistical physics, on the other, are brought together for the first time and shown to complement each other in helping understand and predict the system's behavior. These complementary points of view permit a self-consistent handling of subgrid-scale phenomena as stochastic processes, as well as a unified handling of natural climate variability and forced climate change, along with a treatment of the crucial issues of climate sensitivity, response, and predictability.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A New Mathematical Framework for Atmospheric Blocking Events

We use a simple yet Earth-like atmospheric model to propose a new framework for understanding the mathematics of blocking events. Analysing error growth rates along a very long model trajectory, we show that blockings are associated with conditions of anomalously high instability of the atmosphere. Additionally, the lifetime of a blocking is positively correlated with the intensity of such an anomaly, against intuition. In the case of Atlantic blockings, predictability is especially reduced at the onset and decay of the blocking, while a relative increase of predictability is found in the mature phase, while the opposite holds for Pacific blockings, for which predictability is lowest in the mature phase. We associate blockings to a specific class of unstable periodic orbits (UPOs), natural modes of variability that cover the attractor of the system. The UPOs differ substantially in terms of instability, which explains the diversity of the atmosphere in terms predictability. The UPOs associated to blockings are indeed anomalously unstable, which leads to them being rarely visited. The onset of a blocking takes place when the trajectory of the system hops into the neighbourhood of one of these special UPOs. The decay takes place when the trajectory hops back to the neighbourhood of usual, less unstable UPOs associated with zonal flow. This justifies the classical Markov chains-based analysis of transitions between weather regimes. The existence of UPOs differing in the dimensionality of their unstable manifold indicates a very strong violation of hyperbolicity in the model, which leads to a lack of structural stability. We propose that this is could be a generic feature of atmospheric models and might be a fundamental cause behind difficulties in representing blockings for the current climate and uncertainties in predicting how their statistics will change as a result of climate change.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Response and Sensitivity Using Markov Chains

Dynamical systems are often subject to forcing or changes in their governing parameters and it is of interest to study how this affects their statistical properties. A prominent real-life example of this class of problems is the investigation of climate response to perturbations. In this respect, it is crucial to determine what the linear response of a system is to small perturbations as a quantification of sensitivity. Alongside previous work, here we use the transfer operator formalism to study the response and sensitivity of a dynamical system undergoing perturbations. By projecting the transfer operator onto a suitable finite dimensional vector space, one is able to obtain matrix representations which determine finite Markov processes. Further, using perturbation theory for Markov matrices, it is possible to determine the linear and nonlinear response of the system given a prescribed forcing. Here, we suggest a methodology which puts the scope on the evolution law of densities (the Liouville/Fokker-Planck equation), allowing to effectively calculate the sensitivity and response of two representative dynamical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modelling Complexity: the case of Climate Science

We briefly review some of the scientific challenges and epistemological issues related to climate science. We discuss the formulation and testing of theories and numerical models, which, given the presence of unavoidable uncertainties in observational data, the non-repeatability of world-experiments, and the fact that relevant processes occur in a large variety of spatial and temporal scales, require a rather different approach than in other scientific contexts. A brief discussion of the intrinsic limitations of geo-engineering solutions to global warming is presented, and a framework of investigation based upon non-equilibrium thermodynamics is proposed. We also critically discuss recently proposed perspectives of development of climate science based purely upon massive use of supercomputer and centralized planning of scientific priorities.

0 citations0 reviewsNo code linkedOpen access