3 Sept Irish Math Society
|
|
- Norman Cummings
- 5 years ago
- Views:
Transcription
1 3 Sept Irish Math Society Calculating the Weather: The Mathematics of Atmospheric Modelling Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin
2 Outline of the Lecture Pre-history of NWP Richardson s Forecast The ENIAC Integrations Interlude Data Assimilation Ensemble Prediction Spherical Grids 2
3 Increase in Forecasting Skill 3
4 Relevant Mathematical Areas Partial Differential Equations Numerical Analysis Linear Algebra Variational Methods Dynamical Systems Geometry of the Sphere Something for everyone! 4
5 Ancient Times 5
6 Galileo Galilei ( ) Galileo formulated the basic law of falling bodies, which he verified by careful measurements. He constructed a telescope, with which he studied lunar craters, and discovered four moons revolving around Jupiter. Galileo is credited with the invention of the Thermometer. Thus began quantitative measurements of the atmosphere. 6
7 Evangelista Torricelli Evangelista Torricelli ( ), a student of Galileo, devised the first accurate barometer. Torricelli s Theorem: v = 2gh Torricelli inventing the barometer 7
8 Newton s Law of Motion The rate of change of momentum of a body is equal to the sum of the forces acting on the body. If F is the total applied force, Newton s Second Law gives a differential equation: dp dt = F. The acceleration a is the rate of change of velocity, that is, a = dv/dt. If the mass m is constant, we have F = ma. Force = Mass Acceleration. 8
9 Euler s Equations for Fluid Flow Leonhard Euler Born in Basel in Died 1783 in St Petersburg. Formulated the equations for incompressible, inviscid fluid flow: V t + V V + 1 ρ p = g. V = 0 Partial differential equations. 9
10 Jean Le Rond d Alembert A body moving at constant speed through a gas or a fluid does not experience any resistance (D Alembert 1752). 10
11 George G Stokes, George Gabriel Stokes, founder of modern hydrodynamics. 11
12 ASIDE: Stokes Theorem Γ V dl = Σ V n da. Stokes Theorem was actually discovered by Kelvin in It is of central importance in fluid dynamics. It leads on to Bjerknes Circulation Theorem: dc dt = 1 ρ p da = dp ρ, Σ which generalized Kelvin s Circulation Theorem to baroclinic fluids (ρ varying independently of p), and ushered in the study of Geophysical Fluid Dynamics. Γ 12
13 Resolution of d Alembert s Paradox The minutest amount of viscosity has a profound qualitative impact on the character of the solution. The Navier-Stokes equations incorporate the effect of viscosity. 13
14 The Navier-Stokes Equations Euler s Equations: V t + V V + 1 ρ p = g. The Navier-Stokes Equations V t + V V + 1 ρ p = ν 2 V + g. Motion on the rotating Earth: V t + V V + 2Ω V + 1 ρ p = ν 2 V + g. 14
15 15
16 The Inventors of Thermodynamics 16
17 The Primitive Equations ( du dt f + u tan φ ) v + 1 p a ρ x + F x = 0 ( dv dt + f + u tan φ ) u + 1 p a ρ y + F y = 0 p = RρT p z + gρ = 0 dt dt + (γ 1)T V = Q c p ρ t + ρv = 0 ρ w t + ρ wv = [Sources Sinks] Seven equations; seven variables (u, v, w, p, T, ρ, ρ w ). 17
18 Scientific Weather Forecasting in a Nut-Shell The atmosphere is a physical system Its behaviour is governed by the laws of physics These laws are expressed quantitatively in the form of mathematical equations Using observations, we can specify the atmospheric state at a given initial time: Today s Weather Using the equations, we can calculate how this state will change over time: Tomorrow s Weather The equations are very complicated (non-linear) and a powerful computer is required to do the calculations The accuracy decreases as the range increases; there is an inherent limit of predictibility. 18
19 Richardson s Forecast 19
20 Lewis Fry Richardson, During WWI, Richardson computed by hand the pressure change at a single point. It took him two years! His forecast was a catastrophic failure: p = 145 hpa in 6 hrs But Richardson s method was scientifically sound. 20
21 Tendency of a Noisy Signal x 21
22 Evolution of surface pressure before and after NNMI. (Williamson and Temperton, 1981) 22
23 Initialization of Richardson s Forecast Richardson s Forecast has been repeated on a computer. The atmospheric observations for 20 May, 1910, were recovered from original sources. ORIGINAL: INITIALIZED: dp s = +145 hpa/6 h dt dp s = 0.9 hpa/6 h dt Observations: The barometer was steady! 23
24 Richardson s Forecast Factory c François Schuiten 64,000 Computers: The first Massively Parallel Processor 24
25 The Finite Difference Scheme Let Q be governed by an equation dq dt = F (Q). The time interval under consideration is sliced into a finite number of discrete time steps {0, t, 2 t,..., n t,... }. The time derivitive is approximated by a finite difference: dq dt Q(t + t) Q(t t) 2 t. Thus, a problem in analysis becomes a problem in algebra. 25
26 Reversing History Differential calculus depends upon justifying the limiting process t 0. In approximating a differential equation, we reverse the procedure, and replace derivatives by ratios of increments. We thus... return to the manner in which they did things before the calculus was invented... (Richardson) 26
27 Stepping Forward The time derivative in dq dt = F (Q). is now approximated by a centered difference Q n+1 Q n 1 2 t = F n, Then Q n+1 = Q n t F n. This process of stepping forward is repeated a large number of times, until the desired forecast range is reached. We can discretize space in a similar way, but... 27
28 The Spectral Method The ECMWF Integrated Forecast System (IFS) uses a spectral representation of the meteorological fields. Each field is expanded in spherical harmonics, truncated at a fixed total wavenumber N: N n Q(λ i, φ j, t) = Q m n (t)yn m (λ i, φ j ) n=0 m= n The functions Y m n (λ, φ) are eigensolutions of the Laplacian: 2 Y m n = n(n + 1)Y m n. The coefficients Q m n (t) depend only on time. When the model equations are transformed to spectral space, they become a set of ordinary differential equations for the spectral coefficients Q m n. 28
29 ENIAC Forecast 29
30 The Meteorology Project Project estblished by John von Neumann in Objective of the project: To study the problem of predicting the weather using a digital electronic computer. A Proposal for Funding listed three possibilities : New methods of weather prediction Rational basis for planning observations Step towards influencing the weather! 30
31 The ENIAC 31
32 The ENIAC The ENIAC was the first multi-purpose programmable electronic digital computer. It had: 18,000 vacuum tubes 70,000 resistors 10,000 capacitors 6,000 switches Power: 140 kwatts 32
33 Charney, et al., Tellus, [ ] Absolute Vorticity = [ ] Relative Vorticity + [ ] Planetary Vorticity η = ζ + f. The atmosphere is treated as a single layer. The flow is assumed to be nondivergent. Absolute vorticity is conserved. d(ζ + f) dt = 0. This equation looks deceptively simple. But it is nonlinear: { } ψ 2 ψ t [ 2 ψ] + ψ 2 ψ + β ψ x y y x x = 0, 33
34 ENIAC Forecast for Jan 5,
35 NWP Operations The Joint Numerical Weather Prediction Unit was established on July 1, 1954: Air Weather Service of US Air Force The US Weather Bureau The Naval Weather Service. Operational numerical weather forecasting began in May, 1955, using a three-level quasi-geostrophic model. 35
36 Interlude 36
37 Temperature, Humidity and Climate Change 37
38 Data Assimilation 38
39 Data Assimilation NWP: An initial/boundary value problem Given an estimate of the present state of the atmosphere (initial conditions) appropriate surface and lateral boundary conditions the model forecasts the evolution of the atmosphere. Operational NWP centers produce initial conditions from a statistical combination of observations and short-range forecasts. This is called data assimilation. 39
40 40
41 41
42 42
43 Optimal Interpolation The analysis problem is to find an optimum atmospheric state, x a, given A background field x b (on a regular grid) A set of (irregularly spaced) p observations y o The analysis is cast as background plus increment: x a = x b + W[y o H(x b )] The analysis and the background are vectors of length n. The weights are given by a matrix W of size (n p). 43
44 The Full Set of OI Equations The result of the (least squares) optimization is: x a = x b + W[y o H(x b )] W = BH T (R + HBH T ) 1 P a = (I WH)B All the covariance matrices are modelled using simplifying assumptions. Solution is a formidable computational task: The matrices are huge. Many shortcuts are needed. 44
45 Variational Assimilation Another approach to objective analysis is the variational assimilation technique. Problem: Find the analysis x that minimizes a cost function: J(x) = 1 2 { } (x x b ) T B 1 (x x b ) + [y o H(x)] T R 1 [y o H(x)] the distance between x and the background x b, plus the distance to the observations y o, : Variational assimilation has been shown to yield significant improvements in forecast accuracy. 45
46 The gradient of J with respect to x is J(x) = [B 1 + H T R 1 H](x x b ) H T R 1 {y o H(x b )} To find a minimum of J, we set J(x) = 0. The result is: x = x b + [ B 1 + H T R 1 H ] 1 H T R 1 {y o H(x b )} This is the (formal) solution of the 3-dimensional variational (3D-Var) analysis problem. The matrices are huge: perhaps
47 Minimization In practical 3D-Var, we do not invert a huge matrix. We find the minimum of J(x) by computing the cost function and using an optimization technique. The idea is to proceed downhill as fast as possible: Steepest Descent algorithm, Newton s method, Conjugate Gradient algorithm. 47
48 4D-Variational Assimilation Four-dimensional variational assimilation (4D-Var) is an extension of 3D-Var to allow for observations distributed within a time interval (t 0, t n ). The cost function includes terms for the distance to observations at the time of the observation. J[x(t 0 )] = 1 2 [x(t 0) x b (t 0 )] T B 1 0 [x(t 0) x b (t 0 )] N i=0 [ H(xi ) y o i ] T R 1 i The control variable is the initial state x (t 0 ). [ H(xi ) y o i ] 48
49 Schematic diagram of four dimensional variational assimilation. 49
50 Tangent Linear Model The solution at time t i+1 is computed from the solution at time t i by a (nonlinear) model: x i+1 = M i [x i ]. If we perturb the initial conditions, the solution is x i+1 + δx i+1 = M i [x i + δx i ] The linear tangent model is the (Jacobian) matrix: Then, to first order, [L i ] j,k = [M(x i)] j (x i ) k δx i+1 = L i δx i. 50
51 The Adjoint Model The transpose of the linear tangent model is called the adjoint model. The Gradient of the cost function is: J N [ ] = L T x 0 LT 1 LT i 1 0 i=0 H T i R 1 i d i Every iteration of the 4D-Var minimization requires the computation of the gradient: Compute the observation increments d i during a forward integration Multiply them by H T i R 1 i Integrate these weighted increments backward to the initial time using the adjoint model. 51
52 Atmospheric Normal Modes 52
53 Oscillations of the Atmosphere We treating the atmosphere as a thin single layer: du uv tan φ fv + g h dt a x = 0 dv dt + fu + u2 tan φ + g h a y = 0 ( dh u dt + h x + v y v tan φ ) = 0 a These are the shallow water quations on the sphere. We linearize about a motionless state with depth H: h(λ, φ, t) = Y (φ) exp[i(mλ σt)] After some algebra we get an equation for Y (φ): [( ) ] { [ ] d 1 µ 2 dy 1 m σ 2 + µ 2 dµ σ 2 µ 2 + dµ σ 2 µ 2 σ σ 2 µ 2 m2 1 µ 2 where µ = sin φ and ɛ = (2Ωa) 2 /gh. + ɛ } Y = 0. 53
54 Laplace Tidal Equation Again, the meridional structure is given by [( ) ] { [ ] d 1 µ 2 dy 1 m σ 2 + µ 2 dµ σ 2 µ 2 + dµ σ 2 µ 2 σ σ 2 µ 2 m2 1 µ 2 where µ = sin φ and ɛ = (2Ωa) 2 /gh. + ɛ } Y = 0. The normal modes are determined by the eigensolutions of this second order o.d.e., the Laplace Tidal Equation. Boundary conditions require Y to be regular at the poles. The Laplace Tidal Equation is not in Sturm-Lioville form. 54
55 Mathematical Difficulties The standard form of the Sturm-Liouville equation is ( d p(µ) dy ) + [q(µ) + λr(µ)]y = 0 dµ dµ where p(µ) is regular and has no zeros within the domain. For the Sturm-Liouville Equation: 1. The equation is self-adjoint and the eigenvalues λ are real. 2. The eigenfunctions for different λ are orthogonal. 3. The eigenfunctions form a complete set. 4. There is a denumerable infinity of non-negative eigenvalues with a single limit point at The zeros of the eigenfunctions behave according to the Sturmian oscillation theorems. 55
56 For the LTE, p(µ) = (1 µ 2 )/(σ 2 µ 2 ) blows up at the critical latutudes where µ = ±σ, and the equation is singular. Since the LTE cannot be written in standard Sturm-Liouville form, the five properties may not hold. It has been shown that the eigenvalues ɛ n of the LTE are real and the eigenfunctions form a complete, orthogonal set. Fourth and fifth properties do not hold. For σ < 1 there is a double infinity of eigenvalues, with limit points at both + and. The zeros of the eigenfunctions do not behave in a simple manner like for a regular Sturm-Liouville problem. 56
57 Eigenfrequencies σ of the LTE. 57
58 Atmospheric Predictability and Ensemble Forecasting 58
59 Progress in numerical weather prediction over the past fifty years has been quite dramatic. Forecast skill continues to increase... by one day per decade. However, there is a limit... 59
60 Chaos in Atmospheric Flow In a paper published in 1963, entitled Deterministic Nonperiodic Flow, Edward Lorenz showed that the solutions of the system Edward Lorenz (b. 1917) ẋ = σx + σy ẏ = xz + rx ż = +xy bz are highly sensitive to the initial conditions. 60
61 Identical Twin Experiment 61
62 Ensemble Forecasting In recognition of the chaotic nature of the atmosphere, focus has now shifted to predicting the probability of alternative weather events rather than a single outcome. The mechanism is the Ensemble Prediction System (EPS) and the world leader in this area is the European Centre for Medium-range Weather Forecasts (ECMWF). 62
63 Ensemble remains compact Ensemble spreads out 63
64 Singular Vectors The linear tangent model L i transforms a perturbation at time t i to a perturbation at time t i+1 : δx(t i+1 ) = L i δx(t i ) Perturbation growth is measured by the norm: δx(t i+1 ) 2 = L i δx(t i ) 2 = Lδx, Lδx = δx, L T Lδx. This depends on the eigenvalues of L T L, the singular values. The singular vector corresponding to the maximum singular value gives the component that grows fastest. 64
65 EPS: Ensemble Prediction System We calculate the 25 largest singular values, and the corresponding 25 singular vectors. Fifty perturbed initial states are constructed by adding and subtracting from the analysis. This gives us fifty-one initial states. Fifty-one forecasts are done, starting from these. 65
66 Ensemble of fifty-one 42-hour forecasts. Valid time: 0600 UTC, 26th December,
67 Discretizing the Sphere 67
68 Regular Latitude-Longitude Grid 68
69 69
70 Conformal Stretched Grid 70
71 71
72 Triangulated Icosahedral Grid 72
73 73
74 74
75 75
76 76
77 77
78 Fibonacci Grid Inspired by Sun-flowers and Pineapples 78
79 The ultimate grid remains elusive. Ultimate depends on the application. 79
80 The End Typesetting Software: TEX, Textures, L A TEX, hyperref, texpower, Adobe Acrobat 4.05 Graphics Software: Adobe Illustrator L A TEX Slide Macro Packages: Wendy McKay, Ross Moore
A Century of Numerical Weather Prediction
A Century of Numerical Weather Prediction Peter Lynch School of Mathematical Sciences University College Dublin Royal Meteorological Society, Edinburgh, 10 October, 2008 Outline Prehistory 1890 1920 ENIAC
More informationThe Emergence of Numerical Weather Prediction:
The Emergence of Numerical Weather Prediction: Fulfilment of a Dream & Realization of a Fantasy Peter Lynch School of Mathematics & Statistics University College Dublin Mathematics of Planet Earth Jamboree
More informationOutline. The Spectral Method (MAPH 40260) Part 4: Barotropic Vorticity Equation. Outline. Background. Rossby-Haurwitz Waves.
Outline The Spectral Method (MAPH 40260) Part 4: Barotropic Vorticity Equation Peter Lynch School of Mathematical Sciences Outline The dynamics of non-divergent flow on a rotating sphere are described
More informationThe Spectral Method (MAPH 40260)
The Spectral Method (MAPH 40260) Part 4: Barotropic Vorticity Equation Peter Lynch School of Mathematical Sciences Outline Background Rossby-Haurwitz Waves Interaction Coefficients Transform Method The
More informationIUGG, Perugia July, Energy Spectra from. Entropy Principles. Peter Lynch & Wim Verkely. University College Dublin. Meteorology & Climate Centre
IUGG, Perugia July, 2007 Energy Spectra from Entropy Principles Peter Lynch & Wim Verkely Meteorology & Climate Centre University College Dublin KNMI, De Bilt, Netherlands Introduction The energy distribution
More informationThe Potential Vorticity Equation
The Potential Vorticity Equation The Potential Vorticity Equation The geopotential tendency equation is [ ( )] 2 + f 2 0 Φ t = f p σ p 0 V g + p [ f 2 0 σ V g ( ) 1 2 Φ + f f 0 ( Φ p ) ] The Potential
More informationM.Sc. in Meteorology. Numerical Weather Prediction Prof Peter Lynch
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationOptimal Interpolation ( 5.4) We now generalize the least squares method to obtain the OI equations for vectors of observations and background fields.
Optimal Interpolation ( 5.4) We now generalize the least squares method to obtain the OI equations for vectors of observations and background fields. Optimal Interpolation ( 5.4) We now generalize the
More information4. DATA ASSIMILATION FUNDAMENTALS
4. DATA ASSIMILATION FUNDAMENTALS... [the atmosphere] "is a chaotic system in which errors introduced into the system can grow with time... As a consequence, data assimilation is a struggle between chaotic
More informationRelationship 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 informationRelationship 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 informationFundamentals 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 informationData assimilation; comparison of 4D-Var and LETKF smoothers
Data assimilation; comparison of 4D-Var and LETKF smoothers Eugenia Kalnay and many friends University of Maryland CSCAMM DAS13 June 2013 Contents First part: Forecasting the weather - we are really getting
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Cehtre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. Text for the Course
More informationGFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability
GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What
More informationChapter 2. The continuous equations
Chapter. The continuous equations Fig. 1.: Schematic of a forecast with slowly varying weather-related variations and superimposed high frequency Lamb waves. Note that even though the forecast of the slow
More informationEnsemble 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 informationQuick Recapitulation of Fluid Mechanics
Quick Recapitulation of Fluid Mechanics Amey Joshi 07-Feb-018 1 Equations of ideal fluids onsider a volume element of a fluid of density ρ. If there are no sources or sinks in, the mass in it will change
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationNewton's Laws of Motion
Newton's Laws of Motion Generally, air masses move the way they do because of the forces that act upon them. Definition: A force is something that pushes or pulls on an object, causing it to speed up,
More informationNewton s Second Law of Motion. Isaac Newton (in 1689) The Universe is Regular and Predictable. The Foundation of Science
Key Concepts and Fundamental Theorems Isaac Newton (in 689) of Atmospheric Science This is the earliest portrait of Newton to survive. The artist was Godfrey Kneller, perhaps the greatest portrait painter
More informationBalance in the Atmosphere: Implications for Numerical Weather Prediction
Balance in the Atmosphere: Implications for Numerical Weather Prediction Peter Lynch School of Mathematical Sciences University College Dublin BIRS Summer School, Banff, 10 15 July, 2011 Outline Introduction
More informationt tendency advection convergence twisting baroclinicity
RELATIVE VORTICITY EQUATION Newton s law in a rotating frame in z-coordinate (frictionless): U + U U = 2Ω U Φ α p U + U U 2 + ( U) U = 2Ω U Φ α p Applying to both sides, and noting ω U and using identities
More informationSymmetry methods in dynamic meteorology
Symmetry methods in dynamic meteorology p. 1/12 Symmetry methods in dynamic meteorology Applications of Computer Algebra 2008 Alexander Bihlo alexander.bihlo@univie.ac.at Department of Meteorology and
More informationEnKF 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 informationWeather Forecasting Models in Met Éireann. Eoin Whelan UCD Seminar 3 rd April 2012
Weather Forecasting Models in Met Éireann Eoin Whelan UCD Seminar 3 rd April 2012 Overview Background HIRLAM Models Local Implementation Verification Development work Background Met Éireann Dept of the
More informationComputational challenges in Numerical Weather Prediction
Computational challenges in Numerical Weather Prediction Mike Cullen Oxford 15 September 2008 Contents This presentation covers the following areas Historical background Current challenges Why does it
More informationLecture 1 ATS 601. Thomas Birner, CSU. ATS 601 Lecture 1
Lecture 1 ATS 601 Thomas Birner, CSU About your Instructor: Thomas Birner Assistant Professor, joined CSU 10/2008 M.Sc. Physics (Condensed Matter Theory), U of Leipzig, Germany Ph.D. Atmospheric Science
More informationAttractor of a Shallow Water Equations Model
Thai Journal of Mathematics Volume 5(2007) Number 2 : 299 307 www.math.science.cmu.ac.th/thaijournal Attractor of a Shallow Water Equations Model S. Sornsanam and D. Sukawat Abstract : In this research,
More informationThe spectral transform method
The spectral transform method by Nils Wedi European Centre for Medium-Range Weather Forecasts wedi@ecmwf.int Advanced Numerical Methods for Earth-System Modelling Slide 1 Advanced Numerical Methods for
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationM.Sc. in Meteorology. Numerical Weather Prediction Prof Peter Lynch
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Cehtre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. Text for the Course
More informationObservation 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 information2.5 Shallow water equations, quasigeostrophic filtering, and filtering of inertia-gravity waves
Chapter. The continuous equations φ=gh Φ=gH φ s =gh s Fig..5: Schematic of the shallow water model, a hydrostatic, incompressible fluid with a rigid bottom h s (x,y), a free surface h(x,y,t), and horizontal
More informationIntroduction to Meteorology and Weather Forecasting
Introduction to Meteorology and Weather Forecasting ENVI1400 : Meteorology and Forecasting : lecture 1 2 040909 ENVI1400 : Meteorology and Forecasting : lecture 1 3 040914 ENVI1400 : Meteorology and Forecasting
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationEdward 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 informationLarge-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods
Large-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods Colin Cotter (Imperial College London) & Sebastian Reich (Universität Potsdam) Outline 1. Hydrostatic and semi-geostrophic
More informationPoint Vortex Dynamics in Two Dimensions
Spring School on Fluid Mechanics and Geophysics of Environmental Hazards 9 April to May, 9 Point Vortex Dynamics in Two Dimensions Ruth Musgrave, Mostafa Moghaddami, Victor Avsarkisov, Ruoqian Wang, Wei
More informationCOMPARISON OF THE INFLUENCES OF INITIAL ERRORS AND MODEL PARAMETER ERRORS ON PREDICTABILITY OF NUMERICAL FORECAST
CHINESE JOURNAL OF GEOPHYSICS Vol.51, No.3, 2008, pp: 718 724 COMPARISON OF THE INFLUENCES OF INITIAL ERRORS AND MODEL PARAMETER ERRORS ON PREDICTABILITY OF NUMERICAL FORECAST DING Rui-Qiang, LI Jian-Ping
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationATM 298, Spring 2013 Lecture 2 The Equa;ons of Fluid Mo;on April 3, Paul A. Ullrich (HH 251)
ATM 298, Spring 2013 Lecture 2 The Equa;ons of Fluid Mo;on April 3, 2013 Paul A. Ullrich (HH 251) paullrich@ucdavis.edu Global Atmospheric Modeling Global atmospheric models were originally constructed
More informationChapter 7: Circulation and Vorticity
Chapter 7: Circulation and Vorticity Circulation C = u ds Integration is performed in a counterclockwise direction C is positive for counterclockwise flow!!! Kelvin s Circulation Theorem The rate of change
More informationGeostrophic and Quasi-Geostrophic Balances
Geostrophic and Quasi-Geostrophic Balances Qiyu Xiao June 19, 2018 1 Introduction Understanding how the atmosphere and ocean behave is important to our everyday lives. Techniques such as weather forecasting
More informationM.Sc. in Meteorology. Physical Meteorology Prof Peter Lynch. Mathematical Computation Laboratory Dept. of Maths. Physics, UCD, Belfield.
M.Sc. in Meteorology Physical Meteorology Prof Peter Lynch Mathematical Computation Laboratory Dept. of Maths. Physics, UCD, Belfield. Climate Change???????????????? Tourists run through a swarm of pink
More informationChapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More informationFundamentals of Atmospheric Modelling
M.Sc. in Computational Science Fundamentals of Atmospheric Modelling Peter Lynch, Met Éireann Mathematical Computation Laboratory (Opp. Room 30) Dept. of Maths. Physics, UCD, Belfield. January April, 2004.
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationOverview of the Numerics of the ECMWF. Atmospheric Forecast Model
Overview of the Numerics of the Atmospheric Forecast Model M. Hortal Seminar 6 Sept 2004 Slide 1 Characteristics of the model Hydrostatic shallow-atmosphere approimation Pressure-based hybrid vertical
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationTESTING GEOMETRIC BRED VECTORS WITH A MESOSCALE SHORT-RANGE ENSEMBLE PREDICTION SYSTEM OVER THE WESTERN MEDITERRANEAN
TESTING GEOMETRIC BRED VECTORS WITH A MESOSCALE SHORT-RANGE ENSEMBLE PREDICTION SYSTEM OVER THE WESTERN MEDITERRANEAN Martín, A. (1, V. Homar (1, L. Fita (1, C. Primo (2, M. A. Rodríguez (2 and J. M. Gutiérrez
More informationF11AE1 1. C = ρν r r. r u z r
F11AE1 1 Question 1 20 Marks) Consider an infinite horizontal pipe with circular cross-section of radius a, whose centre line is aligned along the z-axis; see Figure 1. Assume no-slip boundary conditions
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationLecture notes on assimilation algorithms
Lecture notes on assimilation algorithms Elías alur Hólm European Centre for Medium-Range Weather Forecasts Reading, UK April 18, 28 1 Basic concepts 1.1 The analysis In meteorology and other branches
More informationPhase space, Tangent-Linear and Adjoint Models, Singular Vectors, Lyapunov Vectors and Normal Modes
Phase space, Tangent-Linear and Adjoint Models, Singular Vectors, Lyapunov Vectors and Normal Modes Assume a phase space of dimension N where Autonomous governing equations with initial state: = is a state
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationChapter 5. Shallow Water Equations. 5.1 Derivation of shallow water equations
Chapter 5 Shallow Water Equations So far we have concentrated on the dynamics of small-scale disturbances in the atmosphere and ocean with relatively simple background flows. In these analyses we have
More informationContinuum Mechanics Lecture 5 Ideal fluids
Continuum Mechanics Lecture 5 Ideal fluids Prof. http://www.itp.uzh.ch/~teyssier Outline - Helmholtz decomposition - Divergence and curl theorem - Kelvin s circulation theorem - The vorticity equation
More informationIntroduction to Data Assimilation. Reima Eresmaa Finnish Meteorological Institute
Introduction to Data Assimilation Reima Eresmaa Finnish Meteorological Institute 15 June 2006 Outline 1) The purpose of data assimilation 2) The inputs for data assimilation 3) Analysis methods Theoretical
More informationAn Optimal Control Problem Formulation for. the Atmospheric Large-Scale Wave Dynamics
pplied Mathematical Sciences, Vol. 9, 5, no. 8, 875-884 HIKRI Ltd, www.m-hikari.com http://dx.doi.org/.988/ams.5.448 n Optimal Control Problem Formulation for the tmospheric Large-Scale Wave Dynamics Sergei
More informationInitial ensemble perturbations - basic concepts
Initial ensemble perturbations - basic concepts Linus Magnusson Acknowledgements: Erland Källén, Martin Leutbecher,, Slide 1 Introduction Perturbed forecast Optimal = analysis + Perturbation Perturbations
More informationIn the derivation of Optimal Interpolation, we found the optimal weight matrix W that minimizes the total analysis error variance.
hree-dimensional variational assimilation (3D-Var) In the derivation of Optimal Interpolation, we found the optimal weight matrix W that minimizes the total analysis error variance. Lorenc (1986) showed
More informationIntroduction and Vectors Lecture 1
1 Introduction Introduction and Vectors Lecture 1 This is a course on classical Electromagnetism. It is the foundation for more advanced courses in modern physics. All physics of the modern era, from quantum
More informationKey Concepts and Fundamental Theorems of Atmospheric Science
Key Concepts and Fundamental Theorems of Atmospheric Science Part II Global Motion Systems John A. Dutton Meteo 485 Spring 2004 From simple waves to predictability In two lectures a walk through two centuries
More informationMixed Mimetic Spectral Elements for Geophysical Fluid Dynamics
for Geophysical Fluid Dynamics Dave Lee Los Alamos National Laboratory Outline Connection of finite volumes to differential forms Key ideas of differential forms Differential forms for discrete data Construction
More informationA diagnostic interface for ICON Coping with the challenges of high-resolution climate simulations
DLR.de Chart 1 > Challenges of high-resolution climate simulations > Bastian Kern HPCN Workshop, Göttingen > 10.05.2016 A diagnostic interface for ICON Coping with the challenges of high-resolution climate
More information( ) (9.1.1) Chapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation. 9.1 Geostrophy and scaling.
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationMeasurement of Rotation. Circulation. Example. Lecture 4: Circulation and Vorticity 1/31/2017
Lecture 4: Circulation and Vorticity Measurement of Rotation Circulation Bjerknes Circulation Theorem Vorticity Potential Vorticity Conservation of Potential Vorticity Circulation and vorticity are the
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationProf. Simon Tett, Chair of Earth System Dynamics & Modelling: The University of Edinburgh
SAGES Scottish Alliance for Geoscience, Environment & Society Modelling Climate Change Prof. Simon Tett, Chair of Earth System Dynamics & Modelling: The University of Edinburgh Climate Modelling Climate
More informationThe Impact of Background Error on Incomplete Observations for 4D-Var Data Assimilation with the FSU GSM
The Impact of Background Error on Incomplete Observations for 4D-Var Data Assimilation with the FSU GSM I. Michael Navon 1, Dacian N. Daescu 2, and Zhuo Liu 1 1 School of Computational Science and Information
More informationg (z) = 1 (1 + z/a) = 1
1.4.2 Gravitational Force g is the gravitational force. It always points towards the center of mass, and it is proportional to the inverse square of the distance above the center of mass: g (z) = GM (a
More informationLecture 10a: The Hadley Cell
Lecture 10a: The Hadley Cell Geoff Vallis; notes by Jim Thomas and Geoff J. Stanley June 27 In this short lecture we take a look at the general circulation of the atmosphere, and in particular the Hadley
More informationBackground Error Covariance Modelling
Background Error Covariance Modelling Mike Fisher Slide 1 Outline Diagnosing the Statistics of Background Error using Ensembles of Analyses Modelling the Statistics in Spectral Space - Relaxing constraints
More informationA Penalized 4-D Var data assimilation method for reducing forecast error
A Penalized 4-D Var data assimilation method for reducing forecast error M. J. Hossen Department of Mathematics and Natural Sciences BRAC University 66 Mohakhali, Dhaka-1212, Bangladesh I. M. Navon Department
More informationQuantifying observation error correlations in remotely sensed data
Quantifying observation error correlations in remotely sensed data Conference or Workshop Item Published Version Presentation slides Stewart, L., Cameron, J., Dance, S. L., English, S., Eyre, J. and Nichols,
More informationHADLEY CELL EXPANSION IN TODAY S CLIMATE AND PALEOCLIMATES
HADLEY CELL EXPANSION IN TODAY S CLIMATE AND PALEOCLIMATES Bill Langford University Professor Emeritus Department of Mathematics and Statistics University of Guelph, Canada Presented to the BioM&S Symposium
More informationLaplace Transform Integration (ACM 40520)
Laplace Transform Integration (ACM 40520) Peter Lynch School of Mathematical Sciences Outline Basic Theory Residue Theorem Ordinary Differential Equations General Vector NWP Equation Lagrangian Formulation
More informationBalanced Flow on the Spinning Globe
Balanced Flow on the Spinning Globe Peter Lynch School of Mathematical Sciences University College Dublin 14th EMS Annual Meeting Prague, Czech Republic 7 th October 2014 EMS Silver Medal Awardees 2008:
More informationVariational data assimilation
Background and methods NCEO, Dept. of Meteorology, Univ. of Reading 710 March 2018, Univ. of Reading Bayes' Theorem Bayes' Theorem p(x y) = posterior distribution = p(x) p(y x) p(y) prior distribution
More informationThe Hydrostatic Approximation. - Euler Equations in Spherical Coordinates. - The Approximation and the Equations
OUTLINE: The Hydrostatic Approximation - Euler Equations in Spherical Coordinates - The Approximation and the Equations - Critique of Hydrostatic Approximation Inertial Instability - The Phenomenon - The
More informationNumerical Methods for PDE-Constrained Optimization
Numerical Methods for PDE-Constrained Optimization Richard H. Byrd 1 Frank E. Curtis 2 Jorge Nocedal 2 1 University of Colorado at Boulder 2 Northwestern University Courant Institute of Mathematical Sciences,
More informationρ x + fv f 'w + F x ρ y fu + F y Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ Dv Dt + u2 tanφ + vw a a = 1 p Dw Dt u2 + v 2
Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ + uw Dt a a = 1 p ρ x + fv f 'w + F x Dv Dt + u2 tanφ + vw a a = 1 p ρ y fu + F y Dw Dt u2 + v 2 = 1 p a ρ z g + f 'u + F z Dρ Dt + ρ
More informationHow is balance of a forecast ensemble affected by adaptive and non-adaptive localization schemes?
1 How is balance of a forecast ensemble affected by adaptive and non-adaptive localization schemes? Ross Bannister Stefano Migliorini, Laura Baker, Ali Rudd NCEO, DIAMET, ESA University of Reading 2 Outline
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More information7 EQUATIONS OF MOTION FOR AN INVISCID FLUID
7 EQUATIONS OF MOTION FOR AN INISCID FLUID iscosity is a measure of the thickness of a fluid, and its resistance to shearing motions. Honey is difficult to stir because of its high viscosity, whereas water
More informationForecasting and data assimilation
Supported by the National Science Foundation DMS Forecasting and data assimilation Outline Numerical models Kalman Filter Ensembles Douglas Nychka, Thomas Bengtsson, Chris Snyder Geophysical Statistics
More information3.5 Vorticity Equation
.0 - Marine Hydrodynamics, Spring 005 Lecture 9.0 - Marine Hydrodynamics Lecture 9 Lecture 9 is structured as follows: In paragraph 3.5 we return to the full Navier-Stokes equations (unsteady, viscous
More informationCOMPARISON OF FINITE DIFFERENCE- AND PSEUDOSPECTRAL METHODS FOR CONVECTIVE FLOW OVER A SPHERE
COMPARISON OF FINITE DIFFERENCE- AND PSEUDOSPECTRAL METHODS FOR CONVECTIVE FLOW OVER A SPHERE BENGT FORNBERG and DAVID MERRILL Abstract. For modeling convective flows over a sphere or within a spherical
More informationNew Intro (Chapter 1) (plus slides from Lynch) Eugenia Kalnay AOSC
New Intro (Chapter 1) (plus slides from Lynch) Eugenia Kalnay AOSC 614 2011 M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Cehtre School of Mathematical Sciences
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationDust devils, water spouts, tornados
Balanced flow Things we know Primitive equations are very comprehensive, but there may be a number of vast simplifications that may be relevant (e.g., geostrophic balance). Seems that there are things
More informationADJONT-BASED ANALYSIS OF OBSERVATION IMPACT ON TROPICAL CYCLONE INTENSITY FORECASTS
7A.3 ADJONT-BASED ANALYSIS OF OBSERVATION IMPACT ON TROPICAL CYCLONE INTENSITY FORECASTS Brett T. Hoover* and Chris S. Velden Cooperative Institute for Meteorological Satellite Studies, Space Science and
More informationChapter 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 informationParallel Algorithms for Four-Dimensional Variational Data Assimilation
Parallel Algorithms for Four-Dimensional Variational Data Assimilation Mie Fisher ECMWF October 24, 2011 Mie Fisher (ECMWF) Parallel 4D-Var October 24, 2011 1 / 37 Brief Introduction to 4D-Var Four-Dimensional
More informationQuasi-geostrophic system
Quasi-eostrophic system (or, why we love elliptic equations for QGPV) Charney s QG the motion of lare-scale atmospheric disturbances is overned by Laws of conservation of potential temperature and potential
More informationPredictability is the degree to which a correct prediction or forecast of a system's state can be made either qualitatively or quantitatively.
Predictability is the degree to which a correct prediction or forecast of a system's state can be made either qualitatively or quantitatively. The ability to make a skillful forecast requires both that
More informationRelative Merits of 4D-Var and Ensemble Kalman Filter
Relative Merits of 4D-Var and Ensemble Kalman Filter Andrew Lorenc Met Office, Exeter International summer school on Atmospheric and Oceanic Sciences (ISSAOS) "Atmospheric Data Assimilation". August 29
More information