基于CM4模型的中国大陆地区地磁场时空分布特征分析1

李细顺 高登平 刘立申 赵志远 王利兵



基于CM4模型的中国大陆地区地磁场时空分布特征分析1

李细顺 高登平 刘立申 赵志远 王利兵

（河北省地震局红山基准台，河北邢台 055350）

本文利用第四代地磁场综合模型（Comprehensive Model 4，CM4），计算了1982—2001年中国大陆地区同一经度链和同一纬度链上地磁台站的磁层源磁场及其感应场、电离层源磁场及其感应场的地磁北向分量、东向分量、垂直分量的模型值，分析了各场源磁场随时间和空间的变化特征。结果表明：在时间上，经度链和纬度链台站的磁层源磁场及其感应场均呈现出11年和27天周期性变化。电离层源磁场及其感应场具有明显的季节变化，不同年份相同季节变化形态一致但幅度不同。在空间分布上，经度链和纬度链台站磁层源磁场及其感应场的年变化幅度呈现出不同变化特征，电离层源磁场及其感应场在经度链上变化特征不同，而纬度链台站的数值基本一致。日变化分析显示，磁静日和磁扰日期间，模型数据与台站实测数据变化一致性较好，相关性较高。

CM4磁层源磁场 电离层源磁场 实测数据 感应场 时均值 相关系数

引言

利用地磁场模型研究地磁场变化在空间物理领域一直是一个基本问题，国内外均有科学团队在关注这一问题。

地磁场是指在固体地球内部和外部到磁层顶空间范围内所有场源产生的磁场，一般将地磁场分为内源场和外源场，人们对于内源场研究较多（Haines，1985；Alldredge，1987；徐文耀，2002，2009；Wardinski等，2006；Lesur等，2008），对地壳磁异常场也做了很多工作（Hemant等，2005，2007）。地磁模型能对地磁场时空分布进行较为直观的表示。1968年国际地磁和高空大气物理学协会给出了第一代国际地磁参考场模型（IGRF），之后每5年给出一个新的IGRF（杨云存等，2014）。20世纪90年代，美国国家宇航局戈达德飞行中心（NASA/GSFC）和丹麦空间研究所（DSRI）联合开发了一种新的地磁场建模方法——“综合建模”（Comprehensive Modeling，简称CM）（白春华等，2008）。1993年，Sabaka等建立了第一代地磁场综合模型——CM1，1996年，Langel等在CM1的基础上建立了第二代地磁场综合模型——CM2（Langel等，1996），2002年，美国和丹麦空间研究中心的科学家Sabaka和Olsen等人建立了第三代地磁场综合模型——CM3（Sabaka等，2002），它相对于第一代和第二代的提高主要体现在进行了高阶地壳层场的静态表示，能在卫星高度获取大部分地壳异常场。CM3模型内源场最大截断阶数为65，时间跨度为1960—1985年，长期变化通过3次B-spline方法的13阶表示（冯彦等，2011）。

2004年，Sabaka等利用卫星和台站数据在CM3的基础上建立了第四代地磁场综合模型——CM4（Sabaka等，2004），相较于CM3，CM4采用了更多卫星数据，包括POGO和CHAMP卫星标量数据，MAGSAT和Frsted卫星标量和矢量数据。由于卫星位于电离层和磁层间的耦合电流区域（用于区分电离层和磁层的磁场强度），故可获取近地球区域各种场的电流源，并进行各种场的参数化和协同估算，通过迭代重加权最小二乘法进行合理划分，最后顺序计算7种磁场的值（地核场、地壳场、磁层源场、磁层源感应场、电离层源场、电离层源感应场、空间环形磁场）。CM4模型适用于磁静日期间，可对地磁场进行分层计算，相较CM3模型，其适用时间范围更长，为1960—2002年，进一步降低了预处理时的数据噪声。模型1—15阶代表地核场，16—65阶代表地壳场。

