Data Assimilation Research Testbed Tutorial. Section 7: Some Additional Low-Order Models

Size: px
Start display at page:

Download "Data Assimilation Research Testbed Tutorial. Section 7: Some Additional Low-Order Models"

Transcription

1 Data Assimilation Research Testbed Tutorial Section 7: Some Additional Low-Order Models Version.: September, 6 /home/jla/dart_tutorial/dart/tutorial/section7/tut_section7.fm 1 7/13/7

2 Low-order models in DART: Model Size Features lorenz_63 3 Chaotic, nearly integral attractor, bifurcations lorenz_84 3 More complex attractor, not as periodic 9var 9 Transient off-attractor dynamics lorenz_96 4 (variable) forced_ lorenz_96 lorenz_96 _scale 8 (variable) 44 (variable) lorenz_4 variable Higher dimensional system. Attractor dimension 13. Allows assimilation of model parameter (see Section ). Two primary interacting spatial/temporal scales. Multiscale dynamics. /home/jla/dart_tutorial/dart/tutorial/section7/tut_section7.fm 7/13/7

3 Lorenz 84 model: Attractor not sheet-like. Rare significant deviations. 1 1 Trajectories along deviations don t mesh back up with rest of attractor. 4 x This behavior can be challenging for certain filter variants. /home/jla/dart_tutorial/dart/tutorial/section7/tut_section7.fm 3 7/13/7

4 Lorenz 84 model: 3-variables: 1 dx = x dt x 3 ax 1 + af 1 dx = x dt 1 x bx 1 x 3 x + g dx = bx dt 1 x + x 1 x 3 x 3 4 x Parameters: a=.5, b = 4, f = 8, g = 1.5 can be set from model_nml. /home/jla/dart_tutorial/dart/tutorial/section7/tut_section7.fm 4 7/13/7

5 Lorenz 84 model: Run csh workshop_setup.csh in directory models/lorenz_84/work. Each state variable is observed every once every hour. Observational error variance is 1. Use matlab to examine the output. There s a new type of filter challenge represented here. Can you identify it? Can you propose ways to address it with techniques learned to date? /home/jla/dart_tutorial/dart/tutorial/section7/tut_section7.fm 5 7/13/7

6 9 Variable model:. Three groups of variables.1 Variables 1-3: Divergence Variables 4-6: Vorticity. Variables 7-9: Height In general, divergence is small. Height and pressure similar. Height and pressure have attractor similar to Lorenz_ /home/jla/dart_tutorial/dart/tutorial/section7/tut_section7.fm 6 7/13/7

7 9 Variable model: X i = U j U k + V j V k v a i X i + Y i + a i z i = U Y + Y V X v a Y j k j k i i i Y i ż i = U j ( z k h k ) + ( z j h j )V k g X i K a i z i + F i U i = b j x i + cy i (1) () (3) (4) V i = b k x i cy i (5) X i = a i x i (6) Y i = a i y i (7) X => Divergence, Y=>Vorticity, z=> height Parameters can be adjusted from model_nml. /home/jla/dart_tutorial/dart/tutorial/section7/tut_section7.fm 7 7/13/7

8 9 Variable model: When perturbed off the attractor, mimics gravity waves. Transient, high frequency oscillations dominate divergence variables. Can also appear in height and pressure variables. Run csh workshop_setup.csh in directory models/9var/work. Y1, Y, Y3 (the vorticity variables) are observed once every 6 hours Observational error variance is.4. Use matlab to examine the output. How do different filter kinds interact with gravity waves? /home/jla/dart_tutorial/dart/tutorial/section7/tut_section7.fm 8 7/13/7

9 Lorenz_96 (4-variable) model: One dimensional cyclic domain [., 1.]. Acts something like synoptic scale weather around mid-latitude circle Lorenz_96 state varnum 1 Ensemble Members of./prior_diag.nc True State Ensemble Mean Ensemble Members () model time (4 timesteps) Lorenz_96 state varnum Ensemble Members of./prior_diag.nc 1 True State Ensemble Mean Ensemble Members () The Attractor 1 3 True_State.nc Lorenz_96 true state model time (4 timesteps) Lorenz_96 state varnum 3 Ensemble Members of./prior_diag.nc 1 True State Ensemble Mean Ensemble Members () model time (4 timesteps) state variable # state variable # state variable # 1 /home/jla/dart_tutorial/dart/tutorial/section7/tut_section7.fm 9 7/13/7

10 Lorenz_96 (4-variable) model: Attractor dimension 13 by some measures. Start to explore model sizes closer to ensemble size. Can examine possible degeneracy issues with sample covariance. Naive application of small ensembles diverges in many cases. Run csh workshop_setup.csh in directory models/lorenz_96/work. 4 observations, randomly located in space, equally spaced in time. Observed once an hour; Observational error variance is 1.. Use matlab to examine the output. Need new techniques to fix this. /home/jla/dart_tutorial/dart/tutorial/section7/tut_section7.fm 1 7/13/7

