Editorial Type:
Article
Article Type:
Research Article

Four-Dimensional Variational Data Assimilation with Digital Filter Initialization

Saroja Polavarapu Meteorological Service of Canada, Downsview, Ontario, Canada

Search for other papers by Saroja Polavarapu in
Current site
Google Scholar
PubMed Close
,
Monique Tanguay Meteorological Service of Canada, Downsview, Ontario, Canada

Search for other papers by Monique Tanguay in
Current site
Google Scholar
PubMed Close
, and
Luc Fillion Meteorological Service of Canada, Downsview, Ontario, Canada

Search for other papers by Luc Fillion in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Jul 2000
DOI:
https://doi.org/10.1175/1520-0493(2000)128<2491:FDVDAW>2.0.CO;2
Page(s):
2491–2510
Restricted access

Abstract

A four-dimensional variational (4DVAR) data assimilation problem may be constrained so that the solution closely fits the observations but is balanced. In this way, the processes of data analysis and initialization are combined. The method of initialization considered here, digital filtering, is widely used in weather forecasting centers. The digital filter was found to control high-frequency noise when implemented as a strong or as a weak constraint in the context of a global shallow water model. Implementation of a strong constraint did not result in a recovery of small scales although some recovery of intermediate scales did occur. Implementation of a weak constraint as a penalty method with a single fixed value of the penalty parameter resulted in analyses that were smooth, but depended upon the choice of the parameter. With a parameter value that was too large, the divergent kinetic energy spectrum of the analysis was excessively damped in the large scales. The rotational kinetic energy spectrum was also affected by the choice of penalty parameter. Both types of constraint were found to adequately control gravity wave noise although caution is advised in choosing the penalty parameter for the simple penalty term method.

Corresponding author address: Saroja Polavarapu, Data Assimilation and Satellite Meteorology Division, Meteorological Service of Canada, 4905 Dufferin Street, Downsview, ON M3H 4T4, Canada.

Abstract

A four-dimensional variational (4DVAR) data assimilation problem may be constrained so that the solution closely fits the observations but is balanced. In this way, the processes of data analysis and initialization are combined. The method of initialization considered here, digital filtering, is widely used in weather forecasting centers. The digital filter was found to control high-frequency noise when implemented as a strong or as a weak constraint in the context of a global shallow water model. Implementation of a strong constraint did not result in a recovery of small scales although some recovery of intermediate scales did occur. Implementation of a weak constraint as a penalty method with a single fixed value of the penalty parameter resulted in analyses that were smooth, but depended upon the choice of the parameter. With a parameter value that was too large, the divergent kinetic energy spectrum of the analysis was excessively damped in the large scales. The rotational kinetic energy spectrum was also affected by the choice of penalty parameter. Both types of constraint were found to adequately control gravity wave noise although caution is advised in choosing the penalty parameter for the simple penalty term method.

Corresponding author address: Saroja Polavarapu, Data Assimilation and Satellite Meteorology Division, Meteorological Service of Canada, 4905 Dufferin Street, Downsview, ON M3H 4T4, Canada.

