Show simple item record

dc.contributor.authorGopalakrishnan, Ganesh
dc.contributor.authorHoteit, Ibrahim
dc.contributor.authorCornuelle, Bruce D.
dc.contributor.authorRudnick, Daniel L.
dc.date.accessioned2019-04-09T06:53:18Z
dc.date.available2019-04-09T06:53:18Z
dc.date.issued2019-03-31
dc.identifier.citationGopalakrishnan G, Hoteit I, Cornuelle BD, Rudnick DL (2019) Comparison of 4DVAR and EnKF state estimates and forecasts in the Gulf of Mexico. Quarterly Journal of the Royal Meteorological Society. Available: http://dx.doi.org/10.1002/qj.3493.
dc.identifier.issn0035-9009
dc.identifier.issn1477-870X
dc.identifier.doi10.1002/qj.3493
dc.identifier.urihttp://hdl.handle.net/10754/631839
dc.description.abstractAn experiment is conducted to compare four-dimensional variational (4DVAR) and ensemble Kalman filter (EnKF) assimilation systems and their predictability in the Gulf of Mexico (GoM) using the Massachusetts Institute of Technology general circulation model (MITgcm). The quality of the ocean-state estimates, forecasts, and the contribution of ensemble prediction are evaluated. The MITgcm–Estimating the Circulation and Climate of the Ocean (ECCO) 4DVAR (MITgcm–ECCO) and the MITgcm–Data Assimilation Research Testbed (DART) EnKF (MITgcm–DART) systems were used to compute two-month hindcasts (March–April, 2010) by assimilating satellite-derived along-track sea-surface height (SSH) and gridded sea-surface temperature (SST) observations. The estimates from both methods at the end of the hindcast period were then used to initialize forecasts for two months (May–June, 2010). This period was selected because a loop current (LC) eddy (Eddy Franklin: Eddy-F) detachment event occurred at the end of May 2010, immediately after the Deepwater Horizon (DwH) oil spill. Despite some differences between the setups, both systems produce analyses and forecasts of comparable quality and both solutions significantly outperformed model persistence. A reference forecast initialized from the 1/12° Hybrid Coordinate Ocean Model (HYCOM)/NCODA global analysis also performed well. The EnKF experiments for sensitivity to filter parameters showed enhanced predictability when using more ensemble members and stronger covariance localization, but not for larger inflation. The EnKF experiments varying the number of assimilation cycles showed enhanced short-term (long-term) predictability with fewer (more) assimilation cycles. Additional hindcast and forecast experiments at other times of significant LC evolution showed mixed performance of both systems, which depends strongly on the background state of the GoM circulation. The present work demonstrates a practical application of both assimilation methods for the GoM and compares them in a limited number of realizations. The overall conclusion showing improved short-term (long-term) predictability for EnKF (4DVAR) carries an important caveat that the results from this study are specific to a few 4DVAR and EnKF LC eddy separation experiments in the GoM and cannot be generalized to conclude the relative performance of both methods, especially in other applications. However, some of the concepts and methods should carry over to other applications.
dc.description.sponsorshipWe gratefully acknowledge the ECCO consortium, including MIT, JPL, and the University of Hamburg, and the NCAR Data Assimilation Research Section (DAReS). The MITgcm code used in this study is checkpoint 64Y, and was obtained from http://mitgcm.org/. The Ssalto/Duacs altimeter product AVISO is produced and distributed by the Copernicus Marine and Environment Monitoring Service (CMEMS: http://marine.copernicus.eu/). The HYCOM/NCODA 1/12° global analysis data were obtained from the HYCOM consortium (http://hycom.org/dataserver/). The along-track altimetry data were obtained from the Radar Altimeter Database System (RADS: http://rads.tudelft.nl/rads/index.shtml). The SST data were obtained from Remote Sensing Systems Inc. (http://www.remss.com/). The NCEP/NCAR-Reanalysis-1 atmospheric forcings were obtained from http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.html. The MITgcm–DART EnKF and MITgcm–ECCO 4DVAR state estimates, input files including observations and error fields, are available from the authors upon request (email:ggopalakrishnan@ucsd.edu). BDC gratefully acknowledges support from the Office of Naval Research grants N000141512285 and N000141512598. Research reported in this publication was supported by the Gulf Research Program of the National Academy of Sciences under award number 2000006422. The content is solely the responsibility of the authors and does not necessarily represent the official views of the Gulf Research Program or the National Academy of Sciences.
dc.publisherWiley
dc.relation.urlhttps://rmets.onlinelibrary.wiley.com/doi/full/10.1002/qj.3493
dc.rightsArchived with thanks to Quarterly Journal of the Royal Meteorological Society
dc.subjectData Assimilation
dc.subjectEnsemble Forecasting
dc.subjectGeneral Circulation Model
dc.subjectGulf of Mexico
dc.subjectLoop Current
dc.subjectLoop Current Eddy
dc.subjectMITgcm- DART EnKF
dc.subjectMITgcm-ECCO 4DVAR
dc.titleComparison of 4DVAR and EnKF state estimates and forecasts in the Gulf of Mexico
dc.typeArticle
dc.contributor.departmentEarth Fluid Modeling and Prediction Group
dc.contributor.departmentEarth Science and Engineering Program
dc.contributor.departmentPhysical Science and Engineering (PSE) Division
dc.identifier.journalQuarterly Journal of the Royal Meteorological Society
dc.eprint.versionPost-print
dc.contributor.institutionClimate, Atmospheric Science and Physical OceanographyScripps Institution of Oceanography La Jolla California
dc.contributor.institutionMarine Physical LaboratoryScripps Institution of Oceanography La Jolla California
kaust.personHoteit, Ibrahim
refterms.dateFOA2020-03-21T00:00:00Z


Files in this item

Thumbnail
Name:
gom_4dvar_enkf_v20181231.pdf
Size:
5.985Mb
Format:
PDF
Description:
Accepted Manuscript

This item appears in the following Collection(s)

Show simple item record