Source author record

Valerio Lucarini

Valerio Lucarini 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

68works

22topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

preprint2016arXiv

Dynamical Analysis of Blocking Events: Spatial and Temporal Fluctuations of Covariant Lyapunov Vectors

One of the most relevant weather regimes in the mid-latitudes atmosphere is the persistent deviation from the approximately zonally symmetric jet to the emergence of blocking patterns. Such configurations are usually connected to exceptional local stability properties of the flow which come along with an improved local forecast skills during the phenomenon. It is instead extremely hard to predict onset and decay of blockings. Covariant Lyapunov Vectors (CLVs) offer a suitable characterization of the linear stability of a chaotic flow, since they represent the full tangent linear dynamics by a covariant basis which explores linear perturbations at all time scales. Therefore, we assess whether CLVs feature a signature of the blockings. As a first step, we examine the CLVs for a quasi-geostrophic beta-plane 2-layer model in a periodic channel baroclinically driven by a meridional temperature gradient $ΔT$. An orographic forcing enhances the emergence of localized blocked regimes. We detect the blocking events with a Tibaldi-Molteni scheme adapted to the periodic channel. When blocking occurs, the global growth rates of the fastest growing CLVs are significantly higher. Hence, against intuition, the circulation is globally more unstable in blocked phases. Such an increase in the finite time Lyapunov exponents with respect to the long term average is attributed to stronger barotropic and baroclinic conversion in the case of high temperature gradients, while for low values of ΔT, the effect is only due to stronger barotropic instability. In order to determine the localization of the CLVs we compare the meridionally averaged variance of the CLVs during blocked and unblocked phases. We find that on average the variance of the CLVs is clustered around the center of blocking. These results show that the blocked flow affects all time scales and processes described by the CLVs.

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.

preprint2016arXiv

Parametrization of stochastic multiscale triads

We discuss applications of a recently developed method for model reduction based on linear response theory of weakly coupled dynamical systems. We apply the weak coupling method to simple stochastic differential equations with slow and fast degrees of freedom. The weak coupling model reduction method results in general in a non-Markovian system, we therefore discuss the Markovianization of the system to allow for straightforward numerical integration. We compare the applied method to the equations obtained through homogenization in the limit of large time scale separation between slow and fast degrees of freedom. We look at ensemble spread from a fixed initial condition, correlation functions and exit times from domain. The weak coupling method gives better results in all test cases, albeit with a higher numerical cost.

preprint2016arXiv

Predicting Climate Change using Response Theory: Global Averages and Spatial Patterns

The provision of accurate methods for predicting the climate response to anthropogenic and natural forcings is a key contemporary scientific challenge. Using a simplified and efficient open-source general circulation model of the atmosphere featuring O($10^5$) degrees of freedom, we show how it is possible to approach such a problem using nonequilibrium statistical mechanics. Response theory allows one to practically compute the time-dependent measure supported on the pullback attractor of the climate system, whose dynamics is non-autonomous as a result of time-dependent forcings. We propose a simple yet efficient method for predicting - at any lead time and in an ensemble sense - the change in climate properties resulting from increase in the concentration of CO$_2$ using test perturbation model runs. We assess strengths and limitations of the response theory in predicting the changes in the globally averaged values of surface temperature and of the yearly total precipitation, as well as in their spatial patterns. The quality of the predictions obtained for the surface temperature fields is rather good, while in the case of precipitation a good skill is observed only for the global average. We also show how it is possible to define accurately concepts like the the inertia of the climate system or to predict when climate change is detectable given a scenario of forcing. Our analysis can be extended for dealing with more complex portfolios of forcings and can be adapted to treat, in principle, any climate observable. Our conclusion is that climate change is indeed a problem that can be effectively seen through a statistical mechanical lens, and that there is great potential for optimizing the current coordinated modelling exercises run for the preparation of the subsequent reports of the Intergovernmental Panel for Climate Change.

preprint2016arXiv

Statistical and Dynamical Properties of Covariant Lyapunov Vectors in a Coupled Atmosphere-Ocean Model - Multiscale Effects, Geometric Degeneracy, and Error Dynamics

We study a simplified coupled atmosphere-ocean model using the formalism of covariant Lyapunov vectors (CLVs), which link physically-based directions of perturbations to growth/decay rates. The model is obtained via a severe truncation of quasi-geostrophic equations for the two fluids, and includes a simple yet physically meaningful representation of their dynamical/thermodynamical coupling. The model has 36 degrees of freedom, and the parameters are chosen so that a chaotic behaviour is observed. One finds two positive Lyapunov exponents (LEs), sixteen negative LEs, and eighteen near-zero LEs. The presence of many near-zero LEs results from the vast time-scale separation between the characteristic time scales of the two fluids, and leads to nontrivial error growth properties in the tangent space spanned by the corresponding CLVs, which are geometrically very degenerate. Such CLVs correspond to two different classes of ocean/atmosphere coupled modes. The tangent space spanned by the CLVs corresponding to the positive and negative LEs has, instead, a non-pathological behaviour, and one can construct robust large deviations laws for the finite time LEs, thus providing a universal model for assessing predictability on long to ultra-long scales along such directions. It is somewhat surprising to find that the tangent space of the unstable manifold has strong projection on both atmospheric and oceanic components. Our results underline the difficulties in using hyperbolicity as a conceptual framework for multiscale chaotic dynamical systems, whereas the framework of partial hyperbolicity seems better suited, possibly indicating an alternative definition for the chaotic hypothesis. Our results suggest the need for accurate analysis of error dynamics on different time scales and domains and for a careful set-up of assimilation schemes when looking at coupled atmosphere-ocean models.

preprint2015arXiv

Covariant Lyapunov Vectors of a Quasi-geostrophic Baroclinic Model: Analysis of Instabilities and Feedbacks

The classical approach for studying atmospheric variability is based on defining a background state and studying the linear stability of the small fluctuations around such a state. Weakly non-linear theories can be constructed using higher order expansions terms. While these methods have undoubtedly great value for elucidating the relevant physical processes, they are unable to follow the dynamics of a turbulent atmosphere. We provide a first example of extension of the classical stability analysis to a non-linearly evolving quasi-geostrophic flow. The so-called covariant Lyapunov vectors (CLVs) provide a covariant basis describing the directions of exponential expansion and decay of perturbations to the non-linear trajectory of the flow. We use such a formalism to re-examine the basic barotropic and baroclinic processes of the atmosphere with a quasi-geostrophic beta-plane two-layer model in a periodic channel driven by a forced meridional temperature gradient $ΔT$. We explore three settings of $ΔT$, representative of relatively weak turbulence, well-developed turbulence, and intermediate conditions. We construct the Lorenz energy cycle for each CLV describing the energy exchanges with the background state. A positive baroclinic conversion rate is a necessary but not sufficient condition of instability. Barotropic instability is present only for few very unstable CLVs for large values of $ΔT$. Slowly growing and decaying hydrodynamic Lyapunov modes closely mirror the properties of the background flow. Following classical necessary conditions for barotropic/baroclinic instability, we find a clear relationship between the properties of the eddy fluxes of a CLV and its instability. CLVs with positive baroclinic conversion seem to form a set of modes for constructing a reduced model of the atmosphere dynamics.

preprint2015arXiv

Early 21st century snow cover state over the western river basins of the Indus River system

In this paper we assess the snow cover and its dynamics for the western river basins of the Indus River system and their subbasins located in Afghanistan, China, India and Pakistan for the period 2001-2012. First, we validate the Moderate Resolution Imaging Spectroradiometer (MODIS) daily snow products from Terra and Aqua against the Landsat Thematic Mapper/Enhanced Thematic Mapper plus data set, and then improve them for clouds by applying a validated non spectral cloud removal technique. The improved snow product has been analysed on a seasonal and annual basis against different topographic parameters (aspect, elevation and slope). Our results show a decreasing tendency for the annual average snow cover for the westerlies influenced basins (upper Indus basin, Astore, Hunza, Shigar and Shyok) and an increasing tendency for the monsoon influenced basins (Jhelum, Kabul, Swat and Gilgit). Seasonal average snow cover decreases during winter and autumn, and increases during spring and summer, which is consistent with the observed cooling and warming trends during the respective seasons. Subbasins at relatively higher latitudes/altitudes show higher variability than basins at lower latitudes/middle altitudes. Northeastern and northwestern aspects feature greater snow cover. The mean end of summer regional snow line altitude (SLA) zones range from 3000 to 5000 m asl for all basins. Our analysis provides an indication of a descending end of summer regional SLA zone for most of the studied basins, which is significant for the Shyok and Kabul basins, thus indicating a change in their water resources. Such results are consistent with the observed hydroclimatic data, recently collected local perceptions and glacier mass balances for the investigated period within the UIB. Moreover, our analysis shows a significant correlation

