• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    帕米爾—興都庫(kù)什深俯沖殘留體對(duì)410km間斷面起伏形態(tài)的影響

    2015-12-12 08:50:14眭怡周元澤王曉冉
    地球物理學(xué)報(bào) 2015年1期
    關(guān)鍵詞:臺(tái)網(wǎng)震源波形

    眭怡,周元澤*,王曉冉

    1 中國(guó)科學(xué)院計(jì)算地球動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室,北京 100049

    2 中國(guó)科學(xué)院大學(xué)地球科學(xué)學(xué)院,北京 100049

    1 引言

    地震波可以直接帶來(lái)地球內(nèi)部結(jié)構(gòu)的信息,人們利用強(qiáng)體波震相、面波乃至地球自由振蕩等的信息,獲得了一系列可靠的一維徑向分層地球模型,如PREM (Dziewonski and Anderson,1981),IASP91(Kennett and Engdahl,1991),AK135(Kennett et al.,1995).410km間斷面作為地幔轉(zhuǎn)換區(qū)的上界面出現(xiàn)在這三個(gè)速度模型中,其精細(xì)結(jié)構(gòu)對(duì)于認(rèn)識(shí)地球內(nèi)部物質(zhì)組成及地球內(nèi)部運(yùn)移的地球動(dòng)力學(xué)過(guò)程具有重要意義(Deuss,2009).

    人們利用長(zhǎng)周期PP和SS前驅(qū)波(Flanagan and Shearer,1998a,1998b;Deuss,2009)、ScS多次反射波(Revenaugh and Jordan,1991)等,獲得了410間斷面的全球性存在及其大尺度上的起伏形態(tài).結(jié)合高溫高壓礦物物理實(shí)驗(yàn)和理論分析,410km間斷面被認(rèn)為是α相橄欖石到β相橄欖石(也稱(chēng)瓦茲利石)的放熱相變界面(如,Bina and Helffrich,1994;Collier et al.,2001;金振民等,2013).在熱地幔柱地區(qū)410km間斷面相變線表現(xiàn)為下降;而在冷的俯沖板片及其殘留物存在的情況下,410km間斷面相變線表現(xiàn)為上升(Bina and Helffrich,1994;Song et al.,2004;Obayashi et al.2006;Deuss,2009).

    隨著世界上受地震災(zāi)害影響比較嚴(yán)重且經(jīng)濟(jì)發(fā)達(dá)或相對(duì)發(fā)達(dá)的國(guó)家對(duì)于測(cè)震工作的重視,美國(guó)、日本、歐洲部分地區(qū)和中國(guó)均開(kāi)展了高密度、高質(zhì)量固定或流動(dòng)地震臺(tái)網(wǎng)布設(shè)工作,以加強(qiáng)地震活動(dòng)監(jiān)測(cè).近些年來(lái),大量高密度臺(tái)網(wǎng)的高質(zhì)量觀測(cè)波形資料不斷積累,為地球內(nèi)部結(jié)構(gòu)研究提供了觀測(cè)數(shù)據(jù)支撐.基于接收函數(shù)方法提取到相對(duì)低頻的PdS轉(zhuǎn)換震相(d為轉(zhuǎn)換深度),人們對(duì)于區(qū)域尺度上的410km間斷面的存在形態(tài)有了更多的認(rèn)識(shí)(如Zhang et al.,2010;Tian et al.,2010,2011),但相比于受冷或熱物質(zhì)影響的局部大起伏而言,這些結(jié)果顯示的起伏則要小得多(如Li et al.,2000).進(jìn)一步利用相對(duì)短周期波形資料來(lái)研究受冷或熱的局部溫度異常影響下的410km間斷面小尺度上的結(jié)構(gòu)形態(tài)就顯得很重要.

    基于臺(tái)網(wǎng)/臺(tái)陣資料疊加方法提取離源弱次生震相是研究地幔間斷面小尺度起伏形態(tài)的主要方法之一(Rost and Thomas,2002).弱次生轉(zhuǎn)換震相SdP是以S波離源到達(dá)深度為d的間斷面之后轉(zhuǎn)換為P波,并為遠(yuǎn)處地震臺(tái)站所記錄.一般而言,考慮到減少橫向不均勻的影響,常用30~90°震中距范圍內(nèi)的地震記錄來(lái)提取SdP(臧紹先和周元澤,2002;Rost and Thomas,2002).中、深源地震的波形記錄中的SdP因震源深度與附近間斷面距離較小,因而相應(yīng)的地震射線具有較小的Fresnel區(qū),這有助于從觀測(cè)波形的疊加圖像中識(shí)別到該震相,從而能夠確認(rèn)地幔間斷面的小尺度上的存在及起伏形態(tài).

    西太平洋俯沖區(qū)擁有全球最豐富的中、深源地震資源,部分地區(qū)密集地震臺(tái)網(wǎng)可以在合適的震中范圍內(nèi)記錄這些地震.很多研究者基于這些資料,利用傾斜疊加方法來(lái)提取轉(zhuǎn)換震相SdP,很好地確認(rèn)了西太平洋地區(qū)(主要是伊豆—小笠原和湯加—斐濟(jì)地區(qū))冷俯沖物質(zhì)影響下的410km間斷面的相變線深度的變化及復(fù)雜結(jié)構(gòu)(如Collier et al.,2001).陸陸碰撞區(qū)下方的俯沖殘留物質(zhì)較之持續(xù)俯沖過(guò)程中的海洋巖石圈俯沖帶來(lái)滯留物質(zhì)要少(Fukao and Obayashi,2013),溫度異常也會(huì)弱,因而,410km間斷面相變深度所受到的影響也會(huì)小.

    帕米爾—興都庫(kù)什地區(qū)位于印度板塊與歐亞板塊碰撞的西構(gòu)造結(jié)是世界上少有的陸陸碰撞中、深源地震區(qū)(如寧杰遠(yuǎn)和臧紹先,1990;孫文斌等,2009;Negredo et al.,2007;Sippl et al.,2013).基于震源機(jī)制和高精度地震重定位結(jié)果(Pegler and Das,1998;Sippl et al.,2013;Bai and Zhang,2015)表明,興都庫(kù)什地區(qū)板片向北俯沖是重力引起的,而帕米爾地區(qū)板塊向南俯沖則是由上地幔流所決定的,這樣的雙向俯沖(如圖1所示)構(gòu)成了喜馬拉雅西構(gòu)造結(jié)區(qū)域復(fù)雜的高原地貌(張浪平等,2014).由于帕米爾—興都庫(kù)什地區(qū)海洋板塊俯沖進(jìn)入了上地幔深部,直至地幔轉(zhuǎn)換區(qū)處(Negredo,et al.,2007),而大陸巖石圈俯沖則緊隨其后,因而在本區(qū)域具有豐富的中-深源地震最深可以到達(dá)380km左右(Pegler and Das,1998).選擇本區(qū)域研究410km間斷面受俯沖滯留物質(zhì)影響下的起伏形態(tài),可以很好地與海洋巖石圈俯沖板塊對(duì)該間斷面的影響進(jìn)行對(duì)比分析.

    本文基于高密度的日本 Hi-net地震臺(tái)網(wǎng)和Earth Scope計(jì)劃組成部分的USArray流動(dòng)臺(tái)陣(TA:Transportable Array)所記錄的地震波形數(shù)據(jù),利用N次根傾斜疊加方法,提取P尾波中的SdP次生轉(zhuǎn)換震相,分析了帕米爾—興都庫(kù)什地區(qū)俯沖板塊物質(zhì)影響下的410km間斷面的存在及形態(tài).

    2 數(shù)據(jù)處理

    由于震源過(guò)于復(fù)雜或者來(lái)自地表反射的強(qiáng)pP和sP震相及多次波震相會(huì)對(duì)弱次生SdP震相的拾取帶來(lái)很大的干擾,因而我們?cè)谑占卣鸩ㄐ钨Y料時(shí)需要限制地震事件的震級(jí),以保證所使用的波形資料具有相對(duì)簡(jiǎn)單的震源時(shí)間函數(shù),從而有助于提取次生震相.本文中選用的地震事件震級(jí)范圍為Mb5.6~6.4,震源時(shí)間函數(shù)時(shí)長(zhǎng)不超過(guò)4s,且具有足夠高的信噪比.雖然本區(qū)域有較多的中、深源地震,但受限于臺(tái)網(wǎng)記錄時(shí)間以及波形資料選擇的要求,我們收集了2004年以來(lái)的6個(gè)地震的Hi-net臺(tái)網(wǎng)(圖2a)短周期波形資料;同時(shí)還收集了其中一個(gè)地震的TA臺(tái)陣(圖2b)寬頻帶波形資料.地震分布如圖1所示,相關(guān)地震參數(shù)見(jiàn)表1.本文中所用到的地震波形資料來(lái)自于日本防災(zāi)科學(xué)技術(shù)研究所(http://www.hinet.bosai.go.jp)和美國(guó)地震學(xué)聯(lián)合會(huì)(IRIS,http://www.iris.edu).

    表1 本文使用的地震事件相關(guān)參數(shù)Table 1 Paramets of the events used

    圖1 地震分布圖以沙灘球表示的矩張量解來(lái)自http://www.globalcmt.org,沒(méi)有矩張量解的地震事件以大圓點(diǎn)表示,等值線給出的是Wadati-Benioff帶在地幔中存在的深度分布 (Gudmundsson and Sambridge,1998),等值線上的白色數(shù)字為深度,單位為km.Fig.1 The distribution of the events usedThe focal mechanisms labeled with beach balls are the CMT solutions from http://www.globalcmt.org and the solid circle for another one without CMT solution.The contours are for the extension of the Wadati-Benioff slab(Gudmundsson and Sambridge,1998).The white numbers are the depths of the contous with units in km.

    震源深度的確定直接影響到轉(zhuǎn)換點(diǎn)深度反演的精度,因而準(zhǔn)確的震源深度是非常重要的.考慮到本文所使用的局部地區(qū)地震臺(tái)網(wǎng)來(lái)進(jìn)行震源深度確定并不能獲得很好的深度精度,因而本文選用了目前震源深度確定最好的EHB數(shù)據(jù)庫(kù)(http://www.isc.ac.uk/EHB/index.html)給出的2008年前(含2008年)的地震位置及深度參數(shù)(Engdahl et al.,1998),而2009年起的地震參數(shù)則來(lái)自于國(guó)際地震中心(http://www.isc.ac.uk)提供的地震位置及基于全球波形資料的pP震相確定的震源深度.

    文中我們采用0.2~1.0Hz的濾波頻段對(duì)地震波形進(jìn)行濾波.這一濾波頻段可以很好地拾取出與地幔間斷面相關(guān)的SdP轉(zhuǎn)換震相(Castle and Creager,2000),同時(shí)壓制噪聲.在處理中,我們進(jìn)一步人工刪除受脈沖噪聲等影響的非正常波形記錄.

    地幔間斷面比較弱,來(lái)自這樣速度界面的轉(zhuǎn)換震相SdP相對(duì)也是很弱的,一般地是無(wú)法直接觀察到的.圖3所示為Ev.5被 Hi-net子臺(tái)網(wǎng) H6(如圖2a)所記錄的波形資料圖.各記錄以直達(dá)P震相的最強(qiáng)峰值進(jìn)行歸一化處理,并將其到達(dá)時(shí)刻歸零.圖中我們以虛線標(biāo)識(shí)了直達(dá)P震相和后面傾斜疊加結(jié)果中識(shí)讀到的S376P次生轉(zhuǎn)換震相.該次生震相直接從波形圖上則很難觀測(cè)到,甚至應(yīng)該較強(qiáng)的pP震相也很難被清晰地識(shí)讀出.

    為了有效地獲取SdP次生轉(zhuǎn)換震相,我們使用了N次根傾斜疊加方法來(lái)增強(qiáng)次生震相并壓制噪聲(Kanasewich et al.,1973;McFadden et al.,1986;臧紹先和周元澤,2002;Rost and Thomas,2002).由于目標(biāo)震相SdP位于強(qiáng)直達(dá)震相P之后,因而N次根疊加處理中,我們選用直達(dá)P震相作為參考震相.本文所用的地震震源時(shí)間函數(shù)比較簡(jiǎn)單,疊加處理過(guò)程中,我們將直達(dá)P震相的峰值走時(shí)歸為0,并基于峰值對(duì)相應(yīng)的波形資料進(jìn)行幅度歸一化處理(如圖3),因而在疊加圖灰度圖上直達(dá)P震相的走時(shí)差和慢度差均為0.

    圖2 地震臺(tái)網(wǎng)/臺(tái)陣的臺(tái)站分布圖(a)Hi-net臺(tái)網(wǎng)的分布,其中H1—H9為事件Ev.5可用資料對(duì)應(yīng)的各子臺(tái)網(wǎng)臺(tái)站分布圖,圓點(diǎn)為所有的Hi-net臺(tái)站分布;(b)TA的臺(tái)站(三角形)分布圖.Fig.2 Station distributions of the Hi-net seismic network and the TA(a)The station distribution of the Hi-net.The sub-network distributions(H1—H9)of Hi-net which provided the waveform data of the Ev.5are labeled with different color legends,and the small dots are for all others.(b)The station distribution of TA labeled with triangles.

    在具體疊加處理中,TA臺(tái)陣震中范圍相對(duì)較小,且可用臺(tái)站資料比較少,我們對(duì)所有資料直接進(jìn)行4次根疊加處理.日本Hi-net臺(tái)網(wǎng)具有很高的臺(tái)站密度,且質(zhì)量較高.根據(jù)資料的實(shí)際質(zhì)量情況,我們將Hi-net臺(tái)網(wǎng)的信噪比超過(guò)3的波形資料按照震中距和方位角范圍進(jìn)行分組,形成一個(gè)地震為多個(gè)Hi-net子臺(tái)網(wǎng)(如圖2a)所記錄的波形資料組.進(jìn)一步地,我們分別對(duì)各組資料進(jìn)行4次根傾斜疊加處理,獲得各地震-子臺(tái)網(wǎng)的走時(shí)差-慢度差域的相對(duì)幅度圖,從而可以對(duì)相同地震不同臺(tái)網(wǎng)/子臺(tái)網(wǎng)和不同地震相同子臺(tái)網(wǎng)的疊加結(jié)果進(jìn)行對(duì)比,進(jìn)而分析速度界面的起伏,并驗(yàn)證結(jié)果的可靠性.為更好地顯示弱次生震相的存在,進(jìn)一步求取疊加圖的包絡(luò)線并求取20倍的以10為底的對(duì)數(shù)值(dB),據(jù)此繪制出走時(shí)差-慢度差域的灰度圖(Vespegram).參考理論走時(shí)差和慢度差,可以從灰度圖中讀取出慢度差理論值±0.1s/(°)范圍內(nèi)的次生震相SdP;在此基礎(chǔ)上基于IASP91模型反演出相應(yīng)的轉(zhuǎn)換點(diǎn)深度及位置.詳細(xì)的資料分析處理過(guò)程可以參看臧紹先和周元澤(2002).

    3 結(jié)果與討論

    圖3 Hi-net子臺(tái)網(wǎng)H6記錄到的Ev.5地震垂向波形0s時(shí)刻虛線標(biāo)注了歸一化的直達(dá)P震相,標(biāo)有S376P的虛線為對(duì)應(yīng)著來(lái)自376km速度界面的次生轉(zhuǎn)換震相.Fig.3 Waveform data of Ev.5recorded by the sub-network H6of the Hi-netThe dashed line with“P”at 0sis for the peaks of the direct P phases and the dashed line with“S376P”is for the arrival times of the conversion phases from the 410km discontinuities.

    我們以Ev.5為例分析臺(tái)網(wǎng)資料的4次根傾斜疊加結(jié)果.圖4是地震Ev.5為Hi-net各子臺(tái)網(wǎng)記錄波形的4次根傾斜疊加結(jié)果,H6子臺(tái)網(wǎng)記錄的波形資料見(jiàn)圖3.各疊加圖上方的箭頭對(duì)應(yīng)的數(shù)字標(biāo)注出了可能的次生轉(zhuǎn)換震相對(duì)應(yīng)的轉(zhuǎn)換點(diǎn)深度.由于震源深度的限制,較強(qiáng)的直達(dá)震相P和pP之間的次生轉(zhuǎn)換震相均與360~381km深度上速度躍變面有關(guān)系.

    類(lèi)似地,我們將各地震為Hi-net不同子臺(tái)網(wǎng)和TA臺(tái)陣記錄到的波形資料進(jìn)行4次根傾斜疊加處理,共獲得36組能夠清晰地觀察到與410km間斷面相關(guān)SdP次生震相的4次根傾斜疊加灰度圖.相關(guān)的轉(zhuǎn)換點(diǎn)深度見(jiàn)表2.不是每個(gè)子臺(tái)網(wǎng)都能識(shí)讀到具有很好的強(qiáng)度和符合走時(shí)差-慢度差關(guān)系的次生轉(zhuǎn)換震相,因而表中存在著部分子臺(tái)網(wǎng)轉(zhuǎn)換點(diǎn)深度的缺失.

    本文中只有一個(gè)地震(Ev.4)共同為T(mén)A臺(tái)陣和Hi-net臺(tái)網(wǎng)所記錄,并給出了很好的傾斜疊加結(jié)果.TA臺(tái)陣在95°左右的震中距記錄到該地震,而Hi-net則在40°~50°震中距上記錄到該地震.從轉(zhuǎn)換點(diǎn)平均深度上來(lái)看,Hi-net臺(tái)網(wǎng)各子臺(tái)網(wǎng)給出的轉(zhuǎn)換點(diǎn)的平均深度則位于393km上,而TA給出的轉(zhuǎn)換點(diǎn)深度是390km,這說(shuō)明各臺(tái)網(wǎng)/子臺(tái)網(wǎng)疊加結(jié)果具有很好的一致性.

    表2 疊加圖給出的轉(zhuǎn)換點(diǎn)深度列表Table 2 Conversion depths from the vespegrams

    圖4 Hi-net所記錄地震Ev.5的波形資料各子臺(tái)網(wǎng)4次根傾斜疊加灰度圖灰度圖橫坐標(biāo)為相對(duì)于直達(dá)P的走時(shí)差Δt,縱坐標(biāo)為相對(duì)于直達(dá)P的慢度差Δp.箭頭所對(duì)應(yīng)的數(shù)字為轉(zhuǎn)換震相對(duì)應(yīng)的轉(zhuǎn)換深度,單位為km.Fig.4 The vespegrams of the Ev.5derived from the 4-th root slant stacking with the data by Hi-net sub-networks.H1—H9were shown at Fig.2b The numbers above arrows are the conversion depths(in km)related to SdP phases.

    相較于其他地震,Ev.5給出的轉(zhuǎn)換點(diǎn)深度要淺得多.一般地,轉(zhuǎn)換點(diǎn)深度誤差來(lái)源于震源深度誤差(±5km)、相對(duì)走時(shí)識(shí)讀誤差(±5km)以及俯沖板塊速度異常引起的誤差(±1km),總體而言,轉(zhuǎn)換點(diǎn)深度誤差在±11km(Collier and Helffrich,1997).各子臺(tái)網(wǎng)相應(yīng)的傾斜疊加結(jié)果圖上的次生震相對(duì)應(yīng)的轉(zhuǎn)換點(diǎn)深度的差異與局部區(qū)域的速度躍變和走時(shí)差讀取誤差有一定關(guān)系(Zang et al.,2006).因此,該地震給出的轉(zhuǎn)換點(diǎn)深度有些許存疑.

    我們將各轉(zhuǎn)換點(diǎn)投影到地球平面上,可以看到410km間斷面存在的深度位于372~398km附近(圖5),即在俯沖板塊及其殘留物質(zhì)的影響下抬升約20~50km.一般地,410km間斷面附近的克拉伯龍(Clapeyron)斜率為-2.9MPa/K(Bina and Helffrich,1994),那么每抬升8km則約有100K的低溫異常(Flanagan and Shearer,1998b),因而相對(duì)于周?chē)蒯6员疚难芯繀^(qū)域下方俯沖物質(zhì)導(dǎo)致了250~630K的低溫異常.

    圖5 410km間斷面相關(guān)的SdP轉(zhuǎn)換點(diǎn)深度分布三角形為轉(zhuǎn)換點(diǎn)的位置,其附近的數(shù)字為轉(zhuǎn)換點(diǎn)平均深度(表2).其他說(shuō)明同圖1.Fig.5 The SdP conversion points related to the 410km discontinuityTriangles are for the conversion points locations and the related numbers are for the average depths of the conversion points(Table 2).All others are same as Fig.1.

    基于地質(zhì)學(xué)與地球化學(xué)的研究成果,本文研究區(qū)域陸陸碰撞的時(shí)間大致在102~85Ma(Ali et al.,2002),因而海洋巖石圈俯沖殘留應(yīng)該在此之前就已經(jīng)進(jìn)入地球深部了.同時(shí),從地震層析成像結(jié)果來(lái)看(Koulakov and Sobolev,2006;Villase?or et al.,2003,轉(zhuǎn)引自Negredo et al.,2007),帕米爾—興都庫(kù)什俯沖板塊在地球內(nèi)部250km深度處存在著減薄形態(tài),表明俯沖板塊處于拆離過(guò)程中,最深部可能到達(dá)地幔轉(zhuǎn)換區(qū)中.俯沖板塊的殘留體會(huì)給地球深部帶來(lái)不同的物質(zhì),如石榴子石、輝石,甚至水等揮發(fā)分物質(zhì)(Weidner and Wang,2000),這會(huì)導(dǎo)致在α相橄欖石到β相橄欖石相變之外,還有石榴子石和輝石等的相變;同時(shí)水的存在也會(huì)導(dǎo)致部分熔融的發(fā)生,進(jìn)而影響到410km間斷面相變線的起伏形態(tài).

    因而,相比于410km間斷面抬升達(dá)60~70km(如 Collier and Helffrich,1997;Collier et al.,2001)的西太平洋俯沖區(qū)域而言,本研究區(qū)域下方的俯沖殘留體帶來(lái)的低溫異常要小,相應(yīng)的抬升也就要小不少.當(dāng)然俯沖殘留物質(zhì)本身也會(huì)導(dǎo)致該區(qū)域物質(zhì)構(gòu)成的調(diào)整,而這需要進(jìn)一步更精細(xì)的地震學(xué)工作,并在此基礎(chǔ)上,結(jié)合高溫高壓礦物學(xué)成果進(jìn)行分析.

    4 結(jié)論

    本文利用日本Hi-net地震臺(tái)網(wǎng)和美國(guó)TA臺(tái)陣記錄的帕米爾—興都庫(kù)什俯沖區(qū)域的6個(gè)震源深度為154.0~220.9km、震級(jí)為Mb5.6~6.4的中、深源地震的短周期/寬頻帶波形資料,經(jīng)過(guò)4次根傾斜疊加處理,獲得了36組Hi-net子臺(tái)網(wǎng)和TA記錄資料的傾斜疊加灰度圖,從中提取了與410km間斷面相關(guān)的次生轉(zhuǎn)換震相SdP,發(fā)現(xiàn)受俯沖殘留體影響下的410km間斷面的深度位于372~398km.較之持續(xù)俯沖的西太平洋地區(qū)海洋巖石圈對(duì)于410km間斷面的相變線的影響要小得多.

    致謝 感謝兩位評(píng)閱人寶貴意見(jiàn).

    Ali A,Nakamura E,Yamamoto H.2002.Sm-Nd mineral ages of pegmatite veins and their host rocks from Swat area,Chilas complex,northern Pakistan.J.AsianEarthSci.,21(3):331-339.

    Bai L,Zhang T Z.2015.Complex deformation pattern of the Pamir-Hindu Kush region inferred from multi-scale double-difference earthquake relocations.Tectonophysics,638:177-184,doi:10.1016/j.tecto.2014.11.006.

    Bina C R,Helffrich G.1994.Phase transition Clapeyron slopes and transition zone seismic discontinuity topography.J.Geophys.Res.,99(B8):15853-15860.

    Castle J C,Creager K C.2000.Local sharpness and shear wave speed jump across the 660km discontinuity.J.Geophys.Res.,105(B3):6191-6200.

    Collier J D,Helffrich G R.1997.Topography of the“410”and“660”km seismic discontinuities in the Izu-Bonin Subduction Zone.Geophys.Res.Lett.,24(12):1535-1538.

    Collier J D,Helffrich G R,Wood B J.2001.Seismic discontinuities and subduction zones.Phys.EarthPlanet.Inter.,127(1-4):35-49.

    Deuss A.2009.Global observations of mantle discontinuities using SS and PP precursors.Surv.Geophy.,30(4-5):301-326.

    Dziewonski A M,Anderson D L.1981.Preliminary reference Earth model.Phys.EarthPlanet.Inter.,25(4):297-356.

    Engdahl E R,van der Hilst R,Buland R.1998.Global teleseismic earthquake relocation with improved travel times and procedures for depth determination.Bull.Seism.Soc.Am.,88(3):722-743.

    Flanagan M P,Shearer P M.1998a.Global mapping of topography on transition zone velocity discontinuities by stacking SS precursors.J.Geophys.Res.,103(B2):2673-2692.

    Flanagan M P,Shearer P M.1998b.Topography on the 410-km seismic velocity discontinuity near subduction zones from stacking of sS,sP,and pP precursors.J.Geophys.Res.,103(B9):21165-21182.

    Fukao Y,Obayashi M.2013.Subducted slabs stagnant above,penetrating through,and trapped below the 660km discontinuity.J.Geophys.Res.,118(11):5920-5938.

    Gudmundssonó,Sambridge M.1998.A regionalized upper mantle(RUM)seismic model.J.Geophys.Res.,103(B4):7121-7136.

    Jin Z M,Zhang Z J,Wu Y,et al.2013.Mantle transition zone and research methods.//Ding Z L ed.Research Methods in Solid Earth(in Chinese).Beijing:Science Publishing Company:402-434.

    Kanasewich E R,Hemmings C D,Alpaslan T.1973.Nth-root stack nonlinear multichannel filter.Geophysics,38(2):327-338.

    Kennett B L N,Engdahl E R.1991.Traveltimes for global earthquake location and phase identification.Geophys.J.Int.,105(2):429-465.

    Kennett B L N,Engdahl E R,Buland R.1995.Constraints on seismic velocities in the Earth from traveltimes.Geophys.J.Int.,122(1):108-124.

    Koulakov I,Sobolev S V.2006.A tomographic image of Indian lithosphere break-off beneath the Pamir-Hindu Kush region.Geophys.J.Int.,164(2):425-440.

    Li X,Sobolev S V,Kind R,et al.2000.A detailed receiver function image of the upper mantle discontinuities in the Japan subduction zone.EarthPlanet.Sci.Lett.,183(3-4):527-541.

    McFadden P L,Drummond B J,Karvis S.1986.The N-th root stack:theory,application and examples.Geophysics,51(10):1879-1892.

    Negredo A M,Replumaz A,Villase?or A,et al.2007.Modeling the evolution of continental subduction processes in the Pamir-Hindu Kush region.EarthPlanet.Sci.Lett.,259(1-2):212-225.

    Ning J Y,Zang S X.1990.The distribution of Earthquakes and Stress State in Pamir-Hindu Kush Regions.ChineseJ.Geophys.(in Chinese),33(6):657-669.

    Obayashi M,Sugioka H,Yoshimitsu J,et al.2006.High temperature anomalies oceanward of subducting slabs at the 410km discontinuity.EarthPlanet.Sci.Lett.,243(1-2):149-158.

    Pegler G,Das S.1998.An enhanced image of the Pamir-Hindu Kush seismic zone from relocated earthquake hypocentres.Geophys.J.Int.,134(2):573-595.

    Revenaugh J,Jordan T H.1991.Mantle layering from ScS reverberations:1.Waveform inversion of zeroth-order reverberations.J.Geophys.Res.,96(B12):19749-19762.

    Rost S,Thomas C.2002.Array seismology: methods and applications.Rev.Geophys.,40(3):2-1—2-27,doi:10.1029/2000RG000100.

    Sippl C,Schurr B,Yuan X,et al.2013.Geometry of the Pamir-Hindu Kush intermediate-depth earthquake zone from local seismic data.J.Geophys.Res.,118(4):1438-1457.

    Song T R A,Helmberger D V,Grand S P.2004.Low-velocity zone atop the 410-km seismic discontinuity in the northwestern United States.Nature,427(6974):530-533.

    Sun W B,He Y S,Chang Z,et al.2009.The plate subduction and stress state in the Pamir-Hindu Kush region.Seism.Geol.(in Chinese),31(2):207-217.

    Tian X B,Zhao D P,Zhang H S,et al.2010.Mantle transition zone topography and structure beneath the central Tien Shan orogenic belt.J.Geophys.Res.,115(B10):B10308,doi:10.1029/2008JB006229.

    Tian X B,Teng J W,Zhang H S,et al.2011.Structure of crust and upper mantle beneath the Ordos Block and the Yinshan Mountains revealed by receiver function analysis.Phys.Earth Planet.Inter.,184(3-4):186-193.

    Villase?or A,Spakman W,Engdahl E R.2003.Influence of regional travel times in global tomographic models.EGS-AGUEUG Joint Assembly,Abstracts from the meeting held in Nice,F(xiàn)rance,6-11April 2003,abstract 8614.

    Weidner D J,Wang Y B.2000.Phase transformations:implications for mantle structure.//Karaoto S I,F(xiàn)orte A,Liebermann R,et al.,Earth′s Deep Interior:Mineral Physics and Tomography from the Atomic to the Global Scale, Vol. Geophysical Monograph Series 117.American Geophysical Union:215-235.Zang S X,Zhou Y Z.2002.N-th root slant stack and its application in study of mantle discontinuities.ChineseJ.Geophys.(in Chinese),45(3):407-415,doi:10.3321/j.issn:0001-5733.2002.03.012.

    Zang S X,Zhou Y Z,Ning J Y,et al.2006.Multiple discontinuities near 660km beneath Tonga area.Geophys.Res.Lett.,33(20):L20312,doi:10.1029/2006GL027262.

    Zhang L P,Shao Z G,Li Z H.2014.Dynamic action of mutual subduction between the Indian and the Eurasia plates in Hindu Kush-Pamir region.ChineseJ.Geophys.(in Chinese),57(2):459-471,doi:10.6038/cjg20140212.

    Zhang Z J,Yuan X H,Chen Y,et al.2010.Seismic signature of the collision between the east Tibetan escape flow and the Sichuan Basin.EarthPlanet.Sci.Lett.,292(3-4):254-264.

    附中文參考文獻(xiàn)

    金振民,張忠杰,吳耀等.2013.地幔轉(zhuǎn)換帶和研究方法.//丁仲禮主編.固體地球科學(xué)研究方法.北京:科學(xué)出版社:402-434.

    寧杰遠(yuǎn),臧紹先.1990.帕米爾—興都庫(kù)什地區(qū)地震空間分布特征及應(yīng)力場(chǎng)特征.地球物理學(xué)報(bào),33(6):657-669.

    孫文斌,和躍時(shí),常征等.2009.帕米爾—興都庫(kù)什地區(qū)板塊俯沖及其應(yīng)力狀態(tài).地震地質(zhì),31(2):207-217.

    臧紹先,周元澤.2002.N次根傾斜疊加在地幔間斷面研究中的應(yīng)用.地球物理學(xué)報(bào),45(3):407-415,doi:10.3321/j.issn:0001-5733.2002.03.012.

    張浪平,邵志剛,李志海.2014.印度板塊與歐亞板塊在興都庫(kù)什—帕米爾地區(qū)相互俯沖的動(dòng)力作用分析.地球物理學(xué)報(bào),57(2):459-471,doi:10.6038/cjg20140212.

    猜你喜歡
    臺(tái)網(wǎng)震源波形
    地球物理臺(tái)網(wǎng)儀器維修信息管理的研究與實(shí)現(xiàn)
    地震研究(2021年1期)2021-04-13 01:05:30
    對(duì)《壓力容器波形膨脹節(jié)》2018版新標(biāo)準(zhǔn)的理解及分析
    推進(jìn)報(bào)臺(tái)網(wǎng)深度融合 做強(qiáng)區(qū)級(jí)融媒體中心
    基于LFM波形的靈巧干擾效能分析
    震源的高返利起步
    西藏地震應(yīng)急流動(dòng)臺(tái)網(wǎng)淺析
    西藏科技(2016年9期)2016-09-26 12:21:39
    基于ARM的任意波形電源設(shè)計(jì)
    MRG9000媒資衛(wèi)士在市級(jí)電視臺(tái)全臺(tái)網(wǎng)的應(yīng)用
    大連臺(tái)使用CTS-1記錄波形特點(diǎn)
    可控震源地震在張掖盆地南緣逆沖斷裂構(gòu)造勘探中的應(yīng)用
    日本免费a在线| 久99久视频精品免费| 国产成人啪精品午夜网站| 久久久久久久久中文| 国产视频一区二区在线看| 国产又爽黄色视频| 亚洲黑人精品在线| 熟女少妇亚洲综合色aaa.| av电影中文网址| 黄色 视频免费看| 80岁老熟妇乱子伦牲交| 日韩高清综合在线| 国产精品永久免费网站| 国产精品亚洲一级av第二区| 宅男免费午夜| 成人黄色视频免费在线看| 国产片内射在线| 成年版毛片免费区| 久久国产乱子伦精品免费另类| 欧美日韩中文字幕国产精品一区二区三区 | 波多野结衣av一区二区av| 亚洲av日韩精品久久久久久密| 如日韩欧美国产精品一区二区三区| 久久久久精品国产欧美久久久| av电影中文网址| 亚洲伊人色综图| av有码第一页| 国产亚洲欧美98| 久久久久久久久中文| 亚洲全国av大片| 国产一区二区在线av高清观看| a在线观看视频网站| 久久伊人香网站| 亚洲专区字幕在线| av在线播放免费不卡| 91av网站免费观看| 免费久久久久久久精品成人欧美视频| 老熟妇乱子伦视频在线观看| 丰满的人妻完整版| 中文欧美无线码| 国产高清国产精品国产三级| 亚洲 欧美 日韩 在线 免费| 国产一区二区三区综合在线观看| 国产成人av教育| 亚洲精品av麻豆狂野| 亚洲国产欧美日韩在线播放| 精品午夜福利视频在线观看一区| 国产成+人综合+亚洲专区| 大型黄色视频在线免费观看| 久久精品aⅴ一区二区三区四区| 曰老女人黄片| 国产高清videossex| 一级片'在线观看视频| 制服人妻中文乱码| 精品国产乱子伦一区二区三区| 午夜91福利影院| 国产精品98久久久久久宅男小说| 中文亚洲av片在线观看爽| 久久人人97超碰香蕉20202| 日本欧美视频一区| 在线免费观看的www视频| 国产免费现黄频在线看| 亚洲国产欧美一区二区综合| 人人妻人人澡人人看| 免费av中文字幕在线| 亚洲av五月六月丁香网| 国产91精品成人一区二区三区| 18禁裸乳无遮挡免费网站照片 | 99国产极品粉嫩在线观看| www.999成人在线观看| 久久性视频一级片| 99国产精品99久久久久| 欧美一区二区精品小视频在线| 国产精品综合久久久久久久免费 | 亚洲,欧美精品.| 18禁国产床啪视频网站| 交换朋友夫妻互换小说| www国产在线视频色| 一级毛片高清免费大全| 夜夜看夜夜爽夜夜摸 | 国产欧美日韩精品亚洲av| 一个人免费在线观看的高清视频| 亚洲人成77777在线视频| 亚洲五月天丁香| 人人妻人人澡人人看| 法律面前人人平等表现在哪些方面| 一级毛片女人18水好多| 亚洲午夜理论影院| 午夜激情av网站| 91字幕亚洲| 免费av毛片视频| 亚洲人成77777在线视频| 高清毛片免费观看视频网站 | 国产精品乱码一区二三区的特点 | 他把我摸到了高潮在线观看| 精品熟女少妇八av免费久了| 国产精品乱码一区二三区的特点 | 又大又爽又粗| 男男h啪啪无遮挡| 亚洲伊人色综图| 欧美激情高清一区二区三区| 伦理电影免费视频| 成人国产一区最新在线观看| 久久久精品欧美日韩精品| 十分钟在线观看高清视频www| 欧美性长视频在线观看| 日韩欧美一区二区三区在线观看| 69精品国产乱码久久久| 欧美 亚洲 国产 日韩一| 在线观看一区二区三区激情| 99riav亚洲国产免费| 黄片播放在线免费| 久久人人97超碰香蕉20202| av在线天堂中文字幕 | 精品人妻1区二区| 亚洲av第一区精品v没综合| 熟女少妇亚洲综合色aaa.| 亚洲欧美日韩另类电影网站| 国产精品二区激情视频| 十分钟在线观看高清视频www| 一级毛片女人18水好多| 中文字幕最新亚洲高清| 51午夜福利影视在线观看| 大香蕉久久成人网| 黄片播放在线免费| 两性午夜刺激爽爽歪歪视频在线观看 | 两人在一起打扑克的视频| 久久影院123| 国产精品98久久久久久宅男小说| 午夜影院日韩av| 亚洲久久久国产精品| 久久精品国产亚洲av高清一级| 日韩大码丰满熟妇| 欧美成人免费av一区二区三区| 国产成+人综合+亚洲专区| 精品国产亚洲在线| av欧美777| 国产有黄有色有爽视频| 女警被强在线播放| 窝窝影院91人妻| 丰满饥渴人妻一区二区三| 欧美激情 高清一区二区三区| 母亲3免费完整高清在线观看| 欧美日韩亚洲国产一区二区在线观看| 又黄又爽又免费观看的视频| www.999成人在线观看| 天堂√8在线中文| 一级作爱视频免费观看| 欧美人与性动交α欧美精品济南到| 黄色怎么调成土黄色| 国产精品综合久久久久久久免费 | 一进一出好大好爽视频| 大香蕉久久成人网| 黄色毛片三级朝国网站| 最好的美女福利视频网| а√天堂www在线а√下载| 成人三级做爰电影| 国产精品1区2区在线观看.| 1024香蕉在线观看| 一级作爱视频免费观看| 亚洲欧美一区二区三区久久| 亚洲国产精品sss在线观看 | 精品国产乱子伦一区二区三区| 18禁观看日本| 99国产精品免费福利视频| 性色av乱码一区二区三区2| 日韩高清综合在线| 校园春色视频在线观看| 色播在线永久视频| 妹子高潮喷水视频| av视频免费观看在线观看| 18禁美女被吸乳视频| 夜夜夜夜夜久久久久| 怎么达到女性高潮| 麻豆成人av在线观看| 午夜影院日韩av| 视频区欧美日本亚洲| tocl精华| 在线播放国产精品三级| 欧美日韩精品网址| 久久国产精品男人的天堂亚洲| 十八禁人妻一区二区| 看片在线看免费视频| 一级a爱片免费观看的视频| 日韩高清综合在线| 国产亚洲精品一区二区www| 国产av一区二区精品久久| 别揉我奶头~嗯~啊~动态视频| 欧美 亚洲 国产 日韩一| 免费在线观看完整版高清| 日本wwww免费看| 午夜激情av网站| 国产三级黄色录像| 国产高清国产精品国产三级| 久久国产乱子伦精品免费另类| 99riav亚洲国产免费| 91在线观看av| 久久精品影院6| 自拍欧美九色日韩亚洲蝌蚪91| 夜夜看夜夜爽夜夜摸 | 人妻丰满熟妇av一区二区三区| 久久人人精品亚洲av| 一边摸一边抽搐一进一小说| 亚洲专区字幕在线| 深夜精品福利| 多毛熟女@视频| 男人舔女人的私密视频| 久久性视频一级片| 国产精品99久久99久久久不卡| 国产一区二区在线av高清观看| 国产亚洲欧美精品永久| 99国产极品粉嫩在线观看| 国产黄色免费在线视频| 高清欧美精品videossex| 欧美黄色片欧美黄色片| 精品日产1卡2卡| 麻豆久久精品国产亚洲av | 天堂俺去俺来也www色官网| 天天添夜夜摸| 日韩大尺度精品在线看网址 | 男女下面插进去视频免费观看| 人妻丰满熟妇av一区二区三区| 在线观看一区二区三区激情| 极品人妻少妇av视频| 变态另类成人亚洲欧美熟女 | e午夜精品久久久久久久| 久久国产精品影院| 亚洲av日韩精品久久久久久密| 国内久久婷婷六月综合欲色啪| 男人舔女人的私密视频| tocl精华| 亚洲欧美激情在线| 日韩人妻精品一区2区三区| 中文字幕高清在线视频| 亚洲成av片中文字幕在线观看| 可以免费在线观看a视频的电影网站| 动漫黄色视频在线观看| 色精品久久人妻99蜜桃| 欧美人与性动交α欧美软件| 久久精品影院6| 欧美中文综合在线视频| 精品久久久久久,| 久9热在线精品视频| 亚洲一区中文字幕在线| 99久久精品国产亚洲精品| 午夜两性在线视频| 交换朋友夫妻互换小说| 国产免费男女视频| 黑丝袜美女国产一区| 久久国产精品人妻蜜桃| 久久精品国产亚洲av香蕉五月| 国产精品秋霞免费鲁丝片| 国产亚洲精品久久久久久毛片| 极品教师在线免费播放| 亚洲欧美日韩高清在线视频| 一进一出好大好爽视频| www.精华液| 免费看a级黄色片| 两性夫妻黄色片| 手机成人av网站| 精品久久久久久久毛片微露脸| 精品第一国产精品| 久久精品aⅴ一区二区三区四区| 极品人妻少妇av视频| 妹子高潮喷水视频| 桃红色精品国产亚洲av| 水蜜桃什么品种好| 男女下面插进去视频免费观看| 丝袜美足系列| 亚洲精品粉嫩美女一区| 我的亚洲天堂| 国产91精品成人一区二区三区| 91麻豆av在线| 波多野结衣av一区二区av| 国产精品美女特级片免费视频播放器 | 午夜精品久久久久久毛片777| 欧美日韩亚洲国产一区二区在线观看| 欧美国产精品va在线观看不卡| 国产真人三级小视频在线观看| 国产极品粉嫩免费观看在线| 精品卡一卡二卡四卡免费| 国产精品综合久久久久久久免费 | 丰满的人妻完整版| 国产精华一区二区三区| 看免费av毛片| 亚洲伊人色综图| 精品免费久久久久久久清纯| 欧美一级毛片孕妇| 亚洲专区国产一区二区| 国产精品一区二区三区四区久久 | 亚洲一区二区三区色噜噜 | 国产麻豆69| 久久精品亚洲精品国产色婷小说| 黄色女人牲交| 欧美中文综合在线视频| 日韩精品免费视频一区二区三区| 久久草成人影院| 黑人巨大精品欧美一区二区蜜桃| 日韩高清综合在线| 色综合婷婷激情| 日本 av在线| 男女下面插进去视频免费观看| 国产午夜精品久久久久久| 日本撒尿小便嘘嘘汇集6| 亚洲熟女毛片儿| 深夜精品福利| 在线观看免费日韩欧美大片| 人妻丰满熟妇av一区二区三区| 一进一出抽搐gif免费好疼 | 午夜成年电影在线免费观看| 国产99白浆流出| 男女之事视频高清在线观看| 淫妇啪啪啪对白视频| 亚洲自拍偷在线| 久久精品亚洲av国产电影网| 9热在线视频观看99| 91av网站免费观看| 欧美人与性动交α欧美软件| 天堂中文最新版在线下载| 中文字幕高清在线视频| 国产精品爽爽va在线观看网站 | 青草久久国产| 亚洲五月色婷婷综合| 亚洲一码二码三码区别大吗| 欧美激情久久久久久爽电影 | 亚洲第一欧美日韩一区二区三区| 美女大奶头视频| 国产97色在线日韩免费| 在线视频色国产色| 国产成人欧美在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 日本vs欧美在线观看视频| 国产欧美日韩一区二区三区在线| av福利片在线| 91字幕亚洲| av福利片在线| 国产深夜福利视频在线观看| 久久香蕉激情| 欧美亚洲日本最大视频资源| 一级作爱视频免费观看| 久久婷婷成人综合色麻豆| 桃红色精品国产亚洲av| 国产一区二区三区视频了| 又大又爽又粗| 妹子高潮喷水视频| 久久香蕉激情| 嫩草影视91久久| 女警被强在线播放| 婷婷精品国产亚洲av在线| 国产精品秋霞免费鲁丝片| 999久久久国产精品视频| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美一区二区精品小视频在线| 国产熟女午夜一区二区三区| 自线自在国产av| 欧美日韩一级在线毛片| 嫩草影视91久久| 99久久国产精品久久久| 女人高潮潮喷娇喘18禁视频| 老司机午夜十八禁免费视频| 国产三级在线视频| 成年女人毛片免费观看观看9| 精品久久久久久,| 久久久国产成人精品二区 | 国产人伦9x9x在线观看| 美女高潮喷水抽搐中文字幕| 一级毛片女人18水好多| videosex国产| 精品国产国语对白av| 天天添夜夜摸| 日韩免费高清中文字幕av| 久久人妻福利社区极品人妻图片| 两性午夜刺激爽爽歪歪视频在线观看 | 看免费av毛片| 午夜福利免费观看在线| 国产成+人综合+亚洲专区| 免费观看人在逋| 久久精品亚洲熟妇少妇任你| 久久久久久久久久久久大奶| 亚洲成人免费电影在线观看| 欧美人与性动交α欧美精品济南到| 女人被躁到高潮嗷嗷叫费观| 别揉我奶头~嗯~啊~动态视频| 纯流量卡能插随身wifi吗| 欧美日韩一级在线毛片| 国产av在哪里看| 夜夜夜夜夜久久久久| 亚洲成人精品中文字幕电影 | 97碰自拍视频| 亚洲一区高清亚洲精品| 在线观看一区二区三区| 国产高清videossex| 91麻豆精品激情在线观看国产 | 黑丝袜美女国产一区| 亚洲一区二区三区不卡视频| 69av精品久久久久久| 亚洲国产精品合色在线| 男女下面插进去视频免费观看| 成年版毛片免费区| 在线观看免费日韩欧美大片| 久久精品91蜜桃| 国产精品久久视频播放| 可以在线观看毛片的网站| 91大片在线观看| 巨乳人妻的诱惑在线观看| 成人黄色视频免费在线看| 国产精品免费视频内射| 日日干狠狠操夜夜爽| 亚洲avbb在线观看| 天天躁夜夜躁狠狠躁躁| 最近最新免费中文字幕在线| 在线观看舔阴道视频| 精品国产美女av久久久久小说| 制服人妻中文乱码| 国产免费av片在线观看野外av| 天堂影院成人在线观看| 午夜福利,免费看| 天堂影院成人在线观看| 国产黄a三级三级三级人| 午夜日韩欧美国产| 51午夜福利影视在线观看| 欧美日韩av久久| 日韩高清综合在线| ponron亚洲| 免费不卡黄色视频| 日韩 欧美 亚洲 中文字幕| 男女高潮啪啪啪动态图| 久久精品影院6| 久久久久久久久中文| 少妇裸体淫交视频免费看高清 | 亚洲精品久久成人aⅴ小说| 国产片内射在线| 天堂动漫精品| 久久精品国产综合久久久| 99国产综合亚洲精品| www日本在线高清视频| 母亲3免费完整高清在线观看| 亚洲中文字幕日韩| 日本一区二区免费在线视频| 美女大奶头视频| 日本三级黄在线观看| 欧美黄色片欧美黄色片| 国产成人精品久久二区二区免费| 午夜影院日韩av| 国产麻豆69| 亚洲国产毛片av蜜桃av| 免费av中文字幕在线| 中文字幕人妻丝袜一区二区| 亚洲av成人一区二区三| 久久99一区二区三区| 精品少妇一区二区三区视频日本电影| 成人黄色视频免费在线看| 日韩欧美国产一区二区入口| 亚洲人成电影免费在线| 老汉色∧v一级毛片| 国产精品av久久久久免费| 69av精品久久久久久| 丝袜美足系列| 久久香蕉国产精品| 国产在线精品亚洲第一网站| 一进一出好大好爽视频| 国产高清videossex| 国产av一区二区精品久久| 亚洲精品国产一区二区精华液| 久久人人97超碰香蕉20202| 青草久久国产| 亚洲欧美激情在线| 这个男人来自地球电影免费观看| 777久久人妻少妇嫩草av网站| 男女下面插进去视频免费观看| 中文字幕精品免费在线观看视频| 黄色视频,在线免费观看| 啦啦啦在线免费观看视频4| 亚洲,欧美精品.| 国产精品 欧美亚洲| 在线视频色国产色| 人人妻人人澡人人看| 午夜精品国产一区二区电影| 欧美色视频一区免费| 在线观看www视频免费| 亚洲一区二区三区欧美精品| 曰老女人黄片| 日韩av在线大香蕉| 在线十欧美十亚洲十日本专区| 少妇裸体淫交视频免费看高清 | 免费观看精品视频网站| 深夜精品福利| 亚洲情色 制服丝袜| 日韩欧美三级三区| 91大片在线观看| 亚洲成人国产一区在线观看| 91大片在线观看| 欧美老熟妇乱子伦牲交| 亚洲精品国产色婷婷电影| 另类亚洲欧美激情| 成人18禁在线播放| 少妇被粗大的猛进出69影院| 脱女人内裤的视频| 狠狠狠狠99中文字幕| 午夜两性在线视频| 亚洲成a人片在线一区二区| 国产精品久久久人人做人人爽| 久久 成人 亚洲| 大型av网站在线播放| 亚洲精品中文字幕在线视频| 中文字幕色久视频| 亚洲男人天堂网一区| 国产熟女xx| a级毛片在线看网站| 精品电影一区二区在线| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲男人天堂网一区| 美国免费a级毛片| 在线观看午夜福利视频| 午夜免费成人在线视频| www.熟女人妻精品国产| 色综合站精品国产| 黄色成人免费大全| 无人区码免费观看不卡| 自线自在国产av| 两人在一起打扑克的视频| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久成人av| 欧美日本中文国产一区发布| 亚洲一区中文字幕在线| 欧美日韩福利视频一区二区| 国产不卡一卡二| 精品国产乱子伦一区二区三区| 欧美激情极品国产一区二区三区| 中文字幕高清在线视频| 一级,二级,三级黄色视频| 好男人电影高清在线观看| 欧美中文综合在线视频| 国产免费现黄频在线看| 精品少妇一区二区三区视频日本电影| aaaaa片日本免费| 久久精品国产清高在天天线| 精品卡一卡二卡四卡免费| 久久精品成人免费网站| av中文乱码字幕在线| 91大片在线观看| 精品国产乱码久久久久久男人| 亚洲熟女毛片儿| 亚洲免费av在线视频| 在线免费观看的www视频| 国产成+人综合+亚洲专区| 国产色视频综合| 亚洲国产精品999在线| 国产成人av激情在线播放| 日本a在线网址| 欧美日韩亚洲高清精品| 午夜成年电影在线免费观看| 国产成人精品在线电影| 午夜两性在线视频| 少妇的丰满在线观看| 亚洲成av片中文字幕在线观看| 欧美国产精品va在线观看不卡| 久久香蕉激情| 婷婷六月久久综合丁香| 91在线观看av| 成人手机av| 亚洲国产精品sss在线观看 | 宅男免费午夜| 国产免费男女视频| 欧美人与性动交α欧美精品济南到| 亚洲国产看品久久| 看黄色毛片网站| 国产成人影院久久av| 日韩成人在线观看一区二区三区| 一级毛片女人18水好多| 一个人观看的视频www高清免费观看 | 一级a爱视频在线免费观看| 精品久久久久久电影网| 国产精品免费一区二区三区在线| 免费在线观看视频国产中文字幕亚洲| 久久香蕉激情| 国产av一区二区精品久久| 久久青草综合色| 亚洲av电影在线进入| 亚洲国产精品合色在线| 18禁黄网站禁片午夜丰满| 欧美午夜高清在线| 亚洲精品av麻豆狂野| 午夜亚洲福利在线播放| 亚洲五月色婷婷综合| 亚洲一区高清亚洲精品| 久久久水蜜桃国产精品网| av免费在线观看网站| 亚洲免费av在线视频| 国产一区二区三区综合在线观看| 九色亚洲精品在线播放| 久久久久久亚洲精品国产蜜桃av| 老汉色av国产亚洲站长工具| 国产乱人伦免费视频| e午夜精品久久久久久久| 91精品三级在线观看| 手机成人av网站| 97超级碰碰碰精品色视频在线观看| 中文字幕人妻丝袜一区二区| 国产乱人伦免费视频| 国产亚洲欧美98| 国产精品一区二区三区四区久久 | 正在播放国产对白刺激| 中文字幕另类日韩欧美亚洲嫩草| 丝袜在线中文字幕| 在线播放国产精品三级| 男女之事视频高清在线观看| 大型黄色视频在线免费观看| 欧美大码av| 亚洲一码二码三码区别大吗| 99国产极品粉嫩在线观看| 夜夜躁狠狠躁天天躁| 母亲3免费完整高清在线观看| 午夜影院日韩av| 日日摸夜夜添夜夜添小说|