Extension of 3DVAR to 4DVAR: Implementation of 4DVAR at the Meteorological Service of Canada

Pierre Gauthier Meteorological Research Division, Environment Canada, Dorval, Québec, Canada

Search for other papers by Pierre Gauthier in
Current site
Google Scholar
PubMed
Close
,
Monique Tanguay Meteorological Research Division, Environment Canada, Dorval, Québec, Canada

Search for other papers by Monique Tanguay in
Current site
Google Scholar
PubMed
Close
,
Stéphane Laroche Meteorological Research Division, Environment Canada, Dorval, Québec, Canada

Search for other papers by Stéphane Laroche in
Current site
Google Scholar
PubMed
Close
,
Simon Pellerin Meteorological Research Division, Environment Canada, Dorval, Québec, Canada

Search for other papers by Simon Pellerin in
Current site
Google Scholar
PubMed
Close
, and
Josée Morneau Meteorological Service of Canada, Environment Canada, Dorval, Québec, Canada

Search for other papers by Josée Morneau in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

On 15 March 2005, the Meteorological Service of Canada (MSC) proceeded to the implementation of a four-dimensional variational data assimilation (4DVAR) system, which led to significant improvements in the quality of global forecasts. This paper describes the different elements of MSC’s 4DVAR assimilation system, discusses some issues encountered during the development, and reports on the overall results from the 4DVAR implementation tests. The 4DVAR system adopted an incremental approach with two outer iterations. The simplified model used in the minimization has a horizontal resolution of 170 km and its simplified physics includes vertical diffusion, surface drag, orographic blocking, stratiform condensation, and convection. One important element of the design is its modularity, which has permitted continued progress on the three-dimensional variational data assimilation (3DVAR) component (e.g., addition of new observation types) and the model (e.g., computational and numerical changes). This paper discusses some numerical problems that occur in the vicinity of the Poles where the semi-Lagrangian scheme becomes unstable when there is a simultaneous occurrence of converging meridians and strong wind gradients. These could be removed by filtering the winds in the zonal direction before they are used to estimate the upstream position in the semi-Lagrangian scheme. The results show improvements in all aspects of the forecasts over all regions. The impact is particularly significant in the Southern Hemisphere where 4DVAR is able to extract more information from satellite data. In the Northern Hemisphere, 4DVAR accepts more asynoptic data, in particular coming from profilers and aircrafts. The impact noted is also positive and the short-term forecasts are particularly improved over the west coast of North America. Finally, the dynamical consistency of the 4DVAR global analyses leads to a significant impact on regional forecasts. Experimentation has shown that regional forecasts initiated directly from a 4DVAR global analysis are improved with respect to the regional forecasts resulting from the regional 3DVAR analysis.

Corresponding author address: Dr. Pierre Gauthier, Data Assimilation and Satellite Meteorology Section, Meteorological Research Division, Environment Canada, 2121 Trans-Canada Highway, Dorval, QC H9P 1J3, Canada. Email: pierre.gauthier@ec.gc.ca

Abstract

On 15 March 2005, the Meteorological Service of Canada (MSC) proceeded to the implementation of a four-dimensional variational data assimilation (4DVAR) system, which led to significant improvements in the quality of global forecasts. This paper describes the different elements of MSC’s 4DVAR assimilation system, discusses some issues encountered during the development, and reports on the overall results from the 4DVAR implementation tests. The 4DVAR system adopted an incremental approach with two outer iterations. The simplified model used in the minimization has a horizontal resolution of 170 km and its simplified physics includes vertical diffusion, surface drag, orographic blocking, stratiform condensation, and convection. One important element of the design is its modularity, which has permitted continued progress on the three-dimensional variational data assimilation (3DVAR) component (e.g., addition of new observation types) and the model (e.g., computational and numerical changes). This paper discusses some numerical problems that occur in the vicinity of the Poles where the semi-Lagrangian scheme becomes unstable when there is a simultaneous occurrence of converging meridians and strong wind gradients. These could be removed by filtering the winds in the zonal direction before they are used to estimate the upstream position in the semi-Lagrangian scheme. The results show improvements in all aspects of the forecasts over all regions. The impact is particularly significant in the Southern Hemisphere where 4DVAR is able to extract more information from satellite data. In the Northern Hemisphere, 4DVAR accepts more asynoptic data, in particular coming from profilers and aircrafts. The impact noted is also positive and the short-term forecasts are particularly improved over the west coast of North America. Finally, the dynamical consistency of the 4DVAR global analyses leads to a significant impact on regional forecasts. Experimentation has shown that regional forecasts initiated directly from a 4DVAR global analysis are improved with respect to the regional forecasts resulting from the regional 3DVAR analysis.

