西北太平洋的一种潮汐数值同化模型
黄海、渤海TOPEXPOSEDON高度计资料潮汐伴随同化
第27卷 第4期海 洋 学 报V ol.27,No.42005年7月ACTA OCEANOLOGICA SINICAJuly 2005黄海、渤海TOPEX/Poseidon 高度计资料潮汐伴随同化丘仲锋1,2,何宜军1,吕咸青3(11中国科学院海洋研究所,山东青岛266071;21中国科学院研究生院,北京100039;31中国海洋大学,山东青岛266003)收稿日期:2004-04-28;修订日期:2005-04-10.基金项目:物理海洋教育部重点实验室开放课题基金项目(200307);中国科学院海洋研究所领域前沿项目/卫星高度计在海洋动力过程中的应用研究0(KZCX2-202).作者简介:丘仲锋(1979)),男,福建省上杭县人,博士,从事物理海洋和遥感海洋学研究.E -mail:qiuz hongfeng@摘要:首先将大约10a 的TOPEX/Po seidon(T /P)高度计资料沿星下轨迹点做潮汐调和分析,提取得到各分潮的调和常数,利用伴随同化方法,同化到二维非线性潮汐数值模式中,模拟了黄海、渤海区域M 2,S 2,O 1,K 1等4个潮汐分潮,并根据计算结果给出了各分潮的同潮图.将计算值与观测值的进行偏差统计,结果表明计算值与验潮站资料符合良好.研究过程中做了两类试验:一类试验是针对不同的参数进行优化,一类试验是针对不同的资料进行同化.第一类试验表明:将开边界条件和底摩擦系数同时作为模型优化的控制参数,其结果明显优于单独优化开边界条件;第二类试验表明:同时同化高度计资料与验潮站资料,比单独同化其中任一种资料,对模式计算结果都有较好的改进.研究结果表明,采用伴随同化方法,利用T /P 高度计资料和验潮站资料作为同化数据能有效改进模拟结果,用来反演黄海、渤海的潮波系统是可行的.关键词:伴随同化;T OPEX/Poseidon 高度计;潮汐中图分类号:P731.23 文献标识码:A 文章编号:0253-4193(2005)04-0010-091 引言渤海为辽东半岛和山东半岛所环抱,是一个典型的半封闭海湾,平均水深20m.黄海西面与北面为中国大陆,在西北方向与渤海相通,东面为朝鲜半岛,东南与东海相通,是一个半封闭的陆架浅海,最大水深140m.在黄海、渤海区域,潮波运动是主要运动形式之一.渤海、黄海的潮差分布悬殊,潮汐类型复杂,是研究潮波运动的典型海区.早在1980年以前,安希沫[1]等曾对中国邻近海域潮汐进行过研究;1980年以来,许多学者,如沈育疆[2]、夏综万[3]、Fang [4]、丁文兰[5]等,对中国东部海域的潮汐做了大量研究工作.关于东中国东部海域潮汐、潮流的数值模拟,许多学者,如方国洪[6]、赵保仁[7]、叶安乐[8]和万振文等[9,10]等做过大量工作,所得结果与观测资料较为接近.然而,上述文献,观测资料更多用来检验计算结果而不是参与模式运算.传统的潮汐测量仅靠沿岸或沿岛屿的潮汐观测站和几百米深海底的压力仪,所得资料分布较为稀疏.20世纪70年代后期卫星高度计的发展使深海潮汐测量成为可能,特别是90年代升空的T OPEX/Poseido n(以下简称T /P)卫星,提供了大量高精度的测高数据,使潮汐模式获得长足发展,精度越来越高.目前高度计潮汐模式有十几种之多,这些模式用于大洋具有很高的精度,但在浅海和近岸区,不一定适用.Lefevre [11]等将全球潮汐数值模式应用到中国东部海域,其主要分潮均方根偏差M 2为15182cm,S 2为13115cm,K 1为5109cm,O 1为3176cm,结果不如人意.数值模拟是潮波研究的一种有力手段,但在研究中会面临各种具体问题,包括开边界条件的确定、底摩擦系数和耗散系数的选取等.数据同化是解决这些问题的一种途径,即利用有限数量的潮汐观测资料对潮波场进行最优估计.变分同化(伴随同化)是数据同化的一种,变分的目的是设法使观测值尽量接近动力约束的数值解.Lagrange乘子法是一种有效的变分方法.该方法建立在严格的数学基础上,将所要解决的实际问题作为条件最小问题来解,能有效地提高数值模拟的精度.关于变分同化,许多学者作了大量工作,用于解决各种实际问题.比如,Lar der[12]、Seiler[13]及朱江等[14]曾对开边界条件的反演进行过研究.对于中国近海潮汐研究的伴随同化研究,更多集中于方法的探讨,而针对实际海区设计的同化试验较少.吕咸青等[15,16]探讨了实际海区潮汐模型中开边界条件的处理,并利用高度计资料和验潮站调和常数反演了中国东部海域的开边界条件;韩桂军等[17]利用伴随同化,模拟了黄海、东海M2分潮;吴自库等[18]把T/ P资料同化到潮汐潮流数值模式中,模拟了南海的m1,M2分潮.他们虽然对中国近海潮汐的伴随同化作了研究,但都局限于m1,M2分潮.本研究将T/P 高度计资料和验潮站资料同时作为同化数据,模拟了渤海、黄海的M2,S2,O1,K1分潮.本研究采用伴随同化方法,下一节将进行详细描述.第三节介绍了对T/P高度计资料的误差校正,并沿星下轨迹点做潮汐调和分析,提取各分潮的调和常数.最后为计算结果分析和讨论,计算的M2,S2,O1,K1等4个分潮调和常数与实测结果比较,符合程度良好,说明能有效地利用T/P高度计资料和验潮站资料来反演黄海、渤海的潮波系统.2数值模式本研究计算海区(图1所示)为黄海、渤海区域(34b~41b N,117b~127b E)计算网格空间分辨率为(1/6)b,开边界取在34b N.本文采用如下潮汐模式: 9t F=-A-1[9K(H u)+9<(H v cos<),9t u=-A-1u9K u-R-1v9<u+R-1uv tan<+f v-C D H-1(u2+v2)"-u-g A-19K(F-F)+A$U,(1) 9t v=-A-1u9K v-R-1v9<v+R-1u2tan<+f u-C D H-1(u2+v2)"-v-gR-19<(F-F)+A$v,式中,t为时间;K,<为是地球东经和北纬;R为地球半径;g是重力加速度,A=R cos<,f=28sin<,8为地球自转角速度;H=h+F为总水深,h为未扰动水深,F为海面相对于未扰动水深的高度,F为考虑了地潮效应调整后的平衡潮;C D为底摩擦系数;A为水平涡动黏性系数;u和v分别为东向和北向上的流速分量;$为拉普拉斯算子.图1计算海区水深(m)示意图#表示高度计资料点,n表示验潮站资料点为了得到方程组的惟一解,需给定适定的边界条件.这里取闭边界(岸界、岛界)处法向流速为0,即在闭边界上没有入流和出流;与外海相通的开边界处水位由调和常数给定:F=H cos(X t-g),(2)其中,H和g分别为各分潮的振幅和迟角,X为各分潮角频率.在球坐标系下,以$K,$<和$t分别代表经向、纬向和时间步长;采用Ar akaw a C网格,即水位取在网络的中心,而通量则取在网络的边缘,F n i,j=F(K0+i$K,<0+j$<,n$t),u n i,j=u(K0+(i+12)$K,<0+j$<,n$t),v n i,j=v(K0+i$K,<0+(j+12)$<,n$t).(3)对方程组采用如下的差分离散方法:F n+12i,j-F n i,j$t+H i+12,j u n i,j-H i-12,j u n i-1,jA j$K+H i,j+12v n i,j-H i,j-12v n i,j-1R$<=0,(4)v n+12i,j-v n i,j$t+fën i,j+gRF n+12i,j+1-F n+12i,j-F n+12i,j+1+F n+12i,j$<+C Di,j+12H i,j+12u n2i,j+v n2i,j xv n i,j+(1+x)v n+12i,j-114期丘仲锋等:黄海、渤海T O PEX/Poseidon高度计资料潮汐伴随同化A$v n i,j+u n i+1,j(v n i+1,j-v n i-1,j)2A j$K+v n i,j(v n i,j+1-v n i,j-12R$<+u n2i,jRtan<=0,(5)u n+12i,j-u n i,j$t+f v n+12i,j+gA jF n+12i+1,j-F n+12i,j-F n+12i+1,j+F n+12i,j$K+C Di+12,jH i+12,ju n2i,j+v n2i,j x u n i,j+(1-x)u n+12i,j-A$u n i,j+u n i,j(u n i+1,j-u n i-1,j)2A j$K+v n+12i,j(u n i,j+1-u n i,j-1)2R$<-u n i,j v n+12i,jRtan<=0,(6)F n+1i,j-F n+12i,j$t+H i+12,j u n+12i,j-H i-12,j u n+12i-1,jA j$K+H i,j+12v n+12i,j-H i,j-12v n+12i,j-1R$<=0,(7)u n+1i,j-u n+12i,j$t+f v n+1 2i,j+gA jF n+1i+1,j-F n+1i,j-½F n+1i+1,j+F n+1i,j$K+C Di+12,jH i+12,ju(n+12)2i,j+v(n+12)2i,j x u n+12i,j+(1-x)u n+1i,j-A$u n+12i,j+u n+12i,j(u n+12i,j+1-u n+12i,j-1)2A j$K+vn+12i,j(u n+12i,j+1-u n+12i,j-1)2R$<+u n+12 i,j v n+1 2i,jRtan<=0,(8)v n+1i,j-v n+12i,j$t+f u n+1i,j+gRF n+1i,j+1-F n+1i,j-F n+1i,j+1+F n+1i,j$<+C Di,j+12 H i,j+12u(n+12)2i,j+v(n+12)2i,j x v n+12i,j+(1-x)v n+1i,j-A$v n+12i,j+u n+1i,j(v n+12i+1,j-v n+12i-1,j)2A j$K+v n+12i,j(v n+12i,j+1-v n+1 2i,j-1) 2R$<+u(n+1)2i,jRtan<=0,(9)式中,x表示权重系数,这里取x=015,u n i,j=(u n i,j+u n i-1,j+u n i,j+1+u n i-1,j+1)/4,(10) v n i,j=(v n i,j+v n i+1,j+v n i,j-1+v n i+1,j-1)/4.(11)格式计算稳定性要求时间步长满足$tìmin[$x/2gh i,j],(12) $x是最小的网格格距,在这里$x=$K=$<= (1/6)b.计算中初始条件取/冷条件0,即当t=0时,F= u=v=0.当计算稳定后(相邻两周期同一点F之差小于0101m;u,v之差小于0101m/s),再计算一个潮周期,把它作为同化窗口,且对此周期内的时间序列进行调和分析,得出各计算点上潮汐调和常数.利用从T/P资料中提取的沿星下轨迹点的各分潮调和常数,及部分验潮站的调和常数,插值到就近的网格点上,生成观测值:F^=EMi=1H i cos(X i t-g i).(13)式中,{X i B i=1,,,M0}是Fourier模态的角频率, M0是分潮的个数.变分伴随的目的就是通过调整模型控制参数,使模拟值与观测值在最小二乘意义下拟合得最好.因此,定义如下的代价函数:J=12K N E N n=I+1E(i,j)I DN(F n i,j-F^n i,j)2,(14)式中,K N是权重,D N是观测资料所在点的集合(包括高度计资料、验潮站资料).令F*,u*,v*分别为F,u,v的伴随变量,以式(14)为目标函数,以差分方程组为约束的Lag rang e函数为:L=E N n=I+1[E(i,j)I SF(F*n i,j(4)+F*n+12i,j(7))+E(i,j)I Su(u*n+12i,j(6)+u*n+1i,j(8))+E(i,j)I S v(v*n+12i,j(5)+v*n+1i,j(9))]+12K0E k I D{(a0,k-a c0,k)2+E M0l=1[(a l,k-a c l,k)2+(b l,k-b c l,k)2]}+J,(15)其中,(4)~(9)指方程(4)~(9)的左边项,D0为开边界点的集合,{a l,k B l=0,,,M0;b l,k B l=l,,,M0}是开边界点处的Fourier系数,{a c l,k B l=0,,,M0;b c l,k B l=1,,,M0}为相应Fourier系数的先验估计.这样,在差分方程的约束下,求代价函数(14)最小值问题就转化成无约束问题.为使代价函数达到最小值,要求9L9F n i,j=0,9L9F n+12i,j=0,9L9u n i,j=0,9L9u n+12i,j=0,9L9v n i,j=0,9L9v n+12i,j=0,(16)9L9F*n i,j=0,9L9F*n+12i,j=0,9L9u*n i,j=0,9L9u*n+12i,j=0,9L9v*n i,j=0,9L9v*n+12i,j=0,(17)9L9a l,k=0,9L9b l,k=0,(18) 9L9C Di,j=0.(19)式(17)得到原差分形式的模式方程,式(16)导12海洋学报27卷出其伴随模式,式(18)导出代价函数关于开边界的控制梯度表达式,式(19)导出代价函数关于底摩擦系数的控制梯度表达式.这样就构成一个伴随同化系统,主要求解过程如下:(1)正向积分模式差分方程;(2)稳定后再计算一个潮周期作为同化窗口;(3)反向积分伴随模式;(4)计算模型控制变量的梯度;(5)利用最速下降算法解式(14);判断是否最优:是,迭代终止;不是,返回(1).3T/P资料的处理这里分析使用的卫星高度计资料为NA SA提供的T/P资料.T/P卫星在计算区域的轨道分布如图1所示.采用了从1992年10月至2002年7月的高度计数据,对应T/P卫星周期数为2~363,总的时间跨度将近10a.数据处理分为以下两个步骤. 311T/P高度计数据的误差校正处理校正算式为:H=H raw-H i-H d t-H wt-H eb-H ib-E r-H g,(20)式中,H raw是原始海面高度数据;H i是电离层校正量;H dt,H wt分别是干、湿对流层校正量;H eb是电磁偏差校正量,H ib是大气逆压校正量;H g是大地水准面高度;E r是轨道误差.312T/P测高数据潮波分析针对中国近海实际潮汐情况,选择主要分潮进行分析.将沿轨迹任何一点的测高数据系列表示为:h=S0+E M m=1H m f m cos(R m t+u m-H m),(21)式中,H m=g m-v0m,H和g为各分潮调和常数,v0为零时格林威治天文角,R为分潮角频率,S0为多年平均海平面,f m为交点因子,u m为交点订正角.利用最小二乘法可以求出沿轨迹点的潮汐分潮调和常数,并插值到沿迹点附近计算网格点上.4计算结果分析讨论为了对数值结果进行检验与比较,本研究采用了两种偏差计算方法:一种是分别比较振幅和迟角的平均偏差,称为绝均差.即$H mean=1K E Kk=1|H cal,k-H obs,k|,g mean=1KE Kk=1|g cal,k-g obs,k|,(22)式中,下标mean,cal和obs分别代表平均值、计算值和观测值;K为总观测站数.另一种方法与统计学常用的偏差估计方法较一致,称为均方根偏差,R=1KE Kk=1[(a cal,k-a obs,k)2+(b cal,k-b obs,k)2]1/2,(23)式中,a=H cos g,b=H sin g.R也被称为计算值与观测值之间的/距离0,表征了计算值与观测值的偏离程度.411优化不同参数的实验在潮汐数值模拟中,开边界条件的确定和底摩擦系数的选取是一个难点.如果能给定准确的开边界条件和合适的底摩擦系数,数值模拟的精度将得到提高.本文在参数化实验中做了如下两组试验:试验1:只优化开边界条件只将开边界条件作为优化的控制变量,底摩擦系数根据经验在渤海取111@10-3,在黄海取115@10-3.同化资料为T/P高度计资料提取的沿轨迹分潮调和常数,共选取145个观测点.试验2:同时优化开边界条件与底摩擦系数将开边界条件与底摩擦系数同时作为优化的控制变量,初始底摩擦系数取为2@10-3.同化资料与试验1相同.表1和表2给出了两组试验的结果比较.表1给出计算结果与高度计观测值的偏差比较,表2给出计算值与验潮站资料的偏差比较.可以看出,两组试验,计算结果与观测值的偏差都能控制在一个较小的范围之内,试验1只有O1分潮的偏差与试验2较为接近;半日分潮计算结果与观测值的偏差,试验1都比试验2大得多.结果说明,同时优化开边界条件和底摩擦系数能比单独优化开边界条件取得更好的反演结果.因此表明,在潮汐数值伴随同化中,计算结果受模型参数的影响是多方面的,进行多参数的优化,能比单独将某个参数作为优化的控制变量取得更好的结果.表1M2,S2,O1,K1分潮计算值与高度计观测值的偏差比较分潮高度计资料个数振幅绝均差/cm试验1试验2迟角绝均差/(b)试验1试验2均方根偏差/cm试验1试验2M21457172196172151413516S21452161121118318613213O1145113114518419217215K1145116114615412410219134期丘仲锋等:黄海、渤海T O PEX/Poseidon高度计资料潮汐伴随同化表2M2,S2,O1,K1分潮计算值与验潮站观测值的偏差比较分潮高度计资料个数振幅绝均差/cm试验1试验2迟角绝均差/(b)试验1试验2均方根偏差/cm试验1试验2M232141771051051821181714 S23231831114165171112715O131410412515418511514K1292183103612512512418312同化不同资料的试验在试验2的基础上,针对不同资料的同化,做了如下3组试验.试验3:单独同化T/P高度计资料;试验4:单独同化验潮站资料;试验5:同时同化T/ P高度计资料和验潮站资料.试验3的工作与试验2相同,选取145个高度计观测点的资料作为同化数据,进行伴随同化研究;试验4挑选出21个质量好的验潮站资料作为同化数据,模式的其他处理与试验1相同;试验5将高度计资料和验潮站资料一起作为同化数据,共选取了166个观测点的资料.对计算结果的检验采用了与试验2相同的验潮站资料.表3M2,S2,O1,K1分潮计算值与观测值的偏差比较分潮资料个数同化检验振幅绝均差/cm同化检验迟角绝均差/(b)同化检验均方根偏差/cm同化检验M2试验3145322197102155185161714试验4213241041231741710171417试验5166323135133163186161214 S2试验314532112311318517213715试验42132213212316514519716试验516632115215413416310517 O1试验314531114412419418215514试验42131117118319415218310试验516631114315417414217415 K1试验3145291143103412512219418试验42129215219515417419415试验516629116218417418314414表4不同研究者计算和实测偏差值的比较分潮项目赵保仁等[7]叶安乐等[8]万振文等[9]万振文等[10]本研究M2振幅绝均差/cm712612417819513迟角绝均差/(b)614417511713318 S2振幅绝均差/cm318215迟角绝均差/(b)913416 O1振幅绝均差/cm311116213218218迟角绝均差/(b)714419414615418 K1振幅绝均差/cm419315迟角绝均差/(b)71041414海洋学报27卷表5不同研究者M2调和常数计算与实测差值的比较站位北纬度东经度赵保仁等[7]$H/cm$g/(b)万振文等[9]$H/cm$g/(b)本研究$H/cm$g/(b)34b34c126b01c9-81319-1016-1018316 34b45c119b27c2-4316-719017211 35b03c126b05c192711-0181217211 35b23c119b33c2-21711-4111011518 35b21c126b01c0-5-014-3120013 36b05c120b19c-15-6-012-315-913-218 36b23c126b26c-4-14-418-5120-814 36b53c122b25c0-5218-2124144 37b30c122b10c-519316-218112-517 37b44c125b33c1318-2196-912512 38b10c120b45c-59-113-314015017 39b04c123b09c-160-517-013-414-215 39b29c123b05c-418-1194-013 39b42c124b24c-80-512-112-516312 40b43c121b00c-1720-113-315-915-416最大偏差1920171110161018814表3给出了3组试验计算结果与观测值的偏差比较,这里的偏差有两种:一是计算值与用来同化的资料的偏差比较;一是计算值与用来检验的验潮站资料的偏差比较.从表3可以看出:计算结果与同化资料的偏差,试验3的值最小,试验5其次,试验4最大;计算结果与用来检验的验潮站资料的偏差比较,试验5能取得最好的效果,试验4其次,试验3偏差最大.其中O1分潮例外,试验4取得的效果最好,试验3偏差值最大.试验结果表明:以同样的验潮站资料作为计算结果的检验,试验4对改进沿岸潮汐模拟的结果要稍好于试验3,而试验5能取得最好的效果.这是由于试验3仅利用了T/P高度计资料,单独的卫星资料同化对模型在浅水沿岸处的结果改进有限;试验4仅利用了验潮站资料,使得同化的数据只能局限于近岸近岛处,因此尽管试验4的计算结果比试验3有所改进,但改进比较有限;试验5同时同化了T/P资料和验潮站资料,两种资料相互补充,使计算结果相比试验3与试验4有最大的改进.此外,从表3还可以看到,试验4和试验5计算结果与同化数据的偏差都比试验3大.原因可能是试验4和试验5同化的数据里都有验潮站资料,这些资料点位置对应的水深一般较浅,因此要在整体上拟合较为困难;还有一个原因可能是浅水区域潮汐变化快,浅水分潮在潮高误差影响中占较大比重,如果采用包含浅水分潮在内的多分潮伴随同化研究,应该能改善整体的拟合效果.为了更好地说明模式的计算结果,这里将本研究偏差和部分已有数值模式给出的偏差作了比较.虽然不能单纯地以不同研究者提供的计算偏差来判定模式的优劣,但可以作为参考.表4列出了若干已有模式偏差与本研究偏差的比较,表5从相应的文献中找出15个相同的验潮站点的调和常数,进行数值计算结果检验,给出了本研究与其他两个模式M2分潮的逐点偏差比较.因为国内研究者多用振幅和迟角的绝均差表示,所以表4、表5仅给出了绝均差的比较.由于赵保仁等[7]和万振文等[9]对全日潮只计算m1=[(K1+O1)/2]分潮,表4将其偏差换算至K1分潮,即振幅偏差乘112,偏角不变.从表4可以看出,M2分潮与万振文等[9]相比,振幅差相关不大,但迟角差大为改善,与其他3个模式相比,振幅差和迟角差均有显著改进;K1分潮振幅差叶安乐等[8]的最小,为116cm,本研究为218cm,迟角差最小的是万振文等[9]的模式,为414b,本研究为418b;与万振文等[10]4个分潮的偏差相比,本研究计算结果与观测值的拟合程度均有显著提高.从表5的对154期丘仲锋等:黄海、渤海T O PEX/Poseidon高度计资料潮汐伴随同化比可以看出,偏差值小的点,本研究明显多于其他模式;M2分潮的最大振幅差和最大迟角差,本研究均小于用以比较的两个模式.通过比较看出,本研究计算结果与观测值的符合情况较上述模式整体上有提高.对于M2及K1分潮的振幅差与其中某些模式相比改进不大的原因,可能是本研究模式启动时开边界采用的是/静边界0,即当t=0时,a=b=0,而上述模式则是给定开边界值.如果在伴随同化中采用上述模式的开边界值作为初始的开边界先验估计,相信可以进一步改进模式模拟的结果.由试验5得到的M2,S2,O1,K1等4个分潮的等振幅线和同潮时线如图2所示.图2计算等振幅线(虚线)(单位:cm)和同潮时线(实线)[单位:(b)](相对于东八时区)分布31211半日分潮从图2a,b可以看出,研究区域M2,S2分潮有4个无潮点,南黄海无潮点和秦皇岛无潮点位置较稳定,黄河口无潮点较不稳定,方国洪[6]认为黄河口岸形的变迁影响无潮点位置,本研究支持他的观点.和陈国珍[19]给出的同潮图相比,等潮时线的走向较为一致;等振幅线分布基本一致,但在南黄海区域及朝鲜半岛西南近海,振幅相对较大.31212全日分潮从图2c,d可以看出,研究区域O1,K1分潮只有一个无潮点.根据文献显示,O1,K1分潮在34b N附近有一个无潮点.因为本研究中将开边界取在34b N的缘故,虽然在同潮图上看不到明确的无潮点,但在图2c,d上都可以明显地看到34b N附近有封闭的低值区.因此,图2c,d给出的同潮图只有一个无潮点,是因为计算区域不够大,无法将南黄海无潮点完全显示出来.与陈国珍[19]给出的同潮图相比,O1分潮除了北黄海270b等潮时线走向有较大差别,等潮时线吻合得较好,振幅在南黄海区域及朝鲜半岛西南近海偏小;K1分潮则较为一致.16海洋学报27卷以上结果讨论及分析表明,采用伴随同化方法,同时同化T /P 资料与沿岸验潮站资料,能有效地反演黄海、渤海的潮波系统,为高度计资料的浅海应用打下良好的基础.两位审稿专家对本文提出宝贵的修改意见,本文中使用的T/P 卫星资料由美国NASA JPL 物理海洋数据存档中心(PODAAC)提供,谨致谢忱.参考文献:[1] AN H ui Soo.A num erical exp eriment of M 2tide in th e Yellow Sea[J].Jour nal of the Oceanograp hical S ociety of Japan,1977,33(2):103-110.[2] 沈育疆.东中国海潮汐数值计算[J].山东海洋学院学报,1980,10(3):26-35.[3] 夏综万,王钟桾.黄海M 2分潮的数值模拟[J ].黄渤海海洋,1984,2(1):10-18.[4] FANG Gu o -hong.T ide and tidal curren t ch arts for the m arginal seas adjacent to China[J].Chinese Journ al of Oceanology and Limn olo -gy,1986,4(1):1-16.[5] 丁文兰.渤海和黄海潮汐潮流分布的基本特征[J].海洋科学集刊,1985,25:27-39.[6] 方国洪,杨景飞.渤海潮运动的一个二维数值模型[J].海洋与湖沼,1985,16(5):337-346.[7] 赵保仁,方国洪,曹德明.渤黄东海潮汐潮流的数值模拟[J ].海洋学报,1994,16(5):1-10.[8] 叶安乐,梅丽明.渤黄东海潮波数值模拟[J].海洋与湖沼,1995,26(1):64-70.[9] 万振文,乔方利,袁业立.渤、黄、东海三维潮波运动数值模拟[J ].海洋与湖沼,1998,29(6):611-616.[10] 万振文,袁业立,乔方利.东中国海多分潮运动数值模拟[J ].水动力学研究与进展(A 辑),1999,14:82-91.[11] LEFEVRE F,LE PROVOST C,LYARD F H.How can w e imp rove a global ocean tidal model at a regional scale?A test on the YellowS ea and East China Sea[J ].J ou rnal of Geophysical Research,2000,105(C4):8707-8725.[12] LARDER R W.Optim al control of open b oundary con dition s for a numerical tidal model[J].Computer M ethods in Applied M echanicsand Engineerin g,1993,102:367-387.[13] S EILER U.Estimation of open boun dary conditions with th e ad joint m ethod[J].Journ al of Geophysical Research ,1993,98:22855-22870.[14] 朱 江,曾庆存,郭冬建,等.利用伴随算子法从岸边潮位站资料估计近岸模式的开边界条件[J ].中国科学(D 辑),1997,27(5):462-468.[15] 吕咸青,方国洪.渤海开边界潮汐的伴随法反演[J].海洋与湖沼,2002,33(2):113-120.[16] 吕咸青,吴自库,殷忠斌,等.渤、黄、东海潮汐开边界的1种反演方法[J].青岛海洋大学学报,2003,33(2):165-172.[17] 韩桂军,方国洪,马继瑞,等.利用伴随法优化非线性潮汐模型的开边界条件[J].海洋学报,2001,23(2):25-30.[18] 吴自库,田纪伟,吕咸青,等.南海潮汐的伴随同化数值模拟[J ].海洋与湖沼,2003,34(1):101-108.[19] 陈国珍.渤海、黄海、东海海洋图集(水文)[M ].北京:海洋出版社,1992.429-454.Tidal adjoint assimilation with the TOPEX/Poseidon altimetrydata in the Huanghai and Bohai SeasQIU Zhong -feng 1,2,H E Y-i jun 1,L B Xian -qing 3(1.I nstitute of Oceanology ,Chinese A cademy of Sciences ,Qing dao 266071,China;2.Graduate School o f the Chinese A cademy of S ciences ,Beij ing 100039,China;3.Ocean Univers ity o f China,Qing dao 266003,China)Abstract:Tidal co nstituents ,w hich are ex tracted from about 10-y ear T OPEX/Poseidon altimetry data a -long track points,are assimilated into an tw o -dimension non -linear tidal m odel w ith the adjo int assim ilation method.Tidal constituents M 2,S 2,O 1and K 1are simulated in the H uang hai and Bohai Seas,then Co -tidal maps are o btained.The bias betw een calculatio ns and observations are also pro vided and the results tellthat calculations agree w ell w ith observations.Mo reo ver,there are tw o exper im ents in this paper:one is to optim ize differ ent coefficient,the other is to assimilate different data.T he fir st ex perim ent show s that it is better to optim ize the o pen boundary conditions and the bottom friction coefficients at the sam e time than174期 丘仲锋等:黄海、渤海T O PEX/Poseidon 高度计资料潮汐伴随同化18海洋学报27卷to optimize the open boundary conditio ns only.T he second ex periment show s that the results are better w hen assim ilating the altimetr y data and the tidal co nstants from tidal g auges simultaneously than assimila-ting any one separately.T he results also show that it is useful to simulate the tidal system in the H uanghai and Bohai Seas by assim ilating T/P altimetry data and tidal co nstants from tidal g auges w ith adjoint assim-i lation m ethod.Key words:adjoint assimilatio n;T OPEX/Poseidon altim etry;tide。
海洋知识
1. 知识点：船籍港是船舶注册的所在港口，例如天津港。
船籍港所属国家为船旗国。
如果船籍港是天津，则船旗国就是中国。
按规定，远洋货轮到港时应悬挂相应的旗子。
船尾挂船旗国国旗，且每天早晚应该升降旗，进入港口要尽量换全新船旗国国旗。
如果船舶所到港口为非船旗国港口，则在驾驶台顶部的主桅杆右侧升起所抵达港口国国旗以示尊重，大小要与船旗国国旗一致，左侧根据不同情况悬挂检疫旗、引水旗等；如果船舶所到港口是船籍港本国港口，则主桅杆不需要升国旗。
2. 知识点：国际上的1桶油是42加仑，换算成升为158.98升，大约135千克。
3. 知识点：海里是航海上度量距离的单位，它等于地球椭圆子午线上纬度1分所对应的弧长。
由于地球子午圈是一个椭圆，它在不同纬度的曲率是不同的，因此，纬度1分所对应的弧长是不相等的，不同国家的海里标准长度也是不同的。
1海里=1.852千米，这是目前我国和世界上大多数国家采用的1929年国际水文地理学会议通过的海里的标准长度。
1海里=1.851千米是美国标准；1海里=1.853千米是法国标准；1海里=1.854千米是英国标准。
4. 知识点：节（Kn）以前是船员测船速的，每走1海里，船员就在放下的绳子上打个节，后来就用“节”做船速的单位。
1节（kn）=1海里/时=（1852/3600）m/s。
一般来说，VLCC的最高航速的理论值可达到16节。
延伸：“迈”通常是用作汽车的速度单位。
“马赫”则是火箭、航空飞船的专用速度单位。
5. 知识点：世界气象组织定义：中心持续风速在12级至13级（即每秒32.7米至41.4米）的热带气旋为台风（typhoon）。
国际惯例依据其中心附近最大风力分为：强热带风暴：最大风速10～11 级；台风：最大风速12～13级；强台风：最大风速14～15级；超强台风：最大风速≥16级。
6. 知识点：VLCC是超大型油轮“Very Large Crude Carrier”的英文缩写，载重量一般为20至30万吨，相当于200万桶原油的装运量。
卫星海洋学_河海大学中国大学mooc课后章节答案期末考试题库2023年
卫星海洋学_河海大学中国大学mooc课后章节答案期末考试题库2023年1.皮层深度是电场强度衰减为初始值的1/e所在的深度？参考答案:正确2.在可见光和近红外波段范围内，电磁波的穿透深度大小随波长的增加而减小。
参考答案:错误3.大气气溶胶是大气中的微量成分，在大气辐射收支平衡和全球气候模式中扮演着重要角色，在水色遥感的大气校正中也扮演着重要角色。
气溶胶对电磁辐射的影响包括以下哪些方面？参考答案:直接散射和吸收电磁辐射_作为凝结核，在大气中改变云滴的浓度和云滴在大气中存在的时间，通过云滴影响电磁辐射_可以把太阳辐射向太空中散射，造成衰减_可以吸收由地面而来的长波辐射，其作用与温室气体作用相似，形成增益4.如果使用5.3GHz的C波段散射计，当入射角是45°时，水面上波长多大的水波与入射的电磁波共振？参考答案:4 cm5.已知σ0[dB]=10 log10(σ0)，如果σ0增加到原来的100倍，σ0[dB]增加多少？参考答案:206.复折射率包含了实部和虚部，下面哪些参数与复折射率的虚部有关系？参考答案:衰减系数_皮层厚度_穿透深度7.下面哪个深度参量常用于描述热红外和微波辐射在海水表层的电磁波衰减作用？参考答案:皮层深度8.与可见光和热红外波段相比，微波波段的海水发射率相对较（）。
参考答案:低##%_YZPRLFH_%##小9.黑体的发射率等于1，所有非黑体的发射率都小于1。
参考答案:正确10.我国海洋二号（HY-2）系列卫星上未装载以下哪个微波传感器？参考答案:SAR11.热红外波段传感器测量海面亮温的理论依据是遵循哪一个定律？参考答案:普朗克定律12.按照目标的能量来源，遥感可分为哪几类？参考答案:被动遥感_主动遥感13.SAR可以探测以下哪些海洋要素或过程？参考答案:海洋上升流_海面风场_海洋内波14.在水色遥感中，二类水体的光学成分比较复杂，通常是指具有较高的（___）含量的水体。
【国家自然科学基金】_m2分潮_基金支持热词逐年推荐_【万方软件创新助手】_20140801
1 2 3 4 5 6 7 8 9
2014年 科研热词 验潮站 西北太平洋 潮能通量 潮能耗散 大洋潮汐模式 卫星高度计 南海 中国近海 fvcom模式 推荐指数 1 1 1 1 1 1 1 1 1
2012年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
科研热词 潮流 潮汐 趋近法 调和常数 调和分析 西北太平洋 烟墩角 海流 椭圆要素 日本海 季节变化 垂向结构 南海陆架陆坡区 准调和分析 fvcom
推荐指数 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1
2014年 科研热词 推荐指数 序号 潮汐 3 渤海 3 调和常数 2 潮流 2 m2分潮 2 efdc 2 黄海 1 验潮站 1 陆架坡折区 1 长江河口 1 长江口 1 负荷潮 1 西南黄海 1 自吸-负荷潮 1 理论最高潮面 1 理论最低潮面 1 环流结构 1 潮高 1 潮波 1 混合系数 1 深度基准面 1 海洋潮汐 1 泰国湾 1 有限体积法海洋数值模型(fvcom) 1 方差法 1 数值模拟 1 强潮狭长海湾 1 底摩擦 1 底拖曳系数 1 垂直涡黏性系数 1 垂直分辨率 1 地形效应 1 变化趋势 1 双矩形海湾模型 1 历元 1 印尼海 1 南海北部 1 动量平衡法 1 分潮 1 东海 1 三维斜压模型 1 taylor问题 1 nao.99b潮汐模型 1 m2分潮流 1 lagrangian粒子追踪 1 fvcom(finite-volume coastal ocean 1 model) fvcom 1 csr4.0潮汐模型 1
推荐指数 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1
2011年序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
西北太平洋柔鱼资源丰度时空分布的GAM模型分析
２０ ０６年 １ ２月
［ 文章 编号 ］１０ — ４５ ２０ ）４ ０９ ０ ０７ ７０ （０６ ０ — ２５— ６
西 北太 平洋 柔鱼资 源丰度 时空分布 的 Ｇ Ｍ 模 型分析 Ａ
陈新 军 ，田思泉
（ 海水产大 学海洋学院 ，上海 ２０ ９ ） 上 ０ ０ ０
［ 要 ］根据 １９ ２０ 年我 国在西北 太平 洋海 域 柔鱼 生产统 计及相关 数据 ，利 用 Ｇ Ｍ 模 型分析 了 摘 ９ ６～ ０ １ Ａ 表温 、月份 、经纬 度等因子对柔 鱼资 源丰度 Ｃ Ｕ Ｐ Ｅ的影响．研 究认为 ，经纬度 、月份 和表温对 Ｃ Ｕ Ｐ Ｅ时空 分布都有较大 的影 响．１０。 ６ Ｅ以西海域 Ｃ Ｕ Ｐ Ｅ高 ，而 １５。 以东海域 低 ，并主 要集 中在 ４ 。 ６ Ｅ ０ Ｎ～４ Ｎ海 ３。 域．８～１ ０月 Ｃ Ｕ Ｐ Ｅ为最大．不同海域柔鱼 分布的适 宜表 温不 相同，１０。 ５ Ｅ以西海域为 １ ３～１ ８℃ ，１０。 ５ Ｅ
一
１５。 ６ Ｅ海域 为 １ ４～１ Ｃ，１５。 ８ｏ ６ Ｅ～１０。 ８ Ｅ海域为 １ ～１ １ ４ｏ Ｃ．
［ 关键词 ］柔鱼资源 ；影响因子 ；Ｇ Ｍ 模型 ；西北太 平洋 Ａ ［ 中图分类号 ］Ｓ９ １ ３ ［ 文献标识码 ］Ａ
Ｏ 引 言
柔鱼 （ ｍ ａｔｐ ｅ ｂ ｒａ ｉ 广泛分布在北：平洋海域 ，是 目前世界头足类渔业中最重要的种 Ｏ ｍ ｓ ｅｈｓ ａｔ ｍ ） ｒ ｒ 太 类之一．我 国于 １９ 年开发了该资源，１９ ９３ ９４年进行规模性利用 ，目前年产量基本上稳定在 １ １ ０× ０ ｔ 左右 ，已成为我 国鱿钓渔船重要的捕捞对象之一．分析和研究柔鱼资源状况的变动及其影响因素是确
［ 收稿 日期 ］２ ０ ０ ２ ０５— ６— ８ ［ 基金项 目］上海市启 明星计划项 目 （ ２ Ｃ１００ ；上海市重点学科资助项 目 （ １０ ） ０Ｑ ４ ） ４ Ｔ １１ ［ 作者简介 ］陈新军 （ ９ ７一） １６ ，男 ，教授 ，博导 ，从事渔业资源与渔场学及渔业资 源经济学研究．
基于ADCIRC水动力模型的台风风暴潮预报模型验证
个潮汐测站的实测水位过程与模型预报的潮位过程结果。结果表明：①无论是在天文潮期间还是风暴潮
期间，各பைடு நூலகம்站的逐时潮位预报平均绝对误差均在 30 cm 以内；②在对天文潮进行验证时，除崇武站的高潮位
相对误差达到 50 cm 以外，其余各测站的高潮位相对误差均在许可范围±20 cm 以内；③在对风暴潮过程进
洋海域的风暴潮预报模型，并对该模型的合理性进
收稿日期：2018-12-12 基金项目：水利部公益性行业科研专项经费项目（201501014） 作者简介：刘克强，男，高级工程师，主要从事水利前期规划与调度方面的工作。E-mail：15117438@
·27·
2019 年 4 月
水 利 水 电 快 报 EWRHI
为了满足西北太平洋海域复杂岸线和分辨率
的要求，数值模拟采用了目前较为成熟、应用广泛 的 ADCIRC 水动力模型。该模型的构建是基于有限 元的方法，于 1992 年由美国圣母大学的 J.J.Wester⁃ ink 和美国北卡罗来纳大学的 R.A.Luettich 共同研 制。目前，ADCIRC 水动力模型已被广泛应用于河 口 、海 洋 、海 岸 等 区 域 的 潮 汐 、海 流 和 风 暴 潮 预 报 中。本文利用该模型构建了一套适用于西北太平
表 1 2010~2015 年由台风风暴潮灾害导致的经济损失
年份
2010 2011 2012 2013 2014 2015
台风风暴潮 过程次数/次
10 9 13 14 5 6
风暴潮致灾 次数/次 7 5 9 11 5 6
直接经济损失/亿元
65.79 48.81 126.29 152.45 134.69 72.18
2.2 ADCIRC 风暴潮数值预报模型
渤、黄、东海8个主要分潮的数值模拟研究
渤、黄、东海8个主要分潮的数值模拟研究宋泽坤;俞亮亮;向芸芸;施伟勇;许雪峰;杨万康;潘冲【摘要】A high resolution grid model is established for the calculated areas including the Bohai Sea,the Yellow Sea,the East China Sea,the East China Sea continental shelf and the Ryukyu Islands by using the MIKE software and a nonstructural triangular mesh,in which the actual water depth and shoreline are considered and the open bound-ary is resulted from the comprehensive Northwest Pacific tidal model.Based on this model,the tidal wave processes in the China seas are simulated,a harmonic analysis is made for the vertical motion processes of the tidal waves and thus the propagation and distribution characteristics of eight tidal components like M2 ,S2 ,K1 ,O1 ,N2 ,K2 ,P1 and Q1 in the Bohai Sea,the Yellow Sea and the East China Sea are obtained.The re-sults from the model have been validated by using the harmonic constants of 14 tide sta-tions located along the coast of China,indicating that the model is more accurate and re-liable.The research results show that the propagation of the 4 maj or semidiurnal tides (also the diurnal tides)are similar in the above seas,that means,the propagation prop-erty of the tidal waves along the coasts of the above seas is similar to that of the coastal Kelvin Wave.Moreover,4 amphidromic points of the semidiurnal tide component and 2 amphidromic points of the diurnal tide component in the calculated areas are reproduced successfully.The diurnal tidal amplitude is minimal around the amphidromic points and maximal at the free zone of thebay.The tidal range amplitude is large towards the coast and small towards the sea in the East China Sea.It is large along the coasts of Zhej iang and Fuj ian and relatively small in the Yellow Sea.In the Bohai Sea,it is maximal at the tops of the Bohai Bay and the Liaodong Bay and smaller in the area around the two am-phidromic points of the diurnal tide component.%应用MIKE数值模拟软件,采用无结构三角形网格,建立一套计算区域包括整个渤海、黄海、东海以及东海大陆架和琉球群岛的高分辨率数值模型,考虑了实际水深和岸线,外海开边界采用西北太平洋大模型结果的潮位提供,模拟了东中国海潮波的波动过程,对潮波垂直运动过程进行调和分析,得到了渤海、黄海、东海的 M2,S2,K1,O1以及 N2,K2,P1,Q1八个主要分潮的传播和分布特征。
【国家自然科学基金】_数值预报模式_基金支持热词逐年推荐_【万方软件创新助手】_20140801
能见度 背景误差协方差 网格嵌套 综合决策算法 线性和非线性理论 粗糙集 窗宽 空气质量数值预报 积云对流参数化 秦岭山脉的作用 短波辐射效应 短期气候预测 相变 相关分析 相似误差订正 登陆台风 环流结构 环境预报 特大暴雨 热带气旋 灰霾 潮波变化 潮汐预报 潮汐数值模式 漂移物运动轨迹 滑动分区 湿物理参数化 温室效应 海面风场 海面空气动力粗糙度 海浪数值模式 海洋环流模式 海平面上升 海啸 浅水方程模式 波致混合 河南省 沙尘输送 沙尘气溶胶 水物质 水汽通量辐合 水汽输送 水位变化 气候预测 气候模式 气候 正规模 欧拉网格点 模糊神经网络 模拟研究 模拟的水汽场特征 模式模拟 概率预报 核甬数
科研热词 推荐指数 数值模拟 15 资料同化 5 数值天气预报 5 数值预报 4 grapes 4 可预报性 3 青藏高原 2 集合预报 2 暴雨 2 径向风 2 年代际变化 2 对流参数化 2 大气环流模式 2 大暴雨 2 地形 2 台风 2 初始扰动 2 中尺度系统 2 wrf模式 2 mm5模式 2 arps模式 2 飞机颠簸 1 飑线 1 风场分离和重建 1 风切变 1 颠簸指数 1 预条件子 1 预报误差 1 预报因子 1 非静力全可压大气运动动力框架 1 非静力 1 非对称风场订正加强 1 非对称风场 1 静止锋 1 雷达反射率 1 集合卡尔曼滤波(enkf) 1 集合卡尔曼变换(etkf) 1 集合kalman滤波同化 1 集合kalman滤波 1 降雨预报 1 降水预报 1 陆气通量 1 通量 1 递推参数 1 逆温层 1 适应性观测 1 迁入 1 边界层约束 1 起转 1 赫姆霍兹方程 1 资料逐时同化 1 资料 1
2009年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52
温盐流预报讲义
2022年3月7日
目录
概况 温盐流基础知识 全球和区域温盐流时空特征 全球和区域温盐流数值预报 近岸精细化三维温盐流数值预报 海洋中尺度现象诊断预报
1. 概况
1.1 国际温盐流预报概况
系统名称
美国海军 GOFSv3.0
美国NCEP RTOFS
海洋模式 HYCOM
HYCOM
模式区域 全球 全球
3. 全球和区域温盐流时空特征
3.1 全球
➢盐度垂向分布
✓在赤道海区表层较浅的低盐水下，是由南、 北半球副热带海区下沉后向赤道方向扩展 的高盐水，被称为大洋次表层水，具有大 洋垂向上的最高盐度；
✓在高盐次表层下，是由南、北半球中高纬 度表层下沉的大洋（低盐）中层水；
✓在低盐中层水之下，充满了在高纬海区下 沉形成的深层水与底层水，盐度稍高。
度和盐度作误差订正
NEMOVAR (3DVar-FGAT)
14 days
7days
51成员集合预报，18days (每天), 32 days (每周
), 7months (每月)
大气强迫
海军全球环境模式 NAVGEM
NCEP GFS 3-hourly, HWRF (Hurricane WRF)
ECMWF 3hourly 分析场和预报场
全球海洋表层温度年平均分布(WOA13 V2, 1955-2012)
3. 全球和区域温盐流时空特征 3.1 全球
➢一般来说，同一海区夏季的海温要比冬季偏高，特别是在北半球高纬 度地区。
全球海洋表层温度夏季（7-9月）和冬季（1-3月）分布
3. 全球和区域温盐流时空特征
3.1 全球
➢海温深层分布：
目录
概况 温盐流基础知识 全球和区域温盐流时空特征 全球和区域温盐流数值预报 近岸精细化三维温盐流数值预报 海洋中尺度现象诊断预报
我国海洋气象数值预报业务发展与思考
我国海洋气象数值预报业务发展与思考黄彬;阎丽凤;杨超;徐晶【摘要】近年来，我国海洋气象数值预报业务有了一定的进展，国家级和沿海省台建立了海雾、海浪、风暴潮等一些海洋气象数值预报模式。
通过对我国沿海省台业务应用中的海洋气象数值预报模式的调查，概述了我国现行海洋气象数值预报业务模式的发展现状，分析了海洋气象数值预报目前存在的问题和面临的挑战，并探讨了未来的发展方向。
%In recent years, China's marine meteorological numerical prediction has achieved some progress. National-level cities and coastal provinces have established a number of sea-fog, wave, storm-tide-ranging marine meteorological numerical prediction models. Through the investigation of marine meteorological numerical prediction models in China's coastal provinces and their application, this paper summarizes the development status of the Chinese marine meteorological numerical prediction model, analyses the existing marine meteorological numerical prediction problems and challenges, and discusses the future development and suggestion of marine meteorological numerical prediction.【期刊名称】《气象科技进展》【年(卷),期】2014(000)003【总页数】5页(P57-61)【关键词】海洋气象;数值预报;业务发展【作者】黄彬;阎丽凤;杨超;徐晶【作者单位】国家气象中心，中国气象局，北京100081;山东省气象局，济南250000;国家气象中心，中国气象局，北京100081;国家气象中心，中国气象局，北京100081【正文语种】中文我国海洋气象预报业务经过近几年的发展,取得了长足进步,已经初步建立了国家级、海洋区域中心级、省级共3级的海洋气象预报预警服务业务体系[1-4]。
基于FVCOM的泉州湾海域三维潮汐与潮流数值模拟
基于FVCOM的泉州湾海域三维潮汐与潮流数值模拟林作梁;朱学明;鲍献文;刘钦政【摘要】Based on a finite-volume coastal ocean model (FVCOM) , adopting an unstructured triangle grid, a three-dimensional tide and tidal current numerical modeling with high resolution (26 m) is applied to Quanzhou Bay. The simulated results agree well with the observed data from two tide-gauges and three continuing current stations , and reproduce the distribution features of the tide and tidal currents in the Quanzhou Bay famously. The distributions of co-tidal charts and tidal ellipses on the surface layer for four major constituents (M2 , S2, K1 ,O1) are obtained. What's more, the distributions of the maximum probable tidal range and tidal currents velocity and tidal residual currents on the surface and bottom layers are obtained, too. By analyzing, the maximum tidal amplitude and phase-lag range for the four constituents are 219 cm and 19°,85 cm and 25°,26 cm and 12°,26 cm and 9°, respectively. The tidal wave is anti-clockwise standing wave in the east area of Shihu Port, but it is advancing wave in the west area of Shihu Port. The maximum probable tidal range increases from 8. 0 m at the mouth of the bay to 8.8 m inside of the bay. The type of tidal currents is regular semi-diurnal currents inside of the bay, and the maximum velocity of ebbing is larger than flooding. The velocity in the channel of Beiwujiao is stronger than any other area, and the maximum probable tidal-current velocity is 2. 4 m/s on the surface layer. The flow is rotating with anti-clockwise mainly at the mouth of thebay, while rectilinear flow is mainly inside of the bay, such as estuaries and channels. And the directions of major axes are along with channels direction mainly, or paralleling with isobaths and coastlines. The maximum tidal current velocities for the four major constituents are 1. 4 m/s,0. 58m/s, 0. 12 m/s,0. 10 m/s on the surface layer, respectively. Tidal residual currents velocity is closely related to tidal currents, the maximum velocities on the surface layer, the middle layer and the bottom layer are 26 cm/s, 20 cm/ s, 16 cm/s, respectively. All of them are coming into the bay from north and going out of the bay from south.%基于FVCOM海洋数值模式,采用非结构的三角形网格和有限体积法,建立了泉州湾海域高分辨率(26 m)的三维潮汐、潮流数值模型.模拟结果同2个验潮站和3个连续测流站的观测资料符合良好,较好地反映了泉州湾内潮汐、潮流运动的变化状况和分布特征,给出了M2、S2、K1、O14个主要分潮的同潮图、表层潮流椭圆分布,以及模拟区域内最大可能潮差、表层最大可能潮流流速和潮余流分布.分析表明,4个分潮的最大潮汐振幅和迟角差分别为219 cm和19°,85 cm和25°,26 cm和12°,26 cm和9°；石湖港以东海域的潮波为逆时针旋转的驻波,以西海域为前进波；最大可能潮差由湾口的8.0m向湾内增加至8.8m.湾内潮流类型为规则半日潮流,落潮最大流速大于涨潮最大流速,北乌礁水道为强流区,表层最大可能潮流流速为2.4 m/s;湾口潮流运动以逆时针方向的旋转流形式为主,湾内的潮流运动以往复流形式为主,长轴走向主要沿着水道方向,与等深线和海岸线平行；四个分潮流表层最大流速分别为1.4 m/s,0.58m/s,0.12 m/s,0.10 m/s.余流流速大小与潮流强弱有密切的联系,表、中、底层最大余流流速分别为26 cm/s,20 cm/s,16 cm/s,三者在水平方向基本呈北进南出的分布形态.【期刊名称】《海洋学报（中文版）》【年(卷),期】2013(035)001【总页数】10页(P15-24)【关键词】泉州湾;潮汐;潮流;FVCOM;潮余流【作者】林作梁;朱学明;鲍献文;刘钦政【作者单位】泉州市环境监测站,福建泉州362000;国家海洋环境预报中心国家海洋局海洋灾害预报技术研究重点实验室,北京100081;中国海洋大学海洋环境学院,山东青岛266100;国家海洋局海洋减灾中心,北京100194【正文语种】中文【中图分类】P731.231 引言泉州湾地处福建省东南沿海，台湾海峡西侧沿岸的中部，海域总面积211.24 km2，包括围垦面积45.70 k m2，滩涂面积84.84 k m2；海岸线总长229.61 km［1］。
海洋数值预报_国防科技大学中国大学mooc课后章节答案期末考试题库2023年
海洋数值预报_国防科技大学中国大学mooc课后章节答案期末考试题库2023年1.在三维海洋流场模式中，在哪些区域的关键小尺度过程需要做参数化处理？答案:2.以下哪些物理过程不发生在深水条件下答案:3.以下哪个模式是最先提出的第三代海浪数值模式？答案:4.大洋潮汐模式和近海潮汐模式一样。
答案:错误5.在开展海洋环流数值模拟工作时，必须要精确的气象模式提供的预报场资料。
答案:错误6.理想并行效率为1，表明全部处理器都在满负荷工作。
答案:正确7.海浪预报的方法有答案:半经验半理论的方法_海浪数值预报方法_经验统计预报方法8.三波相互作用使得高频能量向低频转移。
答案:错误9.海浪数值预报产品中哪些？答案:平均周期_有效波高_平均波向10.越是复杂的湍封闭模型就是越好的模型。
答案:错误11.在σ坐标系，斜压梯度力误差产生的根本原因是由两个大项代数和去求小项，很容易产生与小项本身量级相当的计算误差。
答案:正确12.采用ADI格式离散二维流场方程组，能够使得数值计算简化，保证_____计算稳定。
答案:重力外波。
13.以下关于矩形网格的说法错误的是答案:矩形网格使得真实的海岸上法向流速为零容易满足。
14.底摩擦在谱峰附近耗散能量最多。
答案:正确15.风暴潮数值预报是“数值天气预报”和“风暴潮数值计算”二者组成的统一整体。
答案:正确16.梯度风原理风场模型是先算出气压，然后通过气压梯度与梯度风的关系得到风速。
答案:正确17.关于矩形网格的说法错误的是答案:18.下列哪些是四维变分同化优点？答案:19.二维流场基本方程中最容易造成计算不稳定的是重力外波项。
答案:20.C网格相比E网格相对较好的原因是答案:C网格相对E网格在边界处理时更加方便。
_C网格相对E网格有更小的计算量。
21.水平湍流摩擦项表达为平均速度梯度和粘性系数的乘积。
答案:正确22.物理过程参数化不会给海洋数值预报带来误差和不确定。
答案:错误23.各类二维网格产生的数值效应均相同。
【国家自然科学基金】_潮差_基金支持热词逐年推荐_【万方软件创新助手】_20140731
推荐指数 3 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2010年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39
2008年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
科研热词 磨刀门水道 盐水入侵 杭州湾 咸界运动规律 鸭绿江口 边界层过程 输沙平衡 辐射沙洲 调和常数 调和分析 设计流速 航道整治 相关分析 盐水 潮流深槽 潮波 潮汐调和分析 潮位 涌潮数学模型 治导线放宽率 河口 沙波 江苏省 水深测量 水位 最大可能潮差 数值模拟 挟沙力公式 成界运动规律 悬沙输运 悬沙浓度(ssc) 小波分析 太平洋 大洋潮波模式 嘉绍大桥 台湾浅滩 再悬浮 topex/poseidon卫星 ions格式 adcp
推荐指数 3 2 2 2 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
2009年 序号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42
西北太平洋的潮能通量与耗散
热带海洋学报JOURNAL OF TROPICAL OCEANOGRAPHY2014年第33卷第1期:1−9西北太平洋的潮能通量与耗散朱学明1,2,宋德海3,鲍献文2,刘桂梅11.国家海洋环境预报中心国家海洋局海洋灾害预报技术研究重点实验室,北京100081;2.中国海洋大学海洋环境学院,山东青岛266100;3.中国海洋大学物理海洋教育部重点实验室,山东青岛266100摘要:利用基于FVCOM(Finite Volume Coastal Ocean Model)模式同化沿岸验潮站数据的高分辨率潮汐数值模型结果,分析研究了包含中国近海、日本海和鄂霍次克海在内的西北太平洋海域全日、半日分潮的潮能通量与耗散。
西北太平洋的潮波能量分3支,分别传入鄂霍次克海、东海和南海。
传入东海的半日潮波能量是传入南海的3倍左右;传入南海的全日潮波能量是传入东海的5倍多。
传入中国东部海域的M2分潮能中,有64.3%耗散在东海,32.4%耗散在黄海,仅有3.2%耗散在渤海;而K1分潮能中分别有48.2%、31.4%及7.1%耗散在东海、黄海及渤海。
进入南海的潮能中,仅有32.7%的M2分潮能和38.3%的K1分潮能耗散在南海的北部海域,另有23.9%的M2分潮能进入并耗散在台湾海峡,其余则进入南海南部。
传入日本海的太平洋潮能很少。
鄂霍次克海消耗的全日潮能是半日潮能的2倍。
关键词:西北太平洋;中国近海;潮能通量;潮能耗散;FVCOM模式中图分类号:P731.23文献标识码:A文章编号:1009-5470(2014)01-0001-09Tidal energy flux and dissipation in the Northwest PacificZHU Xue-ming1,2,SONG De-hai3,BAO Xian-wen2,LIU Gui-mei11.Key Laboratory of Research on Marine Hazards Forecasting,National Marine Environmental Forecasting Center,Beijing100081;2.College of Physical and Environment Oceanography,Ocean University of China,Qingdao266100,China;3.3.Key Laboratory of Physical Oceanography,Ministry of Education,Ocean University of China,Qingdao266100Abstract:A Finite V olume Coastal Ocean Model(FVCOM)that assimilated tidal gauge data was used to study the tidal energy flux and dissipation in the Northwest Pacific,including the China seas,the Japan Sea(JS)and the Sea of Okhotsk(SO). Most of the tidal energy in the Northwest Pacific is transmitted to the East China Sea(ECS),the South China Sea(SCS)and the SO.The semi-diurnal tidal energy transported into the ECS is three times larger than that transported into the SCS;whilst the diurnal tidal energy transported into the SCS is five times of that transported into the ECS.In the waters east of China, 64.3%of M2semi-diurnal tidal energy is dissipated in the ECS,32.4%in the Yellow Sea(YS),and only3.2%in the Bohai Sea (BS).For K1diurnal tide,48.2%,31.4%and7.1%of the tidal energy are dissipated in the ECS,the YS and the BS, respectively.Among the tidal energy transported into the SCS,32.7%of M2and38.3%of K1are dissipated in the northern region,23.9%of M2is transported into and dissipated in the Taiwan Strait,and the rest enters the southern SCS.Little tidal energy is transported from the Pacific Ocean into the JS.The dissipation of diurnal tidal energy in the SO doubles that of the semi-diurnal tidal energy.Key words:Northwest Pacific;China seas;tidal energy flux;tidal energy dissipation;Finite Volume Coastal Ocean Model(FVCOM)有关海洋潮能通量与耗散的最早研究是由Taylor(1919)开始的，他利用潮能底边界耗散公式计收稿日期:2013-01-04;修订日期:2013-03-26。
m海洋科技(知识点)
1．知识点：人类在认知地球的历史过程中使用过三种平台：最早是地面与海面观测平台，其次是空间观测平台，第三个观测平台是海底。
相对于前两个观测平台，海底观测平台技术虽然最难，但是由于海洋占地球面积的70.8%，而深海大洋占据海洋面积的92.4%，且深海海底也是人类认识地球的窗口，不过人类对其又知之甚少，故发展海底观测平台，其科学意义更显重大。
2．知识点：近年来，随着各国对天然气水合物研究的逐步深入，利用地震方法来探测水合物已经成为了大家公认的最为有效的方法。
海洋地震勘探自20 世纪30 年代中期才开始开展。
最初，除所用勘探设备的防水和水密性措施之外，其余设备和方法都是将陆地地震勘探照搬到了海洋上，由于设备和方法的限制，最初的海洋地震勘探主要是集中在濒临陆地的浅水区域。
到了20 世纪50 年代，随着技术的进步，地震勘探接收装置采用了晶体检波器，采用了将光点式地震仪安放在观测船上进行采集的办法。
20 世纪50 年代末期，伴随多次覆盖技术的出现和数据的可重复处理，地震勘探出现了革命性的突破，同时，采集、记录装置的更新以及非炸药震源尤其是电火花震源的出现，使得海洋地震采集技术得到了极大地发展，开始使用漂浮组合电缆在水下进行接收。
极大地提高了勘探的效率和精度。
这期间所使用的地震仪都是模拟地震仪。
20 世纪60 年代中期，伴随电子计算机以及计算机技术的发展，直接促成了20世纪70 年代数字地震仪的出现，伴随地震道数逐步由24 道发展到96 道，对震源能量和激发效率的要求也逐步提高。
20 世纪80 年代以来，海洋地震逐步朝着高分辨率、大激发能量、大接收道数的方向发展。
3．知识点：挑战者号科学考察是1872年至1876年期间，使用英国舰队挑战者号实施的一次科学考察活动，它完成了多项发现，为海洋学的建立奠定了基础。
挑战者号于1872年12月21日从英国的朴茨茅斯启程。
列出了超过4000种先前未知的物种。
经过1606天不间断的考察，713天的海上航行，挑战者号于1876年5月24日回到了汉普郡的斯彼特海德海峡。
海洋温跃层分析方法比较
海洋温跃层分析方法比较江伟;邢博;楼伟;连仁明【摘要】比较了垂直梯度法、曲率极值点法和拟阶梯函数法提取温度跃层信息的异同,结果表明:采用曲率极值点法和拟阶梯函数相结合的方法,能够给出较为准确的跃层上下界面位置,即跃层上界选用曲率极值点法和拟阶梯函数法确定跃层上界的最大值,而跃层下界则选用拟阶梯函数的结果.同时利用再分析资料初步诊断分析了西北太平洋冬季、春季、夏季、秋季温度跃层特征信息分布演变特征.【期刊名称】《海洋预报》【年(卷),期】2016(033)003【总页数】9页(P41-49)【关键词】温跃层;垂直梯度法;曲率极值点法;拟阶梯函数法【作者】江伟;邢博;楼伟;连仁明【作者单位】海军海洋水文气象中心,北京100161;海军海洋水文气象中心,北京100161;海军海洋水文气象中心,北京100161;海军海洋水文气象中心,北京100161【正文语种】中文【中图分类】P731.24跃层是海洋中重要的物理现象，针对所研究的物理量不同，海洋中的跃层可分为温度跃层、密度跃层、盐跃层、声跃层等。
海洋跃层的空间分布和季节变化与水团垂直边界的划定息息相关。
声信号是海洋中重要的通信媒介，声速的铅直分布特征对于水中通讯、水中目标探测具有重要的意义。
而海洋密度场结构直接决定着声速剖面，密度跃层是海洋密度结构的重要且典型的分布特征。
比较海洋跃层诊断分析方法，研究跃层的深度、厚度和强度及其时空演变特征，对于深入研究海洋跃层的形成和演变机理具有重要的科学意义，同时科学合理诊断分析跃层结构特征，有利于水中通讯、水下目标探测活动的开展，对海洋渔业、海上军事活动具有实际应用价值。
国外研究者早在20世纪60年代就开展了海洋温度结构方面的研究工作，比如Turner等［1］利用室内实验与理论分析相结合的方法研究了季节温跃层的形成和维持，Gill等［2］利用实测资料对季节跃层模型进行了检验，着重分析了动力混合和对流混合对上层温度结构的影响。
海洋预报技术及其应用研究
海洋预报技术及其应用研究海洋预报，是指根据海洋环境特征值的历史资料和实时观测结果，运用专门设计的物理模型和数学模型，对一定海域未来时间内的海洋要素、海洋现象、海洋变异及其可能造成的影响，以一定的文字、图表、声像等形式进行描述和发布。
海洋预报包括风暴潮、海浪、海啸、海冰、海流、海温、盐度、潮汐、海平面变化、厄尔尼诺、水质、海岸侵蚀等。
海洋积聚着巨大能量，这些能量一旦以某种突发的方式释放，则可能对人类构成非常大的危害，包括风暴潮、巨大海浪、严重海冰、海啸、赤潮等突发性较强的灾害；海平面变化、海岸侵蚀以及沿岸土地盐渍化等缓发性海洋灾害。
我国沿海地区人口稠密、经济发达，海上各类生产活动蓬勃发展，但频受海洋灾害袭击，每年因海洋灾害造成的直接经济损失在 100亿元以上。
风暴潮、海浪、赤潮、海水入侵和海岸侵蚀等较为严重。
进行海洋预报可以有效地对海洋灾害进行预警，有效经济损失和人员伤亡，在经济社会发展当中发挥了重要的作用。
目前人类社会对海洋分析预报的需求包括：海洋业务化预报、台风预报、季节和气候预测、全球变化研究、海上航运、渔业、目标漂浮物的跟踪预报、海上搜救、溢油预报、海上执法和管理、海上石油和天然气作业、海洋旅游等。
军事方面，海军对海洋分析预报的需求有：水面舰艇的航行安全、作战训练、武器装备性能的发挥、维护国家海上权益、海外救援等。
海洋预报可以对国防安全、海洋经济发展、海洋防灾减灾、海上重大活动、海洋权益维护等起到很好的促进和强有力的支撑，很多海洋预报的产品对军事活动和国防的安全等也有重要的作用。
海洋预报是一切海上活动的基础。
随着海洋数值模拟技术、全球海洋观测系统和海洋数据同化方案技术的发展，以及社会和军事活动对海洋环境信息的迫切需求，驱动着海洋分析预报技术的不断发展。
海洋预报的具体内容有：研制深、浅水海浪数值预报模式，研制适合于海洋环境预报的有限区域海面风场数值预报模式和台风风场数值预报模式，实现中海及邻近海域海浪数值预报。
Nudging
潮 汐现 象 的影 响 。潮 汐学 的研 究还 对其 他方 面如 海洋 污染 物扩 散 、 泥 沙运 动 、 水 体交 换 、 风 暴 潮 等 问题 的研 究 有重 要作 用 。因此 ， 针 对潮 汐进 行研 究具 有重 要 的理论 意义 和实 际应 用价 值 。
摘 要 ： 基于 Ｐ ＯＭ 模 式 ， 以西 北 太 平 洋 海 域 为 试 验 海 区， 针 对 Ｎｕ ｄ ｇ ｉ ｎ ｇ松 弛 项 的 差 分 方 案 及 Ｎ ｕ ｄ ｇ ｉ ｎ ｇ松 弛 系 数 进 行 了数 值 试 验 。结 果 表 明 ， 通 过 引入 Ｎｕ ｄ ｇ ｉ ｎ ｇ松 弛 项 的潮 波 模 拟 结 果 明 显 优 于 未 加 Ｎｕ ｄ ｇ ｉ ｎ ｇ松 弛 项 的模 拟 结 果 ；
中 图分 类 号 ： Ｐ ７ ３ １ ． ２ ３ 文献标识码 ： Ａ 文章编号 ： １ ６ ７ １ — ６ ６ ４ ７ （ ２ ０ １ ５ ） ０ ４ ０ ４ ３ １ － ０ ８
潮 汐是 海 洋最基 本 的运 动形式 之一 。它与人 类 活动关 系 密切， 准确 掌 握潮 汐 的运 动规 律 对 渔业 生 产 和
同化试 验 ， 优 化 了开边 界条件 ， 提 高 了数值 模拟 的精 度 。 吕咸 青和 方 国洪［ ２ 州 采 用伴 随法 ， 由渤 海沿 岸 １ ９个
验 潮 站 的潮汐 调和 常数 来反 演渤 海海 域 的开边 界条 件 ， 数 值模 拟 结果 较 好地 体 现 了渤 海 Ｍ 和 Ｍ。 潮波 的基 本 特 征 。吴 自库等口 阳利用 伴 随 同化 方法 ， 把 Ｔ／ Ｐ卫 星 高度 计 资 料 提 取 得 到 的沿 轨 分 潮 调 和 常数 同化 到 二 维非 线性 潮汐 数值 模式 中去 ， 通 过优 化模 型 中 的开边 界 条 件 和底 摩 擦 系 数 ， 较好 地 模 拟 了南 海 Ｍ 和 Ｍ。 分
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
作 者 简 介 ： 学 明 （９ ４ ， ， 北 省 卢 龙 县 人 ， 理 研 究 员 ， 士 ， 事 海 洋 环 流 和 海 洋 环 境 的 数 值 研 究 与 预 报 。Ｅ ｍａｌｚｕ ｍ＠ ｎ ｆ． 朱 １ ８ 一） 男 河 助 硕 从 — ｉ ｈｘ ： ｍｅｃ
近岸 浅海 区卫 星 高度 计资 料 的精度 并不 十 分令人 满
意， 由此计 算 的潮 汐调 和常 数也 不够 准 确 。
随着 计 算 机 和 电子 通 信 技 术 的 迅 速发 展 ， 值 数
模 式 在海 洋潮 汐 的模拟 和 预报 研究 中得 到 了广泛 应
用Ｕ ９， －］ 因而对 潮 汐 的分布 ， 别是 未 调 查海 区潮 汐 特
基金项 目： 国家 自然 科 学 基 金 （０ ７ ０ １ ； ４ ９ ６ ０ ） 国家 重 点 基 础 研 究 计 划 （ ９ ３ 计 划 ） 目（ ０ １ Ｂ ０ ６ ０ ； 洋 公 益 性 项 目（ ０ ８ ５ ８ ） 国 家 ８ ３ “ ７” 项 ２１Ｃ ４３０）海 ２０００７ ； ６
计 划 （０ ７ ２ ０ ＡＡ０ ２ ０ ） ９ １４ 。
（．国家 海 洋 环 境 预 报 中心 国 家 海 洋 局 海 洋 灾 害预 报 技 术研 究 重 点 实 验 室 ， 京 １０ ８ ；．中 国海 洋 大 学 海 洋 环 １ 北 ００１２
境 学院 ， 山东 青 岛 ２ ６０ ；．中海 石 油 研究 中 心 ， 京 １ ０ ２ ；．中 海 油 能 源 发 展 股 份 有 限 公 司 北 京 分 公 司 ， ６１０３ 北 ００７４ 北 京 １０ ２ ；．国家 海 洋 局 第 三 海 洋 研 究 所 ， 建 厦 门 ３ １ ０ ） ００７５ 福 ６０ ５
关键 词 ： 西北 太平 洋 ； 日本 海 ； 汐 ； 近 法 ； ＶＣ 潮 趋 Ｆ ＯＭ
中 图 分 类 号 ： ７ １ ２ Ｐ ３．３ 文献标志码 ： Ａ 文 章 编 号 ： ２ ３４ ９ （ ０ ２ ０ — ０ １ １ ０ ５ — １ ３ ２ １ ） ２０ ０ — ２
ｌ 引 言
潮 汐 是 中 国大 陆架 最 主 要 的 水 动力 要 素 之 一 ， 对 陆架 海 的海 洋环境 起 着重 要 的控制 作用 。它也 是 塑 造海 岸 带地 貌 ， 运海 底 沉积 物 ， 成 大范 围海 岸 搬 造 侵 蚀 和淤 积 的主要 源动 力之 一 。渤海 、 海 、 黄 东海 的 潮 汐分 布 ， 早是 由 Ｏｇ ｒ［ 在 ２ 最 ｕａ１ ０世 纪 ３ ０年 代 根
摘要： 利用 Ｆ Ｃ Ｍ 海 洋数值模 式， 球坐标 系统 下考虑 非 线性效 应 和天 体 引潮力 的 影响 ， 于 非结 Ｖ Ｏ 在 基 构 的三角形 网格建 立 了包括 中 国近海 、 本 海 、 日 鄂霍 次科海和 部分 西北太平 洋海域 的高 分辨 率海 洋潮 汐数值模 型 ， 并采 用趋近 法 同化 ８ ４个沿岸验 潮站 的观 测资料 。模 拟 结果 与 １ ５个 验潮 站 的实测结 果 ７
第 ３ ４卷
第 ２期
海
洋
学
报
Ｖｏ． １ ３４。 ＮＯ． ２
２１ ０ ２年 ３月
ＡＣＴＡ ＯＣＥＡＮＯＩＯＧＩ ＣＡ ＮＩ ＳＩ ＣＡ
Ｍ ａｒ ｈ ２ ２ ｃ ０１
西 北 太 平 洋 的 一 种 潮 汐 数 值 同化 模 型
朱 学 明 鲍 献 文 黄 必桂。 石 新 刚 ， 春 生 ， ， ， 靖
常数 ， 研究 中国近海 的：Ｏ ＯＯ 一３ ２ １ —４２ ； ２ １ ８ｌ 。
达 到较高 的模 拟精 度 口一 ］ ¨。
中国 近海 的潮 振 动 主要 由两 部 分 组 成 ： 北 太 西
。然 而 ， 在
平洋 潮波 向我 国沿岸 传播 引起 的谐 振潮 及 天体 引 潮
拟合 良好 ， ，２Ｋ Ｏ 四个 主要分潮振 幅和 迟角 的绝对平 均 误差 分别 为 ４ ０ｃ 和 ５ ６ ，． ｍ 和 Ｍ２ Ｓ ， ， １ ． ｍ ． 。２ ４ｃ
７ ５，． ｉ ． 。２ ６ｃ ｎ和 ６ ３，． ｍ 和 ５ ０。依据调 和分 析 结果 给 出 了 ４个 主要分 潮 的 同潮 图分布 ， 到 ８ ． 。１ ５ｃ ．。 得 个半 日分潮和 ５个全 日分潮 的无潮点 ， 实 了宗谷 海峡全 日潮无潮 点的存在 ， 证 首次模 拟 得到 津轻 海峡 的全 日潮无 潮点 ； 还给 出了整个计 算海域 内最大可 能潮差和 潮汐余水 位的分布特 征。
潮波 的伴 随模 式来 反 演 开边 界 和 底 摩 擦 系 数 ， 以求
据 实测 资 料作 出 的 。几 十年 来 我 国学者 分别 采用 不
同 的方 法对 中国陆 架浅 海 的潮 波系统 进 行 了大量 的 观测 和分 析研究 ， 得 了多 项重 要成 果 ， 取 使得 对 我 国 近海 的潮 汐 分 布 已有 了 比 较 好 的认 识 。Ｆ ｎ ｌ 根 ａ ｇ２ ＿ 据 大 陆或 岛屿 沿岸 观测 站 的实 测水位 资 料给 出 了 中 国近 海 Ｍ。 Ｋ 分 潮 的 同潮 图 。近 十 几 年来 ， 着 和 随 卫 星 高度 计 资 料 的 积 累 ， 们 又 采 用 各 种 方 法 从 人 ＴＯＰ ｘ／ ｏ ｅ ｏ Ｅ Ｐ ｓｄ ｎ卫 星高度 计 资 料 中提 取 潮 汐 调 和
分布 的认 识有 了很 大 提高 。但 前人 大都 独 立研究 中
国近 海 的某一 海域 或 海湾口 或 者 将 渤 海 和 黄海 作 ， 为整 体研 究 ¨， 将 渤 海 、 海 、 海 作 为 整 体 研 或 黄 东 究 ＿一 ］ 】 ” 。这 样一 来 模 式 的开 边 界 大 部 分 在 ２ ０１ 。 ０ Ｉ Ｔ 以浅 的海 域 ， 由全 球 潮 汐 模 型 提供 的 潮 位 开边 界 条 件在 这部 分海 区并 不 准 确 。为 了解 决 开 边 界 问 题 ， 人们 又采 用数 据 同化 技术 和伴 随法 ， 立 中 国 近海 建