Ensemble Data Assimilation and Uncertainty Quantification

Ensemble Data Assimilation and Uncertainty Quantification Ensemble Data Assimilation and Uncertainty Quantification Jeff Anderson National Center for Atmospheric Research pg 1 What is Data Assimilation? Observations combined with a Model forecast + to produce

More information

Data Assimilation Research Testbed Tutorial

Data Assimilation Research Testbed Tutorial Data Assimilation Research Testbed Tutorial Section 3: Hierarchical Group Filters and Localization Version 2.: September, 26 Anderson: Ensemble Tutorial 9//6 Ways to deal with regression sampling error:

More information

Data Assimilation Research Testbed Tutorial

Data Assimilation Research Testbed Tutorial Data Assimilation Research Testbed Tutorial Section 2: How should observations of a state variable impact an unobserved state variable? Multivariate assimilation. Single observed variable, single unobserved

More information

Analysis of Dynamical Systems

Analysis of Dynamical Systems 2 YFX1520 Nonlinear phase portrait Mathematical Modelling and Nonlinear Dynamics Coursework Assignments 1 ϕ (t) 0-1 -2-6 -4-2 0 2 4 6 ϕ(t) Dmitri Kartofelev, PhD 2018 Variant 1 Part 1: Liénard type equation

More information

DART Tutorial Sec'on 1: Filtering For a One Variable System

DART Tutorial Sec'on 1: Filtering For a One Variable System DART Tutorial Sec'on 1: Filtering For a One Variable System UCAR The Na'onal Center for Atmospheric Research is sponsored by the Na'onal Science Founda'on. Any opinions, findings and conclusions or recommenda'ons

More information

Ensemble forecasting and flow-dependent estimates of initial uncertainty. Martin Leutbecher

Ensemble forecasting and flow-dependent estimates of initial uncertainty. Martin Leutbecher Ensemble forecasting and flow-dependent estimates of initial uncertainty Martin Leutbecher acknowledgements: Roberto Buizza, Lars Isaksen Flow-dependent aspects of data assimilation, ECMWF 11 13 June 2007

More information

Edward Lorenz: Predictability

Edward Lorenz: Predictability Edward Lorenz: Predictability Master Literature Seminar, speaker: Josef Schröttle Edward Lorenz in 1994, Northern Hemisphere, Lorenz Attractor I) Lorenz, E.N.: Deterministic Nonperiodic Flow (JAS, 1963)

More information

Multiple-scale error growth and data assimilation in convection-resolving models

Multiple-scale error growth and data assimilation in convection-resolving models Numerical Modelling, Predictability and Data Assimilation in Weather, Ocean and Climate A Symposium honouring the legacy of Anna Trevisan Bologna, Italy 17-20 October, 2017 Multiple-scale error growth

More information

6.2 Brief review of fundamental concepts about chaotic systems

6.2 Brief review of fundamental concepts about chaotic systems 6.2 Brief review of fundamental concepts about chaotic systems Lorenz (1963) introduced a 3-variable model that is a prototypical example of chaos theory. These equations were derived as a simplification

More information

DART_LAB Tutorial Section 2: How should observations impact an unobserved state variable? Multivariate assimilation.

DART_LAB Tutorial Section 2: How should observations impact an unobserved state variable? Multivariate assimilation. DART_LAB Tutorial Section 2: How should observations impact an unobserved state variable? Multivariate assimilation. UCAR 2014 The National Center for Atmospheric Research is sponsored by the National

More information

Ensemble Kalman Filter

Ensemble Kalman Filter Ensemble Kalman Filter Geir Evensen and Laurent Bertino Hydro Research Centre, Bergen, Norway, Nansen Environmental and Remote Sensing Center, Bergen, Norway The Ensemble Kalman Filter (EnKF) Represents

More information

Maps and differential equations

Maps and differential equations Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number

More information

xt+1 = 1 ax 2 t + y t y t+1 = bx t (1)

xt+1 = 1 ax 2 t + y t y t+1 = bx t (1) Exercise 2.2: Hénon map In Numerical study of quadratic area-preserving mappings (Commun. Math. Phys. 50, 69-77, 1976), the French astronomer Michel Hénon proposed the following map as a model of the Poincaré

More information

Data assimilation in high dimensions

Data assimilation in high dimensions Data assimilation in high dimensions David Kelly Kody Law Andy Majda Andrew Stuart Xin Tong Courant Institute New York University New York NY www.dtbkelly.com February 3, 2016 DPMMS, University of Cambridge

More information

Model error and parameter estimation

Model error and parameter estimation Model error and parameter estimation Chiara Piccolo and Mike Cullen ECMWF Annual Seminar, 11 September 2018 Summary The application of interest is atmospheric data assimilation focus on EDA; A good ensemble

More information

Reducing the Impact of Sampling Errors in Ensemble Filters

