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

    卡羅琳板塊及其附近地區(qū)的巖石圈有效彈性厚度

    2016-11-16 00:54:57楊安付永濤李安春
    地球物理學(xué)報(bào) 2016年9期
    關(guān)鍵詞:馬里亞納海山卡羅琳

    楊安, 付永濤, 李安春

    1 中國科學(xué)院 海洋研究所, 青島 266071 2 中國科學(xué)院 海洋地質(zhì)與環(huán)境重點(diǎn)實(shí)驗(yàn)室, 青島 266071

    ?

    卡羅琳板塊及其附近地區(qū)的巖石圈有效彈性厚度

    楊安1,2, 付永濤1,2, 李安春1,2

    1 中國科學(xué)院 海洋研究所, 青島 266071 2 中國科學(xué)院 海洋地質(zhì)與環(huán)境重點(diǎn)實(shí)驗(yàn)室, 青島 266071

    本文利用卡羅琳板塊及其附近地區(qū)的自由空氣重力異常和海水深度數(shù)據(jù),結(jié)合滑動(dòng)窗口導(dǎo)納技術(shù)(MWAT),計(jì)算了該地區(qū)的巖石圈有效彈性厚度Te.本文使用Multitaper(多窗譜)方法對(duì)功率譜密度進(jìn)行估計(jì),基于實(shí)際的海底地形,通過模擬計(jì)算得到了MWAT方法較真值的改正,MWAT方法計(jì)算的結(jié)果偏小20%左右.研究結(jié)果顯示卡羅琳板塊及其附近地區(qū)的Te變化范圍為1~34 km.研究區(qū)域包括了海山、海底高原、俯沖帶、擴(kuò)張洋脊等多種構(gòu)造,對(duì)它們的巖石圈強(qiáng)度的研究為認(rèn)識(shí)西太平洋地區(qū)巖石圈的構(gòu)造和演化提供了重要的依據(jù).Te與加載時(shí)的巖石圈年齡、地表熱流相關(guān).Te與海底地殼年齡之間的關(guān)系顯示Te主要位于板塊冷卻模型的450 ℃的等溫線深度以上.西太平洋的Magellan海山和Marcus-Wake Guyots (MWG)地區(qū)的Te主要分布在加載形成時(shí)板塊冷卻模型的200 ℃的等溫線深度附近,較低的等溫線可能受太平洋超級(jí)地幔柱的影響.我們的研究結(jié)果也顯示在研究區(qū)域內(nèi)海洋地殼的熱流與Te之間存在一定的反相關(guān)性.

    有效彈性厚度;撓曲;卡羅琳板塊;巖石圈;滑動(dòng)窗口導(dǎo)納技術(shù)

    1 引言

    卡羅琳板塊位于太平洋西部邊緣,在早期卡羅琳板塊被認(rèn)為是太平洋板塊的一部分,但基于海底地形、重力、地震、板塊運(yùn)動(dòng)等資料的詳細(xì)研究顯示卡羅琳板塊與鄰近的太平洋板塊的運(yùn)動(dòng)歷史以及構(gòu)造特征明顯不一致,卡羅琳板塊被認(rèn)為是一個(gè)獨(dú)立的板塊(Weissel and Anderson, 1978).卡羅琳板塊主要包括東卡羅琳盆地和西卡羅琳盆地,以及分開它們的歐里皮克海隆(圖1).卡羅琳板塊在東和東北邊與太平洋板塊相接,東邊與太平洋板塊的邊界是Mussau海溝,卡羅琳板塊在此俯沖于太平洋板塊之下;東北邊與太平洋板塊的邊界沿著Sorol海槽,其把卡羅琳洋脊分開,呈現(xiàn)出左旋走滑和張裂的特征(Altis, 1999; Bird, 2003).卡羅琳板塊在西和西北邊與菲律賓海板塊相接,在西北邊是Yap海溝,卡羅琳板塊在此以比較緩慢的速度俯沖于菲律賓海板塊之下(Sato et al., 1997);在西邊是Ayu海槽,在此處卡羅琳板塊與菲律賓海板塊沿東西方向緩慢分裂(Fujiwara et al., 1995; Bird, 2003).卡羅琳板塊在南邊與澳大利亞板塊相接,在西南邊俯沖于新幾內(nèi)亞之下,在東南邊俯沖于北俾斯麥板塊之下,形成Manus海溝,其目前還不是一個(gè)成熟的俯沖帶(Bird, 2003).卡羅琳板塊非常年輕(Müller et al., 2008),大約形成于第三紀(jì)中期(Hill and Hegarty, 1987; Bird, 2003).關(guān)于卡羅琳板塊的成因以及構(gòu)造背景仍然存在許多不清楚的地方,其可能是在漸新世時(shí)期形成的弧后盆地,也可能是捕獲的海洋地殼(Hill and Hegarty, 1987; Altis, 1999).卡羅琳洋脊位于卡羅琳-太平洋板塊交界處,其東邊是卡羅琳群島,許多研究認(rèn)為它們形成于熱點(diǎn)的作用(Keating et al., 1984; Altis, 1999),但也有研究認(rèn)為它們的形成與俯沖構(gòu)造環(huán)境下的破裂相關(guān)(Rehman et al., 2013).沿著Yap海溝往北是馬里亞納海溝,太平洋板塊在此處俯沖于菲律賓海板塊之下,這里形成了典型的海溝、島弧、弧后盆地系統(tǒng).西菲律賓海盆、帕里西維拉海盆是由于俯沖帶后退而形成的殘余弧后盆地.馬里亞納海槽正在以大約30 mm/a的速率分裂(Bird, 2003),在馬里亞納海槽和馬里亞納海溝之間形成了馬里亞納板塊.俯沖的太平洋巖石圈在馬里亞納海溝附近經(jīng)歷著很大的加載,有著非常大的彈性撓曲,沿著海溝方向巖石圈的強(qiáng)度也有比較大的變化(Zhang et al., 2014).

    巖石圈的有效彈性厚度(Te)是一個(gè)與巖石圈構(gòu)造背景和演化非常相關(guān)的參數(shù),它表征巖石圈的撓曲強(qiáng)度和動(dòng)力學(xué)響應(yīng),對(duì)認(rèn)識(shí)巖石圈的構(gòu)造背景、力學(xué)性質(zhì)以及演化有重要意義(Forsyth, 1985; Watts et al., 2006; 楊亭等, 2013).海洋巖石圈的有效彈性厚度主要由加載時(shí)(海山、洋脊、海溝等)巖石圈的熱結(jié)構(gòu)控制,近似于板塊冷卻模型的450 ℃的等溫線深度(Watts, 1978).如果已知巖石圈的有效彈性厚度和海底的年齡,就可以根據(jù)這個(gè)關(guān)系,推測(cè)加載形成的年齡.海洋巖石圈有效彈性厚度的研究方法主要包括:通過地震折射數(shù)據(jù)直接測(cè)量巖石圈的彈性撓曲(Watts et al., 1985);重力和地形之間的響應(yīng)函數(shù)分析方法(Kalnins and Watts, 2009);對(duì)比理論和觀測(cè)重力之差的正演法(Calmant et al., 1990; Filmer et al., 1993);對(duì)比理論和觀測(cè)海底地形之差的反演法(Watts et al., 2006);三維黏彈性數(shù)值模擬的方法(Zhong and Watts, 2013).通過以上這些研究方法,國內(nèi)外許多學(xué)者對(duì)海洋巖石圈的有效彈性厚度開展了研究并取得了重要的成果(Watts, 1978; Watts et al., 1985; Calmant et al., 1990; Filmer et al., 1993; Lyons et al., 2000; 付永濤等, 2002; Watts et al., 2006; Kalnins and Watts, 2009; 趙俐紅等, 2010; Zhong and Watts, 2013; 胡敏章等, 2015).

    Watts等(2006)采用海底地形反演的方法計(jì)算了全球超過9000座海山底部巖石圈的有效彈性厚度,他們根據(jù)Te的結(jié)果從小到大把海山類型分為“on-ridge”、“flank-ridge”和“off-ridge”,在太平洋有大量的“on-ridge”型火山活動(dòng),從東南的Foundation海山延伸到西北的Hess 隆起,在這些“on-ridge”型火山之間,也存在著一些“flank-ridge”和“off-ridge”型火山.Te與加載時(shí)的巖石圈年齡之間并不是一個(gè)簡(jiǎn)單的對(duì)應(yīng)關(guān)系,可能是一個(gè)更低或更高的等溫線,其原因包括由于多次噴發(fā)而導(dǎo)致測(cè)量的海山年齡不準(zhǔn)確、計(jì)算巖石圈的有效彈性厚度的方法不準(zhǔn)確、巖石圈的黏彈性應(yīng)力釋放、區(qū)域的等溫線不同等(Watts et al., 2006; Kalnins and Watts, 2009).Kalnins和Watts(2009)用滑動(dòng)窗口導(dǎo)納技術(shù)(MWAT)研究了西太平洋地區(qū)的海洋巖石圈有效彈性厚度,由于MWAT方法使用的是有限窗口的觀測(cè)導(dǎo)納與無限窗口的理論導(dǎo)納之間的對(duì)比,會(huì)導(dǎo)致反演結(jié)果偏小,他們通過海山模型的測(cè)試結(jié)果顯示在Te為30 km的時(shí)候有±5 km的不確定度,并給出了該方法的修正公式.他們的結(jié)果顯示西太平洋地區(qū)的巖石圈有效彈性厚度變化范圍為0~50 km,平均厚度為9.4 km,Te主要由加載時(shí)的巖石圈年齡控制,夏威夷-皇帝海山鏈的Te比較大,萊恩海嶺、赫斯海隆、卡羅琳群島等地方的Te比較小.根據(jù)最佳擬合的Te與加載時(shí)的巖石圈年齡之間的關(guān)系,推測(cè)出這些地形特征大多形成于晚侏羅紀(jì)和早白堊紀(jì),在這兩個(gè)時(shí)期有大量的火山活動(dòng)(Kalnins and Watts, 2009).

    國內(nèi)也有一些關(guān)于海洋巖石圈的有效彈性厚度的研究.趙俐紅等(2010)用正演法對(duì)中西太平洋的麥哲倫海山鏈和中太平洋海山群的巖石圈有效彈性厚度進(jìn)行了研究,推測(cè)它們是在白堊紀(jì)期間形成于現(xiàn)今法屬波利尼亞群島處的大規(guī)模熱點(diǎn)群附近.胡敏章等(2015)采用MWAT方法計(jì)算了西北太平洋巖石圈有效彈性厚度,結(jié)果顯示該區(qū)域的Te平均值為13.2 km,主要分布在150 ℃到450 ℃等溫線深度范圍內(nèi),Te并未隨海山加載時(shí)巖石圈年齡增加而增大.蘇達(dá)權(quán)(2012)比較了海洋巖石圈板塊的響應(yīng)函數(shù)與大陸巖石圈的異同,并分析了海洋巖石圈板塊響應(yīng)函數(shù)的特點(diǎn),他的計(jì)算結(jié)果顯示我國南海南沙海域和南海中央海盆巖石圈板塊有效彈性厚度分別為10 km和6~7 km.

    Altis(1999)研究了卡羅琳板塊北部的卡羅琳洋脊和Sorol海槽的巖石圈有效彈性厚度,為5~7 km,推測(cè)它們形成于洋中脊附近的構(gòu)造環(huán)境,形成時(shí)下部的巖石圈比較年輕并且比較弱,并認(rèn)為卡羅琳洋脊大約是在晚漸新世時(shí)由于熱點(diǎn)的火山作用而形成在漸新世時(shí)的海底之上.卡羅琳群島位于卡羅琳洋脊東部,緊挨著卡羅琳洋脊,其巖石圈有效彈性厚度也比較小(Kalnins and Watts, 2009),由于卡羅琳群島下部的太平洋巖石圈非常古老,而卡羅琳群島比較年輕,較小的Te似乎很難與這樣大的加載形成時(shí)的年齡相一致.Rehman等(2013)的研究認(rèn)為卡羅琳群島不是熱點(diǎn)成因的海山鏈,并認(rèn)為它們的形成與俯沖構(gòu)造環(huán)境下的破裂相關(guān).由于破裂會(huì)導(dǎo)致巖石圈強(qiáng)度降低,這種觀點(diǎn)正好可以解釋卡羅琳群島較小的巖石圈有效彈性厚度.

    海洋巖石圈在俯沖帶受到多種不同的加載,馬里亞納海溝附近有非常大的巖石圈的撓曲,在俯沖帶外圍隆起向海溝方向,巖石圈的有效彈性厚度會(huì)由于巖石圈的破裂屈服而減小(Judge and Mcnutt, 1991).對(duì)馬里亞納海溝附近的彈性撓曲的研究顯示外圍隆起向海溝方向的巖石圈有效彈性厚度為19~40 km,而外圍隆起背離海溝方向的巖石圈有效彈性厚度為45~52 km(Zhang et al., 2014).俯沖帶的加載是一個(gè)正在進(jìn)行中的過程,這里很大的巖石圈有效彈性厚度與此處太平洋板塊非常古老的年齡相一致.Kalnins和Watts(2009)、胡敏章等(2015)的研究結(jié)果也顯示在馬里亞納海溝附近的Te非常大.

    這里我們計(jì)劃對(duì)整個(gè)卡羅琳板塊以及馬里亞納海溝附近地區(qū)的巖石圈有效彈性厚度做一個(gè)詳細(xì)的研究,我們采用重力地形導(dǎo)納法反演巖石圈的有效彈性厚度,類似于Kalnins和Watts(2009)采用的方法,但本文采用了Multitaper方法對(duì)功率譜密度進(jìn)行估計(jì).研究范圍為5°S到25°N,125°E到160°E之間(圖1).海洋重力數(shù)據(jù)采用TOPEX衛(wèi)星重力數(shù)據(jù)topex v23(Sandwell and Smith, 2009),海底地形數(shù)據(jù)采用GEBCO(http:∥www.gebco.net).根據(jù)本文計(jì)算的Te結(jié)果,結(jié)合卡羅琳板塊及其附近地區(qū)的洋殼年齡、海山年齡、海底熱流等資料,將對(duì)該地區(qū)的構(gòu)造背景以及演化提供重要的認(rèn)識(shí).

    2 計(jì)算方法

    2.1 觀測(cè)導(dǎo)納

    已知海底地形和重力異常的時(shí)候,根據(jù)導(dǎo)納的定義,在波數(shù)域中的觀測(cè)自由空氣導(dǎo)納的計(jì)算公式為(Kalnins and Watts, 2009):

    (1)

    其中Zobs(k)為觀測(cè)的自由空氣導(dǎo)納,Sfh(k)和Shh(k)分別表示地形和重力異常的互功率譜密度和地形的自功率譜密度,〈〉表示一個(gè)環(huán)形波數(shù)帶上的平均值.H(k)和ΔGF(k)分別表示地形h和自由空氣重力異常f的傅里葉變換,*表示復(fù)共軛.

    在導(dǎo)納分析前,需要進(jìn)行相關(guān)性的計(jì)算,相關(guān)性的計(jì)算公式為(Kalnins and Watts, 2009):

    (2)

    其中Sff(k)表示重力異常的自功率譜密度.大的相關(guān)系數(shù)表示重力異常主要由地形引起,在進(jìn)行觀測(cè)導(dǎo)納與理論導(dǎo)納的對(duì)比時(shí),需要相關(guān)系數(shù)大于0.4的波數(shù)才采用(Kalnins and Watts, 2009).

    2.2 理論導(dǎo)納

    巖石圈在受到地表地形和內(nèi)部Moho面起伏的垂向加載后,經(jīng)過巖石圈撓曲變形調(diào)整后形成最終的地表地形和最終的Moho面起伏,根據(jù)地表地形和Moho面起伏可以計(jì)算產(chǎn)生的觀測(cè)到的自由空氣重力異常.下面我們首先計(jì)算巖石圈在受到垂向加載后產(chǎn)生的彈性撓曲變形,得到最終的海底地形和Moho面起伏,再計(jì)算它們產(chǎn)生的自由空氣重力異常,最后得到理論的自由空氣導(dǎo)納的計(jì)算公式.

    巖石圈在受到垂向加載q(X)后會(huì)產(chǎn)生撓曲,假設(shè)海洋巖石圈為一個(gè)完全彈性的板塊,其撓曲公式為(Turcotte and Schubert, 2002):

    (3)

    (4)

    其中ρc為地殼密度.對(duì)(4)式進(jìn)行傅里葉變換得到:D|k|4W(k)+(ρm-ρw)gW(k)=-(ρc-ρw)gH(k),

    (6)

    地形加載H、彈性撓曲W和最終地形T之間的關(guān)系為H=T-W,將其代入公式(6)并整理得到彈性撓曲和最終地形之間的關(guān)系為:

    (7)

    根據(jù)Parker(1973)給出的重力異常計(jì)算方法,Moho面撓曲W產(chǎn)生的布格重力異常為:ΔGB(k)=2πG(ρm-ρc)W(k)exp(-k(d+t)),

    (8)其中d,t分別為平均海水深度和平均地殼厚度.自由空氣重力異常為布格重力異常加上地表地形T產(chǎn)生的重力異常:

    ΔGF(k)=ΔGB(k)+2πG(ρc-ρw)T(k)exp(-kd).

    (9)

    綜合公式(7)、(8)、(9)可以得到理論的自由空氣導(dǎo)納Zpred(k)為:

    (10)

    定義

    (11)

    將公式(11)代入公式(10)得到:

    Zpred(k)= 2πG(ρc-ρw)exp(-kd)

    ×[1-Φe(k)exp(-kt)].

    (12)

    2.3 合成測(cè)試

    在計(jì)算觀測(cè)的自由空氣導(dǎo)納時(shí),需要求取海底地形和自由空氣重力異常的互功率譜密度以及海底地形的自功率譜密度.由于Multitaper方法比周期圖法能更精確地確定功率譜密度(Hanssen, 1997),本文使用Multitaper方法進(jìn)行功率譜密度估計(jì).

    對(duì)于巖石圈有效彈性厚度的反演,我們使用MWAT方法.首先選取某大小的一個(gè)滑動(dòng)窗口,根據(jù)公式(1)計(jì)算觀測(cè)導(dǎo)納,同時(shí)根據(jù)公式(12)計(jì)算不同有效彈性厚度和不同地殼密度時(shí)的理論導(dǎo)納,選取與觀測(cè)導(dǎo)納最接近的理論導(dǎo)納對(duì)應(yīng)的巖石圈有效彈性厚度Te作為反演的窗口中央的Te,同時(shí)也反演了地殼密度.在研究區(qū)域內(nèi)移動(dòng)該窗口,就可以獲得Te隨空間變化的情況.由于彈性撓曲的特征波長(zhǎng)和振幅與Te和加載都相關(guān),并沒有一個(gè)標(biāo)準(zhǔn)的反演巖石圈有效彈性厚度的最佳窗口.更小的窗口可以得到更高的空間分辨率,反演Te比較尖銳的變化,但會(huì)丟失長(zhǎng)波長(zhǎng)的撓曲信息;更大的窗口會(huì)得到更高精度的結(jié)果,但結(jié)果會(huì)偏向窗口內(nèi)的最大的撓曲特征,并減小了分辨率.本文選取了比較廣泛的窗口范圍,分別為400 km×400 km、600 km×600 km、800 km×800 km、1000 km×1000 km、1200 km×1200 km、1400 km×1400 km,每個(gè)窗口的權(quán)重為采用的有足夠相關(guān)性的離散波數(shù)的個(gè)數(shù),最終把六個(gè)窗口的結(jié)果加權(quán)平均作為反演的巖石圈有效彈性厚度.這種方法可以得到一個(gè)比較可靠的結(jié)果,同時(shí)還可以獲得較好的分辨率(Kalnins and Watts, 2009).由于MWAT方法使用有限窗口的觀測(cè)導(dǎo)納與無限大窗口的理論導(dǎo)納作對(duì)比,這會(huì)導(dǎo)致MWAT方法反演的巖石圈有效彈性厚度偏小(Pérez-Gussinyé et al., 2004; Kalnins and Watts, 2009).下面我們進(jìn)行MWAT方法的合成測(cè)試,計(jì)算該方法反演的結(jié)果與真實(shí)值的偏差,得到不同Te時(shí)的改正值.

    本測(cè)試基于實(shí)際的海底地形,選取數(shù)據(jù)范圍為10°N到24°N,140°E到154°E,假設(shè)該地形為初始荷載,再設(shè)定平均地殼厚度為10 km,在Te=10 km時(shí),根據(jù)彈性撓曲公式(6)計(jì)算得到巖石圈達(dá)到均衡時(shí)的海底地形(圖2a)和Moho面起伏(圖2b),再根據(jù)公式(9)計(jì)算產(chǎn)生的自由空氣重力異常(圖2c).圖2d給出了Te分別為5、10和20 km的理論導(dǎo)納以及整個(gè)研究區(qū)域的觀測(cè)導(dǎo)納.

    根據(jù)某一大小的Te生成海底地形和重力異常,就可以利用上面提到的方法反演整個(gè)區(qū)域的巖石圈有效彈性厚度.對(duì)不同大小的Te的測(cè)試結(jié)果顯示,當(dāng)窗口比較大的時(shí)候,反演的結(jié)果更接近真實(shí)值,當(dāng)窗口比較小的時(shí)候,反演的結(jié)果偏差更大,標(biāo)準(zhǔn)差也比較大.圖3顯示的是最終六個(gè)窗口加權(quán)平均的結(jié)果,在Te比較小的時(shí)候,偏差較小,總的來看,MWAT方法計(jì)算的結(jié)果大約偏小20%.

    3 結(jié)果

    圖4顯示了研究區(qū)域直接反演的Te和經(jīng)過2.3節(jié)的結(jié)果改正后的Te.由于相關(guān)系數(shù)需要大于0.4,灰色處表示數(shù)據(jù)不足而無法確定Te.圖4b的結(jié)果顯示在馬里亞納海溝附近的Te比較大,最大約為35 km;而卡羅琳板塊、菲律賓海板塊的Te非常小,東部西太平洋板塊的Te主要在8~15 km之間.馬里亞納海溝處大的Te表明俯沖的古老的太平洋巖石圈非常強(qiáng),與前人的研究也相一致(Kalnins and Watts, 2009; Zhang et al., 2014).馬里亞納海溝附近的結(jié)果也顯示了較大的窗口對(duì)Te的平滑作用,由于計(jì)算導(dǎo)納時(shí)窗口的作用,比較小的構(gòu)造的彈性撓曲特征反映并不明顯,導(dǎo)納主要反映了馬里亞納海溝附近的特征.前人用MWAT方法研究的夏威夷附近的巖石圈有效彈性厚度也顯示了類似的效果(Kalnins and Watts, 2009; 胡敏章等, 2015).這里反演的卡羅琳洋脊的Te大于20 km,但該處非常年輕,可能是由于最近的熱柱活動(dòng)才形成,前人的研究也顯示其Te比較小(Altis, 1999),較大的Te是由于窗口的影響,人為導(dǎo)致反演的卡羅琳洋脊、Sorol海槽的Te偏大.

    除了在馬里亞納海溝附近受到窗口平滑作用影響的卡羅琳洋脊和Sorol海槽這些地方外,卡羅琳板塊的有效彈性厚度都比較小,基本都小于15 km.特別是歐里皮克海隆和東卡羅琳盆地,它們的部分地區(qū)的Te小于5 km.歐里皮克海隆非常年輕,并且形成于比較弱的卡羅琳盆地之上,基本處于艾利均衡狀態(tài)(Altis, 1999),它們非常小的巖石圈有效彈性厚度是合理的.卡羅琳板塊南邊的New Guinea海溝的Te比較小,可能與其較小的俯沖巖石圈年齡有關(guān).Manus海溝的Te非常小的原因可能是由于此處處于非彈性均衡導(dǎo)致.卡羅琳板塊西邊的Ayu海槽的Te為10 km左右,符合其類似一個(gè)非常緩慢分裂開的大洋中脊的成因(Fujiwara et al., 1995).卡羅琳群島都比較年輕,其下部的太平洋巖石圈非常古老,如果由加載形成時(shí)的巖石圈的年齡來推斷,其Te應(yīng)該比較大,但這里結(jié)果顯示卡羅琳群島的Te在10 km左右.卡羅琳群島可能并不是熱點(diǎn)成因的海山鏈,如果它們的形成與俯沖構(gòu)造環(huán)境下的破裂相關(guān),那么Te將會(huì)減小(Rehman et al., 2013).

    圖2 Te =10 km時(shí)在彈性均衡時(shí)的地形(a)、Moho面深度(b)、自由空氣重力異常(c)和觀測(cè)導(dǎo)納以及理論導(dǎo)納(d)Fig.2 Bathymetry (a), Moho depth (b) and free-air gravity anomalies (c) for an elastic lithosphere of Te=10 km. Also shown is the observed admittance and theoretical admittance (d)

    在馬里亞納海溝附近地區(qū),結(jié)果顯示Te非常大,表明俯沖的太平洋板塊非常強(qiáng).在外圍隆起向海溝方向,Te并沒有減小,前人研究顯示的在此處由于破裂導(dǎo)致巖石圈有效彈性厚度減小(Zhang et al., 2014)在這里并不明顯.在海溝附近,為了得到更精確的Te值,可能需要考慮在模型中加入底部加載,并且去除表面沉積層以及熱冷卻的影響.

    圖3 反演的最終六個(gè)窗口加權(quán)平均的Te與輸入的Te之間的關(guān)系實(shí)線表示恢復(fù)的Te為輸入的80%.Fig.3 Input Te and mean recovered Te with error bars using a weighted average of all six window sizesThe solid line shows that the recovered Te is 80% of the input Te.

    我們的結(jié)果顯示菲律賓海板塊有一個(gè)非常低的Te,在很多地方小于5 km,甚至接近于0 km,同時(shí)還有部分的區(qū)域無法計(jì)算Te.重力異常和地形之間較低的相關(guān)性可能是由于此處存在底部加載的影響,或者該區(qū)域主要處于地殼均衡狀態(tài),彈性加載理論在這里不成立.菲律賓海板塊大部分地區(qū)非常低的Te與其他研究顯示菲律賓海板塊基本都處于地殼均衡狀態(tài)(Yen et al., 2015)的結(jié)論相一致.

    我們的結(jié)果顯示Magellan海山地區(qū)的Te在10~15 km左右,略微小于趙俐紅等(2010)的剖面研究的結(jié)果.Te在Magellan海山地區(qū)并不是單一的變化,很難用一個(gè)單一的熱點(diǎn)模型來解釋.與Kalnins和Watts(2009)的結(jié)果一致,在MWG地區(qū)Te從南到北逐漸減小,與加載形成時(shí)海底地殼的年齡逐漸減小相一致.Magellan海山和MWG都可以追蹤到形成于太平洋超級(jí)地幔柱附近(Smith et al., 1989),導(dǎo)致其Te較小.在東南角的Ontong-Java高原反演的Te也比較小,這與Kalnins和Watts(2009)的結(jié)果類似,與其可能形成在一個(gè)分裂中心之上相關(guān)(Gladczenko et al., 1997),Ontong-Java高原南邊較大的Te(>15 km)可能與其南邊古老的俯沖帶相關(guān).

    4 討論

    我們反演了卡羅琳板塊及其附近地區(qū)的巖石圈有效彈性厚度,研究區(qū)域包括了海溝、海山、海盆等構(gòu)造.不像大陸巖石圈的結(jié)構(gòu)那么復(fù)雜,海洋巖石圈的熱結(jié)構(gòu)可以用一個(gè)板塊冷卻模型來近似表示(Stein C A and Stein S, 1992),而巖石圈的強(qiáng)度與其熱結(jié)構(gòu)密切相關(guān),隨著離洋中脊距離的增加而增強(qiáng).大量的研究結(jié)果都顯示海洋巖石圈的Te隨著加載時(shí)的巖石圈年齡的增加而增加,Te主要分布在板塊冷卻模型的200~600 ℃的等溫線深度之間(Kalnins and Watts, 2009;Watts, 1978; 胡敏章等, 2015).但最近的一些研究顯示,無論是在海山(Watts et al., 2006),還是海溝(Bry and White,2007)地區(qū),Te與加載時(shí)的巖石圈年齡之間并沒有一個(gè)簡(jiǎn)單的對(duì)應(yīng)關(guān)系.

    圖4 研究區(qū)域直接反演的巖石圈有效彈性厚度Te(a)和經(jīng)過2.3節(jié)的結(jié)果改正后的巖石圈有效彈性厚度Te(b)灰色處表示沒有足夠的相關(guān)性來確定Te.Fig.4 Raw recovered Te (a) and recovered Te after the bias correction from Section 2.3 is applied for the study areaGrey areas in the map indicate regions with insufficient coherence to determine Te.

    圖5顯示了本研究區(qū)域的海底年齡(Müller et al., 2008),對(duì)比圖4可以看出,在很多地方巖石圈的有效彈性厚度與海底年齡之間存在一定的相關(guān)性.菲律賓海板塊和卡羅琳板塊相對(duì)更年輕,其Te比較小,西太平洋的部分地區(qū)更古老,其Te也相對(duì)更大.在馬里亞納海溝有非常大的Te,由于計(jì)算方法的影響,導(dǎo)致其附近區(qū)域計(jì)算的Te都比較大.在海溝附近Te與俯沖巖石圈的年齡之間也存在一定的相關(guān)性,馬里亞納海溝附近的俯沖巖石圈年齡非常大,其Te也最大;其次是菲律賓海溝附近的俯沖巖石圈,其Te為20 km左右.卡羅琳板塊周邊的New guinea海溝、Manus海溝、Yap海溝和Mussau海溝的俯沖巖石圈的年齡都比較小,除Yap海溝以外,它們的Te也都比較小.Yap海溝的Te可能受到計(jì)算方法的影響,需要進(jìn)一步的研究.

    圖5 研究區(qū)域的海底年齡和已知年齡的海山的位置Fig.5 Seafloor ages of the study area and locations of dated seamounts

    圖6a顯示了本研究區(qū)域所有計(jì)算的Te與海底年齡之間的關(guān)系,除了馬里亞納海溝附近地區(qū)部分年輕的海洋巖石圈有較大的Te,數(shù)據(jù)主要都位于板塊冷卻模型的450 ℃的等溫線深度以上.這是因?yàn)楹I叫纬蓵r(shí)的年齡總會(huì)小于海底地殼的年齡,較早的海山加載的巖石圈的Te會(huì)小于最近形成的海山加載的巖石圈的Te.Watts等(2006)對(duì)Te與海底地殼年齡之間關(guān)系的研究也顯示了類似的特征.

    要知道加載形成時(shí)巖石圈的年齡,在已知海底年齡的情況下(Müller et al., 2008),還需要知道海山的年齡.最近這些年對(duì)海山年齡的取樣分析的數(shù)據(jù)有了很大的推進(jìn),本文收集了研究區(qū)域20個(gè)海山的年齡數(shù)據(jù)(Clouard and Bonneville, 2005).圖6b顯示了巖石圈有效彈性厚度與海山加載時(shí)巖石圈年齡之間的關(guān)系,巖石圈的有效彈性厚度主要分布在板塊冷卻模型的200 ℃的等溫線深度附近,小于前人在其他許多地方的研究結(jié)果,Te也并未隨海山加載時(shí)巖石圈年齡的增加而增大.加載形成時(shí)的巖石圈年齡并不是影響Te的唯一因素,本研究的海山的位置主要位于Magellan海山和MWG地區(qū)(圖5),它們都可以追蹤到形成于太平洋超級(jí)地幔柱附近(Smith et al., 1989; 趙俐紅等, 2010),較低的等溫線深度可能是受到太平洋超級(jí)地幔柱的影響,在超級(jí)地幔柱之上形成的巖石圈可能更薄并且更弱(Kalnins and Watts, 2009).另一個(gè)可能的原因就是多次噴發(fā)引起海山測(cè)量年齡的不準(zhǔn)確,多次噴發(fā)的海山的測(cè)量年齡會(huì)偏小,從而引起計(jì)算的加載時(shí)的巖石圈年齡偏大,導(dǎo)致一個(gè)較低的等溫線深度.

    圖6 巖石圈有效彈性厚度與海底年齡之間的關(guān)系(a)和巖石圈有效彈性厚度與海山加載時(shí)巖石圈年齡之間的關(guān)系(b)等溫線基于板塊厚度為95 km的板塊冷卻模型,熱擴(kuò)散系數(shù)為10-6 m·s-2,上下邊界的溫度分別為0 ℃和1300 ℃.Fig.6 Corrected Te versus ages of the oceanic crust (a) and corrected Te versus lithosphere ages at the time of loading (b)Isotherms are based on a cooling plate model with a plate thickness of 95 km, thermal diffusivity of 10-6 m·s-2, and temperatures of 0 ℃ and 1300 ℃ at the upper and lower plate boundaries, respectively.

    地表熱流包含了熱傳導(dǎo)系數(shù)、生熱率、地溫梯度和來自地球內(nèi)部的熱流等重要地球物理信息.對(duì)于大陸巖石圈的研究顯示小的Te與大的地表熱流相關(guān),大的Te對(duì)應(yīng)于較低的地表熱流(楊亭等, 2013; Deng et al., 2014).本文收集了研究區(qū)域的大約600個(gè)地表熱流數(shù)據(jù)(http:∥www.heatflow.und.edu/index2.html).圖7顯示了巖石圈的有效彈性厚度Te與熱流之間的關(guān)系,可以看出在研究區(qū)域內(nèi)海洋地殼的熱流與Te之間存在一定的反相關(guān)性,但當(dāng)Te比較小時(shí)這樣的關(guān)系并不明顯.熱流主要與海洋巖石圈的年齡相關(guān),而Te主要與加載時(shí)的海洋巖石圈的年齡相關(guān),其與海洋巖石圈的年齡之間并不是簡(jiǎn)單的對(duì)應(yīng)關(guān)系(圖6).對(duì)于海洋地區(qū)Te與海底熱流之間的對(duì)應(yīng)關(guān)系,仍然需要更多的研究.

    圖7 巖石圈有效彈性厚度Te與熱流之間的關(guān)系實(shí)心圓表示5 km內(nèi)的平均Te和平均的熱流.Fig.7 Corrected Te versus heat flowSolid circles represent the mean Tewithin 5 km interval versus mean heat flow.

    5 結(jié)論

    本文利用卡羅琳板塊及其附近地區(qū)的自由空氣重力異常和海水深度數(shù)據(jù),結(jié)合滑動(dòng)窗口導(dǎo)納技術(shù),計(jì)算了該地區(qū)的巖石圈有效彈性厚度Te.主要結(jié)論如下:

    (1) 本文使用Multitaper方法對(duì)功率譜密度進(jìn)行估計(jì),基于實(shí)際的海底地形,通過模擬計(jì)算得到了MWAT方法較真值的改正, MWAT方法計(jì)算的結(jié)果大約偏小20%.

    (2) 研究結(jié)果顯示卡羅琳板塊的大部分地區(qū)的Te值都比較低,在歐里皮克海隆的部分地區(qū)Te小于5 km;馬里亞納海溝附近的Te值比較大,最大約為35 km;菲律賓海板塊的Te值也非常低,在很多地方小于5 km;古老的西太平洋板塊地區(qū)的Te值在8~15 km之間.

    (3) 反演的卡羅琳洋脊、Yap海溝和Sorol海槽的Te比較大,主要是在馬里亞納海溝附近受窗口平滑作用影響的原因,由于計(jì)算導(dǎo)納時(shí)窗口的作用,比較小的構(gòu)造的彈性撓曲特征的反映并不明顯,導(dǎo)致了這些地方計(jì)算的Te偏大.

    (4) 反演結(jié)果顯示歐里皮克海隆、菲律賓海板塊局部地區(qū)有著非常低的Te,甚至接近0 km,可能與這些區(qū)域主要處于地殼均衡狀態(tài)相關(guān).

    (5)Te與海底地殼年齡之間的關(guān)系顯示Te主要位于板塊冷卻模型的450 ℃的等溫線深度以上,與前人的結(jié)果類似.西太平洋的Magellan海山和MWG地區(qū)的Te主要分布在加載形成時(shí)板塊冷卻模型的200 ℃的等溫線深度附近,小于前人在其他許多地方的研究結(jié)果.這里較低的等溫線深度可能是受到太平洋超級(jí)地幔柱的影響,在超級(jí)地幔柱之上形成的巖石圈可能更薄并且更弱,另外可能的原因包括多次噴發(fā)引起海山測(cè)量年齡的不準(zhǔn)確.

    (6) 在研究區(qū)域內(nèi)海洋地殼的熱流與Te之間存在一定的反相關(guān)性,但當(dāng)Te比較小時(shí)這樣的關(guān)系并不明顯.

    對(duì)于馬里亞納海溝附近地區(qū),可能存在底部加載的情況,進(jìn)一步研究的時(shí)候需要考慮這個(gè)問題.如果有馬里亞納海溝附近的Yap海溝和卡羅琳洋脊的典型的二維剖面數(shù)據(jù),可以通過剖面的研究減小MWAT方法造成的計(jì)算結(jié)果失真的問題.

    致謝 感謝美國加州理工學(xué)院的楊亭博士提供了多窗譜計(jì)算程序以及在本文寫作過程中的幫助.

    Altis S. 1999. Origin and tectonic evolution of the Caroline Ridge and the Sorol Trough, western tropical Pacific, from admittance and a tectonic modeling analysis.Tectonophysics, 313(3): 271-292.

    Bird P. 2003. An updated digital model of plate boundaries.Geochemistry,Geophysics,Geosystems, 4(3): 1027.

    Bry M, White N. 2007. Reappraising elastic thickness variation at oceanic trenches.J.Geophys.Res., 112: B08414.

    Calmant S, Francheteau J, Cazenave A. 1990. Elastic layer thickening with age of the oceanic lithosphere: a tool for prediction of the age of volcanoes or oceanic crust.Geophys.J.Int., 100(1): 59-67.

    Clouard V, Bonneville A. 2005. Ages of seamounts, islands, and plateaus on the Pacific plate.GSASpecialPapers, 388: 71-90.

    Deng Y F, Zhang Z J, Fan W M, et al. 2014. Multitaper spectral method to estimate the elastic thickness of South China: Implications for intracontinental deformation.GeoscienceFrontiers, 5(2): 193-203. Filmer P E, Mcnutt M K, Wolfe C J. 1993. Elastic thickness of the lithosphere in the Marquesas and Society islands.J.Geophys.Res., 98(B11): 19565-19577.

    Forsyth D W. 1985. Subsurface loading and estimates of the flexural rigidity of continental lithosphere.J.Geophys.Res., 90(B14): 12623-12632.

    Fu Y T, Li A C, Qin Y S. 2002. Effective elastic thickness of the oceanic and continental marginal lithospheres.MarineGeology&QuaternaryGeology(in Chinese), 22(3): 69-75.

    Fujiwara T, Tamaki K, Fujimoto H, et al. 1995. Morphological-studies of the Ayu trough, Philippine-sea Caroline plate boundary.Geophys.Res.Lett., 22(2): 109-112.

    Gladczenko T P, Coffin, M F, Eldholm O. 1997. Crustal structure of the Ontong Java Plateau: Modeling of new gravity and existing seismic data.J.Geophys.Res., 102(B10): 22711-22729.

    Hanssen A. 1997. Multidimensional multitaper spectral estimation.SignalProcessing, 58(3): 327-332.

    Hill K C, Hegarty K A. 1987. New tectonic framework for PNG and the Caroline plate: Implications for cessation of spreading in back-arc basins.∥ Pacific Rim Congress 87, The Australian Institute of Mining and Metallurgy, 179-182.

    Hu M Z, Li J C, Li H, et al. 2015. The lithosphere effective elastic thickness and its tectonic implications in the Northwestern Pacific.ChineseJ.Geophys. (in Chinese), 58(2): 542-555, doi: 10.6038/cjg20150217.

    Judge A V, Mcnutt M K. 1991. The relationship between plate curvature and elastic plate thickness: A study of the Peru-Chile trench.J.Geophys.Res., 96(B10): 16625-16639.

    Kalnins L M, Watts A B. 2009. Spatial variations in effective elastic thickness in the Western Pacific Ocean and their implications for Mesozoic volcanism.EarthandPlanetaryScienceLetters, 286(1-2): 89-100.

    Keating B H, Mattey D P, Helsley C E, et al. 1984. Evidence for a hot spot origin of the Caroline islands.J.Geophys.Res., 89(B12): 9937-9948.

    Lyons S N, Sandwell D T, Smith W H F. 2000. Three-dimensional estimation of elastic thickness under the Louisville Ridge.J.Geophys.Res., 105(B6): 13239-13252.

    Müller R D, Sdrolias M, Gaina C, et al. 2008. Age, spreading rates, and spreading asymmetry of the world's ocean crust.Geochemistry,Geophysics,Geosystems, 9(4): Q04006.

    Pérez-Gussinyé M, Lowry A R, Watts A B, et al. 2004. On the recovery of effective elastic thickness using spectral methods: Examples from synthetic data and from the Fennoscandian Shield.J.Geophys.Res., 109: B10409.

    Parker R L. 1973. The rapid calculation of potential anomalies.Geophys.J.Int., 31(4): 447-455.

    Rehman H U, Nakaya H, Kawai K. 2013. Geological origin of the volcanic islands of the Caroline group in the federated states of Micronesia, western Pacific.SouthPacificStudies, 33(2): 101-118.

    Sandwell D T, Smith W H F. 2009. Global marine gravity from retracked Geosat and ERS-1 altimetry: Ridge segmentation versus spreading rate.J.Geophys.Res., 114: B01411.

    Sato T, Kasahara J, Katao H, et al. 1997. Seismic observations at the Yap Islands and the northern Yap Trench.Tectonophysics, 271(3-4): 285-294.

    Smith W H F, Staudigel H, Watts A B, et al. 1989. The Magellan seamounts: Early Cretaceous record of the South Pacific isotopic and thermal anomaly.J.Geophys.Res., 94(B8): 10501-10523.

    Stein C A, Stein S. 1992. A model for the global variation in oceanic depth and heat flow with lithospheric age.Nature, 359(6391): 123-129.

    Su D Q. 2012. A study of the effective elastic thickness of the oceanic lithosphere.ChineseJ.Geophys. (in Chinese), 55(10): 3259-3265, doi: 10.6038/j.issn.0001-5733.2012.10.008. Turcotte D L, Schubert G. 2002. Geodynamics. Cambridge: Cambridge University Press. Watts A B. 1978. An analysis of isostasy in the world's oceans 1. Hawaiian-Emperor seamount chain.J.Geophys.Res., 83(B12): 5989-6004.

    Watts A B, tem Brind U S, Buhl P, et al. 1985. A multichannel seismic study of lithospheric flexure across the Hawaiian-Emperor seamount chain.Nature, 315(6015): 105-111.

    Watts A B, Sandwell D T, Smith W H F, et al. 2006. Global gravity, bathymetry, and the distribution of submarine volcanism through space and time.J.Geophys.Res., 111: B08408. Weissel J K, Anderson R N. 1978. Is there a Caroline plate?EarthandPlanetaryScienceLetters, 41(2): 143-158.

    Yang T, Fu R S, Huang J S. 2013. Effective elastic thickness of continental lithosphere in China with Moho topography admittance method.ChineseJ.Geophys. (in Chinese), 56(6): 1877-1886, doi: 10.6038/cjg20130610.

    Yen H Y, Lo Y T, Yeh Y L, et al. 2015. The crustal thickness of the Philippine sea plate derived from gravity data.Terr.Atmos.Ocean.Sci., 26(3): 253-259.

    Zhang F, Lin J, Zhan W H. 2014. Variations in oceanic plate bending along the Mariana trench.EarthandPlanetaryScienceLetters, 401: 206-214.

    Zhao L H, Jin X L, Gao J Y, et al. 2010. The effective elastic thickness of lithosphere in the mid-west Pacific and its geological significance.EarthScience(JournalofChinaUniversityofGeosciences) (in Chinese), 35(4): 637-644.

    Zhong S J, Watts A B. 2013. Lithospheric deformation induced by loading of the Hawaiian Islands and its implications for mantle rheology.J.Geophys.Res., 118(11): 6025-6048.

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

    付永濤, 李安春, 秦蘊(yùn)珊. 2002. 大洋和大陸邊緣巖石圈有效彈性厚度的研究意義. 海洋地質(zhì)與第四紀(jì)地質(zhì), 22(3): 69-75.

    胡敏章, 李建成, 李輝等. 2015. 西北太平洋巖石圈有效彈性厚度及其構(gòu)造意義. 地球物理學(xué)報(bào), 58(2): 542-555, doi: 10.6038/cjg20150217.

    蘇達(dá)權(quán). 2012. 海洋巖石圈板塊有效彈性厚度研究. 地球物理學(xué)報(bào), 55(10): 3259-3265, doi: 10.6038/j.issn.0001-5733.2012.10.008.

    楊亭, 傅容珊, 黃金水. 2013. 利用Moho地形導(dǎo)納法(MDDF)反演中國大陸巖石圈有效彈性厚度. 地球物理學(xué)報(bào), 56(6): 1877-1886, doi: 10.6038/cjg20130610.

    趙俐紅, 金翔龍, 高金耀等. 2010. 中西太平洋海山區(qū)的巖石圈有效彈性厚度及其地質(zhì)意義. 地球科學(xué)(中國地質(zhì)大學(xué)學(xué)報(bào)), 35(4): 637-644.

    (本文編輯 何燕)

    The effective elastic thickness of the Caroline plate and its adjacent areas

    YANG An1,2, FU Yong-Tao1,2, LI An-Chun1,2

    1InstituteofOceanology,ChineseAcademyofSciences,Qingdao266071,China2KeyLaboratoryofMarineGeologyandEnvironment,ChineseAcademyofSciences,Qingdao266071,China

    The purpose of this work is to study the effective elastic thickness (Te) of the Caroline plate and its adjacent areas by using free-air gravity anomaly and bathymetric data, in combination with a moving window admittance technique (MWAT). We use the multitaper spectral estimation method to calculate the power spectral density. Correction is obtained through synthetic tests with actual bathymetric data. Results show that the recoveredTeis about 20% lower than the inputTe. The Caroline plate and adjacent areas have aTerange of 1~34 km. The study area includes a wide range of features such as seamounts, plateaux, subduction zones and spreading ridges. The examination of the lithospheric strength for these structures provides important information for the tectonic setting and evolution of the Western Pacific lithosphere.Teis related to the age of the lithosphere at the time of loading, heat flow and other proxy features. Plots ofTeversus ages of the oceanic crust show thatTeis above the depth of 450 ℃ isotherm based on a cooling plate model. For the Magellan seamounts and Marcus-Wake Guyots (MWG) areas, plots ofTeversus the ages of the lithosphere at the time of loading suggest thatTeis at around the depth of 200 ℃ isotherm based on a cooling plate model. The low isotherm could resulted from the effect of the South Pacific supper plume. Our results also show a negative correlation between heat flow andTein the study area.

    Effective elastic thickness; Flexure; Caroline plate; Lithosphere; Moving window admittance technique

    10.6038/cjg20160913.

    中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)(XDA1103010102)和中國科學(xué)院海洋地質(zhì)與環(huán)境重點(diǎn)實(shí)驗(yàn)室開放基金(MGE2015KG07)資助.

    楊安,男,1987年生,博士,研究方向:地幔對(duì)流和巖石圈動(dòng)力學(xué).E-mail:herryangan@163.com

    10.6038/cjg20160913

    P541,P738

    2015-08-02,2016-07-27收修定稿

    楊安, 付永濤, 李安春. 2016. 卡羅琳板塊及其附近地區(qū)的巖石圈有效彈性厚度. 地球物理學(xué)報(bào),59(9):3280-3290,

    Yang A, Fu Y T, Li A C. 2016. The effective elastic thickness of the Caroline plate and its adjacent areas.ChineseJ.Geophys. (in Chinese),59(9):3280-3290,doi:10.6038/cjg20160913.

    猜你喜歡
    馬里亞納海山卡羅琳
    烏克蘭裔“日本小姐”在日引巨大爭(zhēng)議
    深海魅影:馬里亞納海溝里的奇怪生物
    馬里亞納海溝的奇怪生物
    軍事文摘(2021年20期)2021-11-10 01:58:54
    On the heating mechanism of electron cyclotron resonance thruster immerged in a non-uniform magnetic field
    軍器圖說 馬里亞納海戰(zhàn)
    夢(mèng)里鮮花開放
    37°女人(2019年4期)2019-04-20 02:17:42
    我排第一個(gè)
    我排第一個(gè)
    佐藤海山的詩
    作品(2016年10期)2016-12-06 10:41:54
    迎春花
    大江南北(2016年3期)2016-11-22 07:37:20
    三级国产精品片| 国产 一区 欧美 日韩| 日本一二三区视频观看| 少妇裸体淫交视频免费看高清| 久久久久久久久大av| 天天一区二区日本电影三级| 亚洲丝袜综合中文字幕| 嘟嘟电影网在线观看| 青春草国产在线视频| av女优亚洲男人天堂| 亚洲精品久久午夜乱码| 爱豆传媒免费全集在线观看| 交换朋友夫妻互换小说| 久久久久国产网址| 国产毛片在线视频| 午夜精品一区二区三区免费看| 久久精品国产亚洲av天美| 国产乱来视频区| 嫩草影院入口| 我要看日韩黄色一级片| av一本久久久久| 又爽又黄a免费视频| 香蕉精品网在线| 狠狠精品人妻久久久久久综合| 成人特级av手机在线观看| 亚洲精品国产色婷婷电影| 日韩成人av中文字幕在线观看| 美女xxoo啪啪120秒动态图| 欧美性猛交╳xxx乱大交人| 嫩草影院入口| 王馨瑶露胸无遮挡在线观看| 国产在视频线精品| 色视频www国产| 国产91av在线免费观看| 久久久久久久久久人人人人人人| 久久久成人免费电影| 一个人观看的视频www高清免费观看| 国产亚洲一区二区精品| 日本熟妇午夜| 精品一区二区免费观看| 听说在线观看完整版免费高清| 一级毛片我不卡| 中文字幕亚洲精品专区| 精品少妇黑人巨大在线播放| 免费观看a级毛片全部| 嫩草影院精品99| 免费黄色在线免费观看| 国产探花极品一区二区| 少妇丰满av| 久久6这里有精品| 男人狂女人下面高潮的视频| 99久久九九国产精品国产免费| 日产精品乱码卡一卡2卡三| 一级毛片电影观看| 超碰97精品在线观看| 欧美日本视频| 国产 一区精品| 久久人人爽av亚洲精品天堂 | 51国产日韩欧美| 中文字幕免费在线视频6| 伊人久久精品亚洲午夜| 在线观看一区二区三区| 日日摸夜夜添夜夜爱| 亚洲人成网站在线播| 久久久久国产网址| 美女国产视频在线观看| 2021少妇久久久久久久久久久| 国产一区有黄有色的免费视频| 99视频精品全部免费 在线| 国产伦精品一区二区三区视频9| 免费观看无遮挡的男女| 老女人水多毛片| 一级毛片黄色毛片免费观看视频| 爱豆传媒免费全集在线观看| 97精品久久久久久久久久精品| 小蜜桃在线观看免费完整版高清| 91aial.com中文字幕在线观看| 日韩成人av中文字幕在线观看| 91精品国产九色| 日韩,欧美,国产一区二区三区| 免费观看a级毛片全部| 一区二区三区四区激情视频| 亚洲天堂国产精品一区在线| 久久热精品热| 国产高清有码在线观看视频| 亚洲天堂国产精品一区在线| 交换朋友夫妻互换小说| 综合色丁香网| 久久久久久久国产电影| 三级男女做爰猛烈吃奶摸视频| 99热网站在线观看| 国产老妇伦熟女老妇高清| 秋霞伦理黄片| 五月伊人婷婷丁香| 国内精品宾馆在线| 国产毛片在线视频| 亚洲伊人久久精品综合| 欧美日韩国产mv在线观看视频 | 欧美三级亚洲精品| 欧美另类一区| 美女视频免费永久观看网站| 国产熟女欧美一区二区| av一本久久久久| 韩国高清视频一区二区三区| 色综合色国产| 日韩av免费高清视频| 狂野欧美激情性bbbbbb| 黄色一级大片看看| 日日撸夜夜添| 免费黄色在线免费观看| 免费大片黄手机在线观看| 亚洲精品成人久久久久久| 69av精品久久久久久| 久久精品国产鲁丝片午夜精品| 日本免费在线观看一区| 国产黄a三级三级三级人| 一级黄片播放器| 久久久久久久久久成人| 美女国产视频在线观看| 亚洲人与动物交配视频| 国产成人freesex在线| 一区二区av电影网| 亚洲国产精品成人久久小说| 国产探花极品一区二区| av免费在线看不卡| 亚洲精品第二区| 日韩国内少妇激情av| 少妇人妻 视频| 久久国内精品自在自线图片| 欧美变态另类bdsm刘玥| 又大又黄又爽视频免费| 黄色视频在线播放观看不卡| 日韩国内少妇激情av| 九九久久精品国产亚洲av麻豆| 久久久久九九精品影院| 国产精品不卡视频一区二区| 亚洲av男天堂| 人妻系列 视频| 国产毛片a区久久久久| 最近最新中文字幕免费大全7| 老司机影院成人| 免费av观看视频| 亚洲天堂av无毛| av国产久精品久网站免费入址| 男女那种视频在线观看| 2021天堂中文幕一二区在线观| 天美传媒精品一区二区| 真实男女啪啪啪动态图| 久久99蜜桃精品久久| 女人十人毛片免费观看3o分钟| 国产在线一区二区三区精| 午夜福利在线观看免费完整高清在| 大片电影免费在线观看免费| 欧美 日韩 精品 国产| 成人亚洲欧美一区二区av| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲国产精品成人久久小说| 日韩,欧美,国产一区二区三区| 国产成人91sexporn| 禁无遮挡网站| 不卡视频在线观看欧美| 高清av免费在线| 性色av一级| 老司机影院毛片| 国产成人福利小说| 亚洲久久久久久中文字幕| 在线天堂最新版资源| 欧美另类一区| 涩涩av久久男人的天堂| 日本与韩国留学比较| 亚洲婷婷狠狠爱综合网| 日日摸夜夜添夜夜添av毛片| 国产精品一及| 亚洲精品中文字幕在线视频 | 亚洲人成网站高清观看| 亚洲精品成人久久久久久| 一级毛片我不卡| 赤兔流量卡办理| eeuss影院久久| 大香蕉97超碰在线| 国产毛片a区久久久久| av又黄又爽大尺度在线免费看| 免费黄色在线免费观看| 国产精品久久久久久久电影| 久久精品久久精品一区二区三区| av专区在线播放| 五月伊人婷婷丁香| 国产精品一区www在线观看| 日韩国内少妇激情av| 一本色道久久久久久精品综合| 国精品久久久久久国模美| 免费黄频网站在线观看国产| 亚洲精品成人久久久久久| 2021少妇久久久久久久久久久| 日韩在线高清观看一区二区三区| av又黄又爽大尺度在线免费看| 性色av一级| 精品国产露脸久久av麻豆| 国产亚洲午夜精品一区二区久久 | 亚洲一级一片aⅴ在线观看| 人妻系列 视频| 在线a可以看的网站| 成人特级av手机在线观看| 免费黄网站久久成人精品| 三级国产精品欧美在线观看| 一个人观看的视频www高清免费观看| 成年女人在线观看亚洲视频 | 精品久久久噜噜| 欧美日韩视频精品一区| 最近2019中文字幕mv第一页| 插逼视频在线观看| 69av精品久久久久久| 精品视频人人做人人爽| 日本午夜av视频| 国产精品熟女久久久久浪| 在线精品无人区一区二区三 | 久久久久精品性色| 欧美区成人在线视频| 日韩强制内射视频| 亚洲人成网站高清观看| 99热国产这里只有精品6| 久久久久网色| 美女脱内裤让男人舔精品视频| 国产视频首页在线观看| 亚洲人成网站高清观看| 国产精品秋霞免费鲁丝片| 97精品久久久久久久久久精品| 亚洲成色77777| 国产亚洲最大av| av在线app专区| 免费看a级黄色片| 欧美xxxx黑人xx丫x性爽| 亚洲欧美成人精品一区二区| 少妇的逼好多水| 内地一区二区视频在线| 精品人妻视频免费看| 成年人午夜在线观看视频| 少妇高潮的动态图| 成人漫画全彩无遮挡| 国产精品麻豆人妻色哟哟久久| 国产 一区精品| 在线 av 中文字幕| av福利片在线观看| 亚洲,一卡二卡三卡| 亚洲丝袜综合中文字幕| 日韩三级伦理在线观看| 伊人久久国产一区二区| 26uuu在线亚洲综合色| 简卡轻食公司| 一级毛片黄色毛片免费观看视频| 男男h啪啪无遮挡| 亚洲精华国产精华液的使用体验| 大片电影免费在线观看免费| 免费黄频网站在线观看国产| 99久久精品热视频| 狠狠精品人妻久久久久久综合| 91午夜精品亚洲一区二区三区| 欧美日韩国产mv在线观看视频 | 婷婷色综合大香蕉| 日本黄色片子视频| 亚洲精品日韩在线中文字幕| 午夜爱爱视频在线播放| 肉色欧美久久久久久久蜜桃 | 亚洲欧美精品自产自拍| 能在线免费看毛片的网站| 搞女人的毛片| 久久影院123| 国产一区有黄有色的免费视频| 看非洲黑人一级黄片| 亚洲精品第二区| 肉色欧美久久久久久久蜜桃 | 成人毛片60女人毛片免费| av卡一久久| 人人妻人人爽人人添夜夜欢视频 | 国产黄色免费在线视频| 麻豆成人午夜福利视频| 永久网站在线| 春色校园在线视频观看| 亚洲成色77777| 亚洲av一区综合| 日韩欧美一区视频在线观看 | 免费播放大片免费观看视频在线观看| 亚洲精品456在线播放app| 中文字幕av成人在线电影| 一本久久精品| 国产老妇女一区| av国产精品久久久久影院| 少妇的逼水好多| 精品午夜福利在线看| 亚洲图色成人| 搡女人真爽免费视频火全软件| 韩国高清视频一区二区三区| 国产精品国产三级专区第一集| 在线 av 中文字幕| 日本黄大片高清| 97热精品久久久久久| 看免费成人av毛片| 一区二区三区四区激情视频| 国产探花极品一区二区| 校园人妻丝袜中文字幕| 日韩成人伦理影院| 成人特级av手机在线观看| 天堂网av新在线| 亚洲综合色惰| av在线app专区| 久久精品久久久久久噜噜老黄| 成年免费大片在线观看| 久久久午夜欧美精品| 好男人在线观看高清免费视频| 99精国产麻豆久久婷婷| 亚洲国产最新在线播放| 99精国产麻豆久久婷婷| 丰满少妇做爰视频| 精品国产露脸久久av麻豆| 一个人观看的视频www高清免费观看| 高清视频免费观看一区二区| 亚洲av一区综合| kizo精华| 国国产精品蜜臀av免费| 国产欧美日韩精品一区二区| 国产亚洲91精品色在线| 男女无遮挡免费网站观看| 亚洲精品乱久久久久久| 久久久久久久亚洲中文字幕| 大码成人一级视频| 国产成人精品一,二区| 国产精品久久久久久精品电影| 亚洲人成网站在线播| av免费观看日本| 久久99热6这里只有精品| 日韩人妻高清精品专区| 97精品久久久久久久久久精品| 欧美精品一区二区大全| 麻豆成人午夜福利视频| 亚洲精品aⅴ在线观看| 日韩免费高清中文字幕av| 伦理电影大哥的女人| 免费观看性生交大片5| 在线精品无人区一区二区三 | 精华霜和精华液先用哪个| 精品国产一区二区三区久久久樱花 | 免费观看无遮挡的男女| 国产乱人偷精品视频| 日本-黄色视频高清免费观看| 国产男女内射视频| 丝袜脚勾引网站| av.在线天堂| 97热精品久久久久久| 99久久九九国产精品国产免费| 91在线精品国自产拍蜜月| 色5月婷婷丁香| 听说在线观看完整版免费高清| 欧美变态另类bdsm刘玥| 夫妻午夜视频| 97超碰精品成人国产| 久久ye,这里只有精品| 男人添女人高潮全过程视频| 精品久久国产蜜桃| 97超碰精品成人国产| 国产精品爽爽va在线观看网站| 另类亚洲欧美激情| 精品久久国产蜜桃| 人妻 亚洲 视频| 久久久久久国产a免费观看| 日韩视频在线欧美| 久久久久久久午夜电影| 日韩成人伦理影院| 最近中文字幕高清免费大全6| 国产精品99久久久久久久久| 女人十人毛片免费观看3o分钟| 欧美激情国产日韩精品一区| 国产高清有码在线观看视频| 亚洲经典国产精华液单| 亚洲欧美清纯卡通| 国产女主播在线喷水免费视频网站| 亚洲精品日本国产第一区| 成人毛片60女人毛片免费| 精品一区二区三区视频在线| 亚洲精品自拍成人| 超碰97精品在线观看| 少妇人妻一区二区三区视频| 18禁动态无遮挡网站| 精品国产乱码久久久久久小说| 亚洲av不卡在线观看| 久久久久精品久久久久真实原创| av在线app专区| 18禁裸乳无遮挡动漫免费视频 | 精品一区在线观看国产| 亚洲,欧美,日韩| 久久久精品94久久精品| 亚洲,欧美,日韩| 免费观看性生交大片5| 国产成人精品福利久久| 一级毛片我不卡| 成年女人看的毛片在线观看| 下体分泌物呈黄色| 久久99精品国语久久久| 成人亚洲精品av一区二区| 国内揄拍国产精品人妻在线| 伊人久久国产一区二区| 亚洲欧美日韩无卡精品| 最后的刺客免费高清国语| 韩国高清视频一区二区三区| av线在线观看网站| 边亲边吃奶的免费视频| 男男h啪啪无遮挡| 国产精品国产av在线观看| 亚洲av成人精品一二三区| 日韩一区二区三区影片| 国产午夜福利久久久久久| 大片免费播放器 马上看| 另类亚洲欧美激情| 丰满人妻一区二区三区视频av| 国产精品蜜桃在线观看| 国产永久视频网站| 春色校园在线视频观看| 国产精品.久久久| 久久综合国产亚洲精品| 男女无遮挡免费网站观看| 日日撸夜夜添| 亚洲人成网站在线播| av播播在线观看一区| 国产美女午夜福利| 日本-黄色视频高清免费观看| 国产成人a∨麻豆精品| 特级一级黄色大片| 精品一区在线观看国产| 欧美一级a爱片免费观看看| 国产精品秋霞免费鲁丝片| 一级av片app| 亚洲精品久久久久久婷婷小说| 五月天丁香电影| 欧美区成人在线视频| 一二三四中文在线观看免费高清| 亚洲国产精品999| 精品久久国产蜜桃| 一级爰片在线观看| 日日撸夜夜添| 直男gayav资源| 看非洲黑人一级黄片| 又爽又黄无遮挡网站| 国产永久视频网站| 日本欧美国产在线视频| 最近中文字幕2019免费版| 国产免费视频播放在线视频| 欧美激情国产日韩精品一区| 狂野欧美激情性xxxx在线观看| 免费人成在线观看视频色| 黑人高潮一二区| 伊人久久国产一区二区| 少妇丰满av| 欧美 日韩 精品 国产| 黄片无遮挡物在线观看| 亚洲欧美一区二区三区国产| 亚洲丝袜综合中文字幕| 亚洲av男天堂| 九九爱精品视频在线观看| 亚洲精品成人久久久久久| 女人被狂操c到高潮| 欧美亚洲 丝袜 人妻 在线| 免费黄色在线免费观看| 欧美精品人与动牲交sv欧美| 国产人妻一区二区三区在| 大话2 男鬼变身卡| 久久国内精品自在自线图片| 大陆偷拍与自拍| 一区二区三区四区激情视频| 在线精品无人区一区二区三 | 22中文网久久字幕| 老司机影院成人| 国内精品美女久久久久久| 亚洲va在线va天堂va国产| 国产白丝娇喘喷水9色精品| 国产高清有码在线观看视频| 人妻一区二区av| 国产精品三级大全| 夜夜爽夜夜爽视频| 国产成年人精品一区二区| 国产综合懂色| 亚州av有码| 又黄又爽又刺激的免费视频.| 成人国产av品久久久| 日韩三级伦理在线观看| 老女人水多毛片| 国产精品国产三级国产av玫瑰| 伊人久久国产一区二区| 日韩制服骚丝袜av| 国产在线一区二区三区精| 搡女人真爽免费视频火全软件| 亚洲精品国产av蜜桃| 又大又黄又爽视频免费| 国产爽快片一区二区三区| 久久精品久久久久久噜噜老黄| 香蕉精品网在线| 黄片wwwwww| 国产一级毛片在线| 哪个播放器可以免费观看大片| 成人二区视频| 亚洲av成人精品一二三区| 联通29元200g的流量卡| 91在线精品国自产拍蜜月| 成年女人看的毛片在线观看| 麻豆乱淫一区二区| 国产精品秋霞免费鲁丝片| 美女被艹到高潮喷水动态| 在线观看美女被高潮喷水网站| 最近中文字幕2019免费版| 水蜜桃什么品种好| 亚洲成人中文字幕在线播放| tube8黄色片| 欧美性感艳星| 你懂的网址亚洲精品在线观看| 丝袜脚勾引网站| 黄片wwwwww| 国产日韩欧美在线精品| 校园人妻丝袜中文字幕| 日本黄大片高清| 日本黄色片子视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产真实伦视频高清在线观看| 中文在线观看免费www的网站| 日韩强制内射视频| 亚洲经典国产精华液单| 亚洲成人久久爱视频| 最新中文字幕久久久久| 人妻 亚洲 视频| 男男h啪啪无遮挡| 精品人妻偷拍中文字幕| 日本三级黄在线观看| 欧美另类一区| 欧美最新免费一区二区三区| 最近2019中文字幕mv第一页| 国产伦精品一区二区三区视频9| 99久久精品一区二区三区| 一二三四中文在线观看免费高清| 久久久久久九九精品二区国产| 狠狠精品人妻久久久久久综合| 亚洲内射少妇av| 97超碰精品成人国产| 综合色丁香网| 网址你懂的国产日韩在线| 能在线免费看毛片的网站| 国内精品美女久久久久久| 在线免费观看不下载黄p国产| 欧美日韩视频精品一区| 亚洲精品一区蜜桃| 一级片'在线观看视频| av在线亚洲专区| 久久97久久精品| 欧美区成人在线视频| 老司机影院成人| 免费高清在线观看视频在线观看| 欧美日韩亚洲高清精品| 九九久久精品国产亚洲av麻豆| 午夜精品一区二区三区免费看| 激情 狠狠 欧美| 人妻制服诱惑在线中文字幕| 99九九线精品视频在线观看视频| 久久久久国产网址| 亚洲va在线va天堂va国产| 国产永久视频网站| 亚洲成人一二三区av| 国内揄拍国产精品人妻在线| 国产精品麻豆人妻色哟哟久久| 久久久久久国产a免费观看| 99re6热这里在线精品视频| 91精品伊人久久大香线蕉| 天堂中文最新版在线下载 | 小蜜桃在线观看免费完整版高清| 偷拍熟女少妇极品色| 看十八女毛片水多多多| 美女cb高潮喷水在线观看| 2022亚洲国产成人精品| 国产精品久久久久久久久免| 日韩制服骚丝袜av| 亚洲在线观看片| 久久影院123| 欧美激情久久久久久爽电影| 最后的刺客免费高清国语| 91久久精品电影网| 亚洲精品第二区| 欧美日韩精品成人综合77777| 国产在线一区二区三区精| xxx大片免费视频| 国产精品人妻久久久久久| 亚洲,一卡二卡三卡| 久久久久久久久久人人人人人人| 免费观看的影片在线观看| 日韩av在线免费看完整版不卡| 午夜精品国产一区二区电影 | 久久韩国三级中文字幕| 亚洲第一区二区三区不卡| 最近2019中文字幕mv第一页| 97人妻精品一区二区三区麻豆| 亚洲天堂国产精品一区在线| 国产精品三级大全| 国产探花极品一区二区| 日韩欧美精品免费久久| 欧美精品人与动牲交sv欧美| 婷婷色av中文字幕| 精品久久久久久电影网| 一级毛片aaaaaa免费看小| 午夜免费鲁丝| 美女脱内裤让男人舔精品视频| 水蜜桃什么品种好| 日本色播在线视频| tube8黄色片| 激情 狠狠 欧美| 伦理电影大哥的女人| 久久久色成人| 亚洲精品自拍成人| 国产v大片淫在线免费观看| 七月丁香在线播放| 亚洲人成网站在线观看播放| 日日摸夜夜添夜夜爱|