——以湖南蓮花洞LHD5石筍為例"/>
999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

全新世早期弱夏季風事件的精確定位及機制探討
——以湖南蓮花洞LHD5石筍為例

2016-07-15 06:35:14張華生殷建軍程海RLawrenceEdwards林玉石唐偉楊會涂林玲王華潘謀成吳
沉積學報 2016年2期

張華生殷建軍程 海R Lawrence Edwards林玉石唐 偉楊 會涂林玲王 華潘謀成吳 夏

?

全新世早期弱夏季風事件的精確定位及機制探討
——以湖南蓮花洞LHD5石筍為例

張華生1,2殷建軍2程海3,4RLawrence Edwards4林玉石2唐偉2楊會2涂林玲2王 華2潘謀成2吳夏2

(1.西南大學地理科學學院 重慶 400715;2.國土資源部/廣西巖溶動力學重點實驗室國際巖溶研究中心中國地質科學院巖溶地質研究所 桂林 541004;3.西安交通大學全球環境變化研究院 西安 710054;4.Department of Earth Science,University of Minnesota Minneapolis 55455)

摘 要全新世早期是太陽輻射加強、全球溫度上升,并伴隨著冰蓋消融的重要時期,而其間發生的冷事件以及亞洲季風區的弱夏季風事件的成因一直是全新世早期研究的重點。對亞洲季風—海洋—極地聯系研究有著重要的意義。通過分析湖南蓮花洞LHD5石筍28個U/Th年齡和535個氧同位素數據重建了全新世亞洲季風演化特征,其中全新世早期分辨率達8年。LHD5石筍記錄到YD結束時間為11 748±30 a B.P.,全新世開始于11 684±39 a B.P.,轉換時間約為64年,與格陵蘭gicc05記錄在誤差范圍內一致。LHD5石筍記錄到全新世早期6次弱夏季風事件,事件年齡中心點分別為11 461±34 a B.P.、10 354±36 a B.P.、9 957±25 a B.P.、9 062±36 a B.P.、8 744±23 a B.P.、8 144±24 a B.P.,其δ18O值的波動幅度分別為1.08‰、0.94‰、0.66‰、0.90‰、0.55‰、1.02‰,這些弱季風事件在亞洲季風區具有普遍的區域意義。除8.2 ka事件之外,10 ka B.P.之前的弱季風事件除了受到太陽活動的影響,還受到北大西洋IRD事件的影響,而之后更多地受到太陽活動和ITCZ南移的影響。

關鍵詞全新世早期 弱夏季風事件 石筍 δ18O 蓮花洞

0 引言

許多學者通過研究全新世早期各種地質記錄,已經初步認識了全新世早期亞洲季風變化規律[1-18]。全新世早期總體呈現季風增強或強季風期[1-9],并認為全新世早期亞洲季風增強主要受歲差驅動的太陽輻射的影響[10],全新世早期北半球太陽輻射增強推動熱帶輻合(ITCZ)北移,促使亞洲季風增強[11-12]。然而,全新世早期亞洲季風在逐漸增強的趨勢下存在弱季風期[5,8,13-14],多數學者認為全新世早期千年尺度的氣候振蕩與北大西洋冰漂碎屑(IRD)有關[5,11-12,14-15],將亞洲季風區千年尺度的季風變化與北半球高緯地區冰量變化聯系起來,關于短時間尺度的氣候波動則主要受太陽活動影響[3,5,8,16-18]。前人的研究多是覆蓋整個全新世的亞洲季風演變,利用高分辨率石筍記錄討論全新世早期季風演化的研究較少,例如,楊琰等[8]利用衙門洞Y1石筍討論了全新世早期的季風變化并認為短時間尺度的季風變化受太陽活動影響。Cai et al.[18]利用東石崖DSY09石筍討論了中國北方全新世早期東亞季風弱季風事件,認為太陽活動與高緯冰量變化共同影響東亞季風演變。Zhang et al.[12]利用蓮花洞LH2石筍討論了整個全新世的季風演化特征,但LH2石筍底部有暗色層可能存在沉積間斷,導致全新世早期石筍分辨率較低,不足以揭示全新世早期亞洲季風演變的細節特征。另外,Wang et al.[17]利用董哥洞DA石筍證實了全新世9 ka B.P.以來百年—十年尺度亞洲季風受太陽活動的影響,但是卻缺少9 ka B.P.之前的全新世早期,那么全新世早期百年—十年尺度的亞洲季風變化與太陽活動又有怎樣的關系?本文利用受東亞季風和印度季風共同控制的蓮花洞LHD5石筍,討論全新世早期亞洲季風演變特征。LHD5石筍在全新世早期具有平均8年的分辨率,能夠更好地揭示全新世早期亞洲季風的細節特征,年代際尺度的弱季風事件有無區域意義以及弱季風事件的驅動機制,這些都是值得探討的。

1 研究區概況

本文利用LHD5石筍重建古氣候的樣品,采自于湖南省龍山縣蓮花洞(109°32′E,29°29′N,圖1),該洞位于云貴高原和湖南盆地傾斜地帶,其特點為多梯級喀斯特高原。蓮花洞洞口海拔455m,洞長580m,高10~20 m,寬11~36 m,圍巖主要為奧陶系上統(O3)白云質灰巖。洞穴沿斷裂裂隙發育,組成北東、北西走向轉折的洞穴系統,底板由鈣華板組成[19]。蓮花洞所處地區屬亞熱帶大陸性濕潤季風氣候,距離蓮花洞約30 km洛塔鎮現代氣象資料記錄,該地多年平均降雨1 694 mm,多年平均氣溫14.3℃。東亞季風和印度季風共同控制著當地的氣候,每年4—9月主要受濕熱的夏季風控制,降雨豐沛,10月至次年3月主要受西伯利亞干冷冬季風控制,降水較少。

LHD5石筍長1 252mm,呈圓錐形,底部直徑151 mm,頂部直徑75 mm,為文石石筍。

圖1 蓮花洞位置圖。實心正方形代表蓮花洞(109°32′E,29°29′N)位置,實心圓代表其他洞穴位置,箭頭代表現代夏季風包括東亞夏季風和印度夏季風,虛線代表夏季風的北緣。Fig.1 The location of Lianhua Cave.Solid square indicates the location of Lianhua Cave(109°32′E,29°29′N),Solid circles indicate the location of other Caves.Arrows indicate generalized modern summermonsoonal winds including the East Asian and Indian summermonsoons.The dashed line represents the approximate summermonsoon limit at present.

