Predictability of Precipitation from Continental Radar Images. Part IV: Limits to Prediction

Total Page:16

File Type:pdf, Size:1020Kb

Predictability of Precipitation from Continental Radar Images. Part IV: Limits to Prediction 2092 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 63 Predictability of Precipitation from Continental Radar Images. Part IV: Limits to Prediction URS GERMANN MeteoSwiss, Locarno-Monti, Switzerland ISZTAR ZAWADZKI AND BARRY TURNER Atmospheric and Oceanic Sciences, McGill University, Montréal, and J. S. Marshall Weather Radar Observatory, Sainte Anne de Bellevue, Québec, Canada (Manuscript received 6 June 2005, in final form 9 December 2005) ABSTRACT Predictability of precipitation is examined from storm to synoptic scales through an experimental ap- proach using continent-scale radar composite images. The lifetime of radar reflectivity patterns in Eulerian and Lagrangian coordinates is taken as a measure of predictability. The results are stratified according to scale, location, and time in order to determine how predictability depends on these parameters. Three companion papers give a detailed description of the methodology, and present results are obtained for 143 hours of North American warm season rainfall with emphasis on lifetime, scale dependence, optimum smoothing of forecast fields, and predictability in terms of probabilistic rainfall rates. This paper discusses the sources of forecast uncertainty and extends the analysis to a total of 1424 hours of rainfall. In a Lagrangian persistence framework the predictability problem can be separated into a component associated with growth of precipitation and a component associated with changes in the storm motion field. The role of changes in the motion field turned out to be small but not negligible. A stratifi- cation of lifetime according to location reveals the regions with high predictability and significant nonsta- tionary storm motion. This work is of high practical significance for three reasons: First, Lagrangian persistence of radar patterns was proved to have skill for probabilistic precipitation nowcasting. The discussion of the sources of uncertainty provides a guideline for further improvements. Second, a scale- and location-dependent benchmark is obtained against which the progress of other precipitation forecasting techniques can be evaluated. And, third, the experimental approach to predictability presented in this paper is a valuable contribution to the fundamental question of predictability of precipitation. 1. Introduction dimensional system it can be investigated analytically. But, as soon as the system becomes complicated, as is From a purely dynamic point of view predictability is the case for the atmosphere, or even more pronounced related to the sensitive dependence of the trajectory of for precipitation, it is difficult to study its predictability. the atmospheric state in phase space to small perturba- Then, any approach can only give an estimate of the tions in the initial conditions. No matter how small the true predictability. difference between two initial states is, it will grow and There is a variety of studies that allow inferences of eventually be so large that the two atmospheric states the predictability of the atmosphere: some follow a are no more similar than any randomly chosen pair of purely analytical approach, others are more of experi- possible states. This sensitivity is an intrinsic and fun- mental character. Table 1 gives an overview of predict- damental property of any nonlinear system, and is re- ability studies including some selected references. An ferred to as the initial value problem. For a simple low- example of a purely analytical approach is the evalua- tion of the Liouville equation to examine perturbation Corresponding author address: Dr. Urs Germann, MeteoSwiss, growth of a dynamic system. The Liouville equation, CH-6605 Locarno-Monti, Switzerland. introduced in a meteorological context by Gleeson E-mail: [email protected] (1966) and Epstein (1969), describes the evolution of © 2006 American Meteorological Society Unauthenticated | Downloaded 10/03/21 05:02 AM UTC JAS3735 AUGUST 2006 G E R MANN ET AL. 2093 TABLE 1. Predictability studies. Purely experimental approach Studies with observations: Autocorrelation and spectral analysis (Lorenz 1973; Zawadzki 1973) Analogs and periodicity (Lorenz 1993), interdependence (Zawadzki et al. 1981; Lilly 1986) Hovmöller diagram (Carbone et al. 2002; Ahijevych et al. 2004), composite analysis Precursor and predictor analysis (Bocquet 2002) Studies with statistical models: Forecast skill as evaluated against observations: Eulerian and Lagrangian persistence (Zawadzki et al. 1994; Germann and Zawadzki 2002) Dependence on scale (Lorenz 1969; Germann and Zawadzki 2002; Turner et al. 2004) Dependence on location (this paper) and weather Studies with numerical weather prediction models: Precursor analysis (Massacand et al. 1998) Forecast skill as evaluated against observations Sensitivity to initial conditions: singular vectors (Ehrendorfer et al. 1999), ensembles (Palmer 2002; Buizza et al. 1999; Walser et al. 2004) Studies with idealized systems of dynamic equations: Liouville equation (Gleeson 1966; Epstein 1969; Ehrendorfer 1994a,b) Lyapunov exponent (Hergarten 2002) Lorenz model, attractor (Lorenz 1963, 1993; Palmer 1993; Hergarten 2002) Purely analytical approach probability density of the atmospheric state in phase in the model structure, lack of resolution, errors in space and is the analog to the equation of mass conser- scale interactions, inadequate parameterization, param- vation in physical space. It states that realizations of the eter uncertainty, problems in boundary conditions if state of a system cannot spontaneously appear or dis- running a regional model, and numerical and compu- appear. Principally, it can be solved analytically for any tational errors that result, for instance, from finite dynamic system but, in practice, it is only applicable to difference schemes. In practice it is difficult to sepa- low-dimensional problems, such as a one-dimensional rate the problem of predictability into a component Riccati equation (Ehrendorfer 1994a). associated with initial error and a component associ- On the other end of the spectrum of predictability ated with model error (Palmer 2002). The validity of studies are those purely based on observations, such as the results of predictability studies using numerical autocorrelation analysis of the variable of interest weather prediction models can be very sensitive to (Lorenz 1973; Zawadzki 1973) or studies looking for model errors. precursors in the observation space (Bocquet 2002). Zängl (2004) found sensitivity of simulated oro- Between the purely analytical and purely experimental graphic rainfall to model components other than cloud approach there is a variety of studies employing statis- microphysics. He compared the sensitivity to param- tical and numerical weather prediction models that eterization for cloud microphysics with the sensitivity to combine analytical concepts with observations. In a the implementation of horizontal diffusion, the type of model approach predictability can be measured by the vertical coordinates, and boundary layer parameteriza- skill of the model to predict a certain phenomenon as tions. The experiment was conducted with the fifth- evaluated by statistical comparison with observations, generation Pennsylvania State University–National or by the growth of initially small perturbations in Center for Atmospheric Research (NCAR) Mesoscale model phase space using singular vector analysis or en- Model (MM5) and four nested domains with a horizon- sembles of model runs. tal mesh size being 37.8, 12.6, 4.2, and 1.4 km. Precipi- The analytical approach allows one to understand tation of the 1.4-km run accumulated over 36 h and basic concepts of nonlinearity, perturbation growth, averaged over a domain with a 200-km side length var- and sensitive dependence. But, as studies with idealized ied by a factor of 1.4 when using two different imple- systems of dynamic equations or with state-of-the-art mentations of horizontal diffusion. This and other re- numerical weather prediction models both assume the sults indicate that errors in simulated precipitation model to be representative of the true system, the crux fields are not only the result of weaknesses in the cloud lies in all sorts of model errors. These include errors microphysics scheme but may also arise from model Unauthenticated | Downloaded 10/03/21 05:02 AM UTC 2094 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 63 components that are not directly related to precipita- weather in order to determine the dependence of pre- tion such as the diffusion scheme. dictability on these parameters. Another difficulty comes from the dimension and A similar stratification of predictability is not complexity of the weather system, which makes analyti- straightforward when using a numerical weather pre- cal approaches computationally expensive or impracti- diction model because design and implementation of a cal. To obtain a representative estimate of sensitive model both have a significant influence on the perfor- dependence we may need more than just a handful of mance at different scales, locations, and in different ensemble members calculated for a specific weather weather situations. Of course there are problems in the situation, location, and time. But this is often not pos- experimental approaches as well. First there are mea- sible because of limited computation time. surement errors, and, second, there is some difficulty in linking the results to the nonlinear and chaotic nature a. Predictability of precipitation of the system. There are
Recommended publications
  • WHAT IS a CHAOTIC ATTRACTOR? 1. Introduction J. Yorke Coined the Word 'Chaos' As Applied to Deterministic Systems. R. Devane
    WHAT IS A CHAOTIC ATTRACTOR? CLARK ROBINSON Abstract. Devaney gave a mathematical definition of the term chaos, which had earlier been introduced by Yorke. We discuss issues involved in choosing the properties that characterize chaos. We also discuss how this term can be combined with the definition of an attractor. 1. Introduction J. Yorke coined the word `chaos' as applied to deterministic systems. R. Devaney gave the first mathematical definition for a map to be chaotic on the whole space where a map is defined. Since that time, there have been several different definitions of chaos which emphasize different aspects of the map. Some of these are more computable and others are more mathematical. See [9] a comparison of many of these definitions. There is probably no one best or correct definition of chaos. In this paper, we discuss what we feel is one of better mathematical definition. (It may not be as computable as some of the other definitions, e.g., the one by Alligood, Sauer, and Yorke.) Our definition is very similar to the one given by Martelli in [8] and [9]. We also combine the concepts of chaos and attractors and discuss chaotic attractors. 2. Basic definitions We start by giving the basic definitions needed to define a chaotic attractor. We give the definitions for a diffeomorphism (or map), but those for a system of differential equations are similar. The orbit of a point x∗ by F is the set O(x∗; F) = f Fi(x∗) : i 2 Z g. An invariant set for a diffeomorphism F is an set A in the domain such that F(A) = A.
    [Show full text]
  • Attractors of Dynamical Systems in Locally Compact Spaces
    Open Math. 2019; 17:465–471 Open Mathematics Research Article Gang Li* and Yuxia Gao Attractors of dynamical systems in locally compact spaces https://doi.org/10.1515/math-2019-0037 Received February 7, 2019; accepted March 4, 2019 Abstract: In this article the properties of attractors of dynamical systems in locally compact metric space are discussed. Existing conditions of attractors and related results are obtained by the near isolating block which we present. Keywords: Flow, attraction neighborhood, near isolating block MSC: 34C35 54H20 1 Introduction The dynamical system theory studies the rules of changes in the state which depends on time. In the investigation of dynamical systems, one of very interesting topics is the study of attractors (see [1 − 4] and the references therein). In [1 − 3], the authors gave some denitions of attractors and also made further investigations about the properties of them. In [5], the limit set of a neighborhood was used in the denition of an attractor, and in [6] Hale and Waterman also emphasized the importance of the limit set of a set in the analysis of the limiting behavior of a dynamical system. In [7 − 9], the authors dened intertwining attractor of dynamical systems in metric space and obtained some existing conditions of intertwining attractor. In [10], the author studied the properties of limit sets of subsets and attractors in a compact metric space. In [11], the author studied a positively bounded dynamical system in the plane, and obtained the conditions of compactness of the set of all bounded solutions. In [12], the uniform attractor was dened as the minimal compact uniformly pullback attracting random set, several existence criteria for uniform attractors were given and the relationship between uniform and co-cycle attractors was carefully studied.
    [Show full text]
  • Control of Chaos: Methods and Applications
    Automation and Remote Control, Vol. 64, No. 5, 2003, pp. 673{713. Translated from Avtomatika i Telemekhanika, No. 5, 2003, pp. 3{45. Original Russian Text Copyright c 2003 by Andrievskii, Fradkov. REVIEWS Control of Chaos: Methods and Applications. I. Methods1 B. R. Andrievskii and A. L. Fradkov Institute of Mechanical Engineering Problems, Russian Academy of Sciences, St. Petersburg, Russia Received October 15, 2002 Abstract|The problems and methods of control of chaos, which in the last decade was the subject of intensive studies, were reviewed. The three historically earliest and most actively developing directions of research such as the open-loop control based on periodic system ex- citation, the method of Poincar´e map linearization (OGY method), and the method of time- delayed feedback (Pyragas method) were discussed in detail. The basic results obtained within the framework of the traditional linear, nonlinear, and adaptive control, as well as the neural network systems and fuzzy systems were presented. The open problems concerned mostly with support of the methods were formulated. The second part of the review will be devoted to the most interesting applications. 1. INTRODUCTION The term control of chaos is used mostly to denote the area of studies lying at the interfaces between the control theory and the theory of dynamic systems studying the methods of control of deterministic systems with nonregular, chaotic behavior. In the ancient mythology and philosophy, the word \χαωσ" (chaos) meant the disordered state of unformed matter supposed to have existed before the ordered universe. The combination \control of chaos" assumes a paradoxical sense arousing additional interest in the subject.
    [Show full text]
  • ONE HUNDRED YEARS of COMPLEX DYNAMICS the Subject of Complex Dynamics, That Is, the Behaviour of Orbits of Holomorphic Functions
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by University of Liverpool Repository ONE HUNDRED YEARS OF COMPLEX DYNAMICS MARY REES The subject of Complex Dynamics, that is, the behaviour of orbits of holomorphic functions, emerged in the papers produced, independently, by Fatou and Julia, almost 100 years ago. Although the subject of Dynami- cal Systems did not then have a name, the dynamical properties found for holomorphic systems, even in these early researches, were so striking, so unusually comprehensive, and yet so varied, that these systems still attract widespread fascination, 100 years later. The first distinctive feature of iter- ation of a single holomorphic map f is the partition of either the complex plane or the Riemann sphere into two sets which are totally invariant under f: the Julia set | closed, nonempty, perfect, with dynamics which might loosely be called chaotic | and its complement | open, possibly empty, but, if non-empty, then with dynamics which were completely classified by the two pioneering researchers, modulo a few simply stated open questions. Before the subject re-emerged into prominence in the 1980's, the Julia set was alternately called the Fatou set, but Paul Blanchard introduced the idea of calling its complement the Fatou set, and this was immediately universally accepted. Probably the main reason for the remarkable rise in interest in complex dynamics, about thirty-five years ago, was the parallel with the subject of Kleinian groups, and hence with the whole subject of hyperbolic geome- try. A Kleinian group acting on the Riemann sphere is a dynamical system, with the sphere splitting into two disjoint invariant subsets, with the limit set and its complement, the domain of discontinuity, having exactly similar properties to the Julia and Fatou sets.
    [Show full text]
  • Forward Omega Limit Sets of Nonautonomous Dynamical Systems
    DISCRETE AND CONTINUOUS doi:10.3934/dcdss.2020065 DYNAMICAL SYSTEMS SERIES S Volume 13, Number 4, April 2020 pp. 1103{1114 FORWARD OMEGA LIMIT SETS OF NONAUTONOMOUS DYNAMICAL SYSTEMS Hongyong Cui, Peter E. Kloeden and Meihua Yang School of Mathematics and Statistics, and Hubei Key Laboratory of Engineering Modeling and Scientific Computing Huazhong University of Science and Technology Wuhan 430074, China Dedicated to J¨urgen Scheurle on his 65th birthday. Abstract. The forward !-limit set !B of a nonautonomous dynamical system ' with a positively invariant absorbing family B = fB(t); t 2 Rg of closed and bounded subsets of a Banach space X which is asymptotically compact is shown to be asymptotically positive invariant in general and asymptotic negative invariant if ' is also strongly asymptotically compact and eventually continuous in its initial value uniformly on bounded time sets independently of the initial time. In addition, a necessary and sufficient condition for a '- invariant family A = fA(t); t 2 Rg in B of nonempty compact subsets of X to be a forward attractor is generalised to this context. 1. Introduction. Omega limit sets play a fundamental role in characterising the behaviour of dynamical systems. Indeed, the attractor of an autonomous dynamical system is an omega limit set. In such systems the omega limit sets are invariant. The situation is more complicated in nonautonomous dynamical systems because the behaviour now depends on the actual time and not on the the time that has elapsed since starting. This leads to new concepts of nonautonomous attractors that consist of families of compact subsets.
    [Show full text]
  • Toward a Theory of Chaos
    December 3, 2003 12:13 00851 Tutorials and Reviews International Journal of Bifurcation and Chaos, Vol. 13, No. 11 (2003) 3147{3233 c World Scientific Publishing Company TOWARD A THEORY OF CHAOS A. SENGUPTA Department of Mechanical Engineering, Indian Institute of Technology Kanpur, Kanpur 208016, India [email protected] Received February 23, 2001; Revised September 19, 2002 This paper formulates a new approach to the study of chaos in discrete dynamical systems based on the notions of inverse ill-posed problems, set-valued mappings, generalized and multi- valued inverses, graphical convergence of a net of functions in an extended multifunction space [Sengupta & Ray, 2000], and the topological theory of convergence. Order, chaos and complexity are described as distinct components of this unified mathematical structure that can be viewed as an application of the theory of convergence in topological spaces to increasingly nonlinear mappings, with the boundary between order and complexity in the topology of graphical con- vergence being the region in (Multi(X)) that is susceptible to chaos. The paper uses results from the discretized spectral approximation in neutron transport theory [Sengupta, 1988, 1995] and concludes that the numerically exact results obtained by this approximation of the Case singular eigenfunction solution is due to the graphical convergence of the Poisson and conjugate Poisson kernels to the Dirac delta and the principal value multifunctions respectively. In (Multi(X)), the continuous spectrum is shown to reduce to a point spectrum, and we introduce a notion of latent chaotic states to interpret superposition over generalized eigenfunctions. Along with these latent states, spectral theory of nonlinear operators is used to conclude that nature supports complexity to attain efficiently a multiplicity of states that otherwise would remain unavailable to it.
    [Show full text]
  • A Poincaré-Bendixson Theorem for Hybrid Systems
    A POINCARE-BENDIXSON´ THEOREM FOR HYBRID SYSTEMS WILLIAM CLARK, ANTHONY BLOCH, LEONARDO COLOMBO Abstract. The Poincar´e-Bendixsontheorem plays an important role in the study of the qualitative be- havior of dynamical systems on the plane; it describes the structure of limit sets in such systems. We prove a version of the Poincar´e-BendixsonTheorem for two dimensional hybrid dynamical systems and describe a method for computing the derivative of the Poincar´ereturn map, a useful object for the stability analysis of hybrid systems. We also prove a Poincar´e-BendixsonTheorem for a class of one dimensional hybrid dynamical systems. 1. Introduction H. Poincar´econsidered the problem of characterizing the structure of limit sets of trajectories of analytic vector fields on the plane in 1886 [17]. I. Bendixson improved the solution proposed by Poincar´ein 1903 by solving the problem under the weaker hypothesis of C1 vector fields [2]. Since then, the investigation of the asymptotic behavior of dynamical systems has been essential to understanding their behavior. The theory of Poincar´e-Bendixsonstudies so-called limit sets. The classic version of the Poincar´e-BendixsonTheorem states that if a trajectory is bounded and its limit set does not contain any fixed points, then the limt set is a periodic orbit [16]. Therefore, the problem of determining the existence of limit cycles in planar continuous dynamics is well understood. Hybrid systems [6] are non-smooth dynamical systems which exhibit a combination of smooth and discrete dynamics, where the flow evolves continuously on a state space, and a discrete transition occurs when the flow intersects a co-dimension one hypersurface.
    [Show full text]
  • Some Applications of the Poincaré-Bendixson Theorem Robert Roussarie
    Some Applications of the Poincaré-Bendixson Theorem Robert Roussarie To cite this version: Robert Roussarie. Some Applications of the Poincaré-Bendixson Theorem. 2021. hal-03107562v2 HAL Id: hal-03107562 https://hal.archives-ouvertes.fr/hal-03107562v2 Preprint submitted on 19 Jan 2021 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. SOME APPLICATIONS OF THE POINCARE-BENDIXSON´ THEOREM ROBERT ROUSSARIE Abstract. We consider a C1 vector field X defined on an open subset U of the plane, with compact closure. If X has no singular points and if U is simply connected, a weak version of the Poincar´e-BendixsonTheorem says that the limit sets of X in U are empty but that one can defined non empty extended limit sets contained into the boundary of U. We give an elementary proof of this result, independent of the classical Poincar´e-BendixsonTheorem. A trapping triangle T based at p, for a C1 vector field X defined on an open subset U of the plane, is a topological triangle with a corner at a point p located on the boundary @U and a good control of the tranversality of X along the sides.
    [Show full text]
  • Limit Sets of Foliations Annales De L’Institut Fourier, Tome 15, No 2 (1965), P
    ANNALES DE L’INSTITUT FOURIER RICHARD SACKSTEDER ART J. SCHWARTZ Limit sets of foliations Annales de l’institut Fourier, tome 15, no 2 (1965), p. 201-213 <http://www.numdam.org/item?id=AIF_1965__15_2_201_0> © Annales de l’institut Fourier, 1965, tous droits réservés. L’accès aux archives de la revue « Annales de l’institut Fourier » (http://annalif.ujf-grenoble.fr/) implique l’accord avec les conditions gé- nérales d’utilisation (http://www.numdam.org/conditions). Toute utilisa- tion commerciale ou impression systématique est constitutive d’une in- fraction pénale. Toute copie ou impression de ce fichier doit conte- nir la présente mention de copyright. Article numérisé dans le cadre du programme Numérisation de documents anciens mathématiques http://www.numdam.org/ Ann. Inst. Fourier, Grenoble 15, 2 (1965), 201-214. LIMIT SETS OF FOLIATIONS by Richard SACKSTEDER and Arthur J. SCHWARTZ 1. Introduction. Let V be an M-manifold with a foliated structure of co- dimension one. A leaf of the foliation is called proper if its topology as an (n—l)-manifold agrees with its topology as a subset of V. One type of theorem which is proved here asserts that proper leaves behave much like compact leaves with respect to certain stability properties. For example, Theorem 1 can be viewed as an extension of a theorem of Reeb on the behavior of leaves in a neighborhood of compact leaf to non-compact, proper leaves. Other theorems show that leaves whose holonomy is finite in a certain sense have properties like those of non-periodic solutions of differential equations on 2-manifolds.
    [Show full text]
  • On the Julia Set of a Typical Quadratic Polynomial with a Siegel Disk
    Annals of Mathematics, 159 (2004), 1–52 On the Julia set of a typical quadratic polynomial with a Siegel disk By C. L. Petersen and S. Zakeri To the memory of Michael R. Herman (1942–2000) Abstract Let 0 <θ<1 be an irrational number with continued fraction expansion 2πiθ θ =[a1,a2,a3,...], and consider the quadratic polynomial Pθ : z → e z + z2. By performing a trans-quasiconformal surgery on an associated Blaschke product model, we prove that if √ log an = O( n)asn →∞, then the Julia set of Pθ is locally connected and has Lebesgue measure zero. In particular, it follows that for almost every 0 <θ<1, the quadratic Pθ has a Siegel disk whose boundary is a Jordan curve passing through the critical point of Pθ. By standard renormalization theory, these results generalize to the quadratics which have Siegel disks of higher periods. Contents 1. Introduction 2. Preliminaries 3. A Blaschke model 4. Puzzle pieces and a priori area estimates 5. Proofs of Theorems A and B 6. Appendix: A proof of Theorem C References 1. Introduction 2πiθ 2 Consider the quadratic polynomial Pθ : z → e z + z , where 0 <θ<1 is an irrational number. It has an indifferent fixed point at 0 with multiplier 2πiθ − 2πiθ ∞ Pθ(0) = e , and a unique finite critical point located at e /2. Let Aθ( ) be the basin of attraction of infinity, Kθ = C Aθ(∞) be the filled Julia set, 2 C. L. PETERSEN AND S. ZAKERI and Jθ = ∂Kθ be the Julia set of Pθ. The behavior of the sequence of iterates { ◦n} Pθ n≥0 near Jθ is intricate and highly nontrivial.
    [Show full text]
  • Μ-Limit Sets of Cellular Automata from a Computational Complexity
    µ-Limit Sets of Cellular Automata from a Computational Complexity Perspective Laurent Boyer,Martin Delacourt,Victor Poupet,Mathieu Sablik, Guillaume Theyssier Abstract This paper concerns µ-limit sets of cellular automata: sets of configu- rations made of words whose probability to appear does not vanish with time, starting from an initial µ-random configuration. More precisely, we investigate the computational complexity of these sets and of related 0 decision problems. Main results: first, µ-limit sets can have a Σ3-hard language, second, they can contain only α-complex configurations, third, 0 any non-trivial property concerning them is at least Π3-hard. We prove complexity upper bounds, study restrictions of these questions to particu- lar classes of CA, and different types of (non-)convergence of the measure of a word during the evolution. 1 Introduction A cellular automaton (CA) is a complex system defined by a local rule which acts synchronously and uniformly on the configuration space. These simple models have a wide variety of different dynamical behaviors, in particular interesting asymptotic behaviors. In the dynamical systems context, it is natural to study the limit set of a cellular automaton: it is defined as the set of configurations that can appear arbitrarily far in time. This set captures the longterm behavior of the CA and has been widely studied since the end of the 1980s. Given a cellular automaton, arXiv:1309.6730v2 [cs.DM] 22 Jun 2015 it is difficult to determine its limit set. Indeed it is undecidable to know if it contains only one configuration [Kar92] and more generally, any nontrivial property of limit sets is undecidable [Kar94].
    [Show full text]
  • Arxiv:Math/9410221V1
    APPEARED IN BULLETIN OF THE AMERICAN MATHEMATICAL SOCIETY Volume 31, Number 2, October 1994, Pages 155-172 FRONTIERS IN COMPLEX DYNAMICS CURTIS T. MCMULLEN 1. Introduction Rational maps on the Riemann sphere occupy a distinguished niche in the general theory of smooth dynamical systems. First, rational maps are complex-analytic, so a broad spectrum of techniques can contribute to their study (quasiconformal mappings, potential theory, algebraic geometry, etc.). The rational maps of a given degree form a finite-dimensional manifold, so exploration of this parameter space is especially tractable. Finally, some of the conjectures once proposed for smooth dynamical systems (and now known to be false) seem to have a definite chance of holding in the arena of rational maps. In this article we survey a small constellation of such conjectures centering around the density of hyperbolic rational maps — those which are dynamically the best behaved. We discuss some of the evidence and logic underlying these conjectures, and sketch recent progress towards their resolution. Our presentation entails only a brief account of the basics of complex dynamics; a more systematic exposition can be found in the survey articles [Dou1], [Bl], and [EL]; the recent books [Bea] and [CG]; and Milnor’s lecture notes [Mil2]. 2. Hyperbolic rational maps A rational map f : C → C is a holomorphic dynamical system on the Riemann sphere C = C ∪ {∞}. Any such map can be written as a quotient b b d b P (z) a0z + ... + ad f(z) = = d Q(z) b0z + ... + bd of two relative prime polynomials P and Q. The degree of f can be defined topo- logically or algebraically; it is the number of preimages of a typical point z, as arXiv:math/9410221v1 [math.DS] 1 Oct 1994 well as the maximum of the degrees of P and Q.
    [Show full text]