馬青波 解朝娣 向元 王瓊 楊思夢 王英楠

摘要：選取新疆西昆侖地區(qū)的2014年MS7.3和2015年皮山MS6.5地震及其余震序列為研究對象，對地震目錄的完備性進(jìn)行了分析，并計算了主震發(fā)震時刻，主震震源斷層面上的潮汐剪應(yīng)力、正應(yīng)力，得到主震震源斷層面上的潮汐庫侖破裂應(yīng)力；再使用Schuster統(tǒng)計檢驗方法，分析了2個余震序列的潮汐相位角分布，研究新疆西昆侖地區(qū)中強(qiáng)地震及其余震活動與潮汐應(yīng)力之間的觸發(fā)關(guān)系。結(jié)果表明：主震發(fā)生時，皮山震源斷層面上的潮汐庫侖破裂應(yīng)力為正值，而于田震源斷層面上的潮汐庫侖破裂應(yīng)力為負(fù)值，說明潮汐應(yīng)力對皮山主震有明顯的觸發(fā)作用，對于田主震的觸發(fā)作用不明顯。同時，計算結(jié)果表明2個余震序列與潮汐應(yīng)力的統(tǒng)計相關(guān)參數(shù)P值幾乎都小于5%，說明潮汐應(yīng)力對皮山和于田余震序列均有觸發(fā)作用。最后ETAS模型計算結(jié)果也表明了余震序列對地震震級的弱依賴性，可能與外界觸發(fā)因素（潮汐）有關(guān)。

關(guān)鍵詞：2014年于田MS7.3地震；2015年皮山MS6.5地震；Schuster統(tǒng)計檢驗；潮汐應(yīng)力；地震觸發(fā)

中圖分類號：P315.727 文獻(xiàn)標(biāo)識碼：A 文章編號：1000-0666（2018）01-0104-07

0 引言

受印度板塊與歐亞板塊的碰撞作用，新疆塔里木盆地西南邊緣的西昆侖地區(qū)構(gòu)造運(yùn)動十分強(qiáng)烈（潘家偉，2011）。位于該地區(qū)的西昆侖地震帶地震活動頻繁，1985-2011年共發(fā)生MS≥4.0地震近500次，其中4.0≤MS≤5.0地震425次，MS>5.0地震70次。該地區(qū)的地震活動及其發(fā)震機(jī)理研究已取得了不少的進(jìn)展（朱傳慶等，2008；萬永革等，2010；程佳等，2014；姜迪迪等，2014）。吳小平等（2009）結(jié)合構(gòu)造應(yīng)力特點(diǎn)與潮汐應(yīng)力觸發(fā)問題，研究了固體潮與該區(qū)地震活動的關(guān)聯(lián)性，得出該區(qū)的潮汐應(yīng)力觸發(fā)效應(yīng)屬于中等水平；李金等（2014）提出新疆地區(qū)地震活動存在一定的優(yōu)勢潮汐相位。

當(dāng)某區(qū)域或斷層上的構(gòu)造應(yīng)力積累接近臨界狀態(tài)時，若有恰當(dāng)?shù)某毕珣?yīng)力與該構(gòu)造應(yīng)力疊加，總的應(yīng)力疊加效果超過了觸發(fā)閾值，就有可能觸發(fā)地震。已有很多學(xué)者對潮汐觸發(fā)問題進(jìn)行了研究。丁中一等（1983）選取了1957-1983年中國或中國鄰近地區(qū)發(fā)生的70次較大的地震，研究地震與潮汐應(yīng)力的關(guān)系，發(fā)現(xiàn)有61%的地震在一定程度上受到了潮汐應(yīng)力的觸發(fā)。孫長青等（2014）把青藏高原東部及其相鄰地區(qū)劃分成5個構(gòu)造應(yīng)力場分區(qū)，研究潮汐應(yīng)力觸發(fā)效應(yīng)，結(jié)果表明該區(qū)域潮汐觸發(fā)效應(yīng)的強(qiáng)弱與震級大小、震源深度、構(gòu)造應(yīng)力場等有關(guān)。張晶等（2007）選取了1970-2007年多個7級以上強(qiáng)震序列，分析了潮汐應(yīng)力的觸發(fā)效應(yīng)，發(fā)現(xiàn)共有10個強(qiáng)震序列與引潮力相關(guān)性比較好，且在某些特殊區(qū)域其水平分量的方位角在多個方向上存在優(yōu)勢。潮汐觸發(fā)效應(yīng)與地震發(fā)展機(jī)理研究仍是目前地震研究的熱點(diǎn)問題。

