基于迭代PCA的GPS时间序列震后形变估计方法和应用

摘要

GPS时间序列的震后形变分析对于研究区域震后形变机制和岩石圈流变学特性以及维持国际动态地球参考框架具有重要意义.本文在现有参数估计方法的基础上,提出兼顾震后形变衰减特征空间相关性和整体建模的“迭代PCA参数估计方法”,并利用模拟数据证实了新方法可以获取更稳健可靠的震后形变、同震形变和震间速度参数.最后,以37个新西兰GPS连续站坐标时间序列为例,利用迭代PCA方法提取了2016年11月13日Kaikoura地震共性的震后形变时间演化过程和各站点的震后形变,并定量分析了震后形变对地表速度的影响.结果表明各站的震后形变在时间域上以衰减常数τ为4天的对数模型持续松弛;空间域上南岛北东部和北岛最南部震后形变较大,其中,最大震后形变点为cmbl站,截至2017年6月10日NEU方向累计的震后形变分别达到107mm,135mm和187mm,地表速度分别达到133.58 mm·a-1,112.05 mm·a-1和175.58 mm·a-1,仍高于稳定的震间速度.