A three-dimensional numerical internal tidal model involving adjoint method

被引:17
作者
Chen, Haibo [1 ]
Miao, Chunbao [1 ]
Lv, Xianqing [1 ]
机构
[1] Ocean Univ China, Lab Phys Oceanog, Qingdao 266100, Peoples R China
基金
中国国家自然科学基金;
关键词
Parameter optimization; adjoint method; 3D internal tidal model; optimization; marine hydrodynamics; ocean flows; fluid-structure interaction; shallow water; Lagrangian; free surface; OPEN BOUNDARY-CONDITIONS; DETERMINING INITIAL CONDITIONS; VARIATIONAL DATA ASSIMILATION; GENERAL-CIRCULATION MODEL; ATMOSPHERE-OCEAN MODELS; SIMPLE EQUATORIAL MODEL; SHALLOW-WATER FLOW; FREE-SURFACE FLOWS; PARAMETER-ESTIMATION; METEOROLOGICAL OBSERVATIONS;
D O I
10.1002/fld.2650
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
A three-dimensional internal tidal model involving the adjoint method is constructed based on the nonlinear, time-dependent, free-surface hydrodynamic equations in spherical coordinates horizontally, and isopycnic coordinates vertically, subject to the hydrostatic approximations. This model consists of two submodels: the forward model is used for the simulation of internal tides, while the adjoint model is used for optimization of modal parameters. Mode splitting technique is employed in both forward and adjoint models. In this model, the adjoint method is employed to estimate model parameters by assimilating the interior observations. As a preliminary feasibility study, a set of ideal experiments with the model-generated pseudo-observations of surface currents are performed to invert the open boundary conditions (OBCs). In the ideal experiments, 14 kinds of bottom topographies and six kinds of predetermined distributions of OBCs are considered to examine their influence on experiment results. The inversion obtained satisfying results and all the predetermined distributions were successfully inverted. Analysis of results suggests the following: in the case where the spatial variation of the OBC distribution is great or the open boundary is close to a rough topography, the results will be comparatively poor, but still satisfactory; both the tidal elevations and currents can be simulated very accurately with the surface currents at several observation points; the assimilation precision could be reliable and able to reflect both of the inversion and simulation results in the whole field. The performance and results of ideal experiments give a preliminary indication that the construction of this model is successful. Copyright (c) 2011 John Wiley & Sons, Ltd.
引用
收藏
页码:1584 / 1613
页数:30
相关论文
共 109 条
[1]   Comparison of advanced large-scale minimization algorithms for the solution of inverse ill-posed problems [J].
Alekseev, A. K. ;
Navon, I. M. ;
Steward, J. L. .
OPTIMIZATION METHODS & SOFTWARE, 2009, 24 (01) :63-87
[2]   Data assimilation in ocean models [J].
Anderson, DLT ;
Sheinbaum, J ;
Haines, K .
REPORTS ON PROGRESS IN PHYSICS, 1996, 59 (10) :1209-1266
[3]  
[Anonymous], MEM SOC R SCI LIEGE
[4]  
[Anonymous], 2006, Dynamic Data Assimilation: a Least Squares Approach
[5]  
Arakawa A., 1977, Methods in Computational Physics: Advances in Research and Applications, P173, DOI [10.1016/B978-0-12-460817-7.50009-4, DOI 10.1016/B978-0-12-460817-7.50009-4]
[6]   Estimation of boundary values in a North Atlantic circulation model using an adjoint method [J].
Ayoub, N .
OCEAN MODELLING, 2006, 12 (3-4) :319-347
[7]  
Backhaus J.O., 1985, Ocean Dynamics, V38, P165
[8]  
Bennett A.F., 1992, INVERSE METHODS PHYS
[9]  
Bennett A. F., 2002, Inverse modeling of the ocean and atmosphere
[10]  
BENNETT AF, 1982, J PHYS OCEANOGR, V12, P1004, DOI 10.1175/1520-0485(1982)012<1004:OOMAAI>2.0.CO