Reducing the Impact of Sampling Errors in Ensemble Filters Reducing the Impact of Sampling Errors in Ensemble Filters Jeffrey Anderson NCAR Data Assimilation Research Section The National Center for Atmospheric Research is sponsored by the National Science Foundation.

More information

Edward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology

Edward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology The Lorenz system Edward Lorenz Professor of Meteorology at the Massachusetts Institute of Technology In 1963 derived a three dimensional system in efforts to model long range predictions for the weather

More information

Ensemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors

Ensemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors Ensemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors Eugenia Kalnay Lecture 2 Alghero, May 2008 Elements of Ensemble Forecasting

More information

Chapter 6: Ensemble Forecasting and Atmospheric Predictability. Introduction

Chapter 6: Ensemble Forecasting and Atmospheric Predictability. Introduction Chapter 6: Ensemble Forecasting and Atmospheric Predictability Introduction Deterministic Chaos (what!?) In 1951 Charney indicated that forecast skill would break down, but he attributed it to model errors

More information

A variance limiting Kalman filter for data assimilation: I. Sparse observational grids II. Model error

A variance limiting Kalman filter for data assimilation: I. Sparse observational grids II. Model error A variance limiting Kalman filter for data assimilation: I. Sparse observational grids II. Model error Georg Gottwald, Lewis Mitchell, Sebastian Reich* University of Sydney, *Universität Potsdam Durham,

More information

The Role of Mathematics in Understanding the Earth s Climate. Andrew Roberts

The Role of Mathematics in Understanding the Earth s Climate. Andrew Roberts The Role of Mathematics in Understanding the Earth s Climate Andrew Roberts Outline What is climate (change)? History of mathematics in climate science How do we study the climate? Dynamical systems Large-scale

More information

Standing Swells Surveyed Showing Surprisingly Stable Solutions for the Lorenz 96 Model

Standing Swells Surveyed Showing Surprisingly Stable Solutions for the Lorenz 96 Model nternational ournal of Bifurcation and Chaos, Vol., No. () 7 ( pages) c World Scientific Publishing Company DO:./S877 Standing Swells Surveyed Showing Surprisingly Stable Solutions for the Lorenz 9 Model

More information

Relative importance of the tropics, internal variability, and Arctic amplification on midlatitude climate and weather

Relative importance of the tropics, internal variability, and Arctic amplification on midlatitude climate and weather Relative importance of the tropics, internal variability, and Arctic amplification on midlatitude climate and weather Arun Kumar Climate Prediction Center College Park, Maryland, USA Workshop on Arctic

More information

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. EE 5 Homewors Fall 04 Updated: Sunday, October 6, 04 Some homewor assignments refer to the textboos: Slotine and Li, etc. For full credit, show all wor. Some problems require hand calculations. In those

More information

DART Ini)al Condi)ons for a Refined Grid CAM- SE Forecast of Hurricane Katrina. Kevin Raeder (IMAGe) Colin Zarzycki (ASP)

DART Ini)al Condi)ons for a Refined Grid CAM- SE Forecast of Hurricane Katrina. Kevin Raeder (IMAGe) Colin Zarzycki (ASP) DART Ini)al Condi)ons for a Refined Grid CAM- SE Forecast of Hurricane Katrina Kevin Raeder (IMAGe) Colin Zarzycki (ASP) 1 Mo)va)on Thousands of processors on current supercomputers. - > new CAM dynamical

More information

Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF

Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Eugenia Kalnay and Shu-Chih Yang with Alberto Carrasi, Matteo Corazza and Takemasa Miyoshi 4th EnKF Workshop, April 2010 Relationship

More information

Gaussian Process Approximations of Stochastic Differential Equations

Gaussian Process Approximations of Stochastic Differential Equations Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run

More information

Cross-validation methods for quality control, cloud screening, etc.

Cross-validation methods for quality control, cloud screening, etc. Cross-validation methods for quality control, cloud screening, etc. Olaf Stiller, Deutscher Wetterdienst Are observations consistent Sensitivity functions with the other observations? given the background

More information

A Note on the Particle Filter with Posterior Gaussian Resampling

A Note on the Particle Filter with Posterior Gaussian Resampling Tellus (6), 8A, 46 46 Copyright C Blackwell Munksgaard, 6 Printed in Singapore. All rights reserved TELLUS A Note on the Particle Filter with Posterior Gaussian Resampling By X. XIONG 1,I.M.NAVON 1,2 and

More information

Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF

Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Eugenia Kalnay and Shu-Chih Yang with Alberto Carrasi, Matteo Corazza and Takemasa Miyoshi ECODYC10, Dresden 28 January 2010 Relationship

More information

11 Chaos in Continuous Dynamical Systems.

11 Chaos in Continuous Dynamical Systems. 11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization

More information

The integral test and estimates of sums

