CORS站高程非线性速度场及方差波动模型构建方法

本文描述了CORS站高程时间序列采用固定周期项建立线性速度场的基本概念,指出给定周期与实际周期存在偏差的问题。为此,提出了CORS站高程时间序列非线性速度场建模方法。该方法以线性最小二乘解作为迭代初值,利用高斯-牛顿迭代算法求解非线性速度场模型的未知参数,实现了CORS站高程时间序列数据的非线性建模;给出了拟合模型残差平方序列异方差特性的检验方法,阐述了建立GARCH(p,q)模型反映CORS站高程非平稳序列波动情况的基本准则。本文以国内外6个CORS站20余年的高程时间序列为研究对象,建立非线性速度场运动模型。结果表明:CORS站高程运动并不存在严格的整年或半年周期项,近似年周期运动最为明显,近似两年周期运动占比最小,年周期估计偏差达到12%,半年周期偏差18%,两年周期偏差6%,非线性建模精度和效果整体优于线性模型。利用ARCH检验法得到CORS站高程非线性模型的残差平方序列存在异方差性,即残差平方序列具有非平稳特性,引入GARCH(p,q)模型对CORS站高程分量非平稳残差平方序列建模,反映了残差平方序列的非平稳波动。验证了GARCH(p,q)模型对CORS站高程非平稳残差平方序列建模的可行性,为今后利用GARCH(p,q)模型对CORS站高程非平稳噪声序列建模和非线性速度场重构提供了思路。     

基于台州及周边地区CORS站网地表形变时空特征分析

基于台州地区15个连续运行参考站（CORS）2017-2019年的时间序列,对其进行线性拟合,构造三维速度场,并利用功率谱分析及小波变换分析其CORS站周期特征．发现台州地区CORS中基准站水平方向具有平均站速度33.58mm／a,运动方向为东南的水平运动现象,高程方向最大抬升速率3 mm/a,最大沉降速率1 mm/a;从周期上看,CORS站垂直方向相较于水平方向明显很多,具有明显的年周期、半年周期以及季节周期．      

区域CORS网坐标时间序列周期特征分析

通过精密单点定位的方法获取了区域CORS网各基准站的坐标时间序列,并利用最大似然估计法和频谱分析法对其32个基准站坐标时间序列的周期特征进行了分析,获取了各基准站年周期、半年周期的振幅和相位值,以及各基准站的周期功率谱图。结果表明:GNSS基准站坐标时间序列不仅存在线性变化,还存在周期性变化,其中以U方向表征最为明显,与双差定位获得的坐标时间序列周期特性分析结果一致,说明以精密单点定位获取的GNSS基准站坐标时间序列是可靠的,可以用来分析基准站的变化特征。     

利用GNSS数据分析大港验潮站地壳沉降

青岛大港验潮站的地壳沉降关系到该站平均海平面的绝对变化,因而也就关系到我国高程基准面的变化。本文利用青岛GNSS基准站约10年的观测数据对该站的地壳沉降变化进行分析。首先将青岛GNSS基准站纳入由50个国际IGS站和43个国内陆态网络基准站组成的全球网中,进行单日松弛解和单日约束解解算,获得该站坐标时间序列。然后对该站垂向坐标时间序列进行分析,利用粗差探测、偏差探测、趋势项分析、频谱分析等方法对粗差、偏差、趋势项和周期项进行探测、分析,并通过时间序列模型估计获得时间序列中的周期项振幅和偏差估值。分析表明青岛GNSS基准站垂直方向近一段时间未发现存在显著性的地壳沉降变化,但受到比较明显的周年和半周年周期变化影响。结合青岛大港验潮站验潮数据分析结果得出结论:青岛大港验潮站平均海平面的绝对上升速率是1.62mm/a。     

地心运动时间序列的抗差谱分析

本文将抗差估计与谱分析相结合，提出了抗差谱分析。用该方法对全球约120个IGS站5年GPS观测结果求得的地心运动时间序列进行了分析，发现地心运动在各方向有不同的主周期变化。X方向受地球物理交叉影响较大，无法分离出对该方向影响较大的主周期变化；而Y，Z方向都有明显的主周期变化。Y方向有明显的年周期变化、半年周期变化、季节性变化和月周期变化，其中，年周期变化的贡献最大；Z方向的主周期是年周期变化、半年周期变化、51天和609天的周期变化。与X，Y方向相比，Z方向振幅明显偏大。     

南极IGS基准站时间序列极大似然法分析

针对南极大陆IGS的周期特性、噪声特性、趋势特性等问题,该文采用极大似然法和功率谱分析法分析南极地区8个IGS永久观测站坐标的时间序列。研究结果表明,该地区水平和垂直方向上均有较明显的年周期项,垂直方向上还普遍存在半年周期项。闪烁噪声与白噪声组合能更好描述南极基准站3个分量上噪声特性,各方向上的速率不确定度比只考虑白噪声时要大10~15倍。通过估算基准站水平方向的运动速率,发现该方向上具有较规律的运动趋势:南极洲正缓慢向南美洲移动。     

GNSS高程时间序列周期项的经验模态分解提取

针对传统方法——常振幅常相位的周年+半周年谐波模型法不能准确提取GNSS高程时间序列中周期项问题,该文以中国区域10个IGS基准站在ITRF2008框架下2005—2015年高程时间序列为例,采用经验模态分解(EMD)提取各测站高程时间序列的周期项。对传统方法和EMD两种方法提取的周期项做Lomb_Scargle谱分析,用功率谱图分析了这两种方法提取序列周期项的能力。实验结果表明,EMD方法较传统方法更能准确、自适应地提取GNSS高程时间序列的周期项。     

CME对IGS基准站坐标序列噪声模型及速度估计影响分析

