This article presents a new variational data assimilation (VDA) approach for the formal treatment of bias in both model outputs and observations. This approach relies on the Wasserstein metric, stemming from the theory of optimal mass transport, to penalize the distance between the probability histograms of the analysis state and an a priori reference dataset, which is likely to be more uncertain but less biased than both model and observations. Unlike previous bias-aware VDA approaches, the new Wasserstein metric VDA (WM-VDA) treats systematic biases of unknown magnitude and sign dynamically in both model and observations, through assimilation of the reference data in the probability domain, and can recover the probability histogram of the analysis state fully. The performance of WM-VDA is compared with the classic three-dimensional VDA (3D-Var) scheme for first-order linear dynamics and the chaotic Lorenz attractor. Under positive systematic biases in both model and observations, we consistently demonstrate a significant reduction in the forecast bias and unbiased root-mean-squared error.
|Original language||English (US)|
|Number of pages||15|
|Journal||Quarterly Journal of the Royal Meteorological Society|
|State||Published - Jul 1 2020|
Bibliographical noteFunding Information:
The first, second, and fourth authors acknowledge a grant from the National Aeronautics and Space Administration (NASA) Terrestrial Hydrology Program (THP, 80NSSC18K1528) through Dr. J. Entin and the New (Early Career) Investigator Program (NIP, 80NSSC18K0742) through Dr A. Leidner and Dr T. Lee. The third and fourth authors also acknowledge support from the National Science Foundation (NSF, 1830418).
National Aeronautics and Space Administration (NASA) Terrestrial Hydrology Program (THP, 80NSSC18K1528), the New (Early Career) Investigator Program (NIP, 80NSSC18K0742) and National Science Foundation (NSF,1830418). Funding information
© 2020 Royal Meteorological Society
- Wasserstein distance
- bias treatment
- chaotic systems
- optimal mass transport
- variational data assimilation