2 分析方法

從頂部向下沿LHD5石筍生長軸采集石筍年齡樣品共28個(圖2)。年齡測試在美國明尼蘇達大學地球科學系和西安交通大學同位素實驗室完成,采用的儀器是Thermo Finnigan Neptune型多接收電感耦合等離子體質譜(MC-ICP-MS),采用的方法為230Th鈾系不平衡法測定。LHD5石筍共采集535個碳氧同位素樣品,1 252~900 mm每1 mm取一個樣品,900mm至頂部每5 mm取一個樣品。同位素樣品測試在中國地質科學院巖溶地質研究所巖溶地質與資源環境測試中心完成,實驗儀器為Thermo Finnigan MAT253同位素比值質譜儀,外部連接Kiel IV碳酸鹽自動進樣裝置,實驗室標準控制樣多次分析的δ18O值偏差<0.1‰,測試過程中每隔9個樣品插入一個國標GBW 04405進行檢驗,且每個國標后為2個平行樣品,平行樣 δ18O的差值<0.1‰,數據相對 VPDB給出。

3 結果

3.1 年代數據

表1及圖2呈現了LHD5石筍的U、Th同位素含量和28個230Th年齡結果,覆蓋年齡時段為11 918~85 a B.P.,其中11 918~8 030 a B.P.時段有13個年齡控制點,時標精度較高。LHD5石筍樣品238U含量高(1 464~6 610×10-9/g·g-1),232Th含量低(0.2~2.0×10-9/g·g-1),因此測年精度高,測年誤差為2~48年。

3.2 LHD5石筍δ18O記錄

LHD5石筍共采集535個氧同位素樣品,采用線性內插法獲取氧同位素年齡,平均分辨率為22年,其中11 984~9 019 a B.P.時段分辨率為8年。圖3中LHD5δ18O值變化范圍為-7.82‰~-4.67‰,變幅為3.15‰,平均值為-6.61‰,δ18O值變化特征非常明顯。12.0~11.75 ka B.P.δ18O值偏重,平均值為-5.18‰,11.75~11.68 ka B.P.δ18O值從-4.84‰迅速偏輕至-6.69‰,波動幅度達 1.85‰。11.68~9.0 ka B.P.δ18O值持續偏輕,9.0~7.0 ka B.P.δ18O值達到最輕,平均值為-7.26‰。11.68~8.0kaB.P.氧同.)位素持續變輕的時段內,存在6個氧同位素偏重階段。7.0~2.8 ka B.P.氧同位素值持續偏重,2.8~0.05 ka B.P.氧同位素波動較小且偏重,平均值為-5.79‰。

圖2 LHD5 230 Th年齡(a B.P.),Hendy Test取樣點及生長速率Fig.2 Simpling sites of230Th dating and Hendy Test,right figure shows the almost constant growth rates of the whole LHD5

表1 LHD5石筍U,Th同位素組成和230 Th年齡Tab le 1 Isotopic com positions of U,Th and230Th dating results

圖3 蓮花洞LHD5,LH2石筍[12]δ18O值,與65°N與30°N太陽輻射對比圖1,2,3,4,5分別代表YD事件,季風增強期,強季風期,季風減弱期,弱季風期Fig.3 A comparison of the LHD5 and LH2[12]δ18O record to the summer insolation for 65°N and 30°N 1,2,3,4,5 represent Younger Dryas event,Strengtheningmonsoon period,Strongmonsoon period,Weakeningmonsoon period,Weak monsoon period.

圖4 LHD5石筍Hendy檢驗結果4個生長層分別是1 164mm(實心正方形),986 mm(實心圓),490mm(叉號),110 mm(實心三角形)Fig.4 LHD5 stalagm ite Hendy Test resultsfour laminaes at depth of 1 164 mm(solid squares),986 mm(solid circles),490 mm(Cross),110 mm(solid triangles)were analyzed.

3.3 Hendy Test及重現性檢驗

影響石筍氧同位素值的因素,包括氣候因素(溫度,降水)和非氣候因素(動力分餾,蒸發)。Hendy[22]提出在同位素平衡分餾的情況下,石筍同位素受溫度和降雨的影響,即受氣候因素的影響。并提出檢驗石筍同位素平衡分餾的兩個準則:①同一生長層δ18O值相對穩定,向外無富集現象。②同一生長層的δ18O與δ13C不相關。LHD5石筍選取4個不同的生長層,每層7個樣點進行分析(圖2)。從圖4中看出氧同位素值比較恒定,碳同位素波動較大,碳氧的相關性不明顯,表明LHD5石筍在整個全新世達到了同位素平衡分餾。此外,本文又利用近年來較為普遍的氧同位素重現性檢驗方法[23],排除非氣候因子的影響。利用LHD5與同一洞穴LH2石筍作重現性檢驗,由于LH2石筍在全新世早期分辨率較低,本文又利用山寶洞SB10、SB34石筍及董哥洞D4石筍作檢驗,從圖5中可以看出5組氧同位素曲線對應良好,可以確認LHD5氧同位素代表氣候信號。

圖5 LHD5δ18O與LH2[12],山寶洞SB43和SB10[7],董哥洞D4[5]δ18 O值重現性檢驗Fig.5 The replication test ofδ18O record from stalagmite LHD5 and stalagmite LH2,SB43,SB10 from Shanbao Cave and D4 from Dongge cave

4 討論

4.1 12.0 ka B.P.以來季風演變

亞洲季風區大量的方解石石筍記錄表明,石筍δ18O值反映了夏季風的強弱變化,δ18O值偏輕代表夏季風增強,氣候溫暖濕潤,δ18O值偏重代表夏季風減弱,氣候干旱[3,8,10,17-18,24-26]。 前人對于蓮花洞文石石筍δ18O值的研究也已經證實百年尺度上石筍δ18O值代表亞洲季風強度[12,19,27],在年代際及以下尺度則反映區域降水量的變化,指示區域旱澇程度的變化[19]。