本文選取發(fā)生在新疆西昆侖地區(qū)的2014年于田MS7.3、2015年皮山MS6.5地震及其余震序列為研究對象，首先進(jìn)行地震目錄完備性分析，以確定最小完備震級MC；然后計算主震震源斷層面上的潮汐庫侖破裂應(yīng)力值，以此分析2次主震與潮汐應(yīng)力之間的觸發(fā)關(guān)系。在此基礎(chǔ)上，確定序列中每個地震的潮汐相位角，進(jìn)而運(yùn)用Schuster統(tǒng)計檢驗方法分析地震序列與潮汐應(yīng)力的統(tǒng)計相關(guān)性，探討潮汐應(yīng)力對于田和皮山這2個余震序列的觸發(fā)效應(yīng)，進(jìn)一步研究新疆地震活動發(fā)震機(jī)理。

1 研究區(qū)構(gòu)造背景

新疆是地震多發(fā)區(qū)，地震主要集中在昆侖山、塔里木盆地等構(gòu)造活動活躍的地區(qū)，且呈現(xiàn)條帶狀分布（陳建波，2008）。新生代以來，西昆侖地區(qū)由于受到歐亞板塊與印度板塊相互作用所造成的擠壓應(yīng)力，該地區(qū)的構(gòu)造線呈現(xiàn)出北西一南東或接近東西方向（潘家偉，2011）。由于這種擠壓作用的存在，產(chǎn)生了康西瓦、普魯、西昆侖山前等多條較大規(guī)模的斷裂。鄒長橋等（2012）認(rèn)為，位于西昆侖山南邊的西昆侖地震活動區(qū)，由于地殼與巖石圈上地幔的解耦使得區(qū)域構(gòu)造背景更為復(fù)雜。姜迪迪等（2014）在研究中國西部地區(qū)的地殼結(jié)構(gòu)特征與強(qiáng)震活動的關(guān)系中，發(fā)現(xiàn)西昆侖地區(qū)地震的震源深度達(dá)到了莫霍面深度。該地區(qū)特殊的構(gòu)造背景導(dǎo)致大震、強(qiáng)震較多。2014年2月12日，于田縣發(fā)生MS7.3地震，震中位于阿什庫勒斷裂北段（程佳等，2014）；2015年7月3日，皮山縣發(fā)生MS6.5地震，主震位于西昆侖山前斷裂北緣。圖1給出了2次地震序列的空間分布、主震的震源機(jī)制解、研究區(qū)域內(nèi)的主要斷裂分布等情況，其中高程數(shù)據(jù)來源于USGS公布的全球數(shù)字高程模型GTOP030，其它數(shù)據(jù)的來源見表1。

2 地震數(shù)據(jù)選取

基于西昆侖地區(qū)復(fù)雜的構(gòu)造背景和地震空間分布特點(diǎn)，本文采用新疆地震局提供的2014年于田MS7.3和2015年皮山MS6.5地震序列精定位后的震源參數(shù)進(jìn)行分析。這2個序列的主震震級均在6.0級以上，余震序列的個數(shù)均在2000個以上，適合統(tǒng)計分析。表1給出2個地震序列的主要參數(shù)和數(shù)據(jù)來源。

為使潮汐相關(guān)性的統(tǒng)計分析結(jié)果更準(zhǔn)確，本文首先對所選地震序列目錄的完備性進(jìn)行分析，結(jié)果如圖2所示。分析得到皮山地震序列的MC=2.0，序列中M≥2.0地震共有1088個，b值約為0.53；于田地震序列的MC=2.0，M≥2.0地震共1297個，b值約為0.50。

3 研究原理和方法

3.1 潮汐應(yīng)力計算方法

假設(shè)地球內(nèi)部的物理性質(zhì)具有球?qū)ΨQ性，在前人研究的基礎(chǔ)上（蔣駿，張雁濱，1995；Agnew，2007；李智蓉，張曉東，2011；熊子瑤等，2015），地球內(nèi)部或者表面某一點(diǎn)P在球坐標(biāo)下由潮汐作用下產(chǎn)生的應(yīng)變分量eθθ、eλλ、eθλ為：式中：Vn為n階起潮力位；g為地球的平均重力值；φ為P點(diǎn)的地理緯度；r為P點(diǎn)到地心的距離；Hn（r）和Ln（r）為Love數(shù)。表述應(yīng)力與應(yīng)變關(guān)系的胡克定律為：

