Editorial Type:
Article
Article Type:
Research Article

Atmospheric Data Assimilation with an Ensemble Kalman Filter: Results with Real Observations

P. L. Houtekamer Meteorological Service of Canada, Dorval, Québec, Canada

Search for other papers by P. L. Houtekamer in
Current site
Google Scholar
PubMed Close
,
Herschel L. Mitchell Meteorological Service of Canada, Dorval, Québec, Canada

Search for other papers by Herschel L. Mitchell in
Current site
Google Scholar
PubMed Close
,
Gérard Pellerin Meteorological Service of Canada, Dorval, Québec, Canada

Search for other papers by Gérard Pellerin in
Current site
Google Scholar
PubMed Close
,
Mark Buehner Meteorological Service of Canada, Dorval, Québec, Canada

Search for other papers by Mark Buehner in
Current site
Google Scholar
PubMed Close
,
Martin Charron Meteorological Service of Canada, Dorval, Québec, Canada

Search for other papers by Martin Charron in
Current site
Google Scholar
PubMed Close
,
Lubos Spacek Meteorological Service of Canada, Dorval, Québec, Canada

Search for other papers by Lubos Spacek in
Current site
Google Scholar
PubMed Close
, and
Bjarne Hansen Meteorological Service of Canada, Dorval, Québec, Canada

Search for other papers by Bjarne Hansen in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Mar 2005
DOI:
https://doi.org/10.1175/MWR-2864.1
Page(s):
604–620
Restricted access

Abstract

An ensemble Kalman filter (EnKF) has been implemented for atmospheric data assimilation. It assimilates observations from a fairly complete observational network with a forecast model that includes a standard operational set of physical parameterizations. To obtain reasonable results with a limited number of ensemble members, severe horizontal and vertical covariance localizations have been used.

It is observed that the error growth in the data assimilation cycle is mainly due to model error. An isotropic parameterization, similar to the forecast-error parameterization in variational algorithms, is used to represent model error. After some adjustment, it is possible to obtain innovation statistics that agree with the ensemble-based estimate of the innovation amplitudes for winds and temperature. Currently, no model error is added for the humidity variable, and, consequently, the ensemble spread for humidity is too small. After about 5 days of cycling, fairly stable global filter statistics are obtained with no sign of filter divergence.

The quality of the ensemble mean background field, as verified using radiosonde observations, is similar to that obtained using a 3D variational procedure. In part, this is likely due to the form chosen for the parameterized model error. Nevertheless, the degree of similarity is surprising given that the background-error statistics used by the two procedures are rather different, with generally larger background errors being used by the variational scheme.

A set of 5-day integrations has been started from the ensemble of initial conditions provided by the EnKF. For the middle and lower troposphere, the growth rates of the perturbations are somewhat smaller than the growth rate of the actual ensemble mean error. For the upper levels, the perturbation patterns decay for about 3 days as a consequence of diffusive model dynamics. These decaying perturbations tend to severely underestimate the actual error that grows rapidly near the model top.

Corresponding author address: P. L. Houtekamer, Direction de la Recherche en Météorologie, 2121 Route Trans-Canadienne, Dorval, Québec, H9P 1J3, Canada. Email: peter.houtekamer@ec.gc.ca

Abstract

An ensemble Kalman filter (EnKF) has been implemented for atmospheric data assimilation. It assimilates observations from a fairly complete observational network with a forecast model that includes a standard operational set of physical parameterizations. To obtain reasonable results with a limited number of ensemble members, severe horizontal and vertical covariance localizations have been used.

It is observed that the error growth in the data assimilation cycle is mainly due to model error. An isotropic parameterization, similar to the forecast-error parameterization in variational algorithms, is used to represent model error. After some adjustment, it is possible to obtain innovation statistics that agree with the ensemble-based estimate of the innovation amplitudes for winds and temperature. Currently, no model error is added for the humidity variable, and, consequently, the ensemble spread for humidity is too small. After about 5 days of cycling, fairly stable global filter statistics are obtained with no sign of filter divergence.

The quality of the ensemble mean background field, as verified using radiosonde observations, is similar to that obtained using a 3D variational procedure. In part, this is likely due to the form chosen for the parameterized model error. Nevertheless, the degree of similarity is surprising given that the background-error statistics used by the two procedures are rather different, with generally larger background errors being used by the variational scheme.