LHD5δ18O值記錄了12.0 ka B.P.以來亞洲季風強度變化(圖3)。LHD5δ18O值變化趨勢總體上與北半球太陽輻射一致,與65°N太陽輻射的吻合程度強于30°N,主要體現在LHD5δ18O峰值與65°N太陽輻射峰值的年代差較小,尤其是7 ka B.P.以后LHD5 δ18O值的變化幾乎與65°N太陽輻射一致。另外,LHD5δ18O值與亞洲季風區其他石筍δ18O值變化基本一致[3,5,7,11-12,18,28]。LHD5石筍δ18O值所反映的亞洲季風可以分為5個階段:12.0~11.75 ka B.P.為YD事件;11.75~9.0 ka B.P.季風增強期;9.0~7.0 ka B.P.強季風期;7.0~2.8 ka B.P.季風減弱期;2.8~0.05 ka B.P.弱季風期。同時,在季風總體演化趨勢上疊加了多次弱季風事件。同一洞穴LH2石筍討論了12.5 ka B.P.以來全新世亞洲季風演變,由于LH2石筍在全新世早期可能存在沉積間斷不足以揭示季風演變的細節特征,所以本文重點討論全新世早期(12~8.0 ka B.P.)亞洲季風演變。

4.2 全新世早期(12~8.0 ka B.P.)亞洲季風演變

LHD5石筍在12~8.0 ka B.P.期間反映的亞洲季風呈逐漸增強趨勢,具體又可分為3個階段:12.0~11.75 ka B.P.為YD事件結束期;11.75~11.68 ka B.P.為季風快速增強期;11.68~8.0 ka B.P.為季風緩慢增強期。LHD5石筍只記錄到YD事件的結束時間標定在11.75 ka B.P.,本文不對12.0~11.75 ka B.P.時段進行詳細的討論。

4.2.1 季風快速增強期(11.75~11.68 ka B.P.)

11.75~11.68 ka B.P.時段,LHD5石筍δ18O值從-4.84‰迅速偏輕至-6.69‰,波動幅度為1.85‰,為季風快速增強期,或YD事件向全新世轉換期。多數地質記錄把11.5 ka B.P.視為全新世起始點[5,8,18,24],而具有8年分辨率的LHD5石筍記錄到YD結束時間為11 748±30 a B.P.,全新世開始于11 684±39 a B.P.,轉換時間約為64年,與NGRIP冰芯(11 653± 99 a B.P.)極為一致(圖6),與GRIP冰芯(11 500±99 a B.P.)略有差異[29-30],這也充分說明全新世起始點的標定受測年誤差、氣候代用指標或局地環境等影響[24]。另外,董哥洞(108°5′E,25°17′N)、衙門洞(107°5′E,25°29′N)石筍δ18O值增強的趨勢不及LHD5石筍、NGRIP冰芯明顯,這種趨勢與緯度位置有關或受其他因素影響還有待進一步討論。

圖6 YD-全新世轉換期LHD5石筍與格陵蘭NGRIP冰芯的一致變化Fig.6 Synchronous variation ofδ18O record of NGRIP ice core and stalagmite LHD5 from YD to Early Holocence

4.2.2 季風緩慢增強期(11.68~8.0 ka B.P.)

11.68~8.0 ka B.P.時段,LHD5石筍δ18O值的波動范圍為-5.83‰ ~-7.67‰,為季風緩慢增強期。LHD5石筍δ18O值在9.3 ka B.P.開始達到峰值,與亞洲季風區其他石筍記錄一致。利用石筍δ18O值作為季風強度指標的亞洲季風最適宜期并未表現出穿時性,這與An的研究結論有所差異[31],可能的原因是氣候代用指標及其指標建立的基礎不同。該時段LHD5石筍在季風增強的趨勢上疊加了6次弱季風事件,事件中心年齡分別為11 461±34 a B.P.,10 354 ±36 a B.P.,9 957±25 a B.P.,9 062±36 a B.P.,8 744± 23 a B.P.,8 144±24 a B.P.。其中11 461±34 a B.P. 和8 144±24 a B.P.波動幅度超過1‰,其余弱季風事件波動幅度介于0.5‰~1‰之間。

LHD5石筍記錄的11 461±34 a B.P.弱季風事件是全新世早期波動最大的弱季風事件,其波動幅度為1.08‰,與其記錄的8 144±24 a B.P.弱季風事件的波動幅度(1.02‰)相當。11 461±34 a B.P.弱季風事件與NGRIP冰芯記錄的降溫事件(11 440±97 a B.P.)具有一致性(圖7b),同時該時段HGS濃度增加(表示存在冰漂碎屑事件)[32](圖7c),太陽黑子數量減少(代表弱的太陽活動)[33](圖7a)。亞洲季風區衙門洞Y1石筍(圖7f)、東石崖DSY09石筍[18]的弱季風事件分別記錄在11.35 ka B.P.、11.3 ka B.P.,在時間上存在100年的偏差。亞洲季風區其他地質記錄則不存在該弱季風事件,這就需要高分辨率地質記錄進行詳實的討論。LHD5石筍記錄的10 354±36 a B.P.弱季風事件,δ18O值波動幅度為0.94‰,特點為快速進入波動退出,事件持續事件約100年。該時段內太陽黑子數量明顯減少(圖 7a),北大西洋在 10.3 ka B.P.記錄了一次冰漂碎屑事件(圖7c)。亞洲季風區山寶洞石筍在這期間記錄了2次連續弱季風事件(圖7g),衙門洞[8]、東石崖[18]、老母洞[34]、Qunf洞[3]的弱季風事件分別記錄在10.2 ka B.P.、10.3 ka B.P.、10.2 ka B.P.、10.2 ka B.P.,這說明在10.2 ka B.P.左右亞洲季風區普遍存在弱季風事件。

