基于主成分分析方法的大地测量时间序列反演.PDFVIP

翻译者：贺小星，hexiaoxing@ 基于主成分分析方法的大地测量时间序列反演 A. P. Kositsky[1,2]，J.-P. Avouac[1] 收稿日期：2009 年4 月13 日; 修订日期：2009 年8 月24 日;接受日期： 2009 年9 月16 日; 发表日期：2010 年3 月3 日。. [1]全球定位系统（GPS ）系统的发展，使的通过板块边缘分布的台站监测地球表面的形 变，并达到一定的精度成为可能。在滑深板块边界的时空演化，理论上而言，无论是地震形 变或震后滑走，通过时空分布在板块边缘的台站获得的观测数据，都可以采用一些基于弹性 半空间位错理论程序反演得到。笔者讨论和测试了基于主成分分析的反演方法 （PCAIM ），通过地表位移时间序列，采用主成分分析的繁衍方法。结果表明可以通过主 成分分析方法中的各个主分量反演出断层的滑变过程。因为PCAIM 不需要外部强加的时间 滤波，它可以处理任何形式的时间断层滑动的变化。我们把PCAIM 方法应用到大地测量时 间序列，通过此方法来反演一个复杂的滑变过程，包含同震形变，震后形变，非平稳的震间 滑移，以此来测试这种方法。 PCAIM 的滑移模型与标准的反演技术相比，计算更简便。同 时我们也对通过PCAIM 方法和ENIF 方法反演的2008 年8 月6 号尼亚斯地震发生后该地区 的震后形变作比较。我们介绍算法的几个扩展，使得能够应用于多个时间尺度的严谨统计学 中集成多个数据源（例如，GPS 和干涉合成孔径雷达的时间序列）。 PCAIM 可以推广到任 意线性反演算法。 引用：Kositsky ，A,P.,J.-P. Avouac(2010) ，基于主成分分析法的大地测量时间序列反演方法， 地球物理学研究，115， B03401 ，doi ：10.1029/2009JB006535。 1 引言 [2]断层滑动的方式很多，如突然发生地震事件或抗震蠕变的结果。断层滑动速率从时间尺 度上看差异很大，从典型的10 – 100 秒的大地震到几周或者数年的缓慢地震以及震后形变 等。因此监测断层滑动如何随时间变化，成为我们了解滑变过程的关键。深层地表滑变引起 的地表形变可以通过大地测量技术观测 [如Lisowski1991 年，Segall and Davis，1997]古老的 大地测量技术[如Taylor 1987，Sieh 1999]，或遥感技术[如Massonnet ，Feigl ，1998]。因此， 深层断层滑变可通过地表形变模型间接得到。 [3]理论上而言，通过深层的断层滑变来获取地表位移，一般都是基于线性弹性理论[如 Savage1983，Okada1985 年，Cohen，1999 年] 。这重理论是一种线性的、易通过标准的算法 反演。断层滑变的分布一般通过离散化的断层几何学来实现参数化。累积的断层滑动用来解 释发生在两个历元间的位移，可以通过最小二乘反演的大地测量数据来实现。一般而言参数 的数量超过观测值得数量，因此通常需要增加一些约束；例如，滑动分布粗糙度处带来的不 利影响可以增加一个限制约束。一种反演与时间相关的滑变分布的大地测量时间序列包含在 两个观测位移的连续历元中。这种方法是当计算非常密集的时代数量是大的，尤其是当非线 性正规化标准使用。此外，这种方法认为每个划时代表明，所以在不同的时间步长的测量误 差不正确的平衡。当历元数很多时，这种方法计算是很有效的，特别是采用非线性正规化标 化后。而且，该方法单独考虑每一个历元，因此不同历元之间的观测误差不均衡。此外，该 方法还要求每个站点在同一历元获取大地观测的时间序列，这限制了数据融合分析的可能 性，其中可能包括动态数据或干涉合成孔径雷达 （INSAR ）数据。 [4]P.Segall 和他的同事提出了一种改进的逐历元的反演方法，称作延展全网逆推滤波法 (ENIF)[Segall,Mattthews,1997;McGuire,Segall,2003 ]，尤其是在GPS 测量中的应用，ENIF 考 虑了本地基准的随机描述。把滑变速率作为一个非参数的时间函数，估计和去除参考框架误 差，进而使得使用扩展卡尔曼滤波方法去滤去时间相关的噪声。这种方法已应用于多项研 究，并取得巨大的成功[如Miyazaki2004] 。然而，这种方法具有一定的限制性。一是该模型 中包含众多的假设参数，使繁其成为一个繁琐的工具，需要做出一些主观的选择去定义可能 翻译者：贺小星，hexiaoxing@ 的解决方案。另一方面，该方法计算量大