Local Sequential MCMC for Data Assimilation with Applications in Geoscience

Hamza Ruzayqat, Omar Knio
{"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.
查看原文
分享 分享
微信好友 朋友圈 QQ好友 复制链接
本刊更多论文
应用于地球科学数据同化的局部序列 MCMC
本文提出了一种基于序列马尔可夫链蒙特卡洛(SMCMC)数据同化(DA)技术[Ruzayqat et al. 2024]的新数据同化(DA)方案,该方案具有明显的收敛性,最近已被用于滤波,特别是高维非线性和潜在的非高斯状态空间模型的滤波。粒子滤波器被认为是精确的方法,可用于非线性、非高斯模型的滤波,与粒子滤波器不同的是,SMCMC 不给样本/粒子分配权重,因此,在使用相对较少的样本时,该方法不会出现权重退化的问题。我们在 SMCMC 框架内设计了一种本地化方法,该方法侧重于观测值所在的区域,并将状态滤波分布中包含的过渡密度限制在这些区域内。这大大减少了有效自由度,从而提高了效率。我们在高维($d\sim 10^4 - 10^5$)线性高斯模型和具有高斯噪声的非线性浅水模型上测试了新技术,并进行了真实和合成观测。对于其中的两个数值示例,观测数据模仿了美国国家航空航天局(NASA)领导的地表水和海洋地形学(SWOT)任务生成的数据,这是一个在每个同化时间步都会改变位置的海洋高度观测带。我们还使用了一组海洋漂流者的真实观测数据,其中漂流者是根据海洋运动学原理移动的,并假定其在同化时的位置是不确定的。我们的研究表明,当需要更高的精度时,所提出的算法在效率和精度方面都优于其他同类集合方法和原始的 SMCMC 滤波器。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 去求助
来源期刊
自引率
0.00%
发文量
0
期刊最新文献
Poisson approximate likelihood compared to the particle filter Optimising the Trade-Off Between Type I and Type II Errors: A Review and Extensions Bias Reduction in Matched Observational Studies with Continuous Treatments: Calipered Non-Bipartite Matching and Bias-Corrected Estimation and Inference Forecasting age distribution of life-table death counts via α-transformation Probability-scale residuals for event-time data
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
已复制链接
已复制链接
快去分享给好友吧!
我知道了
×
扫码分享
扫码分享
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1