σij=2μeij+λθδij（2）其中：λ，μ為拉梅常量；θ為體應(yīng)變；σij為應(yīng)力分量；eij為應(yīng)變分量。

將（1）式代入（2）式后，便可將P點(diǎn)的潮汐應(yīng)變分量由胡克定律轉(zhuǎn)換為P點(diǎn)的應(yīng)力張量。在此基礎(chǔ)上，根據(jù)主震的震源斷層面參數(shù)，將潮汐應(yīng)力張量通過矩陣旋轉(zhuǎn)投影到震源斷層面上，得到潮汐剪應(yīng)力τ、潮汐正應(yīng)力σn；然后根據(jù)庫侖破裂準(zhǔn)則，計算得出主震斷層面上的潮汐庫侖破裂應(yīng)力：

Tcfs=τ+μ'σn（3）其中：Tcfs為斷層面上因太陽和月亮引力影響所產(chǎn)生的庫侖破裂應(yīng)力，本文中μ'取值0.5，代表斷層的有效摩擦系數(shù)（Harris，1998；解朝娣等，2015）。

3.2 統(tǒng)計檢驗方法

傳統(tǒng)的潮汐相位分析方法是計算諸如Mf、M2等特定周期的潮汐波（解朝娣等，2015）。Schuster（1897）將Rayleigh在1880年提出的時間序列統(tǒng)計檢驗原理應(yīng)用到地震活動與日月相關(guān)性的統(tǒng)計檢驗當(dāng)中。Schuster統(tǒng)計檢驗（Tanaka et al，2006，2010）建立了潮汐現(xiàn)象隨時間變化的函數(shù)，并且這個函數(shù)代表了作用在地震斷層上的潮汐波的總和。該檢驗方法通過計算地震事件所對應(yīng)的潮汐相位角，統(tǒng)計整個地震目錄中地震事件相位角的分布情況，分析某一時間段或某一地區(qū)是否存在優(yōu)勢相位角，進(jìn)而判斷研究區(qū)域內(nèi)的地震活動是否受到了潮汐的觸發(fā)（Tanaka et al，2006，2010；孫長青等，2014）。具體而言，首先將發(fā)震時間轉(zhuǎn)換為國際時，然后采用Gtoic2（Matsumoto et al，2001）按照前文所述的方法計算潮汐應(yīng)變，并確定該地震事件所對應(yīng)的潮汐相位角，再將地震序列中每個地震事件對應(yīng)的潮汐相位角一一標(biāo)出，最后統(tǒng)計分析潮汐相位角是否在特定角度附近集中。

在Schuster檢驗中，規(guī)定第j個地震事件為單位向量ej（方向由該事件對應(yīng)的相位角確定，向量的模為1）。然后將N個地震事件的向量相加得到向量D，表示為：式中：N為地震事件的總數(shù)；θj為第j個地震事件的潮汐相位角。

如果N個地震事件是隨機(jī)發(fā)生的，其潮汐相位角也是隨機(jī)分布的。當(dāng)N足夠大（N>10）時，和向量的長度大于或等于|D|的概率（Tsuruoka etal，1995）為：式中：D2即為向量D模的平方；P值為定量描述潮汐效應(yīng)觸發(fā)地震活動的指標(biāo)，P值范圍為0～1，P值越小，潮汐對地震的觸發(fā)效應(yīng)越明顯（李金等，2014；解朝娣等，2015）。前人通常取P<5%作為判斷潮汐能否觸發(fā)地震的標(biāo)準(zhǔn)（Tanaka，2010；解朝娣等，2015）。

3.3 ETAS模型分析方法

由大森公式發(fā)展得到的ETAS模型（Ogata，2005，2011：Lei el al，2008，2011）可以分析序列中活動性的異常變化，這對于研究地震活動性和外部觸發(fā)因素有幫助。ETAS模型假設(shè)每個地震都有觸發(fā)余震的能力，此能力具有震級依賴性（Leiet al，2013）；并給出總的地震發(fā)生率λ（t）可以由先前地震觸發(fā)的次級余震發(fā)生率和表示“背景”地震活動的發(fā)生率λ0（t）之和來表述：式中：Mi和MC分別表述第i個地震事件的震級和地震序列的最小完備震級。一般情況下，“背景”地震活動的發(fā)生率λ0為常數(shù)，不是時間t的函數(shù)；但如果外部觸發(fā)因素存在時，λ0（t）就包括構(gòu)造背景活動和外部觸發(fā)地震活動，分析λ0（t）就有助于分析外部觸發(fā)地震活動。K0，c，α和P這些模型參數(shù)都是通過最小AIC準(zhǔn)則檢驗給出，合適的平滑化參數(shù)ne也可以通過最小AIC準(zhǔn)則來確定。