地磁台站记录到的地磁变化信号主要受到外源场和内源场两种场源的影响。其中，外源场起源于地表以上的空间电流体系，这些电流体系主要分布在电离层、磁层和行星际空间（李琪等，2015）。我国基于台站和卫星数据对电离层场以及磁层场开展了很多研究。丁鉴海等（2005）利用不同地区电离层台站数据对电离层2层的临界频率02进行了相关研究；姚丽等（2010）通过分析位于向阳面正午两侧的GOES-10和GOES-12卫星观测数据发现，激波作用于磁层时靠近清晨的磁场变化表现出简单压缩效应，而靠近黄昏的磁场变化则显然不同，即B分量减弱，B分量几乎为零，B分量则显著增强；徐文耀（2011）从能量守恒原理出发，讨论了太阳风-磁层-电离层耦合过程的能流路径和能量收支的定量关系问题。本文利用第四代地磁场综合模型（CM4），计算了1982—2001年中国大陆地区经度链和纬度链上地磁台站的磁层源磁场及其感应场、电离层源磁场及其感应场的地磁北向分量、东向分量和垂直分量，分析了各场源磁场随时间和空间的变化特征。对CM4模型计算出的7种磁场值求和，与地磁台站实测数据进行对比，研究了模型数据与实测数据的一致性，为利用CM4模型数据做进一步研究提供依据。

1 数据计算

为了分析研究基于CM4模型的中国大陆地区地磁场时空分布特征，按照经度链和纬度链分别选取4个地磁台站（表1），利用模型计算软件，计算得到了各台站位置处的地核场、地壳场、磁层源磁场及其感应场、电离层源磁场及其感应场、空间环形磁场等7种场源磁场的时均值，由于CM系列模型缺少对观察到的每日磁场变化的表示，本文主要对磁层源磁场及其感应场和电离层源磁场及其感应场的时空分布特征进行研究。

表1 台站分布

2 模型数据各场源时空分布特征分析

2.1 磁层源磁场及其感应场的时空分布特征

CM4模型能够很好地分离地磁场场源，选取表1中固定地磁台1982—2001年间每日世界时第19时磁层源磁场及其感应场时均值进行分析。

2.1.1 经度链数据变化特征

绘制经度链台站地磁场分量磁层源磁场及其感应场每日世界时第19时时均值曲线和F10.7随时间变化曲线图（图1），4个台站时均值曲线图按照台站纬度从高到低的顺序排列，从上到下依次为MZL（a）、LYH（b）、WHN（c）、QZH（d），F10.7（e）是太阳10.7cm波长（2800MHz）射电辐射流量随时间变化曲线，由位于加拿大英属哥伦比亚彭带克顿市的射电天文台每天17时（UTC）对太阳进行观测得到。它的单位是sfu（太阳辐射通量，1sfu＝1×10－22W·m－2·Hz－1）。

从时间上看，4个台站分量磁层源磁场及其感应场变化与太阳活动性指数F10.7的周期性变化相对应，1986年是太阳黑子活动低年，11年之后1997年磁场扰动也较小；相应地，1990年4个台站均出现强扰动，11年之后的2001年同样表现出强烈的磁场变化，存在显著的11年周期。

从空间分布上看，高纬度台站年变化幅度小于低纬度台站，各台站年变化幅度值见表2（其他分量年变化幅度表省略），空间分布特征见图2。

表2 经度链台站X分量年变化幅度

图2显示，随着纬度减小，分量年变化幅度逐渐增大，满洲里台分量年变化幅度最小，泉州台年变化幅度最大，4个台站2001年分量磁场源磁场及其感应场年变化幅度最大，1996年的年变化幅度最小。各台站分量显著相关，且变化幅度随纬度降低而升高，环电流可能是影响其变化的主要因素。

图1 经度链台站磁层源磁场及其感应场X分量时均值随时间变化

