李細(xì)順 高登平 劉立申 趙志遠(yuǎn) 王利兵

?

基于CM4模型的中國(guó)大陸地區(qū)地磁場(chǎng)時(shí)空分布特征分析1

李細(xì)順 高登平 劉立申 趙志遠(yuǎn) 王利兵

（河北省地震局紅山基準(zhǔn)臺(tái)，河北邢臺(tái) 055350）

本文利用第四代地磁場(chǎng)綜合模型（Comprehensive Model 4，CM4），計(jì)算了1982—2001年中國(guó)大陸地區(qū)同一經(jīng)度鏈和同一緯度鏈上地磁臺(tái)站的磁層源磁場(chǎng)及其感應(yīng)場(chǎng)、電離層源磁場(chǎng)及其感應(yīng)場(chǎng)的地磁北向分量、東向分量、垂直分量的模型值，分析了各場(chǎng)源磁場(chǎng)隨時(shí)間和空間的變化特征。結(jié)果表明：在時(shí)間上，經(jīng)度鏈和緯度鏈臺(tái)站的磁層源磁場(chǎng)及其感應(yīng)場(chǎng)均呈現(xiàn)出11年和27天周期性變化。電離層源磁場(chǎng)及其感應(yīng)場(chǎng)具有明顯的季節(jié)變化，不同年份相同季節(jié)變化形態(tài)一致但幅度不同。在空間分布上，經(jīng)度鏈和緯度鏈臺(tái)站磁層源磁場(chǎng)及其感應(yīng)場(chǎng)的年變化幅度呈現(xiàn)出不同變化特征，電離層源磁場(chǎng)及其感應(yīng)場(chǎng)在經(jīng)度鏈上變化特征不同，而緯度鏈臺(tái)站的數(shù)值基本一致。日變化分析顯示，磁靜日和磁擾日期間，模型數(shù)據(jù)與臺(tái)站實(shí)測(cè)數(shù)據(jù)變化一致性較好，相關(guān)性較高。

CM4磁層源磁場(chǎng) 電離層源磁場(chǎng) 實(shí)測(cè)數(shù)據(jù) 感應(yīng)場(chǎng) 時(shí)均值 相關(guān)系數(shù)

引言

利用地磁場(chǎng)模型研究地磁場(chǎng)變化在空間物理領(lǐng)域一直是一個(gè)基本問題，國(guó)內(nèi)外均有科學(xué)團(tuán)隊(duì)在關(guān)注這一問題。

地磁場(chǎng)是指在固體地球內(nèi)部和外部到磁層頂空間范圍內(nèi)所有場(chǎng)源產(chǎn)生的磁場(chǎng)，一般將地磁場(chǎng)分為內(nèi)源場(chǎng)和外源場(chǎng)，人們對(duì)于內(nèi)源場(chǎng)研究較多（Haines，1985；Alldredge，1987；徐文耀，2002，2009；Wardinski等，2006；Lesur等，2008），對(duì)地殼磁異常場(chǎng)也做了很多工作（Hemant等，2005，2007）。地磁模型能對(duì)地磁場(chǎng)時(shí)空分布進(jìn)行較為直觀的表示。1968年國(guó)際地磁和高空大氣物理學(xué)協(xié)會(huì)給出了第一代國(guó)際地磁參考場(chǎng)模型（IGRF），之后每5年給出一個(gè)新的IGRF（楊云存等，2014）。20世紀(jì)90年代，美國(guó)國(guó)家宇航局戈達(dá)德飛行中心（NASA/GSFC）和丹麥空間研究所（DSRI）聯(lián)合開發(fā)了一種新的地磁場(chǎng)建模方法——“綜合建模”（Comprehensive Modeling，簡(jiǎn)稱CM）（白春華等，2008）。1993年，Sabaka等建立了第一代地磁場(chǎng)綜合模型——CM1，1996年，Langel等在CM1的基礎(chǔ)上建立了第二代地磁場(chǎng)綜合模型——CM2（Langel等，1996），2002年，美國(guó)和丹麥空間研究中心的科學(xué)家Sabaka和Olsen等人建立了第三代地磁場(chǎng)綜合模型——CM3（Sabaka等，2002），它相對(duì)于第一代和第二代的提高主要體現(xiàn)在進(jìn)行了高階地殼層場(chǎng)的靜態(tài)表示，能在衛(wèi)星高度獲取大部分地殼異常場(chǎng)。CM3模型內(nèi)源場(chǎng)最大截?cái)嚯A數(shù)為65，時(shí)間跨度為1960—1985年，長(zhǎng)期變化通過3次B-spline方法的13階表示（馮彥等，2011）。