The integral test and estimates of sums The integral test Suppose f is a continuous, positive, decreasing function on [, ) and let a n = f (n). Then the series n= a n is convergent if and only if the improper integral f (x)dx is convergent.

More information

Author: Francine Schevenhoven Utrecht University Supervisors: Dr. ir. Frank Selten (KNMI) Prof. dr. ir. Jason Frank (UU)

Author: Francine Schevenhoven Utrecht University Supervisors: Dr. ir. Frank Selten (KNMI) Prof. dr. ir. Jason Frank (UU) Author: Francine Schevenhoven Utrecht University 3697355 Supervisors: Dr. ir. Frank Selten (KNMI) Prof. dr. ir. Jason Frank (UU) June 3, 26 Abstract Historically, weather and climate forecasting has always

More information

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. EE 533 Homeworks Spring 07 Updated: Saturday, April 08, 07 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. Some homework assignments refer to the textbooks: Slotine

More information

CHAOS -SOME BASIC CONCEPTS

CHAOS -SOME BASIC CONCEPTS CHAOS -SOME BASIC CONCEPTS Anders Ekberg INTRODUCTION This report is my exam of the "Chaos-part" of the course STOCHASTIC VIBRATIONS. I m by no means any expert in the area and may well have misunderstood

More information

Seamless Prediction. Hannah Christensen & Judith Berner. Climate and Global Dynamics Division National Center for Atmospheric Research, Boulder, CO

Seamless Prediction. Hannah Christensen & Judith Berner. Climate and Global Dynamics Division National Center for Atmospheric Research, Boulder, CO Seamless Prediction Can we use shorter timescale forecasts to calibrate climate change projections? Hannah Christensen & Judith Berner Climate and Global Dynamics Division National Center for Atmospheric

More information

Physica D. Aggressive shadowing of a low-dimensional model of atmospheric dynamics

Physica D. Aggressive shadowing of a low-dimensional model of atmospheric dynamics Physica D 241 (2012) 637 648 Contents lists available at SciVerse ScienceDirect Physica D journal homepage: www.elsevier.com/locate/physd Aggressive shadowing of a low-dimensional model of atmospheric

More information

Adaptive ensemble Kalman filtering of nonlinear systems