4 結(jié)果分析

4.1 主震斷層面上的潮汐庫侖破裂應(yīng)力計算結(jié)果

采用上述潮汐應(yīng)力計算方法，首先計算球坐標(biāo)系中于田MS7.3和皮山MS6.5主震的理論潮汐應(yīng)變張量，通過胡克定律轉(zhuǎn)換為應(yīng)力張量，然后結(jié)合表1中列出的主震震源參數(shù)，計算出震源斷層面上前后一天的潮汐剪應(yīng)力τ、潮汐正應(yīng)力σn和潮汐庫侖破裂應(yīng)力Tcfs隨時間變化的關(guān)系。根據(jù)相關(guān)研究（房立華等，2015），2014年于田MS7.3主震的震源斷層面為北東向的節(jié)面，因此采用表1中節(jié)面1的參數(shù)進(jìn)行計算；同理，根據(jù)相關(guān)研究（張廣偉等，2016），2015年皮山MS6.5主震的震源斷層面采用表1中節(jié)面1的參數(shù)進(jìn)行計算。通過計算可以得到于田MS7.3和皮山MS6.5主震發(fā)震時，震源斷層面上的潮汐應(yīng)力，結(jié)果如圖3和表2所示。

結(jié)合表2和圖3可見，在于田主震發(fā)震時刻，于田主震震源斷層面上的潮汐剪應(yīng)力處于波谷值與波峰值之間的上升階段，為-380.78Pa，而潮汐正應(yīng)力和潮汐庫侖破裂應(yīng)力位于波峰值到波谷值的下降階段，分別為-1429.98Pa、-1095.77Pa。在皮山主震發(fā)震時刻，皮山主震震源斷層面上的潮汐剪應(yīng)力和潮汐庫侖破裂應(yīng)力處于波峰值到波谷值的上升階段，分別為1111.25Pa、709.75Pa，潮汐正應(yīng)力處于波峰值到波谷值的下降階段，為-803.00Pa。

4.2 余震序列與潮汐應(yīng)力相關(guān)性統(tǒng)計分析結(jié)果

使用表1中列出的震源參數(shù)，對2個地震序列中震級大于Mc的所有余震事件進(jìn)行潮汐相位分析，確定各事件的潮汐相位角，計算得到潮汐相位角的頻率分布圖，于田和皮山主震的震源斷層面均采用表1中節(jié)面1進(jìn)行分析，結(jié)果如圖4所示。圖中顯示，統(tǒng)計計算得到的P值幾乎都小于5%。再結(jié)合圖3和表2中潮汐庫侖破裂應(yīng)力分析結(jié)果，可以發(fā)現(xiàn)潮汐應(yīng)力對皮山和于田地震的余震序列均有觸發(fā)作用，但是潮汐應(yīng)力對于田主震的觸發(fā)效果不明顯，因為于田主震斷層面上受到的潮汐庫侖破裂應(yīng)力為負(fù)值。

4.3 ETAS模型計算結(jié)果

從圖5可以看出，通過最小AIC準(zhǔn)則得到的2個序列的ne均為14?！氨尘啊钡卣鸹顒勇剩╢.s.）和模型的其他參數(shù)如圖5所示。結(jié)果表明2個余震序列以受迫成分為主（f.s.>80%），說明于田、皮山余震序列中80%以上的事件是來自于外部因素的影響，諸如潮汐應(yīng)力和其他影響因素，其余的是大森型自觸發(fā)地震。于田和皮山地震序列的模型參數(shù)α分別為1.069和0.918，小的α值表明了序列對地震震級的弱依賴性，可以從背景地震活動中區(qū)分開（Ogata，1999；Lei et al，2008）。

5 結(jié)論與討論

本文選取新疆西昆侖地區(qū)的2014年于田MS7.3、2015年皮山MS6.5兩次中強(qiáng)地震及其余震序列進(jìn)行地震潮汐觸發(fā)研究。