2004年，Sabaka等利用衛(wèi)星和臺(tái)站數(shù)據(jù)在CM3的基礎(chǔ)上建立了第四代地磁場(chǎng)綜合模型——CM4（Sabaka等，2004），相較于CM3，CM4采用了更多衛(wèi)星數(shù)據(jù)，包括POGO和CHAMP衛(wèi)星標(biāo)量數(shù)據(jù)，MAGSAT和Frsted衛(wèi)星標(biāo)量和矢量數(shù)據(jù)。由于衛(wèi)星位于電離層和磁層間的耦合電流區(qū)域（用于區(qū)分電離層和磁層的磁場(chǎng)強(qiáng)度），故可獲取近地球區(qū)域各種場(chǎng)的電流源，并進(jìn)行各種場(chǎng)的參數(shù)化和協(xié)同估算，通過迭代重加權(quán)最小二乘法進(jìn)行合理劃分，最后順序計(jì)算7種磁場(chǎng)的值（地核場(chǎng)、地殼場(chǎng)、磁層源場(chǎng)、磁層源感應(yīng)場(chǎng)、電離層源場(chǎng)、電離層源感應(yīng)場(chǎng)、空間環(huán)形磁場(chǎng)）。CM4模型適用于磁靜日期間，可對(duì)地磁場(chǎng)進(jìn)行分層計(jì)算，相較CM3模型，其適用時(shí)間范圍更長(zhǎng)，為1960—2002年，進(jìn)一步降低了預(yù)處理時(shí)的數(shù)據(jù)噪聲。模型1—15階代表地核場(chǎng)，16—65階代表地殼場(chǎng)。

地磁臺(tái)站記錄到的地磁變化信號(hào)主要受到外源場(chǎng)和內(nèi)源場(chǎng)兩種場(chǎng)源的影響。其中，外源場(chǎng)起源于地表以上的空間電流體系，這些電流體系主要分布在電離層、磁層和行星際空間（李琪等，2015）。我國(guó)基于臺(tái)站和衛(wèi)星數(shù)據(jù)對(duì)電離層場(chǎng)以及磁層場(chǎng)開展了很多研究。丁鑒海等（2005）利用不同地區(qū)電離層臺(tái)站數(shù)據(jù)對(duì)電離層2層的臨界頻率02進(jìn)行了相關(guān)研究；姚麗等（2010）通過分析位于向陽(yáng)面正午兩側(cè)的GOES-10和GOES-12衛(wèi)星觀測(cè)數(shù)據(jù)發(fā)現(xiàn)，激波作用于磁層時(shí)靠近清晨的磁場(chǎng)變化表現(xiàn)出簡(jiǎn)單壓縮效應(yīng)，而靠近黃昏的磁場(chǎng)變化則顯然不同，即B分量減弱，B分量幾乎為零，B分量則顯著增強(qiáng)；徐文耀（2011）從能量守恒原理出發(fā)，討論了太陽(yáng)風(fēng)-磁層-電離層耦合過程的能流路徑和能量收支的定量關(guān)系問題。本文利用第四代地磁場(chǎng)綜合模型（CM4），計(jì)算了1982—2001年中國(guó)大陸地區(qū)經(jīng)度鏈和緯度鏈上地磁臺(tái)站的磁層源磁場(chǎng)及其感應(yīng)場(chǎng)、電離層源磁場(chǎng)及其感應(yīng)場(chǎng)的地磁北向分量、東向分量和垂直分量，分析了各場(chǎng)源磁場(chǎng)隨時(shí)間和空間的變化特征。對(duì)CM4模型計(jì)算出的7種磁場(chǎng)值求和，與地磁臺(tái)站實(shí)測(cè)數(shù)據(jù)進(jìn)行對(duì)比，研究了模型數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)的一致性，為利用CM4模型數(shù)據(jù)做進(jìn)一步研究提供依據(jù)。

