第十二章 非Fisher统计.docVIP

第十二章 非Fisher统计 （Lisa Tauxe著，黄宝春译） 建议读物 背景知识：Borradaile (2003)第一章；Efron和Tibshirani (1982)；Tauxe等(1990) 详细了解：Fisher等(1987)第二、三、四、五章 12.1 前言 自上世纪五十年代以来，古地磁学研究人员对单位矢量数据的分析一直依赖于Fisher (1953)提出的特殊的统计方法。基于Fisher统计模型，各种功能强大且实用性很强的工具使得古地磁参数的定量化分析成为可能；如地块的旋转量、平均地球磁场是否为一不依赖于其极性的轴向地心偶极子场。然而，所有这些古地磁学工具，均要求古地磁数据服从Fisher分布。 事实上，在多种情况下，Fisher分布并不能充分地描述古地磁数据，或者地球磁场的方向本身很可能根本不服从Fisher分布（Creer等，1959）。对最近五个百万年以来的全球古地磁方向数据的分析（如：McElhinny和McFadden, 1997）表明，赤道地区的数据分布具显著拉长的特征（图12.1a），只有高纬度地区的数据对称性较好（图12.1c）。 方向数据之所以不符合Fisher统计模型，是因为地球磁场本质上为一个轴向地心偶极子场和其它高阶项产生的“噪声”的矢量和；这样的模型相对应的虚地磁极（VGP）是圆对称分布的。随着采样点靠近赤道地区，圆对称分布的VGP数据转换成方向数据后，其分布变得更加拉长（参见图12.1b, c）。在古地磁数据中，远离磁极的VGP对应于相对较弱的磁场强度；而在地磁场的许多模型中，由于Fisher分布假设了剩磁矢量为单位矢量，因而过分地强调了“外围”的贡献，导致平均倾角比真平均倾角要浅（如：Creer, 1983）。 不服从Fisher分布的另一个情况是地球磁场有两个稳定的极性状态。由于Fisher分布为单峰分布，因此，在对双极性古地磁数据进行平均之前，需要将其分解为两组单一极性进行计算，或者需要将其中一种极性“翻转”成其对称极性进行计算。此外，由多个组分构成的剩磁矢量易于形成条带分布。构造的复杂性（如：褶皱）也有可能造成方向数据呈条带分布。 可见，古地磁数据的确存在非Fisher分布的情况。由此，基于Fisher统计的各种不恰当的检验可能导致错误的解释。在第十一讲中，我们学习了Fisher统计的基本理论以及如何检验数据是否服从Fisher分布。在这一讲中，我们将讨论当古地磁数据不服从Fisher分布时，该如何进行统计与分析。我们将首先介绍几种主要的处理非Fisher分布数据的参数统计方法；然后再讨论如何在古地磁学中应用这些统计方法，如解靴带法12.2 古地磁剩磁矢量的非Fisher分析法 12.2.1 Kent分布 如图12.1a 所示，赤道地区最近五个百万年以来样品的剩磁方向（D′，I′）具有椭圆化的分布，因此并不服从Fisher分布（对称分布）。处理这样的数据，应用Fisher置信椭圆并不适当；而应用一种允许数据在椭圆长轴方向上离散的分布（即椭圆分布）则更贴切。与Fisher分布相对应的椭圆分布是Kent分布（Kent，1982）；其分布密度函数可 图12.1：a) PSVRL数据库内收集的赤道两侧0－5o纬度带内古地磁方向数据（参见McElhinny和McFadden (1997)）。倒转方向用其对趾点表示。所有方向均采用第二讲中介绍的 D′和I′转换方法旋转至以期望方向（五角星：D = 0o, I = 0o）为中心的坐标系中。上半（下半）图面的方向比期望方向倾角浅（陡）；右半（左半）图面的方向位于期望方向的右手（左手）方位。b) 由Tauxe和Kent（2004）统计场模型，在北纬30o（图中由正方形所表示的观测点）处获得的地磁场矢量换算得到的虚地磁极位置（VGP）。三角形代表地理极。黑色圆环代表了纬度为北纬60o的虚地磁极位置。c) 由图 b）中黑色圆环所代表的纬度为北纬60o的虚地磁极位置，转换得到的观测点（图b中正方形）处的剩磁方向。这些方向被投影至以观测点的期望方向（三角形）为中心的坐标系中。值得注意的是，一个以地理极为中心的圆对称圆环上的VGP数据转换成方向数据后，其分布是不对称的，且有显著的倾角浅化现象（引自Tauxe和Kent, 2004）。 表示为： F = c(κ, β)-1 exp (κ cos α + β sin2 α cos 2φ) 其中，α 为给定方向与真平均方向间的夹角（由方向矩阵 T 的最大特征值对应的特征向量 V1（参见第九讲附录C）来估计；φ 是给定方向在垂直于真平均方向的平面内的投影所确定的某一角度，当 φ ＝ 0 时，投影与该平面内的主要特征向量 V2 相平行。κ 是聚集参数，类似于Fisher统计中的精度参数；β 为“椭圆度”参数。c(