preprint2015arXiv

Global instability in the Ghil-Sellers model

The Ghil-Sellers energy balance model of Earth's climate, features -- for a considerable range of the solar intensity -- two stable climate states (a warm and a cold snowball Earth), where the bistability results from the celebrated ice-albedo feedback. The unstable solution is obtained and characterized in this paper. We find such unstable states by applying for the first time in a geophysical context the so-called edge tracking method that has been used for studying multiple coexisting states in shear flows. We examine robustness, efficiency, and accuracy properties of the edge tracking algorithm. We find that the procedure is the most efficient when taking a single bisection per cycle. Due to the strong diffusivity of the system trajectories of transient dynamics, initialized between the stable states with respect to the mean temperature, are confined to the heteroclininc trajectory, one which connects the fixed unstable and stable states, after relatively short transient times. This constraint dictates a functional relationship between observables. We characterize such a relationship between the global average temperature and a descriptor of nonequilibrium thermodynamics, the large scale temperature gradient between low and high latitudes. We find that a maximum of the temperature gradient is realized at the same value of the average temperature, about 270 K, largely independent of the strength of incoming solar radiation. We also examine the structural properties of the system defined by the equilibria depending on the solar strength. We construct new bifurcation diagrams in terms of quantities relevant for describing the thermodynamic disequilibrium, such as the temperature gradient and the material entropy production due to heat transport. We compare our results for the EBM to results for the intermediate complexity GCM PlaSim and find an interesting qualitative agreement.

preprint2015arXiv

Prevailing climatic trends and runoff response from Hindukush-Karakoram-Himalaya, upper Indus basin

We analyze trends in maximum, minimum and mean temperatures (Tx, Tn, and Tavg, respectively), diurnal temperature range (DTR) and precipitation from 18 stations (1250-4500 m asl) for their overlapping period of record (1995-2012), and separately, from six stations of their long term record (1961-2012). We apply Mann-Kendall test on serially independent time series to assess existence of a trend while true slope is estimated using Sen s slope method. Further, we statistically assess the spatial scale (field) significance of local climatic trends within ten identified sub-regions of UIB and analyze whether the spatially significant (field significant) climatic trends qualitatively agree with a trend in discharge out of corresponding sub-region. Over the recent period (1995-2012), we find a well agreed and mostly field significant cooling (warming) during monsoon season i.e. July-October (March-May and November), which is higher in magnitude relative to long term trends (1961-2012). We also find ... The observed hydroclimatic trends, being driven by certain changes in the monsoonal system and westerly disturbances, indicate dominance (suppression) of nival (glacial) runoff regime, altering substantially the overall hydrology of UIB in future. These findings largely contribute to address the hydroclimatic explanation of the Karakoram Anomaly.

preprint2015arXiv

Response Operators for Markov Processes in a Finite State Space: Radius of Convergence and Link to the Response Theory for Axiom A Systems

Using straightforward linear algebra we derive response operators describing the impact of small perturbations to finite state Markov processes. The results can be used for studying empirically constructed - e.g. from observations or through coarse graining of model simulations - finite state approximation of statistical mechanical systems. Recent results concerning the convergence of the statistical properties of finite state Markov approximation of the full asymptotic dynamics on the SRB measure in the limit of finer and finer partitions of the phase space are suggestive of some degree of robustness of the obtained results in the case of Axiom A system. Our findings give closed formulas for the linear and nonlinear response theory at all orders of perturbation and provide matrix expressions that can be directly implemented in any coding language, plus providing bounds on the radius of convergence of the perturbative theory. In particular, we relate the convergence of the response theory to the rate of mixing of the unperturbed system. One can use the formulas obtained for finite state Markov processes to recover previous findings obtained on the response of continuous time Axiom A dynamical systems to perturbations, by considering the generator of time evolution for the measure and for the observables. A very basic, low-tech, and computationally cheap analysis of the response of the Lorenz '63 model to perturbations provides rather encouraging results regarding the possibility of using the approximate representation given by finite state Markov processes to compute the system's response.

preprint2015arXiv

Seasonal cycle of Precipitation over Major River Basins in South and Southeast Asia: A Review of the CMIP5 climate models data for present climate and future climate projections

We review the skill of thirty coupled climate models participating in Coupled Model Intercomparison Project 5 in terms of reproducing properties of the seasonal cycle of precipitation over the major river basins of South and Southeast Asia (Indus, Ganges, Brahmaputra and Mekong) for historical period (1961-2000). We also present projected changes by these models by end of century (2061-2100) under extreme scenario RCP8.5. First, we assess their ability to reproduce observed timings of the monsoon onset and the rate of rapid fractional accumulation (RFA slope) - a measure of seasonality within active monsoon period. Secondly, we apply a threshold-independent seasonality index (SI) - a multiplicative measure of precipitation and extent of its concentration relative to the uniform distribution (relative entropy - RE). We apply SI distinctly for monsoonal precipitation regime (MPR), westerly precipitation regime (WPR) and annual precipitation regime. For present climate, neither any single model nor the multi-model mean performs best in all chosen metrics. Models show overall a modest skill in suggesting right timings of the monsoon onset while the RFA slope is generally underestimated. One third of models fail to capture the monsoon signal over the Indus basin. Mostly, SI estimates for WPR are simulated higher than observed for all basins, while for MPR, it is simulated higher (lower) for Ganges and Brahmaputra (Indus and Mekong) basins, following the pattern of overestimation (underestimation) of precipitation. However, models are biased positive (negative) for RE estimates over Indus and Ganges (Brahmaputra and Mekong) basins, implying the extent of precipitation concentration for MPR and number of dry days within WPR higher (lower) than observed for these basins. Under the RCP8.5 scenario, most of models project a slightly delayed monsoon onset, and a general increase in the RFA slope...

preprint2015arXiv

Seasonality of the hydrological cycle in major South and Southeast Asian river basins as simulated by PCMDI/CMIP3 experiments

In this study, we investigate how PCMDI/CMIP3 general circulation models (GCMs) represent the seasonal properties of the hydrological cycle in four major South and Southeast Asian river basins (Indus, Ganges, Brahmaputra and Mekong). First, we examine the skill of the GCMs by analysing their performance in simulating the 20th century climate (1961-2000 period) using historical forcing (20c3m experiment), and then we analyse the projected changes for the corresponding 21st and 22nd century climates under the SRESA1B scenario. The CMIP3 GCMs show a varying degree of skill in simulating the basic characteristics of the monsoonal precipitation regimes of the Ganges, Brahmaputra and Mekong basins, while the representation of the hydrological cycle over the Indus Basin is poor in most cases, with a few GCMs not capturing the monsoonal signal at all. While the model outputs feature a remarkable spread for the monsoonal precipitation, a satisfactory representation of the western mid-latitude precipitation regime is instead observed. Similarly, most of the models exhibit a satisfactory agreement for the basin-integrated runoff in winter and spring, while their spread is large for the runoff during the monsoon season. For the future climate scenarios, most models foresee a decrease in the winter P - E over all four basins, while agreement is found on the decrease of the spring P - E over the Indus and Ganges basins only. Such decreases in P - E are mainly due to the decrease in precipitation associated with the western mid-latitude disturbances. Consequently, for the Indus and Ganges basins, the runoff drops during the spring season while it rises during the winter season. Such changes indicate a shift from rather glacial and nival to more pluvial runoff regimes, particularly for the Indus Basin. Furthermore, the rise in the projected runoff, along with ...

preprint2014arXiv

A new framework for climate sensitivity and prediction: a modelling perspective

