On the ensemble‐based linearization of numerical models
Abstract
Current parallelization trends in computer technology facilitates development of the algorithms that retrieve linear approximations of the model operators and their adjoints from ensembles of model simulations. In this study we address the problem of obtaining exact linearizations in the presence of semi‐implicit numerics of the parent model under realistic constraints on the ensemble size. The method is based on factorization of the model into a sequence of local and non‐local linear operators and employs prior information on the structure of the respective sparse matrices. The performance of the method is tested using 28 perturbed solutions of the shallow‐water equations with a moderate size () state vector. Numerical experiments have shown feasibility of the approach under relatively general constraints on the structure of the parent model. Because of the substantial expense of the ensemble‐based linearization, special focus is made on the assessment of the optimal frequency of such computations within the time intervals between data injections in typical operational systems.
Document Details
- Document Type
- Pub Defense Publication
- Publication Date
- Jan 01, 2020
- Source ID
- 10.1002/qj.3723
Entities
People
- Dmitri Nechaev
- Max Yaremchuk
- Sergey Frolov
Organizations
- John C. Stennis Space Center
- Office of Naval Research
- University of Colorado
- University of Southern Mississippi