A Simulation Study Using a Local Ensemble Transform Kalman Filter for Data Assimilation in New York Harbor

Ross N. Hoffman Atmospheric and Environmental Research, Inc., Lexington, Massachusetts

Search for other papers by Ross N. Hoffman in
Current site
Google Scholar
PubMed
Close
,
Rui M. Ponte Atmospheric and Environmental Research, Inc., Lexington, Massachusetts

Search for other papers by Rui M. Ponte in
Current site
Google Scholar
PubMed
Close
,
Eric J. Kostelich Department of Mathematics and Statistics, Arizona State University, Tempe, Arizona

Search for other papers by Eric J. Kostelich in
Current site
Google Scholar
PubMed
Close
,
Alan Blumberg Stevens Institute of Technology, Hoboken, New Jersey

Search for other papers by Alan Blumberg in
Current site
Google Scholar
PubMed
Close
,
Istvan Szunyogh University of Maryland, College Park, College Park, Maryland

Search for other papers by Istvan Szunyogh in
Current site
Google Scholar
PubMed
Close
,
Sergey V. Vinogradov Atmospheric and Environmental Research, Inc., Lexington, Massachusetts

Search for other papers by Sergey V. Vinogradov in
Current site
Google Scholar
PubMed
Close
, and
John M. Henderson Atmospheric and Environmental Research, Inc., Lexington, Massachusetts

Search for other papers by John M. Henderson in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

Data assimilation approaches that use ensembles to approximate a Kalman filter have many potential advantages for oceanographic applications. To explore the extent to which this holds, the Estuarine and Coastal Ocean Model (ECOM) is coupled with a modern data assimilation method based on the local ensemble transform Kalman filter (LETKF), and a series of simulation experiments is conducted. In these experiments, a long ECOM “nature” run is taken to be the “truth.” Observations are generated at analysis times by perturbing the nature run at randomly chosen model grid points with errors of known statistics. A diverse collection of model states is used for the initial ensemble. All experiments use the same lateral boundary conditions and external forcing fields as in the nature run. In the data assimilation, the analysis step combines the observations and the ECOM forecasts using the Kalman filter equations. As a control, a free-running forecast (FRF) is made from the initial ensemble mean to check the relative importance of external forcing versus data assimilation on the analysis skill. Results of the assimilation cycle and the FRF are compared to truth to quantify the skill of each.

The LETKF performs well for the cases studied here. After just a few assimilation cycles, the analysis errors are smaller than the observation errors and are much smaller than the errors in the FRF. The assimilation quickly eliminates the domain-averaged bias of the initial ensemble. The filter accurately tracks the truth at all data densities examined, from observations at 50% of the model grid points down to 2% of the model grid points. As the data density increases, the ensemble spread, bias, and error standard deviation decrease. As the ensemble size increases, the ensemble spread increases and the error standard deviation decreases. Increases in the size of the observation error lead to a larger ensemble spread but have a small impact on the analysis accuracy.

Corresponding author address: Dr. Ross N. Hoffman, Atmospheric and Environmental Research, Inc., 131 Hartwell Avenue, Lexington, MA 02421–3126. Email: ross.n.hoffman@aer.com

Abstract

Data assimilation approaches that use ensembles to approximate a Kalman filter have many potential advantages for oceanographic applications. To explore the extent to which this holds, the Estuarine and Coastal Ocean Model (ECOM) is coupled with a modern data assimilation method based on the local ensemble transform Kalman filter (LETKF), and a series of simulation experiments is conducted. In these experiments, a long ECOM “nature” run is taken to be the “truth.” Observations are generated at analysis times by perturbing the nature run at randomly chosen model grid points with errors of known statistics. A diverse collection of model states is used for the initial ensemble. All experiments use the same lateral boundary conditions and external forcing fields as in the nature run. In the data assimilation, the analysis step combines the observations and the ECOM forecasts using the Kalman filter equations. As a control, a free-running forecast (FRF) is made from the initial ensemble mean to check the relative importance of external forcing versus data assimilation on the analysis skill. Results of the assimilation cycle and the FRF are compared to truth to quantify the skill of each.

