太陽黑子、地球自轉與華北地區中強震關系研究

裴東洋，高立新，戴 勇，格 根，梁沙沙，王 磊

(內蒙古自治區地震局，內蒙古 呼和浩特 010010)

0 引言

太陽是距離地球最近的恒星，許多地球物理過程和太陽活動有密切關系[1]，因此研究太陽活動與地震的關系具有重要意義。天文學上一般用太陽黑子相對數作為描述太陽活動強度的指數。地震與太陽黑子相對數的統計關系早已受到研究學者的關注。目前，關于太陽黑子數變化與地震活動頻次的相關性研究認為兩者呈一定程度的負相關，太陽黑子的強烈活動是中強地震發生的重要誘因和影響因子，并且這種影響具有一定的延遲效果。有學者認為地震一般發生在太陽11年活動周期減弱時期[2]。此外，也有學者認為不同構造區域的地震與太陽黑子活動的關系不同，且同一斷裂帶的不同地段也有不同，有的地區集中在太陽活動的下降段發震，有的則在太陽活動上升段發震[3]。

該文主要針對太陽黑子活動周期及地球自轉速率變化與華北地區5級以上中強地震活動之間的關系進行分析研究。研究區的范圍參考華北地區地震活動斷裂邊界[4]，北起陰山褶皺帶，南至河南許昌，西起內蒙古五原，東至遼寧丹東(34°～42°N，108°～124°E)。地震目錄來源于中國地震臺網中心統一編目后的目錄，考慮到7級以上強震造成的余震與地球外部因素關系不大，使用Gardner和Knopoff提出的方法刪除目錄中7級以上強震的余震[5]。太陽黑子相對數來源于比利時皇家天文臺世界數據中心SILSO[6]，時間從1749年至2019年，并進行13個月滑動平均處理得到太陽黑子相對數月滑動平均值。地球自轉速率變化以國際地球自轉服務局(IERS,http://hpiers.obspm.fr/)提供的地球日長時間(LOD)數據來表征地球自轉速度的大小，LOD數據表示地球日長時間實際值與24小時差值。

1 太陽黑子周期特征

太陽黑子相對數存在平均11 a的周期變化規律。按照國際規定，把1755年為太陽黑子活動周第一周的起點，至今已經歷24個完整周期，從2019年3月進入第25個太陽活動周。根據黑子數變化情況，每一個完整的太陽黑子活動周可分為上升段和下降段，通過分析1901年以來太陽黑子相對數變化情況得知，上升段持續約3～4 a，下降段持續約7～8 a。不同活動周期內太陽黑子峰值數及時序曲線形態能夠反映太陽活動的強度。太陽活動強烈的周期內一般表現為單峰值曲線，黑子數峰值均在200以上；太陽活動平穩的周期內一般表現為雙峰值曲線，黑子數峰值在100左右(見第2頁圖1、表1)。

小波變換可以將時間序列分解成時間-頻率形態，使數據在時域與頻域方面均有良好局部性質，可以檢測時間序列中的突變信息，分析它的主要變化周期，并得到這些周期在不同時段的影響大小，因此適合研究太陽黑子相對數時間序列的周期性變化[7]。對太陽黑子相對數作Morlet小波分析，用小波系數來反映全部周期的演化特征及主周期的變化趨勢。

用小波系數的實部等值線圖來反映太陽黑子相對數序列不同時間尺度的周期變化及其在時間域中的分布，判斷在不同時間尺度上，黑子數的未來變化趨勢。從實部等值線圖看出，太陽黑子相對數變化在9～16 a的時間段上顯著(見圖2)。用小波系數模值來反映不同時間尺度變化周期所對應的能量密度在時間域中分布的反映，系數模值愈大，表明其所對應時間尺度的周期性愈強。從模值等值線圖看出，在太陽黑子相對數變化過程中，11～16 a時間尺度模值最大，說明該時間尺度周期變化最明顯(見圖3)。

表1 1901年以來太陽黑子活動周黑子相對數特征表Table 1 Characteristics of relative number of sunspots in the active cycle of sunspots since 1901

根據以上小波系數的分析結果，用小波方差圖來反映太陽黑子相對數時間序列的波動能量隨周期的分布情況，確定黑子數變化過程中的主周期。由小波方差圖看出(見圖4)，在1～32 a尺度上的小波方差檢驗結果中，只有11 a左右的周期顯著，是太陽黑子相對數變化的主周期。因此，對于黑子數變化趨勢的分析應以11 a周期的小波數據為基礎。在主周期小波系數實部圖中看出(見第3頁圖5)，1930-2000年的小波系數實部值較高，反映太陽黑子活動劇烈。2000年以后活動減弱，由主周期小波系數實部曲線變化規律推測，這一過程持續約20～25 a左右。

基于太陽黑子相對數時間序列小波系數、小波方差檢驗結果以及主周期小波系數實部變化曲線的分析，認為2000年以來太陽黑子相對數處于相對低水平。

圖2 太陽黑子小波系數實部等值線圖Fig.2 Contour map of real part of wavelet coefficient of sunspot

圖3 小波系數模值等值線圖Fig.3 Modulus isoline of wavelet coefficient

圖4 太陽黑子數小波方差圖Fig.4 Wavelet variance of sunspot number

圖5 太陽黑子數主周期小波系數圖Fig.5 Main period wavelet coefficients of sunspot number

當前的第25太陽活動周，太陽黑子活動將進一步加強，太陽黑子相對數將處于相對高的水平。

2 華北地區中強震與太陽活動關系

太陽活動對地震活動的作用有兩種方式。一種是通過電磁效應來實現[8]，太陽活動造成的磁暴對地球磁場影響強烈，而地球磁場變化可能通過磁震效應產生觸發地震活動的效果；另一種是通過太陽物質噴發造成日地間引力變化影響地球自轉，進而作用于地震孕育過程[9]。

2.1 太陽黑子相對數與中強震頻次變化關系

按照太陽黑子活動周各變化階段統計華北地區中強震的發生情況看出(見圖6)，太陽活動是誘發中強地震活動的重要因素，且影響效果存在明顯的延遲。華北地區5級以上中強震一般不發生在太陽黑子活動最劇烈的年份，由于空間距離的影響、地殼運動的累積效應，中強震活動強烈的年份一般為太陽黑子數前一個周期波谷與下一個周期波峰之間，太陽黑子數處于上升階段的區間。太陽黑子活動與地震發生之間的時間延遲現象，說明太陽活動并不能直接觸發地震活動。

圖6 太陽黑子相對數時序曲線與華北地區中強地震頻次Fig.6 Time series curve of sunspot relative number and frequency of moderate and strong earthquakes in North China

2.2 地球自轉與太陽黑子、中強地震時序列關系

太陽活動劇烈期間將向外太空拋射大量的物質，導致太陽扁率增加。由于日地間引力關系和引力攝動效應，造成地球自轉速率減小[10]。由于太陽輻射光壓在南、北半球表面上的不平衡分布以及地球軌道變化，太陽輻射光壓可產生地球內核相對地殼地幔的南北震蕩，是形成地球自轉速度周期變化的主要原因[10]。地球自轉速率的變化將導致地球內部圈層力學耦合狀態以及地球物理場的變化，導致地殼，特別是斷裂帶上應力的重新分布[11-12]。

太陽黑子相對數與地球自轉速率呈明顯的負相關關系。太陽黑子相對數上升階段，地球自轉變慢、日長增加，華北地區中強地震活躍。太陽黑子數量達到峰值后約兩年后，地球自轉速度有加快趨勢；太陽黑子數量達到谷值年后約2年，地球自轉速度呈降低趨勢。目前，華北地區地震活動處于第4活動期，已經歷5個活躍幕[13]。其中，第一個強震活躍期大約為1966年到1977年，發生了河北邢臺、遼寧海城和河北唐山3次7級以上強震；第二個活躍期為1988年到1998年，發生了山西大同-陽高、內蒙古包頭和河北張北3次6級以上強震。2次活躍期均處于太陽黑子相對數上升階段及日長數據高值區域(見圖7)。

圖7 太陽黑子、日長數據與華北地區中強地震M-t圖Fig.7 Sunspot, day length data and M-t map of moderate and strong earthquakes in North China

對地球日長日值數據進行Demy小波分解的結果顯示，d3部分存在顯著的約18.6 a的周期變化，其小波分解曲線波谷及其附近2～3 a時間段內華北地區中強震頻次明顯升高，發生該活躍幕內的最強地震(見第4頁圖8)。

地球日長數據中包含有18.6 a周期的變化成分，太陽黑子相對數有約11 a的周期變化規律。兩者的共同作用是影響華北地區中強地震活動的重要因素(見圖7、圖8)，當太陽黑子相對數周期上升段與地球日長18.6 a周期的波谷段疊加時，對華北地震活動的影響達到最大，可能觸發區域中強地震的發生。目前，太陽活動處于黑子相對數的第25周期初始上升階段，地球日長18.6 a周期曲線也處于上升階段，兩者不構成疊加關系，因此誘發華北地區進入中強震多發狀態的可能性不大。

圖8 地球日長數據小波分解與華北地區中強地震M-t圖Fig.8 Wavelet decomposition of earth's daily length data and M-t map of moderate and strong earthquakes in North China

3 結論與討論

(1) 根據Morlet小波分解及小波系數分析結果，顯示太陽黑子相對數時間序列變化的主周期是11 a。主周期小波系數實部曲線的變化規律表明第25太陽活動周太陽黑子相對數將處于相對高水平。

(2) 太陽黑子相對數變化、地球自轉速率變化與華北地區中強地震活動水平有比較明顯的相關關系。太陽黑子相對數上升，地球自轉變慢、日長增大時，華北地區中強地震相對活躍。具體表現為當太陽黑子相對數周期上升段與地球日長18.6 a周期波谷段疊加時，對華北地震活動的影響達到最大，可能觸發區域中強地震的發生。當前處于第25太陽活動周初期，太陽黑子相對數為轉折上升狀態，地球日長小波分解三階細節曲線不處于波谷位置，不構成觸發華北地區中強地震的條件。

(3) 不同地區地殼與地幔的耦合關系不同，自轉速度的變化對不同區域地殼活動的影響也不同。因此，在不同地區自轉速度變化對地震活動的影響需要具體分析。