Adaptive ensemble Kalman filtering of nonlinear systems Adaptive ensemble Kalman filtering of nonlinear systems Tyrus Berry George Mason University June 12, 213 : Problem Setup We consider a system of the form: x k+1 = f (x k ) + ω k+1 ω N (, Q) y k+1 = h(x

More information

DART Tutorial Sec'on 5: Comprehensive Filtering Theory: Non-Iden'ty Observa'ons and the Joint Phase Space

DART Tutorial Sec'on 5: Comprehensive Filtering Theory: Non-Iden'ty Observa'ons and the Joint Phase Space DART Tutorial Sec'on 5: Comprehensive Filtering Theory: Non-Iden'ty Observa'ons and the Joint Phase Space UCAR The Na'onal Center for Atmospheric Research is sponsored by the Na'onal Science Founda'on.

More information

Local Ensemble Transform Kalman Filter: An Efficient Scheme for Assimilating Atmospheric Data

Local Ensemble Transform Kalman Filter: An Efficient Scheme for Assimilating Atmospheric Data Local Ensemble Transform Kalman Filter: An Efficient Scheme for Assimilating Atmospheric Data John Harlim and Brian R. Hunt Department of Mathematics and Institute for Physical Science and Technology University

More information

Exploring and extending the limits of weather predictability? Antje Weisheimer

Exploring and extending the limits of weather predictability? Antje Weisheimer Exploring and extending the limits of weather predictability? Antje Weisheimer Arnt Eliassen s legacy for NWP ECMWF is an independent intergovernmental organisation supported by 34 states. ECMWF produces

More information

Observation Impact Assessment for Dynamic. Data-Driven Coupled Chaotic System

Observation Impact Assessment for Dynamic. Data-Driven Coupled Chaotic System Applied Mathematical Sciences, Vol. 10, 2016, no. 45, 2239-2248 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.65170 Observation Impact Assessment for Dynamic Data-Driven Coupled Chaotic

More information

How does 4D-Var handle Nonlinearity and non-gaussianity?

How does 4D-Var handle Nonlinearity and non-gaussianity? How does 4D-Var handle Nonlinearity and non-gaussianity? Mike Fisher Acknowledgements: Christina Tavolato, Elias Holm, Lars Isaksen, Tavolato, Yannick Tremolet Slide 1 Outline of Talk Non-Gaussian pdf

More information

Question: Total. Points:

Question: Total. Points: MATH 308 May 23, 2011 Final Exam Name: ID: Question: 1 2 3 4 5 6 7 8 9 Total Points: 0 20 20 20 20 20 20 20 20 160 Score: There are 9 problems on 9 pages in this exam (not counting the cover sheet). Make

More information

Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems

Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.

More information

Lecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of:

Lecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Lecture 6 Chaos Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Chaos, Attractors and strange attractors Transient chaos Lorenz Equations

More information

Observability, a Problem in Data Assimilation

Observability, a Problem in Data Assimilation Observability, Data Assimilation with the Extended Kalman Filter 1 Observability, a Problem in Data Assimilation Chris Danforth Department of Applied Mathematics and Scientific Computation, UMD March 10,

More information

DART Tutorial Part IV: Other Updates for an Observed Variable

DART Tutorial Part IV: Other Updates for an Observed Variable DART Tutorial Part IV: Other Updates for an Observed Variable UCAR The Na'onal Center for Atmospheric Research is sponsored by the Na'onal Science Founda'on. Any opinions, findings and conclusions or recommenda'ons

More information

Fundamentals of Data Assimilation

Fundamentals of Data Assimilation National Center for Atmospheric Research, Boulder, CO USA GSI Data Assimilation Tutorial - June 28-30, 2010 Acknowledgments and References WRFDA Overview (WRF Tutorial Lectures, H. Huang and D. Barker)

More information

Data Assimilation: Finding the Initial Conditions in Large Dynamical Systems. Eric Kostelich Data Mining Seminar, Feb. 6, 2006

Data Assimilation: Finding the Initial Conditions in Large Dynamical Systems. Eric Kostelich Data Mining Seminar, Feb. 6, 2006 Data Assimilation: Finding the Initial Conditions in Large Dynamical Systems Eric Kostelich Data Mining Seminar, Feb. 6, 2006 kostelich@asu.edu Co-Workers Istvan Szunyogh, Gyorgyi Gyarmati, Ed Ott, Brian

More information

The ECMWF Hybrid 4D-Var and Ensemble of Data Assimilations

The ECMWF Hybrid 4D-Var and Ensemble of Data Assimilations The Hybrid 4D-Var and Ensemble of Data Assimilations Lars Isaksen, Massimo Bonavita and Elias Holm Data Assimilation Section lars.isaksen@ecmwf.int Acknowledgements to: Mike Fisher and Marta Janiskova

More information

Data assimilation. Polyphemus Training Session. June 9, Introduction 2

Data assimilation. Polyphemus Training Session. June 9, Introduction 2 Data assimilation Polyphemus Training Session June 9, 2009 About Purpose: introduction to optimal interpolation, Kalman filtering and variational data assimilation with Polair3D; application to photochemistry

More information

Why are Discrete Maps Sufficient?

Why are Discrete Maps Sufficient? Why are Discrete Maps Sufficient? Why do dynamical systems specialists study maps of the form x n+ 1 = f ( xn), (time is discrete) when much of the world around us evolves continuously, and is thus well

More information

Short tutorial on data assimilation

Short tutorial on data assimilation Mitglied der Helmholtz-Gemeinschaft Short tutorial on data assimilation 23 June 2015 Wolfgang Kurtz & Harrie-Jan Hendricks Franssen Institute of Bio- and Geosciences IBG-3 (Agrosphere), Forschungszentrum

More information

A model for atmospheric circulation

A model for atmospheric circulation Apeiron, Vol. 19, No. 3, July 2012 264 A model for atmospheric circulation B S Lakshmi JNTU College Of Engineering Hyderabad K L Vasundhara Vidya Jyothi Institute of Technology Hyderabad In this paper

More information

Lecture 20: ODE V - Examples in Physics

Lecture 20: ODE V - Examples in Physics Lecture 20: ODE V - Examples in Physics Helmholtz oscillator The system. A particle of mass is moving in a potential field. Set up the equation of motion. (1.1) (1.2) (1.4) (1.5) Fixed points Linear stability

More information

Chaos. Dr. Dylan McNamara people.uncw.edu/mcnamarad

Chaos. Dr. Dylan McNamara people.uncw.edu/mcnamarad Chaos Dr. Dylan McNamara people.uncw.edu/mcnamarad Discovery of chaos Discovered in early 1960 s by Edward N. Lorenz (in a 3-D continuous-time model) Popularized in 1976 by Sir Robert M. May as an example

More information

TWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations

TWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1

More information

PHY411 Lecture notes Part 5

PHY411 Lecture notes Part 5 PHY411 Lecture notes Part 5 Alice Quillen January 27, 2016 Contents 0.1 Introduction.................................... 1 1 Symbolic Dynamics 2 1.1 The Shift map.................................. 3 1.2

More information

Bred vectors: theory and applications in operational forecasting. Eugenia Kalnay Lecture 3 Alghero, May 2008

Bred vectors: theory and applications in operational forecasting. Eugenia Kalnay Lecture 3 Alghero, May 2008 Bred vectors: theory and applications in operational forecasting. Eugenia Kalnay Lecture 3 Alghero, May 2008 ca. 1974 Central theorem of chaos (Lorenz, 1960s): a) Unstable systems have finite predictability

More information

THREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations

THREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations THREE DIMENSIONAL SYSTEMS Lecture 6: The Lorenz Equations 6. The Lorenz (1963) Equations The Lorenz equations were originally derived by Saltzman (1962) as a minimalist model of thermal convection in a

More information

Dynamical Systems: Lecture 1 Naima Hammoud

Dynamical Systems: Lecture 1 Naima Hammoud Dynamical Systems: Lecture 1 Naima Hammoud Feb 21, 2017 What is dynamics? Dynamics is the study of systems that evolve in time What is dynamics? Dynamics is the study of systems that evolve in time a system

More information

Introduction to ensemble forecasting. Eric J. Kostelich

Introduction to ensemble forecasting. Eric J. Kostelich Introduction to ensemble forecasting Eric J. Kostelich SCHOOL OF MATHEMATICS AND STATISTICS MSRI Climate Change Summer School July 21, 2008 Co-workers: Istvan Szunyogh, Brian Hunt, Edward Ott, Eugenia

More information

EnKF Localization Techniques and Balance

EnKF Localization Techniques and Balance EnKF Localization Techniques and Balance Steven Greybush Eugenia Kalnay, Kayo Ide, Takemasa Miyoshi, and Brian Hunt Weather Chaos Meeting September 21, 2009 Data Assimilation Equation Scalar form: x a

More information

Attractor Learning in Synchronized Chaotic Systems in the Presence of Unresolved Scales

Attractor Learning in Synchronized Chaotic Systems in the Presence of Unresolved Scales Final revised submission to Chaos: An Interdisciplinary Journal of Nonlinear Science Attractor Learning in Synchronized Chaotic Systems in the Presence of Unresolved Scales W. Wiegerinck 1, 2, a) 3, b)

More information

Lecture 1: A Preliminary to Nonlinear Dynamics and Chaos

Lecture 1: A Preliminary to Nonlinear Dynamics and Chaos Lecture 1: A Preliminary to Nonlinear Dynamics and Chaos Autonomous Systems A set of coupled autonomous 1st-order ODEs. Here "autonomous" means that the right hand side of the equations does not explicitly

More information

Experimental Characterization of Nonlinear Dynamics from Chua s Circuit

Experimental Characterization of Nonlinear Dynamics from Chua s Circuit Experimental Characterization of Nonlinear Dynamics from Chua s Circuit John Parker*, 1 Majid Sodagar, 1 Patrick Chang, 1 and Edward Coyle 1 School of Physics, Georgia Institute of Technology, Atlanta,

More information

EnKF Review. P.L. Houtekamer 7th EnKF workshop Introduction to the EnKF. Challenges. The ultimate global EnKF algorithm

EnKF Review. P.L. Houtekamer 7th EnKF workshop Introduction to the EnKF. Challenges. The ultimate global EnKF algorithm Overview 1 2 3 Review of the Ensemble Kalman Filter for Atmospheric Data Assimilation 6th EnKF Purpose EnKF equations localization After the 6th EnKF (2014), I decided with Prof. Zhang to summarize progress

More information

The hybrid ETKF- Variational data assimilation scheme in HIRLAM

The hybrid ETKF- Variational data assimilation scheme in HIRLAM The hybrid ETKF- Variational data assimilation scheme in HIRLAM (current status, problems and further developments) The Hungarian Meteorological Service, Budapest, 24.01.2011 Nils Gustafsson, Jelena Bojarova

More information

What is a Low Order Model?

What is a Low Order Model? What is a Low Order Model? t Ψ = NL(Ψ ), where NL is a nonlinear operator (quadratic nonlinearity) N Ψ (x,y,z,...,t)= Ai (t)φ i (x,y,z,...) i=-n da i = N N cijk A j A k + bij A j + f i v i j;k=-n j=-n

More information

Probability and climate change UK/Russia Climate Workshop, October, 2003

Probability and climate change UK/Russia Climate Workshop, October, 2003 Probability and climate change UK/Russia Climate Workshop, October, 2003 Myles Allen Department of Physics University of Oxford myles.allen@physics.ox.ac.uk South Oxford on January 5 th, 2003 Photo courtesy

More information

... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré

... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré Chapter 2 Dynamical Systems... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré One of the exciting new fields to arise out

More information

Nonlinear atmospheric response to Arctic sea-ice loss under different sea ice scenarios

Nonlinear atmospheric response to Arctic sea-ice loss under different sea ice scenarios Nonlinear atmospheric response to Arctic sea-ice loss under different sea ice scenarios Hans Chen, Fuqing Zhang and Richard Alley Advanced Data Assimilation and Predictability Techniques The Pennsylvania

More information

S56 (5.1) Integration.notebook March 09, 2017

S56 (5.1) Integration.notebook March 09, 2017 Today we will be learning about integration (indefinite integrals) Integration What would you get if you undo the differentiation? Integration is the reverse process of differentiation. It is sometimes

