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

    河南馬溝洞石筍記錄的早中全新世氣候和環(huán)境變化

    2016-03-21 03:17:01毛瑞雪蔡演軍
    地球環(huán)境學(xué)報 2016年3期
    關(guān)鍵詞:石筍季風(fēng)同位素

    毛瑞雪,蔡演軍,馬 樂,成 星

    (1.中國科學(xué)院地球環(huán)境研究所 黃土與第四紀(jì)地質(zhì)國家重點實驗室,西安 710061;2.中國科學(xué)院大學(xué),北京 100049;3.西安交通大學(xué) 全球環(huán)境變化研究院,西安 710054)

    河南馬溝洞石筍記錄的早中全新世氣候和環(huán)境變化

    毛瑞雪1,2,蔡演軍1,3,馬 樂1,成 星1,2

    (1.中國科學(xué)院地球環(huán)境研究所 黃土與第四紀(jì)地質(zhì)國家重點實驗室,西安 710061;2.中國科學(xué)院大學(xué),北京 100049;3.西安交通大學(xué) 全球環(huán)境變化研究院,西安 710054)

    基于河南省禹州市馬溝洞2支石筍(MG-1與MG-40)24個230Th年齡和1988個氧同位素數(shù)據(jù),建立了研究區(qū)13.1 — 4.9 ka BP分辨率為2 — 14 a的石筍氧同位素時間序列。馬溝洞石筍δ18O的時間序列揭示季風(fēng)降水在11.2 — 9.1 ka BP時段在波動中逐漸增加,9.1 — 4.9 ka BP季風(fēng)降水顯著波動但無明顯長期趨勢變化。YD事件、9.3 ka事件、8.2 ka事件記錄與其他石筍δ18O記錄的一致性揭示末次冰消期—早全新世百年—千年尺度氣候突變事件的大范圍存在和共同的驅(qū)動因子。馬溝洞全新世大暖期的δ18O記錄中檢測出的13個數(shù)十年尺度的弱季風(fēng)事件進(jìn)一步證實全新世大暖期氣候的不穩(wěn)定性。與大氣Δ14C記錄、NGRIP冰芯δ18O記錄的對比及周期分析揭示,太陽活動引起的太陽輻射變化和北半球高緯氣候狀況共同影響著亞洲季風(fēng)的變化,ENSO活動及氣候系統(tǒng)內(nèi)部的相互作用也對東亞夏季風(fēng)降水產(chǎn)生重要的影響。

    石筍;δ18O;氣候變化;夏季風(fēng)降水;全新世;馬溝洞

    河南省地處黃河中下游地區(qū),在地形、氣候、水系、植被等地理要素方面具有東西過渡、南北兼具的特征(李中軒等,2013),特殊的地理位置一方面導(dǎo)致這一地區(qū)對氣候敏感響應(yīng),另一方面也成為中原文明的核心區(qū)域。該區(qū)域過去氣候和環(huán)境變化的研究既有助于深刻認(rèn)識東亞季風(fēng)氣候?qū)ξ覈鴸|部平原地區(qū)的影響,也可為查明中原文明演替與氣候變化的關(guān)系提供重要的環(huán)境變化信息。洞穴石筍由于具有可精確定年、高分辨率采樣、沉積時間跨度大、后期沒有擾動等諸多優(yōu)勢,已成為古氣候研究的重要載體之一(馬樂等,2015)。目前有關(guān)河南地區(qū)石筍研究較少,已報道的石筍研究主要集中在豫西山地丘陵區(qū),例如欒川縣的東石崖洞、老母洞、雞冠洞,南陽伏牛山地區(qū)的南陽洞穴。這些研究中有關(guān)洞穴監(jiān)測、洞穴滴水及現(xiàn)代沉積物微量元素記 錄研究較多(馬睿,2013;任小鳳等,2014;趙景耀等,2014),石筍氧、碳同位素研究較少,且研究時段主要集中在早全新世。任曉輝等(2006)通過對南陽洞石筍碳、氧同位素記錄進(jìn)行對比研究,確認(rèn)了研究區(qū)6115 — 4544 a BP期間發(fā)生了3次弱季風(fēng)事件,其中5350 — 5080 a BP、4770 — 4640 a BP弱季風(fēng)事件可與格陵蘭冰芯記錄的降溫事件相對應(yīng),認(rèn)為亞洲季風(fēng)在百年尺度上與高緯極地氣候變化存在遙相關(guān)。Cai et al(2008)研究發(fā)現(xiàn)東石崖洞DSY09記錄的早全新世11.5 — 9.0 ka BP 8次百年 — 十年尺度弱季風(fēng)事件與格陵蘭冷事件在測年誤差內(nèi)一致,進(jìn)而指出東亞季風(fēng)強(qiáng)度與太陽活動相關(guān)。張銀環(huán)等(2015)通過對老母洞研究發(fā)現(xiàn)在測年誤差范圍內(nèi)LM2石筍δ18O與格陵蘭GRIP冰芯δ18O記錄有著較好的一致變化,并記錄了8.2 ka、9.5 ka、10.2 ka、10.9 ka弱季風(fēng)事件,認(rèn)為早全新世亞洲季風(fēng)的演化與太陽活動變化引起的太陽輻射能量的變化和北半球高緯氣候的變化狀況有關(guān)。

    本文基于河南禹州馬溝洞石筍230Th定年結(jié)果和δ18O測試數(shù)據(jù),建立了該地區(qū)13.1 — 4.9 ka BP(相對于1950 a AD)時段石筍δ18O的高分辨率時間序列,結(jié)合已有石筍δ18O記錄,探討早中全新世區(qū)域降水變化的特點和機(jī)理,并在此基礎(chǔ)上討論了文明演替與氣候變化的可能聯(lián)系。

    1 研究區(qū)域和石筍樣品

    馬溝洞(113°23′E,34°19′N,海拔422 m)位于河南省許昌市下轄的禹州市馬溝村(圖1),距離鄭州市區(qū)西南約50 km。馬溝洞地處豫西山地丘陵區(qū)向黃淮海平原區(qū)的過渡地帶,發(fā)育于禹州市北部具茨山系奧陶紀(jì)灰?guī)r地層中,上覆巖層厚約20 m,洞外植被以灌草為主。區(qū)域氣候?qū)儆谂瘻貛Ъ撅L(fēng)氣候,主要受高緯冷高壓南下的冷干氣流(冬季風(fēng))和來自印度-太平洋的暖濕氣流(夏季風(fēng))交替影響。冬季寒冷少雨雪,春季干旱多風(fēng)沙,夏季炎熱而降雨集中,秋季晴和氣爽且日照充足。年平均降雨量為705.6 mm左右,各季節(jié)降水量分布懸殊,6 — 9月雨水最為集中,平均為439.9 mm,占年平均降水量的62.3%(數(shù)據(jù)來自許昌1971 — 2000年資料統(tǒng)計)。歷年年均氣溫在13.0 — 16.0℃,年平均無霜期為218天。該洞穴被采石場采石發(fā)現(xiàn)之前是一個封閉的洞穴,沒有開放的洞口。雖然在洞口炸開后沒有進(jìn)行洞穴環(huán)境監(jiān)測,但可以認(rèn)為該洞穴的溫度變化基本響應(yīng)于地表年均溫的變化,洞內(nèi)濕度接近飽和狀態(tài),應(yīng)十分適合古氣候重建研究。

    所研究石筍樣品編號分別為MG-1、MG-40(圖2)。石筍MG-1總長度約260 mm。該石筍沿生長軸剖開并進(jìn)行剖光后可見清晰的生長紋層,其中距頂部148 — 260 mm段生長軸多次發(fā)生改變,195 — 260 mm紋層致密,生長速度緩慢。MG-1采樣時斷開,分別標(biāo)號為MG-1A、MG-1B。MG-1A頂部缺失,顏色為黃色,紋層清晰,整體上紋層較稀疏,生長速度快。MG-1B整體上淡黃色、棕色相間,66 — 100 mm段和148 — 260 mm段為淡黃色。石筍MG-40總長度約40 mm,整體呈暗黃—棕色,樣品下部18 mm出現(xiàn)多條明顯的沉積間斷,最下部8 mm為基座部分,土黃—白條紋相間,表明碎屑質(zhì)的顯著增多。

    圖1 馬溝洞地理位置示意圖Fig.1 Geographical location of Magou Cave

    圖2 馬溝洞石筍(MG-1、MG-40)剖面圖及年代、測年樣點位置以及年齡-深度模式圖(黑點表示測年點位置;黑色橫線為誤差)Fig.2 The fi gures of the polished sections with position of dating samples and the age-depth model of the two stalagmites (MG-1, MG-40) from Magou Cave

    2 實驗分析方法與結(jié)果

    在石筍拋光面上確定測年點位置后,沿著生長紋層方向用0.5 mm牙鉆共取得24個測年樣品,單個樣品質(zhì)量約為70 — 120 mg。年代樣品在西安交通大學(xué)全球環(huán)境變化研究院U-系年代學(xué)實驗室測試完成,測試儀器為MC-ICP-MS Neptune,鈾釷化學(xué)分離使用Edwards et al(1987)建立的化學(xué)流程,儀器測試使用Cheng et al(2013)建立的測試方法,年齡誤差(2σ)≤1%,測試結(jié)果見表1。年代標(biāo)尺采用線性插值方法,對于測年點以外部分的年代,使用鄰近部分的沉積速率進(jìn)行外推計算。

    表1 馬溝洞石筍230Th測年結(jié)果Tab.1230Th dating results of stalagmites in Magou Cave

    如表1、圖2所示,除MG1B-3S和MG40-1232Th相對較高導(dǎo)致年代誤差較大外(分別為±168、±784),其余樣品點測年誤差均為±28 — ±80。盡管MG1B-3S和MG40-1兩個樣品誤差較大,但所測樣品年齡均按沉積先后序列正常排列,沒有出現(xiàn)年齡倒轉(zhuǎn)。同時對疑似沉積間斷的兩側(cè)做了年代測定,認(rèn)為不存在沉積間斷,或沉積間斷的時間很短,并不影響本文相關(guān)的氣候變化討論。本文通過對相鄰實測年代點進(jìn)行線性等間距內(nèi)插與外延,建立了石筍剖面13.1— 4.9 ka BP年代框架。MG-40生長較慢,13.1— 9.8 ka BP和9.8— 7.1 ka BP生長速率分別為1.5 μm · a-1、6.3 μm · a-1;MG-1總體生長較快,10.3— 8.6 ka BP和8.6— 4.9 ka BP生長速率分別為4.7 μm · a-1、58.6 μm · a-1。圖2為兩個石筍的年代模式圖,顯示兩個石筍在開始沉積時速率較慢,隨后增大并基本保持不變。

    本文石筍碳、氧同位素樣品使用Micro-mill微取樣裝置采取。石筍MG-40的采樣間距為50 μm,MG-1A在0 — 3.6 mm的采樣間距為50 μm,其他樣品的采樣間距均為150 μm。MG-40、MG-1A、MG-1B分別采取439個、381個和1168個子樣品進(jìn)行碳氧同位素組成分析。碳、氧同位素分析測試全部在中國科學(xué)院地球環(huán)境研究所洞穴實驗室完成。分析方法采用磷酸法(McCrea,1950),分析儀器為IsoPrime100穩(wěn)定同位素質(zhì)譜儀,在分析測試時每10個同位素樣品內(nèi)插一個標(biāo)準(zhǔn)NBS-19,使用已知同位素組成的高純CO2氣體為參考?xì)?,?(Rsample/Rstd-1)×1000,其中R=18O /16O,δ值相對于PDB標(biāo)準(zhǔn)給出,分析誤差(2σ)≤± 0.08‰。

    圖3為馬溝洞石筍MG-1和MG-40δ18O變化的時間序列。石筍MG-1涵蓋的5479年內(nèi)共測試碳氧同位素1168對,平均分辨率約4.7年。整個δ18O序列沒有明顯的變化趨勢,δ18O值的變化范圍為-12.7‰ — -9.2‰,平均值為-10.5‰,波動幅度達(dá)3.5‰。石筍MG-40涵蓋的6075年內(nèi)共測試碳氧同位素439對,平均分辨率約14年。δ18O值的變化范圍為-11.9‰ — -8.4‰,平均值為-10.4‰,波動幅度達(dá)3.5‰。與δ18O的變化相比,兩個石筍δ13C的高頻大幅波動變化并不明顯,僅表現(xiàn)出長期變化趨勢。

    圖3 馬溝洞石筍(MG-1、MG-40)氧、碳同位素序列(上圖)及相同時段(10.37 — 7.06 ka BP)兩支石筍氧、碳同位素序列重復(fù)性對比圖(下圖)MG-1氧、碳同位素序列分別用紅色和綠色線條表示,MG-40氧、碳同位素序列分別用藍(lán)色和紫色線條表示。Fig.3 Theδ18O andδ13C time series of the two stalagmites (MG-1, MG-40) of Magou Cave (top panel), and comparisons of theδ18O andδ13C profi les during the overlapped time period (10.37 — 7.06 ka BP) (lower panel)Theδ18O andδ13C time series of the stalagmite MG-1 are shown by the red and green lines, respectively. Theδ18O andδ13C time series of the stalagmite MG-40 are shown by the blue and purple lines, respectively.

    3 討論

    3.1 石筍碳酸鹽δ18O的氣候意義

    應(yīng)用石筍氧同位素組成重建古氣候和古環(huán)境變化的前提條件之一就是石筍碳酸鹽在平衡條件下沉積。重復(fù)性檢驗是目前洞穴石筍研究的主要檢驗準(zhǔn)則。重復(fù)性檢驗的內(nèi)容就是不同石筍樣本相同時段的δ18O如果能夠重復(fù),那么這些石筍樣本沉積時受到動力分餾作用的影響就很小或可以忽略不計,石筍碳酸鹽δ18O就可以用作古環(huán)境重建(Dorale et al,1998;Wang et al,2001),反之,石筍碳酸鹽δ18O的古氣候和古環(huán)境解譯就需慎重。本論文研究的MG-1和MG-40在相同時段內(nèi)(10.37— 7.06 ka BP),盡管因測年誤差和樣點分辨率不同,δ18O記錄在部分細(xì)節(jié)上存在差異,如8.25— 8.0 ka BP、9.2—9.0 ka BP階段,δ18O記錄的峰谷存在些許偏差,但總體上MG-1和MG-40在這一時段能夠較好地對應(yīng)(圖3)。除此之外,石筍MG-40的δ18O也和東石崖洞石筍δ18O記錄在重疊時段具有一致的變化趨勢和特征。這些證據(jù)說明馬溝洞石筍是在同位素平衡分餾狀態(tài)下沉積的,其δ18O變化能夠真實地記錄氣候和環(huán)境變化信息。

    在平衡條件下沉積的石筍碳酸鹽,其δ18O值受滴水的δ18O和洞穴溫度的共同影響。馬溝洞石筍碳酸鹽δ18O由冰消期至早全新世的變化幅度為3.5‰,且早中全新世的波動幅度最大也可達(dá)2.5‰。僅用溫度的變化很難解釋這么大振幅的δ18O的變化,因為洞穴滴水與碳酸鹽之間的氧同位素組成分餾的溫度效應(yīng)很小,只有-0.23‰ ·℃-1(Kim and O'Neil,1997),而早中全新世中國東部的溫度變化也不會超過3℃(溫度變化造成的碳酸鹽δ18O的波動不會超過0.75‰)。因此,石筍碳酸鹽δ18O的變化主要反映滴水的δ18O變化,也就是降水同位素組成的變化。鄰近的鄭州氣象站(1985— 1992年)器測降水資料(圖4)表明,夏季風(fēng)控制期間(6— 9月)δ18O值偏輕,平均值約為-7.4‰,且夏季降水量約占全年降水量的62.3%,夏季風(fēng)強(qiáng)弱變化引起的降水量變化就成為影響區(qū)域降水δ18O組分的主要因素,因為夏季季風(fēng)降水量的變化或者其引起的夏季降水和冬季降水的比值的變化,均能導(dǎo)致年均大氣降水δ18O的顯著變化。盡管溫度是影響冬春季節(jié)大氣降水δ18O的主控因子(Johnson and Ingram,2004),但由于洞穴滴水-碳酸鹽的溫控分餾大體上抵消大氣溫度對降水δ18O的影響(Johnson et al,2006),因此溫度變化的影響較小,這里不做討論。綜合上述的分析,可以認(rèn)為馬溝洞石筍碳酸鹽δ18O變化主要指示了季風(fēng)降水的變化,δ18O值減小,季風(fēng)降水增加,δ18O值增加,季風(fēng)降水減少。

    圖4 鄭州市1985 — 1992年月均降水量、月均溫與大汽降水δ18O值對比Fig.4 Monthly average precipitation, temperature and precipitationδ18O during 1985 — 1992 in Zhengzhou

    3.2 石筍δ18O記錄的早中全新世區(qū)域氣候和環(huán)境變化

    馬溝洞石筍δ18O序列有以下兩個特點:(1)MG-1δ18O序列在百年至數(shù)十年尺度上變化的振幅大于葫蘆洞和董哥洞。這也與東石崖洞δ18O記錄相類似,可能表明馬溝洞臨近干濕區(qū)分界線(南部濕潤區(qū)和北部半干旱區(qū)),且處于豫西山地丘陵區(qū)向黃淮海平原區(qū)的過渡地帶,對夏季風(fēng)強(qiáng)度變化極其敏感。(2)馬溝洞δ18O值由9 ka BP直至5 ka BP一直未出現(xiàn)明顯偏正的變化趨勢,而董哥洞、和尚洞δ18O記錄顯示在8 — 6 ka BP,石筍碳酸鹽δ18O已經(jīng)開始增加,亞洲季風(fēng)逐漸減弱、夏季風(fēng)降水減少,進(jìn)一步支持了較早發(fā)現(xiàn)的全新世中國中東部地區(qū)東亞夏季風(fēng)降水變化具有不同步性,南部季風(fēng)降水量開始減少的時間較早,北部則較晚(Cai et al,2010),但由于缺失晚全新世數(shù)據(jù),本文對此不進(jìn)行深入討論。

    馬溝洞石筍δ18O記錄,提供了河南地區(qū)13.1 — 4.9 ka BP的氣候變化信息。依據(jù)δ18O變化特征,這一記錄大致可以劃分為13.1 — 11.2 ka BP、11.2 — 9.1 ka BP和9.1 — 4.9 ka BP三個時段。其中13.1 — 11.2 ka BP的年代是根據(jù)沉積速率線性外延得到,其δ18O變化特征及與董哥洞、葫蘆洞的對比(圖5)表明該階段可能對應(yīng)于Aller?d暖期和YD事件,但由于分辨率較低和年代標(biāo)尺的可靠性較低,在此不做詳細(xì)討論。

    圖5 馬溝洞石筍δ18O曲線與格陵蘭冰芯δ18O及董哥洞、葫蘆洞、九仙洞、三寶洞和東石崖洞石筍δ18O曲線對比圖Fig.5 The comparison ofδ18O record from Magou Cave (the blue line and the magenta line are 9-point running average) and the NGRIPδ18O record (green, Rasmussen et al, 2006; Vinther et al, 2006) and theδ18O profi les of stalagmites from the Dongge (sea blue, Dykoski et al, 2005; Wang et al, 2005), Hulu (brown, Wang et al, 2001), Jiuxian (light green, Cai et al, 2010), Sanbao (red, Wang et al, 2008), Dongshiya (cyan, Cai et al, 2008) caves in Asian monsoon region

    (1)11.2 — 9.1 ka BP,石筍δ18O值由老至新呈現(xiàn)較明顯的減小趨勢,變化范圍-11.43‰ —-9.7‰。該時段石筍δ18O持續(xù)降低指示的夏季風(fēng)加強(qiáng)過程被幾個百年尺度的弱季風(fēng)事件打斷。MG-40δ18O記錄的弱季風(fēng)事件年代中心點為11.2 ka BP、10.6 ka BP、9.8 ka BP和9.1 ka BP,與Cai et al(2008)在東石崖洞記錄中發(fā)現(xiàn)的弱季風(fēng)事件在誤差范圍內(nèi)有很好的對應(yīng),進(jìn)一步證實了早全新世東亞季風(fēng)氣候的不穩(wěn)定性。

    (2)9.1 — 4.9 ka BP,該階段石筍δ18O變化趨勢總體上比較穩(wěn)定,平均值為-10.51‰,其中9.1 — 8.6 ka BP存在長達(dá)約500年的弱季風(fēng)(夏季降水減少)時期。至5 ka BP,δ18O值一直處于較偏負(fù)水平,5 — 4.9 ka BPδ18O逐漸增加,指示了一個季風(fēng)降水減少的事件。此后該石筍停止生長,可能指示了夏季風(fēng)降水的系統(tǒng)性減少。

    3.3 YD、9.3 ka、8.2 ka氣候突變事件的區(qū)域特征

    已有研究結(jié)果表明(Wang et al,2001;Denton et al,2010),由北大西洋大量淡水的注入引發(fā)的YD事件(12.9 — 11.5 ka BP)是末次冰期向全新世轉(zhuǎn)換、急劇升溫的過程中一次快速降溫事件,該時期東亞地區(qū)出現(xiàn)千年尺度亞洲夏季風(fēng)減弱的狀況。Ma et al(2012)和Liu et al(2013)通過中國中北部石筍(青天洞、苦傈樹洞)230Th和紋層計數(shù)定年基礎(chǔ)上的δ18O記錄,認(rèn)為YD事件開始于約12.9 ka BP,歷經(jīng)約340年,而在11.5 ka BP左右YD事件向早全新世轉(zhuǎn)變在小于20年內(nèi)完成,這在NGRIP冰芯δ18O記錄和葫蘆洞都有體現(xiàn)(圖5)。YD事件時東石崖DSY09石筍出現(xiàn)沉積間斷風(fēng)化層,馬溝洞石筍δ18O值偏正達(dá)-8.5‰,反映了區(qū)域夏季風(fēng)降水減少、氣候較干旱的狀況。

    9.3 ka左右事件是一個百年尺度上的顯著弱季風(fēng)事件,已被多數(shù)學(xué)者證實(Fleitmann et al,2003;Dykoski et al,2005;Cheng et al,2012)。Bond事件中記錄有9.4 ka BP冷事件,NGRIP冰芯記錄到該冷事件的時間中心點大致在9.0 ka BP和9.2 ka BP,三寶洞記錄到該冷事件的時間中心點大致在9.1 ka BP和9.3 ka BP,董哥洞、九仙洞記錄到該冷事件的時間中心點大致都在9.1 ka BP,河南地區(qū)東石崖石筍DSY09記錄存在9.2 ka BP和9.4 ka BP弱季風(fēng)事件,馬溝洞MG-40石筍δ18O值也在9.1 ka BP和9.3 ka BP左右快速偏正,這些記錄說明在早全新世9.3 ka BP左右確實存在一個快速震蕩的弱季風(fēng)時期。這一事件可能主要受融冰淡水脈沖作用的影響,發(fā)生冷事件時,北大西洋表層為融冰淡水控制,抑制了北大西洋深水(NADW)的形成,使北大西洋經(jīng)向環(huán)流(AMO)的強(qiáng)度減弱,向北輸送的熱量大為減少,造成北大西洋地區(qū)氣候的變冷和海冰范圍的擴(kuò)大。北大西洋地區(qū)的變冷一方面使得西風(fēng)環(huán)流向南擺動,另一方面通過西伯利亞高壓的加強(qiáng)而增強(qiáng)東亞冬季風(fēng),進(jìn)而導(dǎo)致東亞夏季風(fēng)的減弱和季風(fēng)降水的減少(Sun et al,2012)。

    來自世界各地區(qū)的記錄顯示8.2 ka事件可能是全新世最強(qiáng)的冷干事件。中國南北方石筍δ18O記錄的8.2 ka BP事件與格陵蘭NGRIP冰芯δ18O記錄的8.2 ka事件相一致(圖5),體現(xiàn)了高緯氣候系統(tǒng)與中低緯季風(fēng)之間存在遙相關(guān)關(guān)系。Dong et al(2015)通過分析山西蓮花洞、湖北和尚洞、貴州董哥洞8.2 ka BP事件δ18O變幅,認(rèn)為自北向南δ18O變幅減小說明自北向南干旱程度減輕,馬溝洞8.2 ka BP事件δ18O變幅約為1.9‰,較董哥洞(~0.8‰)變幅增大一倍多,與這種區(qū)域變化特征一致,但筆者認(rèn)為由于區(qū)域氣候等差異存在,中國南北8.2 ka BP事件中石筍碳酸鹽δ18O值的差異是否就指示了不同區(qū)域干旱程度的不同,還需要更多數(shù)據(jù)支持和證實。

    3.4 石筍δ18O記錄的全新世大暖期內(nèi)部的氣候波動

    施雅風(fēng)等(1992)認(rèn)為中國的大暖期出現(xiàn)于8.5 — 3.0 ka BP,估計大暖期溫度比現(xiàn)今高2℃左右,華北等地可能高3℃,華南僅高1℃。大暖期氣候波動階段可劃分為:(1)8.5 — 7.2 ka BP不穩(wěn)定的溫度波動階段;(2)7.2 — 6 ka BP大暖期中穩(wěn)定的暖濕階段(大暖期的鼎盛期);(3)6 — 5 ka BP氣候波動劇烈、環(huán)境較差階段;(4)5 — 3 ka BP大暖期結(jié)束時的不穩(wěn)定波動階段。大暖期并不是持續(xù)穩(wěn)定的溫暖濕潤期,期間存在多次氣候波動和低溫干旱事件,這一觀點已得到大多數(shù)學(xué)者證實(施雅風(fēng)等,1992;何元慶和章典,2003;彭曉瑩等,2005)。

    圖6顯示全新世大暖期,馬溝洞MG-1石筍δ18O記錄在總體穩(wěn)定背景下,存在多個數(shù)十年尺度的波動變化,最大振幅可達(dá)2‰。根據(jù)波動特征,馬溝洞全新世大暖期的δ18O記錄大致可以分為8.0 — 7.3 ka BP不穩(wěn)定的波動階段、7.3 —6.1 ka BP穩(wěn)定階段、6.1 — 5 ka BP氣候劇烈波動階段等三個階段,大致與施雅風(fēng)等(1992)的劃分相一致。為了凸顯研究區(qū)大暖期出現(xiàn)的氣候事件,將8000 — 4881 a BPδ18O曲線按0.4年進(jìn)行等間距線性內(nèi)插,在此基礎(chǔ)上做50年加權(quán)平均,得到13個弱季風(fēng)時期(從老到新依次命名為WM1—WM13),年代中心點從老到新依次為7.7 ka BP、7.6 ka BP、7.35 ka BP、7.1 ka BP、6.9 ka BP、6.6 ka BP、6.24 ka BP、6.1 ka BP、6.0 ka BP、5.7 ka BP、5.5 ka BP、5.3 ka BP、5.2 ka BP。

    圖6 全新世大暖期MG-1δ18O記錄的13次弱季風(fēng)事件Fig.6 The thirteen weak monsoon events (yellow bars, numbered with WM1—WM13) identifi ed in MG-1δ18O record during the Holocene Megathermal

    這些夏季風(fēng)減弱時段中,WM1(7637 —7605 a BP)和WM2(7707 — 7671 a BP)時期MG1石筍δ18O值出現(xiàn)強(qiáng)烈偏正,比8.2 ka事件時還偏正0.31‰,指示該時段夏季風(fēng)強(qiáng)度/降水在數(shù)十年尺度上的大幅振蕩。祁連山敦德冰芯δ18O在7.8 — 7.1 ka BP出現(xiàn)三次降溫事件,與MG-1石筍δ18O記錄的不同之處是在7.3 ka BP降溫幅度最大(章新平和姚檀棟,1993)。位于6.1 ka BP、6.0 ka BP的WM4和WM5可能是大暖期鼎盛期向氣候大幅波動、夏季風(fēng)減弱變化的轉(zhuǎn)折點,貴州七星洞石筍δ18O值在6.0 ka BP也呈現(xiàn)出顯著下降(蔡演軍等,2001)。河南伏牛山南坡南陽洞石筍δ18O記錄顯示在5830 — 5720 a BP和5350 — 5080 a BP時段發(fā)生弱季風(fēng)事件(任曉輝等,2006);內(nèi)蒙古西部阿拉善高原的湖泊沉積物在7 — 5 ka BP,尤其是6.2 ka BP、5.5 ka BP時段,總體孢粉濃度較低,且95%是耐旱生的白刺植物花粉,說明當(dāng)時湖面較低、氣候異常干旱(陳發(fā)虎等,2004);敦德冰芯δ18O在5.9 ka BP、5.3 ka BP明顯偏負(fù),出現(xiàn)顯著的降溫事件,這些記錄進(jìn)一步印證了6.1 —5 ka BP是氣候波動劇烈、環(huán)境變化相對較大的階段。

    3.5 9.1 — 4.9 ka BP氣候變化的周期特征與驅(qū)動因子

    采用REDFIT38功率譜分析軟件(ofac、n50、iwin參數(shù)分別為3.0、2.0、1)對馬溝洞MG-1石筍頂部至198 mm段的δ18O序列(7881 — 4881 a BP,平均分辨率為2 a)進(jìn)行周期分析,在99%置信度水平下通過紅噪聲檢驗的周期有22 a、19 a、8 a、7 a、6 a,在95%置信度水平下通過紅噪聲檢驗的周期有272 a、83 a、36 a、15 a、11 a、10 a、5 a(圖7)。22 a和11 a周期與樹輪Δ14C(大氣中14C的濃度)檢測出來的周期一致,對應(yīng)于太陽黑子活動的22 a周期(Kurths et al,1993),83 a和太陽活動世紀(jì)周期(70 — 100 a的變化周期,Gleissberg周期)吻合。這些檢測出來的周期與亞洲季風(fēng)區(qū)早中全新世其他高分辨率記錄(Cai et al,2008,2010)檢測出的太陽活動周期相一致,說明馬溝洞石筍記錄的短尺度夏季風(fēng)變化與太陽活動密切相關(guān)。5 — 7 a周期成分與現(xiàn)今ENSO活動周期相對應(yīng)。趙景耀等(2014)通過對河南雞冠洞降水-洞穴滴水-洞穴現(xiàn)生碳酸鈣的δ18O變化特征分析,認(rèn)為欒川地區(qū)大氣降水δ18Op對“環(huán)流效應(yīng)”和水汽來源的強(qiáng)弱變化更加敏感,表現(xiàn)為與ENSO循環(huán)的正相關(guān),可能與降水中不存在顯著的溫度效應(yīng)和雨量效應(yīng)等噪音因素影響有關(guān),鑒于洞穴碳酸鈣對大氣降水δ18Op的信號具有良好的繼承性,本文支持這一觀點。前人研究表明(李崇銀,1992;王紹武,2001)地處夏季風(fēng)鋒面區(qū)域的華北地區(qū),夏季降水與ENSO關(guān)系較為明確,El Ni?o年夏季平均降水量一般較少,El Ni?o事件所引起的熱帶和副熱帶大氣環(huán)流的異常,尤其是西太平洋副熱帶高壓的異常是造成華北地區(qū)汛期降水量異常的重要原因之一。另外,通過7881 — 4900 a BP時段馬溝洞石筍MG-1δ18O序列功率譜與大氣Δ14C記錄、NGRIPδ18O序列功率譜對比發(fā)現(xiàn),MG-1δ18O序列與大氣Δ14C記錄共有的周期有83 a、36 a、22 a,與NGRIPδ18O序列共有的周期為56 a,和大氣Δ14C記錄的59 a周期相近。說明該時段十年至百年尺度上研究區(qū)域夏季風(fēng)變化主要受太陽活動及北半球高緯地區(qū)氣候變化影響,年際尺度上受ENSO影響顯著。

    圖7 馬溝洞石筍MG-1δ18O序列功率譜分析結(jié)果(a)及與相同時段(7881 — 4881 a BP)大氣Δ14C記錄(b)、NGRIPδ18O序列(c)功率譜對比Fig.7 Power spectral analysis of the stalagmite MG-1δ18O sequence (a) from Magou Cave and the comparison with the power spectral analyses of atmosphericΔ14C record (b) and the NGRIPδ18O record (c) during the period of 7881 — 4881 a BP

    大氣Δ14C的變化與太陽活動變化相聯(lián)系,可以作為太陽活動的代用指標(biāo)(Dykoski et al,2005)。圖8中馬溝洞石筍MG-1δ18O與大氣Δ14C去趨勢曲線在百年尺度上存在相似變化,說明在百年乃至幾十年尺度上太陽活動是亞洲季風(fēng)強(qiáng)度變化的驅(qū)動因素之一。然而,對比還發(fā)現(xiàn)7429 — 7137 a BP和6067 — 5924 a BP兩者存在反相位變化,7429 — 7137 a BP時段最為明顯,這很難用年代誤差來解釋兩者的差異,可能說明這兩個時段內(nèi)石筍MG-1δ18O波動變化與其他氣候影響因素(如高緯地區(qū)氣候變化、ENSO)有關(guān)。

    MG-1石筍的δ18O記錄15年加權(quán)平均曲線與格陵蘭NGRIP冰芯δ18O記錄的時間變化序列(圖8),兩者在測年誤差內(nèi)有著較好的一致性,即使在9100 — 7881 a BP石筍MG-1的分辨率較低,但也與冰芯δ18O記錄存在同向變化,說明極地氣溫變化與亞洲夏季風(fēng)之間密切聯(lián)系??赡艿脑蚴菢O地高緯度氣溫的變化可以導(dǎo)致西伯利亞高壓的強(qiáng)度發(fā)生變化,并進(jìn)而影響冬季風(fēng)的強(qiáng)度。而馬溝洞地區(qū)對冬季風(fēng)強(qiáng)度的變化有比較敏感的響應(yīng),強(qiáng)盛的冬季風(fēng)使夏季風(fēng)北界南移,同時持續(xù)時間變短。反之,則夏季風(fēng)持續(xù)時間變長,降水增多。

    圖8 距今9100 — 4900年馬溝洞石筍MG-1δ18O記錄與格陵蘭NGRIP冰芯δ18O記錄、大氣Δ14C去趨勢曲線對比Fig.8 Comparisons of the MG-1δ18O record (red, the bold line is the weighted average of 15 points) with the detrended atmosphericΔ14C record (green, Stuiver and Braziunas, 1993) and the NGRIPδ18O record (blue, Rasmussen et al, 2006; Vinther et al, 2006) during the period of 9100 — 4900 a BP.

    根據(jù)上述周期分析和曲線對比,可知太陽活動引起的太陽輻射變化和北半球高緯氣候狀況共同影響著亞洲季風(fēng)的變化,ENSO活動及氣候系統(tǒng)內(nèi)部的相互作用也對東亞夏季風(fēng)降水產(chǎn)生重要的影響。

    3.6 中原地區(qū)新石器文化發(fā)展與區(qū)域氣候和環(huán)境變化的可能聯(lián)系

    自然氣候和環(huán)境變遷對古文化和社會發(fā)展的影響是不可忽視的。若自然環(huán)境變化的幅度較大,超出了人類適應(yīng)能力的范圍,環(huán)境變化將可能對文化產(chǎn)生一種抑制作用,有時甚至直接導(dǎo)致文化的衰落(吳文祥和葛全勝,2005)。根據(jù)河南地區(qū)現(xiàn)有考古發(fā)掘及相關(guān)研究(魏書亞和宋國定,1997;王幼平,2014;李開封等,2015),河南地區(qū)早中全新世環(huán)境考古時代主要包括:新石器時代的李家溝文化時期(10500 — 8600 cal a BP)、裴李崗文化時期(8500 — 7000 cal a BP)、仰韶文化時期(7000 — 5000 cal a BP)和龍山文化時期(5000 — 4000 cal a BP)(由于文化過渡期問題,不同文化遺存時代上也有交錯現(xiàn)象,本文中不同文化的斷代以李開封等(2015)最新數(shù)據(jù)為準(zhǔn))。

    7.3ka BP處于河南地區(qū)裴李崗文化與仰韶文化的過渡期,氣候由不穩(wěn)定的波動階段轉(zhuǎn)向相對穩(wěn)定的暖濕階段,同時7200 — 6000 a BP仰韶文化在北方地區(qū)雨后春筍般地涌現(xiàn)期(吳文祥和葛全勝,2005)與馬溝洞石筍MG-1δ18O記錄的7.3 —6.1 ka BP穩(wěn)定的暖濕階段相對應(yīng),這可能說明相對穩(wěn)定且適宜的氣候和環(huán)境條件為仰韶文化在中原地區(qū)及中國北方地區(qū)的繁榮和發(fā)展提供了良好的環(huán)境背景,促進(jìn)了仰韶文化的發(fā)展和繁榮。

    4 結(jié)論

    (1)基于河南省禹州市馬溝洞2支石筍(MG-1與MG-40)24個230Th年齡和1988個氧同位素數(shù)據(jù),建立了研究區(qū)13.1 — 4.9 ka BP期間分辨率為2 — 14 a的石筍氧同位素時間序列。MG-1和MG-40在相同時段內(nèi)(10.37 — 7.06 ka BP)變化基本一致,且與周邊洞穴石筍記錄可以對比,滿足重復(fù)性檢驗,表明馬溝洞石筍是在同位素平衡分餾狀態(tài)下沉積的。馬溝洞石筍碳酸鹽δ18O變化主要指示了季風(fēng)降水的變化,溫度變化的影響較小。

    (2)馬溝洞石筍δ18O記錄了11.2 — 9.1 ka BP區(qū)域季風(fēng)降水在波動中逐漸增加,9.1 — 4.9 ka BP季風(fēng)降水顯著波動但無明顯長期趨勢的演變過程。YD事件、9.3 ka事件、8.2 ka事件記錄與其他石筍δ18O記錄的一致性揭示末次冰消期百年—千年尺度氣候突變事件的大范圍存在和共同的驅(qū)動因子。MG-1石筍δ18O記錄的13個數(shù)十年尺度的弱季風(fēng)事件支持了大暖期并不是持續(xù)穩(wěn)定的溫暖濕潤期,期間存在多次降水波動事件這一觀點。

    (3)在7881 — 4900 a BP時段,馬溝洞石筍MG-1δ18O序列與大氣Δ14C記錄共有的周期有83 a、36 a、22 a,與NGRIPδ18O序列共有的周期為56 a,同時MG-1δ18O序列存在5 — 7 a顯著周期,而且9100 — 4881 a BP時間段MG-1δ18O記錄與大氣Δ14C去趨勢曲線變化大體一致、與格陵蘭NGRIP冰芯δ18O記錄的時間變化序列吻合性也較好,這表明太陽活動引起的太陽輻射變化和北半球高緯氣候狀況共同影響著亞洲季風(fēng)的變化,ENSO活動及氣候系統(tǒng)內(nèi)部的相互作用也對東亞夏季風(fēng)降水產(chǎn)生重要的影響。

    (4)仰韶文化繁榮期(7200 a BP — 6000 a BP)與馬溝洞石筍MG-1δ18O記錄的7.3 — 6.1 ka BP穩(wěn)定的暖濕階段相對應(yīng),這可能說明相對穩(wěn)定且適宜的氣候和環(huán)境條件為仰韶文化在中原地區(qū)及中國北方地區(qū)的繁榮和發(fā)展提供了良好的環(huán)境背景,促進(jìn)了仰韶文化的發(fā)展和繁榮。

    致謝:感謝中國科學(xué)院地球環(huán)境研究所周杰同學(xué)在圖形繪制方面給予的幫助。

    蔡演軍, 彭子成, 安芷生, 等. 2001. 貴州七星洞全新世石筍的氧同位素記錄及其指示的季風(fēng)氣候變化[J].科學(xué)通報, 46(16): 1398 – 1402. [Cai Y J, Peng Z C, An Z S, et al. 2001. Theδ18O variation of a stalagmite from Qixing Cave, Guizhou Province and indicated climate change during the Holocene [J].Chinese Science Bulletin, 46(16): 1398 – 1402.]

    陳發(fā)虎, 吳 薇, 朱 艷, 等. 2004. 阿拉善高原中全新世干旱事件的湖泊記錄研究[J].科學(xué)通報, 49(1): 1 – 9. [Chen F H, Wu W, Zhu Y, et al. 2004. Research on the lake record of Holocene drought events in the Alashan Plateau [J].Chinese Science Bulletin, 49(1): 1 – 9.]

    何元慶, 章 典. 2003. 冰芯與其他記錄所揭示的中國全新世大暖期變化特征[J].冰川凍土, 25(1): 11 – 18. [He Y Q, Zhang D. 2003. Climatic Differences in China during the Holocene Indicated by the Various Climatic Proxy Data from Different Parts of China [J].Journal of glaciology and geocryology, 25(1): 11 – 18.]

    李崇銀. 1992. 華北地區(qū)汛期降水的一個分析研究[J].氣象學(xué)報, 50(1): 41 – 49. [Li C Y. 1992. An analytical study on the precipitation in the fl ood period over Huabei area [J].Acta Geographica Sinica, 50(1): 41 – 49.]

    李開封, 馬春梅, 高文華, 等. 2015. 河南地區(qū)全新世環(huán)境考古研究進(jìn)展及未來趨勢[J].地理科學(xué)進(jìn)展, 34(7): 883 – 897. [Li K F, Ma C M, Gao W H, et al. 2015. Progress and trend of Holocene environmental archaeology in Henan Province [J].Progress in Geography, 34(7): 883 – 897.]

    李中軒, 朱 誠, 吳國璽, 等. 2013. 河南省史前人類遺址的時空分布及其驅(qū)動因子[J].地理學(xué)報, 68(11): 1527 – 1537. [Li Z X, Zhu C, Wu G X, et al. 2013. Spatial and temporal distribution of prehistoric human sites and its driving factors in Henan Province [J].Acta Geographica Sincia, 68(11): 1527 – 1537.]

    馬 樂, 蔡演軍, 秦世江, 等. 2015. 貴州七星洞石筍記錄的最近2300年氣候和環(huán)境變化[J].地球環(huán)境學(xué)報, 6(3): 135 – 144. [Ma L, Cai Y J, Q S J, et al. 2015. A high resolution paleoclimate record of the last 2300 years in stalagmite QX-3 from the Qixing Cave, Guizhou Province [J].Journal of Earth Environment, 6(3): 135 – 144.]

    馬 睿. 2013. 河南欒川雞冠洞現(xiàn)代監(jiān)測與巖溶系統(tǒng)碳同位素運移研究[D]. 重慶: 西南大學(xué). [Ma R. 2013. Carbon isotope migration in karst system and Modern monitoring of Jiguan cave, Luanchuan City, Henan Province [D]. Chongqing: Southwest University.]

    彭曉瑩, 鐘 巍, 趙引娟, 等. 2005. 全新世大暖期氣候環(huán)境特征及其機(jī)制的再認(rèn)識[J].華南師范大學(xué)學(xué)報: 自然科學(xué)版, (2): 52 – 60. [Peng X Y, Zhong W, Zhao Y J, et al. 2005. Preliminary research in the Holocene Hypsithermal climate change [J].Journal of South China Normal University: Natural Science Edition, (2): 52 – 60.]

    任曉輝, 董進(jìn)國, 陳昌海, 等. 2006. 中全新世河南南陽石筍記錄的百年季風(fēng)干旱事件[J].中國巖溶, 25(4): 269 – 273. [Ren X H, Dong J G, Chen C H, et al. 2006. Centennial-scale Asian monsoonal aridity events during the Mid-Holocene inferred by stalagmite record from Nanyang Cave, Henan Province [J].Carsologica Sinica, 25(4): 269 – 273.]

    任小鳳, 楊 琰, 彭 濤, 等. 2014. 豫西雞冠洞洞穴水及現(xiàn)代沉積物Mg, Sr和Ba記錄及其意義[J].中國巖溶, 33(1): 57 – 63. [Ren X F, Yang Y, Peng T, et al. 2014. The records and implications of Mg, Sr and Ba in cave water and active speleothems in Jiguan Cave, western Henan Province [J].Carsologica Sinica, 33(1): 57 – 63.]

    施雅風(fēng), 孔昭宸, 王蘇民, 等. 1992. 中國全新世大暖期的氣候波動與重要事件[J].中國科學(xué):化學(xué)·生命科學(xué)·地學(xué), 22(12): 1300 – 1308. [Shi Y F, Kong Z C, Wang S M, et al. 1992. Climate changes and abrupt climate change events during the Holocene Megathermal in China [J].Scientia Sinica Chimica: Chemistry, Life Sciences, Earth Sciences, 22(12): 1300 – 1308.]

    王紹武. 2001.現(xiàn)代氣候?qū)W研究進(jìn)展[M]. 北京:氣象出版社: 161 – 165. [Wang S W. 2001. Progress in modern climatology [M]. Beijing: Meteorological Press: 161 – 165.]

    王幼平. 2014. 新密李家溝遺址研究進(jìn)展及相關(guān)問題[J].中原文物, (1): 20 – 24. [Wang Y P. 2014. On the research progress of Lijiagou site in Xinmi City and other related problems [J].Cultural Relics of Central China, (1): 20 – 24.]

    魏書亞, 宋國定. 1997. 從放射性碳素斷代看中原地區(qū)新石器文化的幾個問題[J].中原文物,(2): 100 – 103, 112. [Wei S Y, Song G D. 1997. Some problems of the Neolithic Culture in Central China from the radiocarbon dating [J].Cultural Relics of Central China, (2): 100 – 103, 112.]

    吳文祥, 葛全勝. 2005. 全新世氣候事件及其對古文化發(fā)展的影響[J].華夏考古, (3): 60 – 67. [Wu W X, Ge Q S. 2005. Climate events of the recent epoch and their influence upon the development of ancient culture [J].Huaxia Archaeology, (3): 60 – 67.]

    章新平, 姚檀棟. 1993. 祁連山敦德冰帽冰芯中氣候記錄的綜述[J].沙漠與綠洲氣象, 16(6): 1 – 6. [Zhang X P, Yao T D. 1993. Review: Climatic records in the Qilian Dunde ice core [J].Desert and Oasis Meteorology, 16(6): 1 – 6.]

    張銀環(huán), 楊 琰, 楊勛林, 等. 2015. 早全新世季風(fēng)演化的高分辨率石筍δ18O記錄研究——以河南老母洞石筍為例[J].沉積學(xué)報, 33(1). 134 – 139. [Zhang Y H, Yang Y, Yang X L, et al. 2015. Early Holocene Monsoon Evolution of High-resolution Stalagmiteδ18O Records: in Henan Laomu Cave [J].Acta Sedimentologica Sinica, 33(1). 134 – 139.]

    趙景耀, 楊 琰, 彭 濤, 等. 2014. 河南雞冠洞降水、滴水和現(xiàn)生碳酸鈣的δ18O變化特征及其環(huán)流意義[J].第四紀(jì)研究, 34(5): 1106 – 1116. [Zhao Y J, Yang Y, Peng T, et al. 2014. Variation ofδ18O values in the precipitation, cave drip water and modern calcite deposition in Jiguan cave, Henan Province and its atmospheric circulation effect [J].Quaternary Sciences, 34(5): 1106 – 1116.]

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

    Cai Y J, Tan L C, An Z S, et al. 2010. The variation of summer monsoon precipitation incentral China since the last deglaciation [J].Earth and Planetary Science Letters, 291: 23 – 31.

    Cheng H, Edwards R L, Shen C C, et al. 2013. Improvements in230Th dating,230Th and234U half-life values, and U-Th isotopic measurements by multi-collector inductively coupled plasma mass spectrometry [J].Earth and Planetary Science Letters, 371: 82 – 91.

    Cheng H, Sinha A, Wang X, et al. 2012. The Global Paleomonsoon as seen through speleothem records from Asia and the Americas [J].Climate Dynamics, 39(5): 1045 – 1062.

    Denton G H, Anderson R F, Toggweiler J R, et al. 2010. The last glacial termination [J].Science, 328(5986): 1652 – 1656.

    Dong J G, Shen C C, Kong X G, et al. 2015. Reconciliation of hydroclimate sequences from the Chinese Loess Plateau and low-latitude East Asian Summer Monsoon regions over the past 14,500 years [J].Palaeogeography, Palaeoclimatology, Palaeoecology, 435: 127 – 135.

    Dorale J A, Edwards R L, Ito E, et al. 1998. Climate and vegetation history of the midcontinent from 75 to 25 ka: a speleothem record from Crevice Cave, Missouri, USA [J].Science, 282(5395): 1871 – 1874.

    Dykoski C A, Edwards R L, Cheng H, et al. 2005. A highresolution, absolute-dated Holocene and deglacial Asian monsoon record from Dongge Cave, China [J].Earth and Planetary Science Letters, 233: 71 – 86.

    Edwards R L, Chen J H, Ku T L, et al. 1987. Precise timing of the last interglacial period from mass spectrometric determination of thorium-230 in corals [J].Science, 236(4808): 1547 – 1553.

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

    Johnson K R, Ingram B L. 2004. Spatial and temporal variability in the stable isotope systematics of modern precipitation in China: implications for paleoclimate reconstructions [J].Earth and Planetary Science Letters, 220(3): 365 – 377.

    Johnson K R, Lynn Ingram B, Sharp W D, et al. 2006. East Asian summermonsoon variability during marine isotope stage 5 based on speleothemδ18O records from Wanxiang Cave, central China [J].Palaeogeography, Palaeoclimatology, Palaeoecology, 236: 5 – 19.

    Kim S T, O'Neil J R. 1997. Equilibrium and nonequilibrium oxygen isotope effects in synthetic carbonates [J].Geochimica et Cosmochimica Acta, 61(16): 3461 – 3475.

    Kurths J, Spiering Ch, Muller-Stoll W, et al. 1993. Search for solar periodicities in Miocene tree ring widths [J].Terra Nova, 5: 359 – 363.

    Liu D, Wang Y, Cheng H, et al. 2013. Centennial-scale Asian monsoon variability during the mid-Younger Dryas from Qingtian Cave, central China [J].Quaternary Research, 80(2): 199 – 206.

    Ma Z B, Cheng H, Tan M, et al. 2012. Timing and structure of the Younger Dryas event in northern China [J].Quaternary Science Reviews, 41: 83 – 93.

    McCrea J M. 1950. On the isotopic chemistry of carbonates and a paleotemperature scale [J].The Journal of Chemical Physics, 18(6): 849 – 857.

    Rasmussen S O, Andersen K K, Svensson A M, et al. 2006. A new Greenland ice core chronology for the last glacial termination [J].Journal of Geophysical Research: Atmospheres (1984—2012), 111, D06102, doi: 10.1029/2005JD006079.

    Stuiver M, Braziunas T F. 1993. Modeling atmospheric14C infl uences and14C ages of marine samples to 10,000 BC [J].Radiocarbon, 35: 137–189.

    Sun Y, Clemens S C, Morrill C, et al. 2012. Infl uence of Atlantic meridional overturning circulation on the East Asian winter monsoon [J].Nature Geoscience, 5(1): 46 – 49.

    Vinther B M, Clausen H B, Johnsen S J, et al. 2006. A synchronized dating of three Greenland ice cores throughout the Holocene [J].Journal of Geophysical Research: Atmospheres (1984—2012), 111(D13), doi: 10.1029/2005JD006921.

    Wang Y J, Cheng H, Edwards R L, et al. 2001. A highresolution absolute-dated Late pleistocene monsoon record from Hulu Cave, China [J].Science, 294(5550): 2345 – 2348.

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

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

    Early to mid-Holocene paleoclimatic changes recorded by the stalagmites from the Magou Cave, Henan Province

    MAO Ruixue1,2, CAI Yanjun1,3, MA Le1, CHENG Xing1,2
    (1. State Key Laboratory of Loess and Quaternary Geology, Institute of Earth Environment, Chinese Academy of Sciences, Xi'an 710061, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. Institute of Global Environmental Change, Xi'an Jiaotong University, Xi'an 710054, China)

    Background, aim, and scopeHenan Province, located in central China, is dominated by monsoonal climate, and has been regarded as the core region of Chinese ancient civilization. Stalagmite has advantages of absolute and high-precision chronology and abundant proxy indexes, and has been widely used to decipher the changes of Asian summer in the past. However, there are few stalagmite records obtained from this region till now and all of them are located in the hilly area of west Henan Province. Here, we present a high-resolution and precisely dated stalagmiteδ18O records attained from two stalagmites (MG-1 and MG-40) from Magou Cave, northern Henan Province, Central China. We reconstruct the history of monsoon precipitation during early to mid Holocene, characterize the climate instability of the Holocene Megathermal and investigate the possible driving forces and dynamics of the Asian Summer Monsoon (ASM) in the study region. We also examine the possible links existed between civilization succession and climate change.Materials and methodsTwo calcite stalagmites, MG-1 and MG-40, have been analyzed in this study. Stalagmite MG-1 is ~26 cm in height and ~9 cm in diameter,and MG-40 is ~4 cm in height and ~8 cm in diameter. We halved these two stalagmites using diamond wire saw. Twenty-four subsamples were drilled parallel to the growth plane using a hand-held carbide dental drill, and were dated with U-series methods on a plasma sourced multi-collector mass spectrometry (MC-ICP-MS) at the Institute of Global Environmental Change, Xi'an Jiaotong University. The chemical procedure is similar to that described in Edwards et al (1987), and the details on instrumental approaches were explained in Cheng et al (2013). Here we use a linear interpolation to establish the chronology, and extrapolate the ages for the portions outside the U-series dates using the growth rate of adjacent segments. Subsamples for stable isotope analysis were collected with a high resolution computer-controlled micro-milling equipment (Micro-mill). A total of 1988 oxygen and carbon isotopic values were obtained on an IsoPrime100 mass spectrometer equipped with the Dual-inlet and connected to the MultiPrep systems in the Institute of Earth Environment, Chinese Academy of Sciences. International standard NBS-19 and inter-laboratory standard TTB1 were run every 10 to 15 samples and arbitrary selected duplicates were run every 10 to 20 samples to check for homogeneity and reproducibility. All oxygen isotopic values are reported inδnotation, the per mil deviation relative to the Vienna Pee Dee Belemnite (VPDB) standard. The standard results show that the precision ofδ18O analysis is better than 0.08‰ (2σ).ResultsAll230Th dates are all in stratigraphic order and it appears that MG-1 and MG-40 grew continuously throughout this period. The calculated growth rate of stalagmites MG-1 and MG-40 varied from 4.7 μm · a-1to 58.6 μm · a-1and from 1.5 μm · a-1to 6.3 μm · a-1respectively. The growth rates of both stalagmites are similar in general, they are relatively low during the early Holocene, then increasing and holding steady. The 1988 isotopic data obtained generate a 2 — 14 years resolution to the oxygen isotope profi le during the 13.1 — 4.9 ka BP. A careful examination of the time-series of speleothemδ18O ratios fi nds that theδ18O of stalagmite MG-1 fluctuated around -10.5‰ (mean value) from 9.1 to 4.9 ka BP, varying from -12.7‰ to -9.2‰. Theδ18O of stalagmite MG-40 fl uctuated around -10.4‰ (mean value) from 11.1 to 7.1 ka BP, varying from -11.9‰ to -8.4‰. The Magouδ18O profile can be divided into three periods: (1) 13.1 — 11.2 ka BP, which possible correspond to Aller?d Warming and Younger Dryas; (2) a gradual increasing precipitation interval between 11.2 — 9.1 ka BP; and (3) a relatively stable interval between 9.1 —4.9 ka BP. Theδ18O value fl uctuated around -10.51‰ (mean value) from 9.1 to 4.9 ka BP, corresponding to the Holocene Megathermal.DiscussionBased on instrumental observations (winter precipitation accounts for only ~37.7% of the amount of annual precipitation in this region) and previous stalagmite studies in this region and southern China, we interpret the stalagmiteδ18O from Magou Cave as a qualitative proxy of ASM precipitation, with heavier values indicating a weak ASM and less summer monsoon precipitation, and vice versa. The Younger Dryas, 9.3 ka and 8.2 ka events can be identified in the Magouδ18O profile and are in good agreement with other speleothem records from Asian monsoon regions. Thirteen signifi cant weak monsoon events on centennial to decadal timescales were singled out within the period of 8.0 —4.9 ka BP. The comparisons between MG-1δ18O record and atmosphericΔ14C record and the NGRIPδ18O record (during the period of 9.1 — 4.9 ka), and the spectrum analysis reveal largely consistent changes among these records.ConclusionsThe occurrence of Younger Dryas, 9.3 ka and 8.2 ka events in the Magouδ18O profile are in good agreement with other speleothem records from Asian monsoon regions, suggesting that these abrupt events have affected the vast monsoon regions at least in East Asian and South Asia. Thirteen signifi cant centennial-decadal weak monsoon events identifi ed within the period of 8.0 — 4.9 ka BP, confi rm the instability of monsoon climate during the Holocene Megathermal. The concordance of MG-1δ18O record and atmosphericΔ14C record and the NGRIPδ18O record (during the period of 9.1 — 4.9 ka), as well as the spectrum analysis, reveal that not only did both changes in solar radiation and climatic changes at high northern latitude affect the Asian summer monsoon variations during the Holocene, bout also ENSO and the internal interaction of climate sub-systems have signifi cant impacts on East Asian summer monsoon changes. We notethat the fl ourishing of the Yangshao Culture at the central China during the 7200 — 6000 a BP was coincident with the relatively stable warm-humid climate within the Holocene Megathermal. It is then not implausible that the climate optimum plays an important role in the development of earlier civilizations (e.g., Yangshao Culture).Recommendations and perspectivesThis paper provides a new high resolution stalagmite record for palaeoenvironment and environmental archaeological research of Henan Province from early to mid-Holocene.

    stalagmite;δ18O ; climate change; summer monsoon precipitation; Holocene; Magou Cave

    CAI Yanjun, E-mail: yanjun_cai@ieecas.cn

    10.7515/JEE201603004

    2016-01-06;錄用日期:2016-04-04

    Received Date:2016-01-06;Accepted Date:2016-04-04

    中國科學(xué)院重點部署項目(KZZD-EW-04);國家自然科學(xué)基金項目(41271229)

    Foundation Item:Key Deployment Project of Chinese Academy of Sciences (KZZD-EW-04); National Natural Science Foundation of China (41271229)

    蔡演軍,E-mail: yanjun_cai@ieecas.cn

    猜你喜歡
    石筍季風(fēng)同位素
    綠水青山圖——海洋季風(fēng)的贊歌
    戶撒刀
    石筍功
    深空探測用同位素電源的研究進(jìn)展
    萬馬奔騰
    中華奇石(2015年9期)2015-07-09 18:32:58
    游石筍山(外一首)
    岷峨詩稿(2014年4期)2014-11-15 00:24:03
    《同位素》(季刊)2015年征訂通知
    同位素(2014年3期)2014-06-13 08:22:28
    硼同位素分離工藝與生產(chǎn)技術(shù)
    同位素(2014年2期)2014-04-16 04:57:15
    穩(wěn)定同位素氘標(biāo)記蘇丹紅I的同位素豐度和化學(xué)純度分析
    同位素(2014年2期)2014-04-16 04:57:12
    客鳥尾石筍奇觀
    成人三级黄色视频| 最近的中文字幕免费完整| 亚洲美女搞黄在线观看| 日本熟妇午夜| av免费观看日本| 18禁在线播放成人免费| 一进一出抽搐动态| 久久人妻av系列| 国语自产精品视频在线第100页| 自拍偷自拍亚洲精品老妇| 亚洲国产欧美人成| 国产老妇伦熟女老妇高清| av黄色大香蕉| 国产白丝娇喘喷水9色精品| 免费看a级黄色片| 国产一区二区在线av高清观看| 99久久精品国产国产毛片| 在线播放国产精品三级| av福利片在线观看| 丰满的人妻完整版| 好男人视频免费观看在线| 老司机福利观看| 亚洲成av人片在线播放无| 国产精品久久久久久久电影| 国内少妇人妻偷人精品xxx网站| av免费观看日本| 一级毛片久久久久久久久女| 又爽又黄无遮挡网站| 97超碰精品成人国产| 中文字幕av成人在线电影| 天堂√8在线中文| 禁无遮挡网站| 国产精品电影一区二区三区| 男女啪啪激烈高潮av片| 九九热线精品视视频播放| 亚洲成人久久爱视频| 精品久久久噜噜| 高清午夜精品一区二区三区 | 一进一出抽搐动态| 成人午夜精彩视频在线观看| 久久久午夜欧美精品| 国产精品日韩av在线免费观看| 97热精品久久久久久| 亚洲精品乱码久久久v下载方式| 国产毛片a区久久久久| 哪个播放器可以免费观看大片| 青春草亚洲视频在线观看| 国产成人aa在线观看| 亚洲三级黄色毛片| 91久久精品电影网| 97在线视频观看| 最近手机中文字幕大全| 国产精品麻豆人妻色哟哟久久 | 成人永久免费在线观看视频| 看非洲黑人一级黄片| 两个人视频免费观看高清| 精品一区二区三区视频在线| 国产成人午夜福利电影在线观看| 久久久午夜欧美精品| 日韩强制内射视频| 桃色一区二区三区在线观看| 国产精品人妻久久久久久| 欧美+亚洲+日韩+国产| 九九热线精品视视频播放| 国产中年淑女户外野战色| a级一级毛片免费在线观看| 1000部很黄的大片| 长腿黑丝高跟| 欧美色欧美亚洲另类二区| 国产伦精品一区二区三区四那| 91久久精品电影网| 乱人视频在线观看| 男人狂女人下面高潮的视频| 一个人看的www免费观看视频| 草草在线视频免费看| 毛片一级片免费看久久久久| 国产乱人偷精品视频| 一级黄色大片毛片| 精品一区二区免费观看| 男人的好看免费观看在线视频| 又爽又黄a免费视频| 国产高清不卡午夜福利| 国产精品人妻久久久影院| 国产精华一区二区三区| 欧美潮喷喷水| 中文欧美无线码| 一级毛片久久久久久久久女| 亚洲欧美精品综合久久99| 成人午夜精彩视频在线观看| av在线亚洲专区| 一进一出抽搐动态| 色尼玛亚洲综合影院| 久久亚洲国产成人精品v| 国产乱人视频| 男女边吃奶边做爰视频| 一边摸一边抽搐一进一小说| 国产精品一二三区在线看| 欧美最黄视频在线播放免费| 亚洲一区二区三区色噜噜| 爱豆传媒免费全集在线观看| 六月丁香七月| 秋霞在线观看毛片| 丝袜喷水一区| 日本成人三级电影网站| 22中文网久久字幕| 色吧在线观看| 免费观看的影片在线观看| av在线老鸭窝| 三级男女做爰猛烈吃奶摸视频| 看非洲黑人一级黄片| 日日摸夜夜添夜夜爱| 欧美色视频一区免费| 国产成人精品一,二区 | 黑人高潮一二区| 亚洲国产精品成人综合色| 国产美女午夜福利| 久久人妻av系列| 欧美性猛交黑人性爽| 国产老妇伦熟女老妇高清| 日本三级黄在线观看| 亚洲,欧美,日韩| 久久国产乱子免费精品| 自拍偷自拍亚洲精品老妇| 欧美日韩在线观看h| 亚洲无线在线观看| 女人被狂操c到高潮| 国产精品一区二区三区四区免费观看| 中文精品一卡2卡3卡4更新| 国产真实乱freesex| 日韩欧美在线乱码| www日本黄色视频网| 精品久久久久久久末码| 青春草视频在线免费观看| 此物有八面人人有两片| 免费搜索国产男女视频| 青春草视频在线免费观看| 日日啪夜夜撸| videossex国产| 国内精品久久久久精免费| 亚洲激情五月婷婷啪啪| 69av精品久久久久久| 午夜久久久久精精品| 国产熟女欧美一区二区| 91久久精品国产一区二区三区| 成人鲁丝片一二三区免费| 国产精品一区二区三区四区免费观看| 麻豆国产97在线/欧美| 桃色一区二区三区在线观看| 国产精品久久久久久久电影| 亚洲精品国产成人久久av| 成人二区视频| 麻豆久久精品国产亚洲av| 有码 亚洲区| 欧美3d第一页| 永久网站在线| 国产黄色小视频在线观看| 亚洲国产精品国产精品| 天天躁夜夜躁狠狠久久av| 欧美性猛交╳xxx乱大交人| 成熟少妇高潮喷水视频| 国产女主播在线喷水免费视频网站 | 亚洲成a人片在线一区二区| 亚洲国产精品合色在线| 波多野结衣高清作品| 99热全是精品| 国产精品一区二区在线观看99 | 寂寞人妻少妇视频99o| 久久久久久大精品| 高清毛片免费看| 日本黄大片高清| 中国国产av一级| 欧美极品一区二区三区四区| 国产爱豆传媒在线观看| 免费人成视频x8x8入口观看| 国产探花在线观看一区二区| 亚洲国产精品sss在线观看| 一级毛片我不卡| 看黄色毛片网站| 午夜a级毛片| 精品一区二区三区视频在线| h日本视频在线播放| 亚洲av成人av| 亚洲精品乱码久久久v下载方式| 亚洲欧美精品自产自拍| 欧美一区二区亚洲| 丝袜喷水一区| 在线观看一区二区三区| 久久精品夜色国产| 亚洲成人av在线免费| 麻豆国产97在线/欧美| 极品教师在线视频| 91久久精品电影网| 亚洲高清免费不卡视频| 中文资源天堂在线| 久久中文看片网| 非洲黑人性xxxx精品又粗又长| 免费无遮挡裸体视频| 国产男人的电影天堂91| 晚上一个人看的免费电影| 国产极品精品免费视频能看的| or卡值多少钱| 国产精品永久免费网站| 自拍偷自拍亚洲精品老妇| 99riav亚洲国产免费| 老师上课跳d突然被开到最大视频| 欧美不卡视频在线免费观看| 午夜福利在线观看免费完整高清在 | 久久这里有精品视频免费| 天天躁夜夜躁狠狠久久av| 亚洲最大成人av| 男插女下体视频免费在线播放| 国产综合懂色| 插逼视频在线观看| 久久久成人免费电影| 男插女下体视频免费在线播放| a级一级毛片免费在线观看| 色哟哟哟哟哟哟| 国产av不卡久久| 亚洲欧美中文字幕日韩二区| 黄色欧美视频在线观看| 欧美人与善性xxx| 村上凉子中文字幕在线| 一边亲一边摸免费视频| 久久午夜亚洲精品久久| 老师上课跳d突然被开到最大视频| 人人妻人人澡人人爽人人夜夜 | 久久99热6这里只有精品| 日韩一区二区视频免费看| 亚洲在线自拍视频| 亚洲最大成人中文| 中文欧美无线码| 欧美性猛交╳xxx乱大交人| 三级毛片av免费| 日韩成人伦理影院| 久久99蜜桃精品久久| 一级二级三级毛片免费看| 乱系列少妇在线播放| 老师上课跳d突然被开到最大视频| 国产熟女欧美一区二区| 69人妻影院| 联通29元200g的流量卡| 亚洲综合色惰| 亚洲av不卡在线观看| 1000部很黄的大片| 99riav亚洲国产免费| 国产伦在线观看视频一区| 亚洲18禁久久av| 中文字幕制服av| 国内少妇人妻偷人精品xxx网站| 国产成人a∨麻豆精品| 在线国产一区二区在线| 日韩三级伦理在线观看| 亚洲图色成人| 少妇被粗大猛烈的视频| 久久久久性生活片| 人妻夜夜爽99麻豆av| 亚洲欧美日韩卡通动漫| 97超碰精品成人国产| 亚洲精品色激情综合| 久久草成人影院| 在线观看av片永久免费下载| 免费观看的影片在线观看| 亚洲最大成人中文| 美女xxoo啪啪120秒动态图| 日本一本二区三区精品| www.色视频.com| 男人狂女人下面高潮的视频| 蜜臀久久99精品久久宅男| videossex国产| 欧美激情国产日韩精品一区| 国产精品爽爽va在线观看网站| 国产色爽女视频免费观看| 精品久久久久久久久亚洲| 一进一出抽搐动态| 亚洲成人久久性| 午夜福利在线观看吧| 欧美+日韩+精品| 成年女人永久免费观看视频| 日韩视频在线欧美| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 两个人的视频大全免费| 午夜福利视频1000在线观看| 亚洲成人av在线免费| 精品免费久久久久久久清纯| 亚洲欧美成人精品一区二区| 欧美一区二区国产精品久久精品| 免费人成在线观看视频色| www.色视频.com| 久久草成人影院| 男人舔奶头视频| 久久九九热精品免费| 免费无遮挡裸体视频| 久久综合国产亚洲精品| 国产精品不卡视频一区二区| 淫秽高清视频在线观看| 精品久久久久久久人妻蜜臀av| 97超视频在线观看视频| 又爽又黄a免费视频| 婷婷色综合大香蕉| 性欧美人与动物交配| 日韩一本色道免费dvd| 国产成人影院久久av| 亚洲欧美日韩东京热| 成年女人永久免费观看视频| av专区在线播放| 色吧在线观看| 国产亚洲av片在线观看秒播厂 | 看十八女毛片水多多多| 狠狠狠狠99中文字幕| 中国美白少妇内射xxxbb| 99在线人妻在线中文字幕| 高清在线视频一区二区三区 | 久久精品影院6| 级片在线观看| 91麻豆精品激情在线观看国产| 人妻久久中文字幕网| 久久人人精品亚洲av| 在线免费观看不下载黄p国产| 色哟哟哟哟哟哟| 一本久久精品| 日本黄色视频三级网站网址| 亚洲欧美清纯卡通| 美女脱内裤让男人舔精品视频 | 色5月婷婷丁香| 成年女人永久免费观看视频| 麻豆一二三区av精品| 亚洲熟妇中文字幕五十中出| 久久久久久久亚洲中文字幕| 久久久久久伊人网av| 亚洲美女搞黄在线观看| 99久久精品热视频| 尤物成人国产欧美一区二区三区| 国产伦精品一区二区三区视频9| 性欧美人与动物交配| av天堂在线播放| 简卡轻食公司| 日韩欧美 国产精品| 最近手机中文字幕大全| 精品熟女少妇av免费看| 久久午夜福利片| 亚洲色图av天堂| 精品免费久久久久久久清纯| 成人高潮视频无遮挡免费网站| 青春草亚洲视频在线观看| 人人妻人人看人人澡| 久久精品综合一区二区三区| 麻豆一二三区av精品| 变态另类丝袜制服| 国产在线男女| 狂野欧美激情性xxxx在线观看| 成人一区二区视频在线观看| 国产av不卡久久| 黄色欧美视频在线观看| 有码 亚洲区| 国产 一区 欧美 日韩| 国产精品综合久久久久久久免费| 成人午夜高清在线视频| 青春草视频在线免费观看| 亚洲av男天堂| 久久久久久久亚洲中文字幕| 国产av不卡久久| 联通29元200g的流量卡| 亚洲国产欧美人成| 97超碰精品成人国产| 噜噜噜噜噜久久久久久91| 小蜜桃在线观看免费完整版高清| 久久精品国产自在天天线| 精品人妻熟女av久视频| 欧美一区二区精品小视频在线| 麻豆成人午夜福利视频| 日韩精品青青久久久久久| 51国产日韩欧美| 日韩精品青青久久久久久| 非洲黑人性xxxx精品又粗又长| 欧美+亚洲+日韩+国产| 中文字幕精品亚洲无线码一区| 亚洲精品日韩在线中文字幕 | 亚洲欧美成人综合另类久久久 | 小说图片视频综合网站| 毛片女人毛片| 悠悠久久av| 成人av在线播放网站| 永久网站在线| 国产精品美女特级片免费视频播放器| 男女啪啪激烈高潮av片| 国产亚洲欧美98| 中文欧美无线码| 91精品国产九色| 黄片wwwwww| 人妻久久中文字幕网| 男人舔女人下体高潮全视频| 精品久久久久久久久久久久久| 国产伦精品一区二区三区视频9| 亚洲第一电影网av| 亚洲精品色激情综合| 中文字幕久久专区| 日韩强制内射视频| 精品少妇黑人巨大在线播放 | 国产精品久久电影中文字幕| 九九热线精品视视频播放| 亚洲欧美日韩高清在线视频| 男女边吃奶边做爰视频| av在线亚洲专区| 亚洲av第一区精品v没综合| 久久精品国产自在天天线| 精品人妻熟女av久视频| 国产欧美日韩精品一区二区| 在线国产一区二区在线| 黄色欧美视频在线观看| 亚洲av免费在线观看| 亚洲美女搞黄在线观看| 亚洲国产精品sss在线观看| 久久久久性生活片| 1000部很黄的大片| 成人无遮挡网站| www.av在线官网国产| 久久精品影院6| 国产亚洲av片在线观看秒播厂 | 久久午夜福利片| 九草在线视频观看| а√天堂www在线а√下载| 尾随美女入室| 黄色视频,在线免费观看| 欧美最黄视频在线播放免费| 久久久久久久久久黄片| 99久久人妻综合| 国产不卡一卡二| 嘟嘟电影网在线观看| 午夜福利高清视频| 亚洲aⅴ乱码一区二区在线播放| av在线亚洲专区| 国产视频首页在线观看| 啦啦啦韩国在线观看视频| 深夜精品福利| 国产久久久一区二区三区| 黄片wwwwww| 欧美色欧美亚洲另类二区| 久久精品人妻少妇| 简卡轻食公司| 悠悠久久av| 一本一本综合久久| 亚洲电影在线观看av| 能在线免费看毛片的网站| www.av在线官网国产| 久久综合国产亚洲精品| 18+在线观看网站| 我的女老师完整版在线观看| 久久精品国产亚洲av天美| 日韩中字成人| 男人狂女人下面高潮的视频| 久久欧美精品欧美久久欧美| 久久久久网色| 免费观看在线日韩| 少妇猛男粗大的猛烈进出视频 | 亚洲精品日韩av片在线观看| 人妻系列 视频| 天天躁夜夜躁狠狠久久av| 午夜福利视频1000在线观看| 久久久久久久久久久免费av| 一级毛片久久久久久久久女| 97超碰精品成人国产| 国产国拍精品亚洲av在线观看| 中文在线观看免费www的网站| 男人舔女人下体高潮全视频| 亚洲精品456在线播放app| 欧美最新免费一区二区三区| 69人妻影院| 久久久久久久久久成人| 亚洲一区高清亚洲精品| 日韩成人伦理影院| 亚洲精品粉嫩美女一区| 老女人水多毛片| 日韩欧美一区二区三区在线观看| 免费看av在线观看网站| 亚洲精品亚洲一区二区| 一本一本综合久久| 亚洲国产高清在线一区二区三| 六月丁香七月| 国产黄片美女视频| 两个人视频免费观看高清| 特级一级黄色大片| 午夜老司机福利剧场| 欧美一区二区精品小视频在线| 麻豆乱淫一区二区| 国内揄拍国产精品人妻在线| 精品久久久噜噜| 欧美性感艳星| 欧美区成人在线视频| 能在线免费观看的黄片| 少妇被粗大猛烈的视频| 啦啦啦观看免费观看视频高清| 一级二级三级毛片免费看| 国产精品一区二区在线观看99 | 色吧在线观看| 亚洲成av人片在线播放无| 色播亚洲综合网| 1000部很黄的大片| 欧美成人精品欧美一级黄| 日韩欧美三级三区| 国产伦一二天堂av在线观看| 久久久久九九精品影院| 只有这里有精品99| 国产在视频线在精品| 久久久久久大精品| 免费av不卡在线播放| av女优亚洲男人天堂| 亚洲国产精品sss在线观看| 观看免费一级毛片| 99九九线精品视频在线观看视频| 国产老妇女一区| 美女 人体艺术 gogo| 一区福利在线观看| 天堂av国产一区二区熟女人妻| 99国产精品一区二区蜜桃av| 国产高潮美女av| 女人被狂操c到高潮| 国产精品免费一区二区三区在线| 久久久久免费精品人妻一区二区| 亚洲va在线va天堂va国产| 欧美+亚洲+日韩+国产| 99久久人妻综合| av一本久久久久| .国产精品久久| 少妇人妻精品综合一区二区| 极品人妻少妇av视频| 国产午夜精品一二区理论片| 王馨瑶露胸无遮挡在线观看| 91午夜精品亚洲一区二区三区| 亚洲美女搞黄在线观看| 激情五月婷婷亚洲| 大陆偷拍与自拍| 久久精品国产亚洲网站| 少妇人妻久久综合中文| 91精品一卡2卡3卡4卡| 国模一区二区三区四区视频| 欧美日韩成人在线一区二区| 亚洲少妇的诱惑av| 丰满迷人的少妇在线观看| 国产熟女欧美一区二区| 超碰97精品在线观看| 日韩中文字幕视频在线看片| 国产在线视频一区二区| 国产精品久久久久久精品电影小说| av福利片在线| 午夜免费观看性视频| 毛片一级片免费看久久久久| 亚洲精品视频女| 国产在线一区二区三区精| 中文字幕最新亚洲高清| 高清午夜精品一区二区三区| 亚洲精品中文字幕在线视频| 又粗又硬又长又爽又黄的视频| 人妻制服诱惑在线中文字幕| 亚洲精品日本国产第一区| 少妇被粗大猛烈的视频| 80岁老熟妇乱子伦牲交| 丰满少妇做爰视频| 亚洲美女黄色视频免费看| 黄色欧美视频在线观看| 国产成人免费无遮挡视频| 日韩av不卡免费在线播放| 乱人伦中国视频| 青春草亚洲视频在线观看| 国产精品久久久久久久电影| 美女视频免费永久观看网站| 亚洲内射少妇av| 五月开心婷婷网| 亚洲精品第二区| 色吧在线观看| 亚洲精品美女久久av网站| 国产精品一区www在线观看| 亚洲怡红院男人天堂| 午夜免费男女啪啪视频观看| 午夜91福利影院| 黄片无遮挡物在线观看| 大又大粗又爽又黄少妇毛片口| 啦啦啦视频在线资源免费观看| 日韩一本色道免费dvd| av在线app专区| 欧美国产精品一级二级三级| 欧美人与善性xxx| 热re99久久国产66热| 80岁老熟妇乱子伦牲交| 在线观看免费日韩欧美大片 | 亚洲一级一片aⅴ在线观看| 欧美日韩亚洲高清精品| 女人精品久久久久毛片| 中文字幕久久专区| 青青草视频在线视频观看| 五月开心婷婷网| av不卡在线播放| 母亲3免费完整高清在线观看 | 色94色欧美一区二区| 中文字幕精品免费在线观看视频 | av在线观看视频网站免费| 亚洲国产欧美日韩在线播放| 精品久久久久久电影网| 日韩电影二区| 9色porny在线观看| 十八禁网站网址无遮挡| 性色avwww在线观看| 成年美女黄网站色视频大全免费 | 熟女av电影| 午夜福利视频精品| 久久久久久久久久人人人人人人| 少妇猛男粗大的猛烈进出视频| 岛国毛片在线播放| 国产精品久久久久成人av| 一级,二级,三级黄色视频| 纯流量卡能插随身wifi吗| xxx大片免费视频| 亚洲国产精品999| av卡一久久| 国产免费一级a男人的天堂| 日韩一区二区三区影片|