经度链台站磁层源磁场及其感应场分量世界时第19时时均值随时间变化特征见图3。

图3所示，地磁场分量的磁层源磁场及其感应场呈现出显著的以1年为周期的一峰两谷规律变化，11年周期性变化不显著。年变化幅度也远远小于分量年变化幅度。

从空间分布上看（图4），相同年份，经度链上4个台站分量年变化幅度基本一致，台站间年变化幅度最大相差不超过3nT。1982年至2001年期间，各台站年变化幅度基本上为9—15nT，但1982年、1990年、2001年除外。台站间分量年变化幅度差别不显著，说明在不同扰动背景下影响分量的电流体系对经度链上台站作用大小是一样的。

图2 经度链台站X分量年变化幅度的空间分布

图3 经度链台站磁层源磁场及其感应场Y分量时均值随时间变化

图4 经度链台站Y分量年变化幅度的空间分布

经度链台站磁层源磁场及其感应场分量世界时第19时时均值随时间变化特征见图5。

图5 经度链台站磁层源磁场及其感应场Z分量时均值随时间变化

如图5所示，磁层源磁场及其感应场的分量与分量具有相同的时间变化特征，但整体变化幅度小于分量。从空间分布看（图6），纬度由高到低，分量年变化幅度从大变小，满洲里台年变化幅度最大，泉州台年变化幅度最小。20年间分量在2001年的年变化幅度最大，1996年最小。

图6 经度链台站Z分量年变化幅度的空间分布

2.1.2 纬度链台站变化特征

采取与分析经度链台站磁层源磁场及其感应场变化特征相同的方法，绘制纬度链台站分量磁层源磁场及其感应场世界时第19时时均值和F10.7随时间变化曲线图（由于篇幅所限，图略）。从时间上来看，纬度链上4个台站分量磁层源磁场及其感应场与经度链上4个台站具有相同变化特征。从空间分布上分析，其变化特征却不同（图7）。

纬度链台站分量磁层源磁场及其感应场年变化幅度基本不随台站的经度改变而变化，2001年各台站年变化幅度相差最大，最大相差4nT。

纬度链台站分量磁层源磁场及其感应场世界时第19时时均值在时间上的变化特征与经度链台站变化特征一致。在空间分布上（图8），除1982年外的其他年份，纬度链台站磁层源磁场及其感应场分量的年变化幅度随台站经度增加而减小，格尔木台年变化幅度最大，泰安台最小。1982年9月6日发生了急始型磁暴，指数最大为9，活动程度为S，红山台（LYH）和泰安台（TAA）世界时第19时时均值变化很大，导致1982年分量磁层源磁场及其感应场年变化幅度变大。格尔木台（GLM）和兰州台（LZH）却没有受到较大影响。

纬度链上台站分量磁层源磁场及其感应场世界时第19时时均值变化特征与经度链台站在时间上的变化特征相同，同样具有显著11年周期变化特征。在空间分布上（图9），分量年变化幅度不受台站经度影响。20年间分量年变化幅度2001年最大，1996年最小。

2.1.3 磁层源磁场及其感应场27天重现性特征

本文选取1997年磁层源磁场及其感应场模型数据时均值，每3个小时的时均值取平均，每天生成8个时均值，与当日的Kp指数相对应，全年每3小时时段均值数据按27天分段，5个27天段对应均值取平均，与对应日期Kp指数进行对比（李细顺等，2015）。研究发现，经度链上4个台站按纬度从高到低顺序排列，依次为MZL（a）、LYH（b）、WHN（c）、QZH（d），分量具有显著的27天周期变化，即地磁活动27天重现性（图10），而分量27天周期性变化则不显著（图11），分量27天周期变化特征与分量一致（图略）。纬度链上4个台站磁层源磁场及其感应场、、三分量具有与经度链台站一样的变化特征。

图7 纬度链台站X分量年变化幅度的空间分布

图8 纬度链台站Y分量年变化幅度的时空分布

