Editorial Type:
Article
Article Type:
Research Article

Four-Dimensional Variational Data Assimilation for Limited-Area Models: Lateral Boundary Conditions, Solution Uniqueness, and Numerical Convergence

Chungu Lu NOAA/OAR Forecast Systems Laboratory, Boulder, Colorado, and Cooperative Institute for Research in the Atmosphere, Colorado State University, Fort Collins, Colorado

Search for other papers by Chungu Lu in
Current site
Google Scholar
PubMed Close
and
Gerald L. Browning NOAA/OAR Forecast Systems Laboratory, Boulder, Colorado, and Cooperative Institute for Research in the Atmosphere, Colorado State University, Fort Collins, Colorado

Search for other papers by Gerald L. Browning in
Current site
Google Scholar
PubMed Close
Print Publication:
01 May 2000
DOI:
https://doi.org/10.1175/1520-0469(2000)057<1341:FDVDAF>2.0.CO;2
Page(s):
1341–1353
Restricted access

Abstract

Mathematical issues arising when applying four-dimensional variational (4DVAR) data assimilation to limited-area problems are studied. The derivation of the adjoint system for the initial-boundary value problem for a general hyperbolic system using the standard variational approach requires that the inflow adjoint variables at an open boundary be zero. However, in general, these “natural” boundary conditions will lead to a different solution than that provided by the global assimilation problem. The impact of using natural boundary conditions when there are errors (on the boundary) in the initial guess on the assimilated initial conditions is discussed.

A proof of the uniqueness of the solution for both forward and adjoint equations in the presence of open boundaries at each iteration of the minimization procedure is provided, along with an assessment of the convergence of numerical solutions.

Numerical experiments with a simple advection equation support the theoretical analyses. Numerical results show that if observational data are perfect, 4DVAR data assimilation using a limited-area model can produce a reasonable initial condition. However, if there are errors in the observational data at the open boundaries and if natural boundary conditions are assumed, boundary errors can contaminate the assimilated solutions.

Corresponding author address: Dr. Chungu Lu, NOAA/OAR Forecast Systems Laboratory, Boulder, CO 80303.

Email: lu@fsl.noaa.gov

Abstract

Mathematical issues arising when applying four-dimensional variational (4DVAR) data assimilation to limited-area problems are studied. The derivation of the adjoint system for the initial-boundary value problem for a general hyperbolic system using the standard variational approach requires that the inflow adjoint variables at an open boundary be zero. However, in general, these “natural” boundary conditions will lead to a different solution than that provided by the global assimilation problem. The impact of using natural boundary conditions when there are errors (on the boundary) in the initial guess on the assimilated initial conditions is discussed.

A proof of the uniqueness of the solution for both forward and adjoint equations in the presence of open boundaries at each iteration of the minimization procedure is provided, along with an assessment of the convergence of numerical solutions.

Numerical experiments with a simple advection equation support the theoretical analyses. Numerical results show that if observational data are perfect, 4DVAR data assimilation using a limited-area model can produce a reasonable initial condition. However, if there are errors in the observational data at the open boundaries and if natural boundary conditions are assumed, boundary errors can contaminate the assimilated solutions.

Corresponding author address: Dr. Chungu Lu, NOAA/OAR Forecast Systems Laboratory, Boulder, CO 80303.

Email: lu@fsl.noaa.gov

Save
Cover Journal of the Atmospheric Sciences
Journal of the Atmospheric Sciences

  • Browning, G. L., and H.-O. Kreiss, 1986: Scaling and computation of smooth atmospheric motions. Tellus,38A, 295–313.

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

  • Errico, R. M., and T. Vukicevic, 1992: Sensitivity analysis using an adjoint of the PSU–NCAR Mesoscale Model. Mon. Wea. Rev.,120, 1644–1660.

  • Gustafsson, B., H.-O. Kreiss, and A. Sundström, 1972: Stability theory of difference approximations for mixed initial boundary value problems. II. Math. Comp.,26, 649–686.

  • Kreiss, H.-O., 1968: Stability theory of difference approximations for mixed initial boundary value problems. I. Math. Comp.,22, 703–714.

  • ——, and E. Lundqvist, 1968: On difference approximations with wrong boundary values. Math. Comp.,22, 1–12.

  • LeDimet, F. X., and O. Talagrand, 1986: Variational algorithms for analysis and assimilation of meteorological observations: Theory aspects. Tellus,38A, 97–110.

  • Lewis, J. M., and J. C. Derber, 1985: The use of adjoint equations to solve a variational adjustment problem with advective constraints. Tellus,37A, 309–322.

  • Lu, C., and G. L. Browning, 1998: The impact of observational and model errors on four-dimensional variational data assimilation. J. Atmos. Sci.,55, 995–1011.

  • ——, and ——, 2000: Discontinuous forcing generating rough initial conditions in 4DVAR data assimilation. J. Atmos. Sci.,57, 1646–1656.

  • Oliger, J., and A. Sundström, 1978: Theoretical and practical aspects of some initial boundary value problems in fluid dynamics. SIAM J. Appl. Math.,35, 419–446.

  • Sirkes, Z., and E. Tziperman, 1997: Finite difference of adjoint or adjoint of finite difference? Mon. Wea. Rev.,125, 3373–3378.

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

  • Xu, Q., 1996: Generalized adjoint for physical processes with parameterized discontinuities. Part I: Basic issues and heuristic examples. J. Atmos. Sci.,53, 1123–1142.

  • Zou, X., and Y.-H. Kuo, 1996: Rainfall assimilation through an optimal control of initial and boundary conditions in a limited-area mesoscale model. Mon. Wea. Rev.,124, 2859–2882.

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

  • Zupanski, M., 1993: Regional four-dimensional variational data assimilation in a quasi-operational forecasting environment. Mon. Wea. Rev.,121, 2396–2408.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 162 27 2
PDF Downloads 39 13 1