1 數(shù)據(jù)計(jì)算

為了分析研究基于CM4模型的中國(guó)大陸地區(qū)地磁場(chǎng)時(shí)空分布特征，按照經(jīng)度鏈和緯度鏈分別選取4個(gè)地磁臺(tái)站（表1），利用模型計(jì)算軟件，計(jì)算得到了各臺(tái)站位置處的地核場(chǎng)、地殼場(chǎng)、磁層源磁場(chǎng)及其感應(yīng)場(chǎng)、電離層源磁場(chǎng)及其感應(yīng)場(chǎng)、空間環(huán)形磁場(chǎng)等7種場(chǎng)源磁場(chǎng)的時(shí)均值，由于CM系列模型缺少對(duì)觀察到的每日磁場(chǎng)變化的表示，本文主要對(duì)磁層源磁場(chǎng)及其感應(yīng)場(chǎng)和電離層源磁場(chǎng)及其感應(yīng)場(chǎng)的時(shí)空分布特征進(jìn)行研究。

表1 臺(tái)站分布

2 模型數(shù)據(jù)各場(chǎng)源時(shí)空分布特征分析

2.1 磁層源磁場(chǎng)及其感應(yīng)場(chǎng)的時(shí)空分布特征

CM4模型能夠很好地分離地磁場(chǎng)場(chǎng)源，選取表1中固定地磁臺(tái)1982—2001年間每日世界時(shí)第19時(shí)磁層源磁場(chǎng)及其感應(yīng)場(chǎng)時(shí)均值進(jìn)行分析。

2.1.1 經(jīng)度鏈數(shù)據(jù)變化特征

繪制經(jīng)度鏈臺(tái)站地磁場(chǎng)分量磁層源磁場(chǎng)及其感應(yīng)場(chǎng)每日世界時(shí)第19時(shí)時(shí)均值曲線和F10.7隨時(shí)間變化曲線圖（圖1），4個(gè)臺(tái)站時(shí)均值曲線圖按照臺(tái)站緯度從高到低的順序排列，從上到下依次為MZL（a）、LYH（b）、WHN（c）、QZH（d），F(xiàn)10.7（e）是太陽(yáng)10.7cm波長(zhǎng)（2800MHz）射電輻射流量隨時(shí)間變化曲線，由位于加拿大英屬哥倫比亞彭帶克頓市的射電天文臺(tái)每天17時(shí)（UTC）對(duì)太陽(yáng)進(jìn)行觀測(cè)得到。它的單位是sfu（太陽(yáng)輻射通量，1sfu＝1×10－22W·m－2·Hz－1）。

從時(shí)間上看，4個(gè)臺(tái)站分量磁層源磁場(chǎng)及其感應(yīng)場(chǎng)變化與太陽(yáng)活動(dòng)性指數(shù)F10.7的周期性變化相對(duì)應(yīng)，1986年是太陽(yáng)黑子活動(dòng)低年，11年之后1997年磁場(chǎng)擾動(dòng)也較小；相應(yīng)地，1990年4個(gè)臺(tái)站均出現(xiàn)強(qiáng)擾動(dòng)，11年之后的2001年同樣表現(xiàn)出強(qiáng)烈的磁場(chǎng)變化，存在顯著的11年周期。

從空間分布上看，高緯度臺(tái)站年變化幅度小于低緯度臺(tái)站，各臺(tái)站年變化幅度值見表2（其他分量年變化幅度表省略），空間分布特征見圖2。

表2 經(jīng)度鏈臺(tái)站X分量年變化幅度

