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

    白云巖風(fēng)化剖面的粒度分布、元素遷移及碳同位素特征——以黔北新蒲剖面為例

    2016-04-13 01:59:08李今今季宏兵
    生態(tài)學(xué)報(bào) 2016年5期
    關(guān)鍵詞:有機(jī)碳

    高 杰,李 銳,李今今,張 莉,季宏兵,2,*

    1 首都師范大學(xué), 北京 100048

    2 中國(guó)科學(xué)院地球化學(xué)研究所, 貴陽(yáng) 550002

    白云巖風(fēng)化剖面的粒度分布、元素遷移及碳同位素特征
    ——以黔北新蒲剖面為例

    高杰1,李銳1,李今今1,張莉1,季宏兵1,2,*

    1 首都師范大學(xué), 北京100048

    2 中國(guó)科學(xué)院地球化學(xué)研究所, 貴陽(yáng)550002

    摘要:研究了影響喀斯特地區(qū)土壤有機(jī)碳(SOC)更替的因素及其與元素地球化學(xué)行為的關(guān)系。選取黔北新蒲地區(qū)的兩個(gè)白云巖風(fēng)化剖面,以粒度分布、有機(jī)碳δ13C和元素質(zhì)量遷移系數(shù)為基礎(chǔ),對(duì)該區(qū)域紅色風(fēng)化殼中有機(jī)碳及δ13C值隨深度變化的趨勢(shì)和元素的地球化學(xué)特征進(jìn)行了探討。結(jié)果顯示,該區(qū)土壤分級(jí)后粉粒(5—50 μm)含量占50%—80%,粘粒次之,砂粒最少;δ13CSOC值介于-26.4‰―-23.6‰,且分級(jí)后的土樣和原樣中δ13C值隨著剖面深度的增加均減小,各粒級(jí)間差異不明顯;草地覆被的剖面中鐵、鉀、鋇、釩和鋅表現(xiàn)出一定的富集,鈉和鈷在風(fēng)化原巖處有一個(gè)富集端點(diǎn),錳在巖粉處富集,鎂、鈣則完全遷出土體。而耕地剖面中,釩和鋅呈現(xiàn)一定程度富集,鐵、鉀在巖土界面富集,鈷在風(fēng)化原巖富集,而鋇在表層微弱富集。

    關(guān)鍵詞:有機(jī)碳;分粒級(jí);元素遷移;碳同位素

    中國(guó)西南喀斯特地區(qū)屬于典型的生態(tài)環(huán)境脆弱區(qū),加之人類不合理地開(kāi)發(fā)利用土地資源,區(qū)域生態(tài)安全面臨諸多威脅[1]。在土壤侵蝕退化與全球碳循環(huán)過(guò)程中,作為土壤重要組成部分和代表土壤主要碳庫(kù)的有機(jī)質(zhì)在生態(tài)系統(tǒng)中起重要作用[2]。對(duì)土壤碳循環(huán)機(jī)制的研究是探討陸地生態(tài)系統(tǒng)碳循環(huán)的重要前提,也是深入理解陸地生態(tài)系統(tǒng)對(duì)全球變化響應(yīng)機(jī)理的重要基礎(chǔ)[3]。SOC 垂直分布和循環(huán)特征及其與土壤中元素的相互關(guān)系是研究土壤地球化學(xué)循環(huán)及氣候變化的重要領(lǐng)域[4]。土壤中Ca,F(xiàn)e,Mg,Cu,Zn,K,S和P等元素循環(huán)和分布特征反映了土壤性質(zhì),同時(shí)這些元素也與土壤生物和微生物及地表植被有緊密聯(lián)系,土壤元素的分布特征也影響著土壤有機(jī)碳的循環(huán)速率[5]。當(dāng)前對(duì)喀斯特山地石漠化與風(fēng)化成土過(guò)程的地球化學(xué)特征的研究較為深入,但從地球化學(xué)角度系統(tǒng)地對(duì)黔北地區(qū)土壤有機(jī)碳與元素地球化學(xué)循環(huán)的關(guān)系及其影響因子的報(bào)道還比較少,開(kāi)展這方面的研究有利于提高人們?cè)谕寥涝氐厍蚧瘜W(xué)循環(huán)對(duì)SOC含量和垂直分布的影響機(jī)理方面的認(rèn)識(shí)。研究選取黔北新蒲兩個(gè)白云巖風(fēng)化剖面,根據(jù)SOC含量及其δ13C 組成差異,結(jié)合土壤pH、顆粒組成及元素地球化學(xué)特征等分析結(jié)果,探討土壤內(nèi)在性質(zhì)對(duì)SOC在風(fēng)化剖面中更替的影響及其與元素地球化學(xué)的關(guān)系,為揭示喀斯特地區(qū)碳循環(huán)和元素地球化學(xué)特征提供基礎(chǔ)資料。

    1剖面特征及分析方法

    1.1研究區(qū)概況

    研究區(qū)位于黔北新蒲巖溶山地背景上的丘陵臺(tái)地,本區(qū)屬中亞熱帶氣候,多年平均氣溫15.1 ℃,年降水量1000—1150,日照少,陰雨多霧,區(qū)內(nèi)地形起伏較大,地貌類型復(fù)雜。自然土壤以黃壤、黃色石灰土為主,成土母質(zhì)為碳酸鹽巖、砂頁(yè)巖。區(qū)內(nèi)分布有馬尾松林、混交林,也有部分灌叢、草被出現(xiàn)。所研究的白云巖剖面是寒武系婁關(guān)山群白云巖經(jīng)自然風(fēng)化形成,基巖面起伏強(qiáng)烈,巖土界面清晰。剖面所處位置地質(zhì)背景及示意圖如圖1,研究區(qū)地質(zhì)圖引自《貴州省1∶50萬(wàn)地質(zhì)圖》。

    2013年7月,在研究區(qū)內(nèi)取兩個(gè)剖面(27°40.174′N,107°00.170′E;27°41.142′N,107°00.26′E),分別編號(hào)為X和Y。剖面X從基巖到地表厚450 cm,垂向上劃分為表土層(A)、土層(B)、巖粉層(C)、風(fēng)化原巖層(D)和基巖層(E)。A層為黑褐色,厚25 cm左右,含大量植物根系。B層為紅褐色土層,是剖面的主體部分,厚約4 m,剖面245—165 cm處發(fā)現(xiàn)少許黑色物質(zhì),125 cm以上部分發(fā)現(xiàn)微量根系殘留物。C層呈灰褐色,厚約15 cm。表層植被為草叢;剖面Y表層為新開(kāi)墾的耕地,厚2.5 m,未見(jiàn)有明顯的巖粉層,風(fēng)化原巖易碎,用手輕搓即成粉砂狀,A層為黃褐色表土層,B層紅褐色,地表植被為玉米。由于邊坡有著近70°的傾角,認(rèn)為剖面近似垂直。

    圖1 采樣區(qū)地質(zhì)背景及剖面示意圖Fig.1 Geological background of the sampling area and profile sketch

    1.2取樣及樣品分析

    為了消除坡積物的影響,采樣過(guò)程中去除剖面表層覆被物,按15 cm間隔自下向上采集5 cm厚新鮮土層,并采集了基巖、風(fēng)化原巖和巖粉,土樣風(fēng)干備用。

    土壤pH值的測(cè)定參考“中華人民共和國(guó)農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)”;取研磨過(guò)100目篩的土樣1 g,用1 mol/L的鹽酸溶液浸泡24 h,以去除土壤中的碳酸鹽,然后用去離子水洗滌至中性,60 ℃烘干后研磨,用TOC分析儀測(cè)定土壤中的有機(jī)碳含量;基于中國(guó)科學(xué)院南京土壤研究所等單位制定的土粒分級(jí)標(biāo)準(zhǔn),以去離子水為分散劑,采用連續(xù)分級(jí)法,結(jié)合使用篩分、沉降、離心及錯(cuò)流超濾(CFF)等技術(shù),將過(guò)60目篩的土樣提取分級(jí)為細(xì)砂粒(50—250 μm)、粉粒(5—50 μm)、粗粘粒(1—5 μm)、細(xì)粘粒(<1 μm)4個(gè)組分[6];根據(jù)中國(guó)地質(zhì)科學(xué)院地質(zhì)實(shí)驗(yàn)測(cè)試中心標(biāo)準(zhǔn)Q/GD001-2002-巖石、土壤、水系沉積物中微量元素的密閉溶樣-電感耦合等離子質(zhì)譜法測(cè)定土壤中常、微量元素,測(cè)試均用標(biāo)樣控制,測(cè)試結(jié)果誤差為±15%;土壤有機(jī)質(zhì)穩(wěn)定碳同位素的測(cè)定參見(jiàn)文獻(xiàn)[7]。將預(yù)處理的土樣送中國(guó)科學(xué)院地球化學(xué)研究所,用MAT252型質(zhì)譜儀測(cè)定樣品有機(jī)質(zhì)的δ13C值(穩(wěn)定碳同位素組成),分析誤差小于±0.1‰。

    2結(jié)果與討論

    2.1剖面的SOC和pH值分布特征

    SOC作為土壤的重要組成部分,在維持土壤的物理、化學(xué)和生物特征中起重要作用,土壤pH值則通過(guò)影響微生物活性而影響有機(jī)質(zhì)降解。自然條件下,表層土壤是植被枯枝落葉的主要接受層,在這一層土壤獲得大量腐殖化、半腐殖化有機(jī)物和新鮮枯枝落葉等有機(jī)碳。隨著土壤深度的加深,SOC的來(lái)源轉(zhuǎn)為植物根系和微生物量,以及表層SOC向下遷移,其來(lái)源的量相對(duì)表層直接接受枯枝落葉所積累的SOC明顯降低[8]。因此,表層SOC含量為整個(gè)剖面最高。

    由圖2,兩剖面中SOC的含量均低于11 g/kg。SOC含量存在一個(gè)最高臨界值為20 g/kg,當(dāng)?shù)陀诖酥禃r(shí),土壤結(jié)構(gòu)將變得不太穩(wěn)定,土地生產(chǎn)力也將降低[9]。說(shuō)明該區(qū)土壤質(zhì)量已經(jīng)下降到對(duì)土壤結(jié)構(gòu)及土地生產(chǎn)不利的程度,這也說(shuō)明該區(qū)土壤侵蝕比較嚴(yán)重,養(yǎng)分流失較多。剖面X中,SOC在0—25 cm范圍急劇增加。以25 cm為分界線,SOC含量急劇減少。這可能是因?yàn)樵撈拭姹韺又脖粸椴輩?,大量的根系殘留?0—25 cm深度的緣故。從25 cm深度到巖粉,SOC含量隨剖面深度的變化呈倒“L”型,即隨剖面深度的加深,SOC含量減少。土體內(nèi)部SOC的累積主要依靠表層SOC溶解于水或吸附于土壤粘粒表面向下遷移,或依靠土體內(nèi)部根系分泌物和微生物量的結(jié)果。SOC在土體內(nèi)部的遷移累積形式非常緩慢,具有層次性,使得剖面X中SOC含量隨深度增加而降低;剖面Y中SOC從表層到65 cm深度急劇減少,65—125 cm范圍SOC略微增加,125—165 cm逐漸減少,165—225 cm范圍略微增加,而到半風(fēng)化原巖又驟然減少。即在125 cm到半風(fēng)化原巖層SOC呈“S”型分布。這可能是因?yàn)槿祟惢顒?dòng)對(duì)土壤的擾動(dòng)(翻耕以及對(duì)作物秸稈的回收等),減少了有機(jī)質(zhì)的輸入,改變了土壤物理性狀,特別是耕作層土壤的通透性,加速了微生物對(duì)這一層次土壤有機(jī)質(zhì)的降解,從而使得表層土壤SOC到65 cm深度急劇減少。而65 cm以下土壤中SOC含量的波動(dòng)則可能是因?yàn)橥鈦?lái)物質(zhì)的加入,擾動(dòng)了SOC在剖面中的遷移轉(zhuǎn)化。

    表層土壤的pH值均呈酸性,且剖面X較Y酸性強(qiáng)(圖2)。隨著深度的增加剖面X的pH值逐漸升高,145—165 cm深度時(shí)突然減小,但很快繼續(xù)升高,一直到巖土界面,pH值接近中性;剖面Y中pH值隨深度增加而升高的趨勢(shì)并不明顯。在0—85 cm,pH值穩(wěn)定在5.2—5.4,之后pH值逐漸升高,一直到風(fēng)化原巖pH值為6.5。剖面表層土壤pH值差異明顯,可能和地表不同植被有關(guān)。

    圖2 剖面中SOC和pH分布特征Fig.2 Distribution characteristics ofsoil organic carbon (SOC) and pH

    2.2剖面的粒度分布特征

    土壤顆粒是構(gòu)成土壤固相骨架的基本顆粒,它們數(shù)目眾多,大小和性狀迥異,按照不同的比例進(jìn)行組合形成不同的土壤質(zhì)地類型,進(jìn)而通過(guò)影響土壤的物理、化學(xué)和生物條件來(lái)改變物質(zhì)循環(huán)和能量流動(dòng)。因此,定量描述土壤的顆粒組成在研究物質(zhì)循環(huán)和能量流動(dòng)過(guò)程中有著重要的現(xiàn)實(shí)意義[10]。

    由圖3,總體上兩個(gè)剖面分級(jí)后的土粒基本集中于粉粒(5—50 μm),最高為Y表層,為86.4%,最低為剖面X125—145 cm處,但是也占到49.1%。剖面X中,各粒級(jí)顆粒含量隨剖面深度的增加波動(dòng)明顯,其中粉粒含量在0—25 cm土壤中明顯下降,在25—145 cm處呈先增加后減少的趨勢(shì),在145—305 cm又逐漸增加,305—425 cm之間基本穩(wěn)定;粘粒含量在0—365 cm深度范圍波動(dòng)明顯,且總體呈下降趨勢(shì),但是之后又逐漸升高,到剖面底部粘粒含量達(dá)到46%。砂粒含量在整個(gè)剖面中變化較粉粒和粘粒稍顯和緩,從表層土壤到65cm處先減少再增加,在65 cm深度達(dá)到20.4%,65—145 cm處,逐漸減少到2.6%,之后一直到325 cm處波動(dòng)明顯,從325 cm深度處一直到巖土界面,砂粒含量一直減少,其中巖土界面最低,為0.3%。

    剖面Y中各粒級(jí)含量隨剖面深度的變化趨勢(shì)較剖面X緩和。表層到85 cm深度,粉粒隨剖面深度的增加而減少,變幅為14%,從85 cm到185 cm處粉粒含量在45%—80%范圍波動(dòng)。之后粉粒含量逐漸增加,到巖土界面達(dá)到最大86.4%。而粘粒的變化趨勢(shì)與粉粒基本相反。砂粒在整個(gè)剖面中變化平緩,從表層到25 cm先減少,25 cm到125 cm逐漸增加,達(dá)到13.9%,從125 cm到185 cm又逐漸減少為0.2%,之后又增加,到半風(fēng)化原巖,砂粒含量達(dá)13.4%。

    圖3 剖面X和Y粒度分布圖Fig.3 Partical size distribution of profile-X and profile-Y

    顆粒組成的深度分布反映了成土淋溶淀積的特點(diǎn),在土壤剖面長(zhǎng)期成土演化的過(guò)程中,表層土粘粒含量最低,隨剖面深度的增加不斷升高,直至出現(xiàn)明顯的粘粒富集層[3]。而據(jù)本文研究,土壤顆粒主要為粉粒,這與張冬青等[11]的研究結(jié)果一致。粘粒含量在整個(gè)剖面中呈現(xiàn)有規(guī)律的波動(dòng),只在較淺層土壤中(X:0—145 cm;Y:0—85 cm)與上述說(shuō)法相符,即隨著剖面深度增加而粘粒含量增加,且兩個(gè)剖面中粘粒含量隨剖面的變化趨勢(shì)與粉粒含量隨剖面的變化趨勢(shì)呈現(xiàn)此消彼長(zhǎng)的態(tài)勢(shì),這可能是因?yàn)樵诔赏吝^(guò)程中,兩個(gè)剖面一方面受成土母質(zhì)的影響,另一方面也說(shuō)明了剖面中的物質(zhì)移動(dòng)或者外來(lái)物質(zhì)的加入,從而使其深度排列變的復(fù)雜。

    對(duì)于均一母巖上發(fā)育的風(fēng)化殼,理想的剖面(假定不發(fā)生淋溶淀積作用)變化趨勢(shì)自下向上為小于5 μm含量(粘粒)逐漸增大[12]。剖面Y中125—85 cm處與半風(fēng)化原巖到205 cm處表現(xiàn)出均一母巖發(fā)育的典型殘積風(fēng)化殼特征,即從下向上,粘粒含量逐漸增加,但是在85 cm到表土層,粘粒含量自下向上呈減少的趨勢(shì);而剖面X在385—145 cm處與65—25 cm處自下向上在整體上保持風(fēng)化殼正向演化方向(即粘粒含量呈增加趨勢(shì))的基礎(chǔ)上,出現(xiàn)一定幅度的波動(dòng)。從巖粉到385 cm處與145—65 cm范圍粘粒含量自下向上逐漸減少。粘粒含量自下向上逐漸減少的現(xiàn)象可能是風(fēng)化殼表層受大氣水垂直下滲的影響,由于機(jī)械淋溶導(dǎo)致粘粒的向下遷移[13]。由裂隙和碎屑礦物組成滲流通道,便于粘粒礦物的遷移,隨著風(fēng)化程度增強(qiáng),剖面逐漸變厚,淀積的深度隨之加深。同時(shí),隨著大量粘粒礦物的形成,原生碎屑變少、變小,滲流管道變窄,變少,從而阻礙了細(xì)粒物質(zhì)的移動(dòng),于是出現(xiàn)了剖面越厚,風(fēng)化程度越強(qiáng),反之亦然[12]。

    剖面中部以下粒度參數(shù)出現(xiàn)的波動(dòng)一方面是受水文條件變化的影響,如地下水位及地下水的溶蝕能力。巖-土界面發(fā)生的風(fēng)化作用代表生物地球化學(xué)過(guò)程,剖面下部發(fā)育時(shí)代較新,風(fēng)化程度低,不穩(wěn)定礦物含量高,抗物理化學(xué)作用弱,對(duì)巖溶地下水條件變化反應(yīng)敏感,所以剖面下部粒度參數(shù)的波動(dòng),可以看作是當(dāng)?shù)厮牡刭|(zhì)條件的指標(biāo);另一方面,也可能是基巖的垂向不均一,導(dǎo)致了上覆風(fēng)化殼粒度分布的差異性[12]。

    2.3剖面的元素分布特征

    成土過(guò)程中,土壤母質(zhì)中的常、微量元素經(jīng)歷一系列物理、化學(xué)和生物化學(xué)反應(yīng),使其形態(tài)及含量發(fā)生轉(zhuǎn)化、移動(dòng)和再分配,構(gòu)成了土壤中元素的循環(huán)。元素的濃度變化不能真實(shí)地反映風(fēng)化成土過(guò)程中元素的地球化學(xué)行為,因?yàn)樵诨瘜W(xué)風(fēng)化過(guò)程中,由于體積的變化可能導(dǎo)致樣品中相對(duì)穩(wěn)定元素含量的增加,從而掩蓋這一過(guò)程中元素活動(dòng)的真實(shí)信息。因此,為了表示風(fēng)化成土過(guò)程中元素的地球化學(xué)行為,目前主要是采用質(zhì)量平衡逼近的方式來(lái)實(shí)現(xiàn),計(jì)算公式見(jiàn)文獻(xiàn)[14]。

    一般認(rèn)為,Ti、Al、Fe、Zr、Th在風(fēng)化過(guò)程中相對(duì)穩(wěn)定,因此常被選作質(zhì)量平衡計(jì)算中的參比元素。但是Zr、Th在巖石中含量較低且分布不均勻,而Ti是常量元素,在樣品中含量高且分布均勻,取樣和測(cè)試分析產(chǎn)生的誤差小,另從其在土層、酸不溶物中的含量來(lái)看,Ti有很強(qiáng)的穩(wěn)定性[15],因此本文選取Ti作為不活化性參比元素,以原巖作對(duì)比基礎(chǔ)。各元素在剖面中的遷移變化如圖4、圖5所示。

    由圖4,剖面X中常量元素鈣、鎂相對(duì)Ti顯著遷出(或稱虧損),并且在巖土界面以上,元素的質(zhì)量遷移系數(shù)基本穩(wěn)定,只在巖土過(guò)渡帶上有一個(gè)驟變。從原巖到巖粉,鈉呈現(xiàn)明顯遷入(或稱富集),而巖粉以上則表現(xiàn)為遷出。鐵除了在剖面表層(0—85 cm虧損程度漸弱)、風(fēng)化原巖和極個(gè)別土層中呈微弱的虧損之外,在整個(gè)剖面呈不同程度的富集態(tài)勢(shì)。鉀從原巖到表層整體呈現(xiàn)遷出的特征。從巖土界面到325 cm處,遷移系數(shù)波動(dòng)明顯,且有兩個(gè)弱富集端點(diǎn),在165—125 cm處,遷出趨勢(shì)減緩,之后恢復(fù)正常;微量元素的變化最明顯的有錳、鈷元素。從原巖到巖粉,鈷呈顯著富集,而錳從原巖到風(fēng)化原巖為虧損,風(fēng)化原巖到巖粉表現(xiàn)為富集。巖粉以上整個(gè)剖面,錳、鈷都呈顯著虧損;釩從原巖到巖粉表現(xiàn)為虧損,從巖土界面到剖面中上部(105 cm處)呈現(xiàn)一定程度的富集,而從剖面中上部到表土層富集程度逐漸減弱,甚至表現(xiàn)為虧損且程度逐漸加強(qiáng);鋇和鋅在整個(gè)剖面上總體表現(xiàn)為虧損,只有鋇在剖面中下部有兩個(gè)富集端點(diǎn)。

    圖4 剖面X中主、微量元素的分布特征Fig.4 Distribution characteritics of major elements and trace elements in profile-X

    從圖5看出,剖面Y中除了鐵和鉀從原巖到巖土界面有微弱富集外,在整個(gè)剖面中鐵、鉀、鈉、鈣、鎂等常量元素均為虧損,且鈣和鎂、鈉和鉀的變化趨勢(shì)相似,而鐵的遷移分布則在整個(gè)剖面有一定的微弱波動(dòng);微量元素中錳和鋇在剖面上整體虧損。其中錳從原巖到巖土界面,虧損趨勢(shì)明顯,巖土界面以上基本完全遷出;鈷在風(fēng)化原巖處有一個(gè)富集端點(diǎn);鋇從原巖到巖土界面虧損趨勢(shì)明顯,而巖土界面以上虧損趨勢(shì)趨于減弱,到表層有微弱富集;釩和鋅在表層和145—205 cm處遷出,而在25—125 cm和風(fēng)化原巖及巖土界面處表現(xiàn)為富集。

    圖5 剖面Y中主、微量元素的分布特征Fig.5 Distribution characteritics of major elements and trace elements in profile-Y

    元素虧損或零獲得正是風(fēng)化殼剖面的基本特征,即部分元素遷出剖面并隨流體進(jìn)入地表或地下水系統(tǒng),部分元素保持穩(wěn)定[16]。剖面中鈣、鎂、鈉的遷移系數(shù)從原巖到巖土界面的土壤中很快就降低,反映了白云石迅速溶解過(guò)程中可溶性離子淋失現(xiàn)象,與大多數(shù)風(fēng)化剖面在巖土界面上發(fā)生較大的元素丟失相近(如Ca, Na, Si等)[17]。剖面中鐵的富集現(xiàn)象則可能是上部富鐵土層(現(xiàn)已被剝蝕)的物理或化學(xué)遷移所致。剖面中鉀的虧損反映出隨著風(fēng)化的進(jìn)行,含鉀礦物不斷的淋失。巖土界面附近出現(xiàn)的波動(dòng)以及弱富集層有3種可能的解釋,一是受當(dāng)?shù)厮牡刭|(zhì)條件及土壤pH值的影響,剖面中鉀的固定與釋放反復(fù)進(jìn)行的結(jié)果;二是剖面演化過(guò)程中環(huán)境變遷的原因,如大氣沉降等;三是風(fēng)化母巖是不均一的,這也從側(cè)面證實(shí)了粒度分析的結(jié)果,在文獻(xiàn)[12]中亦有所體現(xiàn)。

    鋇的主要賦存礦物為長(zhǎng)石、碳酸鈣等,在風(fēng)化剖面演化中它們是持續(xù)溶解的,故Ba元素隨著賦存礦物的淋失而產(chǎn)生虧損。但是兩個(gè)剖面中都表現(xiàn)出了鋇的富集端元,這一方面說(shuō)明了外來(lái)物質(zhì)的加入,另一方面也證明了母巖的不均一性。鈷和釩為變價(jià)元素,它們?cè)诒砩h(huán)境很容易被氧化,釩易被粘土礦物吸附難以遷移,鈷則在酸性氧化環(huán)境溶解,堿性環(huán)境沉淀。當(dāng)原生礦物氧化時(shí),鋅可以進(jìn)入淋濾溶液,但它的遷移受吸附作用控制。在酸性土壤中,錳以兩價(jià)Mn2+為主,中性附近有利于三價(jià)Mn3+的生成,而當(dāng)pH值>8時(shí),大部分向MnO2轉(zhuǎn)化[18]。而兩剖面土壤pH均為酸性到中性,故剖面中錳以遷出為主。

    2.4剖面中δ13CSOC的分布特征

    土壤中的有機(jī)碳主要來(lái)自陸生高等植物,由于SOC分解過(guò)程中的分餾作用很小[19],而且碳同位素分餾產(chǎn)生于早期的植物殘?bào)w分解階段[20],經(jīng)過(guò)長(zhǎng)期分解穩(wěn)定后的土壤總有機(jī)碳的碳同位素組成是穩(wěn)定的,因此SOC中的δ13C值與作為該有機(jī)碳來(lái)源的植被的δ13C值基本一致[21-22],能夠反映局部地區(qū)植被碳固定過(guò)程的同位素組成。土壤有機(jī)碳含量及其δ13C值隨深度變化的趨勢(shì)可以反映植物殘?bào)w的輸入及其在土壤中的分解累積特征,有助于揭示SOC循環(huán)過(guò)程及規(guī)律[23]。

    粒度分級(jí)后的土樣只測(cè)試了原樣和粉粒的全部樣品(圖6)和剖面表層、中部和下部的砂粒與粘粒的δ13C值(表1)。由圖6,兩個(gè)剖面中δ13C值隨著剖面深度的增加而減小??傮w來(lái)看,剖面X中δ13C值隨剖面深度的增加波動(dòng)明顯,剖面表層和巖土界面的δ13C值最大相差2.8‰,在0—25 cm深度δ13C值明顯減小,25—45 cm基本保持不變。45 cm往下,δ13C值隨著剖面深度的增加呈“S”形不規(guī)則波動(dòng),僅在65—85 cm和140—165 cm處δ13C值差異小于0.06‰。粉粒的δ13C值較原樣微弱偏正,但是其隨剖面深度增加的變化趨勢(shì)與原樣基本一致。砂粒和粘粒雖然只測(cè)試了剖面表層、中部和下部的樣品,但是依然可以看出,粘粒較砂粒的δ13C值在表層偏負(fù),而中部和巖土界面則是較砂粒偏正。兩者的δ13C值總體變化趨勢(shì)和原樣基本一致。

    表1 兩剖面表層、中部和巖土界面δ13CSOC值

    Y,X:研究區(qū)兩個(gè)剖面; “—”表示未測(cè)出

    圖6 原樣及粉粒δ13CSOC值的分布特征 Fig.6 δ13CSOC distribution characteristics of original samples and silt samples

    剖面Y原樣風(fēng)化原巖的δ13C值未測(cè)出。整體上δ13C值隨剖面深度增加而減小,變幅在1‰左右。在0—45 cm原樣的δ13C值基本穩(wěn)定。在45 cm以下δ13C值隨剖面深度呈不規(guī)則波動(dòng),但是均集中在-25‰—-26‰之間。粉粒中,在0—25 cm深度δ13C值偏正,25 cm以下δ13C值隨剖面深度的增加緩慢偏負(fù),僅在85 cm和145 cm處微弱偏正。粘粒在剖面表層、中部和下部的δ13C值整體偏負(fù),但是變化不明顯,變幅在0.48‰。砂粒則表現(xiàn)出先偏負(fù)再偏正的趨勢(shì)。而且砂粒和粘粒在表層和中部δ13C值均較原樣和粉粒偏負(fù),而在剖面底部即半風(fēng)化原巖的δ13C值較粉粒偏正。剖面X與剖面Y相比較,X的δ13C值較Y偏負(fù)土壤剖面成土演化對(duì)有機(jī)質(zhì)深度分布具有顯著制約作用[24]。雖然兩個(gè)剖面發(fā)育明顯,土壤厚度較大,但是,在相似的地形條件下,不同層段土壤的成土環(huán)境也可能不相同,從而使土壤有機(jī)質(zhì)分解程度出現(xiàn)差異,這可能是導(dǎo)致兩剖面的δ13CSOC值深度分布不同的主要因素。 表層土壤有機(jī)質(zhì)主要來(lái)源于地表植被的枯枝落葉。剖面X表層土壤δ13C值高于Y,可能是X表層土壤植被為草叢,當(dāng)?shù)貎?yōu)勢(shì)植物種類為C4光合作用的草本植物——芒草,從而使其土壤有機(jī)質(zhì)的δ13C值高于Y表層土壤的δ13C值。隨著深度的增加,剖面X中δ13C值開(kāi)始減小,而剖面Y中δ13C值稍微增加,隨后向減小的方向波動(dòng)。這可能是因?yàn)椋?1)植物新陳代謝產(chǎn)生的分餾作用,使組成土壤有機(jī)質(zhì)不同組分的δ13C值具有一定的變化范圍。相對(duì)于整體植物組織,脂類、木質(zhì)素和纖維素通常虧損13C,而糖類、氨基酸、半纖維素和果膠富集13C[25],由于其相對(duì)化學(xué)穩(wěn)定性不同,在土壤有機(jī)質(zhì)分解早期,富集13C的糖類和氨基酸等物質(zhì)優(yōu)先被降解,虧損13C的脂類和木質(zhì)素可能被優(yōu)先積累,導(dǎo)致它們的含量隨著土層深度的增加而升高,從而使土壤有機(jī)質(zhì)的δ13C值隨著剖面深度的增加而降低。(2)在質(zhì)地黏重的土壤中,淋溶淀積作用強(qiáng)烈,由于毛管作用,下部富13C的有機(jī)質(zhì)可能向上遷移,由于土壤粘粒的聚集導(dǎo)致δ13C值偏負(fù)。(3)據(jù)Guggenberger[26]等的研究在pH值小于5的土壤中,木質(zhì)素的含量明顯偏高,而在pH值大于5的土壤中很少有木質(zhì)素的存在。剖面中δ13C隨剖面深度的增加有偏負(fù)的趨勢(shì),可能是酸性水解反應(yīng)去除了土壤有機(jī)質(zhì)中富集13C的組分,木質(zhì)素等虧損13C的組分在土壤中積累所致。δ13C值在剖面中的不規(guī)則波動(dòng)可能是兩剖面所處的坡度較大,強(qiáng)烈的侵蝕堆積過(guò)程使得剖面內(nèi)部新老有機(jī)質(zhì)混雜所致。這說(shuō)明了土壤有機(jī)質(zhì)δ13C值剖面分布的復(fù)雜性。

    由圖6,兩剖面中粉粒的δ13C值(除了Y下部的幾個(gè)點(diǎn))在整個(gè)剖面中均高于原樣,與粉粒結(jié)合的有機(jī)質(zhì)具有最高的δ13C值。而砂粒δ13C值在剖面表層土壤、中部均較粉粒和原樣偏負(fù)。與砂粒結(jié)合的表層土壤有機(jī)質(zhì)的δ13C分別為X:-24.66 ‰;Y:-25.11 ‰。土壤中與砂粒相結(jié)合的有機(jī)質(zhì)稱為顆粒有機(jī)質(zhì),通常認(rèn)為是由部分分解的植物殘?bào)w所組成,處于有機(jī)質(zhì)分解的早期階段,因而由分解作用所產(chǎn)生的同位素分餾很小[27],它的δ13C值主要由土壤有機(jī)質(zhì)的來(lái)源植物所決定,因此在土壤不同粒徑組分中,與砂粒結(jié)合的有機(jī)質(zhì)的δ13C值是最有可能代表輸入到土壤中地表生物量的有機(jī)碳的平均結(jié)果[28],這說(shuō)明輸入到兩剖面有機(jī)碳的δ13C不一致。

    2.5土壤中元素地球化學(xué)循環(huán)及有機(jī)碳的更替

    SOC及其δ13C受很多因素影響,植被輸入及微生物對(duì)有機(jī)質(zhì)的分解和有機(jī)質(zhì)淋溶和流失等過(guò)程是影響土壤有機(jī)碳循環(huán)的重要因素。土壤質(zhì)地也直接或間接影響著土壤中有機(jī)碳的分解和轉(zhuǎn)化方向。研究土壤剖面的元素分布特征可以從土壤性質(zhì)的角度發(fā)掘影響土壤中δ13C和SOC的垂直分布特征的決定因素。土壤中C含量和元素含量相關(guān),有研究發(fā)現(xiàn)0—50 cm土壤剖面P,Ca,Na和A1等元素都與C含量有很好的相關(guān)性,但是基巖不同,相關(guān)性也不相同[4]。植被的根系等條件也決定著植物從土壤深層吸收元素的能力[29],這些因素也決定著土壤剖面的元素分布特征。表2—表5列出了剖面中SOC含量、δ13C、pH值和粘粒含量及其與各元素質(zhì)量遷移系數(shù)的Pearson相關(guān)分析結(jié)果。

    由表2,剖面X中δ13C與SOC含量相關(guān)(r=0.653),而δ13C與pH值之間呈低度負(fù)相關(guān)。δ13C與所測(cè)得的大多數(shù)常、微量元素(Fe、Li、Na、K、V、Mg、Ba)呈中度負(fù)相關(guān),與Co存在低度正相關(guān);SOC與Fe、K、Li呈低度負(fù)相關(guān),與粘粒含量則存在低度正相關(guān)關(guān)系;土壤pH值與Ba和Fe呈中度正相關(guān),與Mg、Ca和K呈低度正相關(guān),與Co則存在低度負(fù)相關(guān)關(guān)系。以上結(jié)果很好地說(shuō)明了土壤pH和有機(jī)質(zhì)對(duì)成土過(guò)程中元素遷移轉(zhuǎn)化過(guò)程的影響。土壤中的有機(jī)質(zhì)含有多種功能基,這些功能基對(duì)重金屬離子有較強(qiáng)絡(luò)合和富集能力。重金屬離子的存在形態(tài)還受腐殖物質(zhì)的絡(luò)合作用和氧化還原作用的影響,如腐殖物質(zhì)能將V5+還原為V4+、Fe3+還原為Fe2+等。土壤有機(jī)質(zhì)中不同組分的含量與土壤pH值密切相關(guān),低pH值能夠促進(jìn)土壤中酸性水解反應(yīng)的進(jìn)行,從而去除土壤中的蛋白質(zhì)、核酸和多糖等富集13C的組分[30],使得虧損13C的組分在土壤中積累。

    表2 剖面X中δ13C、SOC、pH和粘粒與各元素質(zhì)量遷移系數(shù)的相關(guān)性

    **在 0.01 水平(雙側(cè))上顯著相關(guān);* 在 0.05 水平(雙側(cè))上顯著相關(guān)

    由表3,剖面X中各元素之間呈現(xiàn)良好的相關(guān)性。其中Fe與Li、V和K呈中度正相關(guān);Mg與Ca、Na、Mn和Co呈高度正相關(guān),r均大于0.9;Ca與Mg、Na、Mn和Co呈高度正相關(guān),r均大于0.9;Na與Co、Mg、Ca、Mn呈高度正相關(guān),與K有中度正相關(guān)關(guān)系;K與Fe和Na中度正相關(guān),與Zn中度負(fù)相關(guān)。Mn與Ca、Mg和Na高度正相關(guān)與Co中度正相關(guān);Li與Fe和V中度正相關(guān),與Ba低度正相關(guān),與Co則低度負(fù)相關(guān);Co與Na、Mg、Ca高度正相關(guān),與Mn中度相關(guān),與Li低度負(fù)相關(guān);Ba與Li低度正相關(guān);V與Zn、Fe和Li中度正相關(guān);Zn與K低度負(fù)相關(guān)。元素之間良好的相關(guān)性,說(shuō)明了在成土過(guò)程中各元素之間的協(xié)變關(guān)系。

    表3 剖面X各元素質(zhì)量遷移系數(shù)的相關(guān)性

    由表4,剖面Y中有機(jī)碳的δ13C值與SOC含量呈中度正相關(guān)關(guān)系;土壤的pH值與所測(cè)得的大多數(shù)常、微量元素呈高度正相關(guān)(Mn、Na、Co、Ca、Mg、K)或中度正相關(guān)(V、Fe)關(guān)系。而粘粒含量與各元素質(zhì)量遷移系數(shù)之間不相關(guān)或沒(méi)有線性相關(guān)關(guān)系。說(shuō)明在剖面Y中,元素的遷移轉(zhuǎn)化主要受pH的影響。

    表4 剖面Y中δ13C、SOC、pH和粘粒與各元素質(zhì)量遷移系數(shù)的相關(guān)性

    從表5可以看出剖面Y中各元素之間質(zhì)量遷移系數(shù)的相關(guān)性幾乎均達(dá)到了中度相關(guān)以上。其中Fe尤為突出,與所有被測(cè)元素的遷移系數(shù)均相關(guān)。Ba的遷移系數(shù)除了與Fe、V和Zn呈中度負(fù)相關(guān)外,與其它元素均不相關(guān)或沒(méi)有線性關(guān)系。

    同一氣候條件下,土壤質(zhì)地的不同可能導(dǎo)致兩個(gè)剖面的持水能力、養(yǎng)分狀況、淀積作用的強(qiáng)度、土壤總氧氣的擴(kuò)散能力、土壤穿透阻力、凈初級(jí)生產(chǎn)力等方面的差異,這些因素中的一些將影響輸入到土壤中的有機(jī)質(zhì)的性質(zhì)和數(shù)量,另一些將影響土壤中的物理和化學(xué)過(guò)程[31],從而影響到土壤中元素的遷移富集,從以上的討論可以看出土壤酸堿度對(duì)土壤元素的遷移分布影響較大,因此兩個(gè)剖面粘粒分布的差異及pH值的不同也許是剖面Y中元素之間的相關(guān)性明顯高于剖面X的原因。

    表5 剖面Y各元素質(zhì)量遷移系數(shù)的相關(guān)性

    3結(jié)論

    通過(guò)對(duì)遵義地區(qū)紅色風(fēng)化殼的土壤屬性、元素地球化學(xué)特征和有機(jī)質(zhì)穩(wěn)定碳同位素的分析,得出如下結(jié)論:

    (1)新蒲紅色風(fēng)化殼的粒度分析結(jié)果表明,其成土母巖在垂向上是不均一的,而且當(dāng)?shù)氐乃牡刭|(zhì)條件對(duì)土粒在剖面中的遷移分布有一定的影響。

    (2)總體上兩剖面中有機(jī)質(zhì)的δ13C隨剖面深度的增加均減小,而且δ13C值在各粒級(jí)不同組分中的變化幅度較小。

    (3)對(duì)剖面中常、微量元素遷移分布的研究表明,草地剖面中,鐵、鉀、鋇、釩和鋅表現(xiàn)為一定程度的富集,錳在巖粉處,鈷和鈉在風(fēng)化原巖處富集,而鎂、鈣則幾乎完全遷出土體;農(nóng)耕地剖面中,釩和鋅從原巖到巖土界面是富集的,而巖土界面以上富集趨緩,到表層表現(xiàn)為較明顯的虧損。鐵在巖土界面,鉀和鈷在風(fēng)化原巖處富集,鋇在表層微弱富集,其它元素均虧損。

    (4)通過(guò)對(duì)剖面中δ13C、SOC、pH值和粘粒特征及其與元素垂直分布關(guān)系的研究,發(fā)現(xiàn)草地剖面中δ13C、SOC和pH與元素遷移分布相關(guān)性顯著;而在農(nóng)耕地剖面中各元素在剖面中遷移分布主要與pH值呈顯著正相關(guān),同時(shí)各元素之間也呈現(xiàn)高度相關(guān)性。

    致謝:中國(guó)科學(xué)院地球化學(xué)研究所張風(fēng)雷碩士在野外采樣和同位素測(cè)試給予幫助,特此致謝。

    參考文獻(xiàn)(References):

    [1]王世杰. 喀斯特石漠化——中國(guó)西南最嚴(yán)重的生態(tài)地質(zhì)環(huán)境問(wèn)題. 礦物巖石地球化學(xué)通報(bào), 2003, 22(2): 120- 126.

    [2]Xie X L, Sun B, Zhou H Z, Li A B. Soil organic carbon storage in China. Pedosphere, 2004, 14(4): 491- 500.

    [3]劉叢強(qiáng). 生物地球化學(xué)過(guò)程與地表物質(zhì)循環(huán)——西南喀斯特土壤-植被系統(tǒng)生源要素循環(huán). 北京: 科學(xué)出版社, 2009: 214- 214.

    [4]沈亞婷. EDXRF測(cè)定土壤元素含量及其在有機(jī)碳垂直分布特征研究中的應(yīng)用. 光譜學(xué)與光譜分析, 2012, 32(11): 3117- 3122.

    [6]王楠, 姚佳佳, 高彥征, 劉娟. 黃棕壤中不同粒徑組分的提取分級(jí)與表征. 中國(guó)環(huán)境科學(xué), 2012, 32(12): 2253- 2260.

    [7]朱書(shū)法, 劉叢強(qiáng), 陶發(fā)祥.13C方法在土壤有機(jī)質(zhì)研究中的應(yīng)用. 土壤學(xué)報(bào), 2005, 42(3): 495- 503.

    [8]Chen Q Q, Shen C D, Sun Y M, Peng S L, Li Z A, Jiang M T. Spatial and temporal distribution of carbon isotopes in soil organic matter at the Dinghushan Biosphere Reserve, South China. Plant and Soil, 2005, 273(1/2): 115- 128.

    [9]Loveland P, Webb J. Is there a critical level of organic matter in the agricultural soils of temperate regions: A review. Soil and Tillage Research, 2003, 70(1): 1- 18.

    [10]黃冠華, 詹衛(wèi)華. 土壤顆粒的分形特征及其應(yīng)用. 土壤學(xué)報(bào), 2002, 39(4): 490- 497.

    [11]張冬青, 何騰兵, 林昌虎, 何佳芳. 貴州烏江流域喀斯特土壤物理性質(zhì)的研究. 貴州科學(xué), 2006, 24(2): 39- 43.

    [12]馮志剛, 王世杰, 周德全, 孫承興, 劉秀明. 黔北新蒲紅色風(fēng)化殼剖面的粒度分布特征——兼論巖溶區(qū)上覆土層的物質(zhì)來(lái)源. 地質(zhì)科學(xué), 2004, 39(3): 440- 451.

    [13]馮志剛, 王世杰, 孫承興. 引起紅土表層硅鋁比值增大原因的可能性探討. 地質(zhì)地球化學(xué), 2002, 30(4): 7- 14.

    [14]季宏兵, 歐陽(yáng)自遠(yuǎn), 王世杰, 周德全. 白云巖風(fēng)化剖面的元素地球化學(xué)特征及其對(duì)上陸殼平均化學(xué)組成的意義——以黔北新蒲剖面為例. 中國(guó)科學(xué)(D輯), 1999, 29(6): 504- 513.

    [15]周德全, 王世杰, 劉秀明. 石灰土(碳酸鹽巖風(fēng)化殼)形成地球化學(xué)過(guò)程研究. 地球與環(huán)境, 2005, 33(2): 31- 38.

    [16]Nesbitt H W, Markovics G. Weathering of granodioritic crust, long-term storage of elements in weathering profiles, and petrogenesis of siliciclastic sediments. Geochimica et Cosmochimica Acta, 1997, 61(8): 1653- 1670.

    [17]White A F, Blurn A E, Schulz M S, Vivit D V, Stonestrom D A, Larsen M, Murphy S F, Eberl D. Chemical weathering in a tropical watershed, Luquillo Mountains, Puerto Rico: I. Long-term versus short-term weathering fluxes. Geochimica et Cosmochimica Acta, 1998, 62(2): 209- 226.

    [18]黃昌勇. 土壤學(xué). 北京: 中國(guó)農(nóng)業(yè)出版社, 2000: 209- 212.

    [19]Natelhoffer K J, Fry B. Controls on natural nitrogen-15 and carbon-13 abundances in forest soil organic matter. Soil Science Society of America Journal, 1988, 52(6): 1633- 1640.

    [20]Connin S L, Feng X, Virginia R A. Isotopic discrimination during long-term decomposition in an arid land ecosystem. Soil Biology and Biogeochemistry, 2001, 33(1): 41- 51.

    [21]王國(guó)安, 韓家懋, 周立平, 熊小剛, 譚明, 吳振海, 彭雋. 中國(guó)北方黃土區(qū)C4植物穩(wěn)定碳同位素組成的研究. 中國(guó)科學(xué)(D輯), 2005, 35(12): 1174- 1179.

    [22]樸河春, 劉啟明, 余登利, 郭景恒, 冉景丞. 用天然13C豐度法評(píng)估貴州茂蘭喀斯特森林區(qū)玉米地土壤中有機(jī)碳的來(lái)源. 生態(tài)學(xué)報(bào), 2001, 21(3): 434- 439.

    [23]李龍波, 涂成龍, 趙志琦, 崔麗峰, 劉文景. 黃土高原不同植被覆蓋下土壤有機(jī)碳的分布特征及其同位素組成研究. 地球與環(huán)境, 2011, 39(4): 441- 449.

    [24]陳慶強(qiáng), 沈承德, 孫彥敏, 彭少麟, 易惟熙, 李志安, 姜漫濤. 鼎湖山土壤有機(jī)質(zhì)深度分布的剖面演化機(jī)制. 土壤學(xué)報(bào), 2005, 42(1): 1- 8.

    [25]Boutton T W. Stable carbon isotope ratios of soil organic matter and their use as indicators of vegetation and climate change // Boutton T W, Yamasaki S I, eds. Mass Spectrometry of Soils. New York: Marcel Dekker, 1996: 47- 82.

    [26]Guggenberger G, Christensen B T, Zech W. Land-use effects on the composition of organic matter in particle-size separates of soil: Ⅰ. Lignin and carbohydrate signature. European Journal of Soil Science, 1994, 45(4): 449- 458.

    [27]Krull E S, Skjemstad J O, Burrows W H, Bray S G, Wynn J G, Bol R, Spouncer L, Harms B. Recent vegetation changes in central Queensland, Australia: Evidence fromδ13C and14C analyses of soil organic matter. Geoderma, 2005, 126(3/4): 241- 259.

    [28]Bird M, Kracht O, Derrien D, Zhou Y. The effect of soil texture and roots on the stable carbon isotope composition of soil organic carbon. Australian Journal of Soil Research, 2003, 41(1): 77- 94.

    [29]Berger T W, Swoboda S, Prohaska T, Glatzel G. The role of calcium uptake from deep soils for spruce (Piceaabies) and beech (Fagussylvatica). Forest Ecology and Management, 2006, 229(1/3): 234- 246.

    [30]Krull E S, Skjemstad J O.δ13C andδ15N profiles in14C-dated Oxisol and Vertisols as a function of soil chemistry and mineralogy. Geoderma, 2003, 112(1/2): 1- 29.

    [31]Jobbagy E G and Jackson R B. The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecological Application, 2000, 10: 423- 436.

    Grain size distribution, element migration, and stable carbon isotope characteristics of dolomite weathering profiles in Xinpu, North of Guizhou Province

    GAO Jie1, LI Rui1, LI Jinjin1, ZHANG Li1, JI Hongbing1,2,*

    1CapitalNormalUniversity,Beijing100048,China2InstituteofGeochemistry,ChineseAcademyofScience,Guiyang550002,China

    Abstract:The relation of the vertical distribution and cycle characteristics of soil organic carbon (SOC) with the elements in the soil is an important area of research on the soil geochemical circulation and climate change. At the same time, the distribution characteristics of soil elements can also affect the cycle rate of soil organic carbon. To identify the parameters that may influence the turnover rate of soil organic carbon in karst terrain and their relations with the geochemical behavior of chemical elements, we selected two dolomite weathering profiles for analysis (the surface layer is grass and new cultivated land), which are located in Zunyi, Guizhou Province. On the basis of particle size distribution, of the δ13C value of SOC, and the quality migration coefficient of the elements, we discuss the trends of soil organic carbon and δ13C values with the deepening of the red weathering crust in this area and the geochemical characteristics of the elements. Pearson′s correlation coefficients were calculated to quantify the relations. Our results indicated that SOC concentrations in the grass profile decreased with the deepening of soil layers, but in the farmland profile, SOC content of the surface layer increased rapidly at the depth of 25 cm, and then decreased with deepening of the profile. The pH value of the surface layer was acidic in both profiles, and the acidity of the grass profile was stronger than that in the farmland profile. The pH value increased with the deepening of soil layers, but the rising trend of the grass profile was more obvious: Soil grain size fraction showed that the silt fraction (5—50 μm) made up 50%—80% of the soil, whereas the sand fraction was the least prevalent. The share of each particle size varied with the deepening of the profile. By means of mass balance approximation and by selecting Ti as an inactive reference element, we obtained a quality migration coefficient of each relevant element. The results showed that almost all elements in the profile showed a certain degree of loss or enrichment. In the grass profile, Fe, K, Ba, V, and Zn showed some enrichment; Na and Co were enriched in the weathering protolith; Mn was enriched in the dolomite powder. In contrast, Mg and Ca were almost absent in the profile. In the farmland profile, V and Zn showed certain enrichment. Fe and K were enriched at the interface of the soil and dolomite; Co was enriched in the weathering protolith, and Ba showed weak enrichment at the surface. Analysis of the δ13C value of SOC showed that this metric was between 26.4‰ and 23.6‰. The δ13C value of both the original soil samples and the silt fraction samples decreased with increasing profile depth, and there was no obvious difference among soil grain fractions. In the grass profile, δ13C and pH significantly correlated with most elements; SOC correlated negatively with Fe, K, and Li (P < 0.05); the migration of elements in the farmland profile was mainly affected by pH. Among all of the elements, there was a significant relation in both soil profiles.Thus, it is possible that the basement carbonate in the vertical direction is inhomogeneous in the Xinpu region, and the local hydrogeological conditions have some influence on the distribution of soil particles. The distribution and migration of elements in the profile are related to soil pH, SOC, and clay content in a soil- forming process; there is also a good covariant relation among the elements.

    Key Words:soil organic carbon; grain size distribution; quality migration coefficient; stable carbon isotope

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(41473122, 41073096); 國(guó)家重點(diǎn)研究項(xiàng)目(2013CB956702); 中國(guó)科學(xué)院“百人計(jì)劃”

    收稿日期:2014- 07- 20; 網(wǎng)絡(luò)出版日期:2015- 07- 22

    DOI:10.5846/stxb201407201477

    *通訊作者Corresponding author.E-mail: jih_0000@126.com

    高杰,李銳,李今今,張莉,季宏兵.白云巖風(fēng)化剖面的粒度分布、元素遷移及碳同位素特征——以黔北新蒲剖面為例.生態(tài)學(xué)報(bào),2016,36(5):1409- 1420.

    Gao J, Li R, Li J J, Zhang L, Ji H B.Grain size distribution, element migration, and stable carbon isotope characteristics of dolomite weathering profiles in Xinpu, North of Guizhou Province.Acta Ecologica Sinica,2016,36(5):1409- 1420.

    猜你喜歡
    有機(jī)碳
    貴州黃壤性水稻土不同粒徑有機(jī)碳之間的礦化差異
    長(zhǎng)期施肥下黃壤稻田土壤有機(jī)碳和全氮的演變特征
    差減法測(cè)定土壤中生物炭含量
    不同土地利用方式對(duì)克氏針茅草原表層土壤有機(jī)碳的影響
    綠色科技(2017年14期)2017-08-22 08:36:47
    江漢區(qū)五峰組—龍馬溪組頁(yè)巖氣保存條件分析
    毛烏素沙地典型湖濱帶濕地有機(jī)碳、氮素空間分布特征
    糧草間作培肥處理對(duì)黑土土壤有機(jī)碳含量的影響
    河流水體有機(jī)碳研究綜述
    土壤有機(jī)碳含量測(cè)定方法評(píng)述及最新研究進(jìn)展
    生草栽培技術(shù)對(duì)果園土壤團(tuán)聚體及其有機(jī)碳分布的分析
    国产精品麻豆人妻色哟哟久久| 大片电影免费在线观看免费| 亚洲av国产av综合av卡| 男女午夜视频在线观看| 精品卡一卡二卡四卡免费| 欧美日韩福利视频一区二区| 高清视频免费观看一区二区| 黄色 视频免费看| 亚洲国产欧美网| 亚洲成a人片在线一区二区| av天堂在线播放| 婷婷丁香在线五月| 一级片免费观看大全| 国产不卡av网站在线观看| 欧美亚洲日本最大视频资源| 国产成人av激情在线播放| 国产精品成人在线| 亚洲精品美女久久久久99蜜臀| 日韩欧美国产一区二区入口| 国产在视频线精品| 首页视频小说图片口味搜索| 久久久久久久久久久久大奶| 国产欧美日韩一区二区三区在线| 少妇的丰满在线观看| 高清视频免费观看一区二区| 乱人伦中国视频| 9热在线视频观看99| 欧美另类亚洲清纯唯美| 91国产中文字幕| 国产免费现黄频在线看| 在线观看一区二区三区激情| 又黄又粗又硬又大视频| av一本久久久久| 国产国语露脸激情在线看| 免费女性裸体啪啪无遮挡网站| 少妇精品久久久久久久| 黑人巨大精品欧美一区二区mp4| 国产精品一区二区在线不卡| 久久九九热精品免费| 成年人免费黄色播放视频| 亚洲欧美色中文字幕在线| 女同久久另类99精品国产91| 久久久久久久久免费视频了| 精品久久久久久电影网| 久久精品国产亚洲av高清一级| 精品一区二区三卡| 亚洲精品成人av观看孕妇| 日日爽夜夜爽网站| 国产97色在线日韩免费| 大型av网站在线播放| 欧美黄色片欧美黄色片| 成人国产一区最新在线观看| 午夜福利影视在线免费观看| 欧美日韩亚洲国产一区二区在线观看 | 日韩一卡2卡3卡4卡2021年| 黑丝袜美女国产一区| 亚洲av成人不卡在线观看播放网| 国产精品国产高清国产av | av又黄又爽大尺度在线免费看| 亚洲精华国产精华精| 91成年电影在线观看| 少妇被粗大的猛进出69影院| 老汉色∧v一级毛片| 久久精品国产综合久久久| 久久ye,这里只有精品| 久久久精品国产亚洲av高清涩受| 国产免费福利视频在线观看| 十八禁网站网址无遮挡| 免费在线观看影片大全网站| 又黄又粗又硬又大视频| 一级毛片女人18水好多| 亚洲avbb在线观看| 后天国语完整版免费观看| 欧美激情高清一区二区三区| 亚洲久久久国产精品| 亚洲第一av免费看| 亚洲精品在线美女| 热99国产精品久久久久久7| 50天的宝宝边吃奶边哭怎么回事| 色婷婷av一区二区三区视频| 免费人妻精品一区二区三区视频| 久久人妻福利社区极品人妻图片| 久久中文看片网| 女人高潮潮喷娇喘18禁视频| 日本黄色视频三级网站网址 | 久久国产精品人妻蜜桃| 狂野欧美激情性xxxx| 亚洲色图 男人天堂 中文字幕| 男女之事视频高清在线观看| 男女床上黄色一级片免费看| 国产日韩欧美在线精品| 日韩中文字幕视频在线看片| 最近最新免费中文字幕在线| 亚洲av日韩精品久久久久久密| 午夜福利乱码中文字幕| 日韩免费av在线播放| 一二三四在线观看免费中文在| 精品高清国产在线一区| √禁漫天堂资源中文www| 亚洲伊人久久精品综合| 久久精品成人免费网站| 国产精品自产拍在线观看55亚洲 | 亚洲伊人色综图| 亚洲三区欧美一区| 一级a爱视频在线免费观看| 国产在视频线精品| 看免费av毛片| 国产日韩一区二区三区精品不卡| 精品久久蜜臀av无| 午夜精品久久久久久毛片777| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品久久久久5区| 亚洲熟妇熟女久久| 久久中文看片网| 国产免费av片在线观看野外av| 亚洲精品久久午夜乱码| 手机成人av网站| 涩涩av久久男人的天堂| 一本一本久久a久久精品综合妖精| 青草久久国产| 成人影院久久| 一二三四在线观看免费中文在| 纯流量卡能插随身wifi吗| 一进一出好大好爽视频| 99国产精品一区二区蜜桃av | 亚洲va日本ⅴa欧美va伊人久久| 操出白浆在线播放| 久久久久久免费高清国产稀缺| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品久久午夜乱码| av天堂久久9| 国产精品一区二区免费欧美| 啦啦啦中文免费视频观看日本| 免费一级毛片在线播放高清视频 | 黑人欧美特级aaaaaa片| 久久 成人 亚洲| 变态另类成人亚洲欧美熟女 | 午夜福利乱码中文字幕| 极品少妇高潮喷水抽搐| 少妇精品久久久久久久| 欧美日韩亚洲国产一区二区在线观看 | 国产一区二区在线观看av| 啦啦啦视频在线资源免费观看| 亚洲全国av大片| 岛国毛片在线播放| 日本撒尿小便嘘嘘汇集6| 一本久久精品| 国产精品久久久人人做人人爽| 99riav亚洲国产免费| 制服人妻中文乱码| 久久久久久人人人人人| 国产亚洲精品久久久久5区| 大码成人一级视频| 亚洲熟女毛片儿| 国产在线视频一区二区| 在线av久久热| 女同久久另类99精品国产91| videosex国产| 美女扒开内裤让男人捅视频| 每晚都被弄得嗷嗷叫到高潮| 天堂中文最新版在线下载| 国产成人啪精品午夜网站| 欧美日韩精品网址| 99精品欧美一区二区三区四区| 亚洲一卡2卡3卡4卡5卡精品中文| 91九色精品人成在线观看| 免费少妇av软件| cao死你这个sao货| 高潮久久久久久久久久久不卡| 午夜视频精品福利| 国产伦人伦偷精品视频| 最新美女视频免费是黄的| 999久久久精品免费观看国产| 日本黄色日本黄色录像| 久久久久视频综合| 久久久精品区二区三区| 国产男女超爽视频在线观看| 久久精品国产99精品国产亚洲性色 | 国产成人欧美| 久久中文字幕一级| 国产日韩一区二区三区精品不卡| 黑人巨大精品欧美一区二区蜜桃| 交换朋友夫妻互换小说| 免费在线观看日本一区| 国产在线观看jvid| 欧美精品一区二区免费开放| 香蕉久久夜色| 国产不卡av网站在线观看| 黄片大片在线免费观看| 国产在视频线精品| 久久中文看片网| 麻豆国产av国片精品| 色综合婷婷激情| 日韩成人在线观看一区二区三区| 久久久久国内视频| 久久精品国产综合久久久| av在线播放免费不卡| 1024视频免费在线观看| 亚洲av成人一区二区三| 国产精品国产av在线观看| 欧美日韩视频精品一区| av不卡在线播放| 80岁老熟妇乱子伦牲交| 视频区图区小说| 97在线人人人人妻| 免费女性裸体啪啪无遮挡网站| 免费久久久久久久精品成人欧美视频| 欧美日韩亚洲高清精品| 大型av网站在线播放| 国产高清videossex| 妹子高潮喷水视频| 欧美激情高清一区二区三区| 国产在线精品亚洲第一网站| 丁香六月天网| 成年人黄色毛片网站| 亚洲精品av麻豆狂野| 十八禁高潮呻吟视频| 69av精品久久久久久 | 亚洲性夜色夜夜综合| 午夜成年电影在线免费观看| 成人国产一区最新在线观看| 1024视频免费在线观看| 成人永久免费在线观看视频 | 精品少妇久久久久久888优播| 精品人妻1区二区| 满18在线观看网站| 免费一级毛片在线播放高清视频 | 色尼玛亚洲综合影院| 亚洲国产欧美日韩在线播放| 国产成人精品无人区| 成人影院久久| 波多野结衣av一区二区av| 纵有疾风起免费观看全集完整版| 日日夜夜操网爽| 久久国产精品男人的天堂亚洲| 在线看a的网站| 99riav亚洲国产免费| 99热国产这里只有精品6| 免费不卡黄色视频| 2018国产大陆天天弄谢| 国产97色在线日韩免费| 国产成人免费观看mmmm| 国产精品成人在线| 99香蕉大伊视频| 超碰97精品在线观看| 18禁观看日本| 久久亚洲精品不卡| 成人亚洲精品一区在线观看| 亚洲欧美一区二区三区黑人| 中文字幕人妻熟女乱码| 免费观看a级毛片全部| 多毛熟女@视频| 欧美+亚洲+日韩+国产| 精品福利永久在线观看| 午夜福利影视在线免费观看| 国产在线视频一区二区| av天堂久久9| 无遮挡黄片免费观看| 久久午夜亚洲精品久久| 欧美性长视频在线观看| 亚洲av日韩精品久久久久久密| 午夜久久久在线观看| 日韩免费av在线播放| 69av精品久久久久久 | 亚洲av第一区精品v没综合| 精品国产乱子伦一区二区三区| 天天影视国产精品| 无人区码免费观看不卡 | 欧美日韩精品网址| 久久精品人人爽人人爽视色| 午夜视频精品福利| 国产精品影院久久| 日本av手机在线免费观看| xxxhd国产人妻xxx| 国产成+人综合+亚洲专区| 黑丝袜美女国产一区| 女人精品久久久久毛片| 国产成人欧美| 久久精品亚洲av国产电影网| 亚洲专区国产一区二区| 亚洲国产精品一区二区三区在线| 国产成人精品久久二区二区免费| 99在线人妻在线中文字幕 | 国产99久久九九免费精品| 国产不卡一卡二| 国产国语露脸激情在线看| 侵犯人妻中文字幕一二三四区| 狠狠婷婷综合久久久久久88av| 怎么达到女性高潮| 美女高潮喷水抽搐中文字幕| 丰满人妻熟妇乱又伦精品不卡| 国产男靠女视频免费网站| 如日韩欧美国产精品一区二区三区| 女性被躁到高潮视频| 久久精品人人爽人人爽视色| 啦啦啦在线免费观看视频4| 80岁老熟妇乱子伦牲交| 亚洲欧洲日产国产| av欧美777| 性少妇av在线| 好男人电影高清在线观看| 汤姆久久久久久久影院中文字幕| 9热在线视频观看99| 久久中文字幕人妻熟女| 中文字幕高清在线视频| 黄色视频不卡| 免费在线观看日本一区| 天堂动漫精品| 亚洲欧美日韩高清在线视频 | 久久 成人 亚洲| 久久av网站| 熟女少妇亚洲综合色aaa.| 日本五十路高清| 波多野结衣一区麻豆| 变态另类成人亚洲欧美熟女 | 亚洲精品中文字幕在线视频| 国产黄频视频在线观看| 亚洲欧美激情在线| 1024视频免费在线观看| 侵犯人妻中文字幕一二三四区| 中文字幕人妻丝袜制服| 咕卡用的链子| 国产三级黄色录像| 超碰成人久久| 热99久久久久精品小说推荐| 婷婷成人精品国产| 日韩 欧美 亚洲 中文字幕| 91精品三级在线观看| 久久久久久久精品吃奶| 久久这里只有精品19| 国产aⅴ精品一区二区三区波| 看免费av毛片| 欧美成人午夜精品| 国产xxxxx性猛交| 黄色视频不卡| 国产成人影院久久av| 91麻豆av在线| 久久国产亚洲av麻豆专区| 亚洲欧美精品综合一区二区三区| 亚洲国产毛片av蜜桃av| 女性生殖器流出的白浆| 久久久久久久久久久久大奶| 久久这里只有精品19| 啦啦啦视频在线资源免费观看| 成年版毛片免费区| 国产亚洲精品第一综合不卡| 成人亚洲精品一区在线观看| 丝袜人妻中文字幕| 少妇的丰满在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲午夜精品一区,二区,三区| 男人舔女人的私密视频| 精品少妇黑人巨大在线播放| 国产福利在线免费观看视频| 两人在一起打扑克的视频| 亚洲成人手机| 国产视频一区二区在线看| 美国免费a级毛片| 国产欧美日韩一区二区三区在线| 精品视频人人做人人爽| 肉色欧美久久久久久久蜜桃| 一本—道久久a久久精品蜜桃钙片| 亚洲色图综合在线观看| 女同久久另类99精品国产91| 午夜福利一区二区在线看| 精品国产乱码久久久久久男人| 一进一出抽搐动态| 午夜91福利影院| 美女视频免费永久观看网站| 欧美变态另类bdsm刘玥| 久久久久国内视频| 人妻久久中文字幕网| 精品国产国语对白av| 一个人免费看片子| 欧美在线一区亚洲| 成人国语在线视频| 女人爽到高潮嗷嗷叫在线视频| 国产在线一区二区三区精| 女性生殖器流出的白浆| 18禁黄网站禁片午夜丰满| 国精品久久久久久国模美| 久久久久久久国产电影| 午夜福利免费观看在线| 欧美黄色淫秽网站| 9191精品国产免费久久| 嫩草影视91久久| 一二三四社区在线视频社区8| 18禁黄网站禁片午夜丰满| 黄色怎么调成土黄色| 国产日韩欧美视频二区| 欧美日本中文国产一区发布| 久久久国产成人免费| 99国产精品99久久久久| 亚洲专区中文字幕在线| 亚洲熟妇熟女久久| 日韩一卡2卡3卡4卡2021年| 亚洲色图综合在线观看| 老熟妇仑乱视频hdxx| 91成年电影在线观看| 中文字幕人妻熟女乱码| 亚洲国产欧美网| 丁香六月天网| 久久中文字幕人妻熟女| 人妻久久中文字幕网| 欧美日韩av久久| 男女无遮挡免费网站观看| 母亲3免费完整高清在线观看| 熟女少妇亚洲综合色aaa.| 女性被躁到高潮视频| 人人妻人人爽人人添夜夜欢视频| 丝袜人妻中文字幕| 国产精品一区二区在线不卡| 亚洲专区中文字幕在线| 欧美激情极品国产一区二区三区| 亚洲国产av影院在线观看| 久久精品亚洲精品国产色婷小说| 精品人妻1区二区| 久久久久久久精品吃奶| 黄色丝袜av网址大全| 精品少妇一区二区三区视频日本电影| 一本一本久久a久久精品综合妖精| 亚洲熟女精品中文字幕| 精品国产乱子伦一区二区三区| 亚洲全国av大片| 午夜福利在线免费观看网站| 99久久精品国产亚洲精品| 在线观看66精品国产| 国产不卡av网站在线观看| 两性夫妻黄色片| 50天的宝宝边吃奶边哭怎么回事| 在线观看免费高清a一片| 少妇精品久久久久久久| 亚洲,欧美精品.| 久久久久久久久久久久大奶| 脱女人内裤的视频| 国产深夜福利视频在线观看| 国产精品自产拍在线观看55亚洲 | 久久午夜亚洲精品久久| 50天的宝宝边吃奶边哭怎么回事| 99国产极品粉嫩在线观看| 日韩三级视频一区二区三区| 99精品久久久久人妻精品| 午夜福利影视在线免费观看| 精品视频人人做人人爽| av网站免费在线观看视频| 国产一区二区 视频在线| 亚洲专区中文字幕在线| 老司机影院毛片| 中文字幕最新亚洲高清| 一区二区三区激情视频| 黄色视频不卡| av网站在线播放免费| 亚洲色图av天堂| 成年人免费黄色播放视频| 成人国产av品久久久| 波多野结衣av一区二区av| 午夜福利在线观看吧| 91老司机精品| 日本欧美视频一区| 久久天躁狠狠躁夜夜2o2o| 啦啦啦视频在线资源免费观看| 亚洲欧美色中文字幕在线| 国产精品久久久久久人妻精品电影 | 成人三级做爰电影| 国产亚洲一区二区精品| 久久久久精品人妻al黑| 国产单亲对白刺激| 一二三四在线观看免费中文在| 美女高潮喷水抽搐中文字幕| 最新在线观看一区二区三区| 在线 av 中文字幕| 欧美黄色片欧美黄色片| 国产单亲对白刺激| 黄色成人免费大全| 黄色视频,在线免费观看| 熟女少妇亚洲综合色aaa.| 精品久久久久久久毛片微露脸| 正在播放国产对白刺激| 黄色 视频免费看| 一夜夜www| 国产精品一区二区在线不卡| 天天躁夜夜躁狠狠躁躁| 最近最新免费中文字幕在线| 精品第一国产精品| 国产又爽黄色视频| 国产av又大| 女人爽到高潮嗷嗷叫在线视频| 91精品三级在线观看| 91麻豆av在线| 国产无遮挡羞羞视频在线观看| 亚洲黑人精品在线| 色视频在线一区二区三区| 国产xxxxx性猛交| 成人黄色视频免费在线看| 免费在线观看视频国产中文字幕亚洲| 亚洲精品一二三| 久久久久久亚洲精品国产蜜桃av| 老司机在亚洲福利影院| 国产av精品麻豆| 国产成人啪精品午夜网站| 日韩欧美国产一区二区入口| 亚洲精华国产精华精| 91精品国产国语对白视频| 亚洲精品一二三| 日本wwww免费看| 妹子高潮喷水视频| 91av网站免费观看| 国产日韩一区二区三区精品不卡| 久久国产精品男人的天堂亚洲| 国产精品成人在线| 欧美日韩av久久| 19禁男女啪啪无遮挡网站| 亚洲国产欧美在线一区| www.自偷自拍.com| 电影成人av| www日本在线高清视频| 欧美黄色片欧美黄色片| 2018国产大陆天天弄谢| 天天躁夜夜躁狠狠躁躁| 一级片免费观看大全| 日韩免费高清中文字幕av| 午夜福利一区二区在线看| 一区二区av电影网| av免费在线观看网站| 国产午夜精品久久久久久| 久久热在线av| 日韩成人在线观看一区二区三区| bbb黄色大片| 欧美日韩精品网址| 在线天堂中文资源库| 久久人人爽av亚洲精品天堂| 91成人精品电影| 亚洲熟女精品中文字幕| 国产精品久久久久久精品电影小说| 日本精品一区二区三区蜜桃| 伦理电影免费视频| 男女午夜视频在线观看| 最黄视频免费看| 97人妻天天添夜夜摸| 精品国产乱码久久久久久男人| 国精品久久久久久国模美| 精品人妻在线不人妻| 少妇的丰满在线观看| 悠悠久久av| 中文欧美无线码| 久久久久网色| 少妇裸体淫交视频免费看高清 | 黄频高清免费视频| 一级毛片女人18水好多| 久久精品aⅴ一区二区三区四区| 精品欧美一区二区三区在线| 婷婷成人精品国产| 国产精品成人在线| 亚洲午夜精品一区,二区,三区| 俄罗斯特黄特色一大片| 中文字幕av电影在线播放| 考比视频在线观看| 国产精品欧美亚洲77777| 国产亚洲精品第一综合不卡| 日韩一卡2卡3卡4卡2021年| 女同久久另类99精品国产91| 亚洲七黄色美女视频| 精品国产国语对白av| 精品少妇黑人巨大在线播放| 国产在线精品亚洲第一网站| 国产主播在线观看一区二区| 精品午夜福利视频在线观看一区 | 久热这里只有精品99| 国产精品免费大片| 极品少妇高潮喷水抽搐| 欧美日韩一级在线毛片| 国产精品免费视频内射| 精品免费久久久久久久清纯 | 黄色a级毛片大全视频| 在线观看免费视频网站a站| av网站在线播放免费| 日韩欧美国产一区二区入口| 国产欧美日韩综合在线一区二区| 日本黄色视频三级网站网址 | 国产精品自产拍在线观看55亚洲 | 成人影院久久| 精品视频人人做人人爽| 美女扒开内裤让男人捅视频| 精品亚洲成国产av| 久久精品国产99精品国产亚洲性色 | 一进一出好大好爽视频| 一级毛片电影观看| 丁香欧美五月| 亚洲第一欧美日韩一区二区三区 | 日韩欧美一区二区三区在线观看 | 一本一本久久a久久精品综合妖精| 亚洲精品国产色婷婷电影| avwww免费| 亚洲精品国产一区二区精华液| 色在线成人网| 女人精品久久久久毛片| 午夜免费成人在线视频| 国产欧美日韩综合在线一区二区| 亚洲精品国产区一区二| 亚洲一卡2卡3卡4卡5卡精品中文| 成人特级黄色片久久久久久久 | 国产精品亚洲av一区麻豆| 韩国精品一区二区三区| 国产欧美日韩一区二区精品| 建设人人有责人人尽责人人享有的| 国产av又大| 精品国产亚洲在线| 久9热在线精品视频| 亚洲精品中文字幕在线视频| av福利片在线| 美女视频免费永久观看网站| 最黄视频免费看| 12—13女人毛片做爰片一| 香蕉国产在线看| 久久久久久久大尺度免费视频|