基于卫星高度计资料同化的渤海和黄海M 2 潮流数值模拟*

南京大学学报(自然科学版) ›› 2010, Vol. 46 ›› Issue (6): 654–663.

基于卫星高度计资料同化的渤海和黄海M 2 潮流数值模拟*

张继才 1 ** , 吕咸青 2 , 孙丽艳 3
(

出版日期:2015-04-03 发布日期:2015-04-03
作者简介: 1. 南京大学海岸与海岛开发教育部重点实验室, 南京, 210093; 2. 中国海洋大学物理海洋教育部
重点实验室, 青岛, 266100; 3. 国家海洋局环境监测中心, 大连, 116023)
基金资助:
国家高技术研究发展计划( 2007AA09Z118) , 国家重点基础研究发展计划( 2005CB422308) , 南京大学人才引进启动经费, 国家自然科学基金( 41076006)

Three-dimensional simulation of M 2 tidal current in Bohai and yellow Seas by assimilating satellite altimetry

Zhang Ji Cai 1 , Lu Xian Qing 2 , Sun Li Yan 3

Online:2015-04-03 Published:2015-04-03
About author: ( 1. Ministry of Education Key Laboratory of Coast and Island Development, Nanjing University, Nanjing,
210093, China; 2. Ministry of Education Key Laboratory of Physical Oceanography, Ocean U niversity of
China, Qingdao, 266100, China; 3. National Marine Environmental Monitoring Center, Dalian, 116023, China)

摘要/Abstract

摘要： 将基于数学物理反问题的四维变分同化技术应用到陆架浅海区正压潮流的三维数值模拟中. 孪生实验结果表明, 借助伴随同化技术, 三维正压潮流模型的水位开边界条件可成功反演. 实际实验同
化从 TOPEX/ Poseidon 高度计资料提取的潮汐调和常数, 通过反演水位的第一边值开边界条件, 数值模拟了渤海和黄海的 M 2 潮汐潮流. 同化后代价函数以及代价函数关于开边界 Fourier 系数的梯度的值均
大幅度下降, 证明了伴随同化方法的有效性. 实验过程中潮流资料未参与同化, 但由其表达的代价函数值仍然由1 下降至 0 09, 表明从高度计数据反演陆架浅海区的潮流场具有高度的可行性. 与其它工作相
比, 潮流模拟值与观测值之间的误差较小, 所得同潮图及潮流椭圆分布图准确体现了渤海和黄海 M 2 潮汐潮流的基本特征. 本文工作表明伴随同化方法可有效提高潮流数值模拟的准确度.

Abstract: The four-dimensional variational data assimilation technology based on inverse problem is applied to simulate the tidal currents in marginal shelf seas by assimilating satellite altimetry. T he prescribed open boundary
conditions of a three dimensional barotropic tidal model are successfully inverted in the twin experiments. In the practical experiments, the open boundary conditions are optimized and the M 2 tides and tidal currents in Bohai and
Yellow Seas are simulated by assimilating tidal harmonic constants deduced from the T OPEX/Poseidon altimeter data. During the assimilation, the cost function and the gradients of cost function with respect to the open boundary
conditions are decreased significantly. Although the current observations are not assimilated, the cost function in terms of current data is still decreased from 1 to 0 ? 09, which demonstrates the feasibility of inverting tidal currents
from satellite altimetry or other elevation measurements. T he co?tidal charts and the near surface M 2 tidal current ellipses obtained in the practical experiment are in good agreement with the observed tides and tidal currents in the
Bohai and Yellow Seas. T his work indicates that the adjoint method can significantly improve the accuracy of three dimensional simulation of tidal currents in the marginal shelf seas.