More information

Comparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales

Comparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales Comparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales Meng Zhang and Fuqing Zhang Penn State University Xiang-Yu Huang and Xin Zhang NCAR 4 th EnDA Workshop, Albany, NY

More information

Local Ensemble Transform Kalman Filter

Local Ensemble Transform Kalman Filter Local Ensemble Transform Kalman Filter Brian Hunt 11 June 2013 Review of Notation Forecast model: a known function M on a vector space of model states. Truth: an unknown sequence {x n } of model states

More information

Ensemble Data Assimila.on and Uncertainty Quan.fica.on

Ensemble Data Assimila.on and Uncertainty Quan.fica.on Ensemble Data Assimila.on and Uncertainty Quan.fica.on Jeffrey Anderson, Alicia Karspeck, Tim Hoar, Nancy Collins, Kevin Raeder, Steve Yeager Na.onal Center for Atmospheric Research Ocean Sciences Mee.ng

More information

Supplementary Information Dynamical proxies of North Atlantic predictability and extremes

Supplementary Information Dynamical proxies of North Atlantic predictability and extremes Supplementary Information Dynamical proxies of North Atlantic predictability and extremes Davide Faranda, Gabriele Messori 2 & Pascal Yiou Laboratoire des Sciences du Climat et de l Environnement, LSCE/IPSL,

More information

Representation of model error in a convective-scale ensemble

Representation of model error in a convective-scale ensemble Representation of model error in a convective-scale ensemble Ross Bannister^*, Stefano Migliorini^*, Laura Baker*, Ali Rudd* ^ National Centre for Earth Observation * DIAMET, Dept of Meteorology, University

More information

DART Tutorial Section 18: Lost in Phase Space: The Challenge of Not Knowing the Truth.

DART Tutorial Section 18: Lost in Phase Space: The Challenge of Not Knowing the Truth. DART Tutorial Section 18: Lost in Phase Space: The Challenge of Not Knowing the Truth. UCAR The National Center for Atmospheric Research is sponsored by the National Science Foundation. Any opinions, findings

More information

Mesoscale Predictability of Terrain Induced Flows

Mesoscale Predictability of Terrain Induced Flows Mesoscale Predictability of Terrain Induced Flows Dale R. Durran University of Washington Dept. of Atmospheric Sciences Box 3516 Seattle, WA 98195 phone: (206) 543-74 fax: (206) 543-0308 email: durrand@atmos.washington.edu

More information

APPM 2460 CHAOTIC DYNAMICS

APPM 2460 CHAOTIC DYNAMICS APPM 2460 CHAOTIC DYNAMICS 1. Introduction Today we ll continue our exploration of dynamical systems, focusing in particular upon systems who exhibit a type of strange behavior known as chaos. We will

More information

Nonlinear baroclinic dynamics of surface cyclones crossing a zonal jet

Nonlinear baroclinic dynamics of surface cyclones crossing a zonal jet Nonlinear baroclinic dynamics of surface cyclones crossing a zonal jet Jean-Baptiste GILET, Matthieu Plu and Gwendal Rivière CNRM/GAME (Météo-France, CNRS) 3rd THORPEX International Science Symposium Monterey,

More information

Use of the breeding technique to estimate the structure of the analysis errors of the day

Use of the breeding technique to estimate the structure of the analysis errors of the day Nonlinear Processes in Geophysics (2003) 10: 1 11 Nonlinear Processes in Geophysics c European Geosciences Union 2003 Use of the breeding technique to estimate the structure of the analysis errors of the

More information

Lagrangian data assimilation for point vortex systems

Lagrangian data assimilation for point vortex systems JOT J OURNAL OF TURBULENCE http://jot.iop.org/ Lagrangian data assimilation for point vortex systems Kayo Ide 1, Leonid Kuznetsov 2 and Christopher KRTJones 2 1 Department of Atmospheric Sciences and Institute

More information

Brian J. Etherton University of North Carolina

Brian J. Etherton University of North Carolina Brian J. Etherton University of North Carolina The next 90 minutes of your life Data Assimilation Introit Different methodologies Barnes Analysis in IDV NWP Error Sources 1. Intrinsic Predictability Limitations

More information

Least squares shadowing sensitivity analysis of a modified Kuramoto-Sivashinsky equation

Least squares shadowing sensitivity analysis of a modified Kuramoto-Sivashinsky equation Least squares shadowing sensitivity analysis of a modified Kuramoto-Sivashinsky equation The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.

More information

Igor A. Khovanov,Vadim S. Anishchenko, Astrakhanskaya str. 83, Saratov, Russia

Igor A. Khovanov,Vadim S. Anishchenko, Astrakhanskaya str. 83, Saratov, Russia Noise Induced Escape from Dierent Types of Chaotic Attractor Igor A. Khovanov,Vadim S. Anishchenko, Dmitri G. Luchinsky y,andpeter V.E. McClintock y Department of Physics, Saratov State University, Astrakhanskaya

