Modelling Extreme Rainfall in the Era of Climate Change Concerns: Towards a Consistent Stochastic Methodology

Total Page:16

File Type:pdf, Size:1020Kb

Modelling Extreme Rainfall in the Era of Climate Change Concerns: Towards a Consistent Stochastic Methodology School for Young Scientists: “Modelling and forecasting of river flows and managing hydrological risks: towards a new generation of methods” Moscow, Russia, 22-26 October 2018 Modelling extreme rainfall in the era of climate change concerns: Towards a consistent stochastic methodology One Step Forward, Two Steps Back Demetris Koutsoyiannis Department of Water Resources and Environmental Engineering School of Civil Engineering National Technical University of Athens, Greece ([email protected], http://www.itia.ntua.gr/dk/) Presentation available online: http://www.itia.ntua.gr/1897/ Greetings from the Itia research team http://www.itia.ntua.gr/ D. Koutsoyiannis, Modelling extreme rainfall 1 Giants of the Moscow School of Mathematics D. Koutsoyiannis, Modelling extreme rainfall 2 On names and definitions: A contrast Ἀρχὴ παιδεύσεως ἡ τῶν ὀνομάτων ἐπίσκεψις” What’s in a name? That which (The beginning of education is the inspection of we call a rose, by any other names) name would smell as sweet. Attributed to Socrates by Epictetus, Discourses, Ι.17,12 William Shakespeare, “Romeo and Juliet” (Act 2, scene 2) Each definition is a piece of secret ripped from Let me argue that this Nature by the human spirit. I insist on this: any situation [absence of a complicated thing, being illumined by definitions, definition] ought not create being laid out in them, being broken up into concern and steal time from pieces, will be separated Into pieces completely useful work. Entire fields of transparent even to a child, excluding foggy and mathematics thrive for dark parts that our intuition whispers to us while centuries with a clear but acting, separating into logical pieces, then only evolving self-image, and can we move further, towards new successes due nothing resembling a to definitions . definition Nikolai Luzin (from Graham and Kantor, 2009) Benoit Mandelbrot (1999, p. 14) D. Koutsoyiannis, Modelling extreme rainfall 3 Part A Premises for a stochastic framework about change D. Koutsoyiannis, Modelling extreme rainfall 4 The general framework: Seeking theoretical consistency in analysis of geophysical data (Using stochastics) D. Koutsoyiannis, Modelling extreme rainfall 5 Why stochastics in geophysics and hydrology? Geophysics is the branch of physics that relies most decisively on data. Geophysical data are numbers but to treat them we need to use stochastics, not arithmetic. Stochastics is the mathematics of random variables and stochastic processes. Random variables and stochastic processes are abstract mathematical objects, whose properties distinguish them from typical variables that take on numerical values. It is important to understand these properties before making calculations with data, otherwise the results may be meaningless (not even wrong). The numerical data allow us to estimate (not to determine precisely) expectations. Expectations are defined as integrals of products of functions. For a continuous random variable x with probability density function f(x), the expectation of an arbitrary function g of x (where g(x) ∞ ( ) ( ) . is a random variable per se), the expectation of g(x) is defined as 휃 ≔ E[푔(푥)] ≔ ∫−∞ 푔 푥 푓 푥 d푥 Central among expectations are the moments, in which g(x) is a power of x (or a linear expression of x). To estimate true parameters 휃 from data we need estimators; the estimator 휃̂ of θ is a random variable depending on the stochastic process of interest x(t) and is a model per se, not a number. The estimate 휃̂ is a number, calculated by using the observations and the estimator. Characteristic statistics of the estimator 휃̂ are its bias, E[휃̂] − 휃, and its variance var[휃̂]. When E[휃̂] = 휃 the estimator is called unbiased. Estimation is made possible thanks to two concepts of stochastics: stationarity and ergodicity. Stationarity and ergodicity are not incompatible with, or contradictory to, change. D. Koutsoyiannis, Modelling extreme rainfall 6 Change is crucial in Hydrology: ‘Panta Rhei’—The scientific decade of IAHS 2013-2022 http://iahs.info/Commissions--W-Groups/Working-Groups/Panta-Rhei.do D. Koutsoyiannis, Modelling extreme rainfall 7 ‘Panta Rhei’: © Heraclitus Change and randomness Πάντα ῥεῖ Everything flows (Heraclitus; quoted in Plato’s Cratylus, 339-340) Αἰών παῖς ἐστι παίζων πεσσεύων Time is a child playing, throwing dice (Heraclitus; Fragment 52) D. Koutsoyiannis, Modelling extreme rainfall 8 Change, logic, precision: © Aristotle Μεταβάλλει τῷ χρόνῳ πάντα All is changing in the course of time ( Aristotle; Meteorologica, I.14, 353a 16) Λογική, συλλογισμός, επαγωγή, ορθός λόγος Logic, deduction, induction, (right) reason (Aristotle, Organon & Nicomachean Ethics) …τοσοῦτον τἀκριβὲς ἐπιζητεῖν καθ᾽ ἕκαστον γένος, ἐφ᾽ ὅσον ἡ τοῦ πράγματος φύσις ἐπιδέχεται … look for precision in each class of things just so far as the nature of the subject admits (Aristotle, Nicomachean Ethics 1094b) D. Koutsoyiannis, Modelling extreme rainfall 9 Hurst-Kolmogorov dynamics—Or: Earth’s perpetual change 7 Annual Nile River annual minimum Structured 6 30-year average water level (849 values) random 5 4 3 Minimum water depth (m) 2 1 0 600 700 800 900 1000 1100 1200 1300 1400 1500 Year AD Nilometer data: Koutsoyiannis (2013a) 7 "Annual" Each value is the minimum of m=36 roulette wheel Purely 6 30-"year" average outcomes. The value of m was chosen so that the standard deviation be equal to the Nilometer series random 5 4 3 2 1 Minimum roulette wheel outcome 0 600 700 800 900 1000 1100 1200 1300 1400 1500 "Year" D. Koutsoyiannis, Modelling extreme rainfall 10 The climacogram: A simple statistical tool to quantify change across time scales • Take the Nilometer time series, x1, x2, ..., x849, and calculate the sample estimate of variance γ(1), where the superscript (1) indicates time scale (1 year) • Form a time series at time scale 2 (years): (2) 푥1 + 푥2 (2) 푥3 + 푥4 (2) 푥847 + 푥848 푥 ≔ , 푥 ∶= , . , 푥 ∶= (1) 1 2 2 2 424 2 and calculate the sample estimate of the variance γ(2). • Repeat the same procedure and form a time series at time scale 3, 4, … (years), up to scale 84 (1/10 of the record length) and calculate the variances γ(3), γ(4),… γ(84). • The climacogram is the variance γ (κ) as a function of scale κ; it is visualized as a double logarithmic plot of γ (κ) vs. κ (or alternatively of the standard deviation σ(κ)). • If the time series xτ represented a pure random process, the climacogram would be a straight line with slope –1 (the proof is very easy). • In real world processes, the slope is different from –1, designated as 2H – 2, where H is the so-called Hurst coefficient (0 < H < 1). • The scaling law γ(κ) = γ(1) / κ2 – 2H defines the Hurst-Kolmogorov (HK) process. • High values of H (> 0.5) indicate enhanced change at large scales, else known as long-term persistence, or strong clustering (grouping) of similar values. D. Koutsoyiannis, Modelling extreme rainfall 11 The climacogram of the Nilometer time series 7 • The Hurst-Kolmogorov process Annual 6 30-year average seems consistent with reality. 5 • The Hurst coefficient is H = 0.87 4 (Similar H values are estimated 3 from the simultaneous record of Minimum water depth (m) 2 maximum water levels and from 1 0 the modern, 131-year, flow 600 700 800 900 1000 1100 1200 1300 1400 1500 record of the Nile flows at Year AD Aswan). 1 • The Hurst-Kolmogorov behaviour, seen in the climacogram, indicates that: Variance (m²) (a) long-term changes are more frequent and intense 0.1 than commonly perceived, and (b) future states are much Empirical (from data) more uncertain and Purely random (H = 0.5) Markov unpredictable on long time Hurst-Kolmogorov, theoretical (H = 0.87) Hurst-Kolmogorov adapted for bias horizons than implied by pure 0.01 randomness. 1 10 100 Scale (years) D. Koutsoyiannis, Modelling extreme rainfall 12 Change and predictability Change Koutsoyiannis, 2013a Predictable Unpredictable (regular) (random) Non-periodic Periodic Purely random Structured e.g. acceleration of e.g. daily and e.g. consecutive random outcomes of dice e.g. climatic a falling body annual cycles fluctuations Simple systems – Short time horizons Complex systems – Long time horizons Important but trivial Most interesting D. Koutsoyiannis, Modelling extreme rainfall 13 The cause of change: © Peter Atkins Atkins, 2004 Atkins, 2007 D. Koutsoyiannis, Modelling extreme rainfall 14 Entropy ≡ Uncertainty quantified • Historically entropy was introduced in thermodynamics but later it was given a rigorous definition within probability theory (owing to Boltzmann, Gibbs and Shannon). • Thermodynamic and probabilistic entropy are essentially the same thing (Koutsoyiannis, 2010, 2013b, 2014; but others have different opinion). • Entropy acquires its importance from the principle of maximum entropy (Jaynes, 1957), which postulates that the entropy of a random variable should be at maximum, under the conditions (constraints) which incorporate the available information about this variable. • The tendency of entropy to become maximal explains a spectrum of phenomena from the random outcomes of dice to the Second Law of thermodynamics as the driving force of natural change. • Entropy is a dimensionless measure of uncertainty: Discrete random variable z Continuous random variable z 푓(z) ∞ 푓(푧) Φ[z] := E[–ln P(z)] = ∑푤 푃 ln 푃 Φ[z] := E[– ln ] = – ∫ ln 푓(푧)d푧 푗 = 1 푗 푗 ℎ(z) −∞ ℎ(푧) where P ≔ P{z = z } (probability) where f(z) is probability density and h(z) j j is the density of a background measure D. Koutsoyiannis, Modelling extreme rainfall 15 Memorable moments in the history of stochastics Ludwig Boltzmann George D. Birkhoff Aleksandr Khinchin Andrey N. Kolmogorov (1844 –1906, Universities of (1884 – 1944; (1894 – 1959; Moscow (1903 – 1987; Moscow State Graz and Vienna, Austria, and Princeton, Harvard, State University, University, Russia) Munich, Germany) USA) Russia) 1877 Explanation of the 1931 Discovery of 1933 Purely measure- 1931 Introduction of the terms concept of entropy in the ergodic theoretic proof of the process to describe change of a probability theoretic context.
Recommended publications
  • From Time Series to Stochastics: a Theoretical Framework with Applications on Time Scales Spanning from Microseconds to Megayears
    Orlob International Symposium 2016 University California Davis 20-21 June 2016 From time series to stochastics: A theoretical framework with applications on time scales spanning from microseconds to megayears Demetris Koutsoyiannis & Panayiotis Dimitriadis Department of Water Resources and Environmental Engineering School of Civil Engineering National Technical University of Athens, Greece ([email protected], http://www.itia.ntua.gr/dk/) Presentation available online: http://www.itia.ntua.gr/en/docinfo/1618/ Part 1: On names and definitions D. Koutsoyiannis & P. Dimitriadis, From time series to stochastics 1 Schools on names and definitions The poetic school What’s in a name? That which we call a rose, By any other name would smell as sweet. —William Shakespeare, “Romeo and Juliet” (Act 2, scene 2) The philosophico-epistemic school Ἀρχή σοφίας ὀνομάτων ἐπίσκεψις (The beginning of wisdom is the visiting (inspection) of names) ―Attributed to Antisthenes of Athens, founder of Cynic philosophy Ἀρχὴ παιδεύσεως ἡ τῶν ὀνομάτων ἐπίσκεψις” (The beginning of education is the inspection of names) ―Attributed to Socrates by Epictetus, Discourses, Ι.17,12, The beginning of wisdom is to call things by their proper name. ―Chinese proverb paraphrasing Confucius’s quote “If names be not correct, language is not in accordance with the truth of things.” D. Koutsoyiannis & P. Dimitriadis, From time series to stochastics 2 On names and definitions (contd.) The philosophico-epistemic school (contd.) When I name an object with a word, I thereby assert its existence.” —Andrei Bely, symbolist poet and former mathematics student of Dmitri Egorov, in his essay “The Magic of Words” “Nommer, c’est avoir individu” (to name is to have individuality).
    [Show full text]
  • Biography of N. N. Luzin
    BIOGRAPHY OF N. N. LUZIN http://theor.jinr.ru/~kuzemsky/Luzinbio.html BIOGRAPHY OF N. N. LUZIN (1883-1950) born December 09, 1883, Irkutsk, Russia. died January 28, 1950, Moscow, Russia. Biographic Data of N. N. Luzin: Nikolai Nikolaevich Luzin (also spelled Lusin; Russian: НиколайНиколаевич Лузин) was a Soviet/Russian mathematician known for his work in descriptive set theory and aspects of mathematical analysis with strong connections to point-set topology. He was the co-founder of "Luzitania" (together with professor Dimitrii Egorov), a close group of young Moscow mathematicians of the first half of the 1920s. This group consisted of the higly talented and enthusiastic members which form later the core of the famous Moscow school of mathematics. They adopted his set-theoretic orientation, and went on to apply it in other areas of mathematics. Luzin started studying mathematics in 1901 at Moscow University, where his advisor was professor Dimitrii Egorov (1869-1931). Professor Dimitrii Fedorovihch Egorov was a great scientist and talented teacher; in addition he was a person of very high moral principles. He was a Russian and Soviet mathematician known for significant contributions to the areas of differential geometry and mathematical analysis. Egorov was devoted and openly practicized member of Russian Orthodox Church and active parish worker. This activity was the reason of his conflicts with Soviet authorities after 1917 and finally led him to arrest and exile to Kazan (1930) where he died from heavy cancer. From 1910 to 1914 Luzin studied at Gottingen, where he was influenced by Edmund Landau. He then returned to Moscow and received his Ph.D.
    [Show full text]
  • Yakov Petrovich Terletskii (1912‒1993)
    134 Hiroshi ICHIKAWA Materialist Perestroika of Quantum Dynamics and Soviet Ideology: Yakov Petrovich Terletskii (1912‒1993) Hiroshi ICHIKAWA* Abstract Yakov Petrovich Terletskii (1912‒1993), a professor in theoretical physics at Moscow State University wrote extensively on the “philosophical” issues of modern physics, particularly in the first half of the 1950s. Outside Russia, his name is rarely mentioned. Nevertheless, it is likely that he has been regarded as a “cat’s-paw” of the ideological apparatus of the Marxist state. May we share such a view even now, in spite of the major changes in the historiography of Soviet science we have already seen in the last few decades? Thanks to the progress of the study based on formerly classified documents in various Russian archives, we see that the totalitarian model that had been applied to the understanding of Soviet society for a long time rapidly lost its popularity and was replaced with a more pluralistic view. This study reexamines and reevaluates the thoughts and activities of Yakov Terletskii from the viewpoint of today’s understanding of the historiography of Soviet science. It provides a new analysis, connecting the global quest of physicists for a new approach to quantum dynamics to the local context particular to the Soviet Union. Keywords: Yakov Terletskii, Modern Physics, Quantum Dynamics, the Copenhagen Interpretation, the Soviet Union, Ideology, Dialectical Materialism. Robespierrists! Anti-Robespierrists! For pity’s sake, simply tell us what Robespierre was. [Marc Bloch, The Historian’s Craft. trans. by Peter Putnam (Manchester: Manchester University Press, 1954), p. 140.] 1. Introduction Independently from Ethan Pollock,1 the author of this paper has taken notice of the consequences of an ideologically charged campaign under the guise of a series of academic “discussions” initiated in various fields of science in Postwar Soviet Union, as * Professor at the Graduate School of Integrated Arts and Sciences, Hiroshima University.
    [Show full text]
  • Notes on the Poisson Point Process
    Notes on the Poisson point process Paul Keeler January 5, 2016 This work is licensed under a “CC BY-SA 3.0” license. Abstract The Poisson point process is a type of random object known as a point pro- cess that has been the focus of much study and application. This survey aims to give an accessible but detailed account of the Poisson point process by covering its history, mathematical definitions in a number of settings, and key properties as well detailing various terminology and applications of the process, with rec- ommendations for further reading. 1 Introduction In probability, statistics and related fields, a Poisson point process or a Poisson process or a Poisson point field is a type of random object known as a point pro- cess or point field that consists of randomly positioned points located on some underlying mathematical space [68]. The process has convenient mathematical properties [39], which has led to it being frequently defined in Euclidean space and used as a mathematical model for seemingly random processes in numer- ous disciplines such as astronomy [5], biology [53], ecology [70], geology [16], physics [61], image processing [12], and telecommunications [7][28]. The Poisson point process is often defined on the real line playing an impor- tant role in the field of queueing theory [40] where it is used to model certain random events happening in time such as the arrival of customers at a store or phone calls at an exchange. In the plane, the point process, also known as a spa- tial Poisson process [9], may represent scattered objects such as users in a wireless network [2, 6, 7, 27], particles colliding into a detector, or trees in a forest [68].
    [Show full text]
  • MATH 590: Meshfree Methods Chapter 1 — Part 1: Introduction and a Historical Overview
    MATH 590: Meshfree Methods Chapter 1 — Part 1: Introduction and a Historical Overview Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 [email protected] MATH 590 – Chapter 1 1 Outline 1 Introduction 2 Some Historical Remarks [email protected] MATH 590 – Chapter 1 2 Introduction Outline 1 Introduction 2 Some Historical Remarks [email protected] MATH 590 – Chapter 1 3 Introduction General Meshfree Methods Meshfree Methods have gained much attention in recent years interdisciplinary field many traditional numerical methods (finite differences, finite elements or finite volumes) have trouble with high-dimensional problems meshfree methods can often handle changes in the geometry of the domain of interest (e.g., free surfaces, moving particles and large deformations) better independence from a mesh is a great advantage since mesh generation is one of the most time consuming parts of any mesh-based numerical simulation new generation of numerical tools [email protected] MATH 590 – Chapter 1 4 Introduction General Meshfree Methods Applications Original applications were in geodesy, geophysics, mapping, or meteorology Later, many other application areas numerical solution of PDEs in many engineering applications, computer graphics, optics, artificial intelligence, machine learning or statistical learning (neural networks or SVMs), signal and image processing, sampling theory, statistics (kriging), response surface or surrogate modeling, finance, optimization. [email protected] MATH 590 – Chapter 1 5 Introduction
    [Show full text]
  • Engineering Theory and Mathematics in the Early Development of Information Theory
    1 Engineering Theory and Mathematics in the Early Development of Information Theory Lav R. Varshney School of Electrical and Computer Engineering Cornell University information theory in the American engineering community Abstract— A historical study of the development of information and also touch on the course that information theory took in theory advancements in the first few years after Claude E. Shannon’s seminal work of 1948 is presented. The relationship the American mathematics community. As the early between strictly mathematical advances, information theoretic development of information theory was being made, social evaluations of existing communications systems, and the meanings of the nascent field were just being constructed and development of engineering theory for applications is explored. It incorporated by the two different social groups, mathematical is shown that the contributions of American communications scientists and communications engineering theorists. The engineering theorists are directly tied to the socially constructed meanings constructed had significant influence on the direction meanings of information theory held by members of the group. The role played by these engineering theorists is compared with that future developments took, both in goals and in the role of American mathematical scientists and it is also shown methodology. Mutual interaction between the two social that the early advancement of information is linked to the mutual groups was also highly influenced by these social meanings. interactions between these two social groups. In the half century since its birth, Shannon information theory has paved the way for electronic systems such as compact discs, cellular phones, and modems, all built upon the I. INTRODUCTION developments and social meanings of its first few years.
    [Show full text]
  • English for It-Students Part Ii
    Куліш І.М., Рогульська А.B. ENGLISH FOR IT-STUDENTS PART II МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ ЧЕРКАСЬКИЙ НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ ІМЕНІ БОГДАНА ХМЕЛЬНИЦЬКОГО Куліш І.М., Рогульська А.В. ENGLISH FOR IT- STUDENTS PART II Черкаси - 2015 УДК 811. 111 (075.8) ББК 81. 2Ан – 91 К 90 Рецензенти: І.В. Некоз – кандидат педагогічних наук, доцент кафедри іноземних мов Черкаського національного університету ім. Б. Хмельницького I.Г. Багач . - канд.пед.наук, доц. кафедри романо-германської філології та перекладу Черкаського державного технологічного університету С.Р. Бабушко - канд.філ.наук, доц. кафедри української та іноземних мов Національного університету фізичного виховання і спорту України, м.Київ. Куліш І.М., Рогульська А.В. Англійська мова для студентів ІТ- спеціальностей: Навчальний посібник для студентів ІТ-спеціальностей. – Черкаси: Вид. від ЧНУ імені Богдана Хмельницького, 2015. – 109 с. Навчальний посібник сформований у відповідності до програми з іноземної мови за професійним спрямуванням для студентів 3 та 4 курсу ІТ- спеціальностей. Навчальний матеріал розташовано за тематичним принципом (за розділами). Робота за темою в межах одного розділу організовується на базі декількох основних текстів, зразків усного мовлення та завдань комунікативного характеру, текстів для додаткового читання, системи граматичних вправ репродуктивного та творчого характеру. Посібник вміщує індивідуальні творчі завдання для самостійної роботи студентів та вправи для контролю умінь та навичок студентів. Затверджено до друку на засіданні вченої ради Черкаського Національного університету імені Богдана Хмельницького протокол №3 від 14 грудня 2015 року ISBN 978-966-353-408-4 ЧНУ ім. Б.Хмельницького, 2015 Куліш І.М., Рогульська А.В., 2015 ПЕРЕДМОВА Навчальний посібник призначається для студентів 3 та 4 курсу ІТ- спеціальностей. Мета посібника – формувати у студентів навички читання та розуміння літератури фахової професійної тематики на основі активізації знань загальновживаної лексики та основ граматики.
    [Show full text]
  • Curso: Análisis Real II Clase 3: Aproximación De Funciones
    CONTENIDO REFERENCIAS SUGERIDAS CURSO:ANÁLISIS REAL II CLASE 3: APROXIMACIÓN DE FUNCIONES MEDIBLES POR FUNCIONES CONTINUAS Yamilet Quintana1 1Universidad Simón Bolívar Programa Interuniversitario de Doctorado en Matemáticas Junio, 2020 CONTENIDO REFERENCIAS SUGERIDAS 1 CONTENIDO Aproximación de funciones medibles por funciones continuas p Densidad de Cc(X) en L (X; µ), 1 ≤ p < 1 El teorema de Luzin 2 REFERENCIAS SUGERIDAS CONTENIDO REFERENCIAS SUGERIDAS 1 CONTENIDO Aproximación de funciones medibles por funciones continuas p Densidad de Cc(X) en L (X; µ), 1 ≤ p < 1 El teorema de Luzin 2 REFERENCIAS SUGERIDAS CONTENIDO REFERENCIAS SUGERIDAS APROXIMACIÓN DE FUNCIONES MEDIBLES POR FUNCIONES CONTINUAS Aproximación de funciones medibles por funciones continuas Para finalizar el Tema 1, en esta lección nos ocuparemos de dos resultados sobre aproximaxión de funciones: El teorema de densi- dad de las funciones continuas con soporte compacto en Lp(X; µ), 1 ≤ p < 1, y el teorema de Luzin sobre aproximación de funciones medibles por funciones continuas. Nuestro objetivo será enunciar y demostrar ambos resultados. CONTENIDO REFERENCIAS SUGERIDAS p DENSIDAD DE Cc (X) EN L (X; µ), 1 ≤ p < 1 p Densidad de Cc(X) en L (X; µ), 1 ≤ p < 1 TEOREMA 1-C3 Sean (X; τ) Hausdorff localmente compacto, B(X) su σ-álgebra de Borel, µ : B(X) ! [0; 1] regular, y 1 ≤ p < 1. Entonces Cc (X) es denso en Lp(X; µ). Demostración Seguiremos los enfoques de las referencias [2,10]. p Paso 1: Verifiquemos que Cc (X) es un subconjunto de L (X; µ). En efecto, si f 2 Cc (X) entonces f es acotada y Z p p p kf kLp(X,µ) = jf (x)j dµ(x) ≤ µ(supp(f )) kf k1 : X CONTENIDO REFERENCIAS SUGERIDAS p DENSIDAD DE Cc (X) EN L (X; µ), 1 ≤ p < 1 Demostración (cont.) Paso 2: Sea S el conjunto de las funciones simples que se anulan fuera de conjuntos de medida cero.
    [Show full text]
  • Stochastic Chemical Kinetics
    Stochastic Chemical Kinetics A Special Course on Probability and Stochastic Processes for Physicists, Chemists, and Biologists Summer Term 2011 Version of July 7, 2011 Peter Schuster Institut f¨ur Theoretische Chemie der Universit¨at Wien W¨ahringerstraße 17, A-1090 Wien, Austria Phone: +43 1 4277 527 36 , Fax: +43 1 4277 527 93 E-Mail: pks @ tbi.univie.ac.at Internet: http://www.tbi.univie.ac.at/~pks 2 Peter Schuster Preface The current text contains notes that were prepared first for a course on ‘Stochastic Chemical Kinetics’ held at Vienna University in the winter term 1999/2000 and repeated in the winter term 2006/2007. The current version refers to the summer term 2011 but no claim is made that it is free of errors. The course is addressed to students of chemistry, biochemistry, molecular biology, mathematical and theoretical biology, bioinformatics, and systems biology with particular interests in phenomena observed at small particle numbers. Stochastic kinetics is an interdisciplinary subject and hence the course will contain elements from various disciplines, mainly from probability theory, mathematical statistics, stochastic processes, chemical kinetics, evo- lutionary biology, and computer science. Considerable usage of mathemat- ical language and analytical tools is indispensable, but we have consciously avoided to dwell upon deeper and more formal mathematical topics. This series of lectures will concentrate on principles rather than technical details. At the same time it will be necessary to elaborate tools that allow to treat real problems. Analytical results on stochastic processes are rare and thus it will be unavoidable to deal also with approximation methods and numerical techniques that are able to produce (approximate) results through computer calculations (see, for example, the articles [1–4]).
    [Show full text]
  • Editorial Team European Mathematical Society
    CONTENTS EDITORIAL TEAM EUROPEAN MATHEMATICAL SOCIETY EDITOR-IN-CHIEF ROBIN WILSON Department of Pure Mathematics The Open University Milton Keynes MK7 6AA, UK e-mail: [email protected] ASSOCIATE EDITORS STEEN MARKVORSEN Department of Mathematics Technical University of Denmark Building 303 NEWSLETTER No. 34 DK-2800 Kgs. Lyngby, Denmark e-mail: [email protected] December 1999 KRZYSZTOF CIESIELSKI Mathematics Institute Jagiellonian University Reymonta 4 EMS News : Committee and Agenda .................................................................. 2 30-059 Kraków, Poland Editorial by EMS Secretary David Brannan ..................................................... 3 e-mail: [email protected] KATHLEEN QUINN Introducing the Editorial Team : Part 2 ........................................................... 7 Open University [address as above] EMS Executive meeting in Zurich ....................................................................... 8 e-mail: [email protected] Felix Klein Prize ............................................................................................... 10 SPECIALIST EDITORS 3ecm Barcelona: Second Announcement .......................................................... 12 INTERVIEWS Steen Markvorsen [address as above] EMS Poster Competition .................................................................................. 14 SOCIETIES EMS Council Meeting in Barcelona: Second Announcement ........................... 15 Krzysztof Ciesielski [address as above] Interviews - Jan van Maanen..........................................................................
    [Show full text]
  • Notes on the Poisson Point Process
    Notes on the Poisson point process Paul Keeler March 20, 2018 This work is licensed under a “CC BY-SA 3.0” license. Abstract The Poisson point process is a type of random object in mathematics known as a point process. For over a century this point process has been the focus of much study and application. This survey aims to give an accessible but detailed account of the Poisson point process by covering its history, mathematical defini- tions in a number of settings, and key properties as well detailing various termi- nology and applications of the point process, with recommendations for further reading. 1 Introduction In probability, statistics and related fields, a Poisson point process or a Pois- son process (also called a Poisson random measure or a Poisson point field) is a type of random object known as a point process that consists of randomly po- sitioned points located on some underlying mathematical space [51, 54]. The Poisson point process has interesting mathematical properties [51][54, pages xv- xvi], which has led to it being frequently defined in Euclidean space and used as a mathematical model for seemingly random objects in numerous disciplines such as astronomy [4], biology [62], ecology [84], geology [17], physics [74], im- age processing [10], and telecommunications [33]. For example, defined on the real line, it can be interpreted as a stochastic process 1 called the Poisson pro- cess [83, pages 7–8][2][70, pages 94–95], forming an example of a Levy´ process and a Markov process [73, page ix] [2]. In this setting, the Poisson process plays 1Point processes and stochastic processes are usually considered as different objects, with point processes being associated or related to certain stochastic processes [20, page 194][18, page 3][83, pages 7–8][25, page 580].
    [Show full text]
  • Golden Rotations
    GOLDEN ROTATIONS OLIVER KNILL Abstract. These are expanded preparation notes to a talk given on February 23, 2015 at BU. This was the abstract: "We look Pn k at Birkhoff sums Sn(t)=n = k=1 Xk(t)=n with Xk(t) = g(T t), where T is the irrational golden rotation and where g(t) = cot(πt). Such sums have been studied by number theorists like Hardy and Littlewood or Sinai and Ulcigrain [41] in the context of the curlicue problem. Birkhoff sums can be visualized if the time interval [0; n] is rescaled so that it displays a graph over the interval [0; 1]. While 1 for any L -function g(t) and ergodic T , the sum S[nx](t)=n con- verges almost everywhere to a linear function Mx by Birkhoff's ergodic theorem, there is an interesting phenomenon for Cauchy distributed random variables, where g(t) = cot(πt). The func- tion x ! S[nx]=n on [0; 1] converges for n ! 1 to an explicitly given fractal limiting function, if n is restricted to Fibonacci num- bers F (2n) and if the start point t is 0. The convergence to the \golden graph" shows a truly self-similar random walk. It ex- plains some observations obtained together with John Lesieutre and Folkert Tangerman, where we summed the anti-derivative G of g, which happens to be the Hilbert transform of a piecewise lin- ear periodic function. Recently an observation of [22] was proven by [40]. Birkhoff sums are relevant in KAM contexts, both in an- alytic and smooth situations or in Denjoy-Koksma theory which is a refinement of Birkhoff's ergodic theorem for Diophantine ir- rational rotations.
    [Show full text]