The sensitivity of climate models to increasing CO2 concentration and the climate response at decadal time scales are still major factors of uncertainty for the assessment of the long and short term effects of anthropogenic climate change. While the relative slow progress on these issues is partly due to the inherent inaccuracies of numerical climate models, this also hints at the need for stronger theoretical foundations to the problem of studying climate sensitivity and performing climate change predictions with numerical models. Here we demonstrate that it is possible to use Ruelle's response theory to predict the impact of an arbitrary CO2 forcing scenario on the global surface temperature of a general circulation model. Response theory puts the concept of climate sensitivity on firm theoretical grounds, and addresses rigorously the problem of predictability at different time scales. Conceptually, our results show that performing climate change experiments with general circulation models is a well defined problem from a physical and mathematical point of view. Practically, our results show that considering one single CO2 forcing scenario is enough to construct operators able to predict the response of climatic observables to any other CO2 forcing scenario, without the need to perform additional numerical simulations. We also introduce a general relationship between climate sensitivity and climate response at different time scales, thus providing an explicit definition of the inertia of the system at different time scales. While what we report here refers to the linear response, the general theory allows for treating nonlinear effects as well. Our results pave the way for redesigning and interpreting climate change experiments from a radically new perspective.

preprint2014arXiv

Analysis of rainfall seasonality from observations and climate models

Two new indicators of rainfall seasonality based on information entropy, the relative entropy (RE) and the dimensionless seasonality index (DSI), together with the mean annual rainfall, are evaluated on a global scale for recently updated precipitation gridded datasets and for historical simulations from coupled atmosphere-ocean general circulation models. The RE provides a measure of the number of wet months and, for precipitation regimes featuring one maximum in the monthly rain distribution, it is related to the duration of the wet season. The DSI combines the rainfall intensity with its degree of seasonality and it is an indicator of the extent of the global monsoon region. We show that the RE and the DSI are fairly independent of the time resolution of the precipitation data, thereby allowing objective metrics for model intercomparison and ranking. Regions with different precipitation regimes are classified and characterized in terms of RE and DSI. Comparison of different land observational datasets reveals substantial difference in their local representation of seasonality. It is shown that two-dimensional maps of RE provide an easy way to compare rainfall seasonality from various datasets and to determine areas of interest. CMIP5 models consistently overestimate the RE over tropical Latin America and underestimate it in Western Africa and East Asia. It is demonstrated that positive RE biases in a GCM are associated with simulated monthly precipitation fractions which are too large during the wet months and too small in the months preceding the wet season; negative biases are instead due to an excess of rainfall during the dry months.

preprint2014arXiv

Climate of Earth-like planets with high obliquity and eccentric orbits: implications for habitability conditions

We explore the effects of seasonal variability for the climate of Earth-like planets as determined by the two parameters polar obliquity and orbital eccentricity using a general circulation model of intermediate complexity. In the first part of the paper we examine the consequences of different values of obliquity and eccentricity for the spatio-temporal patterns of radiation and surface temperatures as well as for the main characteristics of the atmospheric circulations. In the second part we analyse the associated implications for the habitability of planets close to the outer edge of the habitable zone (HZ). This part of the paper focuses in particular on the multistability property of climate, i.e. the parallel existence of both an ice-free and an ice-covered climate state. Our results show that seasonal variability affects both the existence of and transitions between the two climate states. Moreover, our experiments reveal that planets with Earth-like atmospheres and high seasonal variability can have ice-free areas at much larger distance from the host star than planets without seasonal variability, which leads to a substantial expansion of the outer edge of the HZ. Sensitivity experiments exploring the role of azimuthal obliquity and surface heat capacity test the robustness of our results. On circular orbits, our findings obtained with a general circulation model agree well with previous studies based on one dimensional energy balance models, whereas significant differences are found on eccentric orbits.

preprint2014arXiv

Equivalence of Non-Equilibrium Ensembles and Representation of Friction in Turbulent Flows: The Lorenz 96 Model

We construct different equivalent non-equilibrium statistical ensembles in a simple yet instructive $N$-degrees of freedom model of atmospheric turbulence, introduced by Lorenz in 1996. The vector field can be decomposed into an energy-conserving, time-reversible part, plus a non-time reversible part, including forcing and dissipation. We construct a modified version of the model where viscosity varies with time, in such a way that energy is conserved, and the resulting dynamics is fully time-reversible. For each value of the forcing, the statistical properties of the irreversible and reversible model are in excellent agreement, if in the latter the energy is kept constant at a value equal to the time-average realized with the irreversible model. In particular, the average contraction rate of the phase space of the time-reversible model agrees with that of the irreversible model, where instead it is constant by construction. We also show that the phase space contraction rate obeys the fluctuation relation, and we relate its finite time corrections to the characteristic time scales of the system. A local version of the fluctuation relation is explored and successfully checked. The equivalence between the two non-equilibrium ensembles extends to dynamical properties such as the Lyapunov exponents, which are shown to obey to a good degree of approximation a pairing rule. These results have relevance in motivating the importance of the chaotic hypothesis. in explaining that we have the freedom to model non-equilibrium systems using different but equivalent approaches, and, in particular, that using a model of a fluid where viscosity is kept constant is just one option, and not necessarily the only option, for describing accurately its statistical and dynamical properties.

preprint2014arXiv

Mathematical and Physical Ideas for Climate Science

The climate is a forced and dissipative nonlinear system featuring non-trivial dynamics of a vast range of spatial and temporal scales. The understanding of the climate's structural and multiscale properties is crucial for the provision of a unifying picture of its dynamics and for the implementation of accurate and efficient numerical models. We present some recent developments at the intersection between climate science, mathematics, and physics, which may prove fruitful in the direction of constructing a more comprehensive account of climate dynamics. We describe the Nambu formulation of fluid dynamics, and the potential of such a theory for constructing sophisticated numerical models of geophysical fluids. Then, we focus on the statistical mechanics of quasi-equilibrium flows in a rotating environment, which seems crucial for constructing a robust theory of geophysical turbulence. We then discuss ideas and methods suited for approaching directly the non-equilibrium nature of the climate system. First, we describe some recent findings on the thermodynamics of climate and characterize its energy and entropy budgets, and discuss related methods for intercomparing climate models and for studying tipping points. These ideas can also create a common ground between geophysics and astrophysics by suggesting general tools for studying exoplanetary atmospheres. We conclude by focusing on non-equilibrium statistical mechanics, which allows for a unified framing of problems as different as the climate response to forcings, the effect of altering the boundary conditions or the coupling between geophysical flows, and the derivation of parametrizations for numerical models.

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

Projected changes of rainfall seasonality and dry spells in a high concentration pathway 21st century scenario

In this diagnostic study we analyze changes of rainfall seasonality and dry spells by the end of the twenty-first century under the most extreme IPCC5 emission scenario (RCP8.5) as projected by twenty-four coupled climate models participating to Coupled Model Intercomparison Project 5. We use estimates of the centroid of the monthly rainfall distribution as an index of the rainfall timing and a threshold-independent, information theory-based quantity such as relative entropy (RE) to quantify the concentration of annual rainfall and the number of dry months and to build a monsoon dimensionless seasonality index (DSI). The RE is projected to increase, with high inter-model agreement over Mediterranean-type regions (southern Europe, northern Africa and southern Australia) and areas of South and Central America, implying an increase in the number of dry days up to one month by the end of the twenty-first century. Positive RE changes are also projected over the monsoon regions of southern Africa and North America, South America. These trends are consistent with a shortening of the wet season associated with a more prolonged pre-monsoonal dry period. The extent of the global monsoon region, characterized by large DSI, is projected to remain substantially unaltered. Centroid analysis shows that most of CMIP5 projections suggest that the monsoonal annual rainfall distribution is expected to change from early to late in the course of the hydrological year by the end of the twenty-first century and particularly after year 2050. This trend is particularly evident over Northern Africa, Southern Africa and western Mexico, where more than 90 % of the models project a delay of the rainfall centroid from a few days up to two weeks. Over the remaining monsoonal regions, there is little inter-model agreement in terms of centroid changes.

preprint2014arXiv

The Impact of Oceanic Heat Transport on the Atmospheric Circulation: a Thermodynamic Perspective