LHD5石筍記錄的9 957±25 a B.P.弱季風事件,其δ18O值波動幅度為0.66‰,事件持續時長約100年。圖7中該時段太陽活動有一次較小的減弱期,NGRIP冰芯在9.84 ka B.P.存在一次較小的降溫事件,東石崖DSY09石筍在9.9 ka B.P.存在一次明顯的弱季風事件[18],其他地質記錄的代用指標波動較小或不存在該事件。LHD5石筍記錄的9 062±36 a B.P.弱季風事件δ18O值波動幅度為0.9‰,是一次比較強烈的弱季風事件,圖7中該時段太陽黑子數量減少,NGRIP冰芯在9 040±64 a B.P.也存在一次波動較小的降溫事件,其顯著降溫事件發生在9 260± 30 a B.P.,時間偏差較大。Rasmussen et al.[35]通過分析格陵蘭DYE-3,GRIP,NGRIP冰芯,得到9.35~9.24 ka B.P.降溫事件。Bond et al.[36]證實北大西洋在9.4 ka B.P.存在冰漂碎屑事件,王寧練在古里雅冰芯中也發現 9.4 ka B.P.冷事件[37]。亞洲季風區董哥洞[5]、衙門洞[8]、山寶洞[7]、東石崖[18]、Qunf洞[3]都在9.3~9.2 ka B.P.存在弱季風事件,足以說明在9.3 ~9.2 ka B.P.期間冷干事件發生區域是大范圍的。且董哥洞D4石筍和Qunf洞Q5石筍的9.3 ka B.P.弱季風事件強度大于其8.2 ka B.P.事件[3,5],并認為9.3 ka B.P.弱季風事件受北大西洋冰漂碎屑事件的影響。LHD5石筍在該時段沒有記錄到明顯的弱季風事件,原因可能是受測年誤差或局地氣候的影響。LHD5石筍記錄的8 744±23 a B.P.弱季風事件δ18O值波動幅度為0.55‰,與NGRIP冰芯降溫事件,弱的太陽活動具有一致性,只是各種指標的波動較小。Cariaco盆地在這期間記錄到一次較強的降水減少事件。亞洲季風區山寶洞石筍及湖光巖瑪珥湖記錄到8.6 ka B.P.弱季風事件[14,38],其他地質記錄在該時間段內不存在弱季風事件。

LHD5石筍記錄的8.2 ka B.P.弱季風事件突變于8 144±24 a B.P.,波動幅度為1.02‰。該時段正好處于弱的太陽活動期(圖7a),NGRIP冰芯降溫事件記錄在8 200±30 a B.P.(圖7b),亞洲季風區大多數地質記錄也識別到8.2 ka弱季風事件[3-5,8,12-14]。關于8.2 ka B.P.事件起止時間,內部結構特征的討論已有一定的認識,Cheng[39]利用亞洲季風區和南美季風區高分辨率的石筍記錄及格陵蘭冰芯記錄證實8.2 ka B.P.事件發生具有同步性,8.2 ka B.P.事件由8.21ka B.P.(持續時間70年)和8.06 ka B.P.(持續時間20年)事件構成,且8.2 ka B.P.事件時亞洲季風表現為弱季風而南美季風為強季風。LHD5石筍δ18O波動幅度為1.02‰,與亞洲季風區的Qunf洞Q5石筍[3],董哥洞DA、D4石筍記錄一致[5,17],但小于和尚洞HS4石筍及老母洞LM2石筍的波動幅度[28,34,40]。

圖7 全新世早期LHD5石筍δ18 O記錄與其他地質記錄對比a.太陽黑子數量[33];b.格陵蘭NGRIPδ18O記錄[29-30];c.北大西洋赤鐵礦濃度(%)[32];d.Cariaco盆地鈦濃[46];e.董哥洞D4δ18O記錄[5];f.衙門洞Y1與和尚洞HS4δ18O記錄[8,28];g.山寶洞SB43 與SB10δ18O記錄[7,14];h.LHD5δ18O記錄。Fig.7 A comparison of the LHD5 record from Lianhua Cave and other record in Early Holocenea.Sunspotnumber[33];b.Greenland NGRIP ice coreδ18O record[29-30];c.Hematite stained grains(%)from The North Atlantic Ocean[32];d.Titanium(%)data from Cariaco Basin[46];e.Dongge cave D4 stalagmite δ18O records[5];f.Yamen cave Y1 stalagmite and Heshang cave HS4δ18O records[8,28];g.Shanbao cave SB43 and SB10 stalagmiteδ18O records[7,14];h.LHD5 stalagmiteδ18O records.

以上討論得出,全新世早期LHD5石筍記錄的弱季風事件與NGRIP冰芯記錄的降溫事件一一對應,整個演變趨勢也極為一致,與弱的太陽活動也基本對應。LHD5石筍記錄的弱季風事件在亞洲季風區其他記錄中基本都存在,說明其記錄的弱季風事件具有區域意義。另外,這些弱季風(降溫)事件在年齡或指標波動幅度上存在差異,例如,圖3中LHD5石筍在全新世早期記錄的6次弱夏季風事件顯示的δ18O值均比現今氣候環境記錄的 δ18O值(-5.8‰~-6.2‰)偏輕0.4‰左后;而圖5中,董哥洞D4石筍在全新世早期記錄的弱夏季風事件所顯示的δ18O值均比現今氣候環境記錄的δ18O值偏重0.2‰~0.4‰。原因可能是蓮花洞現今氣候受太平洋近源水影響較大導致δ18O值偏重[41],而董哥洞在整個全新世受印度洋水汽影響[42]。另外,通過分析亞洲季風區所報道的石筍記錄[3,7,12,14]只有董哥洞D4石筍在全新世早期記錄的弱夏季風事件的δ18O值比現今氣候環境記錄的δ18O值偏重,其原因則可能是受局地氣候或其他因素影響。

5 全新世早期亞洲季風機制討論

全新世早期LHD5石筍記錄的亞洲季風呈逐漸增強的趨勢,與北半球太陽輻射一致,兩者的關系也得到前人的證實[4-5,10,17,24-25]。 然而,LHD5石筍記錄的季風最強盛期出現在9~7 ka B.P.,滯后于太陽輻射峰值約2 ka(圖3),7 ka B.P.以后二者則呈現一致性變化。另外,LHD5石筍與NGRIP冰芯在全新世早期變化趨勢一致,尤其是在YD事件向全新世轉換期間,表明全新世早期高緯冰蓋對亞洲季風變化的影響占主導地位,中晚全新世冰蓋基本消融,太陽輻射對亞洲季風的變化才開始起主導作用[5]。

氣候變化既要研究極地冰蓋對地球表面冷熱變化的影響,也要研究包括季風在內的熱帶過程,探索低緯的水文循環[43]。通過遙感觀測對季風的新認識:季風是熱帶輻合帶 (ITCZ)季節性遷移的表現[44]。亞洲季風具有歲差尺度(19~23 ka)的周期變化[10,41],那么季風的歲差周期也就表現為ITCZ的南北位移[45]。委內瑞拉Cariaco盆地(10°43′N,65° 10′W)是研究ITCZ移動很好的地點,Huag[46]在研究Cariaco盆地沉積物時發現10.5~5.4 ka B.P.時段Ti 和Fe的濃度最高,說明該時段降雨多。LHD5石筍δ18O值與其 Ti濃度變化一致性表明亞洲季風與ITCZ密切相關,即歲差使得全新世早期北半球太陽輻射增強,推動ITCZ北移,導致亞洲夏季風增強。Baker[47]對秘魯南部喀喀湖沉積物研究表明在4.0~2.4 ka B.P.降水多,巴西南部Botuverá洞石筍記錄表明全新世早期降水少,6 ka B.P.以后降水增加[48],LHD5石筍與巴西南部石筍δ18O值呈 “see-saw”關系,證明了在軌道尺度上南北半球低緯季風區共同受到熱帶輻合帶(ITCZ)的影響[7,48],這種關系也得到了亞洲季風區其他地質記錄的證實[5,7]。