圖2顯示，隨著緯度減小，分量年變化幅度逐漸增大，滿洲里臺(tái)分量年變化幅度最小，泉州臺(tái)年變化幅度最大，4個(gè)臺(tái)站2001年分量磁場(chǎng)源磁場(chǎng)及其感應(yīng)場(chǎng)年變化幅度最大，1996年的年變化幅度最小。各臺(tái)站分量顯著相關(guān)，且變化幅度隨緯度降低而升高，環(huán)電流可能是影響其變化的主要因素。

圖1 經(jīng)度鏈臺(tái)站磁層源磁場(chǎng)及其感應(yīng)場(chǎng)X分量時(shí)均值隨時(shí)間變化

經(jīng)度鏈臺(tái)站磁層源磁場(chǎng)及其感應(yīng)場(chǎng)分量世界時(shí)第19時(shí)時(shí)均值隨時(shí)間變化特征見圖3。

圖3所示，地磁場(chǎng)分量的磁層源磁場(chǎng)及其感應(yīng)場(chǎng)呈現(xiàn)出顯著的以1年為周期的一峰兩谷規(guī)律變化，11年周期性變化不顯著。年變化幅度也遠(yuǎn)遠(yuǎn)小于分量年變化幅度。

從空間分布上看（圖4），相同年份，經(jīng)度鏈上4個(gè)臺(tái)站分量年變化幅度基本一致，臺(tái)站間年變化幅度最大相差不超過3nT。1982年至2001年期間，各臺(tái)站年變化幅度基本上為9—15nT，但1982年、1990年、2001年除外。臺(tái)站間分量年變化幅度差別不顯著，說明在不同擾動(dòng)背景下影響分量的電流體系對(duì)經(jīng)度鏈上臺(tái)站作用大小是一樣的。

圖2 經(jīng)度鏈臺(tái)站X分量年變化幅度的空間分布

圖3 經(jīng)度鏈臺(tái)站磁層源磁場(chǎng)及其感應(yīng)場(chǎng)Y分量時(shí)均值隨時(shí)間變化

圖4 經(jīng)度鏈臺(tái)站Y分量年變化幅度的空間分布

經(jīng)度鏈臺(tái)站磁層源磁場(chǎng)及其感應(yīng)場(chǎng)分量世界時(shí)第19時(shí)時(shí)均值隨時(shí)間變化特征見圖5。

圖5 經(jīng)度鏈臺(tái)站磁層源磁場(chǎng)及其感應(yīng)場(chǎng)Z分量時(shí)均值隨時(shí)間變化

如圖5所示，磁層源磁場(chǎng)及其感應(yīng)場(chǎng)的分量與分量具有相同的時(shí)間變化特征，但整體變化幅度小于分量。從空間分布看（圖6），緯度由高到低，分量年變化幅度從大變小，滿洲里臺(tái)年變化幅度最大，泉州臺(tái)年變化幅度最小。20年間分量在2001年的年變化幅度最大，1996年最小。

圖6 經(jīng)度鏈臺(tái)站Z分量年變化幅度的空間分布

2.1.2 緯度鏈臺(tái)站變化特征

采取與分析經(jīng)度鏈臺(tái)站磁層源磁場(chǎng)及其感應(yīng)場(chǎng)變化特征相同的方法，繪制緯度鏈臺(tái)站分量磁層源磁場(chǎng)及其感應(yīng)場(chǎng)世界時(shí)第19時(shí)時(shí)均值和F10.7隨時(shí)間變化曲線圖（由于篇幅所限，圖略）。從時(shí)間上來看，緯度鏈上4個(gè)臺(tái)站分量磁層源磁場(chǎng)及其感應(yīng)場(chǎng)與經(jīng)度鏈上4個(gè)臺(tái)站具有相同變化特征。從空間分布上分析，其變化特征卻不同（圖7）。

緯度鏈臺(tái)站分量磁層源磁場(chǎng)及其感應(yīng)場(chǎng)年變化幅度基本不隨臺(tái)站的經(jīng)度改變而變化，2001年各臺(tái)站年變化幅度相差最大，最大相差4nT。