图9 纬度链台站Z分量年变化幅度的时空分布

图10 1997年经度链台站X分量的27天周期变化

图11 1997年经度链台站Y分量的27天周期变化

磁层源磁场及其感应场27天重现特征主要受太阳27天自转周影响，地磁活动的27天重现性起源于太阳上的特定磁性区，即太阳冕洞，随着太阳的自转，冕洞周期性面对地球，从而引起磁层源磁场及其感应场27天周期性变化，经度链和纬度链上台站变化特征一致。

2.2 电离层源磁场及其感应场变化特征

2.2.1 经度链台站电离层源磁场及其感应场变化特征

根据WDC中心网站提供的Kp指数，每个月找出5个磁静日，相应地，筛选出1991年经度链上4个台站5个磁静日电离层源磁场及其感应场时均值数据，绘制图12，从上到下按台站纬度从高到低顺序排列，依次为MZL（a）、LYH（b）、WHN（c）、QZH（d）。从图中可以看出，分量、分量和分量呈现出显著的季节变化特征，夏季最大，春秋季次之，冬季最小，3个分量中分量变幅最大，分量次之，分量最小。

从空间分布上看，随着纬度从高到低，分量由小变大，越靠近赤道的台站数值越大。分量却随纬度降低而逐渐减小。分量不随纬度变化而变化。

2.2.2 纬度链台站电离层源磁场及其感应场变化特征

同样，选取1991年纬度链上4个台站每个月5天磁静日电离层源磁场及其感应场时均值数据，绘制图13，4个台站按经度从小到大排列，依次为GLM（a）、LZH（b）、LYH（c）、TAA（d）。从图中可以看出，分量、分量和分量呈现出显著的季节变化特征，夏季最大，春秋季次之，冬季最小，3个分量中分量变幅最大，分量次之，分量最小。

图12 经度链台站1991年电离层源磁场及其感应场变化

在空间分布上，分量、分量、分量值不随台站的经度变化而变化，纬度链上的台站各个分量最大值和最小值基本相同。

采用同样的方法选取1996年数据进行分析研究。1996年是地磁活动低年，磁场比较平静，研究结果显示电离层源磁场及其感应场的时空分布特征与1991年的结果一致，只是各分量变化幅度要小的多，由于篇幅所限，在此不再赘述。

2.3 模型数据与实测数据的一致性分析

上述讨论均针对的是CM4模型计算分离出的独立场源。对软件计算出的CM4模型7种场源值求和，求和后的模型值与台站实测值所包含磁场成分较一致，通过分析模型数据与实测数据一致性，对模型可靠性进行了验证。

分别选取1997年7月（磁场平静月份）和2001年3月（磁场扰动月份）经度链和纬度链7个台站模型值和台站实测值的时均值，计算各台站地磁分量（、、）标准差：

图13 纬度链台站1991年电离层源磁场及其感应场变化

表3 实测值与模型值差值的标准差

从表3中看出，7个台站北向分量的标准差均在33nT之内，东向分量的标准差均在13nT之内，垂直分量的标准差均在35nT之内，各台站差值的标准差相差不大。由于通常认为IGRF模型误差为50—100nT（王亶文，2003；张素琴等，2008），参照判断IGRF模型误差水平的标准，上述标准差结果表明台站实测值与CM4模型计算值的差值比较稳定，一致性较好。

去除台站实测数据和模型数据的系统偏差后，对1997年7月的两种数据进行了对比，如图14所示。图中亮蓝线代表模型数据，黑色散点代表实测数据，红色亮线表示相同时间内的Dst指数（对比图形均采用相同的表示方法）。

图14 1997年7月模型数据与实测数据的时均值变化

图14左侧一列台站按经度由小到大而由上到下顺序排列，右侧一列台站按纬度由高到低由上到下顺序排列，在磁场异常平静的时间段内，模型数据与实测数据吻合非常好，分量的相关系数为0.77—0.89，分量为0.89—0.92，分量为0.85—0.90。