张继才 1 ** , 吕咸青 2 , 孙丽艳 3
(. 基于卫星高度计资料同化的渤海和黄海M 2 潮流数值模拟*
[J]. 南京大学学报(自然科学版), 2010, 46(6): 654–663.

Zhang Ji Cai 1 , Lu Xian Qing 2 , Sun Li Yan 3
. Three-dimensional simulation of M 2 tidal current in Bohai and yellow Seas by assimilating satellite altimetry
[J]. Journal of Nanjing University(Natural Sciences), 2010, 46(6): 654–663.

参考文献

[ 1 ] Teague W J, Perkins H T, Hallock Z R, etal. Current and tide observations in the southern Yellow Sea. Journal of Geophysical Research, 1998, 103: 27783~ 27793.
[ 2 ] Munk W H. Once again: Once again-tidal friction. Progress in Oceanography, 1997, 40: 7~ 35.
[ 3 ] Egbert G D, Ray R D. Estimates of M 2 tidal energy dissipation from TOPEX/Poseidon altimeter data. Journal of Geophysical Research, 2001, 106(C10): 22475~ 22502.
[ 4 ] Fang G H. Numerical simulation on the M 2 tidal current in the eastern Yellow Sea and the comparison with observations. The 1 st Workshop on T ide and Sea Level. Tianjin: National
Marine Data and Information Service, 1986, 45~ 56. ( 方国洪. 黄海东部 M 2 潮流的数值计算及其与观测结果的比较. 第一届潮汐与海平面学术讨论会. 天津: 国家海洋局海洋科技情报研究所, 1986, 45~ 56).
[ 5 ] Guo X Y, Yanagi T. T hree dimensional structure of tidal current in the East China Sea and the Yellow Sea. Journal of Oceanography, 1998, 54: 651~ 668.
[ 6 ] Lee J C, Jung K T. Application of eddy viscosity closure models for the M 2 tide and tidal currents in the Yellow Sea and the East China Sea. Continental Shelf Research, 1999, 19: 445~ 475.
[ 7 ] Lefevre F, Le Provost C, Lyard F H. How can we improve a global ocean tide model at a regional scale A test on the Yellow Sea and the East China Sea. Journal of Geophysical Research, 2000, 105: 8707~ 8725.
[ 8 ] Bao X W, Gao G P, Yan J. Three dimensional simulation of tide current characteristics in the East China Sea. Oceanologica Acta, 2001, 24 (2) : 1~ 15.
[ 9 ] Lee H J, Jung K T, So J K, et al. A three dimensional mixed finite difference Galerkin function model for the oceanic circulation in the Yellow Sea and the East China Sea in the presence of M 2 tide. Continental Shelf Research, 2002, 22: 67~ 91.
[ 10] Li P L. The study of the tidal numerical simulation with assimilation and the tidal dissipation in the Bohai Sea, Yellow Sea and East China Sea. Thesis for Doctor% s Degree, Ocean University
of China, Qingdao, 2002. ( 李培良. 渤黄东海潮波同化数值模拟和潮能耗散的研究. 博士学位论文. 青岛: 中国海洋大学海洋环境学院, 2002) .
[ 11] He Y J, Lu X Q, Qiu Z F, etal. Shallow water tidal constituents in the Bohai Sea and the Yellow Sea from a numerical adjoint model with T/ P altimeter data. Continental Shelf Research, 2004, 24: 1521~ 1529.
[ 12] Li L, Du L, Zuo J C, et al. Tidal current simulation for M 2 and K 1 with the finite element method in the Bohai Sea, the Yellow Sea and the East China Sea. Periodical of Ocean Univer
sity of China, 2006, 36(6) : 851~ 858. ( 李磊, 杜凌, 左军成等. 渤、黄、东海 M 2 和 K 1 分潮潮流场的有限元模拟. 中国海洋大学学报, 2006, 36(6): 851~ 858) .
[ 13] Zhang J C, Lu X Q. On the simulation of M 2 tide in the Bohai, Yellow and East China Seas with TOPEX/ Poseidon altimeter data. Terrestrial, Atmospheric and Oceanic Sciences, 2008, 19( 1~ 2): 173~ 181.
[ 14] Ray R D. Inversion of oceanic tidal currents from measured elevations. Journal of Marine System, 2001, 28: 1~ 18.
[ 15] Lardner R W, Das S K. Optimal estimation of eddy viscosity for a quasi three dimensional numerical tidal and storm surge model. Interna tional Journal for Numerical Methods in Fluids, 1994, 18: 295~ 312.
[ 16] Zhang A J, Wei E, Parker B B. Optimal estimation of tidal open boundary conditions using predicted tides and adjoint data assimilation technique. Continental Shelf Research, 2003, 23: 1055~ 1070.
[ 17] Cartwright D E, Ray R D, Sanchez B V. A computer program for predicting oceanic tidal currents. NASA T ech Memo 104578, Goddard Space Flight Center, 1992, 21.
[ 18] Headrick R H, Spiesberger J L, Bushong P J. T idal signals in basin?scale acoustic transmissions. Journal of the Acoustical Society of America, 1993, 93: 790~ 802.
[ 19] Kantha L H. Barotropic tides in the global oceans from a nonlinear tidal model assimilating altimetric tides 1. Model description and re sults. Journal of Geophysical Research, 1995, 101: 25283~ 25308.
[ 20] M cIntosh P C, Bennett A F. Open ocean modeling as an inverse problem: M 2 tides in Bass Strait. Journal of Physical Oceanography, 1984, 14: 601~ 614.
[ 21] Egbert G D, Bennett A F, Foreman M G G. T opex/ Poseidon tides estimated using a global inverse method. Journal of Geophysical Research, 1994, 99: 24821~ 24852.
[ 22] Shum C K, Woodworth P L, Andersen O B, et al. Accuracy assessments of recent global ocean tide models. Journal of Geophysical Research, 1997, 102(C11): 25173~ 25194.
[ 23] Egbert G D, Ray R D. Significant dissipation of tidal energy in the deep ocean inferred from satellite altimeter data. Nature, 2000, 405: 775~ 778.
[ 24] Zhang J C, Lu X Q. Parameter estimation for a three dimensional numerical barotropic tidal model with adjoint method. International Journal for Numerical Methods in Fluids, 2008, 57: 47~ 92.
[ 25] Fang G H, Wang Y G, Wei Z X, etal. Empirical cotidal charts of the Bohai, Yellow, and East China Seas from 10 years of T/ P altimetry. Journal of Geophysical Research, 2004, 109: C11006.
[ 26] Tihonov A N. Regularization of incorrectly posed problems. Soviet mathematics Doklady, 1963, 4: 1624~ 1627.
[ 27] Yeh W W G. Review of parameter identification procedures in groundwater hydrology: The inverse problem. Water Resources Research, 1986, 22: 95~ 108.
[ 28] Smedstad O M, O% Brien J J. Variational data assimilation and parameter estimation in an equatorial Pacific Ocean model. Progress in Oceanography, 1991, 26: 179~ 241.
[ 29] Fang G H, Zheng W Z, Chen Z Y, et al. Analysis and prediction of tides and tidal currents. Beijing: Ocean Press, 1986, 63~ 67. ( 方国洪, 郑文振, 陈宗镛等. 潮汐和潮流的分析与预报. 北京: 海洋出版社, 1986, 63~ 67) .
[ 30] Yang Y, Wang Y P, Gao J H, et al. Hydroolynamics and fire grained sediment resuspension in Changjiang estuary dry seasons. Journal of Nanjing University ( Natural Sci?
ences), 2006, 42( 6) : 643~ 655. (杨旸, 汪亚平, 高建华等. 长江口枯季水动力悬沙特征与再悬浮研究. 南京大学学报(自然科学) , 2006, 42(6): 643~ 655).

Metrics

Viewed

Full text

Abstract

Cited

Shared

Discussed