The present study investigates how global thermodynamic properties of the climate system are affected by the changes in the intensity of the imposed oceanic heat transport in an atmospheric general circulation model in aqua-planet configuration. Increasing the poleward oceanic heat transport results in an overall increase in the surface temperature and a decrease in the equator-to-pole surface temperature difference as a result of the ice-albedo feedback. Following the classical ansatz by Stone, the atmospheric heat transport changes in such a way that the total poleward heat transport remains almost unchanged. We also find that the efficiency of the climate machine, the intensity of the Lorenz energy cycle and the material entropy production of the system decline with increased oceanic heat transport which suggests that the climate system becomes less efficient and turns into a state of reduced entropy production, as the enhanced oceanic transport performs a stronger large-scale mixing between geophysical fluids with different temperature, thus reducing the availability in the climate system and bringing it closer to a state of thermal equilibrium.

preprint2013arXiv

Avalanches, Breathers and Flow Reversal in a Continuous Lorenz-96 Model

For the discrete model suggested by Lorenz in 1996 a one-dimensional long wave approximation with nonlinear excitation and diffusion is derived. The model is energy conserving but non-Hamiltonian. In a low order truncation weak external forcing of the zonal mean flow induces avalanche-like breather solutions which cause reversal of the mean flow by a wave-mean flow interaction. The mechanism is an outburst-recharge process similar to avalanches in a sand pile model.

preprint2013arXiv

Elements of a unified framework for response formulae

We provide a physical interpretation of the first and second order terms occurring in Ruelle's response formalism. We show that entropy fluxes play a major role in determining the response of the system to perturbations. Along this line, we show that our framework allows one to recover a wealth of previous results of response theory in both deterministic and stochastic contexts. In particular, we are able to shed light on the crosslinks between the dynamical systems approach á la Ruelle and large deviations methods.

preprint2013arXiv

Entropy Production and Coarse Graining of the Climate Fields in a General Circulation Model

We extend the analysis of the thermodynamics of the climate system by investigating the role played by processes taking place at various spatial and temporal scales through a procedure of coarse graining. The coarser is the graining of the climatic fields, the lower is the resulting estimate of the material entropy production. In other terms, all the spatial and temporal scales of variability of the thermodynamic fields provide a positive contribution to the material entropy production. At all scales, the temperature fields and the heating fields resulting from the convergence of turbulent fluxes have a negative correlation, while the opposite holds between the temperature fields and the radiative heating fields. Moreover, the latter correlations are stronger, which confirms that radiation acts as primary driver for the climatic processes, while the material fluxes dampen the resulting fluctuations through dissipative processes. Using specific coarse-graining procedures, one can separate the various contributions to the material entropy production coming from the dissipation of kinetic energy, the vertical sensible and latent heat fluxes, and the large scale horizontal fluxes. Most of the entropy production is associated to irreversible exchanges occurring along the vertical direction, and neglecting the horizontal and time variability of the fields has a small impact on the estimate of the material entropy production. The approach presented here seems promising for testing climate models, for assessing the impact of changing their parametrizations and their resolution, and for investigating the atmosphere of exoplanets. Our findings deserve to be explored in a more general context, because they provide a way for understanding the relationship between forced fluctuations and dissipative processes in continuum systems.

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

Habitability and Multistability in Earth-like Planets

We explore the potential multistability of the climate for a planet around the habitable zone. We focus on conditions reminiscent to those of the Earth system, but our investigation aims at presenting a general methodology for dealing with exoplanets. We provide a thorough analysis of the non-equilibrium thermodynamical properties of the climate system and explore, using a a flexible climate model, how such properties depend on the energy input of the parent star, on the infrared atmospheric opacity, and on the rotation rate. It is possible to reproduce the multi-stability properties reminiscent of the paleoclimatologically relevant snowball (SB) - warm (W) conditions. We then study the thermodynamics of the W and SB states, clarifying the role of the hydrological cycle in shaping the irreversibility and the efficiency of the W states, and emphasizing the extreme diversity of the SB states, where dry conditions are realized. Thermodynamics provides the clue for studying the tipping points of the system and leads us to constructing parametrizations where the main thermodynamic properties are expressed as functions of the emission temperature of the planet only. Such functions are rather robust with respect to changing the rotation rate of the planet from the current terrestrial one to half of it. We then explore the dynamical range of slowy rotating and phase locked planets. There is a critical rotation rate below which the multi-stability properties are lost. Such critical rotation rate corresponds roughly to the phase lock 2:1 condition. Therefore, if an Earth-like planet is 1:1 phase locked with respect to the parent star, only one climatic state would be compatible with a given set of astronomical and astrophysical parameters. These results have relevance for the general theory of planetary circulation and for the definition of necessary and sufficient conditions for habitability.

preprint2013arXiv

Hydrological Cycle over South and Southeast Asian River Basins as Simulated by PCMDI/CMIP3 Experiments

We investigate how the climate models contributing to the PCMDI/CMIP3 dataset describe the hydrological cycle over four major South and Southeast Asian river basins (Indus, Ganges, Brahmaputra and Mekong) for the 20th, 21st (13 models) and 22nd (10 models) centuries. For the 20th century, some models do not seem to conserve water at the river basin scale up to a good degree of approximation. The simulated precipitation minus evaporation (P - E), total runoff (R) and precipitation (P) quantities are neither consistent with the observations nor among the models themselves. Most of the models underestimate P - E for all four river basins, which is mainly associated with the underestimation of precipitation. This is in agreement with the recent results on the biases of the representation of monsoonal dynamics by GCMs. Overall, a modest inter-model agreement is found only for the evaporation and inter-annual variability of P - E. For the 21st and 22nd centuries, models agree on the negative (positive) changes of P - E for the Indus basin (Ganges, Brahmaputra and Mekong basins). Most of the models foresee an increase in the inter-annual variability of P - E for the Ganges and Mekong basins, thus suggesting an increase in large low-frequency dry/wet events. Instead, no considerable future change in the inter-annual variability of P - E is found for the Indus and Brahmaputra basins.

preprint2013arXiv

Nonequilibrium thermodynamics of circulation regimes in optically-thin, dry atmospheres

An extensive analysis of an optically-thin, dry atmosphere at different values of the thermal Rossby number Ro and of the Taylor number Ff is per- formed with a general circulation model by varying the rotation rate Ω and the surface drag τ in a wide parametric range. By using nonequilibrium thermodynamics diagnostics such as material entropy production, efficiency, meridional heat transport and kinetic energy dissipation we characterize in a new way the different circulation regimes. Baroclinic circulations feature high mechanical dissipation, meridional heat transport, material entropy pro- duction and are fairly efficient in converting heat into mechanical work. The thermal dissipation associated with the sensible heat flux is found to depend mainly on the surface properties, almost independent from the rotation rate and very low for quasi-barotropic circulations and regimes approaching equa- torial super-rotation. Slowly rotating, axisymmetric circulations have the highest meridional heat transport. At high rotation rates and intermediate- high drag, atmospheric circulations are zonostrohic with very low mechanical dissipation, meridional heat transport and efficiency. When τ is interpreted as a tunable parameter associated with the turbulent boundary layer trans- fer of momentum and sensible heat, our results confirm the possibility of using the Maximum Entropy Production Principle as a tuning guideline in the range of values of Ω. This study suggests the effectiveness of using fun- damental nonequilibrium thermodynamics for investigating the properties of planetary atmospheres and extends our knowledge of the thermodynamics of the atmospheric circulation regimes.

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

Beyond the linear Fluctuation-Dissipation Theorem: the Role of Causality

In this paper we re-examine the traditional problem of connecting the internal fluctuations of a system to its response to external forcings and extend the classical theory in order to be able to encompass also nonlinear processes. With this goal, we try to join on the results by Kubo on statistical mechanical systems close to equilibrium, i.e. whose unperturbed state can be described by a canonical ensemble, the theory of dispersion relations, and the response theory recently developed by Ruelle for non-equilibrium systems equipped with an invariant SRB measure. Our derivations highlight the strong link between causality and the possibility of connecting unambiguously fluctuation and response, both at linear and nonlinear level. We first show in a rather general setting how the formalism of the Ruelle response theory can be used to derive in a novel way Kramers-Kronig relations connecting the real and imaginary part of the linear and nonlinear response to external perturbations. We then provide a formal extension at each order of nonlinearity of the fluctuation-dissipation theorem (FDT) for general systems possessing a smooth invariant measure. Finally, we focus on the physically relevant case of systems close to equilibrium, for which we present explicit fluctuation-dissipation relations linking the susceptibility describing the $n^{th}$ order response of the system with the expectation value of suitably defined correlations of $n+1$ observables taken in the equilibrium ensemble. While the FDT has an especially compact structure in the linear case, in the nonlinear case joining the statistical properties of the fluctuations of the system to its response to external perturbations requires linear changes of variables, simple algebraic sums and multiplications, and a multiple convolution integral. These operations, albeit cumbersome, can be easily implemented numerically.

