An ETKF approach for initial state and parameter estimation in ice sheet modelling
|
|
- Esmond Gibbs
- 6 years ago
- Views:
Transcription
1 An ETKF approach for initial state and parameter estimation in ice sheet modelling Bertrand Bonan University of Reading DARC Seminar, Reading October 30, 2013 DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
2 Context Work supervised by Maëlle Nodet 1,2 and Catherine Ritz 3. Joined work also with Vincent Peyaud 3. 1 : Inria, Laboratoire Jean Kuntzmann (Grenoble) 2 : Université Joseph Fourier Grenoble 1 3 : CNRS, Laboratoire de Glaciologie et Géophysique de l Environnement (Grenoble) Part of project ADAGe funded by French Agence Nationale de la Recherche. DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
3 Outline 1 Problem description 2 Large-scale ice sheet model 3 Data assimilation Twin experiment ETKF and LETKF Initial ensemble Results DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
4 Problem description Outline 1 Problem description 2 Large-scale ice sheet model 3 Data assimilation Twin experiment ETKF and LETKF Initial ensemble Results DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
5 Problem description Sea level change: Antarctica & Greenland contribution Motivations: Computation of the ice discharge of Antarctica and Greenland in the near future, thanks to simulations of polar ice sheet model. Ice discharge: narrow outlets, closely linked to ice velocities, highly sensitive to basal friction parameters, highly sensitive to bedrock topography. Surface ice velocities [Rignot et al. 2011] DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
6 Problem description Basal parameters Typically 2 types of basal parameters are mandatory for simulation: bedrock topography B soc (x) basal friction parameters in basal shear stress relationship { C(x, t) ub (x, t) τ b (x, t) = m(x,t) 1 u b (x, t) if T b (x, t) = T melt 0 if T b (x, t) < T melt with u b basal velocity, T b basal ice temperature, C > 0 and m basal friction parameters. Simplified relationship { β(x)ub (x, t) if T τ b (x, t) = b (x, t) = T melt 0 if T b (x, t) < T melt with β > 0. DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
7 Problem description What is behind β due to sliding of ice over the bed deformation of the bed effective pressure = ice pressure - water pressure (ex: thin water layer between ice and bedrock, water in cavities,...) rock debris sediments... DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
8 Problem description Poorly known basal parameters... Illustration for bedrock topography with Bedmap2 [Fretwell et al. 2013] Bedmap2 estimation of bed topography DARC Seminar (Reading) Estimated uncertainty of bed topography ETKF for ice sheet initialisation October 30, / 50
9 Problem description... Fortunately surface observations Ice velocity Surface elevation (ex.: [Bamber et al, 2009]) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
10 Problem description Overview What we want: Build estimation of actual basal parameters: bedrock topography B soc sliding coefficient β = 10 α Give an estimation of the uncertainty of our estimation. What we have: 20 years of observations of surface elevation and surface velocities Sparse observations of bedrock topography Time dependent ice sheet model What we need: efficient data assimilation system To develop the method: Synthetic experiments along a flowline DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
11 Large-scale ice sheet model Outline 1 Problem description 2 Large-scale ice sheet model 3 Data assimilation Twin experiment ETKF and LETKF Initial ensemble Results DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
12 Large-scale ice sheet model Ice dynamics processes DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
13 Large-scale ice sheet model Simplified physics DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
14 Large-scale ice sheet model Model equations: mass balance Flowline SIA model (1D + time) with a sliding law Mass balance equation: H t = ḃm (UH) x H t=0 = H 0 with x latitude, t time H(x, t) ice thickness, H 0 (x) initial ice thickness U(x, t) ice velocity averaged over ice thickness ḃ m (x, t) surface mass balance rate DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
15 Large-scale ice sheet model Model equations: dynamics (1) Vertically averaged ice velocity is a diagnostic variable = no partial derivative in time involved, computed from geometry at each time step with S(x, t) surface elevation U = u def + u slid S = B soc + H H 0 B soc (x, t) bedrock topography u def ice velocity component due to ice deformation u slid ice velocity component due to sliding DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
16 Large-scale ice sheet model Model equations: dynamics (2) u def is the averaged velocity due to ice deformation. Under SIA approximation S H 2 ( ) S 3 u def = a 1 x 3 a H 4 2 x 5 with a 1, a 2 given coefficients u slid is the velocity component due to basal sliding with β > 0 basal friction coefficient. u slid = 1 β ρ i gh S x DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
17 Large-scale ice sheet model Model equations: mass balance rate ḃ m (x, t) = Acc(x, t) + Abl(x, t) Acc(x, t) 0 ice accumulation rate (mainly snow precipitation) Abl(x, t) 0 ice ablation rate (snow melting) Both are function of surface temperature T s T s (x, t) = F clim (t) + λ x + γs(x, t) (en m/an) Accumulation Ablation Bilan de masse Température en surface (en ) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
18 Outline 1 Problem description 2 Large-scale ice sheet model 3 Data assimilation Twin experiment ETKF and LETKF Initial ensemble Results DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
19 Twin experiment Outline 1 Problem description 2 Large-scale ice sheet model 3 Data assimilation Twin experiment ETKF and LETKF Initial ensemble Results DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
20 Twin experiment Twin experiment What we control: at each grid point Ice thickness H Bedrock topography B soc Basal parameter β What we observe: each year during 20 years (20 observation times) Surface elevation S at each grid point (σ S = 2 m) Surface ice velocity u s at each grid point (σ u = 3 m/yr) Bedrock topography B soc at few grid points (σ B = 20 m) DA system: ETKF [Wang et al., 2004], [Ott et al. 2004] LETKF [Hunt et al., 2007] DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
21 Twin experiment Background and reference states (1) altitude (in m) reference bedrock observed surface backround bedrock x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
22 Twin experiment Background and reference states (2) ! reference! backround! x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
23 Twin experiment Background and reference states (3) velocity (in m/yr) sliding component deformation component surface velocity x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
24 ETKF and LETKF Outline 1 Problem description 2 Large-scale ice sheet model 3 Data assimilation Twin experiment ETKF and LETKF Initial ensemble Results DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
25 ETKF and LETKF Basis of Ensemble Kalman Filter (EnKF) Ensemble formulation: Let x k ensemble of N ens state estimations of model at time t k (in our case, ice thickness, sliding coefficient and bedrock topography) x (i) k ith member of ensemble, i = 1,..., N ens x k ensemble mean P e,k ensemble covariance matrix 2 step algorithm: Forecast: Predict x f k from xa k 1 using the model. Analysis: Correct x f k with observations yo k (in our case, surface velocities, surface elevation and sparse bedrock topography) to produce x a k. DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
26 Data assimilation ETKF and LETKF Sequential data assimilation: EnKF sequence state x f k 2 yk 1 o y o k 2 y o k k 2 k 1 k t DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
27 Data assimilation ETKF and LETKF Sequential data assimilation: EnKF sequence state x f k 2 yk 1 o x a k 2 y o k 2 y o k k 2 k 1 k t DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
28 Data assimilation ETKF and LETKF Sequential data assimilation: EnKF sequence state x f k 2 yk 1 o x a k 2 y o k 2 y x f k o k 1 k 2 k 1 k t DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
29 Data assimilation ETKF and LETKF Sequential data assimilation: EnKF sequence state x f k 2 yk 1 o x a k 2 x a k 1 y o k 2 y x f k o k 1 k 2 k 1 k t DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
30 Data assimilation ETKF and LETKF Sequential data assimilation: EnKF sequence state x f k 2 yk 1 o x a k 2 x a k 1 x f k y o k 2 y x f k o k 1 k 2 k 1 k t DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
31 Data assimilation ETKF and LETKF x Sequential data assimilation: EnKF sequence state x f k 2 yk 1 o y o k 2 x a k 2 x a k 1 x f k 1 x f k k 2 k 1 k t y o k a k DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
32 ETKF and LETKF Ensemble Transform Kalman Filter (ETKF) Detailled in [Hunt et al. 2007], [Harlim and Hunt 2007]. Analysis step } Build an ensemble {x a(i), i = 1... N ens such as x a minimum of J (x) = 1 2 ( x x f ) T P f e 1 ( x x f ) (yo H(x)) T R 1 (y o H(x)) and P a e P a. Attention: rank(p f e) N ens 1 by construction so P f e not invertible. DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
33 ETKF and LETKF First hypothesis Transform Filter w R Nens Assume x = x f + X f w with X f matrix whose ith column is x f (i) x f Why? w R Nens, x R nx and minimise J (x f + X f w) and N ens < n x (generally). ( x x f ) T P f 1 ( e x x f ) = w T X f T P f 1 e X f w (background term) }) if V = Vect ({x f (i) x f, i = 1... N ens, P f e = 1 N 1 Xf X f T invertible as a restriction to V. DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
34 ETKF and LETKF Where are we? J (x f + X f w) = N 1 ( ) w T X f T X f X f T 1 X f w ( y o H(x f + X w)) ( T ) f R 1 y o H(x f + X f w) 2 Attention: J invariant on Ker(X f ) and dim(ker(x f )) 1. = non unicity of w for minimal J. DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
35 ETKF and LETKF Gauge-fixing term Minimise under constraints To avoid this invariance we can minimise J over the w that have a null orthogonal projection on the kernel of X f ( ( ) ) I Nens X f T X f X f T 1 X f w = 0 It is equivalent to add a gauge-fixing term G into J G (w) = N 1 ( ( ) w T I Nens X f T X f X f T 1 X f 2 and minimise the new cost function J. ) w DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
36 ETKF and LETKF Where are we? J (w) = N 1 w T w+ 1 ( y o H(x f + X w)) ( T ) f R 1 y o H(x f + X f w) 2 2 DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
37 ETKF and LETKF Second hypothesis Linear approximation for observation operator ( ) H x f + X f w y f + Y f w with y f (i) = H (x f (i)) y f = 1 N N ens y f (i) i=1 matrix Y f whose ith column is y f (i) y f Finally, we search to minimise a quadratic cost function J (w) = N 1 w T w + 1 ( y o y f Y w) ( T ) f R 1 y o y f Y f w 2 2 DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
38 ETKF and LETKF ETKF analysis step Direct computation of: w a = P ( ) a Y f T R 1 y o y f (minimum of J ) ( ) P a = (N ens 1) I Nens + Y f T 1 ( 1 R 1 Y f (= Hess J (w )) a ) ( ) 1/2 W a = (N ens 1) P a (symetric square root matrix) } Analysis ensemble {x a(i), i = 1... N ens is defined as x a = x f + X f w a (mean) X a = X f W a (anomalies matrix) x a(i) = x a + X a i (with X a i ith column of X a ) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
39 Initial ensemble Outline 1 Problem description 2 Large-scale ice sheet model 3 Data assimilation Twin experiment ETKF and LETKF Initial ensemble Results DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
40 Initial ensemble Initial ensemble: Cooking recipe (1) Generate an ensemble of α = log(β) following a Gaussian law N (α b, B α ) with [B α ] ij = σα 2 exp ( (x i x j ) 2 ) L Generate an ensemble of B soc following a Gaussian law N (B b soc, B b ) Generate an ensemble of S following a Gaussian law N (S obs, σ 2 S I) Compute H = S B soc and correct to avoid negative ice thicknesses. For each member, run the model 1 year in order to obtain more physically balanced ice sheets. Rescale the produced ensemble so that its mean remains equal to the background. Warning: we use observations to produce initial ensemble! DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
41 Initial ensemble Initial ensemble: Cooking recipe (2) B b = ΣCΣ with: Σ square root of diagonal matrix of variances C correlation matrix standart deviation (in m) x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
42 Initial ensemble Initial ensemble: Cooking recipe (3) [C] ij = c 1 exp ( (x i x j ) 2 ) + c 2 exp ( (x i x j ) 2 ) L 1 L 2 1 large scale component short scale component correlation function x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
43 Initial ensemble Initial ensemble 4000 Initial ice sheet ensemble (Ne = 50) altitude (in m) x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
44 Results Outline 1 Problem description 2 Large-scale ice sheet model 3 Data assimilation Twin experiment ETKF and LETKF Initial ensemble Results DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
45 Results ETKF (N ens = 1000) 3000 ETKF result for ice sheet (ensemble size = 1000) altitude (in m) reference background final analysis observed surface x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
46 Results ETKF (N ens = 1000) 800 Standard deviation for bedrock topography (ensemble size = 1000) background final analysis Standard deviation (in m) x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
47 Results ETKF (N ens = 1000) ETKF result for! (ensemble size = 1000) ! reference background final analysis x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
48 Results ETKF (N ens = 1000) 3000 ETKF result for sliding velocity (ensemble size = 1000) velocity (in m/yr) reference background final analysis x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
49 Results ETKF (N ens = 1000) 2500 Standard deviation for sliding velocity (ensemble size = 1000) background final analysis 2000 standard deviation (in m/yr) x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
50 Results LETKF (N ens = 30, 50, 100) LETKF results for ice sheet (ensemble size = 100, 50, 30) reference background final anal. (100 m.) final anal. (50 m.) final anal. (30 m.) observed surface altitude (in m) x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
51 Results LETKF (N ens = 30, 50, 100) 800 Standard deviation for bedrock topography (ensemble size = 100, 50, 30) background final anal. (100 m.) final anal. (50 m.) final anal. (30 m.) standard deviation (in m) x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
52 Results LETKF (N ens = 30, 50, 100) 2500 LETKF results for sliding velocity (ensemble size = 100, 50, 30) sliding velocity (in m/yr) reference background final anal. (100 m.) 2000 final anal. (50 m.) final anal. (30 m.) x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
53 Results LETKF (N ens = 30, 50, 100) Standard deviation for sliding velocity (ensemble size = 100, 50, 30) background final anal. (100 m.) final anal. (50 m) final anal. (30 m.) standard deviation (in m/yr) x axis (in km) DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
54 Conclusion Ontgoing and future works Compare with variational approaches used in glaciology Add more physics: ice shelves, criterion on basal temperature for example. Perform experiments on a real flow line Implement this method into a 2D ice sheet model. DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
55 Conclusion Thank you for your attention! DARC Seminar (Reading) ETKF for ice sheet initialisation October 30, / 50
Asynchronous data assimilation
Ensemble Kalman Filter, lecture 2 Asynchronous data assimilation Pavel Sakov Nansen Environmental and Remote Sensing Center, Norway This talk has been prepared in the course of evita-enkf project funded
More informationSummary for the Greenland ice sheet
Contribution of Greenland and Antarctica to future sea level change Catherine Ritz, Gaël Durand, Fabien Gillet-Chaulet, Olivier Gagliardini, Vincent Peyaud EDGe team, LGGE, CNRS/UJF Grenoble, France Ice
More informationMulti-Modal Flow in a Thermocoupled Model of the Antarctic Ice Sheet, with Verification
Multi-Modal Flow in a Thermocoupled Model of the Antarctic Ice Sheet, with Verification Craig Lingle 1 Jed Brown 2 Ed Bueler 2 1 Geophysical Institute University of Alaska Fairbanks, USA 2 Department of
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 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 informationIce in the climate system. Summary so far. Today. The Cryosphere. 1. Climate history of the Earth. 2. Paleo observations (1)
Ice in the climate system 1. Climate history of the Earth 2. Paleo observations (1) 3. Paleo observations (2) 4. Ice ages 5. Climate sensitivity 6. Ice in the climate system Summary so far Radiation (Milankovitch
More informationA data-driven method for improving the correlation estimation in serial ensemble Kalman filter
A data-driven method for improving the correlation estimation in serial ensemble Kalman filter Michèle De La Chevrotière, 1 John Harlim 2 1 Department of Mathematics, Penn State University, 2 Department
More informationThe Local Ensemble Transform Kalman Filter (LETKF) Eric Kostelich. Main topics
The Local Ensemble Transform Kalman Filter (LETKF) Eric Kostelich Arizona State University Co-workers: Istvan Szunyogh, Brian Hunt, Ed Ott, Eugenia Kalnay, Jim Yorke, and many others http://www.weatherchaos.umd.edu
More informationAn Efficient Ensemble Data Assimilation Approach To Deal With Range Limited Observation
An Efficient Ensemble Data Assimilation Approach To Deal With Range Limited Observation A. Shah 1,2, M. E. Gharamti 1, L. Bertino 1 1 Nansen Environmental and Remote Sensing Center 2 University of Bergen
More informationLocal 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 informationAddressing the nonlinear problem of low order clustering in deterministic filters by using mean-preserving non-symmetric solutions of the ETKF
Addressing the nonlinear problem of low order clustering in deterministic filters by using mean-preserving non-symmetric solutions of the ETKF Javier Amezcua, Dr. Kayo Ide, Dr. Eugenia Kalnay 1 Outline
More informationConvergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit
Convergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit Evan Kwiatkowski, Jan Mandel University of Colorado Denver December 11, 2014 OUTLINE 2 Data Assimilation Bayesian Estimation
More informationSummer School in Glaciology, Fairbanks/McCarthy, Exercises: Glacial geology
Bob Anderson Summer School in Glaciology, Fairbanks/McCarthy, 2010 Exercises: Glacial geology 1. Glacier thickness. We wish to estimate the local thickness of a glacier given only a topographic map of
More informationPractical 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 informationComparison of of Assimilation Schemes for HYCOM
Comparison of of Assimilation Schemes for HYCOM Ashwanth Srinivasan, C. Thacker, Z. Garraffo, E. P. Chassignet, O. M. Smedstad, J. Cummings, F. Counillon, L. Bertino, T. M. Chin, P. Brasseur and C. Lozano
More informationError (E) Covariance Diagnostics. in Ensemble Data Assimilation. Kayo Ide, Daryl Kleist, Adam Kubaryk University of Maryland
Error (E) Covariance Diagnostics in Ensemble Data Assimilation Kayo Ide, Daryl Kleist, Adam Kubaryk University of Maryland 7 th EnKF Data Assimilation Workshop May 2016 Outline Part I. Background: E-diagnostics
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 informationFinite-size ensemble Kalman filters (EnKF-N)
Finite-size ensemble Kalman filters (EnKF-N) Marc Bocquet (1,2), Pavel Sakov (3,4) (1) Université Paris-Est, CEREA, joint lab École des Ponts ParisTech and EdF R&D, France (2) INRIA, Paris-Rocquencourt
More informationIce sheet mass balance from satellite altimetry. Kate Briggs (Mal McMillan)
Ice sheet mass balance from satellite altimetry Kate Briggs (Mal McMillan) Outline Background Recap 25 year altimetry record Recap Measuring surface elevation with altimetry Measuring surface elevation
More informationAssimilation des Observations et Traitement des Incertitudes en Météorologie
Assimilation des Observations et Traitement des Incertitudes en Météorologie Olivier Talagrand Laboratoire de Météorologie Dynamique, Paris 4èmes Rencontres Météo/MathAppli Météo-France, Toulouse, 25 Mars
More informationReview of Covariance Localization in Ensemble Filters
NOAA Earth System Research Laboratory Review of Covariance Localization in Ensemble Filters Tom Hamill NOAA Earth System Research Lab, Boulder, CO tom.hamill@noaa.gov Canonical ensemble Kalman filter update
More informationNew developments in data assimilation in MIKE 21/3 FM Assimilation of along-track altimetry data with correlated measurement errors
New developments in data assimilation in MIKE 21/3 FM Assimilation of along-track altimetry data with correlated measurement errors EnKF Workshop 2016-06-20 Jesper Sandvig Mariegaard Henrik Andersson DHI
More informationOcean Ice Interactions: A
Ocean Ice Interactions: A Cryospheric Perspective Tony Payne a.j.payne@bristol.ac.uk Steph Cornford, Rupert Gladstone and Dan Martin (LLNL) KISS short course Sept. 2013 Slide number 1/37 Outline Evidence
More informationModeled and observed fast flow in the Greenland ice sheet
Modeled and observed fast flow in the Greenland ice sheet Ed Bueler 1 Constantine Khroulev 1 Andy Aschwanden 2 Ian Joughin 3 1 Dept of Mathematics and Statistics, University of Alaska Fairbanks 2 Arctic
More informationPISM, a Parallel Ice Sheet Model: Current status of our Antarctic ice sheet simulation
PISM, a Parallel Ice Sheet Model: Current status of our Antarctic ice sheet simulation Craig Lingle, 1 Ed Bueler, 2 Jed Brown, 1 and David Covey, 1 1 Geophysical Institute, Univ. of Alaska, Fairbanks 2
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 informationDemonstration and Comparison of of Sequential Approaches for Altimeter Data Assimilation in in HYCOM
Demonstration and Comparison of of Sequential Approaches for Altimeter Data Assimilation in in HYCOM A. Srinivasan, E. P. Chassignet, O. M. Smedstad, C. Thacker, L. Bertino, P. Brasseur, T. M. Chin,, F.
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 informationA comparison of sequential data assimilation schemes for. Twin Experiments
A comparison of sequential data assimilation schemes for ocean prediction with HYCOM Twin Experiments A. Srinivasan, University of Miami, Miami, FL E. P. Chassignet, COAPS, Florida State University, Tallahassee,
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 informationEffect of Ocean Warming on West Antarctic Ice Streams and Ice Shelves. Bryan Riel December 4, 2008
Effect of Ocean Warming on West Antarctic Ice Streams and Ice Shelves Bryan Riel December 4, 2008 Ice Sheet Mass Balance/WAIS Dynamics -Mass Balance = (Ice/Snow Accumulation) (Surface melting, ice outflux,
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 informationGlacier (and ice sheet) Mass Balance
Glacier (and ice sheet) Mass Balance The long-term average position of the highest (late summer) firn line is termed the Equilibrium Line Altitude (ELA) Firn is old snow How an ice sheet works (roughly):
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 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 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 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 informationModelling of surface to basal hydrology across the Russell Glacier Catchment
Modelling of surface to basal hydrology across the Russell Glacier Catchment Sam GAP Modelling Workshop, Toronto November 2010 Collaborators Alun Hubbard Centre for Glaciology Institute of Geography and
More informationEstimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations.
Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. May 6, 2009 Motivation Constitutive Equations EnKF algorithm Some results Method Navier Stokes equations
More informationBELISSIMA: BELgian Ice Sheet- Shelf Ice Measurements in Antarctica
BELISSIMA: BELgian Ice Sheet- Shelf Ice Measurements in Antarctica Frank PATTYN 1, Jean-Louis TISON 1, Denis SAMYN 1, Kenichi MATSUOKA², Howard CONWAY², Bryn HUBBARD³ (1)Laboratoire de Glaciologie, DSTE,
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 informationErgodicity in data assimilation methods
Ergodicity in data assimilation methods David Kelly Andy Majda Xin Tong Courant Institute New York University New York NY www.dtbkelly.com April 15, 2016 ETH Zurich David Kelly (CIMS) Data assimilation
More informationModel 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 informationDATA ASSIMILATION OF VIS / NIR REFLECTANCE INTO THE DETAILED SNOWPACK MODEL SURFEX/ISBA- CROCUS
DATA ASSIMILATION OF VIS / NIR REFLECTANCE INTO THE DETAILED SNOWPACK MODEL SURFEX/ISBA- CROCUS Luc CHARROIS 1,2 PhD 2013 / 2016 E. Cosme 2, M. Dumont 1, M. Lafaysse 1, S. Morin 1, Q. Libois 2, G. Picard
More informationParameter Estimation in Reservoir Engineering Models via Data Assimilation Techniques
Parameter Estimation in Reservoir Engineering Models via Data Assimilation Techniques Mariya V. Krymskaya TU Delft July 6, 2007 Ensemble Kalman Filter (EnKF) Iterative Ensemble Kalman Filter (IEnKF) State
More informationFrank PATTYN Laboratoire de Glaciologie Université Libre de Bruxelles Belgium
Frank PATTYN Laboratoire de Glaciologie Université Libre de Bruxelles Belgium 200 x 10 6 people living in coastal floodplains IPCC, 2007 2 x 10 6 km 2 and 10 12 US$ worth of assets lying
More informationImproved analyses and forecasts with AIRS retrievals using the Local Ensemble Transform Kalman Filter
Improved analyses and forecasts with AIRS retrievals using the Local Ensemble Transform Kalman Filter Hong Li, Junjie Liu, and Elana Fertig E. Kalnay I. Szunyogh, E. J. Kostelich Weather and Chaos Group
More informationHow do glaciers form?
Glaciers What is a Glacier? A large mass of moving ice that exists year round is called a glacier. Glaciers are formed when snowfall exceeds snow melt year after year Snow and ice remain on the ground
More informationApport altimétrie pour modélisation hydrologique des grands fleuves
Apport altimétrie pour modélisation hydrologique des grands fleuves C. Emery, V. Pedinotti, S. Biancamaria, A. Boone, A. Paris, S. Calmant, S. Ricci, P.-A. Garambois, B. Decharme, Workshop Niger, 14-15
More informationAccounting for Missing Data in Sparse Wavelet Representation of Observation Error Correlations
Accounting for Missing Data in Sparse Wavelet Representation of Observation Error Correlations Vincent Chabot, Maëlle Nodet, Arthur Vidard To cite this version: Vincent Chabot, Maëlle Nodet, Arthur Vidard.
More informationRobust Ensemble Filtering With Improved Storm Surge Forecasting
Robust Ensemble Filtering With Improved Storm Surge Forecasting U. Altaf, T. Buttler, X. Luo, C. Dawson, T. Mao, I.Hoteit Meteo France, Toulouse, Nov 13, 2012 Project Ensemble data assimilation for storm
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 informationTowards a probabilistic hydrological forecasting and data assimilation system. Henrik Madsen DHI, Denmark
Towards a probabilistic hydrological forecasting and data assimilation system Henrik Madsen DHI, Denmark Outline Hydrological forecasting Data assimilation framework Data assimilation experiments Concluding
More informationIntroduction 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 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 informationIntroduction to Data Assimilation
Introduction to Data Assimilation Alan O Neill Data Assimilation Research Centre University of Reading What is data assimilation? Data assimilation is the technique whereby observational data are combined
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 informationIII Subjective probabilities. III.4. Adaptive Kalman filtering. Probability Course III:4 Bologna 9-13 February 2015
III Subjective probabilities III.4. Adaptive Kalman filtering III.4.1 A 2-dimensional Kalman filter system Obs-Tfc= correction Bias? Corr = A(t) Tfc It appears as if the old bias has abruptly changed into
More informationGlaciology (as opposed to Glacial Geology) Why important? What are glaciers? How do they work?
Glaciology (as opposed to Glacial Geology) Why important? What are glaciers? How do they work? Glaciers are important because of their role in creating glacial landscapes (erosional and depositional features).
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 informationSummary of activities with SURFEX data assimilation at Météo-France. Jean-François MAHFOUF CNRM/GMAP/OBS
Summary of activities with SURFEX data assimilation at Météo-France Jean-François MAHFOUF CNRM/GMAP/OBS Outline Status at Météo-France in 2008 Developments undertaken during 2009-2011 : Extended Kalman
More informationThe Canadian approach to ensemble prediction
The Canadian approach to ensemble prediction ECMWF 2017 Annual seminar: Ensemble prediction : past, present and future. Pieter Houtekamer Montreal, Canada Overview. The Canadian approach. What are the
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 informationSome ideas for Ensemble Kalman Filter
Some ideas for Ensemble Kalman Filter Former students and Eugenia Kalnay UMCP Acknowledgements: UMD Chaos-Weather Group: Brian Hunt, Istvan Szunyogh, Ed Ott and Jim Yorke, Kayo Ide, and students Former
More informationComparisons between 4DEnVar and 4DVar on the Met Office global model
Comparisons between 4DEnVar and 4DVar on the Met Office global model David Fairbairn University of Surrey/Met Office 26 th June 2013 Joint project by David Fairbairn, Stephen Pring, Andrew Lorenc, Neill
More informationSimulating a glacier: application to Engabreen
Utrecht University Institute for Marine and Atmospheric Research Utrecht Bachelor thesis Simulating a glacier: application to Engabreen Supervisors: Carleen Tijm-Reijmer Hans Oerlemans Author: Jesse Reusen
More informationRepresentation of inhomogeneous, non-separable covariances by sparse wavelet-transformed matrices
Representation of inhomogeneous, non-separable covariances by sparse wavelet-transformed matrices Andreas Rhodin, Harald Anlauf German Weather Service (DWD) Workshop on Flow-dependent aspects of data assimilation,
More informationApplication of Mean Recentering Scheme to Improve the Typhoon Track Forecast: A Case Study of Typhoon Nanmadol (2011) Chih-Chien Chang, Shu-Chih Yang
6 th EnKF workshop Application of Mean Recentering Scheme to Improve the Typhoon Track Forecast: A Case Study of Typhoon Nanmadol (2011) Chih-Chien Chang, Shu-Chih Yang National Central University, Taiwan
More informationQuarterly Journal of the Royal Meteorological Society
Quarterly Journal of the Royal Meteorological Society Effects of sequential or simultaneous assimilation of observations and localization methods on the performance of the ensemble Kalman filter Journal:
More informationIce Sheets and Climate Change. William H. Lipscomb Los Alamos National Laboratory
Ice Sheets and Climate Change William H. Lipscomb Los Alamos National Laboratory What is an expert? An expert is somebody who is more than 50 miles from home, has no responsibility for implementing the
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 informationChanging Landscapes: Glaciated Landscapes. How do glaciers move?
Changing Landscapes: Glaciated Landscapes How do glaciers move? What you need to know Differences between cold-and warm-based glaciers, their locations and rates of movement Glacier ice movement including
More informationEnsemble-Based Data Assimilation of GPM/DP R Reflectivity into the Nonhydrostatic Icosahed ral Atmospheric Model NICAM
Ensemble-Based Data Assimilation of GPM/DP R Reflectivity into the Nonhydrostatic Icosahed ral Atmospheric Model NICAM Shunji Kotsuki1, Koji Terasaki1, Shigenori Otsuka1, Kenta Kurosawa1, and Takemasa
More informationData assimilation in the geosciences An overview
Data assimilation in the geosciences An overview Alberto Carrassi 1, Olivier Talagrand 2, Marc Bocquet 3 (1) NERSC, Bergen, Norway (2) LMD, École Normale Supérieure, IPSL, France (3) CEREA, joint lab École
More informationThe Use of a Self-Evolving Additive Inflation in the CNMCA Ensemble Data Assimilation System
The Use of a Self-Evolving Additive Inflation in the CNMCA Ensemble Data Assimilation System Lucio Torrisi and Francesca Marcucci CNMCA, Italian National Met Center Outline Implementation of the LETKF
More informationIce on Earth: An overview and examples on physical properties
Ice on Earth: An overview and examples on physical properties - Ice on Earth during the Pleistocene - Present-day polar and temperate ice masses - Transformation of snow to ice - Mass balance, ice deformation,
More informationConstraining simulated atmospheric states by sparse empirical information
Constraining simulated atmospheric states by sparse empirical information Martin Widmann University of Birmingham Hugues Goosse University of Louvain-la-Neuve NCAR workshop Data assimilation and climate
More informationT. Perron Glaciers 1. Glaciers
T. Perron 12.001 Glaciers 1 Glaciers I. Why study glaciers? [PPT: Perito Moreno glacier, Argentina] Role in freshwater budget o Fraction of earth s water that is fresh (non-saline): 3% o Fraction of earth
More informationLocalization in the ensemble Kalman Filter
Department of Meteorology Localization in the ensemble Kalman Filter Ruth Elizabeth Petrie A dissertation submitted in partial fulfilment of the requirement for the degree of MSc. Atmosphere, Ocean and
More informationThermodynamics of Glaciers
Thermodynamics of Glaciers McCarthy Summer School, June 2010 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA McCarthy Summer School, June 2010 1 / 34 On Notation
More informationCross-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 informationESS 431 Principles of Glaciology ESS 505 The Cryosphere
MID-TERM November 9, 2015 ESS 431 Principles of Glaciology ESS 505 The Cryosphere Instructions: Please answer the following 5 questions. [The actual 5 questions will be selected from these 12 questions
More informationLecture 10 Glaciers and glaciation
Lecture 10 Glaciers and glaciation Outline Importance of ice to people! Basics of glaciers formation, classification, mechanisms of movement Glacial landscapes erosion and deposition by glaciers and the
More informationModel Uncertainty Quantification for Data Assimilation in partially observed Lorenz 96
Model Uncertainty Quantification for Data Assimilation in partially observed Lorenz 96 Sahani Pathiraja, Peter Jan Van Leeuwen Institut für Mathematik Universität Potsdam With thanks: Sebastian Reich,
More informationGlacial Modification of Terrain
Glacial Modification Part I Stupendous glaciers and crystal snowflakes -- every form of animate or inanimate existence leaves its impress upon the soul of man. 1 -Orison Swett Marden Glacial Modification
More informationThe Canadian Precipitation Analysis (CaPA) and the Canadian Land Data Assimilation System (CaLDAS)
The Canadian Precipitation Analysis (CaPA) and the Canadian Land Data Assimilation System (CaLDAS) Vincent Fortin (and many collaborators) Division de la recherche en météorologie Environnement en Changement
More informationAn Ensemble Kalman Filter for NWP based on Variational Data Assimilation: VarEnKF
An Ensemble Kalman Filter for NWP based on Variational Data Assimilation: VarEnKF Blueprints for Next-Generation Data Assimilation Systems Workshop 8-10 March 2016 Mark Buehner Data Assimilation and Satellite
More informationNew Applications and Challenges In Data Assimilation
New Applications and Challenges In Data Assimilation Met Office Nancy Nichols University of Reading 1. Observation Part Errors 1. Applications Coupled Ocean-Atmosphere Ensemble covariances for coupled
More informationEcon 2148, fall 2017 Gaussian process priors, reproducing kernel Hilbert spaces, and Splines
Econ 2148, fall 2017 Gaussian process priors, reproducing kernel Hilbert spaces, and Splines Maximilian Kasy Department of Economics, Harvard University 1 / 37 Agenda 6 equivalent representations of the
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 information"Ice Sheets and Sea Level Rise: How Should IPCC Handle Deep Uncertainty?" Michael Oppenheimer For Inside the IPCC Princeton University 1 April 2008
"Ice Sheets and Sea Level Rise: How Should IPCC Handle Deep Uncertainty?" Michael Oppenheimer For Inside the IPCC Princeton University 1 April 2008 This Talk is about: IPCCs (controversial) assessment
More informationExploring The Polar Connection to Sea Level Rise NGSS Disciplinary Core Ideas Science & Engineering Crosscutting Concepts
Exploring The Polar Connection to Sea Level Rise NGSS Disciplinary Core Ideas Science & Engineering Crosscutting Concepts Practices MS - ESS: Earth & Space Science 1. Ask questions 2. Developing and using
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 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 informationSeminar: Data Assimilation
Seminar: Data Assimilation Jonas Latz, Elisabeth Ullmann Chair of Numerical Mathematics (M2) Technical University of Munich Jonas Latz, Elisabeth Ullmann (TUM) Data Assimilation 1 / 28 Prerequisites Bachelor:
More informationUsing Observations at Different Spatial. Scales in Data Assimilation for. Environmental Prediction. Joanne A. Waller
UNIVERSITY OF READING DEPARTMENT OF MATHEMATICS AND STATISTICS Using Observations at Different Spatial Scales in Data Assimilation for Environmental Prediction Joanne A. Waller Thesis submitted for the
More informationIce Sheet Modeling and Sea Level Rise. William Lipscomb, NCAR CESM Sea Level Session 10 January 2018
Ice Sheet Modeling and Sea Level Rise William Lipscomb, NCAR CESM Sea Level Session 10 January 2018 Ice sheets in IPCC AR4 The IPCC Fourth Assessment Report (AR4) projected 0.18 to 0.59 m of sea level
More informationReport on the Joint SRNWP workshop on DA-EPS Bologna, March. Nils Gustafsson Alex Deckmyn.
Report on the Joint SRNWP workshop on DA-EPS Bologna, 22-24 March Nils Gustafsson Alex Deckmyn http://www.smr.arpa.emr.it/srnwp/ Purpose of the workshop On the one hand, data assimilation techniques require
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 informationAn introduction to data assimilation. Eric Blayo University of Grenoble and INRIA
An introduction to data assimilation Eric Blayo University of Grenoble and INRIA Data assimilation, the science of compromises Context characterizing a (complex) system and/or forecasting its evolution,
More information