Common dynamical properties of business cycle fluctuations are studied in a sample of more than 100 countries that represent economic regions from all around the world. We apply the methodology of multivariate singular spectrum analysis (M-SSA) to identify oscillatory modes and to detect whether these modes are shared by clusters of phase- and frequency-locked oscillators. An extension of the M-SSA approach is introduced to help analyze structural changes in the cluster configuration of synchronization. With this novel technique, we are able to identify a common mode of business cycle activity across our sample, and thus point to the existence of a world business cycle. Superimposed on this mode, we further identify several major events that have markedly influenced the landscape of world economic activity in the postwar era. These findings raise therefore questions about assessments of climate change impacts that are based purely on long-term economic growth models. A key conclusion is the importance of endogenous-dynamics e?ects at the interface between natural climate variability and economic fluctuations.

Watch the video online

Spectral analyses of the North Atlantic temperature field in the Simple Ocean Data Analysis (SODA) reanalysis identify prominent and statistically significant interannual oscillations along the Gulf Stream front and in large regions of the North Atlantic. A 7–8-yr oscillatory mode is characterized by a basin-wide southwest-to-northeast–oriented propagation pattern in the sea surface temperature (SST) field. This pattern is found to be linked to a seesaw in the meridional-dipole structure of the zonal wind stress forcing (TAUX). In the subpolar gyre, the SST and TAUX fields of this mode are shown to be in phase opposition, which suggests a cooling effect of the wind stress on the upper ocean layer. Over all, this mode’s temperature field is characterized by a strong equivalent-barotropic component, as shown by covariations in SSTs and sea surface heights, and by phase-coherent behavior of temperature layers at depth with the SST field. Recent improvements of multivariate singular spectrum analysis (M-SSA) help separate spatio-temporal patterns. This methodology is developed further and applied to studying the ocean’s response to variability in the atmospheric forcing. Statistical evidence is shown to exist for other mechanisms generating oceanic variability of similar 7–8-yr periodicity in the Gulf Stream region; the latter variability is likewise characterized by a strongly equivalent-barotropic component. Two other modes of biennial variability in the Gulf Stream region are also identified, and it is shown that interannual variability in this region cannot be explained by the ocean’s response to similar variability in the atmospheric forcing alone.

The large-scale, near-surface flow of the mid-latitude oceans is dominated by the presence of a larger, anticyclonic and a smaller, cyclonic gyre. The two gyres share the eastward extension of western boundary currents, such as the Gulf Stream or Kuroshio, and are induced by the shear in the winds that cross the respective ocean basins. This physical phenomenology is described mathematically by a hierarchy of systems of nonlinear partial differential equations (PDEs). We study the low-frequency variability of this wind-driven, double-gyre circulation in mid-latitude ocean basins, subject to time-constant, purely periodic and more general forms of time-dependent wind stress. Both analytical and numerical methods of dynamical systems theory are applied to the PDE systems of interest. Recent work has focused on the application of non-autonomous and random forcing to double-gyre models. We discuss the associated pullback and random attractors and the non-uniqueness of the invariant measures that are obtained. The presentation moves from observations of the geophysical phenomena to modeling them and on to a proper mathematical understanding of the models thus obtained. Connections are made with the highly topical issues of climate change and climate sensitivity.

We describe a new approach that allows for systematic causal attribution of weather and climate-related events, in near-real time. The method is designed so as to facilitate its implementation at meteorological centers by relying on data and methods that are routinely available when numerically forecasting the weather. We thus show that causal attribution can be obtained as a by-product of data assimilation procedures run on a daily basis to update numerical weather prediction (NWP) models with new atmospheric observations; hence, the proposed methodology can take advantage of the powerful computational and observational capacity of weather forecasting centers. We explain the theoretical rationale of this approach and sketch the most prominent features of a ``data assimilation–based detection and attribution'' (DADA) procedure. The proposal is illustrated in the context of the classical three-variable Lorenz model with additional forcing. The paper concludes by raising several theoretical and practical questions that need to be addressed to make the proposal operational within NWP centers.