preprint2012arXiv

Bistability of the climate around the habitable zone: a thermodynamic investigation

The goal of this paper is to explore the potential multistability of the climate of a planet around the habitable zone. A thorough investigation of the thermodynamics of the climate system is performed for very diverse conditions of energy input and infrared atmosphere opacity. Using PlaSim, an Earth-like general circulation model, the solar constant S* is modulated between 1160 and 1510 Wm-2 and the CO2 concentration, [CO2], from 90 to 2880 ppm. It is observed that in such a parameter range the climate is bistable, i.e. there are two coexisting attractors, one characterised by warm, moist climates (W) and one by completely frozen sea surface (Snowball Earth, SB). Linear relationships are found for the two transition lines (W\rightarrowSB and SB\rightarrowW) in (S*,[CO2]) between S* and the logarithm of [CO2]. The dynamical and thermodynamical properties - energy fluxes, Lorenz energy cycle, Carnot efficiency, material entropy production - of the W and SB states are very different: W states are dominated by the hydrological cycle and latent heat is prominent in the material entropy production; the SB states are predominantly dry climates where heat transport is realized through sensible heat fluxes and entropy mostly generated by dissipation of kinetic energy. We also show that the Carnot efficiency regularly increases towards each transition between W and SB, with a large decrease in each transition. Finally, we propose well-defined empirical functions allowing for expressing the global non-equilibrium thermodynamical properties of the system in terms of either the mean surface temperature or the mean planetary emission temperature. This paves the way for the possibility of proposing efficient parametrisations of complex non-equilibrium properties and of practically deducing fundamental properties of a planetary system from a relatively simple observable.

preprint2012arXiv

Multi-level Dynamical Systems: Connecting the Ruelle Response Theory and the Mori-Zwanzig Approach

In this paper we consider the problem of deriving approximate autonomous dynamics for a number of variables of a dynamical system, which are weakly coupled to the remaining variables. In a previous paper we have used the Ruelle response theory on such a weakly coupled system to construct a surrogate dynamics, such that the expectation value of any observable agrees, up to second order in the coupling strength, to its expectation evaluated on the full dynamics. We show here that such surrogate dynamics agree up to second order to an expansion of the Mori-Zwanzig projected dynamics. This implies that the parametrizations of unresolved processes suited for prediction and for the representation of long term statistical properties are closely related, if one takes into account, in addition to the widely adopted stochastic forcing, the often neglected memory effects.

preprint2012arXiv

Nambu representation of an extended Lorenz model with viscous heating

We consider the Nambu and Hamiltonian representations of Rayleigh-Benard convection with a nonlinear thermal heating effect proportional to the Eckert number (Ec). The model we use is an extension of the classical Lorenz-63 model with 4 kinematic and 6 thermal degrees of freedom. The conservative parts of the dynamical equations which include all nonlinearities satisfy Liouville's theorem and permit a conserved Hamiltonian H for arbitrary Ec. For Ec=0 two independent conserved Casimir functions exist, one of these is associated with unavailable potential energy and is also present in the Lorenz-63 truncation. This Casimir C is used to construct a Nambu representation of the conserved part of the dynamical system. The thermal heating effect can be represented either by a second canonical Hamiltonian or as a gradient (metric) system using the time derivative of the Casimir. The results demonstrate the impact of viscous heating in the total energy budget and in the Lorenz energy cycle for kinetic and available potential energy.

preprint2011arXiv

Baroclinic stationary waves in aquaplanet models

An aquaplanet model is used to study the nature of the highly persistent low frequency waves that have been observed in models forced by zonally symmetric boundary conditions. Using the Hayashi spectral analysis of the extratropical waves, we find that a quasi-stationary (QS) wave five belongs to a wave packet obeying a well defined dispersion relation with eastward group velocity. The components of the dispersion relation with k>5 baroclinically convert eddy available potential energy into eddy kinetic energy, while those with k<5 are baroclinically neutral. In agreement with the Green's model of baroclinic instability, the wave five is weakly unstable, and the inverse energy cascade, which had been previously proposed as a main forcing for this type of waves, only acts as a positive feedback on its predominantly baroclinic energetics. The QS wave is reinforced by a phase lock to an analogous pattern in the tropical convection, which provides further amplification to the wave. We also find that the Pedlosky bounds on the phase speed of unstable waves provide guidance in explaining the latitudinal structure of the energy conversion, which is shown to be more enhanced where the zonal westerly surface wind is weaker. The wave energy is then trapped in the wave guide created by the upper tropospheric jet stream. In agreement with Green's theory, as the equator to pole SST difference is reduced the stationary marginally stable component shifts toward higher wavenumbers, while the wave five becomes neutral and westward propagating. Some properties of the aquaplanet QS waves are found in interesting agreement with a low frequency wave observed by Salby (1982) in the southern hemisphere DJF, so that this perspective on low frequency variability might be, apart from its value in terms of basic geophysical fluid dynamics, of specific interest for studying the Earth's atmosphere.

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

Climatology of extratropical atmospheric wave packets in the northern hemisphere

Planetary and synoptic scale wave-packets represent one important component of the atmospheric large-scale circulation. These dissipative structures are able to rapidly transport eddy kinetic energy, generated locally (e.g. by baroclinic conversion), downstream along the upper tropospheric flow. The transported energy, moving faster than individual weather systems, will affect the development of the next meteorological system on the leading edge of the wave packet, creating a chain of connections between systems that can be far apart in time and space, with important implications on predictability. In this work we present an automated recognition of atmospheric wave packets which allows the extraction of all the relevant properties, such as location, duration and velocity. Behind this tool lies the need to investigate atmospheric variability in its full complexity, bridging the low-frequency steady-state approach with the storm-tracks lagrangian approach. We have applied the algorithm to the daily analysis (every 12h) from 1958-2010, building an extended climatology of waves packets with different spectral properties. We show that wave packets characteristics over Northern Hemisphere exhibit a strong seasonal dependence, both in their spectral component and in their distribution and localization. The maximum activity is reached in the cold months, from autumn to spring, with a slight weakening in mid-winter and a clear minimum of activity in summer. Preferential areas of genesis are shown to be the Western and Central-Pacific and Western-Atlantic while areas of lyses are the eastern borders of Pacific and Atlantic. We envisage possible applications of this algorithm also for predictability studies and operational activities.

preprint2011arXiv

Disentangling multi-level systems: averaging, correlations and memory

We consider two weakly coupled systems and adopt a perturbative approach based on the Ruelle response theory to study their interaction. We propose a systematic way to parametrize the effect of the coupling as a function of only the variables of a system of interest. Our focus is on describing the impacts of the coupling on the long-term statistics rather than on the finite-time behaviour. By direct calculation, we find that, at first order, the coupling can be surrogated by adding a deterministic perturbation to the autonomous dynamics of the system of interest. At second order, there are additionally two separate and very different contributions. One is a term taking into account the second order contributions of the fluctuations in the coupling, which can be parametrized as a stochastic forcing with given spectral properties. The other one is a memory term, coupling the system of interest to its previous history, through the correlations of the second system. If these correlations are known, this effect can be implemented as a perturbation with memory on the single system. In order to treat this case, we present an extension to Ruelle's response theory able to deal with integral operators. We discuss our results in the context of other methods previously proposed to disentangle the dynamics of two coupled systems. We emphasize that our results do not rely on assuming a time scale separation, and, if such a separation exist, can be used equally well to study the statistics of the slow as well as that of the fast variables. By recursively applying the technique proposed here, we can treat the general case of multilevel systems.

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

Hydrological Cycle in the Danube basin in present-day and XXII century simulations by IPCCAR4 global climate models

