{"title":"Local Sequential MCMC for Data Assimilation with Applications in Geoscience","authors":"Hamza Ruzayqat, Omar Knio","doi":"arxiv-2409.07111","DOIUrl":null,"url":null,"abstract":"This paper presents a new data assimilation (DA) scheme based on a sequential\nMarkov Chain Monte Carlo (SMCMC) DA technique [Ruzayqat et al. 2024] which is\nprovably convergent and has been recently used for filtering, particularly for\nhigh-dimensional non-linear, and potentially, non-Gaussian state-space models.\nUnlike particle filters, which can be considered exact methods and can be used\nfor filtering non-linear, non-Gaussian models, SMCMC does not assign weights to\nthe samples/particles, and therefore, the method does not suffer from the issue\nof weight-degeneracy when a relatively small number of samples is used. We\ndesign a localization approach within the SMCMC framework that focuses on\nregions where observations are located and restricts the transition densities\nincluded in the filtering distribution of the state to these regions. This\nresults in immensely reducing the effective degrees of freedom and thus\nimproving the efficiency. We test the new technique on high-dimensional ($d\n\\sim 10^4 - 10^5$) linear Gaussian model and non-linear shallow water models\nwith Gaussian noise with real and synthetic observations. For two of the\nnumerical examples, the observations mimic the data generated by the Surface\nWater and Ocean Topography (SWOT) mission led by NASA, which is a swath of\nocean height observations that changes location at every assimilation time\nstep. We also use a set of ocean drifters' real observations in which the\ndrifters are moving according the ocean kinematics and assumed to have\nuncertain locations at the time of assimilation. We show that when higher\naccuracy is required, the proposed algorithm is superior in terms of efficiency\nand accuracy over competing ensemble methods and the original SMCMC filter.","PeriodicalId":501425,"journal":{"name":"arXiv - STAT - Methodology","volume":"78 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-09-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - STAT - Methodology","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.07111","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
引用次数: 0
Abstract
This paper presents a new data assimilation (DA) scheme based on a sequential
Markov Chain Monte Carlo (SMCMC) DA technique [Ruzayqat et al. 2024] which is
provably convergent and has been recently used for filtering, particularly for
high-dimensional non-linear, and potentially, non-Gaussian state-space models.
Unlike particle filters, which can be considered exact methods and can be used
for filtering non-linear, non-Gaussian models, SMCMC does not assign weights to
the samples/particles, and therefore, the method does not suffer from the issue
of weight-degeneracy when a relatively small number of samples is used. We
design a localization approach within the SMCMC framework that focuses on
regions where observations are located and restricts the transition densities
included in the filtering distribution of the state to these regions. This
results in immensely reducing the effective degrees of freedom and thus
improving the efficiency. We test the new technique on high-dimensional ($d
\sim 10^4 - 10^5$) linear Gaussian model and non-linear shallow water models
with Gaussian noise with real and synthetic observations. For two of the
numerical examples, the observations mimic the data generated by the Surface
Water and Ocean Topography (SWOT) mission led by NASA, which is a swath of
ocean height observations that changes location at every assimilation time
step. We also use a set of ocean drifters' real observations in which the
drifters are moving according the ocean kinematics and assumed to have
uncertain locations at the time of assimilation. We show that when higher
accuracy is required, the proposed algorithm is superior in terms of efficiency
and accuracy over competing ensemble methods and the original SMCMC filter.