緯度鏈臺(tái)站分量磁層源磁場(chǎng)及其感應(yīng)場(chǎng)世界時(shí)第19時(shí)時(shí)均值在時(shí)間上的變化特征與經(jīng)度鏈臺(tái)站變化特征一致。在空間分布上（圖8），除1982年外的其他年份，緯度鏈臺(tái)站磁層源磁場(chǎng)及其感應(yīng)場(chǎng)分量的年變化幅度隨臺(tái)站經(jīng)度增加而減小，格爾木臺(tái)年變化幅度最大，泰安臺(tái)最小。1982年9月6日發(fā)生了急始型磁暴，指數(shù)最大為9，活動(dòng)程度為S，紅山臺(tái)（LYH）和泰安臺(tái)（TAA）世界時(shí)第19時(shí)時(shí)均值變化很大，導(dǎo)致1982年分量磁層源磁場(chǎng)及其感應(yīng)場(chǎng)年變化幅度變大。格爾木臺(tái)（GLM）和蘭州臺(tái)（LZH）卻沒有受到較大影響。

緯度鏈上臺(tái)站分量磁層源磁場(chǎng)及其感應(yīng)場(chǎng)世界時(shí)第19時(shí)時(shí)均值變化特征與經(jīng)度鏈臺(tái)站在時(shí)間上的變化特征相同，同樣具有顯著11年周期變化特征。在空間分布上（圖9），分量年變化幅度不受臺(tái)站經(jīng)度影響。20年間分量年變化幅度2001年最大，1996年最小。

2.1.3 磁層源磁場(chǎng)及其感應(yīng)場(chǎng)27天重現(xiàn)性特征

本文選取1997年磁層源磁場(chǎng)及其感應(yīng)場(chǎng)模型數(shù)據(jù)時(shí)均值，每3個(gè)小時(shí)的時(shí)均值取平均，每天生成8個(gè)時(shí)均值，與當(dāng)日的Kp指數(shù)相對(duì)應(yīng)，全年每3小時(shí)時(shí)段均值數(shù)據(jù)按27天分段，5個(gè)27天段對(duì)應(yīng)均值取平均，與對(duì)應(yīng)日期Kp指數(shù)進(jìn)行對(duì)比（李細(xì)順等，2015）。研究發(fā)現(xiàn)，經(jīng)度鏈上4個(gè)臺(tái)站按緯度從高到低順序排列，依次為MZL（a）、LYH（b）、WHN（c）、QZH（d），分量具有顯著的27天周期變化，即地磁活動(dòng)27天重現(xiàn)性（圖10），而分量27天周期性變化則不顯著（圖11），分量27天周期變化特征與分量一致（圖略）。緯度鏈上4個(gè)臺(tái)站磁層源磁場(chǎng)及其感應(yīng)場(chǎng)、、三分量具有與經(jīng)度鏈臺(tái)站一樣的變化特征。

圖7 緯度鏈臺(tái)站X分量年變化幅度的空間分布

圖8 緯度鏈臺(tái)站Y分量年變化幅度的時(shí)空分布

圖9 緯度鏈臺(tái)站Z分量年變化幅度的時(shí)空分布

圖10 1997年經(jīng)度鏈臺(tái)站X分量的27天周期變化

圖11 1997年經(jīng)度鏈臺(tái)站Y分量的27天周期變化

磁層源磁場(chǎng)及其感應(yīng)場(chǎng)27天重現(xiàn)特征主要受太陽(yáng)27天自轉(zhuǎn)周影響，地磁活動(dòng)的27天重現(xiàn)性起源于太陽(yáng)上的特定磁性區(qū)，即太陽(yáng)冕洞，隨著太陽(yáng)的自轉(zhuǎn)，冕洞周期性面對(duì)地球，從而引起磁層源磁場(chǎng)及其感應(yīng)場(chǎng)27天周期性變化，經(jīng)度鏈和緯度鏈上臺(tái)站變化特征一致。

2.2 電離層源磁場(chǎng)及其感應(yīng)場(chǎng)變化特征