（1）首先通過計算主震震源斷層面上的潮汐剪應(yīng)力、潮汐正應(yīng)力，得到了主震震源斷層面上的潮汐庫侖破裂應(yīng)力。結(jié)果顯示：于田主震斷層面上的潮汐庫侖破裂應(yīng)力值為負(fù)值，而皮山主震斷層面上的潮汐庫侖破裂應(yīng)力值為正值。依據(jù)地震觸發(fā)的相關(guān)研究結(jié)果（Harris，1998），正的庫侖破裂應(yīng)力變化超過觸發(fā)閾值時，可能起到地震觸發(fā)作用；而負(fù)的庫侖破裂應(yīng)力變化沒有觸發(fā)作用，甚至可能起到抑制后續(xù)地震發(fā)生的作用。但潮汐應(yīng)力的地震觸發(fā)作用除了與潮汐庫侖破裂應(yīng)力值的正負(fù)以及大小（是否超過觸發(fā)閾值）有關(guān)，還與潮汐應(yīng)力的時間震蕩性有關(guān)（解朝娣等，2015）；潮汐觸發(fā)作用的研究不能單一的以潮汐庫侖破裂應(yīng)力的正負(fù)作為唯一判斷標(biāo)準(zhǔn)。因此，本文結(jié)果表明潮汐應(yīng)力對皮山主震有明顯的觸發(fā)作用，而對于田主震的觸發(fā)作用不明顯。

（2）通過Schuster統(tǒng)計檢驗方法，計算2個余震序列潮汐庫侖破裂應(yīng)力的相位角分布發(fā)現(xiàn)，潮汐應(yīng)力對皮山地震和于田地震的余震序列均有觸發(fā)作用。最后結(jié)合ETAS模型分析，也發(fā)現(xiàn)西昆侖地區(qū)的地震活動受到來自于外部因素的觸發(fā)作用，諸如潮汐應(yīng)力和其他影響因素，與大森型自觸發(fā)的余震活動不一致。

新疆西昆侖地區(qū)由于受到歐亞板塊與印度板塊相互作用所造成的擠壓應(yīng)力，產(chǎn)生了多條較大規(guī)模的斷裂，以及一系列的活動斷層，而斷層交匯處，是流體活動的通道，所以該地區(qū)的地?zé)崃黧w活動很豐富的。而流體對潮汐的拉伸作用敏感，可以凸顯潮汐觸發(fā)作用。這可能是該地區(qū)潮汐觸發(fā)作用較為明顯的主要原因。

參考文獻(xiàn)：

陳建波.2008.新疆地震構(gòu)造特征研究[D].蘭州：中國地震局蘭州地震研究所.

程佳，楊文，劉杰，等.2014.2014年2月12日新疆于田MS7.3地震序列及其構(gòu)造背景研究[J].地震學(xué)報，36（3）：350-361.

丁中一，賈晉康，王仁.1983.潮汐應(yīng)力對地震的觸發(fā)作用[J].地震學(xué)報，5（2）：46-58.

房立華，吳建平，王未來，等.200.2014年新疆于田MS7.3級地震序列重定位[J].地球物理學(xué)報，58（3）：802-808.

姜迪迪，江為為，胥頤，等.2014.中國西部地區(qū)地殼結(jié)構(gòu)特征與強(qiáng)震活動相關(guān)性研究[J].地球物理學(xué)報，57（12）：4029-4040.

蔣駿，張雁濱.1995.發(fā)震斷層面上潮汐有效剪切應(yīng)力增量的計算與加卸載的確定[J].中國地震，（1）：72-83.

解朝娣，Lei X L，吳小平，等.2015.潮汐應(yīng)力對2007年雞6.4寧洱地震震源斷層成核失穩(wěn)過程的影響[J].中國科學(xué)：地球科學(xué)，45（9）：1409-1420.

李金，蔣海昆，桂榮，等.2014.新疆喀什—烏恰交匯區(qū)潮汐觸發(fā)地震活動的統(tǒng)計檢驗[J].中國地震，30（1）：64-73.

李金，王瓊.2015.2014年2月12日于田7.3級地震序列震源機(jī)制特征分析[J].中國地震，31（1）：110-120.

李智蓉，張曉東.2011.潮汐應(yīng)力與大震關(guān)系研究[J].地震，31（4）：48-57.

潘家偉.2011.西昆侖構(gòu)造地貌與阿什庫勒地區(qū)活動構(gòu)造研究[D].北京：中國地質(zhì)科學(xué)院.