为了进一步研究模型数据与台站实测数据变化的一致性，本文选取了2001年3月的数据进行分析。3月20日发生了特大型磁暴，图15给出了台站实测数据与模型数据在3月期间的变化，高纬度和中低纬度台站在整个磁暴期间均吻合的很好，分量的相关系数为0.95—0.96，分量为0.80—0.89，Z分量为0.70—0.94。

模型误差（或偏差）与资料源和假设条件密切相关。由于CM4模型使用期限为1960—2002年，在此期限内，模型精度较高，磁静期间模型数据与实测数据的一致性较好，即使在磁扰期间，例如2001年3月磁暴期间，模型数据与实测数据的一致性也很好。

3 结论

本文比较系统地研究了基于CM4模型的地磁场在中国大陆地区的时空分布特征，着重研究了磁层源磁场及其感应场和电离层源磁场及其感应场的时空分布，研究结果显示：

（1）经度链和纬度链台站磁层源磁场及其感应场分量、分量和分量存在显著的11年周期变化，分量变化幅度最大，分量次之，分量最小。分量和分量具有明显的27天太阳自转周变化，即地磁活动的27天重现性，而分量27天周期性变化不显著。从空间分布上看，磁层源磁场及其感应场的分量随着纬度减小年变化幅度逐渐增大，却不随台站所处位置经度的改变而变化，磁场扰动年份的年变化幅度最大；经度链上台站的分量年变化幅度基本一致，纬度链台站的分量年变化幅度除1982年外，其他年份年变化幅度随着台站经度增加而减小。分量年变化幅度随着纬度由高到低而从大变小，基本不随台站经度改变而变化。

（2）经度链和纬度链台站电离层源磁场及其感应场季节变化特征一致，分量、分量、分量均呈现出显著季节变化特征，夏季最大，春秋季次之，冬季最小，3个分量中分量变幅最大，分量次之，分量最小。在空间分布上，随着纬度从高到低，分量值由小变大，越靠近赤道数值越大，分量值却随纬度降低而逐渐减小，分量值不随纬度高低的变化而变化。分量、分量、分量的电离层源磁场及其感应场值不受台站经度变化影响，纬度链上台站各个分量的最大值和最小值基本相同。

（3）无论经度链还是纬度链台站，实测时均值数据与模型时均值数据在磁静和磁扰期间均较好吻合，相关系数较高，CM4模型能够较好地反映台站位置处磁场变化特征。

白春华，徐文耀，康国发，2008．地球主磁场模型．地球物理学进展，23（4）：1045—1057．

丁鉴海，索玉成，余素荣，2005．地磁场与电离层异常现象及其与地震的关系．空间科学学报，25（6）：536—542．

冯彦，安振昌，孙涵等，2011．基于第四代地磁场综合模型（CM4）以改善模型边界效应的研究．地球物理学进展，26（3）：850—857．

李细顺，高登平，李琪等，2015．CM4模型数据与台站实测数据的对比研究．震灾防御技术，10（2）：418—425．

李琪，杨星，蔡绍平，2015．极化方法应用于地磁台阵的震例分析．震灾防御技术，10（2）：412—417．

王亶文，2003．国际地磁参考场在中国大陆地区的误差分析．地球物理学报，46（2）：171—174．

徐文耀，2002．地球主磁场的NOC模型．中国科学（D辑），32（7）：576—587．

徐文耀，2009．地球电磁现象物理学．合肥：中国科学技术大学出版社，18—29．

徐文耀，2011．太阳风-磁层-电离层耦合过程中的能量收支．空间科学学报，31（1）：1—14．

杨云存，高国明，2014．1900—2010年地磁场水平分量梯度的全球变化．震灾防御技术，9（S）：557—571．

