第二章 地球磁场.doc

第二章 地球磁场 （Lisa Tauxe著，常燎译） 建议补充读物 Butler (1992)，3－7页，10－11页。 更多信息可参看：Merrill et al. (1996) 第一、二章。 2.1 地球磁场 2.1.1 地球参考场 这个方程可以改写为： 这个方程的一个解是： 对地球磁场，一般可以写作半径为r，纬度余弦 其中，g和h是高斯系数，可以从特定的年代计算得出，单位为nT，或磁通量（注意，公式中(0由tesla [B]转换到Am-1 [H]）。角标e和i代表外场或者内场的起源，a是地球半径（6.371 х 106 m），(0是自由空间的磁导率（参看第一讲中的表1.1），正比于勒让德多项式，其由传统的施密特多项式归一化而来（可参看建议的读物）。 图2.2展示了三种矢量场的全球倾角分布及对应的面谐函数的：即轴向的（m＝0）偶极子场（l＝1），四极子场（l＝2），以及八极子场（l＝3）。它们的贡献分别由，和确定。相关的多项式（图2.1）为： 如果转动图2.2a中的轴向偶极子场，使其北极指向格林威治子午线，那么它就由系数确定，如果指向90(E，则将由系数而定。所以，总的偶极子贡献将是轴向和两个沿赤道的偶极子项的矢量相加，即。总的四极子贡献（l＝2）由五个系数而定，总的八极子（l＝3）贡献则由七个系数而定。 一般来讲，如果下标（l）和上标（m）的差为奇数（比如，轴向偶极子和八极子），则相应的地球磁场对于赤道是非对称的。然而，如果l和m的差为偶数（如，轴向四极子），则相应的地球磁场是对称的。图2.2a表示，由与现今地磁场方向一致的单一偶极子场产生的倾角。在北半球，倾角都是正的（向下），而在南半球是负的（向上）。相反，由四极子场产生的倾角（图2.2b）是在极区是向下的，在赤道处则是向上的。由轴向八极子场（图2.2c）产生的倾角关于赤道也是非对称的，在两个极区的方向相反，并在中纬度地区被具有相反方向的条带分开。 地球磁场是一个矢量场，所以在每个点都有方向和强度（图2.3）。无论选择怎样的坐标系，三维矢量场都需要三个参数来定义。比如在笛卡儿坐标系下，用x, y, z或x1, x2, x3。对于特定的问题，由于问题本身的对称性，某种坐标系会更合适些。除了笛卡儿坐标系外，也应用其他一些坐标系，但需要在这些坐标系间进行转换。 图2.1：施密特多项式。 水平分量也可以投影到向北（X1）和向东（X2）的轴（一般是测量的方向），即： 公式2.2和2.3对于分解磁场各分量都很有效。 图2.2：全球磁倾角的表面谐函数（附图）及其相应的分布图（图片来源于Tauxe，2005）。a）偶极子；b）四极子；c）八极子。 注意正切函数符号的复杂性。你也许会弄错象限，最后就不得不加180(。 回想第一章（包括附录），一旦标量势已知，地磁场的各分量就可以由计算出，所以，在球坐标系下： 其中，r、(、(分别是半径、余纬（偏离北极的角度）和经度。这里，正向的BV向下，BN向北，和第一章中的Hr和H(相反。注意公式2.1的单位是特斯拉，而不是Am-1。 另外，如果已知磁场矢量场，也可以推导出势场。对于特定参考场，其高斯系数一般是对观测到的地球磁场进行最小二乘法拟合后得出的。为了可靠地估算高斯系数，如果到L＝6则需要至少48个观测点。 图2.3：地磁矢量场B的各分量。BH是矢量场B在地表的切向投影。BH可以分解成向北和向东的分量（BN和BE）。BV是垂直轴向的投影。D是从北向开始0到360度按顺时针方向增加。I是从水平方向开始从－90度到＋90向下增加（因为磁力线也可以指地球的外部）。如果需要，M或H也可以被B代替。 . 为了了解不同高斯系数的作用，表2.1总结了Olsen et al. (2000)估算的前六阶高斯系数。每阶的能量为（Lowes, 1974）(图2.4)。最低阶项（一阶）占主导，几乎占据90％的地磁场。这也正是为什么地球磁场通常被认为是等价于一个处于地心的简单偶极子场。 现在来看看根据1995年的IGRF估算的地磁场参数。对于给定参考场，应用公式2.1和2.5我们可以计算地球任何地方的B，D和I值。图2.5显示地表的磁场是位置的函数。在极区的磁场强度一般可高达~60 (T，而在赤道处最低（~30 (T）。但是与地心轴向偶极子（GAD）产生的磁场不同（图2.6），地球磁场强度的等值线图和纬度不平行。而地心轴向偶极场产生的倾角也规则变化，在极区为－90度和＋90度，在赤道处是0度；等值线图平行于纬度线。虽然图2.5b显示的倾角图类似于GAD模型场，但是还有很多不同之处，这也暗示出地球磁场不能简单地由置于地心的条形磁铁来代表。如果地球磁场是简单的地心轴向偶极子场（GAD场），那么在无论什么地方偏角都为零，实际上显然不是这样的（图2.5c