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

    海南島早白堊世初期A型花崗巖成因: U-Pb年代學、地球化學及Nd-Hf同位素制約

    2016-08-04 09:56:44魏俊浩李艷軍翟玉林
    大地構造與成礦學 2016年3期
    關鍵詞:鋯石A型同位素

    李 翔, 魏俊浩, 李艷軍, 翟玉林

    (中國地質大學(武漢) 資源學院, 湖北 武漢 430074)

    海南島早白堊世初期A型花崗巖成因: U-Pb年代學、地球化學及Nd-Hf同位素制約

    李翔, 魏俊浩*, 李艷軍, 翟玉林

    (中國地質大學(武漢) 資源學院, 湖北 武漢 430074)

    海南島至今未見早白堊世初期 A型花崗巖報道。本文首次在石碌鐵礦床中部鉆孔中發(fā)現(xiàn)了隱伏的 A型花崗巖,巖性為黑云母二長花崗巖。對該隱伏巖體開展了年代學、地球化學和Nd-Hf同位素分析, 探討巖石成因及成巖構造背景。LA-ICP-MS鋯石U-Pb定年結果為133±3 Ma(MSWD=0.33), 為早白堊世初期巖漿活動產物。巖體屬高鉀鈣堿性巖系, 具有高硅(SiO2=75.22%~76.63%)、富堿(K2O+Na2O=8.01%~8.56%)和鐵(FeOT=1.62%~2.01%)、貧鈣(CaO=0.55%~0.69%)和鎂(MgO=0.03%~0.15%)特征。鋁飽和指數(shù) A/CNK=0.99~1.04, 為準鋁–弱過鋁質花崗巖。巖石整體富集 K、Rb等大離子親石元素(LILE)和 U、Th等高場強元素(HFSE), 明顯虧損 Ba、Sr。LREE富集, HREE相對虧損, 且具強烈的 Eu負異常(δEu=0.02~0.08)。巖相學和地球化學顯示其為高分異鋁質A型花崗巖。Nd同位素組成εNd(t)= –6.3~ –7.1, 二階段模式年齡tDM2(Nd)=1.44~1.47 Ga, Hf同位素組成εHf(t)= –18.5~ –21.7, 二階段模式年齡tDM2(Hf)=1.93~2.09 Ga, 揭示該隱伏花崗巖體源區(qū)主要為古老地殼基底, 并混合了少量富集地幔物質。黑云母二長花崗巖具有后碰撞花崗巖的特征, 形成于伸展構造背景, 與早白堊世初期印度板塊北向漂移擠壓及太平洋板塊俯沖作用減弱造成的局部伸展有關。

    A型花崗巖; 早白堊世; 巖石成因; 伸展背景; 海南

    0 引 言

    海南島主體大地構造位置屬于武夷–云開–臺灣造山系海南地塊五指山島弧帶(潘桂棠等, 2009), 受太平洋構造域和特提斯構造域兩大地球動力學系統(tǒng)控制, 表現(xiàn)出復雜的構造格局。區(qū)內巖石以巖漿巖分布面積最廣(圖1a), 占全島面積的53.9%, 其中侵入巖約占 40%, 以花崗質侵入體為主(汪嘯風等, 1991)。A型花崗巖是海南島發(fā)育的重要巖石類型之一, 目前報道的成巖年代主要集中于二疊紀–三疊紀。謝才富等(2006a)首次報道了瓊中地區(qū)發(fā)育的中二疊世的鉀玄質巖體(SHRIMP鋯石U-Pb年齡272± 7 Ma), 認為其與石炭紀–早二疊世華南板塊向印支–南海板塊俯沖后碰撞活動有關。至于三疊紀 A型花崗巖如尖峰嶺超單元、瓊中深堀村、麻山田單元等, SHRIMP和LA-ICP-MS鋯石U-Pb年齡為240~250 Ma (謝才富等, 2005, 2006b; 張小文等, 2009)。這些巖體形成于后造山伸展構造環(huán)境, 由富集地幔分異或殼?;旌先廴谛纬?謝才富等, 2005, 2006b; 周佐民等, 2011)。然而早白堊世初期A型花崗巖至今未見報道。

    近年來華南地區(qū)早白堊世初期巖漿事件屢見報道, 巖石類型有A型花崗巖(王強等, 2005)、火山巖(羅平等, 2009; 張家菁等, 2009a; 邱駿挺等, 2011;蘇慧敏, 2013; 楊水源, 2013)、I型花崗巖(李真等, 2009; Li et al., 2012)、中基性脈巖(Wang et al., 2008)和玄武巖(Chen et al., 2008)等。海南島地區(qū)也發(fā)育該時期巖漿活動事件, 汪嘯風等(1991)和侯威(1996)曾報道海南島131.7 Ma和134.4 Ma的花崗巖漿事件,但局限于黑云母K-Ar法和全巖Rb-Sr等時線法, 該期巖漿活動仍缺乏高精度的年代學制約。另外, 葛小月(2003)曾報道了海南島發(fā)育早白堊世(全巖K-Ar年齡136.4±2.2 Ma)基性巖脈, 認為其形成于伸展構造背景下, 為EMⅡ型富集地幔來源。但關于區(qū)內同時期的A型花崗巖巖石成因和構造背景研究仍缺乏。

    石碌鐵礦床位于海南島西部昌江縣城南側, 礦區(qū)南、北、西面均被侵入巖環(huán)繞, 巖石類型以印支–燕山期花崗巖類為主(汪嘯風等, 1991; 侯威, 1996),并發(fā)育印支期閃長玢巖脈和燕山晚期花崗閃長斑巖、花崗斑巖、石英斑巖、閃長巖、煌斑巖、輝綠巖等巖脈(汪嘯風等, 1991; 侯威, 1996; 王智琳等, 2011; 李艷軍等, 2013; Wang et al., 2015)。近年來海南礦業(yè)股份有限公司探礦過程中, 鉆孔 ZK2302、ZK2303等6個鉆孔均揭露了礦區(qū)中部隱伏的黑云母二長花崗巖。本文對該巖體開展了LA-ICP-MS鋯石U-Pb定年, 精確限定其形成于早白堊世初期, 并進行主量、微量元素、Nd-Hf同位素地球化學研究, 探討其成因類型和構造背景, 以期為研究該時期海南島的構造巖漿事件提供制約和參考。

    1. 中新生代蓋層; 2. 新生代玄武巖; 3. 白堊紀礫巖層; 4. 早古生代火山碎屑巖; 5. 晚古生代火山-沉積巖; 6. 新元古界石碌群; 7. 古中元古界抱板群和片麻狀花崗巖; 8. 新太古代雜巖; 9. 石碌群第1層至第6層; 10. 震旦系石灰頂組; 11. 中-下石炭統(tǒng); 12. 中-下二疊統(tǒng); 13. 變基性巖; 14. 130~90 Ma花崗巖; 15. 270~180 Ma花崗巖; 16. 三疊紀花崗閃長巖; 17. 侏羅紀花崗閃長巖; 18. 隱伏花崗閃長巖; 19. 花崗斑巖脈; 20. 礦體; 21. 向斜; 22. 背斜; 23. 實測及推測斷裂; 24. 地質界線; 25. 鉆孔及編號。圖1 海南島地質略圖(a, 據(jù)許德如等, 2007修改)和石碌鐵礦床地質簡圖(b, 據(jù)許德如等, 2008修改)Fig.1 Tectonic map of Hainan island (a) and geological map of the Shilu iron deposit (b)

    1 地質概況及巖石學特征

    石碌鐵礦床產于近 EW 向區(qū)域性昌江–瓊海深大斷裂與 NE向戈枕韌–脆性斷裂的交匯部位(圖1b)。礦區(qū)中心出露的地層主要為新元古界石碌群(Pt3SL)和石灰頂組(Pt3s), 二者為斷層接觸關系。石碌群為礦區(qū)內主要的賦礦巖系, 是一套(低)綠片巖相(局部可達角閃巖相)變質為主的碎屑沉積巖和碳酸鹽巖建造。礦區(qū)整體被一軸向為NNW–近EW向延伸的石碌復式向斜所控制。斷裂構造發(fā)育, 可劃分為NW向-NNW向、EW向-NEE向及近SN向-NNE向三組。構造體系經歷了加里東中晚期NE向擠壓、海西中期–燕山早期近NE向擠壓、燕山早期NW向擠壓和燕山晚期 NNE向壓扭四次構造作用(李艷軍等, 2013)。礦區(qū)南部、北部為印支–燕山早期斑狀/似斑狀(角閃)黑云母二長花崗巖、花崗閃長巖, 普遍具片麻狀構造, K-Ar和Rb-Sr年齡為190~320 Ma(汪嘯風等, 1991; 侯威, 1996), 但 LA-ICP-MS鋯石U-Pb年齡為262±1 Ma(李艷軍未發(fā)表數(shù)據(jù))。礦區(qū)西部為燕山晚期角閃黑云母二長花崗巖、花崗斑巖,其黑云母 K-Ar同位素和 Rb-Sr同位素年齡分別為131.7 Ma和134.4 Ma(汪嘯風等, 1991; 侯威, 1996)。區(qū)內發(fā)育花崗斑巖、石英斑巖、閃長巖、煌斑巖、輝綠巖等巖脈, K-Ar同位素年齡為97~100 Ma(侯威, 1996)。王智琳等(2011)應用 LA-ICP-MS鋯石 U-Pb法厘定礦區(qū)中部花崗閃長斑巖脈年齡為 93±2 Ma,李艷軍等(2013)獲得了礦區(qū)中部閃長玢巖脈鋯石U-Pb年齡為 248±1 Ma, 指示區(qū)內經歷了多期次的巖漿活動。

    隱伏黑云母二長花崗巖位于石碌鐵礦床中部雞心坳礦段。巖石呈似斑狀結構, 塊狀構造(圖2)。斑晶主要為石英、斜長石、鉀長石、黑云母。石英以它形粒狀結構為主, 粒徑 0.2~2.5 mm, 可見波狀消光, 含量約 15%; 斜長石呈灰白色–白色, 半自形,粒徑 0.5~3 mm, 發(fā)育卡氏雙晶, 偶見環(huán)帶構造, 含量約 15%, 部分斜長石內部見絹云母; 鉀長石呈半自形, 發(fā)育條紋結構, 為條紋長石, 粒徑0.5~3.5 mm, 含量約 10%; 黑云母呈褐色–褐綠色, 多為半自形–自形片狀, 粒徑0.3~1 mm, 含量約5%?;|為細?;◢徑Y構, 主要礦物為石英(25%±)、條紋長石(15%±)、斜長石(15%±), 副礦物組合為磷灰石–鋯石–磁鐵礦。

    (a) 基質細?;◢徑Y構; (b) 似斑狀結構。礦物名稱縮寫: Q. 石英; Pl. 斜長石; Bi. 黑云母; Kfs. 鉀長石。圖2 石碌二長花崗巖體斑晶及基質正交偏光顯微照片F(xiàn)ig.2 Micrographs of the phenocrysts and groundmass from the Shilu monzogranite

    2 測試方法

    用于鋯石U-Pb同位素定年的樣品為ZK2303-1,取于鉆孔ZK2303161 m處。鋯石的挑選在河北省廊坊市區(qū)域地質礦產調查研究所實驗室進行。將利用重磁技術分選的鋯石置于環(huán)氧樹脂中制靶, 磨蝕至鋯石核部出露, 在中國地質大學(武漢)地質過程與礦產資源國家重點實驗室(GPMR)進行陰極發(fā)光(CL)顯微照像, 結合透射光和反射光圖像, 觀察鋯石內部結構。鋯石U-Pb原位微區(qū)測試在GPMR激光剝蝕等離子體質譜儀(LA-ICP-MS)上完成, 激光束斑直徑為24 μm。鋯石U-Th-Pb同位素組成分析以標準鋯石 91500作為外標進行同位素分餾校正; 微量元素含量則利用NIST610作為外標,29Si作為內標的方法進行定量計算, 詳細分析流程和儀器參數(shù)見Liu et al. (2010)。實驗所得數(shù)據(jù)采用ICPMS-DataCal (ver3.0)程序進行處理, 普通 Pb校正見 Andersen (2002), U-Pb諧和年齡圖繪制和加權平均年齡計算用Isoplot3.23(Ludwig, 2003), 分析點表觀年齡采用206Pb/238U年齡。對已做過U-Pb定年的鋯石顆粒, 在GPMR多接收等離子體質譜儀(MC-ICP-MS)上進行Hf同位素原位微區(qū)測試, 激光束斑直徑為 32 μm,具體分析方法及儀器參數(shù)詳見Hu et al. (2012)。用179Hf/177Hf=0.7325及173Yb/171Yb=1.1248 (Blichert-Toft et al., 1997)標定測定樣品的179Hf/177Hf和173Yb/171Yb比值。實驗所得數(shù)據(jù)采用ICPMSDataCal程序(Liu et al., 2010)進行處理。

    室內對巖石樣品進行巖相學鑒定后, 挑選出 7件新鮮、無蝕變或蝕變較弱的樣品, 切取200 g樣品粉碎至 200目以下, 進行主量、微量元素、Sm-Nd同位素測試。主量、微量元素測試在澳實礦物實驗室集團澳實分析檢測(廣州)有限公司完成主量元素用XRF熒光光譜儀分析, 分析精度優(yōu)于5%。微量元素采用 Perkin Elmer Elan 9000等離子質譜儀(ICP-MS)分析, 分析精度優(yōu)于 10%。Sm-Nd同位素比值測試在武漢地質礦產研究所MAT261多接受質譜儀上完成, 全流程本底Nd為2.13×10–11。采用實驗室標準ZK-bzNd(Nd)控制儀器工作狀態(tài),147Sm/144Nd精度優(yōu)于0.5%。

    3 分析結果

    3.1 鋯石U-Pb年代學

    黑云母二長花崗巖(ZK2303-1)鋯石為無色–淺黃色透明狀, 以短軸狀和等軸狀為主, 粒徑多在50~120 μm之間。CL圖像顯示所測鋯石發(fā)育明顯的振蕩環(huán)帶, 基本無增生邊(圖3a)。選擇韻律環(huán)帶明顯的13顆鋯石進行了定年分析, 結果見表1。鋯石的U、Th含量分別為119~721 μg/g, 104~391 μg/g, Th/U比值為0.39~1.06, 屬巖漿成因(Th/U>0.1, Griffin et al., 2004)。輕、重稀土元素分餾明顯, 球粒隕石標準化稀土配分模式圖左傾, 表現(xiàn)為HREE富集、LREE虧損, 具較強烈的 Ce正異常及弱 Eu負異常(圖4),表明所分析的鋯石為典型的巖漿成因鋯石(Hoskin and Schaltegger, 2003)。僅測點ZK2303-1-01 LREE相對富集, 但 HREE與其他測點一致, 估計與后期地質事件擾動時 LREE優(yōu)先進入鋯石晶格有關(Wu and Zheng, 2004), 對應的CL圖像也顯示有后期的干擾。其余大部分數(shù)據(jù)點都位于諧和線上及附近,206Pb/238U年齡介于128±5 Ma~137±6 Ma之間, 加權平均年齡為133±3 Ma (MSWD=0.33) (圖3b), 代表黑云母二長花崗巖體的結晶年齡為早白堊世。

    圖(a)中實線圓圈表示年齡分析點位和編號, 虛線圓圈代表Hf同位素分析點位和編號, 標示年齡為206Pb/238U表面年齡。圖3 石碌二長花崗巖體鋯石陰極發(fā)光圖(a)和U-Pb年齡諧和圖(b)Fig.3 Cathodoluminescence (CL) images of the representative zircon grains (a), and U-Pb concordia diagram (b) for the zircon grains from the Shilu monzogranite

    3.2 巖石地球化學特征

    7件樣品主量和微量元素測試結果列于表2。巖石SiO2含量為75.22%~76.63%, 全堿(K2O+Na2O)含量8.01%~8.56%, Na2O為2.92%~3.5%, K2O/Na2O為1.36~1.74, 在SiO2-K2O圖解上(圖5a), 所有樣品投點位于高鉀鈣堿性巖系列。巖石普遍貧CaO(0.55%~ 0.69%)、P2O5(0.005%~0.016%)、TiO2(0.03%~0.09%) 和MgO(0.03%~0.15%), 富FeOT(1.62%~2.01%)。鋁飽和指數(shù)A/CNK=0.99~1.04, A/NK>1, 為準鋁–弱過鋁質特征(圖5b)。CIPW標準礦物計算顯示部分樣品出現(xiàn)剛玉分子, 堿性長石含量明顯高于鈣長石。巖體總體具高硅、富堿和鐵、貧鈣和鎂的特征。

    樣品ΣREE=137.14~245.11 μg/g (均值為182.2 μg/g), HREE=38.71~56.23 μg/g, LREE=86.94~203.42 μg/g, LREE/HREE為1.55~5.16(表2), 具有強烈的 Eu負異常(δEu=0.02~0.08)。在稀土元素球粒隕石標準化配分圖上(圖 6a)顯示輕稀土富集重稀土相對虧損的特點。(La/Yb)N=0.91~4.57(均值為 2.24), (Gd/Yb)N= 0.79~1.24, 重稀土分布曲線相對平坦。

    表1 石碌二長花崗巖鋯石U-Pb分析結果Table1 Zircon U-Pb isotope results of the Shilu monzogranite

    圖 4 二長花崗巖(ZK2303-2)鋯石球粒隕石標準化稀土配分模式圖解(球粒隕石值標準化值據(jù) Sun and McDonough, 1989)Fig.4 Chondrite normalized REE patterns for the zircon granis from the monzogranite (ZK2303-2)

    微量元素方面, 在原始地幔標準化蛛網圖上所有樣品具顯著的Rb、Th、U、K等元素的正異常和Ba、Sr、P、Ti等元素負異常(圖6b)。樣品的Ga含量較高, 均值為 20.7 μg/g, 104×Ga/Al比值介于3.02~3.38, Zr、Nb、Ce、Y 等元素含量均較高, Zr+Nb+Ce+Y=223.1~ 334.8 μg/g。

    3.3 Nd-Hf同位素

    對4件樣品進行了全巖Sm-Nd同位素分析, 測試結果列于表 3。Nd同位素組成比較均一, 樣品143Nd/144Nd=0.512103~0.512145, εNd(t)變化于–6.3~ –7.1, 表明它們的物質來源具有相似性。對應的二階段Nd模式年齡tDM2(Nd)=1.44~1.47 Ga。

    在鋯石U-Pb定年基礎上, 利用LA-MC-ICP-MS對其中11顆鋯石進行了11個點的Lu-Hf同位素分析, 結果列于表 4。176Lu/177Hf比值為 0.000744~ 0.001332, 顯示鋯石在形成之后具有較低的放射性同位素Hf積累。所有11個點的176Hf/177Hf比值為0.282077~0.282165。以各單顆粒鋯石結晶年齡計算出鋯石εHf(t)= –18.5~ –21.7, 對應的二階段Hf模式年齡tDM2(Hf)=1.93~2.09 Ga。

    4 討 論

    4.1 A型花崗巖的厘定

    圖5 石碌二長花崗巖體SiO2-K2O圖解(a, 底圖引自Collins et al. , 1982)和A/CNK-A/NK圖解(b, 底圖據(jù)Maniar and Piccoli, 1989)Fig.5 SiO2vs. K2O (a) and A/CNK vs. A/NK (b) diagrams for the Shilu monzogranite

    表2 石碌二長花崗巖主量(%)、微量和稀土元素(μg/g)分析結果Table 2 Major (%), trace and rare earth (μg/g) element compositions of the Shilu monzogranite

    續(xù)表2:

    圖 6 石碌二長花崗巖體球粒隕石標準化稀土元素配分模式圖(a)和原始地幔標準化微量元素蛛網圖(b)(標準化值引自Sun and McDonough, 1989)Fig.6 Chondrite normalized REE patterns (a), and primitive mantle-normailized trace element spider diagram (b), for the Shilu monzogranite

    表3 石碌二長花崗巖全巖Sm-Nd同位素分析結果Table 3 Sm-Nd isotopic compositions of the Shilu monzogranite

    表4 石碌二長花崗巖LA-MC-ICPMS鋯石Hf同位素分析結果Table 4 LA-MC-ICPMS ziron Lu-Hf isotopic compositions of the Shilu monzogranite

    目前最常用的花崗巖分類方法是ISMA分類法,其中I、S、M型花崗巖主要是依據(jù)源巖性質劃分。A型花崗巖最早由Loiselle and Wones (1979)提出來,他們將 A 型花崗巖定義為堿性(alkaline)、貧水(anhydrous)、非造山(anorogenic)的花崗巖, 不涉及物質來源。賈小輝等(2009)綜合國內外學者的研究, 對A型花崗巖特征分別從巖石類型、礦物學特征、地球化學成分上進行了總結, 指出: (1)巖石類型上, A型花崗巖不僅包括堿性巖類, 還包括堿鈣性、弱堿–準鋁、弱過鋁甚至強過鋁質巖石, 幾乎囊括除I、S型花崗巖以外的其他花崗巖; (2)礦物學特征上,主要的礦物組合為石英+(富 Fe)鎂鐵質暗色礦物±堿性長石; (3)化學成分方面, A型花崗巖富硅、堿,貧鈣、鎂、鋁, 高(K2O+Na2O)/Al2O3和 FeOT/MgO比值, 富 Rb、Th、Ta、Y, 貧 Sr、Ba, 具有顯著的Eu負異常, Ga/Al值高, 稀土元素配分模式常呈右傾的海鷗型。本文研究的石碌黑云母二長花崗巖屬于高鉀鈣堿性弱過鋁質巖石, 礦物組合上以石英、斜長石、鉀長石、黑云母為主(圖2), 化學成分高Si、Na和K, 低Ca、Mg, 富Rb、Th、Ta、Y等高場強元素, 貧 Sr、Ba, 稀土元素配分曲線呈現(xiàn)輕稀土富集的右傾型配分模式, 具顯著的Eu負異常, 與A型花崗巖具有一致的地球化學特征。樣品的104×Ga/Al 值(3.02~3.38)和 Zr+Nb+Ce+Y值(223.1~334.8 μg/g)較高, 在Whalen et al. (1987)的成因類型判別圖解中(圖7), 均落入A型花崗巖區(qū)域。

    盡管A型花崗巖與高分異的I型、S型花崗巖在地球化學組成上有相似的特點往往很難區(qū)分, 但S型花崗巖通常是過鋁質并且隨著分異程度的增加而更加強烈, 僅根據(jù)特征礦物就可以區(qū)分 S型和 A型花崗巖(King et al., 1997)。此外, S型花崗巖具有高P2O5(均值為0.14%), 低Na2O(均值為2.81%)的特征, 且 P2O5含量與巖石分異程度呈正相關關系(Chappell and White, 1992)。石碌黑云母二長花崗巖中并沒有發(fā)現(xiàn)白云母、石榴石、堇青石等典型富鋁礦物, A/CNK 變化范圍為0.99~1.04, 主體屬準鋁–弱過鋁質。樣品P2O5含量較低(0.005%~ 0.016%), Na2O含量高(2.92%~3.5%), P2O5含量并沒有與 SiO2含量呈現(xiàn)正相關關系, 明顯不同于 S型花崗巖的礦物組合和化學特征, 排除其為高分異的S型花崗巖。

    至于高分異的I型花崗巖與A型花崗巖的區(qū)分,目前普遍采用FeOT含量、形成溫度及判別圖解等判別標準加以判定(Whalen et al., 1987; 王強等, 2000;吳福元等, 2007; 賈小輝等, 2009)。A型花崗巖全鐵(FeOT)含量較高(>1.00%), 且形成溫度高(>800 ℃),而高分異I型花崗巖FeOT含量一般小于1.00%, 形成溫度較低(均值764 ℃)(王強等, 2000; 賈小輝等, 2009; 李小偉等, 2010)。石碌黑云母二長花崗巖體FeOT含量為1.62%~2.01%, 據(jù)Watson et al. (2006)Ti溫度計計算的鋯石結晶溫度均值為 814 ℃(表 1),符合上述A型花崗巖的特征。樣品高硅、高分異指數(shù)(DI最小為93.7%)及強烈的Eu負異常均反映巖體經歷了高程度的分異演化作用。綜上, 石碌隱伏二長花崗巖體應屬于高分異鋁質A型花崗巖。

    圖7 石碌二長花崗巖體成因類型判別圖解(底圖據(jù)Whalen et al., 1987)Fig.7 Petrogenesis discrimination diagrams for the Shilu monzogranite

    4.2 巖石成因

    A型花崗巖成因一直飽受爭議, 提出的成因模式也各式各樣, 主要的爭議焦點是巖漿源區(qū)來源問題。目前關于A型花崗巖的巖漿源區(qū)的觀點歸納起來主要有: (1)幔源巖漿結晶分異(Eby, 1990, 1992; Turner et al., 1992; Frost et al., 1999); (2)幔源巖漿與殼源巖漿混合(邱檢生等, 1996; Yang et al., 2006;李艷軍等, 2014); (3)地殼物質部分熔融(Collins et al., 1982; Whalen et al., 1987)。本文所研究的石碌黑云母二長花崗巖樣品明顯高硅(SiO2含量 75.22%~ 76.63%)、富堿(K2O+Na2O=8.01%~8.56%), 相對富集大離子親石元素和輕稀土元素, 顯示陸殼物質參與了成巖作用。較低的 Nb/Ta比值(6.16~9.53), 高的tDM2(Hf) (1.93~2.09 Ga)和tDM2(Nd)(1.43~1.50 Ga)值,同樣顯示地殼物質對成巖作用的貢獻。

    巖石Hf同位素特征(εHf(t)= –18.5~ –21.7)表明古老的地殼基底對成巖貢獻很大。海南島結晶基底為抱板群, Li et al. (2002)測得其鋯石U-Pb同位素諧和年齡為1.43 Ga, 確定其為中元古代巖系。本次測得的石碌隱伏花崗巖體二階段Nd模式年齡tDM2(Nd)= 1.43~1.50 Ga, 與基底抱板群的形成年齡一致。εNd(t)值介于–6.28和–7.1之間(表3), 在t–εNd(t)圖上(圖8),樣品位于抱板群地殼演化線和原始地幔演化線之間,顯示有幔源物質參與。較高的鋯石結晶溫度(均值為814 ℃, 最高達 952 ℃)也為地幔物質的上侵提供了依據(jù)(李小偉等, 2010)。因此, 我們傾向于認為石碌黑云母二長花崗巖是幔源物質上侵致使原始基底物質重熔形成的, 而在此過程中伴隨有少量幔源物質的混入。葛小月(2003)報道了海南島早白堊世(136.4±2.2 Ma)基性巖脈的存在, 并指出其為 EMⅡ型富集地幔來源, 為海南島這一時期富集地幔物質的上侵提供了佐證。同時期華南新路火山盆地中的A型花崗巖(133~136 Ma, 楊水源, 2013)及武夷地區(qū)天華山盆地早白堊世的火山–侵入巖(137~144 Ma,蘇慧敏, 2013)研究表明其成巖過程也有地幔物質的參與。另外, 武夷地區(qū)同時代的毛斷高演化I型花崗巖(140.0±1.6 Ma, Li et al., 2012)、南嶺地區(qū)基性巖脈(146.2±2.3 Ma, Wang et al., 2008)和何嶺玄武巖(136.8±0.8 Ma, Chen et al., 2008)等研究表明這一時期華南地區(qū)上侵的幔源物質主要為富集地幔。綜上,石碌黑云母二長花崗巖體極有可能是富集地幔上侵, 致使古老基底重熔, 在此過程中伴隨少量富集地幔物質混合, 熔體經過高度分異演化過程而侵位形成的。

    圖8 石碌二長花崗巖體t-εNd(t)圖解(抱板群Nd同位素數(shù)據(jù)引自雷裕紅等, 2005)Fig.8 t vs. εNd(t) diagram for the Shilu monzogranite

    4.3 成巖構造背景

    海南島地處華南地塊與印支地塊之間, 受太平洋構造域和特提斯構造域兩大地球動力學系統(tǒng)的影響, 具有復雜的地質構造演化歷史。然而, 對海南島基底屬性及其構造區(qū)劃目前尚有爭議, 主要有四種觀點: (1)海南島陸殼是華夏古陸的一部分(譚忠福等, 1991; 丁式江, 1995; 張業(yè)明等, 1997); (2)是中生代從親東岡瓦納(澳大利亞)裂解的塊體(殷鴻福等, 1999); (3)以白沙斷裂或九所–陵水斷裂或昌江–瓊海斷裂為界把海南島劃分為不同的塊體(Metcalfe et al., 1993; 雷裕紅等, 2005; Zhang et al., 2011; Xu et al., 2013); (4)以邦溪–晨星蛇綠巖帶作為華南塊體與印支塊體的縫合帶(Li et al., 2002; Cai and Zhang, 2009;陳新躍等, 2011)。第四種觀點目前基本被認可。邦溪–晨星蛇綠巖Sm-Nd等時線年齡為333±12 Ma (Li et al., 2002), 鋯石U-Pb年齡為345±4 Ma陳新躍等, 2013), 表明古特提斯洋東段開始形成于石炭紀。該洋盆在早三疊世消減關閉使華南(包括海南島北部)和印支塊體(包括海南島中南部)碰撞拼合(Li et al., 2002; 陳新躍等, 2013; 溫淑女等, 2013)。自此海南島與華南地塊成為統(tǒng)一的地塊。

    華南地區(qū)中生代巖漿活動發(fā)育, 前人對成巖構造背景研究主要集中于 230~210 Ma、170~150 Ma 和100~90 Ma三個峰值階段, 早白堊世初期巖漿構造背景很少專門論述。然而, 近年來華南的武夷山及南嶺構造帶中有越來越多的早白堊世初期成巖成礦事件的報道, 如蘇村(王強等, 2005)、毛斷(Li et al., 2012)、金山(李真等, 2009)等花崗巖體, 鵝湖嶺組等火山巖(蘇慧敏, 2013), 南嶺地區(qū)基性巖脈(146.2± 2.3 Ma, Wang et al., 2008)和何嶺玄武巖(136.8±0.8 Ma, Chen et al., 2008)等, 以及毛斷(輝鉬礦Re-Os模式年齡 139.0±0.8 Ma, Li et al., 2012)、銅坑嶂(輝鉬礦Re-Os等時線年齡133.8±0.65 Ma, 許建祥等, 2007)、金竹坪(輝鉬礦Re-Os等時線年齡135.5±5.7 Ma, 張家菁等, 2009b)等 Mo-Pb-Zn礦床, 鋸板坑(絹云母Ar-Ar坪年齡139.2±1.5 Ma, 付建明等, 2009)、茅坪(輝鉬礦 Re-Os模式年齡 141.4±2.2 Ma, 曾載淋等, 2009)和大吉山(絹云母 Ar-Ar坪年齡 144.0±0.7 Ma,張文蘭等, 2006)等 W-Sn礦床, 馬坑鐵礦床(輝鉬礦Re-Os模式年齡133.0±0.8 Ma, 張承帥等, 2012)等。另外, 東南沿海火山–侵入巖帶也發(fā)育 145~124 Ma的峰值階段(Guo et al., 2012)。這些地質事件的報道表明早白堊世初期也應為華南地區(qū)一個重要的巖漿活動階段, 甚至是一個重要的成巖成礦時段(Li et al., 2012)。前人也指出該時期是個重要的地球動力學調整的時間段(毛景文等, 2004; 吳根耀, 2006; 薛懷民等, 2009; 毛建仁等, 2014)。對于早白堊世初期的構造背景, 前人普遍認為是伸展背景(王強等,2005; 張岳橋等, 2012)。該時期的 A型花崗巖、A型火山巖、I型花崗巖、基性脈巖及玄武巖的發(fā)現(xiàn),也表明該時期華南地區(qū)處于伸展背景無疑。王強等(2005)總結出139~123 Ma(峰期135±5 Ma)是華南一個重要的 A型花崗巖侵入活動時期。同時, 構造地質學研究也表明, 此時華南處于伸展構造背景(張岳橋等, 2012)。但是, 該時期伸展活動的動力學機制仍存在爭議, 目前主要有大陸伸展和裂谷(Gilder et al., 1996)、古太平洋板塊高角度俯沖(Jahn, 1974; Zhou and Li, 2000; Chen et al., 2014)、古太平洋板塊平板俯沖(Li and Li, 2007; Li et al., 2012)或俯沖角度由緩變陡地殼伸展垮塌引起的區(qū)域性伸展(張岳橋等, 2012)等成因模式。

    然而, 詳細的構造應力場研究表明, 從135 Ma開始, 華南地區(qū)構造主壓應力場是以 NNE-SSW 向近水平擠壓為主要特征, 動力源來自于西南方向(萬天豐, 2004)。Wang et al. (2013)也指出142~132 Ma華夏地塊普遍發(fā)育NW-SE向擠壓變形, 受到SW向動力的擠壓作用。另外, 李艷軍等(2013)通過應力場分析也已指出海南石碌鐵礦區(qū)晚燕山期 NNE向壓扭性斷裂由NNE-SSW向擠壓作用形成。顯然, 上述構造成因模式明顯不能解釋早白堊世初期的應力場。一種可能的解釋是早白堊世古太平洋板塊的俯沖作用可能受到了印度板塊北向漂移(Jaeger et al., 1989; Beck et al., 1995; Gilder et al., 1996, 1999; 譚俊等, 2008; 李艷軍等, 2009, Wang et al., 2013)的影響, 華南地區(qū)處于NE向局部伸展環(huán)境。這一模式明顯符合上述NW向擠壓應力場。受印度板塊快速北移的影響, 古太平洋板塊的俯沖作用可能被逐漸減弱(王強等, 2005), 甚至發(fā)生間歇性后撤, 并出現(xiàn)反轉(Tatsumi et al., 1990; Ren et al., 2002; Zhao et al., 2013), 其中135 Ma是個重要的轉換時間點(薛懷民等, 2009; 周濤發(fā)等, 2012)。在此構造背景下, 形成NE向局部伸展環(huán)境。此時期華南地區(qū)形成了一系列NE向早白堊世斷陷盆地, 巖漿侵入作用和火山作用強烈(張岳橋等, 2012)。正是這一巖漿作用, 形成了蘇村 A型花崗巖(王強等, 2005)、毛斷(Li et al., 2012)和金山(李真, 2009)等I型花崗巖、鵝湖嶺組火山巖(蘇慧敏, 2013)、基性巖脈(Wang et al., 2008)和何嶺玄武巖(Chen et al., 2008)等表征伸展構造背景的地質體。

    石碌隱伏黑云母二長花崗巖體鋯石 U-Pb年齡為 133±3 Ma (MSWD=0.33), 與上述構造巖漿事件時間一致。在104×Ga/Al-R1和Nb-Y-3×Ga圖解中(圖9), 樣品投點于 PA、A2型花崗巖區(qū)域, 表明巖石形成于板塊碰撞后或后造山伸展構造環(huán)境(Eby, 1992; Hong et al., 1996), 微量元素(Y+Nb)-Rb判別圖解(圖10)中樣品均位于火山弧花崗巖、同碰撞花崗巖與板內花崗巖交界部位, 這一范圍也是后碰撞花崗巖的投影區(qū)域(Pearce, 1996; F?rster et al., 1997), 而非板內裂谷等構造環(huán)境的產物。同時, 瓊南地區(qū)136.4±2.2 Ma富集地幔來源基性巖脈(葛小月, 2003)也反映了該時期海南島處于伸展構造背景。因此,我們認為石碌隱伏黑云母二長花崗巖的侵位與早白堊世初期印度板塊北向漂移擠壓及太平洋板塊俯沖作用減弱造成的局部伸展作用有關。

    圖a中: R1=4×Si-11×(Na+K)-2×(Fe+Ti), PA. 造山后A型花崗巖, AA. 非造山A型花崗巖; 圖b中: A1. 大陸裂谷或板內環(huán)境, A2. 陸-陸碰撞或島弧環(huán)境。圖9 石碌二長花崗巖體104×Ga/Al-R1 (a, 據(jù)Hong et al. , 1996)和Nb-Y-3×Ga (b, 據(jù)Eby, 1992)圖解Fig.9 104×Ga/Al vs. R1(a) and Nb-Y-3×Ga (b) diagrams for the Shilu monzogranite

    Syn-COLG. 同碰撞花崗巖; WPG. 板內花崗巖; VAG. 火山弧花崗巖; ORG. 洋脊花崗巖; Post-COLG. 后碰撞花崗巖。圖 10 石碌二長花崗巖體(Y+Nb)-Rb判別圖解(底圖據(jù)Pearce, 1996)Fig.10 (Y+Nb) vs. Rb diagram for the Shilu monzogranite

    5 結 論

    (1) 石碌黑云母二長花崗巖 LA-ICP-MS鋯石U-Pb年齡為133±3 Ma, 形成于早白堊世初期。

    (2) 巖相學、地球化學及 Nd-Hf同位素特征表明石碌黑云母二長花崗巖為高分異鋁質A型花崗巖,是富集地幔上侵致使基底重熔, 并伴隨少量富集地幔物質混合形成的。

    (3) 石碌黑云母二長花崗巖體具有后碰撞花崗巖的特征, 形成于早白堊世初期印度板塊北向漂移擠壓及太平洋板塊俯沖作用減弱造成的局部伸展構造背景。

    致謝: 論文承蒙中山大學王岳軍教授及另外一位匿名審稿專家的審閱, 對他們提出的寶貴意見和建議表示誠摯的謝意。野外工作過程中得到了海南礦業(yè)股份有限公司陳福雄部長、王永輝和梁將科長的幫助, 在此一致表示感謝。

    陳新躍, 王岳軍, 范蔚茗, 張菲菲, 彭頭平, 張玉芝. 2011.海南五指山地區(qū)花崗片麻巖鋯石LA-ICP-MS U-Pb年代學特征及其地質意義. 地球化學, 40(5): 454–463.

    陳新躍, 王岳軍, 張玉芝, 張菲菲, 溫淑女. 2013. 海南晨星安山質火山巖地球化學、年代學特征及其構造意義.大地構造與成礦學, 37(1): 99–108.

    丁式江. 1995. 瓊西戈枕侵入體的地質特征及其發(fā)現(xiàn)意義.大地構造與成礦學, 19(4): 336–343.

    付建明, 李祥能, 程順波, 徐德明, 馬麗艷, 陳希清. 2009.粵北連平地區(qū)鎢錫多金屬礦床成礦時代研究. 中國地質, 36(6): 1331–1339.

    葛小月. 2003. 海南島中生代巖漿作用及其構造意義. 廣州: 中國科學院研究生院博士學位論文.

    侯威. 1996. 海南島大地構造與金成礦學. 北京: 科學出版社: 1–229.

    賈小輝, 王強, 唐功建. 2009. A型花崗巖的研究進展及意義. 大地構造與成礦學, 33(3): 465–480.

    雷裕紅, 丁式江, 馬昌前, 傅楊榮, 楊坤光, 龍文國. 2005.海南島地殼生長和基底性質的Nd同位素制約. 地質科學, 40(3), 439–456.

    李小偉, 莫宣學, 趙志丹, 朱弟成. 2010. 關于A型花崗巖判別過程中若干問題的討論. 地質通報, 29(2-3): 278–285.

    李艷軍, 魏俊浩, 陳華勇, 李歡, 陳沖, 侯本俊. 2014. 義敦島弧帶夏塞早白堊世A型花崗巖成因: 鋯石U-Pb年代學、地球化學及Hf同位素制約. 大地構造與成礦學, 38(4): 939–953.

    李艷軍, 魏俊浩, 伍剛, 譚俊, 石文杰, 趙少卿, 王永輝. 2013. 海南石碌地區(qū)早三疊世閃長玢巖脈U-Pb年代學及構造意義. 地球科學, 38(2): 241–252.

    李艷軍, 魏俊浩, 姚春亮, 鄢云飛, 譚俊, 付樂兵, 潘錦勃, 李偉. 2009. 浙東南石平川花崗巖體LA-ICP-MS鋯石U-Pb年代學及構造意義. 地質論評, 55(5): 673–684.

    李真, 邱檢生, 蔣少涌, 徐夕生, 胡建. 2009. 福建金山花崗質復式巖體的元素和同位素地球化學及其成因研究. 地質學報, 83(4): 515–527.

    羅平, 吳淦國, 張達, 狄永軍, 王長明, 王永慶. 2009. 北武夷生米坑鉛鋅礦床地質地球化學特征與成因探討.地質力學學報, 15(4): 349–362.

    毛建仁, 厲子龍, 葉海敏. 2014. 華南中生代構造-巖漿活動研究: 現(xiàn)狀與前景. 中國科學(D輯), 44(12): 2593–2617.

    毛景文, 謝桂青, 李曉峰, 張長青, 梅燕雄. 2004. 華南地區(qū)中生代大規(guī)模成礦作用與巖石圈多階段伸展. 地學前緣, 11(1): 45–55.

    潘桂棠, 肖慶輝, 陸松年, 鄧晉福, 馮益民, 張克信, 張智勇, 王方國, 邢光福, 郝國杰, 馮艷芳. 2009. 中國大地構造單元劃分. 中國地質, 36(1): 1–16.

    邱檢生, 王德滋, 彭亞鳴, 周金城. 1996. 浙江舟山桃花島堿性花崗巖的巖石學和地球化學特征及成因探討.南京大學學報(自然科學版), 32(1): 82–91.

    邱駿挺, 余心起, 吳淦國, 屈文俊, 狄永軍, 張達, 羅平,杜安道. 2011. 北武夷篁碧礦區(qū)逆沖推覆構造及其與鉬、鉛–鋅成礦作用關系. 地學前緣, 18(5): 243–255.

    蘇慧敏. 2013. 北武夷天華山盆地火山–侵入巖的成因及其與成礦關系的研究. 北京: 中國地質大學(北京)博士學位論文.

    譚俊, 魏俊浩, 李水如, 王忠銘, 付樂兵, 張可清. 2008.廣西昆侖關A型花崗巖地球化學特征及構造意義. 地球科學, 33(6): 743–754.

    譚忠福, 馬國干, 陳明是, 陳平, 黃照先, 蔣大海. 1991.海南抱板群同位素年齡及其地質意義. 中國區(qū)域地質, (3): 241–245.

    萬天豐. 2004. 中國大地構造學綱要. 北京: 地質出版社: 1–387.

    汪嘯風, 馬大銓, 蔣大海. 1991. 海南島地質(二)巖漿巖.北京: 地質出版社: 1–273.

    王強, 趙振華, 簡平, 熊小林, 包志偉, 橦戴謨, 許繼峰,馬金龍. 2005. 華南腹地白堊紀A型花崗巖類或堿性侵入巖年代學及其對華南晚中生代構造演化的制約.巖石學報, 21(3): 795–808.

    王強, 趙振華, 熊小林. 2000. 桐柏–大別造山帶燕山晚期A型花崗巖的厘定. 巖石礦物學雜志, 19(4): 297–306.

    王智琳, 許德如, 張玉泉, 陳福雄, 王力, 吳俊. 2011. 海南石碌鐵礦床花崗閃長斑巖的鋯石LA-ICP-MS U-Pb定年及地質意義. 大地構造與成礦學, 35(2): 292–299.

    溫淑女, 梁新權, 范蔚茗, 王岳軍, 池國祥, 梁細榮, 周云, 蔣英. 2013. 海南島樂東地區(qū)志仲巖體鋯石U-Pb年代學、Hf同位素研究及其構造意義. 大地構造與成礦學, 37(2): 294–307.

    吳福元, 李獻華, 楊進輝, 鄭永飛. 2007. 花崗巖成因研究的若干問題. 巖石學報, 23(6): 1217–1238.

    吳根耀. 2006. 白堊紀: 中國及鄰區(qū)板塊構造演化的一個重要變換期. 中國地質, 33(1): 64–77.

    謝才富, 朱金初, 丁式江, 張業(yè)明, 陳沐龍, 付楊榮, 付太安, 李志宏. 2006b. 海南尖峰嶺花崗巖體的形成時代、成因及其與抱倫金礦的關系. 巖石學報, 22(10): 2493–2508.

    謝才富, 朱金初, 丁式江, 張業(yè)明, 付太安, 李志宏. 2006a. 瓊中海西期鉀玄質侵入巖的厘定及其構造意義. 科學通報, 51(16): 1944–1954.

    謝才富, 朱金初, 趙子杰, 丁式江, 付太安, 李志宏, 張業(yè)明, 徐德明. 2005. 三亞石榴霓輝石正長巖的鋯石SHRIMP U-Pb年齡: 對海南島海西–印支期構造演化的制約. 高校地質學報, 11(1): 47–57.

    許德如, 馬馳, 李鵬春, 夏斌, 張玉泉. 2007. 海南島變碎屑沉積巖鋯石SHRIMP U-Pb年齡及地質意義. 地質學報, 81(3): 381–393.

    許德如, 王力, 肖勇, 劉朝露, 符啟基, 蔡周榮, 黃居銳. 2008. “石碌式”鐵氧化物–銅(金)–鈷礦床成礦模式初探. 礦床地質, 27(6): 681–694.

    許建祥, 曾載淋, 李雪琴, 劉俊生, 陳鄭輝, 劉善寶, 郭春麗, 王成輝. 2007. 江西尋烏銅坑嶂鉬礦床地質特征及其成礦時代. 地質學報, 81(7): 924–928.

    薛懷民, 汪應庚, 馬芳, 汪誠, 王德恩, 左延龍. 2009. 皖南太平–黃山復合巖體的SHRIMP年代學: 由鈣堿性向堿性轉變對揚子克拉通東南部中生代巖石圈減薄時間的約束. 中國科學(D輯), 39(7): 979–993.

    楊水源. 2013. 華南贛杭構造帶含鈾火山盆地巖漿巖的成因機制及動力學背景. 南京: 南京大學博士學位論文.殷鴻福, 吳順寶, 杜遠生, 彭元橋. 1999. 華南是特提斯多島洋體系的一部分. 地球科學, 24(1): 3–14.

    張承帥, 毛景文, 謝桂青, 趙財勝, 于淼, 王金祥, 劉武剛. 2012. 福建馬坑矽卡巖型鐵(鉬)礦床地質特征及輝鉬礦Re-Os同位素年齡. 吉林大學學報(地球科學版), 42(S1): 224–236.

    張家菁, 施光海, 童貴生, 張智宇, 劉海, 吳榮土, 陳磊. 2009a. 浙江徐家墩鵝湖嶺組含銅多金屬礦火山巖的地球化學與年代學. 地質學報, 83(6): 791–799.

    張家菁, 吳木森, 陳鄭輝, 劉善寶, 李立興, 邱良明, 吳斌, 黃安杰, 祝平俊. 2009b. 江西省上饒縣金竹坪鉬多金屬礦床成礦年代學研究. 巖礦測試, 28(3): 228–232.

    張文蘭, 華仁民, 王汝成, 陳培榮, 李惠民. 2006. 贛南大吉山花崗巖成巖與鎢礦成礦年齡的研究. 地質學報, 80(7): 956–962.

    張小文, 向華, 鐘增球, 周漢文, 張利, 楊念, 王婧. 2009.海南尖峰嶺巖體熱液鋯石U-Pb定年及微量元素研究:對熱液作用及抱倫金礦成礦時代的限定. 地球科學, 34(6): 921–930.

    張業(yè)明, 張仁杰, 姚華舟, 馬國干. 1997. 海南島前寒武紀地殼構造演化. 地球科學, 22(4): 59–64.

    張岳橋, 董樹文, 李建華, 崔建軍, 施煒, 蘇金寶, 李勇. 2012. 華南中生代大地構造研究新進展. 地球學報, 33(3): 257–279.

    周濤發(fā), 范裕, 袁峰, 鐘國雄. 2012. 長江中下游成礦帶地質與礦產研究進展. 巖石學報, 28(10): 3051–3066.

    周佐民, 謝才富, 徐倩, 高大飛. 2011. 海南島中三疊世正長巖–花崗巖套的地質地球化學特征與構造意義.地質論評, 57(4): 515–531.

    曾載淋, 張永忠, 朱祥培, 陳鄭輝, 王成輝, 屈文俊. 2009.贛南崇義地區(qū)茅坪鎢錫礦床錸–鋨同位素定年及其地質意義. 巖礦測試, 28(3): 209–214.

    Andersen T. 2002. Correction of common lead in U-Pbanalyses that do not report204Pb. Chemical Geology, 192: 59–79.

    Beck R A, Burbank D W, Sercombe W J, Riley G W, Barndt J K, Berry J R, Afzal J, Khan A M, Jurgen H and Metje J R. 1995. Stratigraphic evidence for an early collision between northwest India and Asia. Nature, 373: 55–58.

    Blichert-Toft J and Albarède F. 1997. The Lu-Hf isotope geochemistry of chondrites and the evolution of the mantle-crust system. Earth and Planetary Science Letters, 148: 243–258.

    Blichert-Toft J, Chauvel C and Albarède F. 1997. Separation of Hf and Lu for high-precision isotope analysis of rock samples by magnetic sector-multiple collector ICP-MS. Contributions to Mineralogy and Petrology, 127: 248–260.

    Cai J X and Zhang K J. 2009. A new model for the Indochina and South China collision during the Late Permian to the Middle Triassic. Tectonophysics, 467: 35–43.

    Chappell B W and White A J R.1992. I- and S-type granites in the Lachlan Fold Belt. Transactions of the Royal Society of Edinburgh: Earth Sciences, 83: 1–26.

    Chen C H, Lee C Y and Shinjo R. 2008. Was there Jurassic paleo-Pacific subduction in South China? Constraints from40Ar/39Ar dating, elemental and Sr-Nd-Pb isotopic geochemistry of the Mesozoic basalts. Lithos, 106: 83–92.

    Chen L, Zhao Z F and Zheng Y F. 2014. Origin of andesitic rocks: Geochemical constraints from Mesozoic volcanics in the Luzong basin, South China. Lithos, 190-191: 220–239.

    Collins W J, Beams S D, White A J R and Chappell B W. 1982. Nature and origin of A-type granites with particular reference to southeastern Australia. Contributions to Mineralogy and Petrology, 80: 189–200.

    DePaolo D J and Wasserburg G J. 1979. Petrogenetic mixing models and Nd-Sr isotopic patterns. Geochimica et Cosmochimica Acta, 43: 615–627.

    Eby G N. 1990. The A-type granitoids: A review of their occurrence and chemical characteristics and speculations on their petrogenesis. Lithos, 26: 115–134.

    Eby G N. 1992. Chemical subdivision of the A-type granitoids: Petrogenetic and tectonic implications. Geology, 20: 641.

    F?rster H J, Tischendorf G and Trumbull R B. 1997. An evaluation of the Rb vs. (Y+Nb) discrimination diagram to infer tectonic setting of silicic igneous rocks. Lithos, 40: 261–293.

    Frost C D, Frost B R, Chamberlain K R and Edwards B R. 1999. Petrogenesis of the 1.43 Ga Sherman batholith, SE Wyoming, USA: A reduced, rapakivi-type anorogenic granite. Journal of Petrology, 40: 1771–1802.

    Gilder S A, Gill J, Coe R S, Zhao X X, Liu Z W, Wang G X, Yuan K, Liu W L, Kuang G and Wu H R. 1996. Isotopic and paleomagnetic constraints on the Mesozoic tectonic evolution of south China. Journal of Geophysical Research: Solid Earth, 101: 16137, 16154.

    Gilder S A, Leloup P H E, Courtillot V, Chen Y, Coe R S, Zhao X X, Xiao W J, Halim N, Cogn Eacute J P and Zhu R X. 1999. Tectonic evolution of the Tancheng-Lujiang (Tan-Lu) fault via Middle Triassic to early Cenozoic paleomagnetic data. Journal of Geophysical Research, 1041: 15365–15390.

    Guo F, Fan W M, Li C W, Zhao L, Li H X and Yang J H. 2012. Multi-stage crust-mantle interaction in SE China: Temporal, thermal and compositional constraints from the Mesozoic felsic volcanic rocks in eastern Guangdong-Fujian provinces. Lithos, 150: 62–84.

    Griffin W L, Belousova E A, Shee S R, Pearson N J and O’Reilly S Y. 2004. Archean crustal evolution in the northern Yilgarn Craton: U-Pb and Hf-isotope evidence from detrital zircons. Precambrian Research, 131: 231–282.

    Griffin W L, Pearson N J, Belousova E, Jackson S E, van Achterbergh E, O’Reilly S Y and Shee S R. 2000. The Hf isotope composition of cratonic mantle: LA-MCICPMS analysis of zircon megacrysts in kimberlites. Geochimica et Cosmochimica Acta, 64: 133–147.

    Hong D W, Wang S G, Han B F and Jin M Y. 1996. Post-orogenic alkaline granites from China and comparisons with anorogenic alkaline granites elsewhere. Journal of Southeast Asian Earth Sciences, 13: 13–27.

    Hoskin P W O and Schaltegger U. 2003. The composition of zircon and igneous and metamorphic petrogenesis. Reviews in Mineralogy and Geochemistry, 53: 27–62.

    Hu Z C, Liu Y S, Gao S, Liu W H, Zhang W, Tong X R, Lin L, Zong K Q, Li M, Chen H, Zhou L and Yang L. 2012. Improved in situ Hf isotope ratio analysis of zircon using newly designed X skimmer cone and jet sample cone in combination with the addition of nitrogen by laser ablation multiple collector ICP-MS. Journal ofAnalytical Atomic Spectrometry, 27: 1391.

    Jaeger J J, Courtillot V and Tapponnier P. 1989. Paleontological view of the ages of the Deccan Traps, the Cretaceous/Tertiary boundary, and the India-Asia collision. Geology, 17: 316–319.

    Jahn B M and Condie K C. 1995. Evolution of the kaapvaal craton as viewed from geochemical and smrnd isotopic analyses of intracratonic pelites. Geochimica et Cosmochimica Acta, 59: 2239–2258.

    Jahn B M. 1974. Mesozoic thermal events in southeast China. Nature, 248: 480-483.

    King P L, White A, Chappell B W and Allen C M. 1997. Characterization and origin of aluminous A-type granites from the Lachlan Fold Belt, southeastern Australia. Journal of Petrology, 38: 371–391.

    Li X H, Zhou H W, Chung S L, Ding S J, Liu Y, Lee C Y, Ge W C, Zhang Y M and Zhang R J. 2002. Geochemical and Sm-Nd isotopic characteristics of metabasites from central Hainan Island, South China and their tectonic significance. The Island Arc, 11: 193–205.

    Li Y J, Wei J H, Chen H Y, Tan J, Fu L B and Wu G. 2012. Origin of the Maoduan Pb-Zn-Mo deposit, eastern Cathaysia Block, China: Geological, geochronological, geochemical, and Sr-Nd-Pb-S isotopic constraints. Mineralium Deposita, 47: 763–780.

    Li Z X and Li X H. 2007. Formation of the 1300-km-wide intracontinental orogen and postorogenic magmatic province in Mesozoic South China: A flat-slab subduction model. Geology, 35: 179.

    Li Z X, Li X H, Chung S, Lo C H, Xu X S and Li W X. 2012. Magmatic switch-on and switch-off along the South China continental margin since the Permian: Transition from an Andean-type to a Western Pacific-type plate boundary. Tectonophysics, 532-535: 271–290.

    Li Z X, Li X H, Zhou H W and Kinny P D. 2002. Grenvillian continental collision in South China: New SHRIMP U-Pb zircon results and implications for the configuration of Rodinia. Geology, 30: 163–166.

    Liu Y S, Gao S, Hu Z C, Gao C G, Zong K Q and Wang D B. 2010. Continental and oceanic crust recycling-induced melt-peridotite interactions in the Trans-North China orogen: U-Pb dating, Hf isotopes and trace elements in zircons from mantle xenoliths. Journal of Petrology, 51: 537–571.

    Loiselle M C and Wones D R. 1979. Characteristics and origin of anorogenic granites. Geological of Society of America, 11: 468.

    Ludwig K R. 2003. User’s manual for Isoplot 3.00: A geochronological toolkit for Microsoft Excel. Berkeley Geochronology Center Special Publication.

    Maniar P D and Piccoli P M. 1989. Tectonic discrimination of granitoids. Geological Society of America Bulletin, 101: 635–643.

    Metcalfe I, Shergold J H and Li Z X. 1993. Gondwana dispersion and Asian accretion: Fieldwork on Hainan Island. Episodes, 16: 443–447.

    Pearce J A. 1996. Sources and settings of granitic rocks. Episodes, 19: 120–125.

    Ren J Y, Tamaki K, Li S and Junxia Z. 2002. Late Mesozoic and Cenozoic rifting and its dynamic setting in Eastern China and adjacent areas. Tectonophysics, 344: 175–205.

    Scherer E, Münker C and Mezger K. 2001. Calibration of the Lutetium-Hafnium Clock. Science, 293: 683–687.

    Steiger R H and J?ger E, 1977. Subcommission on geochronology: Convention on the use of decay constants in geochronology and cosmochronology. Earth and Planetary Science Letters, 36: 359–362.

    Sun S S and Mcdonough W F. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes. Geological Society, London, Special Publications, 42: 313–345.

    Tatsumi Y, Maruyama S and Nohda S. 1990. Mechanism of backarc opening in the Japan Sea: Role of asthenospheric injection. Tectonophysics, 181(1–4): 299–306.

    Turner S P, Foden J D and Morrison R S. 1992. Derivation of some A-type magmas by fractionation of basaltic magma: An example from the Padthaway Ridge, South Australia. Lithos, 28: 151–179.

    Vervoort J D and Blichert-Toft J. 1999. Evolution of the depleted mantle: Hf isotope evidence from juvenile rocks through time. Geochimica et Cosmochimica Acta, 63: 533–556.

    Wang Y J, Fan W M, Cawood P A and Li S Z. 2008. Sr-Nd-Pb isotopic constraints on multiple mantle domains for Mesozoic mafic rocks beneath the South China Block hinterland. Lithos, 106: 297–308.

    Wang Y J, Fan W M, Zhang G W and Zhang Y H. 2013. Phanerozoic tectonics of the South China Block: Key observations and controversies. Gondwana Research, 23: 1273–1305.

    Wang Z L, Xu D R, Zhang Z C, Zou F H, Wang L, Yu L L and Hu M Y. 2015. Mineralogy and trace element geochemistry of the Co-and Cu-bearing sulfides from the Shilu Fe-Co-Cu ore district in Hainan Province of South China. Journal of Asian Earth Sciences, 113:980-997.

    Watson E B, Wark D A and Thomas J B. 2006. Crystallization thermometers for zircon and rutile. Contributions to Mineralogy and Petrology, 151: 413–433.

    Whalen J B, Currie K L and Chappell B W. 1987. A-type granites: Geochemical characteristics, discrimination and petrogenesis. Contributions to Mineralogy and Petrology, 95: 407–419.

    Wu Y H and Zheng Y F. 2004. Genesis of zircon and its constraints on interpretation of U-Pb age. Chinese Science Bulletin, 49: 1554–1569.

    Xu D R, Wang Z L, Cai J X, Wu C J, Bakun-Czubarow N, Wang L, Chen H Y, Baker M J and Kusiak M A. 2013. Geological characteristics and metallogenesis of the shilu Fe-ore deposit in Hainan Province, South China. Ore Geology Reviews, 53: 318–342.

    Yang J H, Wu F Y, Chung S L, Wilde S A and Chu M F. 2006. A hybrid origin for the Qianshan A-type granite, northeast China: Geochemical and Sr-Nd-Hf isotopic evidence. Lithos, 89: 89–106.

    Zhang F F, Wang Y J, Chen X Y, Fan W M, Zhang Y H, Zhang G W and Zhang A M. 2011. Triassic high-strain shear zones in Hainan Island (South China) and their implications on the amalgamation of the Indochina and South China Blocks: Kinematic and40Ar/39Ar geochronological constraints. Gondwana Research, 19: 910–925.

    Zhao Z F, Dai L Q and Zheng Y F. 2013. Postcollisional mafic igneous rocks record crust-mantle interaction during continental deep subduction. Scientific Reports, 3(12): 3413–3413.

    Zhou X M and Li W X. 2000. Origin of Late Mesozoic igneous rocks in Southeastern China: Implications for lithosphere subduction and underplating of mafic magmas. Tectonophysics, 326(3–4): 269–287.

    Zircon U-Pb Age and Geochemical and Nd-Hf Isotopic Constraints on the Origin of the Early Cretaceous Shilu A-type Granite in Hainan Island, South China

    LI Xiang, WEI Junhao*, LI Yanjun and ZHAI Yulin
    (Faculty of Earth Resources, China University of Geosciences (Wuhan), Wuhan 430074, Hubei, China)

    Mesozoic granites widely distribute in Hainan island, South China. However, U-Pb age, geochemistry and Nd-Hf isotopes of the Early Cretaceous A-type granites have rarely been reported so far. This paper reports a concealed A-type granite, mainly composed of biotite monzogranite, found in drill cores from the Shilu iron deposit (it is called the Shilu A-type granite). LA-ICP-MS zircon U-Pb dating of the A-type granite yields a weighted206Pb/238U mean age of 133±3 Ma (MSWD=0.33), indicating that it was emplaced in the Early Cretaceous. Petrography and geochemical compositions show that the granite belongs to high-K cal-alkaline series and is characterized by high silicon (SiO2=75.22%–76.63%), alkali (K2O+Na2O=8.01%–8.56%) and iron (FeOT=1.62%–2.01%), but low calcium (CaO=0.55%–0.69%) and magnesium (MgO=0.03%–0.15%). The granite is metaluminous to weak peraluminous with A/CNK values in the range of 0.99 to 1.04. The granitic rocks are enriched in large ion lithophile element (e.g. K and Rb) and high field strength elements (e.g. U and Th) but depleted in Ba, Sr. The chondrite normalized REE patterns of the granitic rocks are right dipping with significant negative Eu anomalies (δEu=0.02- 0.08). The evidence from mineralogy and geochemistry suggests that the concealed granite is a highly fractionated aluminous A-type granite. The εNd(t) values range from –6.3 to –7.1 with two-stage model ages tDM2(Nd) from 1.44 to 1.47 Ga. The εHf(t) values of zircon from the A-type granite range from –18.5 to –21.7 and the two-stage model ages tDM2(Hf) range from 1.93 Ga to 2.09 Ga, which suggests that the granite was dominantly derived from the basement with minor enriched mantle materials. The trace element features show that the A-type granite was formed in an extensional tectonic setting in the Early Cretaceous, which was related to the northward movement of the Indian plate and the atteunuation of the subduction of the Pacific plate.

    A-type granite; Early Cretaceous; petrogenesis; extensional tectonic setting; Hainan

    P595; P597

    A

    1001-1552(2016)03-0587-016

    2015-07-22; 改回日期: 2015-10-18
    項目資助: 國家自然科學基金項目(41202054)、中國地質大學(武漢)中央高校基本科研業(yè)務費專項基金(CUG120702)和全國危機礦山“桂東粵西地區(qū)鉛鋅金等礦床成礦規(guī)律總結研究”項目(20089946)聯(lián)合資助。

    李翔(1990–), 男, 碩士研究生, 礦產普查與勘探專業(yè)。Email: lixiang9007@163.com

    魏俊浩(1961–), 教授, 博士生導師, 從事礦床地球化學、成礦規(guī)律與成礦預測研究。Email: junhaow@163.com

    猜你喜歡
    鋯石A型同位素
    鋯石成因礦物學及Lu-Hf同位素的應用
    俄成功試射“鋯石”高超音速巡航導彈
    軍事文摘(2020年24期)2020-02-06 05:56:36
    鋯石微區(qū)原位U-Pb定年的測定位置選擇方法
    華北地質(2015年3期)2015-12-07 05:13:23
    深空探測用同位素電源的研究進展
    電源技術(2015年7期)2015-08-22 08:49:04
    DF100A型發(fā)射機馬達電源板改進
    新聞傳播(2015年6期)2015-07-18 11:13:15
    A型肉毒素在注射面部皺紋中的應用及體會
    A型肉毒毒素聯(lián)合減張壓迫法在面部整形切口的應用
    AZA型號磨齒機工件主軸的改造
    鉆石與鋯石的區(qū)別知多少?
    《同位素》(季刊)2015年征訂通知
    同位素(2014年3期)2014-06-13 08:22:28
    久久青草综合色| 国产成人一区二区在线| 91精品伊人久久大香线蕉| av天堂久久9| 人人澡人人妻人| 精品少妇内射三级| 99久久精品国产国产毛片| 亚洲精品成人av观看孕妇| 国产日韩欧美亚洲二区| 久久av网站| 三级国产精品欧美在线观看| 亚洲欧美精品专区久久| √禁漫天堂资源中文www| 国精品久久久久久国模美| 久久久久国产精品人妻一区二区| 日日啪夜夜爽| 黑丝袜美女国产一区| 深夜a级毛片| 国产男女内射视频| 久久久久精品性色| 欧美一级a爱片免费观看看| 五月开心婷婷网| 日产精品乱码卡一卡2卡三| 亚洲av男天堂| 午夜福利视频精品| 老熟女久久久| 大香蕉久久网| 日韩不卡一区二区三区视频在线| 毛片一级片免费看久久久久| 蜜桃久久精品国产亚洲av| 色5月婷婷丁香| 欧美精品人与动牲交sv欧美| 男女边吃奶边做爰视频| 人妻少妇偷人精品九色| 成年女人在线观看亚洲视频| 亚洲av日韩在线播放| 亚洲国产av新网站| 丰满少妇做爰视频| 国产亚洲av片在线观看秒播厂| 99re6热这里在线精品视频| 国产精品久久久久久久电影| 丝袜喷水一区| 99九九在线精品视频 | av免费观看日本| 午夜激情福利司机影院| 精品久久国产蜜桃| 国产av码专区亚洲av| 日本91视频免费播放| videossex国产| 久久久久久久久大av| 人人妻人人澡人人爽人人夜夜| 一个人免费看片子| 日韩 亚洲 欧美在线| 国产成人aa在线观看| 欧美精品高潮呻吟av久久| 黄色日韩在线| 亚洲av男天堂| 色婷婷久久久亚洲欧美| 美女中出高潮动态图| 国产精品久久久久久久电影| 国产亚洲91精品色在线| 欧美日韩av久久| 成人美女网站在线观看视频| 男的添女的下面高潮视频| 大码成人一级视频| 熟女人妻精品中文字幕| 少妇的逼水好多| 中文天堂在线官网| 精品一区在线观看国产| 天堂俺去俺来也www色官网| 少妇的逼水好多| 精品人妻偷拍中文字幕| 在线观看国产h片| 亚洲欧美一区二区三区国产| 亚洲无线观看免费| 国产熟女午夜一区二区三区 | 少妇人妻精品综合一区二区| a级片在线免费高清观看视频| 制服丝袜香蕉在线| 国产精品三级大全| 久久99热6这里只有精品| 99久久精品一区二区三区| 国产精品一区二区性色av| 色吧在线观看| 精品一区二区三区视频在线| 91精品一卡2卡3卡4卡| freevideosex欧美| 免费观看性生交大片5| 在线观看www视频免费| 国产精品国产av在线观看| 国内少妇人妻偷人精品xxx网站| 国产男女内射视频| 青春草国产在线视频| 免费看不卡的av| 日韩一本色道免费dvd| 精品久久久久久久久av| a级毛色黄片| 色5月婷婷丁香| av播播在线观看一区| 久久久久人妻精品一区果冻| 国产精品久久久久久久久免| 久久国产精品男人的天堂亚洲 | 日本黄色日本黄色录像| 99热国产这里只有精品6| 久久久精品94久久精品| 韩国高清视频一区二区三区| 大香蕉久久网| av又黄又爽大尺度在线免费看| 国产亚洲精品久久久com| 精品久久久噜噜| 夜夜看夜夜爽夜夜摸| av免费在线看不卡| 国产免费一区二区三区四区乱码| 黄色欧美视频在线观看| 久久热精品热| 日韩成人av中文字幕在线观看| 另类精品久久| 欧美日韩综合久久久久久| 亚洲综合色惰| 久久久久精品性色| 不卡视频在线观看欧美| 我的女老师完整版在线观看| 国产精品一区www在线观看| 精品一品国产午夜福利视频| 国产精品无大码| 日本爱情动作片www.在线观看| 特大巨黑吊av在线直播| 国产有黄有色有爽视频| 欧美另类一区| 狠狠精品人妻久久久久久综合| 成人免费观看视频高清| 国产亚洲91精品色在线| 欧美丝袜亚洲另类| 久久国产精品大桥未久av | 人人妻人人看人人澡| 男人舔奶头视频| 美女xxoo啪啪120秒动态图| 免费看av在线观看网站| 老司机影院毛片| 亚洲精品色激情综合| 色婷婷久久久亚洲欧美| 内射极品少妇av片p| 少妇被粗大猛烈的视频| 2021少妇久久久久久久久久久| 黄色欧美视频在线观看| 亚洲熟女精品中文字幕| 我要看黄色一级片免费的| 久久精品熟女亚洲av麻豆精品| 国产永久视频网站| 精品少妇久久久久久888优播| 欧美日韩亚洲高清精品| 免费观看在线日韩| 久久久久人妻精品一区果冻| 十八禁高潮呻吟视频 | 另类精品久久| 久久人妻熟女aⅴ| 国产精品欧美亚洲77777| 97精品久久久久久久久久精品| 赤兔流量卡办理| 国产精品福利在线免费观看| 能在线免费看毛片的网站| 亚洲第一av免费看| 亚洲精品日韩av片在线观看| 国内精品宾馆在线| 内射极品少妇av片p| 国产一区二区在线观看av| 夜夜骑夜夜射夜夜干| 久久婷婷青草| 国产极品粉嫩免费观看在线 | 国产日韩欧美视频二区| 亚洲中文av在线| 99热网站在线观看| 青春草国产在线视频| 少妇高潮的动态图| 亚洲内射少妇av| 欧美3d第一页| 插逼视频在线观看| 少妇裸体淫交视频免费看高清| 新久久久久国产一级毛片| 天天躁夜夜躁狠狠久久av| 久久久久久久精品精品| 免费黄频网站在线观看国产| 一区在线观看完整版| 中文字幕精品免费在线观看视频 | 男人添女人高潮全过程视频| 少妇丰满av| 一级a做视频免费观看| 色视频在线一区二区三区| 亚洲欧美中文字幕日韩二区| 蜜桃在线观看..| 亚洲美女视频黄频| 麻豆精品久久久久久蜜桃| 三级经典国产精品| 国产精品一区二区在线不卡| 国产色爽女视频免费观看| 亚洲真实伦在线观看| 精品卡一卡二卡四卡免费| 精品一品国产午夜福利视频| 九九久久精品国产亚洲av麻豆| 黄色视频在线播放观看不卡| 美女脱内裤让男人舔精品视频| 国国产精品蜜臀av免费| av黄色大香蕉| 国产熟女午夜一区二区三区 | 大又大粗又爽又黄少妇毛片口| 国产精品成人在线| 色婷婷久久久亚洲欧美| 国产高清不卡午夜福利| 亚洲情色 制服丝袜| 又爽又黄a免费视频| freevideosex欧美| 一本大道久久a久久精品| 日日摸夜夜添夜夜添av毛片| 大片电影免费在线观看免费| 极品少妇高潮喷水抽搐| 五月伊人婷婷丁香| 国产精品一区www在线观看| 国产在线免费精品| 亚洲av电影在线观看一区二区三区| 亚洲精品乱码久久久v下载方式| 91午夜精品亚洲一区二区三区| 91精品国产国语对白视频| 成年女人在线观看亚洲视频| 99久久精品热视频| 国产精品久久久久久久久免| 国产精品国产三级专区第一集| 男男h啪啪无遮挡| av视频免费观看在线观看| 日本wwww免费看| 亚洲人与动物交配视频| 久久人妻熟女aⅴ| 久久久久久久精品精品| 亚洲欧美成人精品一区二区| av又黄又爽大尺度在线免费看| 晚上一个人看的免费电影| 国产成人精品无人区| 亚洲欧美日韩卡通动漫| 中国美白少妇内射xxxbb| 久久久久视频综合| 亚洲精品乱码久久久v下载方式| 男人舔奶头视频| 性色avwww在线观看| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品一区蜜桃| 日本色播在线视频| 校园人妻丝袜中文字幕| 九九爱精品视频在线观看| 伊人亚洲综合成人网| videossex国产| 高清视频免费观看一区二区| a级毛片免费高清观看在线播放| 亚洲精品国产av蜜桃| av网站免费在线观看视频| 王馨瑶露胸无遮挡在线观看| www.av在线官网国产| 国国产精品蜜臀av免费| 久久狼人影院| 日韩不卡一区二区三区视频在线| 香蕉精品网在线| 久久久亚洲精品成人影院| 一级毛片我不卡| 中文资源天堂在线| 777米奇影视久久| 亚洲人成网站在线播| 男女啪啪激烈高潮av片| 亚洲精品日本国产第一区| 日韩熟女老妇一区二区性免费视频| 噜噜噜噜噜久久久久久91| 久久久a久久爽久久v久久| 日本午夜av视频| 最近的中文字幕免费完整| 一级毛片黄色毛片免费观看视频| 色视频在线一区二区三区| 午夜免费观看性视频| 日日撸夜夜添| 亚洲国产精品国产精品| 一区二区av电影网| 亚洲国产成人一精品久久久| 欧美精品一区二区免费开放| 国产黄色视频一区二区在线观看| 高清黄色对白视频在线免费看 | av.在线天堂| 9色porny在线观看| 欧美高清成人免费视频www| 精品午夜福利在线看| 99久久人妻综合| 亚洲熟女精品中文字幕| 国产无遮挡羞羞视频在线观看| av网站免费在线观看视频| 最近手机中文字幕大全| 色网站视频免费| 欧美日韩视频高清一区二区三区二| 亚洲丝袜综合中文字幕| 国产成人精品婷婷| videos熟女内射| 26uuu在线亚洲综合色| 人妻系列 视频| 纵有疾风起免费观看全集完整版| 五月天丁香电影| 国产精品一区二区在线观看99| 天美传媒精品一区二区| 中文资源天堂在线| 亚洲国产精品一区三区| 日日摸夜夜添夜夜添av毛片| 青春草亚洲视频在线观看| 欧美97在线视频| 七月丁香在线播放| 丁香六月天网| 亚洲精品乱久久久久久| 欧美精品高潮呻吟av久久| 亚洲美女视频黄频| 99久久人妻综合| 欧美bdsm另类| 纯流量卡能插随身wifi吗| 91aial.com中文字幕在线观看| 国产 精品1| av.在线天堂| 精品国产国语对白av| 国产伦在线观看视频一区| 国产视频首页在线观看| 97超碰精品成人国产| 国产男女内射视频| 99久久精品一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品国产成人久久av| 大又大粗又爽又黄少妇毛片口| 久久久久视频综合| 夜夜爽夜夜爽视频| 亚洲av福利一区| 国产色爽女视频免费观看| www.色视频.com| 亚洲国产成人一精品久久久| 九九在线视频观看精品| 国产成人免费观看mmmm| 好男人视频免费观看在线| 国产91av在线免费观看| 免费大片黄手机在线观看| 最近最新中文字幕免费大全7| 男人爽女人下面视频在线观看| 看免费成人av毛片| 日日爽夜夜爽网站| 亚洲中文av在线| 啦啦啦中文免费视频观看日本| 91精品一卡2卡3卡4卡| 国国产精品蜜臀av免费| 亚洲人成网站在线观看播放| 国产精品欧美亚洲77777| 欧美97在线视频| 亚洲精品色激情综合| 乱码一卡2卡4卡精品| 99久久精品热视频| 极品少妇高潮喷水抽搐| 成人国产av品久久久| 97超碰精品成人国产| 成人特级av手机在线观看| 国产伦精品一区二区三区四那| 高清av免费在线| 男人添女人高潮全过程视频| 欧美国产精品一级二级三级 | 亚洲精品久久午夜乱码| 欧美+日韩+精品| 久久人人爽av亚洲精品天堂| 国产一区二区三区综合在线观看 | 亚洲精品国产色婷婷电影| 另类亚洲欧美激情| 亚洲av电影在线观看一区二区三区| 亚洲精华国产精华液的使用体验| 9色porny在线观看| 免费大片黄手机在线观看| av在线播放精品| 九草在线视频观看| 美女福利国产在线| 全区人妻精品视频| 男人爽女人下面视频在线观看| 国产免费一级a男人的天堂| 在现免费观看毛片| 男女啪啪激烈高潮av片| 建设人人有责人人尽责人人享有的| 国产亚洲5aaaaa淫片| 韩国av在线不卡| 国产亚洲5aaaaa淫片| 久久精品久久久久久久性| 日产精品乱码卡一卡2卡三| 国产精品国产三级国产专区5o| 只有这里有精品99| 女的被弄到高潮叫床怎么办| 日韩不卡一区二区三区视频在线| 中文在线观看免费www的网站| 成人国产av品久久久| 黄色日韩在线| 国产精品久久久久成人av| 欧美xxⅹ黑人| 51国产日韩欧美| 亚洲一区二区三区欧美精品| 男人狂女人下面高潮的视频| 日本与韩国留学比较| 亚洲精品中文字幕在线视频 | 久久久久国产精品人妻一区二区| 国产精品成人在线| 国产爽快片一区二区三区| 视频区图区小说| 亚洲图色成人| 日韩三级伦理在线观看| 青春草视频在线免费观看| 亚洲欧洲日产国产| 日日摸夜夜添夜夜添av毛片| 亚洲欧美日韩卡通动漫| 九九在线视频观看精品| 国产精品久久久久久精品电影小说| 久久免费观看电影| 国产亚洲一区二区精品| 久久精品久久久久久噜噜老黄| 自线自在国产av| 久久久久久久久大av| 国产高清三级在线| 国产精品三级大全| 欧美激情国产日韩精品一区| 大片免费播放器 马上看| 国产精品久久久久久久电影| 日韩 亚洲 欧美在线| 一区二区三区精品91| 欧美 日韩 精品 国产| 黄色一级大片看看| 亚洲欧美清纯卡通| 国内少妇人妻偷人精品xxx网站| 少妇熟女欧美另类| 精品国产乱码久久久久久小说| 久久精品国产a三级三级三级| 一本—道久久a久久精品蜜桃钙片| 国产成人精品一,二区| 久久久久视频综合| 国产伦精品一区二区三区四那| 亚洲av电影在线观看一区二区三区| 国产成人freesex在线| 亚洲成人手机| 亚洲欧美一区二区三区黑人 | 免费播放大片免费观看视频在线观看| 观看免费一级毛片| 亚洲欧洲日产国产| 一个人看视频在线观看www免费| 偷拍熟女少妇极品色| 黄色一级大片看看| 国产亚洲精品久久久com| a级毛片在线看网站| 欧美3d第一页| 美女国产视频在线观看| 精品一品国产午夜福利视频| 视频中文字幕在线观看| 色网站视频免费| 麻豆成人午夜福利视频| 免费高清在线观看视频在线观看| 晚上一个人看的免费电影| 亚洲综合色惰| 成人午夜精彩视频在线观看| 成人漫画全彩无遮挡| 黄色配什么色好看| 99久国产av精品国产电影| 亚洲精品自拍成人| 亚洲人成网站在线播| 一本色道久久久久久精品综合| 男女无遮挡免费网站观看| www.色视频.com| 又爽又黄a免费视频| 一级片'在线观看视频| 亚洲欧美成人精品一区二区| 少妇猛男粗大的猛烈进出视频| 亚洲电影在线观看av| 熟女人妻精品中文字幕| 最新的欧美精品一区二区| 麻豆成人午夜福利视频| 午夜福利网站1000一区二区三区| 少妇熟女欧美另类| 精品国产一区二区久久| 一级毛片aaaaaa免费看小| 青青草视频在线视频观看| 欧美最新免费一区二区三区| 国产 精品1| 日本色播在线视频| 国产一区有黄有色的免费视频| 多毛熟女@视频| 日本黄色日本黄色录像| 国产精品久久久久成人av| 最新的欧美精品一区二区| 纯流量卡能插随身wifi吗| 亚洲精品乱久久久久久| 日韩成人伦理影院| 成人无遮挡网站| 亚洲精品aⅴ在线观看| 性色av一级| 亚洲国产精品专区欧美| 少妇被粗大的猛进出69影院 | 久久久久精品久久久久真实原创| 成人综合一区亚洲| 色视频在线一区二区三区| 最近中文字幕2019免费版| 日韩 亚洲 欧美在线| 69精品国产乱码久久久| 亚洲第一av免费看| 黄色怎么调成土黄色| 全区人妻精品视频| 久久99热6这里只有精品| 精品久久久久久电影网| 成人午夜精彩视频在线观看| 在线观看免费高清a一片| 在线看a的网站| 日本黄色片子视频| 美女国产视频在线观看| 蜜桃在线观看..| 王馨瑶露胸无遮挡在线观看| 亚洲精品亚洲一区二区| 久久精品国产亚洲av天美| 精品久久久久久久久av| www.av在线官网国产| 日韩亚洲欧美综合| 国产真实伦视频高清在线观看| 91成人精品电影| 中国国产av一级| 久久久久久久久久人人人人人人| 亚洲欧美日韩另类电影网站| www.色视频.com| 免费看不卡的av| 久久午夜福利片| 亚洲av中文av极速乱| 乱人伦中国视频| 激情五月婷婷亚洲| 国产黄片美女视频| 免费黄色在线免费观看| 亚洲国产精品专区欧美| 97超碰精品成人国产| 丰满少妇做爰视频| 自拍偷自拍亚洲精品老妇| 国产日韩欧美亚洲二区| 成人18禁高潮啪啪吃奶动态图 | 日日啪夜夜撸| 涩涩av久久男人的天堂| 丰满迷人的少妇在线观看| 国产av码专区亚洲av| 天堂俺去俺来也www色官网| 成人国产av品久久久| 精品人妻偷拍中文字幕| 亚洲精品乱码久久久v下载方式| 欧美国产精品一级二级三级 | 国产美女午夜福利| 国内揄拍国产精品人妻在线| 少妇的逼水好多| 久久午夜福利片| 免费高清在线观看视频在线观看| 老女人水多毛片| 丰满迷人的少妇在线观看| 夜夜看夜夜爽夜夜摸| 亚洲av综合色区一区| 亚州av有码| 欧美+日韩+精品| 在线 av 中文字幕| 视频区图区小说| 午夜免费男女啪啪视频观看| 久久午夜福利片| 一区二区三区精品91| 国产高清三级在线| 久久97久久精品| 日韩人妻高清精品专区| 黑人巨大精品欧美一区二区蜜桃 | 一级片'在线观看视频| 国产成人精品福利久久| 少妇人妻精品综合一区二区| av免费在线看不卡| 一个人看视频在线观看www免费| 日韩亚洲欧美综合| 最新中文字幕久久久久| 两个人免费观看高清视频 | 精品人妻偷拍中文字幕| 亚洲国产精品国产精品| 人人澡人人妻人| 国产女主播在线喷水免费视频网站| 熟妇人妻不卡中文字幕| 日本猛色少妇xxxxx猛交久久| 亚洲av二区三区四区| 青春草国产在线视频| 国产精品人妻久久久影院| 国产在线免费精品| 秋霞在线观看毛片| 亚洲精品色激情综合| 国产在线免费精品| 免费在线观看成人毛片| 日韩大片免费观看网站| 国产在线视频一区二区| 亚洲人与动物交配视频| 十分钟在线观看高清视频www | 欧美少妇被猛烈插入视频| 欧美区成人在线视频| 嘟嘟电影网在线观看| 国产成人精品婷婷| 国产精品一区二区在线观看99| 丰满饥渴人妻一区二区三| 又粗又硬又长又爽又黄的视频| 国产黄片美女视频| 久久精品久久久久久久性| 国产国拍精品亚洲av在线观看| .国产精品久久| 啦啦啦啦在线视频资源| 免费久久久久久久精品成人欧美视频 | 成人国产av品久久久| 午夜av观看不卡| 国产淫语在线视频| 成人无遮挡网站| av免费观看日本| 国产又色又爽无遮挡免| 日本欧美国产在线视频| 日韩不卡一区二区三区视频在线| 亚洲无线观看免费| 麻豆成人午夜福利视频| 美女国产视频在线观看| 人妻制服诱惑在线中文字幕| 少妇熟女欧美另类| 久久人人爽人人爽人人片va|