2.2.1 經(jīng)度鏈臺(tái)站電離層源磁場(chǎng)及其感應(yīng)場(chǎng)變化特征

根據(jù)WDC中心網(wǎng)站提供的Kp指數(shù)，每個(gè)月找出5個(gè)磁靜日，相應(yīng)地，篩選出1991年經(jīng)度鏈上4個(gè)臺(tái)站5個(gè)磁靜日電離層源磁場(chǎng)及其感應(yīng)場(chǎng)時(shí)均值數(shù)據(jù)，繪制圖12，從上到下按臺(tái)站緯度從高到低順序排列，依次為MZL（a）、LYH（b）、WHN（c）、QZH（d）。從圖中可以看出，分量、分量和分量呈現(xiàn)出顯著的季節(jié)變化特征，夏季最大，春秋季次之，冬季最小，3個(gè)分量中分量變幅最大，分量次之，分量最小。

從空間分布上看，隨著緯度從高到低，分量由小變大，越靠近赤道的臺(tái)站數(shù)值越大。分量卻隨緯度降低而逐漸減小。分量不隨緯度變化而變化。

2.2.2 緯度鏈臺(tái)站電離層源磁場(chǎng)及其感應(yīng)場(chǎng)變化特征

同樣，選取1991年緯度鏈上4個(gè)臺(tái)站每個(gè)月5天磁靜日電離層源磁場(chǎng)及其感應(yīng)場(chǎng)時(shí)均值數(shù)據(jù)，繪制圖13，4個(gè)臺(tái)站按經(jīng)度從小到大排列，依次為GLM（a）、LZH（b）、LYH（c）、TAA（d）。從圖中可以看出，分量、分量和分量呈現(xiàn)出顯著的季節(jié)變化特征，夏季最大，春秋季次之，冬季最小，3個(gè)分量中分量變幅最大，分量次之，分量最小。

圖12 經(jīng)度鏈臺(tái)站1991年電離層源磁場(chǎng)及其感應(yīng)場(chǎng)變化

在空間分布上，分量、分量、分量值不隨臺(tái)站的經(jīng)度變化而變化，緯度鏈上的臺(tái)站各個(gè)分量最大值和最小值基本相同。

采用同樣的方法選取1996年數(shù)據(jù)進(jìn)行分析研究。1996年是地磁活動(dòng)低年，磁場(chǎng)比較平靜，研究結(jié)果顯示電離層源磁場(chǎng)及其感應(yīng)場(chǎng)的時(shí)空分布特征與1991年的結(jié)果一致，只是各分量變化幅度要小的多，由于篇幅所限，在此不再贅述。

2.3 模型數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)的一致性分析

上述討論均針對(duì)的是CM4模型計(jì)算分離出的獨(dú)立場(chǎng)源。對(duì)軟件計(jì)算出的CM4模型7種場(chǎng)源值求和，求和后的模型值與臺(tái)站實(shí)測(cè)值所包含磁場(chǎng)成分較一致，通過分析模型數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)一致性，對(duì)模型可靠性進(jìn)行了驗(yàn)證。

分別選取1997年7月（磁場(chǎng)平靜月份）和2001年3月（磁場(chǎng)擾動(dòng)月份）經(jīng)度鏈和緯度鏈7個(gè)臺(tái)站模型值和臺(tái)站實(shí)測(cè)值的時(shí)均值，計(jì)算各臺(tái)站地磁分量（、、）標(biāo)準(zhǔn)差：

圖13 緯度鏈臺(tái)站1991年電離層源磁場(chǎng)及其感應(yīng)場(chǎng)變化

表3 實(shí)測(cè)值與模型值差值的標(biāo)準(zhǔn)差