全新世早期LHD5石筍與NGRIP冰芯不僅在總趨勢上一致,在谷峰值的波動上也一一對應,說明LHD5石筍記錄的弱季風事件在百年—十年尺度上受高緯冰量變化的影響。Teller et al.[49]發現全新世早期Agassiz湖發生10次淡水爆發事件,LHD5記錄的10 354±36 a B.P.,9 957±25 a B.P.,9 062±36 a B.P.弱季風事件與其中的3次淡水爆發事件對應,LHD5支持了全新世早期氣候波動與淡水爆發有關[50-51]。LHD5記錄的8 144±24 a B.P.弱季風事件非常接近8.2 ka B.P.事件,關于驅動8.2 ka B.P.事件的機制有學者認為受太陽活動影響[16-17],然而太陽活動并不能很好的解釋亞洲季風與南美季風反相位的關系,Cheng[39]認為8.2 ka B.P.事件由Agassiz和Ojibway湖淡水釋放造成,大量淡水通過哈德遜灣注入北大西洋導致北大西洋經向環流(AMOC)減弱,推動熱帶輻合帶的南移,減弱亞洲季風強度[48]。此外,由于8.2 ka B.P.事件亦處于弱的太陽活動期(圖7a),本文認為LHD5石筍8.2 ka B.P.事件受太陽活動和Agassiz 和Ojibway湖淡水釋放共同影響。然而,Rohling[52]指出由于北大西洋淡水注入要早于亞洲季風區8.2 ka B.P.事件,并不支持淡水注入是造成8.2 ka B.P.事件的原因。Liu[40]則指出亞洲季風區8.2 ka B.P.事件與北大西洋地區的遙相關關系,是由快速的大氣傳輸造成的。Magny[53]在研究全新世早期的氣候波動時將太陽活動和淡水爆發結合起來進行討論,為研究全新世早期的氣候振蕩提供了新思路。Solanki[33]指出太陽黑子反映太陽活動強度,太陽黑子數量多,太陽活動強,反之亦然。圖7中LHD5石筍在10 354± 36 a B.P.,9 062±36 a B.P.,8 144±24 a B.P.發生弱季風事件時都對應著太陽黑子的低值,即與弱的太陽活動相對應,說明全新早期LHD5石筍記錄的弱季風事件在一定程度上受到太陽活動的影響。

為進一步分析太陽活動對亞洲季風的影響,采用REDFIT38軟件對LHD5石筍δ18O值進行功率譜分析(圖8)。在95%置信度水平下通過紅噪聲檢驗的周期主要有88 a,80 a,59 a,31 a,27 a,23 a,20 a,19 a,17 a。88 a周期為太陽活動的Gleissberg周期,LHD5石筍檢驗出的80 a周期與董歌洞D4石筍的81 a周期非常接近[5]。59 a周期亦反映太陽活動周期[3],Agnihotri et al.[54]的研究表明,60 a的太陽周期直接控制了印度夏季風強度的變化。31 a周期接近樹輪記錄的14C活動周期,在先前石筍研究中也有發現[8,17]。23 a,20 a,19 a,17 a周期接近太陽太陽活動22 a的雙海爾周期[8,55],Cai et al.[56]在研究九仙洞石筍也發現明顯的23 a周期,Cai etal.對東石崖洞DSY09石筍[18]及Dykoski對董歌洞D4石筍[5]的研究中分別發現20 a和19 a周期。從以上分析可知LHD5石筍所反映的亞洲季風不僅在軌道尺度受太陽輻射控制,在百年—十年尺度也要受到太陽活動的影響。

圖8 LHD5石筍功率譜分析圖Fig.8 Power spectral analysis of LHD5 stalagmite

6 結論

通過分析湖南龍山縣蓮花洞LHD5石筍δ18O記錄得到以下結論:

(1)LHD5石筍記錄的全新世早期亞洲季風呈逐漸增強的趨勢,對比研究表明亞洲季風在全新世早期相位具有一致性,主要受北半球太陽輻射的增強推動熱帶輻合帶(ITCZ)北移影響。LHD5石筍記錄的YD結束時間為11 748±30 a B.P.,全新世開始于11 684±39 a B.P.,與亞洲季風區大多地質記錄略有差異,轉換時間約為64年,與格陵蘭gicc05記錄在誤差范圍內一致。全新世早期LHD5石筍記錄的弱季風事件與NGRIP冰芯記錄的降溫事件一一對應,整個演變趨勢也極為一致,表明亞洲季風與高緯氣候變化密切相關。

(2)全新世早期LHD5石筍記錄的6次弱季風事件,11 461±34 a B.P.(1.08‰)弱季風事件δ18O值波動幅度與8 144±24 a B.P.(1.02‰)弱季風事件相當,10 354±36 a B.P.、9 957±25 a B.P.、9 062±36 a B.P.弱季風事件時長約100年,且10 354±36 a B.P.弱季風事件的特點為快速進入,波動退出。8 744±23 a B.P.弱季風事件δ18O值波動幅度最小,且只有少數地質記錄識別到該弱季風事件。通過對比發現這些弱季風事件在測年誤差范圍內具有區域意義,除8.2 ka事件之外,10 ka B.P.之前的弱季風事件除了受到太陽活動的影響,還受到北大西洋IRD事件的影響,而之后更多地受到太陽活動和ITCZ南移的影響。

參考文獻(References)

1 Overpeck J,Anderson D,Trumbore S,et al.The southwest Indian Monsoon over the last 18 000 years[J].Climate Dynamics,1996,12 (3):213-225.

2 Gupta A K,Anderson D M,Overpeck JT.Abrupt changes in the Asian southwestmonsoon during the Holocene and their links to the North Atlantic Ocean[J].Nature,2003,421(6921):354-357.

