A comparison of sequential data assimilation schemes for. Twin Experiments
|
|
- Arron Silas Knight
- 5 years ago
- Views:
Transcription
1 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, FL O. M. Smedstad, QinetiQ North America, Stennis Space Center, MS, USA T. M. Chin, University of Miami, Miami, FL F. Counillon, Nansen Center, Bergen, Norway J. M. Brankart and P. Brasseur LEGI, Grenoble, France W. C. Thacker, AOML, NOAA, Miami, FL J. A. Cummings, NRL Monterey, CA LOM-29 - June 2 nd, 29 p.1/31
2 Assimilation Schemes for HYCOM Multivariate Optimal Interpolation (MVOI) J. A. Cummings; Lorenc 1981, Daley 1991, Cummings 25 Ensemble Optimal Interpolation (EnOI) F. Counillon; Evensen 23, Oke 22 Ensemble Reduced Order Information Filter (EnROIF) T. M. Chin; Chin etal 1999,1 Fixed basis variant of the SEEK filter (SEEK) J. M. Brankart and P. Brasseur; Pham et al Objective: To assess the performance of these schemes in identically configured twin experiments assimilating synthetic altimeter and surface temperature obs. Domain & Model Configuration: 1/12 Gulf of Mexico nested in the 1/12 North Atlantic LOM-29 - June 2 nd, 29 p.2/31
3 Data Update Equation Common linear formula for updating the model-forecast x f to obtain data-analysis x a : x a = x f + K(y Hx f ) (1) y is the data to be assimilated H is the observation operator K is an optimization parameter often called the gain matrix K = P f H T (HP f H T + R) 1 (2) P f is the forecast error covariance R is the observation error covariance Main difference among the four methods is in the numerical representation of P f LOM-29 - June 2 nd, 29 p.3/31
4 MVOI P f is separated into several variances and correlations C h = (1 + s h )exp( s h ) (3) C v = (1 + s v )exp( s v ) (4) C f = (1 + s f )exp( s f ) (5) C b = C f C h C v (6) where C h and C v are the horizontal and vertical correlations and s h and s v are scaled distances. analysis on z levels (42 zlevels between -25 m) state variables: U, V, T, S and intf. Pressure dynamical method (Cooper-Haines 1996) for vertical projection simultaneous 3D analysis obs errors uncorrelated LOM-29 - June 2 nd, 29 p.4/31
5 EnOI In EnOI, the forecast covariance matrix is essentially the sample covariance of an ensemble of model states P f EnOI = 1 M 1 M m=1 (x f m x f ) ( x f m x f ) T (7) where x f m is the m th sample of the forecast ensemble, x f is the ensemble mean, and M is the number of samples. The analysis update is computed as: x a = x f + αp f H T (αhp f H T + R) 1 (y Hx f ) (8) where α (, 1] is a parameter introduced to allow different weights for the ensemble and measurements. State Varaibles: UT, VT, UB, VB, PB, DP, T, S, SSH, SST Vertical Projection: static correlations; analysis in observation space LOM-29 - June 2 nd, 29 p.5/31
6 EnROIF EnROIF uses a Markov random field (MRF) to model the forecast error process as e(i, j) = ( i, j) N γ(i, j, i, j) e(i i, j j) + δ(i, j) (9) where N specifies a set of local grid locations, γ is the regression coefficient (small matrix), and δ(i, j) is a white noise with unit variance. Computation for the state analysis x a in ROIF is performed as L a (x a x f ) = H T R 1 (y Hx f ) (1) where the sparsely banded analysis information matrix L a = L + H T R 1 H is numerically inverted. State Varaibles: DP, SSH, SST Vertical Projection: static correlations; 2D vertically decoupled analysis LOM-29 - June 2 nd, 29 p.6/31
7 Fixed basis SEEK In SEEK filter, P f, is assumed to be of low rank, M and is represented by dominant modes of empirical orthogonal functions (EOFs) as P f SEEK = 1 M 1 M m=1 S f ms ft m (11) where S f m, m = 1,..., M, are the M most dominant EOF modes. In the SEEK analysis the Kalman Gain is rewritten using the Sherman-Morrrison-Woodberry matrix identity as K = S f [I + ( HS f ) R 1 ( HS f ) T ] 1 ( HS f ) R 1 (12) State Varaibles: UT, VT, UB, VB, PB, DP, T, S, SSH, SST Vertical Projection: static correlations; analysis in reduced space LOM-29 - June 2 nd, 29 p.7/31
8 Twin Experiments: Initial States 3 o N Truth State Aug SSH Assimilation Run State Aug SSH 3 o N Latitude Latitude 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W temp. sect N temp. sect N depth [m] 5 1 depth [m] Longitude Longitude LOM-29 - June 2 nd, 29 p.8/31
9 Twin Experiments: Data 3 o N altimeter 3 Aug. 9 Sep o N MCSST 3 Aug Latitude Latitude 96 o W 92 o W 88 o W 84 o W 8 o W Longitude 96 o W 92 o W 88 o W 84 o W 8 o W 76 o W Longitude altimeter data Aug Dec MCSST data Aug Dec data count 5 data count Days 5 1 Days LOM-29 - June 2 nd, 29 p.9/31
10 Twin Experiments: SSH SSH (Forecast) RMS Error MVOI EnOI SEEK EnROIF.16 RMS error [m] Days LOM-29 - June 2 nd, 29 p.1/31
11 Twin Experiments: SSH Sea Surface Elevation (m) (Forecast Truth) Oct 18, 1999 (day 5) MVOI EnOI 3 o N 3 o N 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W SEEK EnROIF 3 o N 3 o N 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W LOM-29 - June 2 nd, 29 p.11/31
12 Twin Experiments: SSH Sea Surface Elevation (m) (Forecast Truth) Dec 7, 1999 (day 1) MVOI EnOI 3 o N 3 o N 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W SEEK EnROIF 3 o N 3 o N 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W LOM-29 - June 2 nd, 29 p.12/31
13 Twin Experiments: SST.8.7 SST (Forecast) RMS Error MVOI EnOI SEEK EnROIF.6 RMS error [degc] Days LOM-29 - June 2 nd, 29 p.13/31
14 Twin Experiments: SST Sea Surface Temperature (degc) (Forecast Truth) Oct 18, 1999 (day 5) MVOI EnOI 3 o N 3 o N 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W SEEK EnROIF 3 o N 3 o N 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W LOM-29 - June 2 nd, 29 p.14/31
15 Twin Experiments: Unobserved Variables [ u,v, t,s] RMS error [m/s] U velocity (Forecast) RMS Error RMS error [m/s] V velocity (Forecast) RMS Error MVOI EnOI SEEK EnROIF Temperature (Forecast) RMS Error 2 Salinity (Forecast) RMS Error.35 RMS error [degc] Days RMS error [psu] Days LOM-29 - June 2 nd, 29 p.15/31
16 Twin Experiments: Temperature(degC) section 25.4 N (Forecast Truth) Aug 31, 1999 (day 2) MVOI EnOI depth [m] Longitude Longitude SEEK EnROIF depth [m] LOM June 2 nd, 29 p.16/31
17 Twin Experiments: Temperature(degC) section 25.4 N (Forecast Truth) Oct 18, 1999 (day 5) MVOI EnOI depth [m] Longitude Longitude SEEK EnROIF depth [m] LOM June 2 nd, 29 p.17/31
18 Twin Experiments: depth of 2 deg isotherm (m) (Forecast Truth) Aug 31, 1999 (day 2) MVOI EnOI 3 o N 3 o N 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W SEEK EnROIF 3 o N 3 o N 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W LOM June 2 nd, 29 p.18/31
19 Twin Experiments: depth of 2 deg isotherm (m) (Forecast Truth) Oct 18, 1999 (day 5) MVOI EnOI 3 o N 3 o N 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W SEEK EnROIF 3 o N 3 o N 96 o W 92 o W 88 o W 84 o W 8 o W 96 o W 92 o W 88 o W 84 o W 8 o W LOM June 2 nd, 29 p.19/31
20 Twin Experiments: 3D U&V u velocity day 2 u velocity day 5 depth [m] MVOI 4 EnOI 6 SEEK EnROIF v velocity day 2 v velocity day 5 depth [m] rms error [m/s] rms error [m/s] LOM-29 - June 2 nd, 29 p.2/31
21 Twin Experiments: 3D T&S Temperature day 2 Temperature day 5 depth [m] RMS error [degc] Salinity day 2 2 MVOI 4 EnOI 6 SEEK EnROIF RMS error [degc] Salinity day 2 depth [m] RMS error [psu] depth [m] RMS error [psu] LOM-29 - June 2 nd, 29 p.21/31
22 Twin Experiments: Layer Thickness Mean Layer Thickness 3 Layer 5 6 Layer 6 Truth MVOI Layer Layer EnOI SEEK EnROIF depth [m] Layer Layer Layer Days Layer Days LOM-29 - June 2 nd, 29 p.22/31
23 Twin Experiments: Layer Thickness Mean Layer Thickness 55 Layer Layer 14 Truth MVOI Layer Layer EnOI SEEK EnROIF depth [m] Layer Layer Layer Layer Days Days LOM-29 - June 2 nd, 29 p.23/31
24 Twin Experiments: Net Transports.6.4 Net transports across 26 N MVOI EnOI SEEK EnROIF.2 Transport (SV) Days LOM-29 - June 2 nd, 29 p.24/31
25 Twin Experiments: State Variables Impact of State Variables EnOI_142 SEEK_154 SEEK_155 SEEK_156 SEEK_157 MVOI_12 MVOI_116 RMS error [m] Days LOM-29 - June 2 nd, 29 p.25/31
26 Twin Experiments: Ensemble Size and Covariance Rank Impact of Ensemble Size.2.18 EnOI_146 EnOI_142 SEEK_155 SEEK_162 SEEK_163 SEEK_ RMS error [m] Days LOM-29 - June 2 nd, 29 p.26/31
27 Twin Experiments: Localization Impact of Localization.2.18 EnOI_142 EnOI_144 EnOI_145 SEEK_153 SEEK_ RMS error [m] Days LOM-29 - June 2 nd, 29 p.27/31
28 Twin Experiments: Data Window Impact of Data Window.2.18 MVOI_12 MVOI_121 EnOI_142 SEEK_155 SEEK_ RMS error [m] Days LOM-29 - June 2 nd, 29 p.28/31
29 Twin Experiments: Assimilation Strength Impact of Assimilation Strength.2 EnOI_142 EnOI_143 EnOI_ RMS error [m] Days LOM-29 - June 2 nd, 29 p.29/31
30 Computational Costs Scheme Scalability Memory Req. Typical Wall Clock/update MVOI O(m 3 /V ) 22 min EnOI O(m 3 ) nn 31 min EnROIF sn 8 min SEEK O(r 3 ) nn 14 min n ==> no of members in the ensemble N ==> size of the forecast model s ==> size of the MRF neighborhood r ==> covariance rank (== no of members in our experiments) v == no of analysis volumes in MVOI all experiments were done on 8 core/2.6 GHz Intel/12GB RAM) LOM-29 - June 2 nd, 29 p.3/31
31 Summary All schemes are able to reduce errors in both observed and unobserved variables. 3D multivariate (MVOI, EnOI, SEEK) are clearly better than decoupled 2D scheme used currently in EnROIF Best performance is obtained when all state variables are used in the estimation space Subsurface correction based on correlations(enoi and SEEK) are better balanced and seem more effective than the dynamical method(mvoi) The nearly identical performance of EnOI and SEEK is both expected and reassuring Several parameters (ensemble size, covariance rank, radius of localization etc.) have to be explored for a given problem to find a good solution that is also cost effective. LOM-29 - June 2 nd, 29 p.31/31
Demonstration 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 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 informationDemonstration and Comparison of Sequential Approaches to Data Assimilation in HYCOM
Demonstration and Comparison of Sequential Approaches to Data Assimilation in HYCOM A. Srinivasan 1, E. Chassignet, O.M. Smedstad, W. C. Thacker, L. Bertino, P. Brasseur, T.M. Chin, F. Counillon, and J.
More informationImplementation of the SEEK filter in HYCOM
Implementation of the SEEK filter in HYCOM P. Brasseur, J. Verron, J.M. Brankart LEGI/CNRS, Grenoble, France HYCOM model SST, SSH, ocean colour Assimilation SEEK filter Y H x f In situ, XBT K x a Real-time
More informationEnhancing predictability of the Loop Current variability using Gulf of Mexico Hycom
Enhancing predictability of the Loop Current variability using Gulf of Mexico Hycom Matthieu Le Hénaff (1) Villy Kourafalou (1) Ashwanth Srinivasan (1) Collaborators: O. M. Smedstad (2), P. Hogan (2),
More informationImpact of frontal eddy dynamics on the Loop Current variability during free and data assimilative HYCOM simulations
Impact of frontal eddy dynamics on the Loop Current variability during free and data assimilative HYCOM simulations Matthieu Le Hénaff (1) Villy H. Kourafalou (1) Ashwanth Srinivasan (1) George R. Halliwell
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 informationOverview of the NRL Coupled Ocean Data Assimilation System
Overview of the NRL Coupled Ocean Data Assimilation System James Cummings and Rich Hodur Naval Research Laboratory, Monterey, CA HYCOM Ocean Model Workshop NCEP August 2003 Report Documentation Page Form
More informationNCODA Implementation with re-layerization
NCODA Implementation with re-layerization HeeSook Kang CIMAS/RSMAS/U. Miami with W. Carlisle Thacker NOAA/AOML HYCOM meeting December 6 2005 1 GULF OF MEXICO MODEL CONFIGURATION: Horizontal grid: 1/12
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 informationOOPC-GODAE workshop on OSE/OSSEs Paris, IOCUNESCO, November 5-7, 2007
OOPC-GODAE workshop on OSE/OSSEs Paris, IOCUNESCO, November 5-7, 2007 Design of ocean observing systems: strengths and weaknesses of approaches based on assimilative systems Pierre Brasseur CNRS / LEGI
More informationAlexander Barth, Aida Alvera-Azc. Azcárate, Robert H. Weisberg, University of South Florida. George Halliwell RSMAS, University of Miami
Ensemble-based based Assimilation of HF-Radar Surface Currents in a West Florida Shelf ROMS Nested into HYCOM and filtering of spurious surface gravity waves. Alexander Barth, Aida Alvera-Azc Azcárate,
More informationJohn Kindle. Sergio derada Igor Shulman Ole Martin Smedstad Stephanie Anderson. Data Assimilation in Coastal Modeing April
John Kindle Sergio derada Igor Shulman Ole Martin Smedstad Stephanie Anderson Data Assimilation in Coastal Modeing April 3 2007 MODELS Motivation: Global->Coastal Real-Time Regional Coastal Models Global
More informationOcean data assimilation for reanalysis
Ocean data assimilation for reanalysis Matt Martin. ERA-CLIM2 Symposium, University of Bern, 14 th December 2017. Contents Introduction. On-going developments to improve ocean data assimilation for reanalysis.
More informationAssimilation of satellite altimetry referenced to the new GRACE geoid estimate
GEOPHYSICAL RESEARCH LETTERS, VOL. 32, L06601, doi:10.1029/2004gl021329, 2005 Assimilation of satellite altimetry referenced to the new GRACE geoid estimate F. Birol, 1 J. M. Brankart, 1 J. M. Lemoine,
More informationData Assimilation with the Ensemble Kalman Filter and the SEIK Filter applied to a Finite Element Model of the North Atlantic
Data Assimilation with the Ensemble Kalman Filter and the SEIK Filter applied to a Finite Element Model of the North Atlantic L. Nerger S. Danilov, G. Kivman, W. Hiller, and J. Schröter Alfred Wegener
More informationEarly results and plans for the future. Robert Atlas
Observing System Simulation Experiments: Methodology, Early results and plans for the future Robert Atlas National Oceanic and Atmospheric Administration Atlantic Oceanographic and Meteorological Laboratory
More informationPerformance of a 23 years TOPAZ reanalysis
Performance of a 23 years TOPAZ reanalysis L. Bertino, F. Counillon, J. Xie,, NERSC LOM meeting, Copenhagen, 2 nd -4 th June 2015 Outline Presentation of the TOPAZ4 system Choice of modeling and assimilation
More informationLOM Layered Ocean Model Workshop, 2-4 June 2015, Copenhagen
LOM Layered Ocean Model Workshop, 2-4 June 2015, Copenhagen 2 HYCOM Preliminary operational system Grid 1/3 (22 levels) Spin-up COADS 30 years (USP-UFBA) 2007/2008 NCEP atmospheric forcing HYCOM Grid nesting
More informationO.M Smedstad 1, E.J. Metzger 2, R.A. Allard 2, R. Broome 1, D.S. Franklin 1 and A.J. Wallcraft 2. QinetiQ North America 2. Naval Research Laboratory
An eddy-resolving ocean reanalysis using the 1/12 global HYbrid Coordinate Ocean Model (HYCOM) and the Navy Coupled Ocean Data Assimilation (NCODA) scheme O.M Smedstad 1, E.J. Metzger 2, R.A. Allard 2,
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 informationSensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts
Sensitivity of West Florida Shelf Simulations to Initial and Boundary Conditions Provided by HYCOM Data-Assimilative Ocean Hindcasts George Halliwell, MPO/RSMAS, University of Miami Alexander Barth, University
More informationPreliminary Test of Glider Data Assimilation Along the Labrador Sea Shelf Break
Preliminary Test of Glider Data Assimilation Along the Labrador Sea Shelf Break Third Annual VITALS Science Meeting October 19, 2015 Changheng Chen, K. Andrea Scott Department of Systems Design Engineering
More informationAn Update on the 1/12 Global HYCOM Effort
An Update on the 1/12 Global HYCOM Effort E.J. Metzger 1, O.M. Smedstad 2, A.J. Wallcraft 1, P.G. Posey 1, and D.S. Franklin 2 1 Naval Research Laboratory 2 Qinetiq North America 2013 Layered Ocean Model
More informationGlobal HYCOM and Advanced Data Assimilation
Global HYCOM and Advanced Data Assimilation Harley E. Hurlburt, Robert C. Rhodes, Gregg A. Jacobs Naval Research Laboratory Stennis Space Center, MS 39529-5004 phone: (228) 688-4626 fax:(228) 688-4759
More informationStrongly coupled data assimilation: Could scatterometer winds improve ocean analyses?
Strongly coupled data assimilation: Could scatterometer winds improve ocean analyses? Sergey Frolov 1 Craig H. Bishop 2 ; Teddy Holt 2 ; James Cummings 3 ; David Kuhl 4 1 UCAR, Visiting Scientist 2 Naval
More informationMultigrid state vector for data assimilation in. a two-way nested model of the Ligurian Sea
Multigrid state vector for data assimilation in a two-way nested model of the Ligurian Sea A. Barth a,,1, A. Alvera-Azcárate a,1, J.-M. Beckers a M. Rixen b, L. Vandenbulcke a a University of Liége, GHER,
More informationJOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. C3, 3074, doi: /2001jc001198, 2003
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. C3, 3074, doi:10.1029/2001jc001198, 2003 Implementation of a multivariate data assimilation scheme for isopycnic coordinate ocean models: Application to a
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 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 informationInterannual Variability of the Gulf of Mexico Loop Current
Interannual Variability of the Gulf of Mexico Loop Current Dmitry Dukhovskoy (COAPS FSU) Eric Chassignet (COAPS FSU) Robert Leben (UC) Acknowledgements: O. M. Smedstad (Planning System Inc.) J. Metzger
More informationData Assimilation Research Testbed Tutorial
Data Assimilation Research Testbed Tutorial Section 2: How should observations of a state variable impact an unobserved state variable? Multivariate assimilation. Single observed variable, single unobserved
More 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 informationWP2 task 2.2 SST assimilation
WP2 task 2.2 SST assimilation Matt Martin, Dan Lea, James While, Rob King ERA-CLIM2 General Assembly, Dec 2015. Contents Impact of SST assimilation in coupled DA SST bias correction developments Assimilation
More informationInitial Results of Altimetry Assimilation in POP2 Ocean Model
Initial Results of Altimetry Assimilation in POP2 Ocean Model Svetlana Karol and Alicia R. Karspeck With thanks to the DART Group, CESM Software Development Group, Climate Modeling Development Group POP-DART
More informationDART_LAB Tutorial Section 2: How should observations impact an unobserved state variable? Multivariate assimilation.
DART_LAB Tutorial Section 2: How should observations impact an unobserved state variable? Multivariate assimilation. UCAR 2014 The National Center for Atmospheric Research is sponsored by the National
More informationInhomogeneous and Nonstationary Feature Analysis: Melding of Oceanic Variability and Structure (INFAMOVS)
Inhomogeneous and Nonstationary Feature Analysis: Melding of Oceanic Variability and Structure (INFAMOVS) Dr. Arthur J. Mariano Rosenstiel School of Marine and Atmospheric Science Division of Meteorology
More informationHYCOM in the South Atlantic Bight: Performance and Client Applications
HYCOM in the South Atlantic Bight: Performance and Client Applications Brian Blanton, Alfredo Aretxabaleta Department of Marine Sciences UNC-Chapel Hill UNC Group Activities HYCOM/GODAE NOPP Provide SEACOOS
More informationThe Bureau of Meteorology Coupled Data Assimilation System for ACCESS-S
The Bureau of Meteorology Coupled Data Assimilation System for ACCESS-S Yonghong Yin, Angus Gray-Weale, Oscar Alves, Pavel Sakov, Debra Hudson, Xiaobing Zhou, Hailing Yan, Mei Zhao Research and Development
More informationThe Center for Ocean Atmospheric Prediction Studies, FSU, Tallahassee, FL 32306
Global Ocean Prediction Using HYCOM E. Joseph Metzger 1 joe.metzger@nrlssc.navy.mil Eric P. Chassignet 2 echassignet@coaps.fsu.edu James A. Cummings 3 james.cummings@nrlmry.navy.mil Harley E. Hurlburt
More informationModel and observation bias correction in altimeter ocean data assimilation in FOAM
Model and observation bias correction in altimeter ocean data assimilation in FOAM Daniel Lea 1, Keith Haines 2, Matt Martin 1 1 Met Office, Exeter, UK 2 ESSC, Reading University, UK Abstract We implement
More informationDefense Technical Information Center Compilation Part Notice
UNCLASSIFIED Defense Technical Information Center Compilation Part Notice ADP023750 TITLE: Global Ocean Prediction Using HYCOM DISTRIBUTION: Approved for public release, distribution unlimited This paper
More informationUniversity of Miami/RSMAS
Observing System Simulation Experiments to Evaluate the Potential Impact of Proposed Observing Systems on Hurricane Prediction: R. Atlas, T. Vukicevic, L.Bucci, B. Annane, A. Aksoy, NOAA Atlantic Oceanographic
More informationBalanced multivariate model errors of an intermediate coupled model for ensemble Kalman filter data assimilation
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 113,, doi:10.1029/2007jc004621, 2008 Balanced multivariate model errors of an intermediate coupled model for ensemble Kalman filter data
More informationArgo data assimilation into HYCOM with an EnOI method in the Atlantic Ocean
doi:10.5194/os-11-195-2015 Author(s) 2015. CC Attribution 3.0 License. Argo data assimilation into HYCOM with an EnOI method in the Atlantic Ocean D. Mignac 1,2, C. A. S. Tanajura 2,3,4, A. N. Santana
More informationHYCOM Caspian Sea Modeling. Part I: An Overview of the Model and Coastal Upwelling. Naval Research Laboratory, Stennis Space Center, USA
HYCOM Caspian Sea Modeling. Part I: An Overview of the Model and Coastal Upwelling By BIROL KARA, ALAN WALLCRAFT AND JOE METZGER Naval Research Laboratory, Stennis Space Center, USA MURAT GUNDUZ Institute
More informationAn Overview of Nested Regions Using HYCOM
An Overview of Nested Regions Using HYCOM Patrick Hogan Alan Wallcraft Luis Zamudio Sergio DeRada Prasad Thoppil Naval Research Laboratory Stennis Space Center, MS 10 th HYCOM Consortium Meeting COAPS,
More informationA Modeling Study on Flows in the Strait of Hormuz (SOH)
A Modeling Study on Flows in the Strait of Hormuz (SOH) Peter C Chu & Travis Clem Naval Postgraduate School Monterey, CA 93943, USA IUGG 2007: PS005 Flows and Waves in Straits. July 5-6, Perugia, Italy
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 informationApplication and improvement of Ensemble Optimal Interpolation on Regional Ocean Modeling System (ROMS)
Application and improvement of Ensemble Optimal Interpolation on Regional Ocean Modeling System (ROMS) Zhaoyi Wang Guokun Lyu, Hui Wang, Jiang Zhu, Guimei Liu et al. National Marine Environmental Forecasting
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 informationSensitivity of Satellite Altimetry Data Assimilation on a Naval Anti-Submarine Warfare Weapon System
Sensitivity of Satellite Altimetry Data Assimilation on a Naval Anti-Submarine Warfare Weapon System Steven Mancini LCDR, USN Advisor: Prof. Peter Chu Second Readers: Dr. Charlie Barron Eric Gottshall,
More informationOverview of data assimilation in oceanography or how best to initialize the ocean?
Overview of data assimilation in oceanography or how best to initialize the ocean? T. Janjic Alfred Wegener Institute for Polar and Marine Research Bremerhaven, Germany Outline Ocean observing system Ocean
More informationModelling forecast error statistics in the Mercator ocean and sea-ice reanalysis system.
Modelling forecast error statistics in the Mercator ocean and sea-ice reanalysis system. C.E Testut 1, G. Ruggiero 1, L. Parent 1, J.M. Lellouche 1, O. Legalloudec 1, C. Bricaud 1, J. Chanut 1, G. Smith
More informationstatistical methods for tailoring seasonal climate forecasts Andrew W. Robertson, IRI
statistical methods for tailoring seasonal climate forecasts Andrew W. Robertson, IRI tailored seasonal forecasts why do we make probabilistic forecasts? to reduce our uncertainty about the (unknown) future
More informationDevelopments to the assimilation of sea surface temperature
Developments to the assimilation of sea surface temperature James While, Daniel Lea, Matthew Martin ERA-CLIM2 General Assembly, January 2017 Contents Introduction Improvements to SST bias correction Development
More informationA Dressed Ensemble Kalman Filter Using the Hybrid Coordinate Ocean Model in the Pacific
ADVANCES IN ATMOSPHERIC SCIENCES, VOL. 26, NO. 5, 2009, 1042 1052 A Dressed Ensemble Kalman Filter Using the Hybrid Coordinate Ocean Model in the Pacific WAN Liying 1 (!"#), ZHU Jiang 2 ($%), WANG Hui
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 informationObservations assimilated
Observations assimilated R-factor to adjust model error variances and K-Factor to adjust observation error variances. ETKF employs SST and SLA bias correction via AR() model fit. TABLE Summary of ocean
More informationNested Gulf of Mexico Modeling with HYCOM
Nested Gulf of Mexico Modeling with HYCOM Patrick J. Hogan Alan J. Wallcraft Naval Research Laboratory Stennis Space Center, MS HYCOM Meeting December 6-8, 2005 University of Miami, Miami, FL 1/25 (~4
More informationComparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales
Comparing Variational, Ensemble-based and Hybrid Data Assimilations at Regional Scales Meng Zhang and Fuqing Zhang Penn State University Xiang-Yu Huang and Xin Zhang NCAR 4 th EnDA Workshop, Albany, NY
More informationDynamical Evaluation of Gulf Stream Simulations in Models with High Vertical Resolution
Dynamical Evaluation of Gulf Stream Simulations in Models with High Vertical Resolution Harley E. Hurlburt 1, Eric P. Chassignet 2, E. Joseph Metzger 1, James G. Richman 1, William J. Schmitz, Jr. 3, Jay
More informationArctic Cap Nowcast Forecast System (ACNFS) end of summer 2013 Ice Extent Projection July Report
Arctic Cap Nowcast Forecast System (ACNFS) end of summer 2013 Ice Extent Projection July Report Naval Research Laboratory, Stennis Space Center, MS The NRL Ice Team consists of: Pamela Posey 1, E. Joseph
More informationAssimilation Impact of Physical Data on the California Coastal Ocean Circulation and Biogeochemistry
Assimilation Impact of Physical Data on the California Coastal Ocean Circulation and Biogeochemistry Yi Chao, Remote Sensing Solutions (RSS)/UCLA; John D. Farrara, RSS; Fei Chai, University of Maine; Hongchun
More informationCoupled Ocean-Atmosphere Assimilation
Coupled Ocean-Atmosphere Assimilation Shu-Chih Yang 1, Eugenia Kalnay 2, Joaquim Ballabrera 3, Malaquias Peña 4 1:Department of Atmospheric Sciences, National Central University 2: Department of Atmospheric
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 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 informationApplication of 3-D ensemble variational data assimilation to a Baltic Sea reanalysis
PUBLISHED BY THE INTERNATIONAL METEOROLOGICAL INSTITUTE IN STOCKHOLM SERIES A DYNAMIC METEOROLOGY AND OCEANOGRAPHY Application of 3-D ensemble variational data assimilation to a Baltic Sea reanalysis 989
More informationCoastal Data Assimilation: progress and challenges in state estimation for circulation and BGC models
Coastal Data Assimilation: progress and challenges in state estimation for circulation and BGC models Emlyn Jones, Roger Scott, Mark Baird, Frank Colberg, Paul Sandery, Pavel Sakov, Gary Brassington, Mathieu
More informationConstrained data assimilation. W. Carlisle Thacker Atlantic Oceanographic and Meteorological Laboratory Miami, Florida USA
Constrained data assimilation W. Carlisle Thacker Atlantic Oceanographic and Meteorological Laboratory Miami, Florida 33149 USA Plan Range constraints: : HYCOM layers have minimum thickness. Optimal interpolation:
More informationT2.3: Use of ensemble information in ocean analysis and development of efficient 4D-Var
T2.3: Use of ensemble information in ocean analysis and development of efficient 4D-Var A. Weaver 1,B.Ménétrier 1,J.Tshimanga 1 and A. Vidard 2 1 CERFACS, Toulouse 2 INRIA/LJK, Grenoble December 10, 2015
More informationApplication and validation of polynomial chaos methods to quantify uncertainties in simulating the Gulf of Mexico circulation using HYCOM.
Application and validation of polynomial chaos methods to quantify uncertainties in simulating the Gulf of Mexico circulation using HYCOM. Mohamed Iskandarani Matthieu Le Hénaff Carlisle Thacker University
More informationActivityoftheJapanCoastalOcean PredictabilityExperiment -RecentprogressoftheKuroshio forecastandadownscalingstudy- TakashiKagimoto(FRCGC/JAMSTEC)
ActivityoftheJapanCoastalOcean PredictabilityExperiment -RecentprogressoftheKuroshio forecastandadownscalingstudy- TakashiKagimoto(FRCGC/JAMSTEC) JCOPEmember YasumasaMiyazawa(FRCGC/JAMSTEC) XinyuGuo(CMES,EhimeUniv.)
More informationGFDL, NCEP, & SODA Upper Ocean Assimilation Systems
GFDL, NCEP, & SODA Upper Ocean Assimilation Systems Jim Carton (UMD) With help from Gennady Chepurin, Ben Giese (TAMU), David Behringer (NCEP), Matt Harrison & Tony Rosati (GFDL) Description Goals Products
More informationStatus of 1/25 Global HYCOM Simulations
Status of 1/25 Global HYCOM Simulations Luis Zamudio 1 & Joe Metzger 2 1 Center for Ocean-Atmospheric Prediction Studies, Florida State University 2 Naval Research Laboratory, Stennis Space Center, Mississippi
More informationMERSEA IP WP 5 Integrated System Design and Assessment. Internal Metrics for the MERSEA Global Ocean: Specifications for Implementation
MERSEA IP WP 5 Integrated System Design and Assessment Internal Metrics for the MERSEA Global Ocean: Specifications for Implementation Planning Workshop CSIRO Hobart, Australia 17-18 March 7 Participants
More informationUsing sea level data to constrain a finite-element primitive-equation ocean model with a local SEIK filter
Noname manuscript No. (will be inserted by the editor) Using sea level data to constrain a finite-element primitive-equation ocean model with a local SEIK filter L. Nerger(1) S. Danilov W. Hiller J. Schröter
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 informationHorizon Marine Inc. Numerical Ocean Model Nowcast and Forecast Skill Assessment in the Gulf of Mexico
Horizon Marine Inc. Numerical Ocean Model Nowcast and Forecast Skill Assessment in the Gulf of Mexico Robert Leben and Kevin Corcoran Colorado Center for Astrodynamics Research University of Colorado,
More informationThe Ensemble Kalman Filter: Theoretical Formulation and Practical Implementation
Noname manuscript No. (will be inserted by the editor) The Ensemble Kalman Filter: Theoretical Formulation and Practical Implementation Geir Evensen Nansen Environmental and Remote Sensing Center, Bergen
More informationPreparation of the SWOT Mission
Preparation of the SWOT Mission M.Benkiran, E. Greiner, E. Rémy, P.Y. Le Traon and the Mercator Ocean team. Study done in the framework of a CNES/Mercator Ocean convention, in collaboration with CLS. GODAE
More informationApplication of a hybrid EnKF-OI to ocean forecasting
Ocean Sci., 5, 389 401, 2009 Author(s) 2009. This work is distributed under the Creative Commons Attribution 3.0 License. Ocean Science Application of a hybrid EnKF-OI to ocean forecasting F. Counillon,
More informationGlobal Ocean Reanalysis Simulations at Mercator Océan GLORYS1: the Argo years
Global Ocean Reanalysis Simulations at Mercator Océan GLORYS1: the Argo years 2002-2008 Laurent Parent, Nicolas Ferry and the Mercator Océan team Mercator-Océan, Toulouse, France: http://www.mercator-ocean.fr
More informationOptimal Spectral Decomposition (OSD) for GTSPP Data Analysis
Optimal Spectral Decomposition (OSD) for GTSPP Data Analysis Peter C Chu (1),Charles Sun (2), & Chenwu Fan (1) (1) Naval Postgraduate School, Monterey, CA 93943 pcchu@nps.edu, http://faculty.nps.edu/pcchu/
More informationNorwegian Climate Prediction Model (NorCPM) getting ready for CMIP6 DCPP
Norwegian Climate Prediction Model (NorCPM) getting ready for CMIP6 DCPP Francois Counillon, Noel Keenlyside, Mats Bentsen, Ingo Bethke, Laurent Bertino, Teferi Demissie, Tor Eldevik, Shunya Koseki, Camille
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 informationStatus and future of data assimilation in operational oceanography
Status and future of data assimilation in operational oceanography MJ Martin, Met Office, Exeter, UK. M Balmaseda, ECMWF, Reading, UK L Bertino, NERSC, Bergen, Norway P Brasseur, LEGI, Grenoble, France
More informationAn ETKF approach for initial state and parameter estimation in ice sheet modelling
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
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 informationChallenges and Advances of Regional Ocean Data Assimilation. Andy Moore Department of Ocean Science University of California Santa Cruz
Challenges and Advances of Regional Ocean Data Assimilation Andy Moore Department of Ocean Science University of California Santa Cruz The Large Scale Ocean Circulation Diverse Approach Global versus Regional
More informationAssimilating Altimetry Data into a HYCOM Model of the Pacific: Ensemble Optimal Interpolation versus Ensemble Kalman Filter
APRIL 2010 W A N E T A L. 753 Assimilating Altimetry Data into a HYCOM Model of the Pacific: Ensemble Optimal Interpolation versus Ensemble Kalman Filter LIYING WAN National Marine Environmental Forecasting
More informationA demonstration of ensemble-based assimilation methods with a layered OGCM from the perspective of operational ocean forecasting systems
Journal of Marine Systems 40 41 (2003) 253 289 www.elsevier.com/locate/jmarsys A demonstration of ensemble-based assimilation methods with a layered OGCM from the perspective of operational ocean forecasting
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 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 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 informationOcean Data Assimilation: a case for ensemble optimal interpolation
Ocean Data Assimilation: a case for ensemble optimal interpolation Peter R. Oke 1,2, Gary B. Brassington 1,3, David A. Griffin 1,2, Andreas Schiller 1,2 1 Centre for Australian Weather and Climate Research
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 informationDevelopment and Application of the Chinese Operational Hydrological Forecasting System
Development and Application of the Chinese Operational Hydrological Forecasting System Guimei LIU National Marine Environmental Forecasting Center, Beijing www.nmefc.gov.cn Contents Background Operational
More informationThe TOPAZ3 forecasting system. L. Bertino, K.A. Lisæter, Mohn-Sverdrup Center/NERSC
The TOPAZ3 forecasting system L. Bertino, K.A. Lisæter, Mohn-Sverdrup Center/NERSC LOM meeting, Bergen Beach, 20 th Aug. 2007 Introduction Main objective of data assimilation Estimate the most likely
More informationConcurrent simulation of the eddying general circulation and tides in a global ocean model
Concurrent simulation of the eddying general circulation and tides in a global ocean model Brian K. Arbic 1 E. Joseph Metzger 2 Alan J. Wallcraft 2 1 Department of Oceanography and Center for Ocean-Atmospheric
More information