Corresponding author address: Dr. Pierre Gauthier, Data Assimilation and Satellite Meteorology Section, Meteorological Research Division, Environment Canada, 2121 Trans-Canada Highway, Dorval, QC H9P 1J3, Canada. Email: pierre.gauthier@ec.gc.ca

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

  • Aparicio, J., and G. Deblonde, 2007: Evaluation of GPS radio occultation assimilation. Mon. Wea. Rev., in press.

  • Bélair, S., J. Mailhot, C. Girard, and P. Vaillancourt, 2005: Boundary layer and shallow cumulus clouds in a medium-range forecast of a large-scale weather system. Mon. Wea. Rev., 133 , 19381960.

    • Search Google Scholar
    • Export Citation
  • 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.

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

    • Search Google Scholar
    • Export Citation
  • Côté, J., J. G. Desmarais, S. Gravel, A. Méthot, A. Patoine, M. Roch, and A. N. Staniforth, 1998b: The operational CMC-MRB Global Environmental Multiscale (GEM) model. Part II: Results. Mon. Wea. Rev., 126 , 13971418.

    • Search Google Scholar
    • Export Citation
  • Courtier, P., 1997: Dual formulation of four-dimensional variational data assimilation. Quart. J. Roy. Meteor. Soc., 123 , 24492461.

  • Courtier, P., J. N. Thépaut, and A. Hollingsworth, 1994: A strategy for operational implementation of 4D-Var, using an incremental approach. Quart. J. Roy. Meteor. Soc., 120 , 13671387.

    • Search Google Scholar
    • Export Citation
  • Deblonde, G., 1999: Variational assimilation of SSM/I total precipitable water retrievals in the CMC analysis system. Mon. Wea. Rev., 127 , 14581476.

    • Search Google Scholar
    • Export Citation
  • Desgagné, M., S. Thomas, and M. Valin, 2000: Performance of MC2 and the ECMWF IFS forecast model on the Fujitsu VPP700 and NEC SX 4M. Sci. Program., 8 , 2330.

    • Search Google Scholar
    • Export Citation
  • Durran, D. R., 1999: Numerical Methods for Wave Equations in Geophysical Fluid Dynamics. Springer Verlag, 465 pp.

  • Eyre, J. R., 1989: Inversion of cloudy satellite sounding radiances by nonlinear optimal estimation. I: Theory and simulation of TOVS. Quart. J. Roy. Meteor. Soc., 115 , 10011026.

    • Search Google Scholar
    • Export Citation
  • Fillion, L., H. L. Mitchell, H. R. Ritchie, and A. N. Staniforth, 1995: The impact of a digital filter finalization technique in a global data assimilation system. Tellus, 47A , 304323.

    • Search Google Scholar
    • Export Citation
  • Garand, L., A. Beaulne, and N. Wagneur, 2005: Assimilation of AIRS hyperspectral radiances at the Meteorological Service of Canada. Preprints, 14th Conf. on Satellite Meteorology and Oceanography, Atlanta, GA, Amer. Meteor. Soc., CD-ROM, P5.13.

  • Gauthier, P., 2003: Operational implementation of variational data assimilation. Data Assimilation for the Earth System, R. Swinbank, V. Shutyaev, and W. A. Lahoz, Eds., NATO Science Series, IV: Earth and Environmental Sciences, Vol. 26, Kluwer Academic, 167–176.

  • Gauthier, P., and J-N. Thépaut, 2001: Impact of the digital filter as a weak constraint in the preoperational 4DVAR assimilation system of Météo-France. Mon. Wea. Rev., 129 , 20892102.

    • Search Google Scholar
    • Export Citation
  • Gauthier, P., C. Charette, L. Fillion, P. Koclas, and S. Laroche, 1999a: Implementation of a 3D variational data assimilation system at the Canadian Meteorological Centre. Part I: The global analysis. Atmos.–Ocean, 37 , 103156.

    • Search Google Scholar
    • Export Citation
  • Gauthier, P., M. Buehner, and L. Fillion, 1999b: 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.

  • Gauthier, P., C. Chouinard, and B. Brasnett, 2003: Quality control: Methodology and applications. Data Assimilation for the Earth System, R. Swinbank, V. Shutyaev, and W. A. Lahoz, Eds., NATO Science Series, IV: Earth and Environmental Sciences, Vol. 26, Kluwer Academic, 177–187.

  • Gilbert, J. C., and C. LeMaréchal, 1989: Some numerical experiments with variable-storage quasi-Newton algorithms. Math. Program., 45 , 407435.

    • Search Google Scholar
    • Export Citation
  • Ingleby, N. B., and A. C. Lorenc, 1993: Bayesian quality control using multivariate normal distributions. Quart. J. Roy. Meteor. Soc., 119 , 11951225.

    • Search Google Scholar
    • Export Citation
  • Lagarde, T., A. Piacentini, and O. Thual, 2001: A new representation of data assimilation methods: The PALM flow-charting approach. Quart. J. Roy. Meteor. Soc., 127 , 189207.

    • Search Google Scholar
    • Export Citation
  • Laroche, S., and P. Gauthier, 1998: A validation of the incremental formulation of 4D variational data assimilation in a nonlinear barotropic flow. Tellus, 50A , 557572.

    • Search Google Scholar
    • Export Citation
  • Laroche, S., P. Gauthier, J. St. James, and J. Morneau, 1999: Implementation of a 3D variational data assimilation system at the Canadian Meteorological Centre. Part II: The regional analysis. Atmos.– Ocean, 37 , 281307.

    • Search Google Scholar
    • Export Citation
  • Laroche, S., M. Tanguay, and Y. Delage, 2002a: Linearization of a simplified planetary boundary layer parameterization. Mon. Wea. Rev., 130 , 20742087.

    • Search Google Scholar
    • Export Citation
  • Laroche, S., M. Tanguay, A. Zadra, and J. Morneau, 2002b: Use of adjoint sensitivity analysis to diagnose the CMC global analysis performance: A case study. Atmos.–Ocean, 40 , 423443.

    • Search Google Scholar
    • Export Citation
  • Laroche, S., P. Gauthier, M. Tanguay, S. Pellerin, and J. Morneau, 2007: Impact of the different components of 4DVAR on the global forecast system of the Meteorological Service of Canada. Mon. Wea. Rev., 135 , 23552364.

    • Search Google Scholar
    • Export Citation
  • Liu, Z. Q., and F. Rabier, 2003: The potential of high-density observations for numerical weather prediction: A study with simulated observations. Quart. J. Roy. Meteor. Soc., 129 , 30133035.

    • Search Google Scholar
    • Export Citation
  • Lorenc, A. C., 1986: Analysis methods for numerical weather prediction. Quart. J. Roy. Meteor. Soc., 112 , 11771194.

  • Mahfouf, J-F., 1999: Influence of physical processes on the tangent-linear approximation. Tellus, 51A , 147166.

  • Mahfouf, J-F., 2005: Linearization of a simple moist convection scheme for large-scale NWP models. Mon. Wea. Rev., 133 , 16551670.

  • Ménard, R., 2005: Advances in atmospheric chemistry and dynamics research by development of coupled chemistry-dynamics data assimilation. ESTEC Contract 18560/04/NL/FF, 116 pp.

  • Pellerin, S., P. Gauthier, M. Tanguay, S. Laroche, and J. Morneau, 2005: A modular implementation of the operational 4D-Var data assimilation at the Meteorological Service of Canada. Proc. Fourth WMO Symp. on Assimilation of Observations in Meteorology and Oceanography, Prague, Czech Republic, WMO, WWRP 9, WMO-TD 1316, CD-ROM.

  • Polavarapu, S., and M. Tanguay, 1998: Linearizing iterative processes for four-dimensional data assimilation schemes. Quart. J. Roy. Meteor. Soc., 124 , 17151742.

    • Search Google Scholar
    • Export Citation
  • Polavarapu, S., M. Tanguay, R. Ménard, and A. Staniforth, 1996: The tangent linear model for semi-Lagrangian schemes: Linearizing the process of interpolation. Tellus, 48A , 7495.

    • Search Google Scholar
    • Export Citation
  • Polavarapu, S., S. Ren, Y. Rochon, D. Sankey, N. Ek, J. Koshyk, and D. Tarasick, 2005: Data assimilation with the Canadian Middle Atmosphere Model. Atmos.–Ocean, 44 , 77100.

    • Search Google Scholar
    • Export Citation
  • Pudykiewicz, J., and A. Staniforth, 1984: Some properties and comparative performance of the semi-Lagrangian method of Robert in the solution of the advection-diffusion equation. Atmos.–Ocean, 22 , 283308.

    • Search Google Scholar
    • Export Citation
  • Qaddouri, A., J. Côté, and M. Valin, 1999: A parallel direct 3D elliptic solver. High Performance Computing Systems and Applications, A. Pollard, D. Mewhort, and D. Weaver, Eds., Kluwer Academic, 429–442.

    • Search Google Scholar
    • Export Citation
  • Rabier, F., H. Järvinen, E. Klinker, J-F. Mahfouf, and A. Simmons, 2000: The ECMWF operational implementation of four dimensional variational assimilation. Part I: Experimental results with simplified physics. Quart. J. Roy. Meteor. Soc., 126 , 11431170.

    • Search Google Scholar
    • Export Citation
  • Smolarkiewicz, P., and J. Pudykiewicz, 1992: A class of semi-Lagrangian approximations for fluids. J. Atmos. Sci., 49 , 20822096.

  • Talagrand, O., and P. Courtier, 1987: Variational assimilation of meteorological observations with the adjoint vorticity equation. I—Theory. Quart. J. Roy. Meteor. Soc., 113 , 13111328.

    • Search Google Scholar
    • Export Citation
  • Tanguay, M., and S. Polavarapu, 1999: The adjoint of the semi-Lagrangian treatment of the passive tracer equation. Mon. Wea. Rev., 127 , 551564.

    • Search Google Scholar
    • Export Citation
  • Tanguay, M., S. Polavarapu, and P. Gauthier, 1997: Temporal accumulation of first-order linearization error for semi-Lagrangian passive advection. Mon. Wea. Rev., 125 , 12961311.

    • Search Google Scholar
    • Export Citation
  • Trémolet, Y., 2004: Diagnostics of linear and incremental approximations in 4D-Var. Quart. J. Roy. Meteor. Soc., 130 , 22332251.

  • Zadra, A., M. Buehner, S. Laroche, and J-F. Mahfouf, 2004: Impact of the GEM model simplified physics on the extratropical singular vectors. Quart. J. Roy. Meteor. Soc., 130 , 25412569.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 921 260 22
PDF Downloads 565 135 17