The present work applies singular spectrum analysis (SSA) to the study of macroeconomic fluctuations in three European countries: Italy, The Netherlands, and the United Kingdom. This advanced spectral method provides valuable spatial and frequency information for multivariate data sets and goes far beyond the classical forms of time domain analysis. In particular, SSA enables us to identify dominant cycles that characterize the deterministic behavior of each time series separately, as well as their shared behavior. We demonstrate its usefulness by analyzing several fundamental indicators of the three countries' real aggregate economy in a univariate, as well as a multivariate setting. Since business cycles are international phenomena, which show common characteristics across countries, our aim is to uncover supranational behavior within the set of representative European economies selected herein. Finally, the analysis is extended to include several indicators from the U.S. economy, in order to examine its influence on the European economies under study and their interrelationships.

This paper addresses the effect of interannual variability in jet stream orientation on weather systems over the North Atlantic basin (NAB). The observational analysis relies on 65 yr of NCEP–NCAR reanalysis (1948–2012). The total daily kinetic energy of the geostrophic wind (GTKE) is taken as a measure of storm activity over the North Atlantic. The NAB is partitioned into four rectangular regions, and the winter average of GTKE is calculated for each quadrant. The spatial GTKE average over all four quadrants shows striking year-to-year variability and is strongly correlated with the North Atlantic Oscillation (NAO).The GTKE strength in the northeast quadrant is closely related to the diffluence angle of the jet stream in the northwest quadrant. To gain insight into the relationship between the diffluence angle and its downstream impact, a quasigeostrophic baroclinic model is used. The results show that an initially zonal jet persists at its initial latitude over 30 days or longer, while a tilted jet propagates meridionally according to the Rossby wave group velocity, unless kept stationary by external forcing.A Gulf Stream–like narrow sea surface temperature (SST) front provides the requisite forcing for an analytical steady-state solution to this problem. This SST front influences the atmospheric jet in the northwest quadrant: it both strengthens the jet and tilts it northward at higher levels, while its effect is opposite at lower levels. Reanalysis data confirm these effects, which are consistent with thermal wind balance. The results suggest that the interannual variability found in the GTKE may be caused by intrinsic variability of the thermal Gulf Stream front.

The recent development of dense and continuously operating Global Navigation Satellite System (GNSS) networks worldwide has led to a significant increase in geodetic data sets that sometimes capture transient-deformation signals. It is challenging, however, to extract such transients of geophysical origin from the background noise inherent to GNSS time series and, even more so, to separate them from other signals, such as seasonal redistributions of geophysical fluid mass loads. In addition, because of the very large number of continuously recording GNSS stations now available, it has become impossible to systematically inspect each time series and visually compare them at all neighboring sites. Here we show that Multichannel Singular Spectrum Analysis (M-SSA), a method derived from the analysis of dynamical systems, can be used to extract transient deformations, seasonal oscillations, and background noise present in GNSS time series. M-SSA is a multivariate, nonparametric, statistical method that simultaneously exploits the spatial and temporal correlations of geophysical fields. The method allows for the extraction of common modes of variability, such as trends with nonconstant slopes and oscillations shared across time series, without a priori hypotheses about their spatiotemporal structure or their noise characteristics. We illustrate this method using synthetic examples and show applications to actual GPS data from Alaska to detect seasonal signals and microdeformation at the Akutan active volcano. The geophysically coherent spatiotemporal patterns of uplift and subsidence thus detected are compared to the results of an idealized model of such processes in the presence of a magma chamber source.