We present an intercomparison and verification analysis of 20 GCMs included in the 4th IPCC assessment report regarding their representation of the hydrological cycle on the Danube river basin for 1961-2000 and for the 2161-2200 SRESA1B scenario runs. The basin-scale properties of the hydrological cycle are computed by spatially integrating the precipitation, evaporation, and runoff fields using the Voronoi-Thiessen tessellation formalism. The span of the model simulated mean annual water balances is of the same order of magnitude of the observed Danube discharge of the Delta; the true value is within the range simulated by the models. Some land components seem to have deficiencies since there are cases of violation of water conservation when annual means are considered. The overall performance and the degree of agreement of the GCMs are comparable to those of the RCMs analyzed in a previous work, in spite of the much higher resolution and common nesting of the RCMs. The reanalyses are shown to feature several inconsistencies and cannot be used as a verification benchmark for the hydrological cycle in the Danubian region. In the scenario runs, for basically all models the water balance decreases, whereas its interannual variability increases. Changes in the strength of the hydrological cycle are not consistent among models: it is confirmed that capturing the impact of climate change on the hydrological cycle is not an easy task over land areas. Moreover, in several cases we find that qualitatively different behaviours emerge among the models: the ensemble mean does not represent any sort of average model, and often it falls between the models' clusters.

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.

preprint2011arXiv

New Results on the Thermodynamical Properties of the Climate System

In this paper we exploit two equivalent formulations of the average rate of material entropy production in a planetary system to propose an approximate splitting between contributions due to vertical and eminently horizontal processes. Our approach is based only upon 2D radiative fields at the surface and at the top of atmosphere of a general planetary body. Using 2D fields at the top of atmosphere alone, we derive lower bounds to the rate of material entropy production and to the intensity of the Lorenz energy cycle. By introducing a measure of the efficiency of the planetary system with respect to horizontal thermodynamical processes, we provide insight on a previous intuition on the possibility of defining a baroclinic heat engine extracting work from the meridional heat flux. The approximate formula of the material entropy production is verified and used for studying the global thermodynamic properties of climate models (CMs) included in the PCMDI/CMIP3 dataset in pre-industrial climate conditions. It is found that about 90% of the material entropy production is due to vertical processes such as convection, whereas the large scale meridional heat transport contributes only about 10%. The total material entropy production is typically 55 mWK-1m-2, with discrepancies of the order of 5% and CMs' baroclinic efficiencies are clustered around 0.055. When looking at the variability and co-variability of the considered thermodynamical quantities, the agreement among CMs is worse, suggesting that the description of feedbacks is more uncertain.

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

Stochastic perturbations to dynamical systems: a response theory approach

We study the impact of stochastic perturbations to deterministic dynamical systems using the formalism of the Ruelle response theory and explore how stochastic noise can be used to explore the properties of the underlying deterministic dynamics of a system. We find the expression for the change in the expectation value of a general observable when a white noise forcing is introduced in the system, both in the case of additive and multiplicative noise. We also show that the difference between the expectation value of the power spectrum of an observable in the stochastically perturbed case and of the same observable in the unperturbed case is equal to the variance of the noise times the square of the modulus of the susceptibility describing the frequency-dependent response of the system to perturbations with the same spatial patterns as the considered stochastic forcing. Using Kramers-Kronig theory, it is then possible to derive the susceptibility and thus deduce the Green function of the system for any desired observable. We then extend our results to rather general patterns of random forcing, from the case of several white noise forcings, to noise terms with memory, up to the case of a space-time random field. Explicit formulas are provided for each relevant case. As a general result, we find, using an argument of positive-definiteness, that the power spectrum of the stochastically perturbed system is larger at all frequencies than the power spectrum of the unperturbed system. We provide a example of application of our results by considering the Lorenz 96 model. These results clarify the property of stochastic stability of SRB measures in Axiom A flows, provide tools for analysing stochastic parameterisations and related closure ansatz to be implemented in modelling studies, and introduce new ways to study the response of a system to external perturbations.

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.

preprint2010arXiv

A Statistical Mechanical Approach for the Computation of the Climatic Response to General Forcings

The climate belongs to the class of non-equilibrium forced and dissipative systems, for which most results of quasi-equilibrium statistical mechanics, including the fluctuation-dissipation theorem, do not apply. We show for the first time how the Ruelle linear response theory, developed for studying rigorously the impact of perturbations on general observables of non-equilibrium statistical mechanical systems, can be applied to analyze the climatic response. We choose as test bed the Lorenz 96 model, which has a well-recognized prototypical value. We recapitulate the main aspects of the response theory and propose some new results. We then analyze the frequency dependence of the response of both local and global observables to perturbations with localized as well as global spatial patterns. We derive analytically the asymptotic behaviour, validity of Kramers-Kronig relations, and sum rules for the susceptibilities, and related them to parameters describing the unperturbed properties of the system. We verify the theoretical predictions from the outputs of the simulations with great precision. The theory is used to explain differences in the response of local and global observables, in defining the intensive properties of the system and in generalizing the concept of climate sensitivity to all time scales. We also show how to reconstruct the linear Green function, which maps perturbations of general time patterns into changes in the expectation value of the considered observable. Finally, we propose a general methodology to study Climate Change problems by resorting to few, well selected simulations and discuss the specific case of surface temperature response to changes of the $CO_2$ concentration. This approach may provide a radically new perspective to study rigorously the problem of climate sensitivity and climate change.

preprint2010arXiv

Energetics of PCMDI/CMIP3 Climate Models: Energy Budget and Meridional Enthalpy Transport

We analyze the PCMDI/CMIP3 simulations performed by climate models (CMs) using pre-industrial and SRESA1B scenarios. Relatively large biases are present for most CMs when global energy budgets and when the atmospheric, oceanic, and land budgets are considered. Apparently, the biases do not result from transient effects, but depend on the imperfect closure of the energy cycle in the fluid components and on inconsistencies over land. Therefore, the planetary emission temperature is underestimated. This may explain the CMs' cold bias. In the pre-industrial scenario, CMs agree on the location in the mid-latitudes of the peaks of the meridional atmospheric enthalpy transport, while large discrepancies exist on the intensity. Disagreements on the location and intensity of the oceanic transport peaks are serious. With increased $CO_2$ concentration, a small poleward shift of the peak and an increase in the intensity of the atmospheric transport of up to 10% are detected in both hemispheres. Instead, most CMs feature a decrease in the oceanic transport intensity in the northern hemisphere and an equatorward shift of the peak in both hemispheres. The Bjerkens compensation mechanism is active both on climatological and interannual time scales. The peak of the total meridional transport is typically around $35^\circ$ in both hemispheres and scenarios, whereas disagreements on the intensity are relevant. With increased $CO_2$ concentration, the total transport increases by up to 10%, thus contributing to polar amplification. Advances in the representation of physical processes are definitely needed for providing a self-consistent representation of climate as a non-equilibrium thermodynamical system.

preprint2009arXiv

Thermodynamic Analysis of Snowball Earth Hysteresis Experiment: Efficiency, Entropy Production, and Irreversibility

We present an extensive thermodynamic analysis of a hysteresis experiment performed on a simplified yet Earth-like climate model. We slowly vary the solar constant by 20% around the present value and detect a substantial bistability: for a large range of values the realization of snowball (SB) or of warm (W) climate conditions depend on the history of the system. Using recent results on the global climate thermodynamics, we show that that the two regimes feature radically different properties. The efficiency of the climate machine increases with decreasing solar constant in W climate conditions, whereas the opposite takes place in the SB regime. Instead, entropy production is increasing with the solar constant in both branches of climate conditions, and its value is about 4 times as large in the W branch than in the corresponding C state. Finally, the degree of irreversibility of the system is much higher in the W conditions, with an explosive growth in the upper range of the considered values of solar constants. Whereas in the SB regime a dominating role is played by changes in the meridional albedo contrast, in the W climate regime changes in the intensity of latent heat fluxes are crucial for determining the observed properties. This substantiates the importance of addressing correctly the variations of the hydrological cycle in a changing climate. An interpretation of the transitions at the boundary of the bistable region based upon macro-scale thermodynamic properties is also proposed. Our results support the adoption of a new generation of diagnostic tools based on the 2nd law of thermodynamics for auditing climate model and outline a set of parameterizations to be used in conceptual and intermediate complexity models or for the reconstruction of the past climate conditions.

preprint2009arXiv

Thermodynamic Efficiency and Entropy Production in the Climate System

