opexhowidonp卫星高度计数据的潮汐分析.docxVIP

opexhowidonp卫星高度计数据的潮汐分析 1 国外海洋潮普组合测高模型的研究 海洋卫星研究技术的发展为海洋大陆地区的潮汐观测提供了新的方法。卫星探测高数据与沿海和岛货潮站数据的结合极大地提高了人们对全球海上潮汐分布的认识。随着卫星探测高距离定轨技术的发展，测量的精度不断提高，相应的全球海上潮汐模型也不断改进。1992年8月，主场景卫星t.p发射。测量表面高度的精度为5厘米，明显优于早期的海拔高度卫星。自t.p卫星发射以来，国外使用t.p测量数据进行了大量的研究，并发表了基于t.p卫星高度计算数据建立的潮汐模型文章。有些模式是基于和谐分析，而其他模式是基于反应法。 2 交叉点位置的确定 本研究采用法国空间研究中心(CNES)提供的综合地球物理数据记录(MGDR,Merged Geophysical Data Records). 从T/P卫星的CYCLE 11(1993年1月)到CYCLE 249(1999年6月),共有6 a多的资料. CNES的MGDR资料也同时提供了每个CYCLE的交叉点资料. 为进行调和分析,交叉点资料与沿轨资料需先经过一系列改正,并形成时间序列. 卫星高度计资料包含大量的大气和海洋信息,例如对流层、电离层、大地水准面、潮汐、波浪和海流等. 为有效确定潮汐,必须对卫星测高数据进行一系列改正,滤除与海潮无关的信息. 本文在交叉点处进行下列改正: ζ=Sat—Alt-Η—Alt-Dry—Corr-Wet—Corr-Ιono—Corr-Ιnv—Bar-SSB—Corr-Η—ΜSS-Η—Set-Η—Ρol-Η—Lt?(1)ζ=Sat—Alt?H—Alt?Dry—Corr?Wet—Corr?Iono—Corr?Inv—Bar?SSB—Corr?H—MSS?H—Set?H—Pol?H—Lt?(1) 其中,ζ是改正后的海面高度;Sat—Alt是相对于参考椭球的卫星高度;H—Alt是卫星测高值;Dry—Corr和Wet—Corr分别是干、湿大气订正;Iono—Corr是电离层订正;Inv—Bar是反变气压计订正;SSB—Corr是海况订正;H—MSS是相对于参考椭球的平均海面高度;H—Set和H—Pol分别是固体潮和极潮;H—Lt是负载潮. 上述各项改正值皆来自于CNES的MGDR数据光盘. 对交叉点数据,先经过上述一系列改正,然后形成本文交叉点处海面高度时间序列. 沿轨资料的处理与交叉点资料的处理类似,但其形成时间序列过程要复杂一些. 因为对于某一CYCLE的观测点,其下一个CYCLE并非都是精确重复的. T/P卫星地面轨道的摆动范围约为±1 km;沿轨方向的最大偏差可达一个测距,即约6.2 km. 因此进行共线处理是必要的:选一个参考CYCLE,将其他CYCLE的测值都处理到参考CYCLE的测点上,从而形成时间序列. 显然,此参考CYCLE须满足以下两个条件:一是观测点尽量多,二是其地面轨迹处于摆动范围中间. 本文取CYCLE 91作为参考CYCLE,经共线处理后形成时间序列. 关于沿轨资料的改正及共线内插形成时间序列的方法和过程可参考董晓军等的文章. 3 潮化过程和调和过程的求解 T/P重复取样周期是9.915 6 d,本文所用T/P资料取样长度为6 a多. 因此存在两类频率成分混淆问题:一类是由于采样间隔引起的高频折叠,即“aliasing”问题,对此李立等已作了详细的讨论;另一类,由于取样长度不够长,即短于某两个分潮的会合周期,也将导致这两个分潮的混淆,本文将主要讨论这一类混淆问题. Schrama和Ray指出:严格应用Rayleigh判据进行判断,分离M2和S2需要将近3 a资料,分离K2和P1需要9 a资料,分离K1和SSA也需要至少9 a的资料. 而本文所用T/P沿轨资料序列的长度不足9 a,直接应用调和分析方法显然难以分离混淆较严重的分潮,如K2和P1,K1和SSA. 我们注意到,交叉点资料包括升轨资料和降轨资料,其资料量比沿轨点资料量多1倍. 据此先应用调和分析方法对交叉点资料进行分析. 然后将交叉点处分潮之间的关系引入到沿轨的相邻点处,对沿轨点的时间序列进行调和分析,实现对某些混淆分潮的分离. 潮汐的调和分析方法是根据最小二乘法将观测的海面高度拟合到主要的分潮,从而得出各分潮的调和常数. 把交叉点或沿轨点的实测水位ζ作为多个调和分潮的叠加,即: ζ=Ζ0+Ν∑n=1fnΗncos(ωnt+Vn-gn)?(2)ζ=Z0+∑n=1NfnHncos(ωnt+Vn?gn)?(2) 其中,Z0为平均海面高度;ωn为第n个分潮的角速率;Hn,gn分别为第n个分潮的振幅和迟角;fn和Vn为相应的交点因子和天文初位相;N为分潮个数. 在每一个观测点,要求其观测次数大于潮汐调和常数的个数. 这样形成超定方程组