This paper has two interrelated foci: (i) obtaining stable and efficient data-driven closure models by using a multivariate time series of partial observations from a large-dimensional system; and (ii) comparing these closure models with the optimal closures predicted by the Mori–Zwanzig (MZ) formalism of statistical physics. Multilayer stochastic models (MSMs) are introduced as both a generalization and a time-continuous limit of existing multilevel, regression-based approaches to closure in a data-driven setting; these approaches include empirical model reduction (EMR), as well as more recent multi-layer modeling. It is shown that the multilayer structure of MSMs can provide a natural Markov approximation to the generalized Langevin equation (GLE) of the MZ formalism. A simple correlation-based stopping criterion for an EMR–MSM model is derived to assess how well it approximates the GLE solution. Sufficient conditions are derived on the structure of the nonlinear cross-interactions between the constitutive layers of a given MSM to guarantee the existence of a global random attractor. This existence ensures that no blow-up can occur for a broad class of MSM applications, a class that includes non-polynomial predictors and nonlinearities that do not necessarily preserve quadratic energy invariants. The EMR–MSM methodology is first applied to a conceptual, nonlinear, stochastic climate model of coupled slow and fast variables, in which only slow variables are observed. It is shown that the resulting closure model with energy-conserving nonlinearities efficiently captures the main statistical features of the slow variables, even when there is no formal scale separation and the fast variables are quite energetic. Second, an MSM is shown to successfully reproduce the statistics of a partially observed, generalized Lotka–Volterra model of population dynamics in its chaotic regime. The challenges here include the rarity of strange attractors in the model’s parameter space and the existence of multiple attractor basins with fractal boundaries. The positivity constraint on the solutions’ components replaces here the quadratic-energy–preserving constraint of fluid-flow problems and it successfully prevents blow-up.
The monograph covers the fundamentals and the consequences of extreme geophysical phenomena like asteroid impacts, climatic change, earthquakes, tsunamis, hurricanes, landslides, volcanic eruptions, flooding, and space weather. This monograph also addresses their associated, local and worldwide socio-economic impacts. The understanding and modeling of these phenomena is critical to the development of timely worldwide strategies for the prediction of natural and anthropogenic extreme events, in order to mitigate their adverse consequences. This monograph is unique in as much as it is dedicated to recent theoretical, numerical and empirical developments that aim to improve: (i) the understanding, modeling and prediction of extreme events in the geosciences, and, (ii) the quantitative evaluation of their economic consequences. The emphasis is on coupled, integrative assessment of the physical phenomena and their socio-economic impacts. With its overarching theme, Extreme Events: Observations, Modeling and Economics will be relevant to and become an important tool for researchers and practitioners in the fields of hazard and risk analysis in general, as well as to those with a special interest in climate change, atmospheric and oceanic sciences, seismo-tectonics, hydrology, and space weather.
Groth, Andreas, Patrice Dumas, Michael Ghil, and Stéphane Hallegatte. “
Impacts of natural disasters on a dynamic economy.” In
Extreme Events : Observations, Modeling, and Economics,
edited by Eric Chavez, Michael Ghil, and Jaime Urrutia-Fucugauchi, 343–360. American Geophysical Union and Wiley-Blackwell, 2015.
Abstract This chapter presents a modeling framework for macroeconomic growth dynamics; it is motivated by recent attempts to formulate and study “integrated models” of the coupling between natural and socioeconomic phe nomena. The challenge is to describe the interfaces between human activities and the functioning of the earth system. We examine the way in which this interface works in the presence of endogenous business cycle dynam ics, based on a nonequilibrium dynamic model. Recent findings about the macroeconomic response to natural disasters in such a nonequilibrium setting have shown a more severe response to natural disasters during expan sions than during recessions. These findings raise questions about the assessment of climate change damages or natural disaster losses that are based purely on long-term growth models. In order to compare the theoretical findings with observational data, we analyze cyclic behavior in the U.S. economy, based on multivariate singular spectrum analysis. We analyze a total of nine aggregate indicators in a 52 year interval (1954–2005) and demon strate that the behavior of the U.S. economy changes significantly between intervals of growth and recession, with higher volatility during expansions.
PDFRecent estimates of climate evolution over the coming century still differ by several degrees. This uncertainty motivates the work presented here. There are two basic approaches to apprehend the complexity of climate change: deterministically nonlinear and stochastically linear, i.e., the Lorenz and the Hasselmann approach. The grand unification of these two approaches relies on the theory of random dynamical systems. We apply this theory to study the random attractors of nonlinear, stochastically perturbed climate models. Doing so allows one to examine the interaction of internal climate variability with the forcing, whether natural or anthropogenic, and to take into account the climate system's non-equilibrium behavior in determining climate sensitivity. This non-equilibrium behavior is due to a combination of nonlinear and random effects. We give here a unified treatment of such effects from the point of view of the theory of dynamical systems and of their bifurcations. Energy balance models are used to illustrate multiple equilibria, while multi-decadal oscillations in the thermohaline circulation illustrate the transition from steady states to periodic behavior. Random effects are introduced in the setting of random dynamical systems, which permit a unified treatment of both nonlinearity and stochasticity. The combined treatment of nonlinear and random effects is applied to a stochastically perturbed version of the classical Lorenz convection model. Climate sensitivity is then defined mathematically as the derivative of an appropriate functional or other function of the system’s state with respect to the bifurcation parameter. This definition is illustrated by using numerical results for a model of the El Niño–Southern Oscillation. The concept of a hierarchy of models is the thread that runs across this chapter, and the robustness of elementary bifurcations across such a hierarchy is emphasized.
PDFSingular spectrum analysis (SSA) along with its multivariate extension (M-SSA) provides an efficient way to identify weak oscillatory behavior in high-dimensional data. To prevent the misinterpretation of stochastic fluctuations in short time series as oscillations, Monte Carlo (MC)–type hypothesis tests provide objective criteria for the statistical significance of the oscillatory behavior. Procrustes target rotation is introduced here as a key method for refining previously available MC tests. The proposed modification helps reduce the risk of type-I errors, and it is shown to improve the test’s discriminating power. The reliability of the proposed methodology is examined in an idealized setting for a cluster of harmonic oscillators immersed in red noise. Furthermore, the common method of data compression into a few leading principal components, prior to M-SSA, is reexamined, and its possibly negative effects are discussed. Finally, the generalized Procrustes test is applied to the analysis of interannual variability in the North Atlantic’s sea surface temperature and sea level pressure fields. The results of this analysis provide further evidence for shared mechanisms of variability between the Gulf Stream and the North Atlantic Oscillation in the interannual frequency band.
PDFWe apply multivariate singular spectrum analysis to the study of U.S. business cycle dynamics. This method provides a robust way to identify and reconstruct oscillations, whether intermittent or modulated. We show such oscillations to be associated with comovements across the entire economy. The problem of spurious cycles generated by the use of detrending filters is addressed and we present a Monte Carlo test to extract significant oscillations. The behavior of the U.S. economy is shown to change significantly from one phase of the business cycle to another: the recession phase is dominated by a five-year mode, while the expansion phase exhibits more complex dynamics, with higher-frequency modes coming into play. We show that the variations so identified cannot be generated by random shocks alone, as assumed in ‘real’ business-cycle models, and that endogenous, deterministically generated variability has to be involved.
PDF
In this second volume, a general approach is developed to provide approximate parameterizations of the "small" scales by the "large" ones for a broad class of stochastic partial differential equations (SPDEs). This is accomplished via the concept of parameterizing manifolds (PMs), which are stochastic manifolds that improve, for a given realization of the noise, in mean square error the partial knowledge of the full SPDE solution when compared to its projection onto some resolved modes. Backward-forward systems are designed to give access to such PMs in practice. The key idea consists of representing the modes with high wave numbers as a pullback limit depending on the time-history of the modes with low wave numbers. Non-Markovian stochastic reduced systems are then derived based on such a PM approach. The reduced systems take the form of stochastic differential equations involving random coefficients that convey memory effects. The theory is illustrated on a stochastic Burgers-type equation.
The El-Nino/Southern-Oscillation (ENSO) phenomenon is the most prominent signal of seasonal-to-interannual climate variability. The past 30 years of research have shown that ENSO dynamics is governed, by and large, by the interplay of the nonlinear mechanisms, and that their simplest version can be studied in autonomous or forced delay differential equation (DDE) models. This chapter briefly reviews the results of Ghil et al., Zaliapin and Ghil, and Ghil and Zaliapin and pursues their DDE model analysis by focusing on multiple model solutions for the same parameter values and the dynamics of local extrema. It first introduces the DDE model of ENSO variability, reviews the main theoretical results concerning its solutions, and comments on the appropriate numerical integration methods. Novel results on multiple solutions and their extrema are reported and illustrated. After discussing the model's pullback attractor, the chapter explores parameter dependence in the model over its entire 3D parameter space.
PDF