On Sensitivity Analysis within the 4DVAR Framework

Max Yaremchuk Naval Research Laboratory, Stennis Space Center, Mississippi

Search for other papers by Max Yaremchuk in
Current site
Google Scholar
PubMed
Close
and
Paul Martin Naval Research Laboratory, Stennis Space Center, Mississippi

Search for other papers by Paul Martin in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

The sensitivity of model forecasts to uncertainties in control variables is evaluated using the adjoint technique and the ensemble generated by the reduced-order four-dimensional variational data assimilation (R4DVAR) algorithm within the framework of twin-data experiments with a quasigeostrophic model. To simulate real applications where the true state is unknown, the sensitivities were estimated using model solutions that were obtained after assimilating sparse observations extracted from the true solutions. The numerical experiments were conducted in the linear, weakly nonlinear, and strongly nonlinear (NL) regimes with special emphasis on the NL case characterized by the instability of the tangent linear model. It is shown that the ensemble-based R4DVAR method provides better sensitivity estimates in the NL case, primarily due to the better accuracy of the optimized solutions. The concept of sensitivity in the NL case is also considered within the statistical framework. Using analytical arguments and numerical experimentation, averaging the adjoint sensitivity estimates over an ensemble of model trajectories generated by finite perturbations of the optimal control is shown to provide an estimate similar to that obtained with the adjoint model stabilized by enhanced dissipation. This observation allows for evaluation of the sensitivities of strongly nonlinear optimal solutions by using both the adjoint (4DVAR) and ensemble (R4DVAR) optimization algorithms.

Naval Research Laboratory Contribution Number 7320-13-1728.

Corresponding author address: Max Yaremchuk, Naval Research Laboratory, Bldg. 1009, Stennis Space Center, MS 39529. E-mail: max.yaremchuk@nrlssc.navy.mil

Abstract

The sensitivity of model forecasts to uncertainties in control variables is evaluated using the adjoint technique and the ensemble generated by the reduced-order four-dimensional variational data assimilation (R4DVAR) algorithm within the framework of twin-data experiments with a quasigeostrophic model. To simulate real applications where the true state is unknown, the sensitivities were estimated using model solutions that were obtained after assimilating sparse observations extracted from the true solutions. The numerical experiments were conducted in the linear, weakly nonlinear, and strongly nonlinear (NL) regimes with special emphasis on the NL case characterized by the instability of the tangent linear model. It is shown that the ensemble-based R4DVAR method provides better sensitivity estimates in the NL case, primarily due to the better accuracy of the optimized solutions. The concept of sensitivity in the NL case is also considered within the statistical framework. Using analytical arguments and numerical experimentation, averaging the adjoint sensitivity estimates over an ensemble of model trajectories generated by finite perturbations of the optimal control is shown to provide an estimate similar to that obtained with the adjoint model stabilized by enhanced dissipation. This observation allows for evaluation of the sensitivities of strongly nonlinear optimal solutions by using both the adjoint (4DVAR) and ensemble (R4DVAR) optimization algorithms.

Naval Research Laboratory Contribution Number 7320-13-1728.