從表3中看出，7個(gè)臺(tái)站北向分量的標(biāo)準(zhǔn)差均在33nT之內(nèi)，東向分量的標(biāo)準(zhǔn)差均在13nT之內(nèi)，垂直分量的標(biāo)準(zhǔn)差均在35nT之內(nèi)，各臺(tái)站差值的標(biāo)準(zhǔn)差相差不大。由于通常認(rèn)為IGRF模型誤差為50—100nT（王亶文，2003；張素琴等，2008），參照判斷IGRF模型誤差水平的標(biāo)準(zhǔn)，上述標(biāo)準(zhǔn)差結(jié)果表明臺(tái)站實(shí)測(cè)值與CM4模型計(jì)算值的差值比較穩(wěn)定，一致性較好。

去除臺(tái)站實(shí)測(cè)數(shù)據(jù)和模型數(shù)據(jù)的系統(tǒng)偏差后，對(duì)1997年7月的兩種數(shù)據(jù)進(jìn)行了對(duì)比，如圖14所示。圖中亮藍(lán)線代表模型數(shù)據(jù)，黑色散點(diǎn)代表實(shí)測(cè)數(shù)據(jù)，紅色亮線表示相同時(shí)間內(nèi)的Dst指數(shù)（對(duì)比圖形均采用相同的表示方法）。

圖14 1997年7月模型數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)的時(shí)均值變化

圖14左側(cè)一列臺(tái)站按經(jīng)度由小到大而由上到下順序排列，右側(cè)一列臺(tái)站按緯度由高到低由上到下順序排列，在磁場(chǎng)異常平靜的時(shí)間段內(nèi)，模型數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)吻合非常好，分量的相關(guān)系數(shù)為0.77—0.89，分量為0.89—0.92，分量為0.85—0.90。

為了進(jìn)一步研究模型數(shù)據(jù)與臺(tái)站實(shí)測(cè)數(shù)據(jù)變化的一致性，本文選取了2001年3月的數(shù)據(jù)進(jìn)行分析。3月20日發(fā)生了特大型磁暴，圖15給出了臺(tái)站實(shí)測(cè)數(shù)據(jù)與模型數(shù)據(jù)在3月期間的變化，高緯度和中低緯度臺(tái)站在整個(gè)磁暴期間均吻合的很好，分量的相關(guān)系數(shù)為0.95—0.96，分量為0.80—0.89，Z分量為0.70—0.94。

模型誤差（或偏差）與資料源和假設(shè)條件密切相關(guān)。由于CM4模型使用期限為1960—2002年，在此期限內(nèi)，模型精度較高，磁靜期間模型數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)的一致性較好，即使在磁擾期間，例如2001年3月磁暴期間，模型數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)的一致性也很好。

3 結(jié)論

本文比較系統(tǒng)地研究了基于CM4模型的地磁場(chǎng)在中國(guó)大陸地區(qū)的時(shí)空分布特征，著重研究了磁層源磁場(chǎng)及其感應(yīng)場(chǎng)和電離層源磁場(chǎng)及其感應(yīng)場(chǎng)的時(shí)空分布，研究結(jié)果顯示：

（1）經(jīng)度鏈和緯度鏈臺(tái)站磁層源磁場(chǎng)及其感應(yīng)場(chǎng)分量、分量和分量存在顯著的11年周期變化，分量變化幅度最大，分量次之，分量最小。分量和分量具有明顯的27天太陽(yáng)自轉(zhuǎn)周變化，即地磁活動(dòng)的27天重現(xiàn)性，而分量27天周期性變化不顯著。從空間分布上看，磁層源磁場(chǎng)及其感應(yīng)場(chǎng)的分量隨著緯度減小年變化幅度逐漸增大，卻不隨臺(tái)站所處位置經(jīng)度的改變而變化，磁場(chǎng)擾動(dòng)年份的年變化幅度最大；經(jīng)度鏈上臺(tái)站的分量年變化幅度基本一致，緯度鏈臺(tái)站的分量年變化幅度除1982年外，其他年份年變化幅度隨著臺(tái)站經(jīng)度增加而減小。分量年變化幅度隨著緯度由高到低而從大變小，基本不隨臺(tái)站經(jīng)度改變而變化。