姚丽，刘振兴，左平兵等，2010．行星际激波对地球磁层的压缩效应分析．空间科学学报，30（2）：113—120．

张素琴，杨冬梅，李琪等，2008．中国部分地磁台站年均值与IGRF模型一致性分析．地震地磁观测与研究，29（2）：42—49．

Alldredge L. R., 1987. On regional magnetic charts. Journal of Geomagnetism and Geoelectricity, 39: 723—738.

Haines G. V., 1985. Spherical cap harmonic analysis. Journal of Geophysical Research: Solid Earth, 90(B3): 2583—2591.

Hemant K., Maus S., 2005. Geological modeling of the new CHAMP magnetic anomaly maps using a geographical information system technique. Journal of Geophysical Research: Solid Earth, 110(B12): B12103.

Hemant K., Thébault E., Mandea M., et al., 2007. Magnetic anomaly map of the world: merging satellite, airborne, marine and ground-based magnetic data sets. Earth and Planetary Science Letters, 260(1—2): 56—71.

Langel R. A., Sabaka T. J., Baldwin R. T., et al., 1996. The near-earth magnetic field from magnetospheric and quiet-day ionospheric sources and how it is modeled. Physics of the Earth and Planetary Interiors, 98(3—4): 235—267.

Lesur V., Wardinski I., Rother M., et al., 2008. GRIMM: the GFZ reference internal magnetic model based on vector satellite and observatory data. Geophysical Journal International, 173(2): 382—394.

Sabaka T. J., Olsen N., Langel R. A., 2002. A comprehensive model of the quiet—time, near—Earth magnetic field: phase 3. Geophysical Journal International, 151(1): 32—68.

Sabaka T. J., Olsen N., Purucker M. E., 2004. Extending comprehensive models of the Earth´s magnetic field with ørsted and CHAMP data. Geophysical Journal International, 159(2): 521—547.

Wardinski I., Holme R., 2006. A time—dependent model of the earth’s magnetic field and its secular variation for the period 1980—2000. Journal of Geophysical Research: Solid Earth, 111(B12): B12101.

Characteristics of Spatial and Temporal Distribution of Geomagnetic Field in Chinese Mainland Based on CM4 Model

Li Xishun, Gao Dengping, Liu Lishen, Zhao Zhiyuan and Wang Libing

(Hongshan Seismological Observatory, Hebei Earthquake Agency, Xingtai 055350, Hebei, China)

In this paper, based on the comprehensive model4 (CM4), we calculated the latitude and longitude chain stations magnetosphere and its induction magnetic field, ionosphere and its induction magnetic field of the,,component from 1982 to 2001 in chinese mainland and also analyzed the change characteristics of the magnetic field with time and space. The results show that the magnetosphere and its induction magnetic field of the longitude chain and latitude chain show periodic variation in 11 years and 27 days. The ionosphere and its induction magnetic field have shown seasonal variation, morphological changes in the same season in which the amplitude is different in the spatial distribution. The annual variation range of the magnetosphere and its induction magnetic field in the longitude chain and latitude chain showed different characteristics. The ionosphere field and its induction magnetic field vary spatially in the longitude chain, but they are about the same of latitude chain station. Diurnal variation analysis showed that the model data is in good consistency with the measured data of the station on the magnetically quiet and disturbed date, and they are highly correlated.

CM4；Magnetosphere field；Ionosphere field；Observatory data；Induction field；The mean time value；Correlation coefficient

李细顺，高登平，刘立申，赵志远，王利兵，2018．基于CM4模型的中国大陆地区地磁场时空分布特征分析．震灾防御技术，13（1）：98—113．

10.11899/zzfy20180109

中国地震局监测、预报、科研三结合课题（160301）

2017-06-16

李细顺，女，生于1976年。高级工程师。主要从事地磁观测与研究。E-mail：270886468@qq.com

高登平，男，生于1973年。高级工程师。主要从事地震监测与研究。E-mail：276221635@qq.com