The LETKF performs well for the cases studied here. After just a few assimilation cycles, the analysis errors are smaller than the observation errors and are much smaller than the errors in the FRF. The assimilation quickly eliminates the domain-averaged bias of the initial ensemble. The filter accurately tracks the truth at all data densities examined, from observations at 50% of the model grid points down to 2% of the model grid points. As the data density increases, the ensemble spread, bias, and error standard deviation decrease. As the ensemble size increases, the ensemble spread increases and the error standard deviation decreases. Increases in the size of the observation error lead to a larger ensemble spread but have a small impact on the analysis accuracy.

Corresponding author address: Dr. Ross N. Hoffman, Atmospheric and Environmental Research, Inc., 131 Hartwell Avenue, Lexington, MA 02421–3126. Email: ross.n.hoffman@aer.com

Save
  • Ahsan, Q., and Blumberg A. , 1999: Three-dimensional hydrothermal model of Onondaga Lake, New York. J. Hydrol. Eng., 125 , 912923.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Allen, J. S., Newberger P. A. , and Federiuk J. , 1995: Upwelling circulation on the Oregon continental shelf. Part I: Response to idealized forcing. J. Phys. Oceanogr., 25 , 18431866.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Baek, S-J., Hunt B. R. , Kalnay E. , Ott E. , and Szunyogh I. , 2006: Local ensemble Kalman filtering in the presence of model bias. Tellus, 58A , 293306.

    • Search Google Scholar
    • Export Citation
  • Bishop, C. H., Etherton B. J. , and Majumdar S. J. , 2001: Adaptive sampling with the ensemble transform Kalman filter. Part I: Theoretical aspects. Mon. Wea. Rev., 129 , 420436.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Blumberg, A. F., 1996: An estuarine and coastal ocean version of POM. Proc. Princeton Ocean Model Users Meeting (POM96), Princeton, NJ, Princeton University, 9.

  • Blumberg, A. F., and Mellor G. L. , 1987: A description of a three-dimensional coastal ocean circulation model. Three-Dimensional Coastal Ocean Models, N. Heaps, Ed., American Geophysical Union, 1–16.

    • Search Google Scholar
    • Export Citation
  • Blumberg, A. F., and Goodrich D. M. , 1990: Modeling of wind-induced destratification in Chesapeake Bay. Estuaries Coasts, 13 , 236249.

  • Blumberg, A. F., Khan L. A. , and St. John J. P. , 1999: Three-dimensional hydrodynamic simulations of the New York Harbor, Long Island Sound, and the New York Bight. J. Hydraul. Eng., 125 , 799816.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Blumberg, A. F., Ahsan Q. , and Lewis J. K. , 2000: Modeling hydrodynamics of the Mississippi Sound and adjoining rivers, bays, and shelf waters. Proc. OCEANS 2000 MTS/IEEE Conf. and Exhibition, Vol. 3, Providence, RI, IEEE, 1983–1989.

  • Bruno, M. S., Blumberg A. F. , and Herrington T. O. , 2006: The urban ocean observatory—Coastal ocean observations and forecasting in the New York Bight. J. Mar. Sci. Environ., C4 , 19.

    • Search Google Scholar
    • Export Citation
  • Courtier, P., 1997: Variational methods. J. Meteor. Soc. Japan, 75 , 211218.

  • Daley, R., 1991: Atmospheric Data Analysis. Cambridge University Press, 457 pp.

  • Evensen, G., 1994: Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics. J. Geophys. Res., 99 , C5. 1014310162.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Evensen, G., 2006: Data Assimilation: The Ensemble Kalman Filter. Springer, 280 pp.

  • Fan, S., Oey L-Y. , and Hamilton P. , 2004: Assimilation of drifter and satellite data in a model of the northeastern Gulf of Mexico. Cont. Shelf Res., 24 , 10011013.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Fukumori, I., 2002: A partitioned Kalman filter and smoother. Mon. Wea. Rev., 130 , 13701383.

  • Fukumori, I., and Malanotte-Rizzoli P. , 1995: An approximate Kalman filter for ocean data assimilation: An example with an idealized Gulf Stream model. J. Geophys. Res., 100 , C4. 67776793.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Galperin, B., and Mellor G. L. , 1990: A time-dependent, three-dimensional model of the Delaware Bay and River system. Part 1: Description of the model and tidal analysis. Estuarine Coastal Shelf Sci., 31 , 231253.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hunt, B., Kostelich E. J. , and Szunyogh I. , 2007: Efficient data assimilation for spatiotemporal chaos: A local ensemble transform Kalman filter. Physica D, 230 , 112126.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Jazwinski, A. H., 1970: Stochastic Processes and Filtering Theory. Academic Press, 376 pp.

  • Kalnay, E., 2002: Atmospheric Modeling, Data Assimilation, and Predictability. Cambridge University Press, 364 pp.

  • Keppenne, C., and Rienecker M. , 2002: Initial testing of a massively parallel ensemble Kalman filter with the Poseidon isopycnal ocean general circulation model. Mon. Wea. Rev., 130 , 29512965.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lermusiaux, P. F. J., and Coauthors, 2006: Quantifying uncertainties in ocean predictions. Oceanography, 19 , 90103.

  • Lin, S-J., and Rood R. B. , 1996: Multidimensional flux-form semi-Lagrangian transport schemes. Mon. Wea. Rev., 124 , 20462070.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Lorenc, A. C., 1997: Development of an operational variational assimilation scheme. J. Meteor. Soc. Japan, 75 , 339346.

  • Mellor, G. L., and Yamada T. , 1982: Development of a turbulence closure model for geophysical fluid problems. Rev. Geophys. Space Phys., 20 , 851875.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Mellor, G. L., and Ezer T. , 1991: A Gulf Stream model and an altimetry assimilation scheme. J. Geophys. Res., 96 , 87798795.

  • Ott, E., and Coauthors, 2004: A local ensemble Kalman filter for atmospheric data assimilation. Tellus, 56A , 415428.

  • Phillips, N. A., 1957: A coordinate system having some special advantages for numerical forecasting. J. Meteor., 14 , 184185.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Signell, R. P., Jenter H. L. , and Blumberg A. F. , 1994: Modeling the seasonal circulation in Massachusetts Bay. Proc. Third Int. Conf. on Estuarine and Coastal Modeling, ASCE, 578–590.

    • Search Google Scholar
    • Export Citation
  • Stammer, D., and Coauthors, 2002: Global ocean circulation during 1992–1997, estimated from ocean observations and a general circulation model. J. Geophys. Res., 107 .3118, doi:10.1029/2001JC000888.

    • Search Google Scholar
    • Export Citation
  • Szunyogh, I., Kostelich E. J. , Gyarmati G. , Patil D. J. , Hunt B. R. , Kalnay E. , Ott E. , and Yorke J. A. , 2005: Assessing a local ensemble Kalman filter: Perfect model experiments with the National Centers for Environmental Prediction global model. Tellus, 57A , 528545.

    • Search Google Scholar
    • Export Citation
  • Szunyogh, I., Kostelich E. J. , Gyarmati G. , Kalnay E. , Hunt B. R. , Ott E. , Satterfield E. , and Yorke J. A. , 2008: A local ensemble transform Kalman filter data assimilation system for the NCEP global model. Tellus, 60A , 113130. doi:10.1111/j.1600-0870.2007.00274.x.

    • Search Google Scholar
    • Export Citation
  • Vinogradova, N., Vinogradov S. , Nechaev D. , Kamenkovich V. , Blumberg A. F. , Ahsan Q. , and Li H. , 2005: Evaluation of the northern Gulf of Mexico littoral initiative (NGLI) model based on the observed temperature and salinity in the Mississippi Bight. Mar. Tech. Soc. J., 39 , 2538.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Whitaker, J. S., Compo G. P. , Wei X. , and Hamill T. M. , 2004: Reanalysis without radiosondes using ensemble data assimilation. Mon. Wea. Rev., 132 , 11901200.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Whitaker, J. S., Hamill T. M. , Wei X. , Song Y. , and Toth Z. , 2007: Ensemble data assimilation with the NCEP global forecast system. Mon. Wea. Rev., 136 , 463482.

    • Search Google Scholar
    • Export Citation
  • Wilks, D. S., 2006: Statistical Methods in the Atmospheric Sciences. 2nd ed. Academic Press, 648 pp.

  • Wunsch, C., and Heimbach P. , 2007: Practical global ocean state estimation. Physica D, 230 , 197208. doi:10.1016/j.physd.2006.09.040.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 307 121 4
PDF Downloads 44 27 4