Miscellaneous Thoughts on Ocean Data Assimilation
|
|
- Sherilyn Sylvia Stevens
- 5 years ago
- Views:
Transcription
1 Miscellaneous Thoughts on Ocean Data Assimilation Mike Dowd Dept of Mathematics & Statistics (and Dept of Oceanography) Dalhousie University, Halifax, Canada STATMOS Summer School in Data Assimila5on, Boulder, Colorado, May
2 DA for Oceanography (Physical and Biological) What are the problems? What are the key features? What are some Data Assimilation (DA) issues and approaches used or being considered (and that we have not really covered here)?
3 The State Space Framework for DA x t = d(x t 1,θ) + e t x t ~ [x t,θ x ] or t 1 y t = h(x t,φ) + v t y t x t ~ [y t,φ x t ] dynamics observa)ons t = 1,,T General Features: - A spa)o- temporal problem. High dimensional state space. Actual degrees of freedom? - State vector indexes space as defined on model grid, but could switch space index with )me index and run sequen)ally (later ). - Dynamics relies on numerical ocean models (d operator). Markovian - > sequen)al - Observa)ons indirectly related to the ocean state (h operator), are ojen themselves models (with complex error distribu)ons and sampling protocols) Goal is to es5mate the ocean state x t and the model parameters θ.
4 Ocean Physics/Circulation: Governing Equations Solved numerically via community ocean model code (finite difference, finite element, finite volume) Applied to global, basin-scale, regional, and coastal domains Grid size are O(1-10km), time steps are O(minute)
5 Global Ocean Circulation Model: SST snapshot
6 Coastal Circulation Model
7 Physical Oceanographic DA Features: Nonlinear. Mul)ple interac)ng scales, complex geometry, and boundary condi)ons. Goals Online state es)ma)on for opera)onal oceanogrpahy, i.e. the filtering/predic)on problem. BUT also many scien)fically mo)vated retrospec)ve studies, i.e. reconstruc)on of the past ocean state. i.e. the smoothing problem. Sampling or observing array design of emerging interest (OSSEs)
8 Biological Oceanography: Biogeochemical Models Systems of ODEs: (complex dynamics for biology) dx dt = f (X,θ) + e(t) O- D model Systems of PDEs: (biology within ocean circulation) (complex dynamics and high dimensionality) 1- D or 3- D models X i t + u! X i (K X i ) = f i (X 1,..., X m,θ) + e(t), i = 1,...,m
9 Biogeochemical Dynamics dp dt = g { N;k N }γ P λ P P g P;k P dz dt = εg P;k P dn dt dd dt { } IZ λ Z Z + n Z { } IZ + n P = φd + βg { P;k P } IZ g { N;k N }γ P + νλ Z Z + n N = φd + λ P P + (1 ε β)g{ P;k P } IZ + (1 ν)λ Z Z + n D
10 Snapshot of Ocean Chlorophyll (Phytoplankton)
11 Biological Oceanographic DA Features: Built on physical circula)on models as interac)ng non- conserva)ve tracers. Uncertainty in the biological model structure, and appropriate level of complexity. Primary Goals Parameter es)ma)on, state es)ma)on, as well as system iden)fica)on to aid in scien)fic understanding. Some interest in opera)onal systems è Emerging area, and hence a variety of approaches being considered with respect to DA and dynamical models used. TRADITIONALLY, DEVELOPMENT OF OCEAN DA HAS BEEN DRIVEN MORE BY THE DYNAMICAL MODELS, AND LESS BY THE OBSERVATIONS
12 Measurements of the Physical Ocean
13 Measurements of Ocean Biology 5me series (moored instruments - plankton, nutrients) spa5al data (satellite chorophyll) 5me- space series (glider robots - plankton) OBSERVATION REVOLUTION IN OCEAN SCIENCES, STILL A PARTIALLY OBSERVED SYSTEM AND UNDERSAMPLED COMPLEX SPATIO- TEMPORAL SAMPLING PROTOCOLS
14 Variational DA Approaches Strong Constraint (assumes no model error ): minimize J(θ) = T 2 y t h(x t ) 1 Σv w.r.t θ t=0 subject to x t = d(x t 1,θ) Sta5s5cally, nonlinear (generalized least- squares) regression, with observa)on errors assumed Gaussian. NOTE: For DA, parameters are ojen ini)al or boundary condi)ons. Covariance matrices contain dynamical balances (e.g. geostrophy) The bread and bu[er of prac)cal data assimila)on (3DVAR, 4DVAR). For opera)onal DA, embedded with sequen)al DA cycle. Key aspect is Compu)ng minimiza)on θ J using adjoint equa)ons to faciliate
15 Given the cost function: 4DVAR and Adjoints Consider the prototypical DA problem of es)ma)ng the ini)al condi)ons x 0 J(x 0 ) = 1 2 (x 0 x 0 p ) T Σ 0 1 (x 0 x 0 p ) T (y t Hx t ) T Σ 1 v ( y t Hx t ) t=1 subject to x t = Dx t 1 We define the Lagrange function as: L = J + T t=1 λ T t 1 (x t Dx t 1 ) L λ t = 0 x t = Dx t 1 + e t L = 0 λ t 1 = D T λ t + H T Σ 1 v (y t Hx t ) x t L = 0 x 0 = x p 0 + Σ 0 Dλ 0 x 0 Forward Model Adjoint Model (backward in )me) Parameter es)ma)on
16 A harder varia)onal DA problem.. incorpora)ng dynamical model error Weak Constraint: Sta5s5cally: min J(θ) = T t=1 2 y t h(x t ) 1 Σv wrt x t and θ + x t d(x t 1,θ) Σe 1 2 If linear and Gaussian, random- beta regression (Duncan and Horne 1972, JASA), and the same as the Kalman smoother. Nonlinear case is a much more complex problem can use representers (i.e. solve in observa)on space), or op)miza)on based on calculus of varia)ons (i.e. inverse problems).
17 PROS Remarks on Variational DA Varia)onal DA works for large scale GFD systems. Flexible. Allows for joint es)ma)on for parameters and state Dynamical balances readily incorporated Can use approximate gradients - > incremental approach CONS Loss/Cost func)ons are user defined (usually least squares based, or otherwise assume MV normality) Li[le emphasis on error es)mates Uniqueness of solu)on? ISSUES: ill- posedness, iden)fiability, need for regulariza)on (smoothness priors) Deriva)on of tangent linear and adjoint equa)ons is model- specific, and maintenance of adjoint equa)ons for complex code is non- trivial.
18 Sequential Methods for DA Single stage transition of system from time t-1 to time t prediction measurement x t = d(x t 1,θ) +ν t y t nowcast [x t 1,θ y 1:t 1 ] forecast nowcast [x t,θ y 1:t 1 ] [x t,θ y 1:t ] time = t-1 time = t time - Recursive estimation of system state through time - Forecast and Measurement steps
19 Nudging ˆx t t = ˆx t t 1 + K(y t Hˆx t t 1 ) Many ad- hoc methods used for nudging (= adjus)ng forecast toward measurements). Major dynamical Issues! Now used mainly to correct for model drij from climatology. For linear Gaussian state space model, Kalman filter provides op)mal K (and state covariance propaga)on) - Nonlinear extensions available, e.g. extended Kalman filter. - NonGaussian extensions available, e.g. Gaussian sum filter. Problems with nonlinear dynamics (e.g. extended KF has instabili)es in covariance propaga)on, Evensen s PhD) Ensemble Kalman filter has supplanted these methods.
20 The Ensemble Kalman Filter!x (i) (i) t t = x t t 1 + K(y (i) (i) t Hx t t 1 ), i = 1,,n Works for large systems. Easy to implement Gaussian Approxima)on uses KF upda)ng equa)on NEEDS (spa)al) localiza)on due to spurious long range correla)ons (e.g. the LETKF, or parametric approxima)on/distance based filtering of covariance matrices) Some problems with highly nonlinear dynamics, nonlinear observa)on operators, plus dynamical balances, mass conserva)on
21 Ensemble Mean PON A Comparison Ensemble Std Dev PON depth Stochastic Simulation Stochastic Simulation depth Ens Kalman Filter Ens Kalman Filter depth Particle Filter Particle Filter time Both PF and enkf underperformed but for different reasons
22 Ensembles approaches rely on stochashc dynamic ensemble predichon Key Issue: small ensemble must cover the part of large state space
23 How to do Stochastic Simulation? 1. Dynamical noise: disturbance terms, or random forcing, added to each of the state variables (act as sources/sinks) 2. Stochastic parameters: Model parameters can follow a probability distribution or a stochastic process. 3. Initial and Boundary Conditions: Random initial states and/ or boundary forcing may be incorporated (NWP ensemble forecasting). Like treating these as parameters. 4. Multimodel Ensembles: incorporates structural model errors via using multiple models each with multiple configurations (used in climate modelling) Treatment of Model error is a neglected area?
24 Example: Stochastic Parameters Realization of O-D ODE based PZND model max growth parameter follows AR(1) Frequent transitioning across bifurcation point aperiodic/episodic dynamical dependencies maintained Represent equatorial Pacific
25 We want to generate large large ensembles, GFD models are computahonally expensive What about cheap evaluators for the dynamics, i.e. model surrogates?
26 Idea: The Transition Density and Copulas We want: x t ~ [x t x t 1,θ] - predic)ve/transi)on density We have: x t = d(x t 1,θ) + e t - a numerical model to generate samples Rather than a sample- based representa5on of model errors, what about a parametric representa5on? Must be: accurate, flexible, easy to sample from, high dimensions Idea: create mul)variate distribu)on for [x t x t 1,θ] using copulas - control the marginal distribu)ons, and dependence separately - construc)on and sampling from copulas is straighnorward - can be derived from ensembles and their predic)on errors using method of moments.
27 5 20 De Det Chl Residual Residual Residual 1.0 P(x1,, xn ) = C(P(x1 ),, P(xn )) Oxy Residual 0.15 Based on 1D PZND model of BATS site DIN 5 Depth= 250m Density: Oxy 6.5. SPECIFICATION OF THE STATE-SPACE MODEL Chl Phy Det DIN Depth= 150m Phy Density: MARGINALS Residual e 6.15: kernel density estimates for residuals at -50, -150 and -250 metres deep. read and magnitude of the residuals is quite different at different layers and for t state variables, requiring all 1, 750 marginal distributions to be parametrised ely. Also, nearby residuals are highly correlated. Finally, the residuals are not Samples from predic)ve density, or model errors For example the residuals imply a negative model bias for deep DET and shallow symmetrical about zero implying that the h(xt xt 1 ) is biased for some states and ates and a positive model bias for deep PHY. It is possible to centre h t (xt xt 1 ) at 1 ) b to account for the bias b, but this was not done in this application because CORRELATION mperfection of the data available for estimating b. Figure 6.16: Contour plot for the empirical residual correlation.
28 Idea: Emulator Based DA IDEA: Replace complex dynamical model with an emulator to facilitate computationally efficient data assimilation. Estimate seasonal cycle of two key biological parameters (C:chl ratio, Zooplankton grazing). APPROACH Approximate 3-D biological ocean model with a polynomial chaos based emulator Estimate time-varying parameters by minimizing a discrepancy (metric: adaptive grey block distance) between SeaWiFS satellite chlorophyll images and its model predicted field.
29 Most of our sequenhal eshmahon procedures come from Hme series. AdaptaHons for spaho- temporal problems?
30 Idea: A Location Particle Smoother for Spatio-Temporal Systems IDEA: run par)cle filter through )me, but apply par)cle smoother in the spa>al domain at every )me step. - At each )me, the )me- domain predic)ve distribu)on is used as proposal distribu)on, followed by sequen)al importance resampling with loca)on filter/smoother sweeping through space. Successful for mul)variate state es)ma)on for 1- D biological model (Bermuda Atlan)c Time Series site) Surface Depth Oxygen OXY LPS Deep year one year two
31 Performance of LocaHon ParHcle Smoother Prediction PDF Filter PDF: LPS Simula)on study: Lorenz96 system 10 true states generated 10 replicates per state Cauchy observa)on error (plus different error types, mis- specifica)ons wrt variance and tail density) Frequency σ 2 =.08 σ 2 =.015 Frequency proportion locations covered proportion locations covered Filter PDF: Particle Filter Filter PDF: EnKF Performance Metrics: Percent coverage of true state by 95% credible intervals Average ensemble variance (σ 2 ) Frequency σ 2 è collapse σ 2 =.025 Frequency proportion locations covered proportion locations covered PROS: outperforms EnKF for non- Gaussian likelihood CONS: less computa)onally efficient, no be[er for Gaussian errors
32 What about the (stahc) parameter eshmahon?
33 Approaches for Parameter Estimation 1. State Augmentation: append parameters to the state x t = x t θ t Can use standard sequential MC methods with iterative filtering 2. Likelihood Methods: Use sample based likelihoods T [y 1:T θ] = L(θ y 1:T ) = [y t x t,θ][x t y 1:t 1,θ]dx t t=1 3. Bayesian Hierarchical: particle MCMC, SMC^2
34 Sample Based Likelihood Surface L(θ y 1:T ) 1 n T t=1 ( n [y t x (i),θ] ) i=1 t t 1
35 Concluding Remarks Ocean DA has used a variety of approaches. Successful methods require consideration of problem specific features. Plug and play is unlikely to work (though modularity of methods are important). Physical DA and Biological DA have different goals (i.e. it is not all just filtering). Challenge is to incorporate full suite of observations, large-scale nonlinear models, characterize uncertainty, allow for sampling design, model building/ selection, and multi-model inference. Emerging interest in (hierarchical) Bayesian framework for estimation in dynamical systems. Useful prior information exists, and counters data paucity. But computational issues remain due to expensive GFD models (can cheap, approximate evaluators help?)
Introduction to Particle Filters for Data Assimilation
Introduction to Particle Filters for Data Assimilation Mike Dowd Dept of Mathematics & Statistics (and Dept of Oceanography Dalhousie University, Halifax, Canada STATMOS Summer School in Data Assimila5on,
More informationFundamentals of Data Assimila1on
014 GSI Community Tutorial NCAR Foothills Campus, Boulder, CO July 14-16, 014 Fundamentals of Data Assimila1on Milija Zupanski Cooperative Institute for Research in the Atmosphere Colorado State University
More informationPar$cle Filters Part I: Theory. Peter Jan van Leeuwen Data- Assimila$on Research Centre DARC University of Reading
Par$cle Filters Part I: Theory Peter Jan van Leeuwen Data- Assimila$on Research Centre DARC University of Reading Reading July 2013 Why Data Assimila$on Predic$on Model improvement: - Parameter es$ma$on
More informationFundamentals of Data Assimila1on
2015 GSI Community Tutorial NCAR Foothills Campus, Boulder, CO August 11-14, 2015 Fundamentals of Data Assimila1on Milija Zupanski Cooperative Institute for Research in the Atmosphere Colorado State University
More informationNonlinear ensemble data assimila/on in high- dimensional spaces. Peter Jan van Leeuwen Javier Amezcua, Mengbin Zhu, Melanie Ades
Nonlinear ensemble data assimila/on in high- dimensional spaces Peter Jan van Leeuwen Javier Amezcua, Mengbin Zhu, Melanie Ades Data assimila/on: general formula/on Bayes theorem: The solu/on is a pdf!
More informationQuan&fying Uncertainty. Sai Ravela Massachuse7s Ins&tute of Technology
Quan&fying Uncertainty Sai Ravela Massachuse7s Ins&tute of Technology 1 the many sources of uncertainty! 2 Two days ago 3 Quan&fying Indefinite Delay 4 Finally 5 Quan&fying Indefinite Delay P(X=delay M=
More informationLeast Squares Parameter Es.ma.on
Least Squares Parameter Es.ma.on Alun L. Lloyd Department of Mathema.cs Biomathema.cs Graduate Program North Carolina State University Aims of this Lecture 1. Model fifng using least squares 2. Quan.fica.on
More informationEnsemble 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 informationThe Ensemble Kalman Filter:
p.1 The Ensemble Kalman Filter: Theoretical formulation and practical implementation Geir Evensen Norsk Hydro Research Centre, Bergen, Norway Based on Evensen 23, Ocean Dynamics, Vol 53, No 4 p.2 The Ensemble
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 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 informationSequential Monte Carlo Methods for Bayesian Computation
Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter
More informationSmoothers: Types and Benchmarks
Smoothers: Types and Benchmarks Patrick N. Raanes Oxford University, NERSC 8th International EnKF Workshop May 27, 2013 Chris Farmer, Irene Moroz Laurent Bertino NERSC Geir Evensen Abstract Talk builds
More informationLeast Square Es?ma?on, Filtering, and Predic?on: ECE 5/639 Sta?s?cal Signal Processing II: Linear Es?ma?on
Least Square Es?ma?on, Filtering, and Predic?on: Sta?s?cal Signal Processing II: Linear Es?ma?on Eric Wan, Ph.D. Fall 2015 1 Mo?va?ons If the second-order sta?s?cs are known, the op?mum es?mator is given
More informationGaussian Filtering Strategies for Nonlinear Systems
Gaussian Filtering Strategies for Nonlinear Systems Canonical Nonlinear Filtering Problem ~u m+1 = ~ f (~u m )+~ m+1 ~v m+1 = ~g(~u m+1 )+~ o m+1 I ~ f and ~g are nonlinear & deterministic I Noise/Errors
More informationEnsemble Data Assimila.on for Climate System Component Models
Ensemble Data Assimila.on for Climate System Component Models Jeffrey Anderson Na.onal Center for Atmospheric Research In collabora.on with: Alicia Karspeck, Kevin Raeder, Tim Hoar, Nancy Collins IMA 11
More informationDART 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 informationCS 6140: Machine Learning Spring What We Learned Last Week. Survey 2/26/16. VS. Model
Logis@cs CS 6140: Machine Learning Spring 2016 Instructor: Lu Wang College of Computer and Informa@on Science Northeastern University Webpage: www.ccs.neu.edu/home/luwang Email: luwang@ccs.neu.edu Assignment
More informationCS 6140: Machine Learning Spring 2016
CS 6140: Machine Learning Spring 2016 Instructor: Lu Wang College of Computer and Informa?on Science Northeastern University Webpage: www.ccs.neu.edu/home/luwang Email: luwang@ccs.neu.edu Logis?cs Assignment
More informationModeling radiocarbon in the Earth System. Radiocarbon summer school 2012
Modeling radiocarbon in the Earth System Radiocarbon summer school 2012 Outline Brief introduc9on to mathema9cal modeling Single pool models Mul9ple pool models Model implementa9on Parameter es9ma9on What
More informationLocal 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 informationDART 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 informationEnsemble of Climate Models
Ensemble of Climate Models Claudia Tebaldi Climate Central and Department of Sta7s7cs, UBC Reto Knu>, Reinhard Furrer, Richard Smith, Bruno Sanso Outline Mul7 model ensembles (MMEs) a descrip7on at face
More information4DEnVar. Four-Dimensional Ensemble-Variational Data Assimilation. Colloque National sur l'assimilation de données
Four-Dimensional Ensemble-Variational Data Assimilation 4DEnVar Colloque National sur l'assimilation de données Andrew Lorenc, Toulouse France. 1-3 décembre 2014 Crown copyright Met Office 4DEnVar: Topics
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationShort 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 informationDART Tutorial Part II: How should observa'ons impact an unobserved state variable? Mul'variate assimila'on.
DART Tutorial Part II: How should observa'ons impact an unobserved state variable? Mul'variate assimila'on. UCAR The Na'onal Center for Atmospheric Research is sponsored by the Na'onal Science Founda'on.
More informationEnsemble 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 informationEnsembles and Particle Filters for Ocean Data Assimilation
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Ensembles and Particle Filters for Ocean Data Assimilation Robert N. Miller College of Oceanic and Atmospheric Sciences
More informationEnsemble 4DVAR and observa3on impact study with the GSIbased hybrid ensemble varia3onal data assimila3on system. for the GFS
Ensemble 4DVAR and observa3on impact study with the GSIbased hybrid ensemble varia3onal data assimila3on system for the GFS Xuguang Wang University of Oklahoma, Norman, OK xuguang.wang@ou.edu Ting Lei,
More informationCOMP 562: Introduction to Machine Learning
COMP 562: Introduction to Machine Learning Lecture 20 : Support Vector Machines, Kernels Mahmoud Mostapha 1 Department of Computer Science University of North Carolina at Chapel Hill mahmoudm@cs.unc.edu
More information(Extended) Kalman Filter
(Extended) Kalman Filter Brian Hunt 7 June 2013 Goals of Data Assimilation (DA) Estimate the state of a system based on both current and all past observations of the system, using a model for the system
More informationDART Tutorial Sec'on 9: More on Dealing with Error: Infla'on
DART Tutorial Sec'on 9: More on Dealing with Error: Infla'on 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 informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond Department of Biomedical Engineering and Computational Science Aalto University January 26, 2012 Contents 1 Batch and Recursive Estimation
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 informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More informationAdaptive 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 informationOrganization. I MCMC discussion. I project talks. I Lecture.
Organization I MCMC discussion I project talks. I Lecture. Content I Uncertainty Propagation Overview I Forward-Backward with an Ensemble I Model Reduction (Intro) Uncertainty Propagation in Causal Systems
More informationDART 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 informationEstimation of positive sum-to-one constrained parameters with ensemble-based Kalman filters: application to an ocean ecosystem model
Estimation of positive sum-to-one constrained parameters with ensemble-based Kalman filters: application to an ocean ecosystem model Ehouarn Simon 1, Annette Samuelsen 1, Laurent Bertino 1 Dany Dumont
More informationLagrangian Data Assimilation and Manifold Detection for a Point-Vortex Model. David Darmon, AMSC Kayo Ide, AOSC, IPST, CSCAMM, ESSIC
Lagrangian Data Assimilation and Manifold Detection for a Point-Vortex Model David Darmon, AMSC Kayo Ide, AOSC, IPST, CSCAMM, ESSIC Background Data Assimilation Iterative process Forecast Analysis Background
More informationLagrangian Data Assimilation and Its Application to Geophysical Fluid Flows
Lagrangian Data Assimilation and Its Application to Geophysical Fluid Flows Laura Slivinski June, 3 Laura Slivinski (Brown University) Lagrangian Data Assimilation June, 3 / 3 Data Assimilation Setup:
More informationData assimilation with and without a model
Data assimilation with and without a model Tim Sauer George Mason University Parameter estimation and UQ U. Pittsburgh Mar. 5, 2017 Partially supported by NSF Most of this work is due to: Tyrus Berry,
More informationData assimilation with and without a model
Data assimilation with and without a model Tyrus Berry George Mason University NJIT Feb. 28, 2017 Postdoc supported by NSF This work is in collaboration with: Tim Sauer, GMU Franz Hamilton, Postdoc, NCSU
More informationData 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 informationGaussian 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 informationCSE 473: Ar+ficial Intelligence. Probability Recap. Markov Models - II. Condi+onal probability. Product rule. Chain rule.
CSE 473: Ar+ficial Intelligence Markov Models - II Daniel S. Weld - - - University of Washington [Most slides were created by Dan Klein and Pieter Abbeel for CS188 Intro to AI at UC Berkeley. All CS188
More informationDATA ASSIMILATION FOR FLOOD FORECASTING
DATA ASSIMILATION FOR FLOOD FORECASTING Arnold Heemin Delft University of Technology 09/16/14 1 Data assimilation is the incorporation of measurement into a numerical model to improve the model results
More informationState 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 informationParallelizing Gaussian Process Calcula1ons in R
Parallelizing Gaussian Process Calcula1ons in R Christopher Paciorek UC Berkeley Sta1s1cs Joint work with: Benjamin Lipshitz Wei Zhuo Prabhat Cari Kaufman Rollin Thomas UC Berkeley EECS (formerly) IBM
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Centre for Computational Statistics and Machine Learning University College London c.archambeau@cs.ucl.ac.uk CSML
More informationThe analog data assimilation: method, applications and implementation
The analog data assimilation: method, applications and implementation Pierre Tandeo pierre.tandeo@imt-atlantique.fr IMT-Atlantique, Brest, France 16 th of February 2018 Before starting Works in collaboration
More informationApplied Time Series Analysis FISH 507. Eric Ward Mark Scheuerell Eli Holmes
Applied Time Series Analysis FISH 507 Eric Ward Mark Scheuerell Eli Holmes Introduc;ons Who are we? Who & why you re here? What are you looking to get from this class? Days and Times Lectures When: Tues
More informationParameter Estimation in a Moving Horizon Perspective
Parameter Estimation in a Moving Horizon Perspective State and Parameter Estimation in Dynamical Systems Reglerteknik, ISY, Linköpings Universitet State and Parameter Estimation in Dynamical Systems OUTLINE
More informationAspects of the practical application of ensemble-based Kalman filters
Aspects of the practical application of ensemble-based Kalman filters Lars Nerger Alfred Wegener Institute for Polar and Marine Research Bremerhaven, Germany and Bremen Supercomputing Competence Center
More informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
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 informationA new Hierarchical Bayes approach to ensemble-variational data assimilation
A new Hierarchical Bayes approach to ensemble-variational data assimilation Michael Tsyrulnikov and Alexander Rakitko HydroMetCenter of Russia College Park, 20 Oct 2014 Michael Tsyrulnikov and Alexander
More informationL09. PARTICLE FILTERING. NA568 Mobile Robotics: Methods & Algorithms
L09. PARTICLE FILTERING NA568 Mobile Robotics: Methods & Algorithms Particle Filters Different approach to state estimation Instead of parametric description of state (and uncertainty), use a set of state
More informationThe Ensemble Kalman Filter:
p.1 The Ensemble Kalman Filter: Theoretical formulation and practical implementation Geir Evensen Norsk Hydro Research Centre, Bergen, Norway Based on Evensen, Ocean Dynamics, Vol 5, No p. The Ensemble
More informationConsider the joint probability, P(x,y), shown as the contours in the figure above. P(x) is given by the integral of P(x,y) over all values of y.
ATMO/OPTI 656b Spring 009 Bayesian Retrievals Note: This follows the discussion in Chapter of Rogers (000) As we have seen, the problem with the nadir viewing emission measurements is they do not contain
More informationA Spectral Approach to Linear Bayesian Updating
A Spectral Approach to Linear Bayesian Updating Oliver Pajonk 1,2, Bojana V. Rosic 1, Alexander Litvinenko 1, and Hermann G. Matthies 1 1 Institute of Scientific Computing, TU Braunschweig, Germany 2 SPT
More informationFeb 21 and 25: Local weighted least squares: Quadratic loess smoother
Feb 1 and 5: Local weighted least squares: Quadratic loess smoother An example of weighted least squares fitting of data to a simple model for the purposes of simultaneous smoothing and interpolation is
More informationMul$- model ensemble challenge ini$al/model uncertain$es
Mul$- model ensemble challenge ini$al/model uncertain$es Yuejian Zhu Ensemble team leader Environmental Modeling Center NCEP/NWS/NOAA Acknowledgments: EMC ensemble team staffs Presenta$on for WMO/WWRP
More informationCan hybrid-4denvar match hybrid-4dvar?
Comparing ensemble-variational assimilation methods for NWP: Can hybrid-4denvar match hybrid-4dvar? WWOSC, Montreal, August 2014. Andrew Lorenc, Neill Bowler, Adam Clayton, David Fairbairn and Stephen
More informationSensor Tasking and Control
Sensor Tasking and Control Sensing Networking Leonidas Guibas Stanford University Computation CS428 Sensor systems are about sensing, after all... System State Continuous and Discrete Variables The quantities
More informationApplication of the Ensemble Kalman Filter to History Matching
Application of the Ensemble Kalman Filter to History Matching Presented at Texas A&M, November 16,2010 Outline Philosophy EnKF for Data Assimilation Field History Match Using EnKF with Covariance Localization
More informationA statistical overview and perspectives on data assimilation for marine biogeochemical models
Special Issue Paper Received: 6 October 2013, Revised: 22 January 2014, Accepted: 24 January 2014, Published online in Wiley Online Library (wileyonlinelibrary.com) DOI: 10.1002/env.2264 A statistical
More informationEnKF-based particle filters
EnKF-based particle filters Jana de Wiljes, Sebastian Reich, Wilhelm Stannat, Walter Acevedo June 20, 2017 Filtering Problem Signal dx t = f (X t )dt + 2CdW t Observations dy t = h(x t )dt + R 1/2 dv t.
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
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 informationMachine Learning and Data Mining. Linear regression. Prof. Alexander Ihler
+ Machine Learning and Data Mining Linear regression Prof. Alexander Ihler Supervised learning Notation Features x Targets y Predictions ŷ Parameters θ Learning algorithm Program ( Learner ) Change µ Improve
More informationGEOG 401 Climate Change
GEOG 401 Climate Change Climate Downscaling GCMs have coarse resolu/on Spa
More informationCSE 473: Ar+ficial Intelligence. Example. Par+cle Filters for HMMs. An HMM is defined by: Ini+al distribu+on: Transi+ons: Emissions:
CSE 473: Ar+ficial Intelligence Par+cle Filters for HMMs Daniel S. Weld - - - University of Washington [Most slides were created by Dan Klein and Pieter Abbeel for CS188 Intro to AI at UC Berkeley. All
More informationEnhancing information transfer from observations to unobserved state variables for mesoscale radar data assimilation
Enhancing information transfer from observations to unobserved state variables for mesoscale radar data assimilation Weiguang Chang and Isztar Zawadzki Department of Atmospheric and Oceanic Sciences Faculty
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 informationMethods of Data Assimilation and Comparisons for Lagrangian Data
Methods of Data Assimilation and Comparisons for Lagrangian Data Chris Jones, Warwick and UNC-CH Kayo Ide, UCLA Andrew Stuart, Jochen Voss, Warwick Guillaume Vernieres, UNC-CH Amarjit Budiraja, UNC-CH
More informationCosmological N-Body Simulations and Galaxy Surveys
Cosmological N-Body Simulations and Galaxy Surveys Adrian Pope, High Energy Physics, Argonne Na3onal Laboratory, apope@anl.gov CScADS: Scien3fic Data and Analy3cs for Extreme- scale Compu3ng, 30 July 2012
More informationDecadal Hindcasts and Forecasts at GFDL
Decadal Hindcasts and Forecasts at GFDL Tony Rosati T. Delworth, R. Gudgel, F. Zang, S. Zhang Key Ques;ons What seasonal decadal predictability exists in the climate system, and what are the mechanisms
More informationKalman Filter and Ensemble Kalman Filter
Kalman Filter and Ensemble Kalman Filter 1 Motivation Ensemble forecasting : Provides flow-dependent estimate of uncertainty of the forecast. Data assimilation : requires information about uncertainty
More informationX t = a t + r t, (7.1)
Chapter 7 State Space Models 71 Introduction State Space models, developed over the past 10 20 years, are alternative models for time series They include both the ARIMA models of Chapters 3 6 and the Classical
More informationBayesian statistical data assimilation for ecosystem models using Markov Chain Monte Carlo
Available online at www.sciencedirect.com Journal of Marine Systems 68 (2007) 439 456 www.elsevier.com/locate/jmarsys Bayesian statistical data assimilation for ecosystem models using Markov Chain Monte
More informationGaussian with mean ( µ ) and standard deviation ( σ)
Slide from Pieter Abbeel Gaussian with mean ( µ ) and standard deviation ( σ) 10/6/16 CSE-571: Robotics X ~ N( µ, σ ) Y ~ N( aµ + b, a σ ) Y = ax + b + + + + 1 1 1 1 1 1 1 1 1 1, ~ ) ( ) ( ), ( ~ ), (
More informationAdaptive Data Assimilation and Multi-Model Fusion
Adaptive Data Assimilation and Multi-Model Fusion Pierre F.J. Lermusiaux, Oleg G. Logoutov and Patrick J. Haley Jr. Mechanical Engineering and Ocean Science and Engineering, MIT We thank: Allan R. Robinson
More informationIntroduc)on to Ar)ficial Intelligence
Introduc)on to Ar)ficial Intelligence Lecture 13 Approximate Inference CS/CNS/EE 154 Andreas Krause Bayesian networks! Compact representa)on of distribu)ons over large number of variables! (OQen) allows
More informationLatent Dirichlet Alloca/on
Latent Dirichlet Alloca/on Blei, Ng and Jordan ( 2002 ) Presented by Deepak Santhanam What is Latent Dirichlet Alloca/on? Genera/ve Model for collec/ons of discrete data Data generated by parameters which
More informationApplications of an ensemble Kalman Filter to regional ocean modeling associated with the western boundary currents variations
Applications of an ensemble Kalman Filter to regional ocean modeling associated with the western boundary currents variations Miyazawa, Yasumasa (JAMSTEC) Collaboration with Princeton University AICS Data
More informationA 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 informationExpectation Propagation in Dynamical Systems
Expectation Propagation in Dynamical Systems Marc Peter Deisenroth Joint Work with Shakir Mohamed (UBC) August 10, 2012 Marc Deisenroth (TU Darmstadt) EP in Dynamical Systems 1 Motivation Figure : Complex
More informationStatistical Inference and Methods
Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 31st January 2006 Part VI Session 6: Filtering and Time to Event Data Session 6: Filtering and
More informationBayesian Calibration of Simulators with Structured Discretization Uncertainty
Bayesian Calibration of Simulators with Structured Discretization Uncertainty Oksana A. Chkrebtii Department of Statistics, The Ohio State University Joint work with Matthew T. Pratola (Statistics, The
More information6 Sequential Data Assimilation for Nonlinear Dynamics: The Ensemble Kalman Filter
6 Sequential Data Assimilation for Nonlinear Dynamics: The Ensemble Kalman Filter GEIR EVENSEN Nansen Environmental and Remote Sensing Center, Bergen, Norway 6.1 Introduction Sequential data assimilation
More informationBiogeochemical modelling and data assimilation: status in Australia
Biogeochemical modelling and data assimilation: status in Australia Richard Matear, Andrew Lenton, Matt Chamberlain, Mathieu Mongin, Emlyn Jones, Mark Baird www.cmar.csiro.au/staff/oke/ Biogeochemical
More informationThe minimisation gives a set of linear equations for optimal weights w:
4. Interpolation onto a regular grid 4.1 Optimal interpolation method The optimal interpolation method was used to compute climatological property distributions of the selected standard levels on a regular
More informationPart 2. Representation Learning Algorithms
53 Part 2 Representation Learning Algorithms 54 A neural network = running several logistic regressions at the same time If we feed a vector of inputs through a bunch of logis;c regression func;ons, then
More informationCS 6140: Machine Learning Spring What We Learned Last Week 2/26/16
Logis@cs CS 6140: Machine Learning Spring 2016 Instructor: Lu Wang College of Computer and Informa@on Science Northeastern University Webpage: www.ccs.neu.edu/home/luwang Email: luwang@ccs.neu.edu Sign
More informationExperimental Designs for Planning Efficient Accelerated Life Tests
Experimental Designs for Planning Efficient Accelerated Life Tests Kangwon Seo and Rong Pan School of Compu@ng, Informa@cs, and Decision Systems Engineering Arizona State University ASTR 2015, Sep 9-11,
More informationBayesian Nonlinear Assimilation of Eulerian and Lagrangian Coastal Flow Data
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Bayesian Nonlinear Assimilation of Eulerian and Lagrangian Coastal Flow Data Dr. Pierre F.J. Lermusiaux Department of Mechanical
More informationA State Space Model for Wind Forecast Correction
A State Space Model for Wind Forecast Correction Valrie Monbe, Pierre Ailliot 2, and Anne Cuzol 1 1 Lab-STICC, Université Européenne de Bretagne, France (e-mail: valerie.monbet@univ-ubs.fr, anne.cuzol@univ-ubs.fr)
More information