3 Fleitmann D,Burns S J,Mudelsee M,et al.Holocene forcing of the Indian monsoon recorded in a stalagmite from southern Oman[J].Science,2003,300(5626):1737-1739.

4 Hong Y T,Hong B,Lin Q H,etal.Correlation between Indian Ocean summermonsoon and North Atlantic climate during the Holocene[J]. Earth and Planetary Science Letters,2003,211(3/4):371-380.

5 Dykoski C A,Edwards R L,Cheng H,et al.A high-resolution,absolute-dated Holocene and deglacial Asian monsoon record from Dongge Cave,China[J].Earth and Planetary Science Letters,2005,233(1/ 2):71-86.

6 Staubwasser M.An overview of Holocene South Asianmonsoon records-Monsoon domains and regional contrasts[J].Journal of the Geological Society of India,2006,68(3):433-446.

7 Dong JG,Wang Y J,Cheng H,etal.A high-resolution stalagmite record of the Holocene East Asianmonsoon from Mt Shennongjia,central China[J].The Holocene,2010,20(2):257-264.

8 楊琰,袁道先,程海,等.末次冰消期亞洲季風突變事件的精確定年:以貴州衙門洞石筍為例[J].中國科學(D輯):地球科學,2010,40(2):199-210.[Yang Yan,Yuan Daoxian,Cheng Hai,etal. Precise dating of abruptshifts in the Asian Monsoon during the lastdeglaciation based on stalagmite data from Yamen Cave,Guizhou province,China[J].Science China(Seri.D):Earth Sciences,2010,40 (2):199-210.]

9 吳江瀅,汪永進,董進國.全新世東亞夏季風演化的遼寧暖和洞石筍δ18O記錄[J].第四紀研究,2011,31(6):990-998.[Wu Jiangying,Wang Yongjin,Dong Jinguo.Change in East Asian summermonsoon during the Holocene recorded by stalagmiteδ18O records from Liaoning province[J].Quaternary Sciences,2011,31(6):990-998.]

10 Yuan D X,Cheng H,Edwards R L,et al.Timing,duration,and transitions of the last interglacial Asianmonsoon[J].Science,2004,304(5670):575-578.

11 Fleitmann D,Burns SJ,ManginiA,etal.Holocene ITCZ and Indian monsoon dynamics recorded in stalagmites from Oman and Yemen (Socotra)[J].Quaternary Science Reviews,2007,26(1/2):170-188.

12 Zhang H L,Yu K F,Zhao JX,et al.East Asian Summer Monsoon variations in the past12.5 ka:High-resolutionδ18O record from a precisely dated aragonite stalagmite in central China[J].Journal of Asian Earth Sciences,2013,73:162-175.

13 覃嘉銘,袁道先,程海,等.新仙女木及全新世早中期氣候突變事件:貴州茂蘭石筍氧同位素記錄[J].中國科學(D輯):地球科學,2004,34(1):69-74.[Qin Jiaming,Yuan Daoxian,Cheng Hai,etal.The Y.D.and climate abrupt events in the early and Middle Holocene:Stalagmite oxygen isotope record from Maolan,Guizhou,China[J].Science China(Seri.D):Earth Sciences,2004,34(1):69-74.]

14 邵曉華,汪永進,程海,等.全新世季風氣候演化與干旱事件的湖北神農架石筍記錄[J].科學通報,2006,51(1):80-86.[Shao Xiaohua,Wang Yongjin,Cheng Hai,et al.Long-term trend and abrupt events of the Holocene Asian monsoon inferred from a stalagmiteδ10O record from Shennongjia in central China[J].Chinese Science Bulletin,2006,51(1):80-86.]

15 Teller,JT,Leverington,DW.Glacial Lake Agassiz:a5000 yr history of change and its relationship to the record of Greenland[J].Geological Society of America Bulletin,2004,116(5):729-742.

16 Neff U,Burns S J,Mangini A,et al.Strong coherence between solar variability and themonsoon in Oman between 9 and 6 kyr ago[J].Nature,2001,411(6835):290-293.

17 Wang Y J,Cheng H,Edwards R L,et al.The Holocene Asianmonsoon:links to solar changes and North Atlantic climate[J].Science,2005,308(5723):854-857.

18 Cai B G,Edwards R L,Cheng H,et al.A dry episode during the Younger Dryas and centennial-scale weak monsoon events during the early Holocene:A high-resolution stalagmite record from southeast of the Loess Plateau,China[J].Geophysical Research Letters,2008,35(2):L02075,doi:10.1029/2007GL030986.

19 殷建軍.湖南湘西龍山蓮花洞近2000年來氣候變化石筍記錄研究[D].重慶:西南大學,2013.[Yin Jianjun.Climate change in the past 2000 years revealed by stalagmites from Lianhua Cave in Northwest Hunan,China[D].Chongqing:Southwest University,2013.]

20 Jaffey A H,Flynn K F,Glendenin L E,etal.Precision measurement of half-lives and specific activities of235U and238U[J].Physical Review C,1971,4(5):1889-1906.

21 Cheng H,Edwards R L,Shen CC,et al.Improvements in230Th dating,230Th and234U half-life values,and U-Th isotopic measurements bymulti-collector inductively coupled plasmamass spectrometry[J]. Earth and Planetary Science Letters,2013,371-372:82-91.

22 Hendy C H.The isotopic geochemistry of speleothems—Ι.The calculation of the effectsof differentmodes of formation on the isotopic composition of speleothems and their applicability as palaeoclimatic indicators[J].Geochimica et Cosmochimica Acta,1971,35(8):801-824.

23 Dorale JA,Liu Z.Limitations of Hendy test criteria in judging the paleoclimatic suitability of speleothems and the need for replication[J]. Journal of Cave and Karst Studies,2009,71(1):73-80.

24 Wang Y J,Cheng H,Edwards R L,etal.A high-resolution absolutedated late Pleistocene monsoon record from Hulu Cave,China[J]. Science,2001,294(5550):2345-2348.

25 Wang Y J,Cheng H,Edwards R L,etal.Millennial-and orbital-scale changes in the East Asian monsoon over the past 224,000 years[J]. Nature,2008,451(7182):1090-1093.

26 程海,艾思本,王先鋒,等.中國南方石筍氧同位素記錄的重要意義[J].第四紀研究,2005,25(2):157-163.[Cheng Hai,Edwards R L,Wang Xianfeng,etal.Oxygen isotope recordsofstalagmites from southern China[J].Quaternary Sciences,2005,25(2):157-163.]