孫長青，閻春恒，吳小平，等.2014.青藏高原東部及鄰區(qū)地震斷層面上的潮汐應(yīng)力觸發(fā)效應(yīng)[J].地球物理學(xué)報，57（7）：2054-2064.

萬永革，沈正康，盛書中，等.2010.2008年新疆于田7.3級地震對周圍斷層的影響及其正斷層機(jī)制的區(qū)域構(gòu)造解釋[J].地球物理學(xué)報，53（2）：280-289.

吳小平，冒蔚，黃雍，等.2009.基于不同構(gòu)造分區(qū)中國地震的潮汐應(yīng)力觸發(fā)效應(yīng)及相關(guān)天文特征[J].中國科學(xué)：物理學(xué)力學(xué)天文學(xué)，39（6）：901-912.

熊子瑤，解朝娣，雷興林，等.2015.云南大姚、姚安及附近地區(qū)地震活動與潮汐應(yīng)力的觸發(fā)關(guān)系[J]地震研究，38（1）：66-74.

張廣偉，張洪艷，孫長青.2016.2015年新疆皮山MS6.5地震震源機(jī)制及余震序列定位[J].地震地質(zhì)，38（3）：711-720.

張晶，郗欽文，楊林章，等.2007.引潮力與潮汐應(yīng)力對強(qiáng)震觸發(fā)的研究[J].地球物理學(xué)報，50（2）；448-454.

朱傳慶，楊書江，李同彬，等.2008.塔里木盆地主要斷裂的地球物理特征及與天然地震的關(guān)系[J].地球物理學(xué)進(jìn)展，23（6）：1685-1691.

鄒長橋，賀日政，張智.2012.藏北高原地震活動性特征及其大地構(gòu)造意義[J].地球物理學(xué)進(jìn)展，27（2）；429-440.

Agnew D C.2007.Earth Tides，in Treatise on Geophysics：Geodesy，T AHerring，ed[M].New York：Elsevier，163-195.

Harris R A.1998.Introduction to Special Section；Stress Triggers，StressShadows，and Implications for Seismic Hazard[J].J Geophys Res，103（B10）：24347-24358.

Lei X，Shengli M A，Jinrong S U，et al.2013.Inelastic triggering of the2013 MW6.6 Lushan earthquake by the 2008 MW7.9 Wenchuanearthquake[J].Seismology And Geology，35（2）：411-422.

Lei X，Xie C，F(xiàn)u B.2011.Remotely triggered seismicity in Yunnan，south-western China，following the 2004 MW9.3 Sumatra earthquake[J].JGeophys Res，116（B8）：114-123.

Lei X，Yu G，Ma S，et al.2008.Earthquakes induced by water injection at～3 km depth within the Rongchang gas field，Chongqing，China[J].J Geophys Res，113（B10：1-12），doi：10.1029/2008JB005604.

Matsumoto K，Takanezawa T，Ooe M.2001.GOTIC2：A program for com-putation of oceanic tidal loading effect[J].J Geod Soc Japan，47（1），243-248.

Ogata Y.1999.Seismicity Analysis through Point-process Modeling：AReview[J].Pure and Applied Geophysics，155（2）：471-507.

Ogata Y.2005.Detection of anomalous seismicity as a stress change sen-sor[J].J Geophys Res，110（B5）：767-782.

Ogata Y.2011.Pre-seismic anomalies in seismicity and crustal deforma-tion：case studies of the 2007 Noto Hanto earthquake of M6.9 andthe 2007 Chuetsu-oki earthquake of M6.8 after the 2004 Chuetsuearthquake of M6.8[J].Geophys J Int，186（1）：331-348.

Schuster A.1897.On Lunar and Solar Periodicities of Earthquakes[J].Nature，61（1449）：321-321.

Tanaka S，Sato H，Matsumura S，et al.2006.Tidal triggering of earth-quakes in the subducting Philippine Sea plate beneath the lockedzone of the plate interface in the Tokai region，Japan[J].Tectono-physics，417（1-2）：69-80.

Tanaka S.2010.Tidal triggering of earthquakes precursory to the recentSumatra megathrust earthquakes of 26 December 2004（MW9.0），28March 2005（MW8.6），and 12 September 2007（MW8.5）[J].Geo-phys Res Lett，37（37）：195-205.

Tsuruoka H，Ohtake M，Sato H.1995.Statistical test of the tidal triggeringof earthquakes：contribution of the ocean tide loading effect[J].Geo-phys J Int，122（1）：183-194.