（2）經(jīng)度鏈和緯度鏈臺(tái)站電離層源磁場(chǎng)及其感應(yīng)場(chǎng)季節(jié)變化特征一致，分量、分量、分量均呈現(xiàn)出顯著季節(jié)變化特征，夏季最大，春秋季次之，冬季最小，3個(gè)分量中分量變幅最大，分量次之，分量最小。在空間分布上，隨著緯度從高到低，分量值由小變大，越靠近赤道數(shù)值越大，分量值卻隨緯度降低而逐漸減小，分量值不隨緯度高低的變化而變化。分量、分量、分量的電離層源磁場(chǎng)及其感應(yīng)場(chǎng)值不受臺(tái)站經(jīng)度變化影響，緯度鏈上臺(tái)站各個(gè)分量的最大值和最小值基本相同。

（3）無論經(jīng)度鏈還是緯度鏈臺(tái)站，實(shí)測(cè)時(shí)均值數(shù)據(jù)與模型時(shí)均值數(shù)據(jù)在磁靜和磁擾期間均較好吻合，相關(guān)系數(shù)較高，CM4模型能夠較好地反映臺(tái)站位置處磁場(chǎng)變化特征。

白春華，徐文耀，康國(guó)發(fā)，2008．地球主磁場(chǎng)模型．地球物理學(xué)進(jìn)展，23（4）：1045—1057．

丁鑒海，索玉成，余素榮，2005．地磁場(chǎng)與電離層異常現(xiàn)象及其與地震的關(guān)系．空間科學(xué)學(xué)報(bào)，25（6）：536—542．

馮彥，安振昌，孫涵等，2011．基于第四代地磁場(chǎng)綜合模型（CM4）以改善模型邊界效應(yīng)的研究．地球物理學(xué)進(jìn)展，26（3）：850—857．

李細(xì)順，高登平，李琪等，2015．CM4模型數(shù)據(jù)與臺(tái)站實(shí)測(cè)數(shù)據(jù)的對(duì)比研究．震災(zāi)防御技術(shù)，10（2）：418—425．

李琪，楊星，蔡紹平，2015．極化方法應(yīng)用于地磁臺(tái)陣的震例分析．震災(zāi)防御技術(shù)，10（2）：412—417．

王亶文，2003．國(guó)際地磁參考場(chǎng)在中國(guó)大陸地區(qū)的誤差分析．地球物理學(xué)報(bào)，46（2）：171—174．

徐文耀，2002．地球主磁場(chǎng)的NOC模型．中國(guó)科學(xué)（D輯），32（7）：576—587．

徐文耀，2009．地球電磁現(xiàn)象物理學(xué)．合肥：中國(guó)科學(xué)技術(shù)大學(xué)出版社，18—29．

徐文耀，2011．太陽(yáng)風(fēng)-磁層-電離層耦合過程中的能量收支．空間科學(xué)學(xué)報(bào)，31（1）：1—14．

楊云存，高國(guó)明，2014．1900—2010年地磁場(chǎng)水平分量梯度的全球變化．震災(zāi)防御技術(shù)，9（S）：557—571．

姚麗，劉振興，左平兵等，2010．行星際激波對(duì)地球磁層的壓縮效應(yīng)分析．空間科學(xué)學(xué)報(bào)，30（2）：113—120．

張素琴，楊冬梅，李琪等，2008．中國(guó)部分地磁臺(tái)站年均值與IGRF模型一致性分析．地震地磁觀測(cè)與研究，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

李細(xì)順，高登平，劉立申，趙志遠(yuǎn)，王利兵，2018．基于CM4模型的中國(guó)大陸地區(qū)地磁場(chǎng)時(shí)空分布特征分析．震災(zāi)防御技術(shù)，13（1）：98—113．

10.11899/zzfy20180109

中國(guó)地震局監(jiān)測(cè)、預(yù)報(bào)、科研三結(jié)合課題（160301）

2017-06-16

李細(xì)順，女，生于1976年。高級(jí)工程師。主要從事地磁觀測(cè)與研究。E-mail：270886468@qq.com

高登平，男，生于1973年。高級(jí)工程師。主要從事地震監(jiān)測(cè)與研究。E-mail：276221635@qq.com