27 Cosford J,Qing H R,Eglington B,et al.East Asian monsoon variability since the Mid-Holocene recorded in a high-resolution,absolutedated aragonite speleothem from eastern China[J].Earth and Planetary Science Letters,2008,275(3/4):296-307.

28 Hu C Y,Henderson G M,Huang JH,et al.Quantification of Holocene Asian monsoon rainfall from spatially separated cave records [J].Earth and Planetary Science Letters,2008,266(3/4):221-232.

29 Vinther B M,Clausen H B,Johnsen SJ,etal.A synchronized dating of three Greenland ice cores throughout the Holocene[J].Journal of Geophysical Research,2006,111(D13):D13102,doi:10.1029/ 2005JD006921.

30 Rasmussen SO,Andersen K K,Svensson A M,et al.A new Greenland ice core chronology for the last glacial termination[J].Journalof Geophysical Research,2006,111(D6):D06102,doi:10.1029/ 2005JD006079.

31 An Z S,Porter SC,Kutzbach JE,etal.Asynchronous Holocene optimum of the East Asian monsoon[J].Quaternary Science Reviews,2000,19(8):743-762.

32 Bond G,Kromer B,Beer J,et al.Persistent solar influence on North Atlantic climate during the Holocene[J].Science,2001,294 (5549):2130-2136.

33 Solanki SK,Usoskin IG,Kromer B,et al.Unusual activity of the Sun during recent decades compared to the previous 11,000 years [J].Nature,2004,431(7012):1084-1087.

34 張銀環,楊琰,楊勛林,等.早全新世季風演化的高分辨率石筍δ18O記錄研究——以河南老母洞石筍為例[J].沉積學報,2015,33(1):134-141.[Zhang Yinhuan,Yang Yan,Yang Xunlin,et al. Early Holocenemonsoon evolution of high-resolution stalagmiteδ18O records:in Henan Laomu Cave[J].Acta Sedimentologica Sinica,2015,33(1):134-141.]

35 Rasmussen SO,Vinther B M,Clausen H B,et al.Early Holocene climate oscillations recorded in three Greenland ice cores[J].Quaternary Science Reviews,2007,26(15/16):1907-1914.

36 Bond G,ShowersW,Cheseby M,et al.A pervasivemillennial-scale cycle in North Atlantic Holocene and glacial climates[J].Science,1997,278(5341):1257-1266.

37 王寧練,姚檀棟,Thompson L G,等.全新世早期強降溫事件的古里雅冰芯記錄證據[J].科學通報,2002,47(11):818-823.[Wang Ninglian,Yao Tandong,Thompson L G,et al.Evidence for cold events in the Early Holocene from the Guliya Ice Core,Tibetan Plateau,China[J].Chinese Science Bulletin,2002,47(11):818-823.]

38 劉嘉麒,呂厚遠,Negendank J,等.湖光巖瑪珥湖全新世氣候波動的周期性[J].科學通報,2000,45(11):1190-1195.[Liu Jiaqi,Lü Houyuan,Negendank J,et al.Periodicity of Holocene climatic variations in the Huguangyan Maar Lake[J].Chinese Science Bulletin,2000,45(11):1190-1195.]

39 Cheng H,Fleitmann D,Edwards R L,et al.Timing and structure of the 8.2kyr B.P.event inferred fromδ18O records of stalagmites from China,Oman,and Brazil[J].Geology,2009,37(11):1007-1010.

40 Liu Y H,Henderson GM,Hu C Y,etal.Links between the East A-sian monsoon and North Atlantic climate during the 8,200 year event [J].Nature Geoscience,2013,6(2):117-120.

41 Cheng H,Edwards R L,Broecker W S,et al.Ice age terminations [J].Science,2009,326(5950):248-252.

42 Yang X L,Liu JB,Liang F Y,et al.Holocene stalagmiteδ18O records in the East Asianmonsoon region and their correlation with those in the Indian monsoon region[J].The Holocene,2014,24(12):1657-1664,doi:10.1177/0959683614551222.

43 汪品先.低緯過程的軌道驅動[J].第四紀研究,2006,26(5):694-701.[Wang Pinxian.Orbital forcing of the low-latitude processes[J]. Quaternary Sciences,2006,26(5):694-701.]

44 Gadgil S.The Indian monsoon and its variability[J].Annual Review of Earth and Planetary Sciences,2003,31(1):429-467.

45 汪品先.全球季風的地質演變[J].科學通報,2009,54(5):535-556.[Wang Pinxian.Globalmonsoon in a geological perspective[J]. Chinese Science Bulletin,2009,54(5):535-556.]

46 Haug G H,Hughen K A,Sigman D M,et al.Southward migration of the intertropical convergence zone through the Holocene[J].Science,2001,293(5533):1304-1308.

47 Baker P A,Seltzer G O,Fritz SC,etal.The history of South American tropical precipitation for the past 25,000 years[J].Science,2001,291(5504):640-643.

48 Wang X F,Auler A S,Edwards R L,etal.Millennial-scale precipitation changes in southern Brazil over the past90,000 years[J].Geophysical Research Letters,2007,34(23):L23701,doi:10.1029/ 2007GL031149.

49 Teller JT,Leverington DW,Mann JD.Freshwater outbursts to the oceans from glacial Lake Agassiz and their role in climate change during the last deglaciation[J].Quaternary Science Reviews,2002,21 (8/9):879-887.

50 Mayewski P A,Rohling E E,Stager JC,etal.Holocene climate variability[J].Quaternary Research,2004,62(3):243-255.

51 Hong Y T,Hong B,Lin Q H,et al.Synchronous climate anomalies in the western North Pacific and North Atlantic regions during the last 14,000 years[J].Quaternary Science Reviews,2009,28(9/10):840-849.

52 Rohling E J,P?like H.Centennial-scale climate coolingwith a sudden cold event around 8,200 years ago[J].Nature,2005,434(7036):975-979.

53 Magny M,Vannière B,De Beaulieu J L,et al.Early-Holocene climatic oscillations recorded by lake-level fluctuations in west-central Europe and in central Italy[J].Quaternary Science Reviews,2007,26(15/16):1951-1964.

54 Mavromichalaki H,BelehakiA,Rafios X,et al.Hale-cycle effects in cosmic-ray intensity during the last four cycles[J].Astrophysics and Space Science,1996,246(1):7-14.

