A delay differential model of ENSO variability, Part 2: Phase locking, multiple solutions, and dynamics of extrema

We consider a highly idealized model for El Nino/Southern Oscillation (ENSO) variability, as introduced in an earlier paper. The model is governed by a delay differential equation for sea surface temperature in the Tropical Pacific, and it combines two key mechanisms that participate in ENSO dynamics: delayed negative feedback and seasonal forcing. We perform a theoretical and numerical study of the model in the three-dimensional space of its physically relevant parameters: propagation period of oceanic waves across the Tropical Pacific, atmosphere-ocean coupling, and strength of seasonal forcing. Phase locking of model solutions to the periodic forcing is prevalent: the local maxima and minima of the solutions tend to occur at the same position within the seasonal cycle. Such phase locking is a key feature of the observed El Nino (warm) and La Nina (cold) events. The phasing of the extrema within the seasonal cycle depends sensitively on model parameters when forcing is weak. We also study co-existence of multiple solutions for fixed model parameters and describe the basins of attraction of the stable solutions in a one-dimensional space of constant initial model histories.Comment: Nonlin. Proc. Geophys., 2010, accepte

A delay differential model of ENSO variability: Parametric instability and the distribution of extremes

We consider a delay differential equation (DDE) model for El-Nino Southern Oscillation (ENSO) variability. The model combines two key mechanisms that participate in ENSO dynamics: delayed negative feedback and seasonal forcing. We perform stability analyses of the model in the three-dimensional space of its physically relevant parameters. Our results illustrate the role of these three parameters: strength of seasonal forcing $b$, atmosphere-ocean coupling $\kappa$, and propagation period $\tau$ of oceanic waves across the Tropical Pacific. Two regimes of variability, stable and unstable, are separated by a sharp neutral curve in the $(b,\tau)$ plane at constant $\kappa$. The detailed structure of the neutral curve becomes very irregular and possibly fractal, while individual trajectories within the unstable region become highly complex and possibly chaotic, as the atmosphere-ocean coupling $\kappa$ increases. In the unstable regime, spontaneous transitions occur in the mean ``temperature'' ({\it i.e.}, thermocline depth), period, and extreme annual values, for purely periodic, seasonal forcing. The model reproduces the Devil's bleachers characterizing other ENSO models, such as nonlinear, coupled systems of partial differential equations; some of the features of this behavior have been documented in general circulation models, as well as in observations. We expect, therefore, similar behavior in much more detailed and realistic models, where it is harder to describe its causes as completely.Comment: 22 pages, 9 figure

Transport on river networks: A dynamical approach

This study is motivated by problems related to environmental transport on river networks. We establish statistical properties of a flow along a directed branching network and suggest its compact parameterization. The downstream network transport is treated as a particular case of nearest-neighbor hierarchical aggregation with respect to the metric induced by the branching structure of the river network. We describe the static geometric structure of a drainage network by a tree, referred to as the static tree, and introduce an associated dynamic tree that describes the transport along the static tree. It is well known that the static branching structure of river networks can be described by self-similar trees (SSTs); we demonstrate that the corresponding dynamic trees are also self-similar. We report an unexpected phase transition in the dynamics of three river networks, one from California and two from Italy, demonstrate the universal features of this transition, and seek to interpret it in hydrological terms.Comment: 38 pages, 15 figure

Predictability of extreme events in a branching diffusion model

We propose a framework for studying predictability of extreme events in complex systems. Major conceptual elements -- hierarchical structure, spatial dynamics, and external driving -- are combined in a classical branching diffusion with immigration. New elements -- observation space and observed events -- are introduced in order to formulate a prediction problem patterned after the geophysical and environmental applications. The problem consists of estimating the likelihood of occurrence of an extreme event given the observations of smaller events while the complete internal dynamics of the system is unknown. We look for premonitory patterns that emerge as an extreme event approaches; those patterns are deviations from the long-term system's averages. We have found a single control parameter that governs multiple spatio-temporal premonitory patterns. For that purpose, we derive i) complete analytic description of time- and space-dependent size distribution of particles generated by a single immigrant; ii) the steady-state moments that correspond to multiple immigrants; and iii) size- and space-based asymptotic for the particle size distribution. Our results suggest a mechanism for universal premonitory patterns and provide a natural framework for their theoretical and empirical study