Save
Cover Monthly Weather Review
Monthly Weather Review

  • 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.

  • ——, 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.

  • Courtier, P., and O. Talagrand, 1990: Variational assimilation of meteorological observations with the direct and adjoint shallow-water equations. Tellus,42A, 531–549.

  • ——, 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, 1367–1387.

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

  • Derber, J., D. F. Parrish, W. S. Wu, Z. Pu, and S. R. H. Rizvi, 1994:Improvements to the operational SSI global analysis system. Preprints, 10th Conf. on Numerical Weather Prediction, Portland, OR, Amer. Meteor. Soc., 149–150.

  • Errico, R. M., 2000: The dynamical balance of singular vectors in a primitive equation model. Quart. J. Roy. Meteor. Soc., in press.

  • ——, and K. D. Raeder, 1999: An examination of the accuracy of the linearization of a mesoscale model with moist physics. Quart. J. Roy. Meteor. Soc.,125, 169–195.

  • Fillion, L., H. L. Mitchell, H. Ritchie, and A. Staniforth, 1995: The impact of a digital filter finalization technique in a global data assimilation system. Tellus,47A, 304–323.

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

  • Gill, P. E., W. Murray, and M. H. Wright, 1981: Practical Optimization. Academic Press, 401 pp.

  • Gustafsson, N., 1993: Use of a digital filter as weak constraint in variational data assimilation. Proc. Workshop on Variational Assimilation, with Special Emphasis on Three-Dimensional Aspects, Reading, United Kingdom, ECMWF, 327–338.

  • Huang, X.-Y., and P. Lynch, 1993: Diabatic digital-filtering initialization: Application to the HIRLAM model. Mon. Wea. Rev.,121, 589–603.

  • ——, and H. Sundqvist, 1993: Initialization of cloud water content and cloud cover for numerical prediction models. Mon. Wea. Rev.,121, 2719–2726.

  • ——, A. Cederskov, and E. Källén, 1994: A comparison between digital filtering initialization and nonlinear normal-mode initialization in a data assimilation system. Mon. Wea. Rev.,122, 1001–1015.

  • Janiskova, M., J. N. Thepaut, and J. F. Geleyn, 1999: Simplified and regular physical parametrizations for incremental four-dimensional variational assimilation. Mon. Wea. Rev.,127, 26–45.

  • Laroche, S., and P. Gauthier, 1998: A validation of the incremental formulation of 4D variational data assimilation in a nonlinear barotropic flow. Tellus,50A, 557–572.

  • Lynch, P., and X.-Y. Huang, 1992: Initialization of the HIRLAM model using a digital filter. Mon. Wea. Rev.,120, 1019–1034.

  • ——, and ——, 1994: Diabatic initialization using recursive filters. Tellus,46A, 583–597.

  • ——, D. Giard, and V. Ivanovici, 1997: Improving the efficiency of a digital filtering scheme for diabatic initialization. Mon. Wea. Rev.,125, 1976–1982.

  • ——, D. Majewski, and P. Prohl, 1998: Digital filtering initialization of the global icosahedral-hexagonal grid NWP model (GME) of the Deutscher Wetterdienst. CAS/JSC Working Group on Numerical Experimentation Rep. 27 on Research Activities in Atmospheric and Oceanic Modelling, A. Staniforth, Ed., WMO, 459 pp.

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

  • Mitchell, H. L., C. Charette, C. Chouinard, and B. Brasnett, 1990: Revised interpolation statistics for the Canadian data assimilation procedure: Their derivation and application. Mon. Wea. Rev.,118, 1591–1614.

  • Moorthi, S., R. W. Higgins, and J. R. Bates, 1995: A global multilevel atmospheric model using a vector semi-Lagrangian finite-difference scheme. Part II: Version with physics. Mon. Wea. Rev.,123, 1523–1541.

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

  • Puri, K., G. S. Dietachmayer, G. A. Mills, N. E. Davidson, R. A. Bowen, and L. W. Logan, 1998: The new BMRC Limited Area Prediction System, LAPS. Aust. Meteor. Mag.,47, 203–223.

  • Sasaki, Y., 1970: Some basic formalisms in numerical variational analysis. Mon. Wea. Rev.,98, 875–883.

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

  • Tanguay, M., P. Bartello, and P. Gauthier, 1995: Four-dimensional data assimilation with a wide range of scales. Tellus,47A, 974–997.

  • Thépaut, J.-N., and P. Courtier, 1991: Four-dimensional data assimilation using the adjoint of a multilevel primitive equation model. Quart. J. Roy. Meteor. Soc.,117, 1225–1254.

  • ——, D. Vasiljevic, P. Courtier, and J. Pailleux, 1993: Variational assimilation of conventional meteorological observations with a multilevel primitive-equation model. Quart. J. Roy. Meteor. Soc.,119, 153–186.

  • Zou, X., I. M. Navon, and F. X. Le Dimet, 1992: Incomplete observations and control of gravity waves in variational data assimilation. Tellus,44A, 273–296.

  • ——, ——, and J. Sela, 1993: Control of gravitational oscillations in variational data assimilation. Mon. Wea. Rev.,121, 272–289.

  • Zupanski, M., and D. Zupanski, 1995: Recent developments of the NMC’s regional four-dimensional variational data assimilation system. Proc. Second Int. Symp. on Assimilation of Observations in Meteorology and Oceanography, Tokyo, Japan, WMO/TD-No. 651, 367–372.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 752 571 17
PDF Downloads 171 47 1