We present a new outlook on the climate system thermodynamics, studying some of its macroscopic properties in terms of the 1st and 2nd laws of thermodynamics. We review and clarify the notion of efficiency of the climate system by constructing formally an equivalent Carnot engine with efficiency eta, and show how the Lorenz energy cycle can be framed in a macro-scale thermodynamic context. Then, by exploiting the 2nd law, we prove that the lower bound to the entropy production is eta times the integrated absolute value of the internal entropy fluctuations. An exergetic interpretation is also proposed. Finally, the controversial maximum entropy production principle is re-interpreted as requiring the joint optimization of heat transport and mechanical work production. These results provide new tools for climate change analysis and for climate models' validation.

preprint2008arXiv

Evidence of dispersion relations for the nonlinear response of the Lorenz 63 system

Along the lines of the nonlinear response theory developed by Ruelle, in a previous paper we have proved under rather general conditions that Kramers-Kronig dispersion relations and sum rules apply for a class of susceptibilities describing at any order of perturbation the response of Axiom A non equilibrium steady state systems to weak monochromatic forcings. We present here the first evidence of the validity of these integral relations for the linear and the second harmonic response for the perturbed Lorenz 63 system, by showing that numerical simulations agree up to high degree of accuracy with the theoretical predictions. Some new theoretical results, showing how to obtain recursively harmonic generation susceptibilities for general observables, are also presented. Our findings confirm the conceptual validity of the nonlinear response theory, suggest that the theory can be extended for more general non equilibrium steady state systems, and shed new light on the applicability of very general tools, based only upon the principle of causality, for diagnosing the behavior of perturbed chaotic systems and reconstructing their output signals, in situations where the fluctuation-dissipation relation is not of great help.

preprint2008arXiv

Symmetry-break, mixing, instability, and low frequency variability in a minimal Lorenz-like system

Starting from the classical Saltzman 2D convection equations, we derive via a severe spectral truncation a minimal 10 ODE system which includes the thermal effect of viscous dissipation. Neglecting this process leads to a dynamical system which includes a decoupled (generalized) Lorenz system. The consideration of this process breaks an important symmetry, couples the dynamics of fast and slow variables, ensuing modifications of the structural properties of the attractor and of the spectral features. When the relevant nondimensional number (Eckert number Ec) is different from zero, the system is ergodic and hyperbolic, the slow variables feature long term memory with 1/f 3/2 power spectra, and the fast variables feature amplitude modulation on time scale of 1/Ec. Increasing the strength of the thermal-viscous feedback has a stabilizing effect, as both the metric entropy and the Kaplan-Yorke attractor dimension decrease monotonically with Ec. The analyzed system features very rich dynamics: it overcomes some of the limitations of the Lorenz system and might have prototypical value in relevant processes in complex systems dynamics, such as the interaction between slow and fast variables. the presence of long term memory and the associated extreme value statistics. Analysis shows how, neglecting the coupling of slow and fast variables only on the basis of scale analysis can be catastrophic. In fact, this leads to spurious invariances that affect essential dynamical properties (ergodicity, hyperbolicity) and that cause the model losing ability in describing intrinsically multiscale processes.

preprint2008arXiv

Three-dimensional Random Voronoi Tessellations: From Cubic Crystal Lattices to Poisson Point Processes

We perturb the SC, BCC, and FCC crystal structures with a spatial Gaussian noise whose adimensional strength is controlled by the parameter a, and analyze the topological and metrical properties of the resulting Voronoi Tessellations (VT). The topological properties of the VT of the SC and FCC crystals are unstable with respect to the introduction of noise, because the corresponding polyhedra are geometrically degenerate, whereas the tessellation of the BCC crystal is topologically stable even against noise of small but finite intensity. For weak noise, the mean area of the perturbed BCC and FCC crystals VT increases quadratically with a. In the case of perturbed SCC crystals, there is an optimal amount of noise that minimizes the mean area of the cells. Already for a moderate noise (a>0.5), the properties of the three perturbed VT are indistinguishable, and for intense noise (a>2), results converge to the Poisson-VT limit. Notably, 2-parameter gamma distributions are an excellent model for the empirical of of all considered properties. The VT of the perturbed BCC and FCC structures are local maxima for the isoperimetric quotient, which measures the degre of sphericity of the cells, among space filling VT. In the BCC case, this suggests a weaker form of the recentluy disproved Kelvin conjecture. Due to the fluctuations of the shape of the cells, anomalous scalings with exponents >3/2 is observed between the area and the volumes of the cells, and, except for the FCC case, also for a->0. In the Poisson-VT limit, the exponent is about 1.67. As the number of faces is positively correlated with the sphericity of the cells, the anomalous scaling is heavily reduced when we perform powerlaw fits separately on cells with a specific number of faces.

preprint2007arXiv

From symmetry break to Poisson point process in 2D Voronoi tessellations: the generic nature of hexagons

We bridge the properties of the regular square and honeycomb Voronoi tessellations of the plane to those of the Poisson-Voronoi case, thus analyzing in a common framework symmetry-break processes and the approach to uniformly random distributions of tessellation-generating points. We consider ensemble simulations of tessellations generated by points whose regular positions are perturbed through a Gaussian noise controlled by the parameter alpha. We study the number of sides, the area, and the perimeter of the Voronoi cells. For alpha>0, hexagons are the most common class of cells, and 2-parameter gamma distributions describe well the statistics of the geometrical characteristics. The symmetry break due to noise destroys the square tessellation, whereas the honeycomb hexagonal tessellation is very stable and all Voronoi cells are hexagon for small but finite noise with alpha<0.1. For a moderate amount of Gaussian noise, memory of the specific unperturbed tessellation is lost, because the statistics of the two perturbed tessellations is indistinguishable. When alpha>2, results converge to those of Poisson-Voronoi tessellations. The properties of n-sided cells change with alpha until the Poisson-Voronoi limit is reached for alpha>2. The Desch law for perimeters is confirmed to be not valid and a square root dependence on n is established. The ensemble mean of the cells area and perimeter restricted to the hexagonal cells coincides with the full ensemble mean; this might imply that the number of sides acts as a thermodynamic state variable fluctuating about n=6; this reinforces the idea that hexagons, beyond their ubiquitous numerical prominence, can be taken as generic polygons in 2D Voronoi tessellations.

preprint2007arXiv

Response Theory for Equilibrium and Non-Equilibrium Statistical Mechanics: Causality and Generalized Kramers-Kronig relations

We consider the general response theory proposed by Ruelle for describing the impact of small perturbations to the non-equilibrium steady states resulting from Axiom A dynamical systems. We show that the causality of the response functions allows for writing a set of Kramers-Kronig relations for the corresponding susceptibilities at all orders of nonlinearity. Nonetheless, only a special class of observable susceptibilities obey Kramers-Kronig relations. Specific results are provided for arbitrary order harmonic response, which allows for a very comprehensive Kramers-Kronig analysis and the establishment of sum rules connecting the asymptotic behavior of the susceptibility to the short-time response of the system. These results generalize previous findings on optical Hamiltonian systems and simple mechanical models, and shed light on the general impact of considering the principle of causality for testing self-consistency: the described dispersion relations constitute unavoidable benchmarks for any experimental and model generated dataset. In order to connect the response theory for equilibrium and non equilibrium systems, we rewrite the classical results by Kubo so that response functions formally identical to those proposed by Ruelle, apart from the measure involved in the phase space integration, are obtained. We briefly discuss how these results, taking into account the chaotic hypothesis, might be relevant for climate research. In particular, whereas the fluctuation-dissipation theorem does not work for non-equilibrium systems, because of the non-equivalence between internal and external fluctuations, Kramers-Kronig relations might be more robust tools for the definition of a self-consistent theory of climate change.

preprint2006arXiv

Does the Danube exist? Versions of reality given by various regional climate models and climatological datasets