以50个IGS基准站坐标时间序列为研究对象,采用赤池信息量准则（BIC）模型估计准则对四种组合噪声特性进行估计分析,探讨共模噪声（CME）对IGS基准站坐标时间序列噪声模型及站速度影响. 结果表明CME会导致IGS基准站坐标序列噪声模型的有偏估计,并影响站速度的确定精度,准确估计基准站速度参数时应对CME进行修正;经CME修正后IGS基准站周年运动周年项振幅有所减小,高程方向更为明显,表明滤波后IGS基准站周年运动相对稳定.      

华北区域GNSS基准站高程时间序列的非线性变化分析

基于陆态网络全球卫星导航系统（GNSS）观测成果,采用功率谱分析法和最小二乘方法,以华北地区为例,研究了区域基准站高程时间序列的非线性变化特征,并分析了不同环境负载的影响.结果表明,GNSS基准站高程方向存在显著的周年和半年周期特征,且周年特征要显著于半周年特征.位于不同地区的基准站的振幅和相位存在差异,华北平原南部地区的周年振幅要大于北部地区,整体上华北地区周年变化在秋季时节振幅达到最大.不同环境负载效应对华北GNSS高程位移的影响不一致,利用三种环境负载修正GNSS序列后,水文负载的修正效果最好,非潮汐大气负载次之,非潮汐海洋负载修正结果不理想.      

利用CEEMD分析中国沿海GNSS站高程时序变化

利用完备经验模态分解方法（CEEMD）对我国沿海地区6个GNSS基准站（2010—2018）的高程时序数据进行了处理分析。结果表明：CEEMD在高程时间序列分析中具有一定的优越性,可准确分解出各GNSS站高程时序中存在的周、月、季节、年等变化周期项,其中周年运动是主要贡献项,各站高程时间序列的短周期变化与潮汐变化周期具有密切关联性;沿海GNSS站的地面沉降既具有区域的一致性,又存在区域间差异性,其中D区DBJO、DZJJ站呈现先下降后上升的趋势,N区NZUH、NWZU站呈下降趋势,B区的BZMW呈上升趋势,而同海区的BLHT站则呈显著的下降趋势。     

GPS参考站垂直位移时间序列的主成分分析

利用了主成分分析方法,分析了中国地壳运动观测网络主要基准站垂直位移时间序列的频谱特征,得出了站点的垂直位移变化具有很强的周年变化特征。     

香港CORS站坐标时间序列分析研究

针对不同地区连续运行参考站(CORS)坐标时间序列周期特性与噪声特性存在的差异性问题,本文采用功率谱分析法和极大似然法分析香港卫星参考站网(SatRef)的坐标时间序列.研究结果表明,参考站点的三个坐标分量上都存在明显的年周期项与半年周期项.各参考站坐标分量上的噪声特性存在多样化的特征,白噪声(WN)加闪烁噪声(FN)是主要噪声模型,仅考虑WN而忽视有色噪声(CN)的做法会大大低估参数估计的不确定度.根据最优噪声模型估值得出的SatRef站速度场在水平方向上有整体向东南方向运动的趋势,与华南块体的运动结果基本一致.     

汶川地震前后四川盆地CORS站运动特性分析

为探讨四川盆地在汶川地震前后地壳运动的规律,本文对2006年3月至2012年9月四川12个连续运行参考站（Continuously Operating Reference Stations, CORS）6年多的观测资料进行了详细的处理和分析。通过比较震前震后CORS站运动的相似性和差异性发现：由于震后形变的影响,震后第一主成分(Principle Component First, PC1)的贡献率下降20%~40%;震后多站的速度场发生了明显的变化,突出表现在PIXI,CHDU,MYAN等站,YAAN,QLAI站处于闭锁状态;盆地垂直分量的周年振幅在震前一年偏大,震后逐年减小。以上发现表明汶川地震可能改变了四川盆地原有的运动特性。     

北斗卫星伪距偏差模型估计及其对精密定位的影响

现阶段北斗卫星导航系统（BeiDou navigation satellite system,BDS）的同步地球轨道（geostationary orbits,GEO）卫星、中倾斜地球同步轨道（inclined geo-synchronous orbits,IGSO）卫星和中圆地球轨道（medium earth orbit,MEO）卫星均存在伪距偏差,该伪距偏差的存在对精密定位的研究及其应用产生了较大的影响。根据北斗IGSO和MEO卫星的伪距偏差与高度角和频率相关的误差特性,本文分析了测站数目及分布,以及观测时长对建模的影响,选择18个测站2015年全年的数据作为MEO卫星的建模数据,其中可以连续观测到全弧段IGSO卫星的4个测站用于IGSO卫星的建模,采用加权分段线性拟合联合抗差估计的方法建立了北斗卫星伪距偏差改正模型。模型改正后,北斗IGSO和MEO卫星的伪距偏差得到明显的削弱,相比于传统的伪距偏差改正模型,精密单点定位（precise point positioning,PPP）的定位精度和收敛时间均得到提升。     

我国IGS站的周期变化特征

利用谱方法和小波技术分析了我国IGS站坐标时间序列,阐述了其周期变化特征,即绝大部分台站在N、E、U方向都具有半年周期性、年周期性,U方向周期特性相对明显,同时台站的主要周期成分和时频局部特性又具有一定区域性.     

利用GPS和GRACE研究澳大利亚地壳垂向季节性变化

为探究重力场恢复与气候实验(gravity recovery and climate experiment,GRACE)卫星与全球定位系统(global positioning system,GPS)两种独立技术获取的因陆地水储量变化引起的地壳垂向季节性位移的一致性,选取澳大利亚27个GPS站点5～10 a的高程时间序...