More information

Data assimilation in high dimensions

Data assimilation in high dimensions Data assimilation in high dimensions David Kelly Courant Institute New York University New York NY www.dtbkelly.com February 12, 2015 Graduate seminar, CIMS David Kelly (CIMS) Data assimilation February

More information

The Application of Contraction Theory in Synchronization of Coupled Chen Systems

The Application of Contraction Theory in Synchronization of Coupled Chen Systems ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.9(2010) No.1,pp.72-77 The Application of Contraction Theory in Synchronization of Coupled Chen Systems Hongxing

More information

State and Parameter Estimation in Stochastic Dynamical Models

State and Parameter Estimation in Stochastic Dynamical Models State and Parameter Estimation in Stochastic Dynamical Models Timothy DelSole George Mason University, Fairfax, Va and Center for Ocean-Land-Atmosphere Studies, Calverton, MD June 21, 2011 1 1 collaboration

More information

THE IMPACT OF SATELLITE-DERIVED WINDS ON GFDL HURRICANE MODEL FORECASTS

THE IMPACT OF SATELLITE-DERIVED WINDS ON GFDL HURRICANE MODEL FORECASTS THE IMPACT OF SATELLITE-DERIVED WINDS ON GFDL HURRICANE MODEL FORECASTS Brian J. Soden 1 and Christopher S. Velden 2 1) Geophysical Fluid Dynamics Laboratory National Oceanic and Atmospheric Administration

More information

4.5 Integration of Rational Functions by Partial Fractions

4.5 Integration of Rational Functions by Partial Fractions 4.5 Integration of Rational Functions by Partial Fractions From algebra, we learned how to find common denominators so we can do something like this, 2 x + 1 + 3 x 3 = 2(x 3) (x + 1)(x 3) + 3(x + 1) (x

More information

HOW TO IMPROVE FORECASTING WITH BRED VECTORS BY USING THE GEOMETRIC NORM

HOW TO IMPROVE FORECASTING WITH BRED VECTORS BY USING THE GEOMETRIC NORM HOW TO IMPROVE FORECASTING WITH BRED VECTORS BY USING THE GEOMETRIC NORM Diego Pazó Instituto de Física de Cantabria (IFCA), CSIC-UC, Santander (Spain) Santander OUTLINE 1) Lorenz '96 model. 2) Definition

More information

Four-Dimensional Ensemble Kalman Filtering

Four-Dimensional Ensemble Kalman Filtering Four-Dimensional Ensemble Kalman Filtering B.R. Hunt, E. Kalnay, E.J. Kostelich, E. Ott, D.J. Patil, T. Sauer, I. Szunyogh, J.A. Yorke, A.V. Zimin University of Maryland, College Park, MD 20742, USA Ensemble

More information

Practical Aspects of Ensemble-based Kalman Filters

Practical Aspects of Ensemble-based Kalman Filters Practical Aspects of Ensemble-based Kalman Filters Lars Nerger Alfred Wegener Institute for Polar and Marine Research Bremerhaven, Germany and Bremen Supercomputing Competence Center BremHLR Bremen, Germany

More information

Comparison of ensemble and NMC type of background error statistics for the ALADIN/HU model

Comparison of ensemble and NMC type of background error statistics for the ALADIN/HU model Comparison of ensemble and NMC type of background error statistics for the ALADIN/HU model Kristian Horvath horvath@cirus.dhz.hr Croatian Meteorological and Hydrological Service supervised by Bölöni Gergely

More information

Data assimilation with Lorenz 3-variable model. Prepared by Shu-Chih Yang Modified by Juan Ruiz.

Data assimilation with Lorenz 3-variable model. Prepared by Shu-Chih Yang Modified by Juan Ruiz. Data assimilation with Lorenz 3-variable model Prepared by Shu-Chih Yang Modified by Juan Ruiz. Governing equations dx dt dy dt dz dt = σ(y x) = rx y xz = xy bz Lorenz, E. N, 1963: Deterministic nonperiodic

More information

Updated: January 16, 2016 Calculus II 7.4. Math 230. Calculus II. Brian Veitch Fall 2015 Northern Illinois University

Updated: January 16, 2016 Calculus II 7.4. Math 230. Calculus II. Brian Veitch Fall 2015 Northern Illinois University Math 30 Calculus II Brian Veitch Fall 015 Northern Illinois University Integration of Rational Functions by Partial Fractions From algebra, we learned how to find common denominators so we can do something

More information

Controlling a Novel Chaotic Attractor using Linear Feedback

Controlling a Novel Chaotic Attractor using Linear Feedback ISSN 746-7659, England, UK Journal of Information and Computing Science Vol 5, No,, pp 7-4 Controlling a Novel Chaotic Attractor using Linear Feedback Lin Pan,, Daoyun Xu 3, and Wuneng Zhou College of

More information