We present an intercomparison and verification analysis of several regional climate models (RCMs) nested into the same run of the same Atmospheric Global Circulation Model (AGCM) regarding their representation of the statistical properties of the hydrological balance of the Danube river basin for 1961-1990. We also consider the datasets produced by the driving AGCM, from the ECMWF and NCEP-NCAR reanalyses. The hydrological balance is computed by integrating the precipitation and evaporation fields over the area of interest. Large discrepancies exist among RCMs for the monthly climatology as well as for the mean and variability of the annual balances, and only few datasets are consistent with the observed discharge values of the Danube at its Delta, even if the driving AGCM provides itself an excellent estimate. Since the considered approach relies on the mass conservation principle and bypasses the details of the air-land interface modeling, we propose that the atmospheric components of RCMs still face difficulties in representing the water balance even on a relatively large scale. Their reliability on smaller river basins may be even more problematic. Moreover, since for some models the hydrological balance estimates obtained with the runoff fields do not agree with those obtained via precipitation and evaporation, some deficiencies of the land models are also apparent. NCEP-NCAR and ERA-40 reanalyses result to be largely inadequate for representing the hydrology of the Danube river basin, both for the reconstruction of the long-term averages and of the seasonal cycle, and cannot in any sense be used as verification. We suggest that these results should be carefully considered in the perspective of auditing climate models and assessing their ability to simulate future climate changes.

preprint2006arXiv

Intercomparison of the northern hemisphere winter mid-latitude atmospheric variability of the IPCC models

We compare, for the overlapping time frame 1962-2000, the estimate of the northern hemisphere (NH) mid-latitude winter atmospheric variability within the XX century simulations of 17 global climate models (GCMs) included in the IPCC-4AR with the NCEP and ECMWF reanalyses. We compute the Hayashi spectra of the 500hPa geopotential height fields and introduce an integral measure of the variability observed in the NH on different spectral sub-domains. Only two high-resolution GCMs have a good agreement with reanalyses. Large biases, in most cases larger than 20%, are found between the wave climatologies of most GCMs and the reanalyses, with a relative span of around 50%. The travelling baroclinic waves are usually overestimated, while the planetary waves are usually underestimated, in agreement with previous studies performed on global weather forecasting models. When comparing the results of various versions of similar GCMs, it is clear that in some cases the vertical resolution of the atmosphere and, somewhat unexpectedly, of the adopted ocean model seem to be critical in determining the agreement with the reanalyses. The GCMs ensemble is biased with respect to the reanalyses but is comparable to the best 5 GCMs. This study suggests serious caveats with respect to the ability of most of the presently available GCMs in representing the statistics of the global scale atmospheric dynamics of the present climate and, a fortiori, in the perspective of modelling climate change.

preprint2005arXiv

Physical and Mathematical Properties of a Quasi-Geostrophic Model of Intermediate Complexity of the Mid-Latitudes Atmospheric Circulation

A quasi-geostrophic intermediate complexity model is considered, providing a schematic representation of the baroclinic conversion processes which characterize the physics of the mid-latitudes atmospheric circulation. The model is relaxed towards a given latitudinal temperature profile, which acts as baroclinic forcing, controlled by a parameter TE determining the forced equator-to-pole temperature gradient. As TE increases, a transition takes place from a stationary regime to a periodic regime, and eventually to an earth-like chaotic regime where evolution takes place on a strange attractor. The dependence of the attractor dimension, metric entropy, and bounding box volume in phase space is studied by varying both TE and model resolution. The statistical properties of observables having physical relevance, namely the total energy of the system and the latitudinally averaged zonal wind, are also examined. It is emphasized that while the attractor's properties are quite sensitive to model resolution, the global physical observables depend less critically on it. For more detailed physical observables, such as the latitudinal profiles of the zonal wind, model resolution again may be critical: the effectiveness of the zonal wind convergence, acting as barotropic stabilization of the baroclinic waves, heavily relies on the details of the latitudinal structure of the fields. The necessity and complementarity of both the dynamical systems and physical approach is underlined.

Valerio Lucarini

What is connected

Connect this record

See the researcher in context

Building this map preview

68 published item(s)

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

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

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

Robustness of competing climatic states

Typicality of the 2021 Western North America Summer Heatwave

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

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

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

Introduction to the Special Issue on the Statistical Mechanics of Climate

The Physics of Climate Variability and Climate Change

A New Mathematical Framework for Atmospheric Blocking Events

Response and Sensitivity Using Markov Chains

Dynamical Analysis of Blocking Events: Spatial and Temporal Fluctuations of Covariant Lyapunov Vectors

Extremes and Recurrence in Dynamical Systems

Parametrization of stochastic multiscale triads

Predicting Climate Change using Response Theory: Global Averages and Spatial Patterns

Statistical and Dynamical Properties of Covariant Lyapunov Vectors in a Coupled Atmosphere-Ocean Model - Multiscale Effects, Geometric Degeneracy, and Error Dynamics

Covariant Lyapunov Vectors of a Quasi-geostrophic Baroclinic Model: Analysis of Instabilities and Feedbacks

Early 21st century snow cover state over the western river basins of the Indus River system

Global instability in the Ghil-Sellers model

Prevailing climatic trends and runoff response from Hindukush-Karakoram-Himalaya, upper Indus basin

Response Operators for Markov Processes in a Finite State Space: Radius of Convergence and Link to the Response Theory for Axiom A Systems

Seasonal cycle of Precipitation over Major River Basins in South and Southeast Asia: A Review of the CMIP5 climate models data for present climate and future climate projections

Seasonality of the hydrological cycle in major South and Southeast Asian river basins as simulated by PCMDI/CMIP3 experiments

A new framework for climate sensitivity and prediction: a modelling perspective

Analysis of rainfall seasonality from observations and climate models

Climate of Earth-like planets with high obliquity and eccentric orbits: implications for habitability conditions

Equivalence of Non-Equilibrium Ensembles and Representation of Friction in Turbulent Flows: The Lorenz 96 Model

Mathematical and Physical Ideas for Climate Science

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

Projected changes of rainfall seasonality and dry spells in a high concentration pathway 21st century scenario

The Impact of Oceanic Heat Transport on the Atmospheric Circulation: a Thermodynamic Perspective

Avalanches, Breathers and Flow Reversal in a Continuous Lorenz-96 Model

Elements of a unified framework for response formulae

Entropy Production and Coarse Graining of the Climate Fields in a General Circulation Model

Extreme value statistics for dynamical systems with noise

Habitability and Multistability in Earth-like Planets

Hydrological Cycle over South and Southeast Asian River Basins as Simulated by PCMDI/CMIP3 Experiments

Nonequilibrium thermodynamics of circulation regimes in optically-thin, dry atmospheres

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

Beyond the linear Fluctuation-Dissipation Theorem: the Role of Causality

Bistability of the climate around the habitable zone: a thermodynamic investigation

Multi-level Dynamical Systems: Connecting the Ruelle Response Theory and the Mori-Zwanzig Approach

Nambu representation of an extended Lorenz model with viscous heating

Baroclinic stationary waves in aquaplanet models

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

Climatology of extratropical atmospheric wave packets in the northern hemisphere

Disentangling multi-level systems: averaging, correlations and memory

Extreme Value distribution for singular measures

Hydrological Cycle in the Danube basin in present-day and XXII century simulations by IPCCAR4 global climate models

Modelling Complexity: the case of Climate Science

New Results on the Thermodynamical Properties of the Climate System

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

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

Stochastic perturbations to dynamical systems: a response theory approach

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

A Statistical Mechanical Approach for the Computation of the Climatic Response to General Forcings

Energetics of PCMDI/CMIP3 Climate Models: Energy Budget and Meridional Enthalpy Transport

Thermodynamic Analysis of Snowball Earth Hysteresis Experiment: Efficiency, Entropy Production, and Irreversibility

Thermodynamic Efficiency and Entropy Production in the Climate System

Evidence of dispersion relations for the nonlinear response of the Lorenz 63 system

Symmetry-break, mixing, instability, and low frequency variability in a minimal Lorenz-like system

Three-dimensional Random Voronoi Tessellations: From Cubic Crystal Lattices to Poisson Point Processes

From symmetry break to Poisson point process in 2D Voronoi tessellations: the generic nature of hexagons

Response Theory for Equilibrium and Non-Equilibrium Statistical Mechanics: Causality and Generalized Kramers-Kronig relations

Does the Danube exist? Versions of reality given by various regional climate models and climatological datasets

Intercomparison of the northern hemisphere winter mid-latitude atmospheric variability of the IPCC models

Physical and Mathematical Properties of a Quasi-Geostrophic Model of Intermediate Complexity of the Mid-Latitudes Atmospheric Circulation