On the convergence of data assimilation for the one-dimensional shallow water equations with sparse observations

被引:0
作者
N. K.-R. Kevlahan
R. Khan
B. Protas
机构
[1] McMaster University,Department of Mathematics and Statistics
来源
Advances in Computational Mathematics | 2019年 / 45卷
关键词
Shallow water equations; Data assimilation; Tsunami; 35L05; 35Q35; 35Q93; 65M06;
D O I
暂无
中图分类号
学科分类号
摘要
The shallow water equations (SWE) are a widely used model for the propagation of surface waves on the oceans. We consider the problem of optimally determining the initial conditions for the one-dimensional SWE in an unbounded domain from a small set of observations of the sea surface height. In the linear case, we prove a theorem that gives sufficient conditions for convergence to the true initial conditions. At least two observation points must be used and at least one pair of observation points must be spaced more closely than half the effective minimum wavelength of the energy spectrum of the initial conditions. This result also applies to the linear wave equation. Our analysis is confirmed by numerical experiments for both the linear and nonlinear SWE data assimilation problems. These results show that convergence rates improve with increasing numbers of observation points and that at least three observation points are required for practically useful results. Better results are obtained for the nonlinear equations provided more than two observation points are used. This paper is a first step towards understanding the conditions for observability of the SWE for small numbers of observation points in more physically realistic settings.
引用
收藏
页码:3195 / 3216
页数:21
相关论文
共 63 条
[1]  
Lacasta A(2003)Calibration of the 1d shallow water equations: a comparison of Monte Carlo and gradient-based optimization methods J Hydroinformatics 19 282-298
[2]  
Morales-Hernández M(2003)Data assimilation (4d-VAR) for shallow-water flow: the case of the Chicoutimi river Vis. Geosci. 8 1-17
[3]  
Burguete J(2017)Determination of the bottom deformation from space- and time-resolved water wave measurements J. Fluid Mech. 835 301-326
[4]  
Brufau P(2010)Exact boundary observability for nonautonomous quasilinear wave equations J. Math. Anal. Appl. 364 41-50
[5]  
García-Navarro P(1997)Unified notation for data assimilation: operational, sequential and variational J. Meteorol. Soc. Jpn. 75 181-189
[6]  
Bélanger E(2015)Adaptive wavelet simulation of global ocean dynamics using a new Brinkman volume penalization Geosci. Model Dev. 8 3891-3909
[7]  
Vincent A(2006)Exact boundary observability for 1-D quasilinear wave equations Math. Method Appl. Sci. 29 1543-1553
[8]  
Fortin A(2015)Successive estimation of a tsunami wavefield without earthquake source data: a data assimilation approach toward real-time tsunami forecasting Geophys. Res. Lett. 42 7923-7932
[9]  
Cobelli PJ(1985)Surface deformation due to shear and tensile faults in a half-space Bull. Seism. Soc. Amer. 75 1135-1154
[10]  
Petitjeans P(2001)Tsunami waveform inversion by adjoint methods J. Geophys. Res.: Oceans 106 19773-19796