55 AgnihotriR,Dutta K,Bhushan R,etal.Evidence for solar forcing on the Indianmonsoon during the lastmillennium[J].Earth and Planetary Science Letters,2002,198(3/4):521-527.

56 Cai Y J,Tan L C,Cheng H,et al.The variation of summermonsoon precipitation in central China since the lastdeglaciation[J].Earth and Planetary Science Letters,2010,291(1/2/3/4):21-31.

Discussion about the M echanism of the W eak Summer M onsoon Events during the Early Holocence:A case study of precisely dated stalagm ite record from Lianhua Cave,Hunan province,China

ZHANG HuaSheng1,2YIN JianJun2CHENG Hai3,4R Lawrence Edwards4LIN YuShi2TANGWei2YANG Hui2TU LinLing2WANG Hua2PAN MouCheng2WU Xia2
(1.School of Geographical Sciences,Southwest University,Chongqing 400715,China;2.Key Laboratory of Karst Dynam ics,M LR&Guangxi;International Research Center on Karst;Institute of Karst Geology,CAGS,Guilin 541004,China;3.Institute of G lobal Environmental Change,Xi'an Jiaotong University,Xi'an 710054,China;4.Department of Earth Science,University of M innesota,M inneapolis 55455,USA)

Abstract:The Early Holocene is an important period of Solar radiation strengthening,globalwarming and along with ice sheetsmelting.The cause of the colder events and the weaker summermonsoon in the Asian monsoon region during the Early Holocene has always been the focus of research.Ithas important significance for the study of the linksamong the Asianmonsoon,ocean and polar regions.Based on 28 U/Th dates and 535 oxygen isotopic data of stalagmite LHD5 from Lianhua Cave,Hunan province,China,a Holocene Asian monsoon evolution record was reconstructed,and the average resolution is8 year in the Early Holocene.From the LHD5 stalagmite record,the end of Younger Dryas was at 11 748±30 a B.P.,the startof the Holocene was at11 684±39 a B.P.,and the conversion time is about 64 years,which is consistentwith the recordsof the gicc05 ice core from Greenland in the error range.The centralage of the six weak summermonsoon eventswere at11 461±34 a B.P.(1.08‰),10 354±36 a B.P.(0.94‰),9 957±25 a B.P.(0.66‰),9 062±36 a B.P.(0.90‰),8 744±23 a B.P.(0.55‰)and 8 144±24 a B.P.(1.02‰)in the record of the stalagmite LHD5 during the Early Holocene,respectively.These weaker summermonsoon events in the A-sianmonsoon region have a common regional significance.In addition to the 8.2 ka event,the weaker summermonsoon events before 10 ka B.P.were notonly affected by the summer insolation but also by the ice-rafted debris(IRD)events in the North Atlantic,but the following weak summermonsoon events were mostly influenced by solar activity and the southward migration of themean position of the Inter-Tropical Convergence Zone.

Key words:Early Holocene;weak summermonsoon events;stalagmite;δ18O;Lianhua cave

第一作者簡介張華生 男 1988年出生 碩士 石筍古氣候研究 E-mail:peanutzhs@163.com

通訊作者殷建軍 男 助理研究員 E-mail:david1985_2005@163.com

中圖分類號P532

文獻標識碼A

文章編號:1000-0550(2016)02-0281-11

doi:10.14027/j.cnki.cjxb.2016.02.007

收稿日期:2015-07-06;收修改稿日期:2015-08-03

基金項目:中國地質科學院基本科研業務費項目(YWF201414);巖溶地質研究所基本科研業務費項目(2014025,201319);中國地質調查局地調項目(12120113006700)[Foundation:Fundamental Research Fund of Chinese Academy of Geological Science,No.YWF201414;Fundamental Research Fund of Institute of Karst Geology,China Geological Survey,No.2014025,201319;China Geological Survey Project,No. 12120113006700]

主站蜘蛛池模板: 国产一在线| 91亚洲免费视频| 91精品国产情侣高潮露脸| 亚洲天堂啪啪| 久久综合色天堂av| 日韩黄色精品| 欧美高清三区| 无码中文字幕精品推荐| 福利一区三区| 波多野结衣一级毛片| 国产精品亚洲日韩AⅤ在线观看| 免费一级毛片在线播放傲雪网| 国产一级无码不卡视频| 拍国产真实乱人偷精品| 人妻熟妇日韩AV在线播放| 日韩国产一区二区三区无码| 极品性荡少妇一区二区色欲| 在线观看视频99| 波多野结衣爽到高潮漏水大喷| 大乳丰满人妻中文字幕日本| 亚洲自偷自拍另类小说| 片在线无码观看| 国产剧情无码视频在线观看| 91免费国产在线观看尤物| 亚洲欧洲日韩国产综合在线二区| 91在线日韩在线播放| 激情综合网激情综合| 激情爆乳一区二区| 欧美日韩免费| 国产精品九九视频| 亚洲日本韩在线观看| 中文字幕欧美成人免费| 在线观看av永久| 国产视频欧美| 91精品啪在线观看国产91九色| 国禁国产you女视频网站| 国产成人亚洲综合a∨婷婷| 国产性生大片免费观看性欧美| 中文字幕不卡免费高清视频| 毛片大全免费观看| 久久99国产乱子伦精品免| 无码网站免费观看| 欧美视频免费一区二区三区| 免费不卡在线观看av| 国产日韩精品欧美一区喷| 亚洲成人一区在线| 精品国产欧美精品v| 大乳丰满人妻中文字幕日本| 亚洲福利视频一区二区| 1024你懂的国产精品| 91啪在线| 久久精品免费国产大片| 亚洲最大综合网| 999国内精品视频免费| 日本一区二区三区精品视频| 91黄视频在线观看| 色婷婷电影网| 四虎免费视频网站| 国产毛片不卡| 亚洲AV无码不卡无码| 亚洲无码在线午夜电影| 国产精品不卡片视频免费观看| 国产视频你懂得| a亚洲视频| 2020国产在线视精品在| 三级国产在线观看| 丝袜国产一区| 免费一级毛片在线观看| 国产在线欧美| 乱人伦99久久| 99无码熟妇丰满人妻啪啪| 成人免费黄色小视频| 中文字幕在线视频免费| a级毛片在线免费观看| 色综合五月婷婷| 国产午夜一级毛片| 国产精品香蕉在线观看不卡| 亚洲性视频网站| 亚洲精品色AV无码看| 国产高潮流白浆视频| 国产电话自拍伊人| 亚洲日韩第九十九页|