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

    遼南地區(qū)金州西團(tuán)瓢湖沼相沉積物粒度化學(xué)特征及環(huán)境意義

    2020-07-02 03:21:48張威成然馬瑞豐麻浩男葛潤澤孫波趙賀
    關(guān)鍵詞:粉砂沉積物平均值

    張威,成然,馬瑞豐,麻浩男,葛潤澤,孫波,趙賀

    遼寧師范大學(xué)地理科學(xué)學(xué)院,大連 116029

    遼南地區(qū)作為海陸交互地帶,受來自陸地和海洋的雙重氣候變化影響。該地區(qū)的湖泊沉積物序列連續(xù),是研究歷史時(shí)期氣候變化的理想載體。20世紀(jì)70年代,中國科學(xué)院貴陽地球化學(xué)研究所對(duì)遼南地區(qū)普蘭店湖相沉積進(jìn)行研究,通過孢粉分析與14C測年討論了全新世以來的氣候變化,早全新世(10.0~8.0 kaBP)氣候總體溫和干燥;中全新世(8.0~2.5 kaBP)前期氣候溫暖濕潤,后期溫暖較干燥;晚全新世(2.5 kaBP至今)氣溫偏涼,前期濕潤,后期干燥[1]。20世紀(jì)90年代,李雪銘對(duì)遼南大蓮泡湖相沉積進(jìn)行研究,通過孢粉與粒度分析恢復(fù)了區(qū)內(nèi)全新世以來的氣候演化,所得到的結(jié)論與中國科學(xué)院貴陽地化所基本一致[2]。雖然不同的學(xué)者對(duì)研究區(qū)全新世的氣候變化特征進(jìn)行了初步研究,但是時(shí)間尺度比較短。對(duì)于晚更新世以來的氣候環(huán)境變化研究成果較少,且絕大多數(shù)根據(jù)陸架黃土,采用一些比較傳統(tǒng)的環(huán)境替代性指標(biāo)對(duì)環(huán)境變化進(jìn)行研究[3-4],分辨率并不高,尤其是對(duì)粒度與化學(xué)數(shù)據(jù)的處理方面,分析方法和手段比較單一。因此,本文在詳細(xì)野外調(diào)查的基礎(chǔ)上,對(duì)遼南地區(qū)保存完整的西團(tuán)瓢湖相沉積剖面進(jìn)行高分辨率采樣,采用端元模型對(duì)湖相沉積物的粒度變化特征進(jìn)行深入分析,并結(jié)合粒度參數(shù)以及地球化學(xué)元素CIA、燒失量、電導(dǎo)率等指標(biāo),探討研究區(qū)晚更新世末期以來的環(huán)境演化特征。

    1 研究區(qū)概況

    研究區(qū)位于遼東半島南部,地貌表現(xiàn)為中央高、東西兩側(cè)低(圖1)。區(qū)內(nèi)由泡子崖向斜、付家村-七頂山背斜和西團(tuán)瓢向斜組成一個(gè)北西向復(fù)向斜構(gòu)造。研究區(qū)出露地層有震旦系灰?guī)r、泥質(zhì)灰?guī)r、石英巖,古生界寒武系云母砂巖、頁巖、灰?guī)r,奧陶系硅質(zhì)灰?guī)r、白云質(zhì)灰?guī)r。本區(qū)屬于暖溫帶亞濕潤氣候區(qū),西臨渤海,受海洋性氣候影響大,年平均氣溫10.1 ℃,年降水量 660 mm[5]。西團(tuán)瓢剖面位于遼寧省大連市金州區(qū)西團(tuán)瓢村,處于由大西山、七頂山及大東山圍成的盆地中(圖2)。盆地東側(cè)為三十里堡河,該河從東南流向西北,穿盆而過,注入渤海灣。盆地南部剝蝕發(fā)育了多條溝谷,出山口后均匯入湖盆[2]。剖面距離渤海約10.36 km,距離黃海約36.29 km。大蓮泡位于西團(tuán)瓢西南部2 km處,同處于山間湖盆中,發(fā)育了約3 m厚的松散沉積物。下層(1.65~3 m),為黃棕色亞黏土層,以坡面溝谷沖積為主,大蓮泡湖沼處于初始階段,孢粉數(shù)量較少,且多為耐寒冷干旱的花粉,說明當(dāng)時(shí)氣候惡劣,植被稀疏,雨水沖刷坡面攜帶碎屑物進(jìn)入湖盆;中層(0.85~1.65 m),為黑褐色黏土層,富含有機(jī)質(zhì),早期坡面溝谷沖積作用減少,以湖泊沉積為主,環(huán)境穩(wěn)定,沉積了大量細(xì)顆粒物質(zhì),有大量闊葉樹種及水生植物的花粉,推斷此階段處于全新世溫暖濕潤期,但有一段出現(xiàn)大量鹽生蒿屬的孢粉,與海水侵入湖盆有關(guān);上層(地表至0.85 m處)為黃棕色粉砂層,下段為草炭層,14C年代測定為667±100a,以河流作用為主,草本植物花粉增加,水生植物減少,標(biāo)志氣候又逐漸惡劣,湖盆變淺,最終干涸消失[2]。

    2 資料和方法

    2.1 剖面特征

    西團(tuán)瓢剖面(39°14′54.32″N、121°43′39.61″E)位于遼寧省大連市金州區(qū)西團(tuán)瓢村(圖3)。剖面厚354 cm,根據(jù)巖性特征將剖面劃為8個(gè)地層單元,剖面巖性特征描述如下:

    (1)0~24 cm,現(xiàn)代生草層,富含植物根系;

    (2)24~74 cm,粉—中砂層,由深棕色過渡至淺棕色,夾有灰白色黏土條帶,較疏松;

    (3)74~136 cm,泥炭層,由棕黑色過度到黑色,緊實(shí),富含植物殘?bào)w;

    圖1 研究區(qū)和采樣點(diǎn)位置圖Fig.1 Sampling site location in study area

    圖2 研究區(qū)地形圖Fig.2 Topographic map of study area

    圖3 西團(tuán)瓢剖面照片和巖性柱狀圖Fig.3 Photographs and lithologic histogram of the Xi Tuanpiao profile

    (4)136~176 cm,粉—中砂層,深灰色,夾淺灰色黏土條帶,較緊實(shí);

    (5)176~222 cm,粉—細(xì)砂層,淺灰色,212 cm以上可見螺殼;

    (6)222~232 cm,中砂層,灰白色,有黃色銹斑,較疏松;

    (7)232~276 cm,粉—中砂層,棕灰色,夾雜黃色銹斑及黏土條,偶有白色鈣結(jié)核,較緊實(shí);

    (8)276~354 cm,粉砂—黏土層,棕黃色,有黃色銹斑,夾有灰色黏土條;下部為灰色大石塊。

    在剖面每個(gè)層位進(jìn)行樣品采集,每隔2 cm采集一個(gè)樣品,共采集樣品177個(gè)。

    2.2 實(shí)驗(yàn)方法

    粒度測試使用LS13320激光衍射粒度儀,粒級(jí)以Krumbein粒度(φ)表示,轉(zhuǎn)換公式為 φ=?log2d,其中d為粒級(jí)(mm)。

    化學(xué)元素測試使用日本理學(xué)公司生產(chǎn)的ZSX PrinmusⅡ型X射線熒光光譜儀。測試了樣品中元素所對(duì)應(yīng)的氧化物的含量。

    燒失量測量:先將樣品在110 ℃下烘干1 h以去除水分,再用精度為0.000 1 g的電子秤稱取1 g樣品,使用550 ℃馬沸爐燃燒4 h,冷卻后稱重。

    黏土混濁水電導(dǎo)率測量:先將樣品在110 ℃下烘干 48 h,再將樣品研磨至 31 目以下,稱取 10 g 放入燒杯中,加入 120 mL 蒸餾水。攪拌 3 min,靜置 1 h后開始測量。儀器為DDS-320型電導(dǎo)率儀。以上實(shí)驗(yàn)均在遼寧師范大學(xué)實(shí)驗(yàn)中心完成。

    14C年代樣品為96 cm深處的植物殘?bào)w和186 cm深處的螺殼,利用加速器質(zhì)譜儀測定,在美國Beta實(shí)驗(yàn)室完成。測定的年代分別為 1 170±30 cal.aBP、11 510±40 cal.aBP。0~96 cm 的沉積速率為 8.20 mm/a,96~186 cm 的沉積速率為 0.87 mm/a。底部 354 cm處的年代推測為30 808 cal.aBP。其他各層位年代如圖3所示。

    2.3 分析方法

    2.3.1 粒度特征

    沉積物粒度特征通過粒度組成特征、頻率曲線特征、粒度統(tǒng)計(jì)參數(shù)特征來表達(dá),其中統(tǒng)計(jì)參數(shù)包括平均粒徑、中值粒徑、標(biāo)準(zhǔn)偏差、偏態(tài)和峰態(tài),以上參數(shù)計(jì)算采用圖解法[6],對(duì)各參數(shù)的分級(jí)以及定性描述術(shù)語參照盧連戰(zhàn)的分級(jí)標(biāo)準(zhǔn)[7]。

    2.3.2 端元分析

    沉積物是由多種沉積動(dòng)力因素共同作用的結(jié)果,因此將不同沉積動(dòng)力過程分解出來可以更有效地反演物質(zhì)來源及沉積過程。端元分析模型主要用于對(duì)沉積物粒度數(shù)據(jù)進(jìn)行組分分離[8],在識(shí)別和分析沉積物各動(dòng)力組分中得到廣泛應(yīng)用[9-10]。Weltje結(jié)合前人的研究,歸納出了端元分析模型,并提出了動(dòng)力組分(Dynamic Population)這一概念。端元分析模型的目的就在于從沉積物的粒度數(shù)據(jù)中提取出與動(dòng)力組分相對(duì)應(yīng)的端元(End Member,簡稱EM)。

    湖泊沉積物一般由多個(gè)不同時(shí)期的動(dòng)力組分混合而成,根據(jù)端元分析模型,粒度數(shù)據(jù)集X可以用多個(gè)端元成分的數(shù)據(jù)集B的線性組合來表示:

    其中,X(m*n)為沉積物樣品粒度的矩陣,M(m*k)為端元含量的混合矩陣,B(k*n)為端元粒度的矩陣(m:樣品總數(shù),n:粒級(jí)總數(shù),k:端元總數(shù))

    由于沉積物搬運(yùn)過程的特殊性以及在采集樣品和進(jìn)行實(shí)驗(yàn)時(shí)可能出現(xiàn)的誤差,一般很難得到符合公式(1)的完美混合矩陣,因此矩陣由兩部分組成,即完美混合矩陣X’和誤差矩陣E[11]:

    本文采用Peterson改進(jìn)后的端元分析模型[12],在Matlab軟件中運(yùn)行端元計(jì)算程序,導(dǎo)入粒級(jí)數(shù)據(jù),并假定端元數(shù)在1~10之間。端元數(shù)越大,粒級(jí)復(fù)相關(guān)系數(shù)R2越接近于1,R2在0.8以上表示端元數(shù)基本滿足擬合的要求。為了便于分離出沉積物的動(dòng)力組分,端元分析盡量以端元數(shù)為少的原則進(jìn)行,因此本文最終選擇3個(gè)端元對(duì)該組粒度數(shù)據(jù)進(jìn)行分析(圖4)。

    3 結(jié)果分析

    3.1 粒度組成特征

    按照吳正的粒度分級(jí)標(biāo)準(zhǔn)[13],將沉積物粒徑分為 6 個(gè)等級(jí):黏土:0~4 μm,粉砂:4~63 μm,極細(xì)砂:63~125 μm,細(xì)砂:125~250 μm,中砂:250~500 μm,粗砂:>500 μm。按照這個(gè)標(biāo)準(zhǔn),對(duì)西團(tuán)瓢湖沼相沉積物粒度組成特征進(jìn)行分析,結(jié)果如圖5所示。

    從粒度組成圖中可以看出,整個(gè)剖面粗砂含量最少,平均值僅為4.4%,極細(xì)砂含量次之,粉砂含量最高,平均值達(dá)到43.4%。不同層位沉積物的粒級(jí)組成存在一定差異。0~24 cm的現(xiàn)代生草層以中砂和細(xì)砂含量為主,平均含量分別為34.98%和28.24%;24~74 cm沉積物以中砂和粉砂為主,平均含量分別為29.97%和25.48%,該層沉積物粒度組分隨深度變化的幅度較大。與上一層相比,74~136 cm的泥炭層的粉砂和黏土明顯增加,粉砂含量由25.48%增加到53.77%,黏土含量由5.51%增加到14.26%。136~176 cm的沉積物以粉砂和中砂為主,平均含量分別為46.56%和15.33%。黏土與粗砂含量比泥炭層偏低。176~222 cm沉積物以粉砂和極細(xì)砂為主,平均含量分別為59.54%和15.75%,其中粉砂的平均含量為整個(gè)剖面的最大值。222~232 cm的沉積物以中砂和細(xì)砂為主,平均含量分別為35.30%和24.70%,黏土含量僅為0.50%,是整個(gè)剖面的最低值。232~276 cm的沉積物與上一層位相比顆粒較細(xì),以粉砂和中砂為主,平均含量分別為34.42%和18.30%。276~354 cm的沉積物以粉砂和黏土為主,平均含量分別為46.73%和23.00%,其中黏土組分的平均含量比其他層位的黏土含量都高,粗砂含量極少,僅為0.61%。

    圖4 端元線性相關(guān)(a)和角度離差(b)Fig.4 Linear correlation of end members (a) and angular dispersion (b)

    圖5 西團(tuán)瓢剖面粒度姐成Fig.5 The composition of the grain size of the Xi Tuanpiao profile

    3.2 頻率曲線特征

    頻率曲線的橫坐標(biāo)一般用沉積物的粒徑區(qū)間來表示,而縱坐標(biāo)用不同粒徑區(qū)間的體積百分比表示。頻率曲線呈現(xiàn)的粒度特征可以直接反映當(dāng)時(shí)的沉積環(huán)境。峰值個(gè)數(shù)越多,代表著沉積環(huán)境越復(fù)雜,曲線只有一個(gè)峰意味著沉積環(huán)境比較穩(wěn)定,水動(dòng)力條件單一。曲線出現(xiàn)兩個(gè)及以上峰表示沉積環(huán)境不穩(wěn)定。沉積時(shí)期的水動(dòng)力條件會(huì)關(guān)系到峰值分布的區(qū)間,如果峰值處于較粗的粒度區(qū)間,則表示水動(dòng)力條件較強(qiáng),如果峰值處于較細(xì)的粒度區(qū)間,則表示水動(dòng)力條件較弱。

    如圖6所示,剖面沉積物的粒度頻率分布曲線整體以單峰和雙峰為主,峰值主要出現(xiàn)在2Ф附近,另有一部分雙峰曲線的峰值出現(xiàn)在2Ф、5Ф和7Ф附近,可見沉積環(huán)境不穩(wěn)定。由圖5所示,0~24、24~74、222~232、232~276 cm 的沉積物粒度頻率分布曲線以單峰為主,且峰值出現(xiàn)在2Ф附近,說明在這些階段湖泊沉積環(huán)境穩(wěn)定,且水動(dòng)力較強(qiáng),其中24~74和232~276 cm的頻率曲線出現(xiàn)弱峰,說明在這兩個(gè)階段沉積環(huán)境出現(xiàn)一定波動(dòng)。176~222 cm的頻率曲線呈現(xiàn)單雙峰交雜的情況,峰值主要集中在5Ф附近,還有部分雙峰峰值集中在2Ф和7Ф附近,說明此時(shí)期沉積環(huán)境不穩(wěn)定,且水動(dòng)力較弱,呈現(xiàn)低能。74~136、136~176和276~354 cm的頻率曲線以雙峰為主,且交錯(cuò)繁復(fù),峰值大致出現(xiàn)在2Ф和5~7Ф附近,說明這兩個(gè)階段沉積環(huán)境復(fù)雜,水動(dòng)力強(qiáng)弱交替。

    圖6 西團(tuán)瓢剖面粒度頻率分布曲線Fig.6 Drain-size distribution of the Xi Tuanpiao profile

    3.3 粒度統(tǒng)計(jì)參數(shù)特征

    沉積物粒度特征除了通過頻率分布曲線表達(dá),還可用平均粒徑、標(biāo)準(zhǔn)偏差、偏態(tài)和峰態(tài)等參數(shù)進(jìn)行分析。不同參數(shù)表示的沉積意義各不相同,同時(shí)各參數(shù)之間相互關(guān)聯(lián),因此通過對(duì)粒度參數(shù)進(jìn)行分析可以更好地推斷沉積環(huán)境。

    如圖7所示,金州西團(tuán)瓢剖面的平均粒徑為0.70~7.30Ф,且波動(dòng)明顯?,F(xiàn)代生草層(0~24 cm)平均粒徑最大,平均值為2.53Ф;中砂層(222~232 cm)次之,平均值為 2.81Ф;粉—中砂層(24~74 cm)居于其后,平均值為3.14Ф;之后平均粒徑由大到小的層位依次是粉—中砂層(232~276 cm)、粉—中砂層(136~176 cm)、粉—細(xì)砂層(176~222 cm)、泥炭層(74~136 cm),平均粒徑的平均值分別為:4.32、4.43、4.86、4.99Ф。粉砂—黏土層(276~354 cm)最小,平均值為 5.88Ф。

    剖面沉積物粒度的分選系數(shù)范圍為1.15~3.11,數(shù)值越大,代表分選效果越差。根據(jù)圖解法對(duì)分選系數(shù)等級(jí)的定性描述可知,該剖面整體處于分選較差到分選差的程度,但是不同層位的分選程度不同,可進(jìn)行剖面內(nèi)部比較。現(xiàn)代生草層(0~24 cm)的分選系數(shù)平均值為1.89,與其他層位相比,分選效果最好;粉—細(xì)砂層(176~222 cm)次之,分選系數(shù)平均值為 1.95,;中砂層(222~232 cm)居于其后,分選系數(shù)平均值為2.01;之后分選效果由好變差的層位依次是粉—中砂層(24~74 cm)、粉—中砂層(136~176 cm)、粉砂—黏土層(276~354 cm)、泥炭層(74~136 cm),分選系數(shù)依次是 2.05、2.46、2.48、2.54;粉—中砂層(232~276 cm)分選效果最差,分選系數(shù)平均值為2.79。

    剖面整體偏態(tài)值為-0.29~0.74,平均值為0.13,剖面整體表現(xiàn)為正偏態(tài),沉積物顆粒較粗,分選較差,沉積時(shí)期水動(dòng)力較強(qiáng)且不穩(wěn)定;不同層位偏態(tài)值不同,可進(jìn)行內(nèi)部比較。中砂層(222~232 cm)偏態(tài)值最大,現(xiàn)代生草層(0~24 cm)次之,粉—中砂層(24~74 cm)位于其后,平均值分別為 0.49、0.48、0.39,這3個(gè)層位的偏態(tài)特征都定性描述為極正偏。粉—中砂層(232~276 cm)、粉—細(xì)砂層(176~222 cm)的偏態(tài)值的平均值分別為 0.26、0.1,這兩個(gè)層位的偏態(tài)特征描述為正偏。泥炭層(74~136 cm)、粉—中砂層(136~176 cm)、粉砂—黏土層(276~354 cm)的偏態(tài)值平均值分別為-0.07、0.09、-0.08,這3個(gè)層位的偏態(tài)特征描述為近對(duì)稱。

    整個(gè)剖面的峰態(tài)值為0.63~2.15,平均值為1.06,西團(tuán)瓢剖面的峰態(tài)多表現(xiàn)為中等寬度,說明沉積物粒徑較為集中,沉積時(shí)期的水動(dòng)力條件總體表現(xiàn)為中能。粉—中砂層(136~176 cm)峰態(tài)值最小,其平均值為0.82,定性描述為平坦;泥炭層(74~136 cm)、粉砂—黏土層(276~354 cm)、粉—細(xì)砂層(176~222 cm)峰 態(tài)值的平 均值分別 為 0.91、0.99、1.1,這3個(gè)層位的峰態(tài)特征為中等寬度;粉—中砂層(24~74 cm)峰態(tài)值的平均值為 1.23,峰態(tài)特征描述為尖銳;中砂層(222~232 cm)、現(xiàn)代生草層(0~24 cm)、粉—中砂層(232~276 cm)峰態(tài)值的平均值分別為1.65、1.74、1.85,峰態(tài)特征描述為很尖銳。

    圖7 西團(tuán)瓢剖面粒度特征曲線對(duì)比Fig.7 Comparison diagram of profile particle size characteristic curve of the Xi Tuanpiao profile

    3.4 端元分析結(jié)果

    通過端元分析分離出的3個(gè)端元的粒級(jí)頻率分布曲線圖可以看出,各端元呈現(xiàn)的不同粒級(jí)的頻率分布曲線均為單峰。根據(jù)各自的曲線形態(tài)特征,結(jié)合各深度地層中各端元組分的含量(圖8)以及其他氣候指標(biāo)綜合判斷各端元所代表的沉積動(dòng)力特征。從端元1到端元3主體粒級(jí)粒徑由小變大,其中端元1的平均粒徑為 7.49 μm,眾數(shù)為 9 μm左右,屬于粉砂,分選性最差,峰態(tài)最平坦,可能代表風(fēng)成沉積和流水?dāng)y帶的入湖細(xì)顆粒;端元2的平均粒徑為 38.02 μm,眾數(shù)為 40 μm 左右,屬于粉砂,分選水平最好,接近正態(tài)分布,與端元1相比峰態(tài)略尖銳,可能代表弱水流沉積;端元3的平均粒徑最大,為195.23 μm,眾數(shù)為 300 μm 左右,屬于中砂,分選水平居中,峰度值最大(表1),可能屬于沖洪積相沉積,水動(dòng)力較強(qiáng)。從端元組分的含量上可以看出,端元1與端元2、3大致呈相反的趨勢,隨著地層的加深,端元1的含量呈現(xiàn)波動(dòng)上升的趨勢,端元2與端元3組分含量隨著地層的加深都呈現(xiàn)波動(dòng)下降的趨勢。由圖9可以看出,端元1與化學(xué)蝕變指數(shù)(CIA)較為對(duì)應(yīng)?;瘜W(xué)蝕變指數(shù)(CIA)是利用常量元素判斷化學(xué)風(fēng)化程度的指標(biāo),可以有效指示長石風(fēng)化成黏土礦物的程度,CIA值越大意味著風(fēng)化程度越強(qiáng),CIA值在50~65之間表示干冷的氣候條件導(dǎo)致的低等化學(xué)風(fēng)化程度;CIA值在65~85之間表示暖濕的氣候條件導(dǎo)致的中等化學(xué)風(fēng)化程度。金州西團(tuán)瓢剖面沉積物CIA值的平均值為63.31,屬于中等偏低的化學(xué)風(fēng)化程度[14-16]。

    圖8 西團(tuán)瓢剖面各端元頻率分布曲線Fig.8 Frequency distribution curves of each end member in the Xi Tuanpiao profile

    表1 西團(tuán)瓢剖面各端元粒度參數(shù)特征Table 1 The end member parameters characteristic of The Xi Tuanpiao profile

    燒失量表示在高溫條件下,樣品損失的部分占總質(zhì)量的百分比。燒失量值越大表示樣品中有機(jī)質(zhì)含量越高。湖泊沉積物中的有機(jī)質(zhì)含量可以有效地指示湖區(qū)氣候環(huán)境特征,濕潤的環(huán)境是有機(jī)質(zhì)富集的重要條件,高溫和寒冷都有助于有機(jī)質(zhì)的積累[17]。金州西團(tuán)瓢剖面有機(jī)質(zhì)含量的平均值為13.61%,剖面內(nèi)部各層位之間又存在差別,可進(jìn)一步對(duì)比。

    端元1是細(xì)顆粒組分,在現(xiàn)代生草層(0~24 cm)及向泥炭過渡的地層(24~74 cm)中含量較少,在泥炭層(74~136 cm)含量較多,在含有螺殼的湖泊沉積層(136~222 cm)中逐漸減少至最低。但隨著地層的加深(222~354 cm),在 3.08~1.68 kaBP 端元1的含量逐漸增加至最高,此階段氣溫低,降水少,風(fēng)力強(qiáng)勁,由東亞季風(fēng)帶來的西北內(nèi)陸地區(qū)的細(xì)顆粒組分較多。同時(shí)少量大氣降水形成動(dòng)力較弱的坡面徑流,攜帶山區(qū)細(xì)顆粒進(jìn)入湖泊,因此端元1可能指示風(fēng)成沉積及坡面流水?dāng)y帶的細(xì)顆粒沉積,由于不同礦物抗風(fēng)化能力的差異性,使得沉積物總體上粒徑變小,差異性變大,分選性變大;端元2可能指示弱水流沉積,沉積物粒徑較粗,分選較好。弱水流沉積指代由坡面徑流與溝谷水流攜帶流入湖泊的偏粗顆粒沉積物。隨著地層的加深(0~232 cm),端元 2 含量逐漸增加,在 186~232 cm(1.68~1.15 kaBP)富集。此階段為氣候由干冷向暖濕的過渡時(shí)期,降水量逐漸增大,西團(tuán)瓢剖面屬于山前平原,位于七頂山東南部,歷史時(shí)期由于坡面流水作用,接收來自山區(qū)的偏粗顆粒沉積物;端元3可能指示水動(dòng)力較強(qiáng)的沖洪積相沉積,分選較差。

    3.5 剖面沉積環(huán)境

    遼南大蓮泡與西團(tuán)瓢同處于由大西山、七頂山、大東山組成的山間小盆地,三十里堡河居于盆地東側(cè),由東南向西北穿盆而過,注入渤海灣。通過孢粉分析發(fā)現(xiàn),大蓮泡1.10~1.70 m的湖沼相沉積物中的淡水水生植物孢粉被鹽生蒿屬植物孢粉代替,但該段沉積物中還存在淡水螺殼,說明此層沉積物形成時(shí)可能受到海水的影響,但是沒有大規(guī)模海水入侵。從區(qū)域地理環(huán)境看,研究區(qū)三面環(huán)山,有入海河流流經(jīng)此地,在全新世高海面時(shí)期,海水漲潮時(shí)可能溯古河道東進(jìn)影響研究區(qū)淡水湖體的環(huán)境[2]。全新世以來,遼東半島的海陸分布有過明顯的變遷,在約 10 000~5 000 aBP 前是海侵時(shí)期,北黃海海岸線向陸地延伸 5~15 km,5 000 aBP以來是海退時(shí)期,海岸線呈波浪式地向海洋方向移動(dòng)[1、18]。

    黏土混濁水電導(dǎo)率值(以下簡稱EC)與溶液中的硫酸根離子及硫酸鹽的含量成正比[19],沉積環(huán)境從海相過渡到陸相,沉積物中硫化物的含量會(huì)逐漸減少。日本學(xué)者將黏土混濁水電導(dǎo)率與硅藻組合進(jìn)行對(duì)比,總結(jié)出陸相沉積物的EC值小于0.6 mS/cm,海陸過渡相沉積物的EC值為0.6~1.5 mS/cm,海相沉積物的EC值大于1.5 mS/cm[20]。在海岸帶沉積物形成初期,水體中的硫酸根離子由細(xì)菌活動(dòng)被還原成硫化氫,并與鐵結(jié)合形成黃鐵礦。海岸線退縮,黃鐵礦依然被作為穩(wěn)定的物質(zhì)保存在沉積物中[21-23]。近年來,多位學(xué)者通過黏土混濁水電導(dǎo)率劃分遼東半島的大洋河平原、八岔溝平原、大清河平原以及下遼河平原的海陸相沉積[24-28]。金州西團(tuán)瓢沉積物測得的電導(dǎo)率值最大為0.451 mS/cm,遠(yuǎn)小于0.6 mS/cm(見圖9)。結(jié)合前人的研究,在西團(tuán)瓢剖面沉積時(shí)期,海水曾由三十里堡河上溯,但影響范圍與持續(xù)時(shí)間都不甚長遠(yuǎn),因此該剖面為可能短暫受海水影響的淡水湖沉積物。

    圖9 西團(tuán)飄剖面粒度端元組分、化學(xué)蝕變指數(shù)、燒失量及電導(dǎo)率的垂向變化特征Fig.9 The change of each end member content,CIA,loss on nlgnition and electric conductivity in depth of the Xi Tuanpiao profile

    4 討論

    4.1 沉積物粒度與化學(xué)指標(biāo)縱向變化及其環(huán)境意義

    0~24 cm(0.29 kaBP 至今),顆粒物相對(duì)較粗,以中砂和細(xì)砂含量為主,平均含量分別為34.98%和28.24%;粗砂的平均含量為9.28%。平均粒徑范圍為2.13~2.91Ф,平均值為2.53Ф。分選系數(shù)為1.13~2.49,平均值為1.89,在整個(gè)剖面中屬于分選較好的一層。偏態(tài)值為0.38~0.57,平均值約為0.48,為極正偏,說明沉積物以粗粒級(jí)為主。峰態(tài)值的平均值為1.74,很尖銳。粒度頻率分布曲線以單峰為主,而且峰值出現(xiàn)在2Ф附近??偟膩碚f,此階段沉積環(huán)境穩(wěn)定。CIA指數(shù)為64.32%,指示中等偏低的化學(xué)風(fēng)化程度。燒失量為14.27%,高于剖面整體的平均值,由于地表生長著大片蘆葦,根系很深,使得該層有機(jī)質(zhì)含量偏高。此階段氣候偏涼偏干。

    24~74 cm(0.29~0.90 kaBP)沉積物以中砂和粉砂為主,平均含量分別為29.97%和25.48%,該層沉積物粒度組分隨深度變化的幅度較大。平均粒徑為1.85~4.42Ф,平均值為3.14Ф。分選系數(shù)為1.15~2.70,平均值為2.05,分選效果屬于中等水平。偏態(tài)值為-0.02~0.60,平均值約為0.39,極正偏。表面沉積物大多為粗顆粒,分選水平較差。峰態(tài)值為0.70~2.15,平均值為1.23,尖銳。粒度頻率分布曲線以單峰為主,且峰值出現(xiàn)在2Ф附近,但部分樣品出現(xiàn)弱峰,說明此階段沉積環(huán)境不穩(wěn)定。CIA指數(shù)為66.43%,表示中等水平的化學(xué)風(fēng)化程度;燒失量為12.53%,略低于剖面整體平均值,說明該時(shí)期氣溫偏涼,降水偏少。

    與上一層相比,74~136 cm(0.90~5.77 kaBP)的泥炭層的粉砂和黏土明顯增加,粉砂含量由25.48%增加到53.77%,黏土含量由5.51%增加到14.26%。分選系數(shù)為1.80~3.11,平均值為2.54,分選水平較差。峰態(tài)值的平均值為0.91,中等寬度(正態(tài))。頻率曲線以雙峰為主,且交錯(cuò)繁復(fù),峰值大致出現(xiàn)在2Ф和5~7Ф附近,總體來看,這個(gè)階段沉積環(huán)境復(fù)雜,水動(dòng)力強(qiáng)弱交替。CIA指數(shù)為71.67%,表示中等偏高的化學(xué)風(fēng)化程度;燒失量為17.37%,該層是泥炭層,由植物腐化形成,說明沉積時(shí)期的氣候溫暖濕潤,適宜于植物的生長。

    136~176 cm(5.77~10.36 kaBP)的沉積物大多為粉砂和中砂,平均含量分別為46.56%和15.33%。黏土與粗砂含量比泥炭層偏低。平均粒徑范圍為3.27~5.03Ф,平均值為4.43Ф。分選系數(shù)為1.97~2.76,平均值為2.46,分選水平較差。峰態(tài)值為0.70~1.27,平均值為0.82,峰態(tài)平坦。與上層相似,頻率曲線以雙峰為主,且交錯(cuò)繁復(fù),峰值大致出現(xiàn)在2Ф和5~7Ф附近,說明此階段沉積環(huán)境復(fù)雜,水動(dòng)力強(qiáng)弱交替。CIA指數(shù)最大為72.26%,表示中等偏高的化學(xué)風(fēng)化程度,燒失量為5.74%,說明此階段氣候溫暖但降水量偏低。

    176~222 cm(10.36~15.64 kaBP)沉積物多為粉砂和極細(xì)砂,平均含量分別為59.54%和15.75%,其中粉砂的平均含量為整個(gè)剖面的最大值。平均粒徑為3.99~5.69Ф,平均值為4.86Ф。分選系數(shù)為1.56~2.56,平均值為1.95。分選效果要比上兩層好。偏態(tài)值為-0.14~0.29,平均值約為0.10,為正偏,分選程度略差。峰態(tài)值為0.77~1.31,平均值為1.10,峰態(tài)為中等(正態(tài))。頻率曲線呈現(xiàn)單雙峰交雜的情況,峰值主要集中在5Ф附近,還有部分雙峰峰值集中在2和7Ф附近??偟膩砜?,此時(shí)期沉積環(huán)境不穩(wěn)定。CIA指數(shù)最小為37.37%,屬于低等化學(xué)風(fēng)化程度;燒失量為21.60%,這是因?yàn)樵搶映练e了大量螺殼,導(dǎo)致該層有機(jī)質(zhì)含量高。由此說明此階段氣候寒冷干旱。

    222~232 cm(15.64~16.79 kaBP)的沉積物多為中砂和細(xì)砂,平均含量分別為35.30%和24.70%,黏土含量僅為0.50%,是整個(gè)剖面的最低值。平均粒徑為2.12~4.32Ф。分選系數(shù)平均值為2.01。偏態(tài)值為0.31~0.58,平均值約為0.49,極正偏,說明沉積物以粗顆粒為主,分選水平較差。峰態(tài)值為0.83~2.02,平均值為1.65,峰態(tài)很尖銳。沉積物粒度頻率分布曲線以單峰為主,且峰值出現(xiàn)在2Ф附近。綜上所述,該階段湖泊沉積環(huán)境穩(wěn)定。CIA指數(shù)的平均值為68.85%,指示中等化學(xué)風(fēng)化水平;燒失量為2.22%,該層有機(jī)質(zhì)含量為整個(gè)剖面的最低值。說明此階段氣候條件不適宜動(dòng)植物的生長,呈現(xiàn)寒冷干旱的狀態(tài)。

    232~276 cm(16.79~21.85 kaBP)的沉積物與上一層位相比顆粒較細(xì),以粉砂和中砂為主,平均含量分別為34.42%和18.30%。平均粒徑為0.72~5.72Ф,平均值為4.32Ф。分選系數(shù)為1.84~3.06,平均值為2.79,此階段分選系數(shù)跨度大,且平均值比其他層位都大,分選效果差。偏態(tài)值為-0.12~0.74,平均值為0.26,正偏,與上層相比,沉積物顆??傮w偏細(xì),分選水平較差。峰態(tài)值為0.69~1.96,平均值為1.85,峰態(tài)很尖銳。沉積物粒度頻率分布曲線以單峰為主,但少量樣品出現(xiàn)弱峰,說明在這個(gè)階段沉積環(huán)境不穩(wěn)定,出現(xiàn)波動(dòng)。CIA指數(shù)為65.95%,表示中等化學(xué)風(fēng)化程度;燒失量為12.87%,低于剖面整體平均值,說明此階段氣候條件不適宜動(dòng)植物的生長,呈現(xiàn)寒冷干旱的狀態(tài)。

    276~354 cm(21.85~30.81 kaBP)的沉積物仍以細(xì)顆粒(粉砂和黏土)組分為主,平均含量分別為46.73%和23.00%,其中黏土組分的平均含量比其他層位都高,粗砂含量極少,僅為0.61%。平均粒徑為4.60~7.30Ф,平均值為5.88Ф。分選系數(shù)為1.75~3.03,平均值為2.48。偏態(tài)值為-0.27~0.22,平均值約為-0.08,近對(duì)稱,與上層相比,分選水平較好。峰態(tài)值的平均值為0.99,峰態(tài)中等(正態(tài))。頻率曲線以雙峰為主,且交錯(cuò)繁復(fù),峰值大致出現(xiàn)在2和5~7Ф附近。綜上說明此階段沉積環(huán)境復(fù)雜,水動(dòng)力強(qiáng)弱交替。CIA指數(shù)64.68%,指示中等偏低的化學(xué)風(fēng)化水平;燒失量為12.31%,數(shù)值偏低,表示氣候條件不適于生物生存,呈現(xiàn)寒冷干旱的狀態(tài)。

    4.2 三個(gè)端元所指示的環(huán)境意義

    端元1組分在地層中的變化趨勢與CIA指數(shù)的變化趨勢相近,說明端元1的組分特征可能代表了東亞季風(fēng)的動(dòng)力強(qiáng)度,由季風(fēng)帶來的溫濕變化導(dǎo)致沉積物風(fēng)化程度的變化,使得礦物顆粒總體變小。此外,少量大氣降水形成動(dòng)力較弱的坡面徑流,攜帶山區(qū)細(xì)顆粒進(jìn)入湖泊。276~354 cm沉積層端元1含量最大,平均值為63.00%,CIA指數(shù)的平均值為64.68%,說明這個(gè)階段氣候較干旱。

    端元2組分在各地層中所占的比重偏少,整體平均值為29.99%。該組分的平均粒徑較端元1偏大,為38.02 μm,眾數(shù)值也在此附近,分選效果最好,端元頻率分布曲線接近正態(tài)分布。176~222 cm沉積層中端元2組分含量最大,平均值為64.09%,CIA值最小,平均值為37.37%。該層沉積物顆粒較細(xì),以粉砂為主,并有螺殼出現(xiàn)。結(jié)合上文分析可知,該層水動(dòng)力較弱,呈現(xiàn)低能。由此說明,端元2可能指示弱水流沉積,以此組分為主導(dǎo)的沉積時(shí)期氣候偏寒冷干燥,化學(xué)風(fēng)化能力弱。

    端元3組分的平均粒徑最大,為195.23 μm,在24~76 cm 和 232~276 cm 沉積層中含量較多,平均值分別為66.43%和49.36%,這兩層的CIA平均值分別為67.19%和65.95%,說明氣候偏溫暖濕潤,降水量豐富。由上文分析可知,這兩層沉積物顆粒以中砂為主,沉積時(shí)期水動(dòng)力較強(qiáng)。由此說明,端元3可能指示沖洪積物,以此組分為主導(dǎo)的沉積時(shí)期氣候溫暖濕潤,化學(xué)風(fēng)化能力較強(qiáng)。

    4.3 剖面年代序列

    西團(tuán)瓢剖面中存在典型的泥炭層和完整的螺殼,有機(jī)質(zhì)含量很高,因此采用14C測年。測年材料分別選取96 cm深處的植物殘?bào)w與186 cm深處的螺殼,測得年齡分別為 1 170±30、11 510±40 cal.aBP,利用線性內(nèi)插、外推的方法計(jì)算出各層位年代。分別為 292、901、5 765、10 360、15 644、16 792、21 847、30 808 cal.aBP。0~96 cm 的沉積速率為 8.20 mm/a,96~354 cm 的沉積速率為 0.87 mm/a。雖然頂部(24~74 cm)與底部(232 cm 以下)的巖性較為一致,但是其沉積過程不同。頂部沉積物的組分以端元3為主,沖洪積過程水動(dòng)力大,攜帶物質(zhì)多,單位時(shí)間內(nèi)沉積層厚,因此沉積速率大。底部沉積物的組分以端元1為主,風(fēng)成沉積過程動(dòng)力小,攜帶物質(zhì)少且顆粒細(xì)小,單位時(shí)間內(nèi)沉積層薄,沉積速率小。

    4.4 區(qū)域?qū)Ρ?/h3>

    遼東半島位于東部季風(fēng)氣候區(qū),前人對(duì)中國東北地區(qū)晚更新世以來古氣候研究多以泥炭孢粉研究為主,并結(jié)合TOC、TN、干密度、含水量及干燥度等指標(biāo)獲得了晚更新世以來的氣候演化信息。吉林二龍灣瑪珥湖記錄的氣候演化過程自晚更新世以來可分為8個(gè)階段[29],遼寧省南部普蘭店剖面記錄的氣候自全新世以來劃分了3個(gè)階段[1],遼寧省南部大蓮泡剖面記錄了全新世以來該區(qū)的沉積環(huán)境演變以及氣候演化過程[2]。本文將西團(tuán)瓢剖面記錄的氣候演化過程劃分了5個(gè)階段,現(xiàn)將3個(gè)剖面記錄的氣候演化過程作以下對(duì)比,如圖10。

    15.64 ~30.81 kaBP,西團(tuán)瓢剖面沉積層以細(xì)顆粒為主,CIA指數(shù)多指示中等的化學(xué)風(fēng)化水平,燒失量整體偏低,表示有機(jī)質(zhì)含量較少,說明此階段氣溫低,降水量少,沉積動(dòng)力較弱,氣候條件整體呈現(xiàn)寒冷和干旱的狀態(tài),但是,在16.79~15.64 kaBP,沉積物顆粒變粗,說明經(jīng)歷了短暫的降水激增的過程。二龍灣植被以冷杉、云杉等針葉林和草本植物為主,沉積物顆粒細(xì),TOC、TN含量下降,含水量低、干密度大,氣候寒冷干旱,不利于植物生長;在植被稀疏的狀態(tài)下,少量的雨水沖刷坡面攜帶較細(xì)的碎屑物進(jìn)入湖盆,同時(shí)伴有東亞季風(fēng)帶來的西北內(nèi)陸地區(qū)的細(xì)顆粒組分,這一時(shí)期湖盆內(nèi)的沉積物整體偏細(xì)。

    圖10 西團(tuán)瓢與二龍灣[29]、普蘭店[1]、大蓮泡[2] 環(huán)境指標(biāo)變化曲線對(duì)比圖Fig.10 Comparision of environmental index change curves of Xi tuanpiao,Er longwan[29],Pu landian[1] and Da lianpao[2]

    10.36 ~15.64 kaBP,西團(tuán)瓢沉積層以細(xì)顆粒為主,CIA指數(shù)指示低等化學(xué)風(fēng)化程度,燒失量為21.60%,這是因?yàn)樵搶映练e了大量螺殼,導(dǎo)致該層有機(jī)質(zhì)含量高。由此說明此階段氣候寒冷干燥。二龍灣在16.1~11.2 kaBP植被仍以針葉林和草本植物為主,氣候寒冷干旱。這一時(shí)期出現(xiàn)淡水螺殼且沉積物顆粒偏細(xì),說明坡面徑流作用減少,湖泊處于相對(duì)穩(wěn)定的環(huán)境。

    5.77 ~10.36 kaBP,西團(tuán)瓢的沉積物以細(xì)顆粒為主,CIA指數(shù)表示中等偏高的化學(xué)風(fēng)化程度,燒失量低于剖面平均水平。說明此階段氣候溫暖但降水量偏少。普蘭店剖面在10.0~8.0 kaBP榆和樺等闊葉樹種發(fā)展起來,植被為針闊混交林,年均溫度較低,濕潤度稍差,氣候總體溫和干燥。與西團(tuán)瓢剖面較為一致。同一時(shí)期大蓮泡剖面的沉積物也為灰色粉砂層,從孢粉分析看,植物孢粉種類少,且多為耐寒耐旱的植物花粉。據(jù)此可推測湖盆當(dāng)時(shí)依然受末次冰期的影響,溫度較低,集中性的降水沖刷山體,形成溝谷泥石流攜帶碎屑物進(jìn)入湖盆。此階段海岸線向內(nèi)陸延伸,海水上溯三十里堡河影響湖盆水質(zhì),致使植物種類減少,多生長鹽生蒿屬。

    0.90 ~5.77 kaBP,西團(tuán)瓢沉積物以細(xì)顆粒為主,CIA指數(shù)指示中等偏高的化學(xué)風(fēng)化程度;燒失量為剖面各層位的最高值,該層是泥炭層,由植物腐化形成,說明沉積時(shí)期的氣候溫暖濕潤,適宜于植物的生長。二龍灣的植被類型以落葉闊葉林為主,TOC、TOC/TN處于高值,含水量上升,氣候整體表現(xiàn)為溫暖濕潤。此時(shí),湖盆內(nèi)水草豐美,植被繁茂,水體平靜,坡面沖刷逐漸消失,沉積物主要由河流攜帶進(jìn)入湖盆內(nèi),粗顆粒物質(zhì)來源減少,湖內(nèi)沉積了大量細(xì)粒物質(zhì)。

    0.90 kaBP至今,西團(tuán)瓢沉積物顆粒相對(duì)較粗,CIA指數(shù)表示中等偏低的化學(xué)風(fēng)化程度,燒失量略低于剖面整體平均值,而且目前湖泊處于干涸的狀態(tài)。說明該時(shí)期氣溫偏涼,降水偏少。二龍灣在1.7 kaBP至今,區(qū)內(nèi)再次出現(xiàn)針葉林,草本植物增加,孢粉濃度下降,氣候向涼、干的方向發(fā)展。普蘭店在晚全新世時(shí)期生長針闊混交林,0.7 kaBP至今,氣候呈溫涼偏干的狀態(tài)。此階段,早期大蓮泡剖面中水生植物孢粉減少,草本植物孢粉增加,說明水體逐漸減少,湖盆處于萎縮的狀態(tài)。晚期,草本植物與水生植物孢粉都幾近消失,植被稀疏,坡面沖刷作用又得到加強(qiáng),沖洪積物比重增多。后期加上人類活動(dòng)的影響,湖體逐漸萎縮干涸。

    5 結(jié)論

    (1)15.64~30.81 kaBP,222~354 cm 沉積層以細(xì)顆粒為主,CIA指數(shù)多指示中等的化學(xué)風(fēng)化水平,燒失量整體偏低,表示有機(jī)質(zhì)含量較少,說明此階段氣溫低,降水量少,沉積動(dòng)力較小,氣候條件整體呈現(xiàn)寒冷和較干旱的狀態(tài)。

    (2)10.36~15.64 kaBP,176~222 cm 沉積層以細(xì)顆粒為主,CIA指數(shù)指示低等化學(xué)風(fēng)化程度;燒失量為21.60%,這是因?yàn)樵搶映练e了大量螺殼,導(dǎo)致該層有機(jī)質(zhì)含量高。由此說明此階段氣候寒冷干燥。

    (3)5.77~10.36 kaBP,136~176 cm 的沉積物大多為粉砂和中砂,CIA指數(shù)表示中等偏高的化學(xué)風(fēng)化程度,燒失量低于剖面平均水平。說明此階段氣候溫暖但降水量偏少。

    (4)0.90~5.77 kaBP,74~136 cm 沉積物以細(xì)顆粒為主,CIA指數(shù)指示中等偏高的化學(xué)風(fēng)化程度;燒失量為剖面各層位的最高值,該層是泥炭層,由植物腐化形成,說明沉積時(shí)期的氣候溫暖濕潤,適宜于植物的生長。

    (5)0.90 kaBP 至今,西團(tuán)瓢沉積物顆粒相對(duì)較粗,CIA指數(shù)表示中等偏低的化學(xué)風(fēng)化程度,燒失量略低于剖面整體平均值,而且目前湖泊處于干涸的狀態(tài)。說明該時(shí)期氣溫偏涼,降水偏少。

    端元分析模型根據(jù)粒度數(shù)據(jù)能夠有效地提取出沉積物動(dòng)力組分對(duì)應(yīng)的端元。本文將西團(tuán)瓢剖面沉積物粒度組分分為3個(gè)端元:端元1可能指示風(fēng)成沉積及坡面流水?dāng)y帶的細(xì)顆粒沉積,端元2可能代表動(dòng)力較低的弱水流沉積物,端元3則可能代表動(dòng)力較強(qiáng)的沖洪積物。由此可見,沉積物質(zhì)在沉積過程中往往是多種動(dòng)力因素綜合的結(jié)果,通過對(duì)沉積物中多動(dòng)力組分的分離,可以有效地判斷各組分所代表的沉積動(dòng)力的變化特征。

    猜你喜歡
    粉砂沉積物平均值
    “平均值代換”法在數(shù)學(xué)解題中的應(yīng)用
    晚更新世以來南黃海陸架沉積物源分析
    渤海油田某FPSO污水艙沉積物的分散處理
    海洋石油(2021年3期)2021-11-05 07:43:12
    水體表層沉積物對(duì)磷的吸收及釋放研究進(jìn)展
    典型粉砂地層盾構(gòu)選型及施工參數(shù)研究
    躍進(jìn)總干渠粉砂地基上節(jié)制閘抗滑穩(wěn)定性分析
    討論用ICP-AES測定土壤和沉積物時(shí)鈦對(duì)鈷的干擾
    平面圖形中構(gòu)造調(diào)和平均值幾例
    基于電流平均值的改進(jìn)無功檢測法
    電測與儀表(2014年6期)2014-04-04 11:59:46
    地下連續(xù)墻在富水粉砂層中降水加固成槽施工技術(shù)
    69av精品久久久久久| 精品久久久久久久末码| 激情 狠狠 欧美| 91久久精品国产一区二区成人| 亚洲欧美日韩卡通动漫| 中国美女看黄片| 乱人视频在线观看| 成人午夜高清在线视频| 老司机影院成人| 简卡轻食公司| 亚洲最大成人手机在线| 亚洲成人久久爱视频| 一卡2卡三卡四卡精品乱码亚洲| 性欧美人与动物交配| 非洲黑人性xxxx精品又粗又长| 非洲黑人性xxxx精品又粗又长| 色尼玛亚洲综合影院| 国产精品久久电影中文字幕| 观看免费一级毛片| 国产精品永久免费网站| 校园春色视频在线观看| 久久久久久久久久久丰满| 国产在视频线在精品| 一本精品99久久精品77| 国产欧美日韩精品一区二区| 看片在线看免费视频| av黄色大香蕉| 久久久a久久爽久久v久久| 又黄又爽又刺激的免费视频.| 乱码一卡2卡4卡精品| 又爽又黄a免费视频| 国内久久婷婷六月综合欲色啪| av黄色大香蕉| 日韩成人伦理影院| 黄色配什么色好看| 国产视频一区二区在线看| 神马国产精品三级电影在线观看| 我要看日韩黄色一级片| 99热6这里只有精品| 亚洲人成网站在线观看播放| 国产男靠女视频免费网站| 九九热线精品视视频播放| 国产真实乱freesex| 啦啦啦观看免费观看视频高清| 亚洲专区国产一区二区| 99久久精品热视频| 久久国内精品自在自线图片| 久久人妻av系列| 黄色日韩在线| 亚洲精品色激情综合| 久久久久九九精品影院| 亚洲三级黄色毛片| 嫩草影视91久久| 久久久久久久午夜电影| 国产单亲对白刺激| 人妻制服诱惑在线中文字幕| 精品久久久久久久久av| 免费在线观看成人毛片| 久久久久久久久久成人| 精华霜和精华液先用哪个| 亚洲av一区综合| 最新在线观看一区二区三区| 日韩,欧美,国产一区二区三区 | 久久精品国产鲁丝片午夜精品| 亚洲av五月六月丁香网| 亚洲一级一片aⅴ在线观看| 亚洲欧美清纯卡通| 美女被艹到高潮喷水动态| 精品一区二区三区人妻视频| a级毛片a级免费在线| 97超级碰碰碰精品色视频在线观看| 久久久欧美国产精品| 97碰自拍视频| 午夜亚洲福利在线播放| 在线播放国产精品三级| 成人毛片a级毛片在线播放| 九色成人免费人妻av| 成人av在线播放网站| a级一级毛片免费在线观看| 91午夜精品亚洲一区二区三区| 国产精品日韩av在线免费观看| 51国产日韩欧美| 午夜日韩欧美国产| 少妇人妻一区二区三区视频| 亚洲婷婷狠狠爱综合网| 日韩大尺度精品在线看网址| 亚洲精品一区av在线观看| 乱人视频在线观看| 蜜桃亚洲精品一区二区三区| 免费一级毛片在线播放高清视频| 国产精品久久久久久久电影| 久久欧美精品欧美久久欧美| 久久久久久久午夜电影| 成人二区视频| 婷婷亚洲欧美| 波多野结衣高清作品| 亚洲欧美中文字幕日韩二区| 日韩欧美在线乱码| 搡女人真爽免费视频火全软件 | 中国美白少妇内射xxxbb| 亚洲欧美日韩无卡精品| 最近中文字幕高清免费大全6| 老司机福利观看| .国产精品久久| 免费电影在线观看免费观看| 午夜激情欧美在线| 精品福利观看| 床上黄色一级片| 少妇人妻精品综合一区二区 | 在线a可以看的网站| 91麻豆精品激情在线观看国产| 午夜免费激情av| 中文字幕久久专区| 搡老岳熟女国产| 99久久无色码亚洲精品果冻| 成人一区二区视频在线观看| 亚洲乱码一区二区免费版| 精品免费久久久久久久清纯| 国产成人a区在线观看| 国产女主播在线喷水免费视频网站 | 久99久视频精品免费| 在线国产一区二区在线| 久久这里只有精品中国| 成人av在线播放网站| 欧美日本视频| 色吧在线观看| 少妇熟女欧美另类| 亚洲成人av在线免费| 97超视频在线观看视频| 最新在线观看一区二区三区| 婷婷精品国产亚洲av| 丰满人妻一区二区三区视频av| 国产真实乱freesex| 亚洲国产日韩欧美精品在线观看| www.色视频.com| 国产亚洲精品久久久com| 国产精品电影一区二区三区| av国产免费在线观看| 18禁在线无遮挡免费观看视频 | 综合色av麻豆| 久久草成人影院| 成人亚洲欧美一区二区av| 五月玫瑰六月丁香| 久久草成人影院| 国产精品亚洲一级av第二区| 嫩草影院精品99| 精品久久久久久成人av| 欧美+亚洲+日韩+国产| 一个人免费在线观看电影| 亚洲国产精品久久男人天堂| 男女啪啪激烈高潮av片| 中国美女看黄片| 亚洲丝袜综合中文字幕| 亚洲美女视频黄频| 日本撒尿小便嘘嘘汇集6| 黄色配什么色好看| 成人特级av手机在线观看| 成人二区视频| 国产女主播在线喷水免费视频网站 | 国产色婷婷99| 黄片wwwwww| 国产亚洲精品久久久com| 在线国产一区二区在线| 色吧在线观看| 99久国产av精品| 国产熟女欧美一区二区| 精品免费久久久久久久清纯| 观看免费一级毛片| 亚洲欧美日韩无卡精品| 99在线视频只有这里精品首页| 91在线观看av| 深夜精品福利| 亚洲av电影不卡..在线观看| 精品一区二区三区视频在线观看免费| 黄色视频,在线免费观看| 日本在线视频免费播放| av天堂在线播放| 卡戴珊不雅视频在线播放| 黄色欧美视频在线观看| 男人舔女人下体高潮全视频| 狂野欧美白嫩少妇大欣赏| 又黄又爽又免费观看的视频| 女生性感内裤真人,穿戴方法视频| 性欧美人与动物交配| 国产精品亚洲一级av第二区| 嫩草影院新地址| av在线老鸭窝| 一a级毛片在线观看| 精品欧美国产一区二区三| 日韩成人伦理影院| 国产成年人精品一区二区| 久久99热6这里只有精品| 亚洲av二区三区四区| 国产黄a三级三级三级人| 亚洲精品久久国产高清桃花| 看黄色毛片网站| 男女下面进入的视频免费午夜| 亚洲av免费高清在线观看| 国产伦在线观看视频一区| 色综合亚洲欧美另类图片| 亚洲成人久久爱视频| 老熟妇乱子伦视频在线观看| 欧美性猛交╳xxx乱大交人| 在线播放国产精品三级| 亚洲精品在线观看二区| 97热精品久久久久久| 波多野结衣巨乳人妻| 一区福利在线观看| 国产亚洲欧美98| 小说图片视频综合网站| 国产黄色小视频在线观看| 久久久a久久爽久久v久久| 午夜影院日韩av| 亚洲最大成人手机在线| 亚洲经典国产精华液单| 免费看美女性在线毛片视频| 女的被弄到高潮叫床怎么办| 嫩草影视91久久| videossex国产| 日韩成人伦理影院| 精品久久国产蜜桃| 蜜桃久久精品国产亚洲av| 大型黄色视频在线免费观看| 亚洲av免费高清在线观看| 不卡视频在线观看欧美| 久久人人爽人人爽人人片va| 尤物成人国产欧美一区二区三区| 蜜臀久久99精品久久宅男| 97热精品久久久久久| 欧美日韩国产亚洲二区| 日韩三级伦理在线观看| 欧美极品一区二区三区四区| 日产精品乱码卡一卡2卡三| videossex国产| 尾随美女入室| 97在线视频观看| 级片在线观看| 精品国产三级普通话版| 91精品国产九色| 1000部很黄的大片| 日本一二三区视频观看| 婷婷亚洲欧美| 精品免费久久久久久久清纯| 少妇的逼水好多| 九九久久精品国产亚洲av麻豆| 秋霞在线观看毛片| 成人av一区二区三区在线看| 免费观看的影片在线观看| 一个人看视频在线观看www免费| 亚洲自偷自拍三级| 中文在线观看免费www的网站| 亚洲乱码一区二区免费版| 1024手机看黄色片| 九九在线视频观看精品| 欧美在线一区亚洲| 伊人久久精品亚洲午夜| 听说在线观看完整版免费高清| 国产又黄又爽又无遮挡在线| 亚洲精品成人久久久久久| 久99久视频精品免费| 小说图片视频综合网站| av卡一久久| 国产精品人妻久久久久久| 亚州av有码| 极品教师在线视频| 日产精品乱码卡一卡2卡三| 亚洲最大成人手机在线| 国产亚洲精品久久久久久毛片| 久久综合国产亚洲精品| 国产精品不卡视频一区二区| 国国产精品蜜臀av免费| 中文字幕久久专区| 菩萨蛮人人尽说江南好唐韦庄 | www.色视频.com| 国产黄色视频一区二区在线观看 | 久久亚洲国产成人精品v| 精品午夜福利在线看| 国产精品久久电影中文字幕| 日韩中字成人| 蜜桃亚洲精品一区二区三区| 亚洲成人精品中文字幕电影| 可以在线观看的亚洲视频| 三级男女做爰猛烈吃奶摸视频| 欧美日韩乱码在线| 亚洲国产日韩欧美精品在线观看| 有码 亚洲区| 国产大屁股一区二区在线视频| 久久久久久九九精品二区国产| 男女之事视频高清在线观看| 蜜桃亚洲精品一区二区三区| 日韩一区二区视频免费看| 人人妻人人澡欧美一区二区| 成人无遮挡网站| 午夜视频国产福利| 少妇丰满av| 毛片女人毛片| 日韩在线高清观看一区二区三区| 欧美日本亚洲视频在线播放| 一区二区三区四区激情视频 | 中国美白少妇内射xxxbb| 久久这里只有精品中国| av卡一久久| 美女被艹到高潮喷水动态| 精华霜和精华液先用哪个| 国产亚洲av嫩草精品影院| 99精品在免费线老司机午夜| 亚洲精品一区av在线观看| 18禁在线播放成人免费| 中文字幕人妻熟人妻熟丝袜美| 国产视频一区二区在线看| 亚洲无线在线观看| 国产久久久一区二区三区| 久久精品人妻少妇| 美女 人体艺术 gogo| 99久久久亚洲精品蜜臀av| 午夜福利视频1000在线观看| 欧美最新免费一区二区三区| 六月丁香七月| 三级经典国产精品| 18禁在线无遮挡免费观看视频 | 亚洲四区av| 亚洲专区国产一区二区| 一级毛片久久久久久久久女| 国产精品人妻久久久久久| 女生性感内裤真人,穿戴方法视频| 精品一区二区三区av网在线观看| 免费不卡的大黄色大毛片视频在线观看 | 午夜福利在线在线| 久久久久性生活片| 男人舔奶头视频| 国产真实乱freesex| 亚洲欧美精品自产自拍| 国内精品宾馆在线| 我的老师免费观看完整版| 99久久精品国产国产毛片| 插逼视频在线观看| 给我免费播放毛片高清在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲乱码一区二区免费版| 在线免费观看不下载黄p国产| 1000部很黄的大片| 乱系列少妇在线播放| av在线蜜桃| 国模一区二区三区四区视频| 你懂的网址亚洲精品在线观看 | 久久6这里有精品| 久久精品国产亚洲av涩爱 | 精品无人区乱码1区二区| 午夜福利高清视频| 欧美xxxx性猛交bbbb| 日本撒尿小便嘘嘘汇集6| 亚洲成人久久性| 给我免费播放毛片高清在线观看| 我的老师免费观看完整版| 久久久色成人| 国产女主播在线喷水免费视频网站 | 国产伦在线观看视频一区| 听说在线观看完整版免费高清| 成人精品一区二区免费| a级毛片免费高清观看在线播放| 日日摸夜夜添夜夜添小说| 有码 亚洲区| 成人欧美大片| 亚洲国产欧美人成| 欧美潮喷喷水| 国产一区二区激情短视频| 欧美一级a爱片免费观看看| 免费电影在线观看免费观看| 天堂网av新在线| 免费高清视频大片| 小蜜桃在线观看免费完整版高清| 噜噜噜噜噜久久久久久91| 少妇的逼水好多| 天堂√8在线中文| 又爽又黄a免费视频| 国产精品av视频在线免费观看| 国产91av在线免费观看| 少妇丰满av| 亚洲精品456在线播放app| 国产午夜精品久久久久久一区二区三区 | 麻豆一二三区av精品| 亚洲无线在线观看| 伦精品一区二区三区| 少妇裸体淫交视频免费看高清| 色5月婷婷丁香| 欧美日韩乱码在线| 国产白丝娇喘喷水9色精品| 国产成人a∨麻豆精品| h日本视频在线播放| 欧美xxxx黑人xx丫x性爽| 成人亚洲欧美一区二区av| 亚洲婷婷狠狠爱综合网| 99视频精品全部免费 在线| 男人的好看免费观看在线视频| 夜夜看夜夜爽夜夜摸| 亚洲最大成人手机在线| 在线观看av片永久免费下载| 夜夜夜夜夜久久久久| 国产 一区精品| 亚洲av不卡在线观看| 成人国产麻豆网| 99久久成人亚洲精品观看| 美女黄网站色视频| av专区在线播放| 久久精品综合一区二区三区| 欧美+日韩+精品| 嫩草影院新地址| 伦理电影大哥的女人| 97热精品久久久久久| 欧美日韩一区二区视频在线观看视频在线 | 欧美极品一区二区三区四区| 亚洲精品日韩av片在线观看| 国产视频内射| 成人三级黄色视频| 亚洲国产精品国产精品| 国产亚洲91精品色在线| 国产在线男女| 欧美日韩乱码在线| 老熟妇乱子伦视频在线观看| 色播亚洲综合网| 国内精品宾馆在线| 国产av麻豆久久久久久久| 少妇猛男粗大的猛烈进出视频 | 国产69精品久久久久777片| 亚洲中文字幕日韩| 深爱激情五月婷婷| 91狼人影院| 综合色av麻豆| 日韩在线高清观看一区二区三区| 国产精品爽爽va在线观看网站| 99热精品在线国产| 国产欧美日韩精品亚洲av| 啦啦啦韩国在线观看视频| 亚洲精品粉嫩美女一区| 日韩在线高清观看一区二区三区| 观看美女的网站| 久久久久国产网址| 欧美xxxx黑人xx丫x性爽| 插逼视频在线观看| 99在线视频只有这里精品首页| 国产精品综合久久久久久久免费| 少妇熟女欧美另类| 精华霜和精华液先用哪个| 成年av动漫网址| 五月伊人婷婷丁香| 久久精品国产鲁丝片午夜精品| av免费在线看不卡| 91精品国产九色| 成人毛片a级毛片在线播放| 国产国拍精品亚洲av在线观看| 国产欧美日韩一区二区精品| 国产色爽女视频免费观看| 国内少妇人妻偷人精品xxx网站| 亚洲婷婷狠狠爱综合网| 国产高清三级在线| 久久久久久九九精品二区国产| 国产精品爽爽va在线观看网站| 成年女人毛片免费观看观看9| 国产精品亚洲美女久久久| 最新中文字幕久久久久| 秋霞在线观看毛片| 最新在线观看一区二区三区| 国产午夜精品论理片| 日本精品一区二区三区蜜桃| 蜜桃久久精品国产亚洲av| 国内精品美女久久久久久| 欧美性猛交黑人性爽| 老师上课跳d突然被开到最大视频| 亚洲国产高清在线一区二区三| 听说在线观看完整版免费高清| 日本一本二区三区精品| 97人妻精品一区二区三区麻豆| 成人性生交大片免费视频hd| 亚洲国产精品合色在线| a级毛片免费高清观看在线播放| 国产aⅴ精品一区二区三区波| 国产不卡一卡二| 桃色一区二区三区在线观看| 精品熟女少妇av免费看| 国产精品美女特级片免费视频播放器| 村上凉子中文字幕在线| 国产欧美日韩一区二区精品| 黄色日韩在线| 亚洲无线在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产真实伦视频高清在线观看| 中文字幕精品亚洲无线码一区| 日韩,欧美,国产一区二区三区 | 伦精品一区二区三区| 欧美性猛交╳xxx乱大交人| 麻豆一二三区av精品| 日日啪夜夜撸| 亚洲国产精品成人久久小说 | 三级男女做爰猛烈吃奶摸视频| 久久6这里有精品| 免费看光身美女| 日日干狠狠操夜夜爽| 大香蕉久久网| 欧美在线一区亚洲| 久久久久免费精品人妻一区二区| 国产av不卡久久| 免费人成在线观看视频色| 国产精品国产三级国产av玫瑰| 免费看日本二区| 亚洲国产高清在线一区二区三| 一级毛片久久久久久久久女| 亚洲av二区三区四区| 老司机影院成人| 最好的美女福利视频网| 国内久久婷婷六月综合欲色啪| 哪里可以看免费的av片| 欧美激情国产日韩精品一区| 99riav亚洲国产免费| 又粗又爽又猛毛片免费看| 特大巨黑吊av在线直播| 国产一区亚洲一区在线观看| av在线亚洲专区| 18禁裸乳无遮挡免费网站照片| 精品人妻偷拍中文字幕| 欧美成人免费av一区二区三区| 国产高清视频在线观看网站| 国产伦在线观看视频一区| 一区二区三区免费毛片| 日日啪夜夜撸| av黄色大香蕉| 露出奶头的视频| 亚洲色图av天堂| av女优亚洲男人天堂| 精品人妻视频免费看| 一本精品99久久精品77| 中文字幕熟女人妻在线| av天堂中文字幕网| 亚洲人与动物交配视频| 成人国产麻豆网| 成人欧美大片| 女人十人毛片免费观看3o分钟| 精品一区二区三区视频在线观看免费| 精品人妻熟女av久视频| 国产午夜福利久久久久久| 欧美成人一区二区免费高清观看| 搡老熟女国产l中国老女人| 国产一区二区亚洲精品在线观看| 亚洲av成人精品一区久久| 亚洲成人中文字幕在线播放| 国产淫片久久久久久久久| 亚洲精品一卡2卡三卡4卡5卡| 超碰av人人做人人爽久久| 男女那种视频在线观看| 国产男靠女视频免费网站| 日本一二三区视频观看| 国产极品精品免费视频能看的| 色视频www国产| 99久久中文字幕三级久久日本| 尾随美女入室| 最近手机中文字幕大全| 成人三级黄色视频| 精品久久久噜噜| 一夜夜www| 色哟哟哟哟哟哟| 国产精品久久久久久久久免| 十八禁网站免费在线| 黄色一级大片看看| 国产真实乱freesex| 一a级毛片在线观看| 午夜视频国产福利| 久久午夜福利片| 91午夜精品亚洲一区二区三区| 成年av动漫网址| 露出奶头的视频| 成人特级黄色片久久久久久久| 色视频www国产| 欧美精品国产亚洲| 久久久久国产网址| 狂野欧美白嫩少妇大欣赏| 婷婷六月久久综合丁香| 尾随美女入室| 99在线视频只有这里精品首页| 男女下面进入的视频免费午夜| 国产精品野战在线观看| 成人午夜高清在线视频| av在线蜜桃| 欧美极品一区二区三区四区| 91久久精品国产一区二区三区| 高清午夜精品一区二区三区 | av专区在线播放| 国产成人a∨麻豆精品| 亚洲激情五月婷婷啪啪| 亚洲精品乱码久久久v下载方式| 欧美性猛交黑人性爽| 高清午夜精品一区二区三区 | 白带黄色成豆腐渣| 国产麻豆成人av免费视频| av中文乱码字幕在线| 中文字幕人妻熟人妻熟丝袜美| 日韩成人av中文字幕在线观看 | 国产精品一区www在线观看| 久久亚洲精品不卡| 伦精品一区二区三区| 国产片特级美女逼逼视频| 男女做爰动态图高潮gif福利片| 亚洲在线自拍视频| 毛片一级片免费看久久久久| 免费看a级黄色片| 精品久久久噜噜| 美女内射精品一级片tv| 亚洲自拍偷在线| 免费观看精品视频网站| 毛片一级片免费看久久久久| 久久人人精品亚洲av| 淫妇啪啪啪对白视频| av卡一久久| 深夜a级毛片| 赤兔流量卡办理| 最近手机中文字幕大全| 在线观看一区二区三区|