Corresponding author address: Max Yaremchuk, Naval Research Laboratory, Bldg. 1009, Stennis Space Center, MS 39529. E-mail: max.yaremchuk@nrlssc.navy.mil
Save
  • Ancell, B., and G. J. Hakim, 2007: Comparing adjoint- and ensemble-sensitivity analysis with applications to observation targeting. Mon. Wea. Rev., 135, 41174134.

    • Search Google Scholar
    • Export Citation
  • Baker, N., and R. Daley, 2000: Observation and background adjoint sensitivity in the adaptive observation-targeting problem. Quart. J. Roy. Meteor. Soc., 126, 14311454.

    • Search Google Scholar
    • Export Citation
  • Cao, Y., J. Zhu, I. M. Navon, and Z. Luo, 2007: A reduced-order approach to four-dimensional variational data assimilation using proper orthogonal decomposition. Int. J. Numer. Methods Fluids, 53, 15711583.

    • Search Google Scholar
    • Export Citation
  • Clayton, A. M., A. C. Lorenc, and D. M. Barker, 2013: Operational implementation of a hybrid ensemble/4D-Var global data assimilation system at the Met Office. Quart. J. Roy. Meteor. Soc., 139, 1445–1461, doi:10.1002/qj.2054.

    • Search Google Scholar
    • Export Citation
  • Farrell, B. F., and P. J. Ioannou, 2005: Approximating optimal state estimation. Predictability of Weather and Climate, T. N. Palmer and R. Hagedorn, Eds., Cambridge University Press, 181–216.

  • Gilbert, J., and C. Lemarechal, 1989: Some numerical experiments with variable-storage quasi-Newton algorithms. Math. Programming, 45, 407435.

    • Search Google Scholar
    • Export Citation
  • Hoteit, I., B. Cornuelle, A. Köhl, and D. Stammer, 2005: Treating strong adjoint sensitivities in tropical eddy-permitting variational data assimilation. Quart. J. Roy. Meteor. Soc., 131, 36593682.

    • Search Google Scholar
    • Export Citation
  • Hoteit, I., B. Cornuelle, and P. Heimbach, 2010: An eddy-permitting, dynamically consistent adjoint-based assimilation system for the tropical Pacific: Hindcast experiments in 2000. J. Geophys. Res., 115, C03001, doi:10.1029/2009JC005437.

    • Search Google Scholar
    • Export Citation
  • Lee, T., B. Cheng, and R. Giering, 2001: Adjoint sensitivity of Indonesian Throughflow transport to wind stress: Application to interannual variability. Jet Propulsion Laboratory Publ. 01-11, 34 pp.

  • Liu, C., Q. Xiao, and B. Wang, 2008: An ensemble-based four-dimensional variational data assimilation scheme. Part I: Technical formulation and preliminary test. Mon. Wea. Rev., 136, 33633373.

    • Search Google Scholar
    • Export Citation
  • Losch, M., and P. Heimbach, 2007: Adjoint sensitivity of an ocean general circulation model to bottom topography. J. Phys. Oceanogr., 37, 377393.

    • Search Google Scholar
    • Export Citation
  • Mirouse, I., and A. Weaver, 2010: Representation of correlation functions in variational data assimilation using an implicit diffusion operator. Quart. J. Roy. Meteor. Soc., 136, 14211443.

    • Search Google Scholar
    • Export Citation
  • Moore, A. M., H. G. Arango, E. D. Lorenzo, A. J. Miller, and B. D. Cornuelle, 2009: An adjoint sensitivity analysis of the Southern California Current circulation and ecosystem. J. Phys. Oceanogr., 39, 702720.

    • Search Google Scholar
    • Export Citation
  • Ngodock, H., S. Smith, and G. Jacobs, 2009: Cycling the representer method with non-linear models. Data Assimilation for Atmospheric, Oceanic and Hydrological Applications, S. K. Park and L. Xu, Eds., Vol. 1, Springer, 321–340.

  • Purser, R. J., W. Wu, D. F. Parrish, and N. M. Roberts, 2003: Numerical aspects of the application of recursive filters to variational statistical analysis. Part II: Spatially inhomogeneous and anisotropic covariances. Mon. Wea. Rev., 131, 15361548.

    • Search Google Scholar
    • Export Citation
  • Torn, R. D., and G. J. Hakim, 2008: Ensemble-based sensitivity analysis. Mon. Wea. Rev., 136, 663677.

  • Torn, R. D., and G. J. Hakim, 2009: Initial condition sensitivity of western Pacific extratropical transitions determined using ensemble-based sensitivity analysis. Mon. Wea. Rev., 137, 33883406.

    • Search Google Scholar
    • Export Citation
  • Veneziani, V., C. A. Edwards, and A. M. Moore, 2009: A central California coastal ocean modeling study: Adjoint sensitivity to local and remote forcing mechanisms. J. Geophys. Res., 114, C04020, doi:10.1029/2008JC004775.

    • Search Google Scholar
    • Export Citation
  • Xu, L., and R. Daley, 2000: Towards a true 4-dimensional data assimilation algorithm: Application of a cycling representer algorithm to a simple transport problem. Tellus, 52A, 109–128.

    • Search Google Scholar
    • Export Citation
  • Xu, Q., 2005: Representation of the inverse covariances by differential operators. Adv. Atmos. Sci., 22, 181194.

  • Yaremchuk, M., D. Nechaev, and G. G. Panteleev, 2009: A method of successive corrections of the control subspace in the reduced-order variational data assimilation. Mon. Wea. Rev., 137, 29662978.

    • Search Google Scholar
    • Export Citation
  • Yaremchuk, M., M. Carrier, S. Smith, and G. Jacobs, 2013: Background error correlation modeling with diffusion operators. Data Assimilation for Atmospheric, Oceanic and Hydrological Applications, S. K. Park and L. Xu, Eds., Vol. 2, Springer, 177–203.

  • Zack, J., E. Natenberg, S. Young, G. V. Knowe, K. Waight, J. Manobianco, and C. Kamath, 2010: Application of ensemble sensitivity analysis to observation targeting for short-term wind speed forecasting in the Washington–Oregon region. LLNL Tech. Rep. 458086, Lawrence Livermore National Laboratory, Livermore, CA, 65 pp.

  • Zhang, X., B. Cornuelle, and D. Roemmich, 2011: Adjoint sensitivity of the Niño-3 surface temperature to wind forcing. J. Climate, 24, 44804493.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 99 35 4
PDF Downloads 41 12 1