A low-order quasigeostrophic double-gyre ocean model is subjected to an aperiodic forcing that mimics time dependence dominated by interdecadal variability. This model is used as a prototype of an unstable and nonlinear dynamical system with time-dependent forcing to explore basic features of climate change in the presence of natural variability. The study relies on the theoretical framework of nonautonomous dynamical systems and of their pullback attractors (PBAs), that is, of the time-dependent invariant sets attracting all trajectories initialized in the remote past. The existence of a global PBA is rigorously demonstrated for this weakly dissipative nonlinear model. Ensemble simulations are carried out and the convergence to PBAs is assessed by computing the probability density function (PDF) of localization of the trajectories. A sensitivity analysis with respect to forcing amplitude shows that the PBAs experience large modifications if the underlying autonomous system is dominated by small-amplitude limit cycles, while less dramatic changes occur in a regime characterized by large-amplitude relaxation oscillations. The dependence of the attracting sets on the choice of the ensemble of initial states is then analyzed. Two types of basins of attraction coexist for certain parameter ranges; they contain chaotic and nonchaotic trajectories, respectively. The statistics of the former does not depend on the initial states whereas the trajectories in the latter converge to small portions of the global PBA. This complex scenario requires separate PDFs for chaotic and nonchaotic trajectories. General implications for climate predictability are finally discussed.

This article revisits the approximation problem of systems of nonlinear delay differential equations (DDEs) by a set of ordinary differential equations (ODEs). We work in Hilbert spaces endowed with a natural inner product including a point mass, and introduce polynomials orthogonal with respect to such an inner product that live in the domain of the linear operator associated with the underlying DDE. These polynomials are then used to design a general Galerkin scheme for which we derive rigorous convergence results and show that it can be numerically implemented via simple analytic formulas. The scheme so obtained is applied to three nonlinear DDEs, two autonomous and one forced: (i) a simple DDE with distributed delays whose solutions recall Brownian motion; (ii) a DDE with a discrete delay that exhibits bimodal and chaotic dynamics; and (iii) a periodically forced DDE with two discrete delays arising in climate dynamics. In all three cases, the Galerkin scheme introduced in this article provides a good approximation by low-dimensional ODE systems of the DDE's strange attractor, as well as of the statistical features that characterize its nonlinear dynamics.

Evaluating the effects of climate variation on ecosystems is of paramount importance for our ability to forecast and mitigate the consequences of global change. However, the ways in which complex food webs respond to climate variations remain poorly understood. Here, we use long-term time series to investigate the effects of temperature variation on the intraguild-predation (IGP) system of Windermere (UK), a lake where pike (Esox lucius, top predator) feed on small-sized perch (Perca fluviatilis) but compete with large-sized perch for the same food sources. Spectral analyses of time series reveal that pike recruitment dynamics are temperature controlled. In 1976, expansion of a size-truncating perch pathogen into the lake severely impacted large perch and favoured pike as the IGP-dominant species. This pathogen-induced regime shift to a pike-dominated IGP apparently triggered a temperature-controlled trophic cascade passing through pike down to dissolved nutrients. In simple food chains, warming is predicted to strengthen top–down control by accelerating metabolic rates in ectothermic consumers, while pathogens of top consumers are predicted to dampen this top–down control. In contrast, the local IGP structure in Windermere made warming and pathogens synergistic in their top–down effects on ecosystem functioning. More generally, our results point to top predators as major mediators of community response to global change, and show that size-selective agents (e.g. pathogens, fishers or hunters) may change the topological architecture of food webs and alter whole ecosystem sensitivity to climate variation.

Abstract The Rosenzweig–MacArthur model is a set of ordinary differential equations (ODEs) that provides an aggregate description of the dynamics of a predator–prey system. When including an Allee effect on the prey, this model exhibits bistability and contains a pitchfork bifurcation, a Hopf bifurcation and a heteroclinic bifurcation. We develop an agent-based model (ABM) on a two-dimensional, square lattice that encompasses the key assumptions of the aggregate model. Although the two modelling approaches – \ODE\ and \ABM\ – differ, both models exhibit similar bifurcation patterns. The \ABM\ model's behaviour is richer and it is analysed using advanced statistical methods. In particular, singular spectrum analysis is used to robustly locate the transition between apparently random, small-amplitude fluctuations around a fixed point and stable, large-amplitude oscillations. Critical slowing down of model trajectories anticipates the heteroclinic bifurcation. Systematic comparison between the \ABM\ and the \ODE\ models’ behaviour helps one understand the predator–prey system better; it provides guidance in model exploration and allows one to draw more robust conclusions on the nature of predator–prey interactions.