A set of 5-day integrations has been started from the ensemble of initial conditions provided by the EnKF. For the middle and lower troposphere, the growth rates of the perturbations are somewhat smaller than the growth rate of the actual ensemble mean error. For the upper levels, the perturbation patterns decay for about 3 days as a consequence of diffusive model dynamics. These decaying perturbations tend to severely underestimate the actual error that grows rapidly near the model top.

Corresponding author address: P. L. Houtekamer, Direction de la Recherche en Météorologie, 2121 Route Trans-Canadienne, Dorval, Québec, H9P 1J3, Canada. Email: peter.houtekamer@ec.gc.ca

Save
Cover Monthly Weather Review
Monthly Weather Review

  • Andersson, E., and H. Järvinen, 1999: Variational quality control. Quart. J. Roy. Meteor. Soc., 125 , 697–722.

  • Chouinard, C., C. Charette, J. Hallé, P. Gauthier, J. Morneau, and R. Sarrazin, 2001: The Canadian 3D-VAR analysis scheme on model vertical coordinate. Preprints, 14th Conf. on Numerical Weather Prediction, Fort Lauderdale, FL, Amer. Meteor. Soc., 14–18.

  • Chouinard, C., J. Hallé, C. Charette, and R. Sarrazin, 2002: Recent improvements in the use of TOVS satellite radiances in the unified 3D-VAR system of the Canadian Meteorological Centre. Proc. 12th Int. ATOVS Study Conf., Lorne, Australia, Bureau of Meteorology Research Centre, 38–44.

  • Côté, J., J-G. Desmarais, S. Gravel, A. Méthot, A. Patoine, M. Roch, and A. Staniforth, 1998a: The operational CMC–MRB Global Environmental Multiscale (GEM) model. Part II: Results. Mon. Wea. Rev., 126 , 1397–1418.

    • Search Google Scholar
    • Export Citation

  • Côté, J., S. Gravel, A. Méthot, A. Patoine, M. Roch, and A. Staniforth, 1998b: The operational CMC–MRB Global Environmental Multiscale (GEM) model. Part I: Design considerations and formulation. Mon. Wea. Rev., 126 , 1373–1395.

    • Search Google Scholar
    • Export Citation

  • Dee, D. P., 1995: On-line estimation of error covariance parameters for atmospheric data assimilation. Mon. Wea. Rev., 123 , 1128–1145.

    • Search Google Scholar
    • Export Citation

  • 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. 10143–10162.

    • Search Google Scholar
    • Export Citation

  • Evensen, G., 2003: The ensemble Kalman filter: Theoretical formulation and practical implementation. Ocean Dyn., 53 , 343–367.

  • Garand, L., 1993: A pattern recognition technique for retrieving humidity profiles from Meteosat or GOES imagery. J. Appl. Meteor., 32 , 1592–1607.

    • Search Google Scholar
    • Export Citation

  • Gaspari, G., and S. E. Cohn, 1999: Construction of correlation functions in two and three dimensions. Quart. J. Roy. Meteor. Soc., 125 , 723–757.

    • Search Google Scholar
    • Export Citation

  • Gauthier, P., M. Buehner, and L. Fillion, 1999a: Background-error statistics modelling in a 3D variational data assimilation scheme: Estimation and impact on the analyses. Proc. ECMWF Workshop on Diagnosis of Data Assimilation Systems, Reading, United Kingdom, ECMWF, 131–145. [Available from ECMWF, Shinfield Park, Reading, Berkshire RG2 9AX, United Kingdom.].

  • Gauthier, P., C. Charette, L. Fillion, P. Koclas, and S. Laroche, 1999b: Implementation of a 3d variational data assimilation system at the Canadian Meteorological Centre. Part I: The global analysis. Atmos.–Ocean, 37 , 103–156.

    • Search Google Scholar
    • Export Citation

  • Hamill, T. M., and C. Snyder, 2000: A hybrid ensemble Kalman filter–3D-variational analysis scheme. Mon. Wea. Rev., 128 , 2905–2919.

  • Houtekamer, P. L., and H. L. Mitchell, 1998: Data assimilation using an ensemble Kalman filter technique. Mon. Wea. Rev., 126 , 796–811.

    • Search Google Scholar
    • Export Citation

  • Houtekamer, P. L., and H. L. Mitchell, 2001: A sequential ensemble Kalman filter for atmospheric data assimilation. Mon. Wea. Rev., 129 , 123–137.

    • Search Google Scholar
    • Export Citation

  • Houtekamer, P. L., L. Lefaivre, J. Derome, H. Ritchie, and H. L. Mitchell, 1996: A system simulation approach to ensemble prediction. Mon. Wea. Rev., 124 , 1225–1242.

    • Search Google Scholar
    • Export Citation

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

    • Search Google Scholar
    • Export Citation

  • Lawson, W. G., and J. A. Hansen, 2004: Implications of stochastic and deterministic filters as ensemble-based data assimilation methods in varying regimes of error growth. Mon. Wea. Rev., 132 , 1966–1981.

    • Search Google Scholar
    • Export Citation

  • Lönnberg, P., and D. Shaw, 1987: ECMWF data assimilation scientific documentation. ECMWF, 96 pp. [Available from ECMWF, Shinfield Park, Reading, Berkshire RG2 9AX, United Kingdom.].

  • Lorenc, A. C., 2003: The potential of the ensemble Kalman filter for NWP—A comparison with 4D-Var. Quart. J. Roy. Meteor. Soc., 129 , 3183–3203.

    • Search Google Scholar
    • Export Citation

  • Mitchell, H. L., and P. L. Houtekamer, 2000: An adaptive ensemble Kalman filter. Mon. Wea. Rev., 128 , 416–433.

  • Mitchell, H. L., C. Chouinard, C. Charette, R. Hogue, and S. J. Lambert, 1996: Impact of a revised analysis algorithm on an operational data assimilation system. Mon. Wea. Rev., 124 , 1243–1255.

    • Search Google Scholar
    • Export Citation

  • Mitchell, H. L., P. L. Houtekamer, and G. Pellerin, 2002: Ensemble size, balance, and model-error representation in an ensemble Kalman filter. Mon. Wea. Rev., 130 , 2791–2808.

    • Search Google Scholar
    • Export Citation

  • Moghaddamjoo, A. R., and R. Kirlin, 1993: Robust adaptive Kalman filtering. Approximate Kalman Filtering, G. Chen, Ed., World Scientific, 65–85.

    • Search Google Scholar
    • Export Citation

  • Molteni, F., R. Buizza, T. N. Palmer, and T. Petroliagis, 1996: The ECMWF ensemble prediction system: Methodology and validation. Quart. J. Roy. Meteor. Soc., 122 , 73–119.

    • Search Google Scholar
    • Export Citation

  • Orrell, D., L. Smith, J. Barkmeijer, and T. N. Palmer, 2001: Model error in weather forecasting. Nonlinear Processes Geophys., 8 , 357–371.

    • Search Google Scholar
    • Export Citation

  • Pellerin, G., L. Lefaivre, P. Houtekamer, and C. Girard, 2003: Increasing the horizontal resolution of ensemble forecasts at CMC. Nonlinear Processes Geophys., 10 , 463–468.

    • Search Google Scholar
    • Export Citation

  • Sarrazin, R., and B. Brasnett, 2002: Modifications in the operational use of satellite atmospheric motion winds at CMC. Proc. Sixth Int. Winds Workshop, Madison, WI, EUMETSAT Publication 35, 147–154.

  • Saunders, R. W., 2000: RTTOV-6: Science and validation report. EUMETSAT, 31 pp. [Available from EUMETSAT Satellite Application Facility on NWP, The Met Office, London Road, Bracknell, Berkshire RG12 2SZ, United Kingdom.].

  • Saunders, R. W., M. Matricardi, and P. Brunel, 1999: An improved fast radiative transfer model for assimilation of satellite radiance observations. Quart. J. Roy. Meteor. Soc., 125 , 1407–1425.

    • Search Google Scholar
    • Export Citation

  • Tippett, M. K., J. L. Anderson, C. H. Bishop, T. M. Hamill, and J. S. Whitaker, 2003: Ensemble square root filters. Mon. Wea. Rev., 131 , 1485–1490.

    • Search Google Scholar
    • Export Citation

  • Toth, Z., and E. Kalnay, 1993: Ensemble forecasting at NMC: The generation of perturbations. Bull. Amer. Meteor. Soc., 74 , 2317–2330.

  • Whitaker, J. S., G. P. Compo, X. Wei, and T. M. Hamill, 2004: Reanalysis without radiosondes using ensemble data assimilation. Mon. Wea. Rev., 132 , 1190–1200.

    • Search Google Scholar
    • Export Citation

  • Zadra, A., M. Roch, S. Laroche, and M. Charron, 2003: The subgrid-scale orographic blocking parameterization of the GEM model. Atmos.–Ocean, 41 , 155–170.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 4572 2516 910
PDF Downloads 1154 302 52