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

    瑪曲高寒退化草甸土壤重金屬污染風險評價

    2024-04-28 05:47:56孫建財白紅桃熊增連張杰雪全小龍鄧得婷史惠蘭
    生態(tài)科學 2024年1期
    關(guān)鍵詞:瑪曲草甸金屬元素

    孫建財, 白紅桃, 熊增連, 張杰雪, 全小龍, 鄧得婷, 史惠蘭,*

    瑪曲高寒退化草甸土壤重金屬污染風險評價

    孫建財1, 白紅桃1, 熊增連3, 張杰雪2, 全小龍2, 鄧得婷1, 史惠蘭1,*

    1. 青海大學生態(tài)環(huán)境工程學院, 西寧 810016 2. 青海大學省部共建三江源生態(tài)和高原農(nóng)牧業(yè)國家重點實驗室, 西寧 810016 3. 青海省自然資源綜合調(diào)查監(jiān)測院, 西寧 810016

    為了分析瑪曲高寒草甸土壤重金屬含量分布及重金屬污染風險,采集瑪曲高寒退化草甸土壤樣品, 分析研究區(qū)表層土壤(0—20 cm)中Cr、Ni、Cu、Zn、As、Cd、Hg、Pb等八種重金屬含量, 描述瑪曲高寒草甸土壤重金屬的空間分布特征, 采用單因子污染指數(shù)、內(nèi)梅羅污染指數(shù)、地累積指數(shù)評價重金屬污染程度并分析重金屬污染潛在生態(tài)風險, 以相關(guān)性分析、主成分分析和聚類分析表征同源性重金屬。結(jié)果表明: (1)研究區(qū)內(nèi)As、Cd的質(zhì)量分數(shù)高于背景值, Ni的質(zhì)量分數(shù)偏低, 變異系數(shù)大小為Cr> Hg> Cd> Cu> As> Ni> Zn> Pb;(2)Cd單因子污染指數(shù)為輕微污染, 地累積指數(shù)反映As屬于未受污染與輕微污染之間, 各重金屬元素均無潛在生態(tài)風險; 以國標和土壤背景值為標準, 內(nèi)梅羅綜合污染指數(shù)表明瑪曲未受到重金屬污染; (3)主成分和相關(guān)性分析表明瑪曲土壤中Ni、Cu、Zn、As、Cd、Pb等元素具有相同或相似的來源, Ni和Pb的來源優(yōu)先順序最大, Cr和Hg含量呈現(xiàn)弱相關(guān)。

    高寒草甸; 重金屬污染; 污染指數(shù); 潛在生態(tài)風險; 主成分分析

    0 前言

    土壤是巖石圈、水圈、大氣圈和生物圈物質(zhì)循環(huán)和能量流動的緊密過渡地帶, 是人類賴以生存和發(fā)展的重要資源之一, 是陸地生態(tài)系統(tǒng)的重要組成部分, 具有重要的生態(tài)意義[1]。重金屬作為典型的累積性污染物, 其污染效果具有持久性、隱蔽性和難降解性, 近年來土壤生態(tài)不斷受到重金屬污染導致土壤穩(wěn)定性降低、生產(chǎn)力下降, 造成地表植被覆蓋率降低, 水土流失、環(huán)境惡化等一系列問題, 并且重金屬易進入食物鏈, 影響農(nóng)作物質(zhì)量和安全, 危害人體健康, 引起了社會廣泛關(guān)注[2]。

    瑪曲縣地處青藏高原東北緣, 土壤以亞高寒草甸土和泥炭土為主, 主要植被類型為亞高寒草甸, 平均海拔3700 m。退化高寒草甸是草地生態(tài)學研究的重點, 草地退化最直接的表現(xiàn)是草地群落類型的轉(zhuǎn)變, 造成雜草入侵, 生物量下降等。由于獨特的地理位置和有限的工農(nóng)業(yè)活動, 瑪曲退化草地重金屬污染輕, 但隨著自然資源的開發(fā)與利用, 青藏高原土壤生態(tài)系統(tǒng)已遭受一定程度的重金屬污染[3,4]。有研究表明青藏高原東北部湖泊汞含量顯著增加[7], 青藏高原路基兩側(cè)受重金屬污染最為嚴重[5,6]。由于重金屬污染的隱蔽性和不可逆性, 預防輸入成為控制重金屬污染的主要手段[8]。本研究對8種重金屬(Cr、Ni、Cu、Zn、As、Cd、Hg、Pb)進行分析, 探討瑪曲高寒退化草甸重金屬污染狀況, 評價重金屬污染潛在生態(tài)風險, 以期為瑪曲高寒退化草甸重金屬污染防控提供科學依據(jù)。

    1 材料和方法

    1.1 研究區(qū)概況

    瑪曲縣地處青藏高原東北緣, 位于甘肅省甘南藏族自治州西南部, 甘、青、川三省交界處, 黃河第一彎曲部, 地處東經(jīng)100°45′45″—102°29′00″, 北緯33°06′30″—34°30′15″之間?,斍h總面積10190 km2, 平均海拔3700 m。該區(qū)域冬春季漫長, 嚴寒多風; 夏秋季短, 多雨潮濕; 年均氣溫1.2 ℃, 年均降雨量 600 mm, 年際降水量呈下降趨勢, 年蒸發(fā)量1352.4 mm, 年均風速2.5 m·s–1, 地表部分被沙土覆蓋, 土壤以亞高寒草甸土和泥炭土為主, 草地類型為亞高寒草甸。

    1.2 樣品采集與分析

    于2017年7月赴甘肅省瑪曲縣歐拉鄉(xiāng), 該鄉(xiāng)為退化草地生態(tài)修復試驗示范地, 主要針對禿斑地采用人工撒播、施肥和輕耙處理, 補播牧草種子包括同德垂穗披堿草、青海草地早熟禾、青海中華羊茅三種, 針對退化草地生態(tài)修復測定重金屬污染水平, 布設16個樣條, 調(diào)查草地植物群落特征, 采集植物生物量, 同時, 取深度為0—20 cm的土壤樣品, 帶回實驗室晾干, 挑出植物根莖等, 將土研磨后過100目篩, 稱取0.1000 g樣品, 加入7 mLHNO3, 1mL HF及0.5 mLHClO4, 進行微波消解[9], 同時做試劑空白; 消解完成后, 將樣品轉(zhuǎn)移至坩堝中, 在電熱板上100 ℃條件下趕酸, 完全冷卻后, 以超純水定容至50 mL, 用ICP-MS[10]測定重金屬含量。

    1.3 重金屬污染評價方法

    1.3.1 單因子指數(shù)法

    單因子污染指數(shù)是指土壤中某一重金屬的含量與土壤環(huán)境質(zhì)量評價標準之比, 計算方式為[11]:

    P=C/S

    式中:P為污染因子的污染指數(shù),C為污染因子的實測濃度,S為污染物的標準值。

    1.3.2 內(nèi)梅羅污染指數(shù)法

    內(nèi)梅羅指數(shù)法涉及單因子污染指數(shù)的平均值和最大值, 可以反映各重金屬對土壤的綜合作用, 計算公式如下[12]:

    綜合=[(平均+Pmax)/2]1/2

    式中:綜合為內(nèi)梅羅綜合污染指數(shù),平均為所有單項污染指數(shù)的平均值,Pmax為各單項污染指數(shù)中的最大值。

    1.3.3 地累積指數(shù)法

    地累積指數(shù)(geo)又稱 Muller 指數(shù), 反映了單一重金屬元素的污染水平, 在計算過程中加入了巖石地質(zhì)及其他因素的修正指數(shù), 充分考慮到人為活動對重金屬污染的影響, 其公式為[13]:

    表1 采樣點位置

    式中,CB分別為第種重金屬的實測值和土壤中該重金屬的地化學背景值。

    1.3.4 潛在生態(tài)風險指數(shù)法

    潛在生態(tài)風險指數(shù)法由瑞典科學家Hakanson在1980年創(chuàng)立, 適用于對沉積物或土壤重金屬元素危害的程度以及生態(tài)風險進行評價。該方法是應用沉積學原理, 并結(jié)合其他生物、化學領(lǐng)域的理論知識, 對重金屬危害程度進行定量的劃分, 計算公式如下[14]:

    式中,為土壤中重金屬的潛在生態(tài)風險指數(shù),E為重金屬的潛在生態(tài)風險系數(shù),T為重金屬的毒性系數(shù),C為重金屬實測值,C為重金屬的背景值。

    1.4 數(shù)據(jù)分析

    采用SPSS 25.0軟件和Excle2016對所測數(shù)據(jù)進行統(tǒng)計分析; 采用單因素方差分析(One-Way ANOVA)和最小顯著差異法(LSD)比較不同數(shù)據(jù)組間的差異, 采用Spearman相關(guān)系數(shù)對各指標間相關(guān)性進行分析, 顯著性采用雙尾檢驗; 采用Origin2018作圖。

    2 結(jié)果與分析

    2.1 瑪曲高寒退化草甸植物群落特征

    對植物群落結(jié)構(gòu)調(diào)查結(jié)果顯示研究區(qū)草地有22個物種, 隸屬13科21屬。禾本科(Gramineae)和菊科(Compositae)最多, 各有4種植物, 各占總種數(shù)的17%。垂穗披堿草優(yōu)勢度為10.94, 生物量為80.56g·m–2, 矮嵩草(優(yōu)勢度為1.71, 生物量為11.45g·m–2, 雜草類生物量占112.40g·m–2。

    高寒草甸發(fā)生的主要演替過程可以分為禾草-矮嵩草群落、矮嵩草群落、小嵩草群落、雜類草“黑土型”次生裸地4個典型過程[15]。本研究區(qū)中, 群落優(yōu)勢種為露蕊烏頭(優(yōu)勢度為16.67%, 科優(yōu)勢度為雜類草>禾本科>菊科>莎草科>豆科, 生物量大小為雜類草>禾本科>莎草科, 雜類草的優(yōu)勢度、生物量大于禾草, 且雜草種類多, 研究區(qū)草地呈現(xiàn)一定程度的退化。

    2.2 土壤重金屬含量描述性統(tǒng)計

    瑪曲高寒草甸重金屬質(zhì)量分數(shù)的描述性統(tǒng)計見表4。背景值采用青海省土壤環(huán)境背景值[16], 與土壤背景值相比, 瑪曲As、Cd的質(zhì)量分數(shù)高于背景值, 均值與背景值之比分別是5.38和1.82, 表明瑪曲As和Cd的含量高, Cr、Cu、Zn、Hg、Pb質(zhì)量分數(shù)相近, 而Ni的質(zhì)量分數(shù)偏低。土壤重金屬的變異系數(shù)由大到小的順序為Cr> Hg> Cd> Cu> As> Ni> Zn> Pb, 其中Ni、Cu、Zn、As、Cd、Pb屬于弱變異(CV≤0.1), 說明其含量在空間上變化小, 受人為活動影響弱, Cr和Hg屬中等變異(0.1

    表2 污染指數(shù)評價等級

    表3 重金屬毒性系數(shù)

    2.3 重金屬污染狀況評價

    以瑪曲土壤重金屬背景值分析, 瑪曲高寒草甸單因子污染指數(shù)由大到小為Cd(1.78)>As(1.50)>Pb (1.02)>Hg(1.00)>Ni(0.97)>Cu(0.94)>Zn(0.91)>Cr(0.40)(圖1), As和Cd含量如前文所述, 含量較高, 使用單因子評價方法As和Cd為輕度污染, 其他重金屬元素均為無污染清潔水平。以國標分析, 單因子指數(shù)污染大小為Cd(1.19)>Ni(0.77)>As(0.75)> Zn(0.36)>Cu(0.33)>Hg(0.30)>Pb(0.12)>Cr(–0.02), Cd輕度污染外其余重金屬元素均處于正常水平。

    內(nèi)梅羅指數(shù)法是一種綜合考慮各種污染元素極值的計權(quán)型多因子環(huán)境質(zhì)量評價指數(shù), 涉及單因子污染指數(shù)的平均值和最大值, 可以反映各重金屬對土壤的綜合作用, 以瑪曲土壤重金屬背景值和國標分析, 瑪曲重金屬內(nèi)梅羅綜合污染指數(shù)分別是為0.37和0.25(P<0.7), 表明瑪曲未受到重金屬污染, 屬清潔水平(表5)。

    圖2反映了重金屬地累積指數(shù)范圍和平均值, 與瑪曲土壤背景值相比, 重金屬污染強度為AS>Cd>Cr>Pb>Zn> Hg>Cu>Ni, 土壤中As的累積程度高, 平均值為1.84, 屬于中度污染, 其余重金屬都屬于未污染, 其中Cd的平均值為0.27, 其余重金屬geo指數(shù)都為負值, 表現(xiàn)出健康無污染水平, 與單因子指數(shù)評價結(jié)果一致。與國標相比, As的平均累積指數(shù)為0.84, 屬于未受污染與輕度污染之間, 瑪曲高寒草地重金屬累積程度程度均小于0, 表現(xiàn)為未受污染。

    2.4 潛在生態(tài)風險分析

    各重金屬的潛在生態(tài)風險系數(shù)計算結(jié)果見圖3(a)。以土壤背景值為標準, 重金屬潛在生態(tài)風險程度大小順序為Cd>Hg>As>Pb>Ni>Cu>Zn>Cr, 土壤Cd的E平均值為53.53, 屬中等風險, 其余重金屬元素潛在生態(tài)風險系數(shù)均小于40, 生態(tài)風險等級較低, 屬輕微風險。和國標相比, 各元素風險大小為Cd>Hg>As>Ni>Cu> Pb> Zn> Cr, 其Eir值均小于40, 表現(xiàn)為無風險或輕微風險。

    圖3(b)是各采樣點綜合風險指數(shù), 以土壤背景值為標準, 除14、16號點位綜合生態(tài)風險指數(shù)在145以上外, 其余點位綜合生態(tài)風險指數(shù)均為低風險水平, 在12、15、4、3、19點位綜合生態(tài)風險指數(shù)低,

    表4 土壤重金屬含量描述性統(tǒng)計

    注: 背景值采用青海省土壤環(huán)境背景值[16]。

    均在111以下。和國標相比, 各樣點綜合生態(tài)風險指數(shù)均小于80, 其生態(tài)風險指數(shù)低, 表現(xiàn)為無風險或輕微風險。

    2.5 土壤重金屬同源性分析

    2.5.1 土壤重金屬相關(guān)性分析

    土壤重金屬元素相關(guān)性分析見表6。其中Cr與Zn、Ni、Hg、Pb、Cd、AS無顯著相關(guān), Hg與Cu、Ni、Pb、Cd、As、Zn無顯著相關(guān), 但是Cr與Zn, Cd與As顯著相關(guān)(< 0.05), 其余重金屬兩兩之間呈現(xiàn)極顯著相關(guān)(< 0.01)。顯著相關(guān)的重金屬, 可能具有共同或相似的來源, 重金屬之間存在不顯著相關(guān), 則可能存在多種重金屬來源。

    圖1 土壤重金屬背景值(a)與國標(b)單因子污染評價

    Figure 1 Soil heavy metal background value(a) and national standard(b) single factor pollution assessment

    表5 土壤重金屬單因子和內(nèi)梅羅污染指數(shù)

    圖2 土壤重金屬地累計指數(shù)

    Figure 2 Accumulation index of soil heavy metals

    2.5.2 主成分分析(PCA)

    主成分分析和聚類分析是判別重金屬來源的有效方法。利用SPSS25對土壤重金屬含量進行KMO檢測, 得到的統(tǒng)計量值為0.661(0.661>0.5), 說明數(shù)據(jù)適用于主成分分析。對Kaiser標準化后的因子進行Varimax正交旋轉(zhuǎn), 得到了2個特征值大于1的主成分, 累計方差貢獻率為79.508%, 可解釋土壤重金屬元素的大部分信息(表7)。

    第一主成分(F1)的貢獻率是63.296%, 載荷較高的元素是Ni、Cu、Zn、As、Cd、Pb, 前文分析已知, Ni、Cu、Zn、As、Pb的變異系數(shù)小, 屬弱變異, 受人為活動因素的影響小, 且Ni、Cu、Zn、As、Cd、Pb之間存在極顯著正相關(guān)(<0.01), 推斷這6種重金屬具有相同或相似的來源, 有較強的同源性。第二主成分(F2)的貢獻率是16.212%, 載荷較高的元素是Hg、Cr, Hg和Cr屬中等變異, 均值與土壤背景值之比高于1, 受到一定的人為因素干擾。Hg和Cr的相關(guān)系數(shù)為0.243, 無顯著相關(guān), 可能有不同的來源。

    主成分分析僅能反映各個重金屬元素歸屬于某一主成分因素, 而不能說明各重金屬元素對主成分的貢獻率。各重金屬元素對兩個主成分的貢獻率如圖4, 對F1來說, 貢獻率大小順序為Ni> Cu >Cd> Zn> Pb> As> Cr> Hg, 在F2中, 重金屬貢獻率為Hg> Cr> As> Pb> Zn> Cd> Cu> Ni, 其中Ni對F1的貢獻率最高, 貢獻率82%, Hg對F2貢獻率最大, 達到98%, As對F1和F2的貢獻率均高, 分別是36%和64%。

    圖3 土壤重金屬潛在生態(tài)風險(a)與各樣點綜合生態(tài)風險(b)評價

    Figure 3 Assessment of potential ecological risk of soil heavy metals(a) and comprehensive ecological risk of various points(b)

    表6 土壤重金屬相關(guān)性矩陣

    注: *表示在0.05水平(雙尾)相關(guān)性顯著; **表示在0.01水平(雙尾)相關(guān)性顯著。

    表7 土壤重金屬主成分分析

    以8種重金屬元素均值標準化后的質(zhì)量分數(shù)作為分類指標, 對樣本點進行聚類分析, 結(jié)果如圖5所示, 當組間距離為1時, Ni和Pb聚為一類, 距離為2時Zn加入, 距離為4時Cu加入, 聚類為5時As加入, 聚類為6時Cd加入, Ni、Pb、Zn、Cu、As、Cd聚為一類; 當組間距離為16時Cr和Hg聚為一類。8種元素共分為兩大類, 結(jié)果與主成分分析中2個主成分的元素構(gòu)成一致, 聚類分析在主成分分析結(jié)果的基礎上, 進一步提供了同源污染物的優(yōu)先順序。

    3 討論

    3.1 瑪曲高寒退化草甸重金屬含量分布

    瑪曲退化高寒草甸重金屬含量整體處于較低水平, 重金屬含量較低, 本研究區(qū)不同采樣點重金屬含量存在顯著差異?,斍嘶莸橥寥繟s、Cd質(zhì)量分數(shù)高于背景值, Ni含量偏低, 土壤Cr和Hg屬中等變異(0.1< CV< 1), 可能存在點源污染。

    As含量較高可能是由于青藏高原廣泛分布著富含As的頁巖, 導致As的質(zhì)量分數(shù)大于土壤背景值[19], 而且使用含As肥料, 過量使用殺蟲劑和除草劑等, 都會引As含量增加[20]。江濤等[21]研究發(fā)現(xiàn), 重金屬元素遷移能力為: Cd> Pb> Zn> Cu> Ni> Cr, 其中Cd的遷移能力強, Ni的遷移能力弱, 重金屬含量可能受到不同元素遷移能力的影響, 遷移能力強的元素會大量進入土壤造成富集, 遷移能力弱的元素難以進入土壤導致其含量低于土壤背景值。青藏高原廣布由超鎂鐵巖發(fā)育而來的土壤, Cr和Ni含量豐富[4], Ni多以無機和有機絡合物的形式溶解于水流失, 所以Ni含量會降低, Cr、Cd不易被植物被吸收或在土壤中沉淀[22]。

    圖4 重金屬對主成分的貢獻率

    Figure 4 Contribution rate of heavy metals to principal components

    圖5 聚類分析譜系圖

    Figure 5 Genealogy diagram of cluster analysis

    有研究表明青藏高原表層土壤中Cd含量較高主要受交通因素的影響[18], 通過汽車燃料的不完全燃燒及輪胎的磨損、撕裂等釋放Cd[3,4], Zhang等[17]研究表明青藏公路沿線表層土壤中Cd、Pb和Zn含量較高。Cr由含超鎂鐵巖的土壤發(fā)育而來[4], 含量豐富, 研究區(qū)西側(cè)為瑪曲縣X417西扎段公路, 重金屬含量可能受到交通運輸因素的影響。本研究中潛在生態(tài)風險分析表明Cd屬于中等風險, 結(jié)果與國內(nèi)其他研究的結(jié)果一致[23]。

    3.2 瑪曲高寒退化草地重金屬可能來源分析

    一般而言, Cr和Ni在土壤中的含量主要受成土母質(zhì)的含量影響[24], 這兩種重金屬元素在表生地球化學過程中均為親鐵元素, 在土壤中的聚集與含量極有可能由成土過程主導, 并且Cu含量受成土母質(zhì)的影響大, Cu在F1的載荷高達0.9以上, 在本研究中Cr在F1和F2上均有較高載荷。土壤中Zn的累積主要與有色金屬加工等工業(yè)及牲畜糞便有關(guān)[25], 研究區(qū)以高寒草甸牦牛養(yǎng)殖為主, 對草地使用農(nóng)藥、化肥、牦牛產(chǎn)生的糞便等也可能成為Zn、As元素污染的重要來源。

    相比其他重金屬, Hg的人為污染來源相對集中, 主要來源于燃煤燃燒、城市生活垃圾焚燒、建筑材料生產(chǎn)等[18]。Hg主要是通過燃煤、電廠和焦化廠等向大氣中排放, 大氣中的Hg很穩(wěn)定, 通過大氣干濕沉降的方式進入土壤并累積[26]。土壤重金屬富集的重要途徑是大氣沉降[27], 隨著顆粒物的遠距離飄散, Hg, Cr可分布整個研究區(qū), 造成區(qū)域污染。研究發(fā)現(xiàn)[28], 畜禽糞便是農(nóng)田土壤Hg的主要來源。由此可見, Hg, Cr的累積可能來自于大氣沉降。

    化肥農(nóng)藥不合理施用也是造成土壤污染的途徑之一[29], 化肥農(nóng)藥中通常含有微量Cd、Zn、Cu、As和Hg等重金屬元素, 即使是符合國家標準的化肥農(nóng)藥, 在長期施用的過程中農(nóng)田土壤也會不斷富集重金屬。有機肥料中重金屬的殘留的增加, 在土壤中會累積形成不同的污染情況, 導致植被污染或不能正常生長[30]?,斍吆莸橹饕糜谀撩穹拍羀31], 草地施加了大量的氮磷等營養(yǎng)物質(zhì), 導致人為因素影響重金屬含量。研究顯示[32]中國土壤 Cd、Cu、Pb和Zn的年大氣輸入量分別為0.4、10.67、20和64.08 mg·m–2, 汽車尾氣排放、垃圾焚燒以及建筑材料粉塵等因素產(chǎn)生的重金屬微粒通過特定途徑沉降到土壤環(huán)境中, 進一步加劇土壤重金屬富集和污染。

    由于自然條件的難控制性和人為活動干擾的多樣性, 重金屬元素在環(huán)境中的遷移、富集和轉(zhuǎn)化較為復雜, 某一個地區(qū)的重金屬元素可能受到多種因素的共同作用, 包括本研究未涉及的其他重金屬元素, 同時, 重金屬元素之間也存在相互影響, 今后對土壤重金屬含量應當予以長期監(jiān)測, 進一步探究各元素之間的相互關(guān)系, 為某一地區(qū)重金屬污染來源進行系統(tǒng)評價和分析。

    4 結(jié)論

    1.地累積指數(shù)反映As屬于未受污染與輕微污染之間, 各重金屬元素均無潛在生態(tài)風險; 以國標和土壤背景值為標準, 內(nèi)梅羅綜合污染指數(shù)表明瑪曲未受到重金屬污染。

    2. 土壤中Ni、Cu、Zn、As、Cd、Pb等元素具有相同或相似的來源, Ni和Pb的來源優(yōu)先順序最大, Cr和Hg含量呈現(xiàn)弱相關(guān)。

    [1] 陳世寶, 王萌, 李杉杉, 等. 中國農(nóng)田土壤重金屬污染防治現(xiàn)狀與問題思考[J]. 地學前緣, 2019, 26(6): 35–41.

    [2] LUO X S, JING D, BO X, et al. Incorporating bioaccessibility into human health risk assessments of heavy metals in urban park soils[J]. Science of the Total Environment,2012, 424: 88–96.

    [3] SHENG J J, WANG X P, GONG P, et al. Heavy metals of the Tibetan top soils[J]. Environmental Science&Pollution Research, 2012, 19(8): 3362–3370

    [4] WU J, DUAN D, JIAN L, et al. Inorganic pollution around the Qinghai-Tibet Plateau: An overview of the current observations[J]. Science of The Total Environment, 2016, 550: 628–636.

    [5] 高興川, 曹小曙, 李濤, 等. 1976-2016年青藏高原地區(qū)通達性空間格局演變[J]. 地理學報, 2019, 74(6): 1190–1204.

    [6] 王冠星, 閆學東, 張凡, 等. 青藏高原路側(cè)土壤重金屬含量分布規(guī)律及影響因素研究[J]. 環(huán)境科學學報, 2014, 34(2): 431–438.

    [7] WANG X P, YANG H D, GONG P, et al. One century sedimentary records of polycyclic aromatic hydrocarbons, mercury and trace elements in the Qinghai Lake, Tibetan Plateau. [J]. Environmental Pollution, 2010, 158(10): 3065–3070.

    [8] 楊啟良, 武振中, 陳金陵, 等. 植物修復重金屬污染土壤的研究現(xiàn)狀及其水肥調(diào)控技術(shù)展望[J]. 生態(tài)環(huán)境學報, 2015, 24(6): 1075–1084.

    [9] 楊曉輝. 長江下游典型區(qū)土壤—作物系統(tǒng)重金屬富集特征[D]. 南京: 南京大學, 2013.

    [10] 徐聰, 趙婷, 池海濤, 等. 微波消解ICP-MS法測定土壤及耕作物小麥中的8種重金屬元素[J]. 中國測試, 2019, 45(5): 85–92.

    [11] 周長松, 鄒勝章, 李錄娟, 等. 幾種土壤重金屬污染評價方法的對比[J]. 地球與環(huán)境, 2015, 43(6): 709–713.

    [12] 張金婷, 孫華. 內(nèi)梅羅指數(shù)法和模糊綜合評價法在土壤重金屬污染評價應用中的差異分析[J]. 環(huán)境監(jiān)測管理與技術(shù), 2016, 28(4): 27–31.

    [13] YAQIN J I, FENG Y, JIANHUI W U, et al. Using geoaccumulation index to study source profiles of soil dust in China[J]. Journal of Environmental Sciences, 2008, 20: 571–578.

    [14] JIE F, ZHAO C, LUO Y, et al. Heavy metals in surface sediments of the Jialu River, China: Their relations to environmental factors[J]. Journal of Hazardous Materials, 2014, 270: 102–109.

    [15] 王建兵, 張德罡, 曹廣民, 等. 青藏高原高寒草甸退化演替的分區(qū)特征[J]. 草業(yè)學報, 2013, 22(2): 1–10.

    [16] 郭立本, 郭競世.青海省土壤環(huán)境背景值研究實驗室質(zhì)量控制[J]. 青海環(huán)境, 1995(1): 40–42.

    [17] JIANG F, REN B, HURSTHOUSE A, et al. Distribution, source identification, and ecological-health risks of potentially toxic elements (PTEs) in soil of thallium mine area southwestern Guizhou, China[J]. Environmental Science and Pollution Research, 2019, doi: 10.1007/ s11356-019-04997-3.

    [18] ZHANG H, WANG Z, ZHANG Y, et al. Identification of traffic-related metals and the effects of different environments on their enrichment in roadside soils along the Qinghai–Tibet highway[J]. Science of The Total Environment, 2015(521/522): 160–172.

    [19] LI C, KANG S, CHEN P, et al. Geothermal spring causes arsenic contamination in river waters of the southern Tibetan Plateau, China[J]. Environmental Earth Sciences, 2014, 71(9): 4143–4148.

    [20] 楊輝, 陳國光, 劉紅櫻, 等. 長江下游主要湖泊沉積物重金屬污染及潛在生態(tài)風險評價[J]. 地球與環(huán)境, 2013, 41(2): 160–165.

    [21] 江濤, 林偉穩(wěn), 曹英杰, 等. 梅江流域清涼山水庫沉積物重金屬污染,生態(tài)風險評價及來源解析[J]. 環(huán)境科學, 2020, 41(12): 5410–5418.

    [22] 周會程, 姚玉嬌, 梁婷, 等. 天祝不同退化梯度高寒草甸土壤重金屬污染及風險評價[J]. 生態(tài)環(huán)境學報, 2020, 29(10): 2102–2109.

    [23] 方曉波, 史堅, 廖欣峰, 等.臨安市雷竹林土壤重金屬污染特征及生態(tài)風險評價[J]. 應用生態(tài)學報, 2015, 26(6): 1883–1891.

    [24] 李春芳,王菲,曹文濤, 等. 龍口市污水灌溉區(qū)農(nóng)田重金屬來源、空間分布及污染評價[J]. 環(huán)境科學, 2017, 38(3): 1018–1027.

    [25] 李金峰, 聶兆君, 趙鵬, 等. 畜禽糞便配施對冬小麥產(chǎn)量及Cu、Zn、As在植株累積和土壤中垂直分布的影響[J]. 江蘇農(nóng)業(yè)科學, 2016, 44(4): 137–140.

    [26] 孔絲紡,劉陽生,曾輝, 等. 焚燒源揮發(fā)性重金屬的排放對周邊土壤和植被污染的研究進展[J]. 生態(tài)環(huán)境學報, 2010, 19(4): 985–990.

    [27] 楊潔, 周曉芳, 王新富, 等. 基于PCA/APCS的草原區(qū)煤礦土壤重金屬環(huán)境質(zhì)量評價與定量溯源研究[J]. 草地學報, 2021, 29(3): 537–544.

    [28] SHI T, MA J, WU F, et al. Mass balance-based inventory of heavy metals inputs to and outputs from agricultural soils in Zhejiang Province, China[J]. Science of The Total Environment, 2018, 649: 1269–1280.

    [29] 李福燕, 李許明, 楊帆, 等. 海南島農(nóng)田土壤重金屬污染的評價及其來源分析[J]. 海南大學學報(自然科學版), 2013, 31(3): 211–217.

    [30] 陸泗進, 王業(yè)耀, 何立環(huán). 風險評價代碼法對農(nóng)田土壤重金屬生態(tài)風險的評價[J].環(huán)境化學, 2014, 33(11): 1857–1863.

    [31] 馬青青, 柴林榮, 馬海玲, 等. 瑪曲高寒草甸放牧強度的遙感監(jiān)測[J]. 草業(yè)科學, 2018, 35(5): 941–948.

    [32] LEI L, MA Y, ZHANG S, et al. An inventory of trace element inputs to agricultural soils in China[J]. Journal of Environmental Management, 2009, 90(8): 2524–2530.

    Risk assessment of heavy metal pollution in degraded alpine meadow in Maqu County

    Sun Jiancai1, BAI Hongtao1, XIONG Zenglian3, ZHANG Jiexue2, QUAN Xiaolong2, DENG Deting1, SHI Huilan1,*

    1. College of Eco-Environmental Engineering, Qinghai University, Xining 810016, China 2. The State Key Laboratory of Plateau Ecology and Agriculture, Qinghai University, Xining 810016, China 3.Qinghai Provincial Natural Resources Survey and Monitoring Institute, Qinghai 810016, China

    To analyze the risk assessment of heavy metal pollution in maqu alpine meadow, we collected the surface soil samples (0-20 cm) of degraded alpine meadow in maqu county to determining the contents of eight heavy metals in soil, including Cr, Ni, Cu, Zn, As, Cd, Hg and Pb. Single factor pollution index, Nemerow pollution index and geoaccumulation index were used to evaluate the pollution degree of heavy metals and analyze the potential ecological risk. We used correlation analysis, principal component analysis and cluster analysis to characterize the homologous heavy metals. The results showed that :(1) Therewas greater mass fraction of As and Cd than the background value in soil, while there was lower mass fraction of Ni than the background value in soil. The ranking of variation coefficient was Cr> Hg> Cd> Cu> As> Ni> Zn> Pb. (2) The single factor pollution index of Cdwas slight pollution, the ground accumulation index reflects that As was between unpolluted and slight pollution. Each heavy metal element was no potential ecological risk. Based on the national standard and soil background value, Nemerow composite pollution index showed that maqu was not contaminated by heavy metals. (3) Principal component analysis and correlation analysis showed that Ni, Cu, Zn, As, Cd, Pbin Maqu soil had the same or similar sources.The Ni and Pb had the highest priority order, while Cr and Hg contents showed a weak correlation.

    Alpine meadow; Heavy metal pollution; Pollution index; Potential ecological risk; Principal component analysis

    10.14108/j.cnki.1008-8873.2024.01.021

    X825

    A

    1008-8873(2024)01-177-09

    2021-05-25;

    2021-07-13

    青海省科技廳基礎研究項目(2021-ZJ-715); 國家自然科學基金項目(C170201); 國家重點研發(fā)項目(2017YFC0504800); 2021年促進與加澳新 及拉美地區(qū)科研合作與高層次人才培養(yǎng)項目(留金美[2021]109號)

    孫建財(1997—), 男, 青海樂都人, 碩士研究生, 主要從事草地生態(tài)學研究. E-mail: 2267290531@qq.com

    通信作者:史惠蘭(1977—), 女, 在職博士, 副教授, 主要從事草地生態(tài)學研究. E-mail: hlshi7701@126.com

    孫建財, 白紅桃, 熊增連,等. 瑪曲高寒退化草甸土壤重金屬污染風險評價[J]. 生態(tài)科學, 2024, 43(1): 177–185.

    Sun Jiancai, BAI Hongtao, Xiong Zenglian, et al. Risk assessment of heavy metal pollution in degraded alpine meadow soil in Maqu[J]. Ecological Science, 2024, 43(1): 177–185.

    猜你喜歡
    瑪曲草甸金屬元素
    高山草甸
    The world's narrowest river
    山地草甸
    固體廢物中金屬元素檢測方法研究進展
    云南化工(2021年8期)2021-12-21 06:37:08
    瑪曲:美麗的黃河首曲大草原
    學生天地(2020年28期)2020-06-01 02:18:48
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    微波消解-ICP-MS法同時測定牛蒡子中8種重金屬元素
    中成藥(2018年11期)2018-11-24 02:57:28
    瑪曲行吟
    中國三峽(2016年10期)2017-01-15 13:59:44
    在瑪曲
    紅豆(2016年9期)2016-05-14 11:11:49
    環(huán)境監(jiān)測中重金屬元素分析方法探析
    河南科技(2014年8期)2014-02-27 14:08:07
    亚洲av.av天堂| 久久99热6这里只有精品| 国产成人精品一,二区| 日韩,欧美,国产一区二区三区| 国产成人精品久久久久久| 最新中文字幕久久久久| 亚洲欧美日韩卡通动漫| 三级国产精品片| 巨乳人妻的诱惑在线观看| 精品酒店卫生间| 在线观看免费日韩欧美大片| 国产精品欧美亚洲77777| 欧美日韩国产mv在线观看视频| 亚洲精品456在线播放app| 国产成人午夜福利电影在线观看| 日韩av在线免费看完整版不卡| 国产成人免费观看mmmm| 大香蕉97超碰在线| 91成人精品电影| 成人免费观看视频高清| 亚洲美女视频黄频| 日日爽夜夜爽网站| 久久精品国产亚洲av天美| 不卡视频在线观看欧美| 蜜臀久久99精品久久宅男| √禁漫天堂资源中文www| 大香蕉97超碰在线| 亚洲精品久久成人aⅴ小说| www.av在线官网国产| 国产成人精品无人区| 亚洲精品aⅴ在线观看| 欧美97在线视频| 在线观看免费视频网站a站| 插逼视频在线观看| 18在线观看网站| 一本—道久久a久久精品蜜桃钙片| 色视频在线一区二区三区| 一级毛片我不卡| 亚洲熟女精品中文字幕| 国产一区有黄有色的免费视频| 亚洲精品一二三| av有码第一页| 90打野战视频偷拍视频| 免费女性裸体啪啪无遮挡网站| 亚洲激情五月婷婷啪啪| 在线观看免费日韩欧美大片| 全区人妻精品视频| kizo精华| 免费黄频网站在线观看国产| 久久久久久人妻| av线在线观看网站| 国产免费又黄又爽又色| 午夜免费男女啪啪视频观看| 欧美激情国产日韩精品一区| 性高湖久久久久久久久免费观看| 在线看a的网站| 亚洲熟女精品中文字幕| 精品卡一卡二卡四卡免费| 丝瓜视频免费看黄片| 国产成人精品福利久久| 日韩大片免费观看网站| 最近2019中文字幕mv第一页| 久久av网站| 国产黄频视频在线观看| 久久 成人 亚洲| 自线自在国产av| 欧美最新免费一区二区三区| 大码成人一级视频| 少妇 在线观看| 亚洲人成网站在线观看播放| 日本av免费视频播放| 黑人欧美特级aaaaaa片| 熟妇人妻不卡中文字幕| 制服丝袜香蕉在线| 丰满乱子伦码专区| 如何舔出高潮| 亚洲成av片中文字幕在线观看 | 亚洲人成网站在线观看播放| 日韩中字成人| 国产成人a∨麻豆精品| 国产爽快片一区二区三区| 性高湖久久久久久久久免费观看| 五月玫瑰六月丁香| 黄色视频在线播放观看不卡| 精品少妇黑人巨大在线播放| 亚洲熟女精品中文字幕| 好男人视频免费观看在线| 国产乱人偷精品视频| 九九在线视频观看精品| 国国产精品蜜臀av免费| 国产一级毛片在线| 国产av码专区亚洲av| 色吧在线观看| 亚洲人与动物交配视频| 成年人免费黄色播放视频| 男人操女人黄网站| 热re99久久国产66热| 成年av动漫网址| 老司机影院毛片| 色婷婷av一区二区三区视频| 久久久久精品性色| 在线精品无人区一区二区三| 日韩大片免费观看网站| 人妻系列 视频| 伦精品一区二区三区| 午夜av观看不卡| 啦啦啦中文免费视频观看日本| 熟妇人妻不卡中文字幕| 寂寞人妻少妇视频99o| 色哟哟·www| 秋霞在线观看毛片| 久久精品国产亚洲av涩爱| 国产一区二区在线观看日韩| 国产精品嫩草影院av在线观看| h视频一区二区三区| 亚洲精品视频女| 日韩成人av中文字幕在线观看| 午夜精品国产一区二区电影| 亚洲欧美一区二区三区黑人 | 女性被躁到高潮视频| 午夜视频国产福利| 街头女战士在线观看网站| 精品少妇久久久久久888优播| 内地一区二区视频在线| 国产永久视频网站| 啦啦啦视频在线资源免费观看| 九九爱精品视频在线观看| 一级毛片黄色毛片免费观看视频| 亚洲av电影在线观看一区二区三区| 国产无遮挡羞羞视频在线观看| 伦理电影大哥的女人| 亚洲精品自拍成人| 日韩伦理黄色片| 精品久久久精品久久久| 亚洲欧美成人精品一区二区| 草草在线视频免费看| 亚洲综合精品二区| 色5月婷婷丁香| 久久久久精品久久久久真实原创| 欧美 亚洲 国产 日韩一| 在线观看免费视频网站a站| 一区二区三区精品91| 久久午夜综合久久蜜桃| 十八禁高潮呻吟视频| 国产成人精品无人区| 亚洲欧洲精品一区二区精品久久久 | 成人毛片a级毛片在线播放| 欧美精品av麻豆av| 伊人久久国产一区二区| 婷婷色麻豆天堂久久| 最近中文字幕高清免费大全6| 如何舔出高潮| 国产在线视频一区二区| 搡老乐熟女国产| 99久久精品国产国产毛片| 欧美亚洲日本最大视频资源| 午夜免费男女啪啪视频观看| av在线app专区| 女性被躁到高潮视频| 久久精品国产鲁丝片午夜精品| 国产毛片在线视频| 18在线观看网站| 欧美国产精品一级二级三级| 女性被躁到高潮视频| 极品少妇高潮喷水抽搐| 国产一区二区三区综合在线观看 | 国产日韩一区二区三区精品不卡| 久久久久国产精品人妻一区二区| 亚洲三级黄色毛片| 99国产综合亚洲精品| 欧美bdsm另类| 在线观看美女被高潮喷水网站| 飞空精品影院首页| 青青草视频在线视频观看| 国产男人的电影天堂91| 中文字幕av电影在线播放| 999精品在线视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 天堂俺去俺来也www色官网| 亚洲精华国产精华液的使用体验| 国产永久视频网站| 少妇的丰满在线观看| 又黄又粗又硬又大视频| 在线观看一区二区三区激情| 99国产综合亚洲精品| 国产精品秋霞免费鲁丝片| 咕卡用的链子| 亚洲欧美色中文字幕在线| 亚洲激情五月婷婷啪啪| 人人澡人人妻人| 久久久精品区二区三区| 久久这里有精品视频免费| 91午夜精品亚洲一区二区三区| 成人国语在线视频| 国产av国产精品国产| 国产在线一区二区三区精| 午夜福利在线观看免费完整高清在| 看非洲黑人一级黄片| 亚洲成人手机| 激情视频va一区二区三区| 午夜影院在线不卡| 狠狠婷婷综合久久久久久88av| 欧美精品人与动牲交sv欧美| av片东京热男人的天堂| av在线播放精品| 美女视频免费永久观看网站| 国产高清三级在线| 春色校园在线视频观看| 中文字幕免费在线视频6| 精品酒店卫生间| 国产欧美日韩综合在线一区二区| 高清不卡的av网站| 亚洲精品第二区| 国产综合精华液| 久久精品国产亚洲av天美| 亚洲在久久综合| 国产精品 国内视频| 亚洲综合色网址| 五月天丁香电影| 欧美变态另类bdsm刘玥| 五月开心婷婷网| 又大又黄又爽视频免费| 一本久久精品| 亚洲丝袜综合中文字幕| 性高湖久久久久久久久免费观看| 午夜福利在线观看免费完整高清在| 国产xxxxx性猛交| 国产成人精品在线电影| 国产一区二区在线观看av| 午夜福利视频精品| 99久久人妻综合| 成人亚洲精品一区在线观看| 日韩大片免费观看网站| 免费黄色在线免费观看| 免费看不卡的av| 日韩精品有码人妻一区| 亚洲综合精品二区| 久久国产精品男人的天堂亚洲 | 亚洲一区二区三区欧美精品| 极品人妻少妇av视频| 一级毛片电影观看| 亚洲国产色片| 少妇精品久久久久久久| 国产在线视频一区二区| 久久99一区二区三区| 91aial.com中文字幕在线观看| 免费不卡的大黄色大毛片视频在线观看| 美国免费a级毛片| 成人二区视频| 日韩成人伦理影院| 最新的欧美精品一区二区| a级毛片在线看网站| 91久久精品国产一区二区三区| 成年人午夜在线观看视频| 欧美精品高潮呻吟av久久| 你懂的网址亚洲精品在线观看| 两性夫妻黄色片 | 老司机亚洲免费影院| 欧美精品高潮呻吟av久久| 欧美xxⅹ黑人| 久久久国产一区二区| 男人爽女人下面视频在线观看| 国产成人午夜福利电影在线观看| 国产一区二区三区av在线| 精品少妇久久久久久888优播| 国产深夜福利视频在线观看| 国产成人精品无人区| 一级片'在线观看视频| 精品国产国语对白av| 久久精品国产亚洲av天美| 久久精品人人爽人人爽视色| 日本wwww免费看| 日韩欧美精品免费久久| 一级爰片在线观看| 国产福利在线免费观看视频| 久久久精品94久久精品| 精品99又大又爽又粗少妇毛片| 亚洲av日韩在线播放| 久久热在线av| 韩国精品一区二区三区 | 日韩 亚洲 欧美在线| 黄色配什么色好看| 国产一区亚洲一区在线观看| 欧美 日韩 精品 国产| 插逼视频在线观看| 夫妻午夜视频| 高清在线视频一区二区三区| 国产女主播在线喷水免费视频网站| 人人澡人人妻人| 伊人久久国产一区二区| 精品人妻偷拍中文字幕| 免费大片18禁| 老女人水多毛片| 亚洲精品乱码久久久久久按摩| 日韩成人av中文字幕在线观看| 丝袜喷水一区| 91午夜精品亚洲一区二区三区| 国产成人免费无遮挡视频| 新久久久久国产一级毛片| 母亲3免费完整高清在线观看 | 黑人高潮一二区| 免费看av在线观看网站| 精品酒店卫生间| xxx大片免费视频| 如何舔出高潮| 黑丝袜美女国产一区| 欧美bdsm另类| 午夜福利影视在线免费观看| 亚洲欧美色中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 99热网站在线观看| 欧美日韩av久久| 99国产精品免费福利视频| 亚洲美女视频黄频| 成年人午夜在线观看视频| 人人妻人人爽人人添夜夜欢视频| 一区二区av电影网| 欧美人与性动交α欧美软件 | 亚洲成国产人片在线观看| 51国产日韩欧美| 亚洲国产日韩一区二区| 日韩大片免费观看网站| 久久99热6这里只有精品| 99精国产麻豆久久婷婷| 久久久久人妻精品一区果冻| av在线老鸭窝| 91aial.com中文字幕在线观看| 夜夜骑夜夜射夜夜干| av有码第一页| 精品少妇内射三级| 国产男女超爽视频在线观看| 亚洲高清免费不卡视频| 国产av精品麻豆| 国产1区2区3区精品| 日本黄色日本黄色录像| 啦啦啦视频在线资源免费观看| 亚洲在久久综合| 久久人人爽av亚洲精品天堂| 亚洲,欧美,日韩| 狂野欧美激情性xxxx在线观看| 久热这里只有精品99| 十分钟在线观看高清视频www| 亚洲av在线观看美女高潮| 亚洲欧美成人综合另类久久久| 亚洲精品456在线播放app| 国产精品成人在线| 超碰97精品在线观看| 国产精品成人在线| 性色av一级| 国产精品偷伦视频观看了| 日韩一区二区三区影片| 街头女战士在线观看网站| 青春草国产在线视频| 亚洲精品久久久久久婷婷小说| 少妇人妻精品综合一区二区| 一区二区三区精品91| 麻豆乱淫一区二区| 亚洲精品,欧美精品| 99热这里只有是精品在线观看| 国产精品国产三级专区第一集| 亚洲五月色婷婷综合| 美女大奶头黄色视频| 老司机亚洲免费影院| 欧美xxⅹ黑人| 国产精品一区www在线观看| 男的添女的下面高潮视频| 9热在线视频观看99| 边亲边吃奶的免费视频| 97精品久久久久久久久久精品| 成人毛片a级毛片在线播放| 9热在线视频观看99| av网站免费在线观看视频| 久久久久网色| 日韩一本色道免费dvd| 一二三四在线观看免费中文在 | 亚洲欧美清纯卡通| 国产免费一级a男人的天堂| 免费日韩欧美在线观看| 交换朋友夫妻互换小说| 欧美xxxx性猛交bbbb| 国产精品 国内视频| 国产黄色视频一区二区在线观看| 国产黄色免费在线视频| 婷婷色av中文字幕| 国产国语露脸激情在线看| 精品99又大又爽又粗少妇毛片| 日韩精品有码人妻一区| 乱码一卡2卡4卡精品| 少妇的丰满在线观看| 日本欧美视频一区| 欧美国产精品一级二级三级| 中国美白少妇内射xxxbb| 精品第一国产精品| 国产成人精品福利久久| 国语对白做爰xxxⅹ性视频网站| 国产有黄有色有爽视频| 成年人午夜在线观看视频| 国产永久视频网站| 成年av动漫网址| 国产欧美日韩综合在线一区二区| 极品少妇高潮喷水抽搐| 久久人人爽人人爽人人片va| 80岁老熟妇乱子伦牲交| 欧美性感艳星| 五月开心婷婷网| 精品国产国语对白av| 在线天堂最新版资源| 国产精品嫩草影院av在线观看| 亚洲av欧美aⅴ国产| 亚洲欧美一区二区三区国产| 91aial.com中文字幕在线观看| 丝袜在线中文字幕| 亚洲欧美成人综合另类久久久| 中文字幕亚洲精品专区| 女人精品久久久久毛片| 美女中出高潮动态图| 欧美国产精品一级二级三级| 深夜精品福利| 久久国产亚洲av麻豆专区| 久久精品人人爽人人爽视色| 欧美 亚洲 国产 日韩一| 中文天堂在线官网| 国产69精品久久久久777片| 婷婷色综合大香蕉| 乱码一卡2卡4卡精品| 在线观看三级黄色| 巨乳人妻的诱惑在线观看| 亚洲欧洲日产国产| 亚洲成av片中文字幕在线观看 | 久久国内精品自在自线图片| 欧美精品人与动牲交sv欧美| 啦啦啦啦在线视频资源| 搡女人真爽免费视频火全软件| 老女人水多毛片| 国产国拍精品亚洲av在线观看| 国产av一区二区精品久久| 久久ye,这里只有精品| 十八禁网站网址无遮挡| 美女内射精品一级片tv| 纵有疾风起免费观看全集完整版| 一边摸一边做爽爽视频免费| 观看av在线不卡| 精品人妻一区二区三区麻豆| 精品久久国产蜜桃| 精品人妻熟女毛片av久久网站| 热re99久久国产66热| 久久人妻熟女aⅴ| 曰老女人黄片| 蜜臀久久99精品久久宅男| 18禁动态无遮挡网站| 国产成人91sexporn| 男人添女人高潮全过程视频| 亚洲情色 制服丝袜| 久久久久国产网址| 中文字幕精品免费在线观看视频 | 国产免费现黄频在线看| 美女福利国产在线| 国产在线视频一区二区| 欧美日韩视频高清一区二区三区二| 国产av一区二区精品久久| 自拍欧美九色日韩亚洲蝌蚪91| 久久人人97超碰香蕉20202| 日本与韩国留学比较| 欧美精品亚洲一区二区| 在线看a的网站| 日韩免费高清中文字幕av| 亚洲成色77777| 热re99久久国产66热| 免费高清在线观看日韩| 校园人妻丝袜中文字幕| 国产精品不卡视频一区二区| 狂野欧美激情性xxxx在线观看| 少妇被粗大的猛进出69影院 | 一二三四中文在线观看免费高清| 久久精品国产亚洲av涩爱| 两个人看的免费小视频| 少妇被粗大的猛进出69影院 | 男人爽女人下面视频在线观看| 大话2 男鬼变身卡| 国内精品宾馆在线| 制服丝袜香蕉在线| 在现免费观看毛片| 国产福利在线免费观看视频| 色婷婷久久久亚洲欧美| 亚洲av日韩在线播放| 在线 av 中文字幕| 日本免费在线观看一区| 亚洲精品国产av成人精品| 大陆偷拍与自拍| 国产av国产精品国产| 国产一区二区三区综合在线观看 | 国产成人午夜福利电影在线观看| 国产一区亚洲一区在线观看| 免费黄频网站在线观看国产| 一本色道久久久久久精品综合| 男女午夜视频在线观看 | 伊人久久国产一区二区| 精品一区二区三区四区五区乱码 | 日韩电影二区| 高清av免费在线| 久久久久久久亚洲中文字幕| 日韩中文字幕视频在线看片| 天天影视国产精品| 国产欧美另类精品又又久久亚洲欧美| 国产成人精品一,二区| 国产男女内射视频| 80岁老熟妇乱子伦牲交| 国产精品一区二区在线不卡| 欧美老熟妇乱子伦牲交| 成人18禁高潮啪啪吃奶动态图| 久热这里只有精品99| 有码 亚洲区| 人妻 亚洲 视频| 免费av中文字幕在线| 久久精品aⅴ一区二区三区四区 | 如日韩欧美国产精品一区二区三区| 在线免费观看不下载黄p国产| 亚洲av国产av综合av卡| 一级黄片播放器| 免费人妻精品一区二区三区视频| 在线免费观看不下载黄p国产| 丝袜人妻中文字幕| 99热这里只有是精品在线观看| 大香蕉久久成人网| 国产无遮挡羞羞视频在线观看| 亚洲熟女精品中文字幕| 亚洲国产欧美日韩在线播放| 中文字幕人妻熟女乱码| 国产欧美日韩一区二区三区在线| 免费在线观看黄色视频的| 久久精品国产亚洲av天美| 国产一区二区在线观看av| 亚洲少妇的诱惑av| 少妇的丰满在线观看| 欧美97在线视频| 国产精品三级大全| 亚洲人成网站在线观看播放| 国产综合精华液| 全区人妻精品视频| 热re99久久精品国产66热6| 爱豆传媒免费全集在线观看| 亚洲在久久综合| 国产精品嫩草影院av在线观看| 日本色播在线视频| 成人综合一区亚洲| 国产1区2区3区精品| 一级片免费观看大全| 丝瓜视频免费看黄片| 精品一区在线观看国产| 99国产综合亚洲精品| 国产亚洲一区二区精品| 各种免费的搞黄视频| 亚洲成av片中文字幕在线观看 | 三级国产精品片| 美女国产视频在线观看| 人妻一区二区av| 热99国产精品久久久久久7| 69精品国产乱码久久久| 亚洲三级黄色毛片| 久久国产亚洲av麻豆专区| 一级,二级,三级黄色视频| 伊人久久国产一区二区| 9色porny在线观看| 在线看a的网站| 中文字幕免费在线视频6| 欧美少妇被猛烈插入视频| 亚洲人成77777在线视频| 国产乱来视频区| 99久久中文字幕三级久久日本| a级片在线免费高清观看视频| 亚洲国产最新在线播放| 黄片无遮挡物在线观看| 在线观看一区二区三区激情| 免费观看无遮挡的男女| 亚洲av国产av综合av卡| 少妇精品久久久久久久| av在线播放精品| a级毛片在线看网站| 新久久久久国产一级毛片| 国产成人av激情在线播放| 成人毛片60女人毛片免费| 日韩制服丝袜自拍偷拍| 丰满迷人的少妇在线观看| av.在线天堂| 久久久久久人人人人人| av一本久久久久| 国产精品久久久久久久电影| 如何舔出高潮| 亚洲,一卡二卡三卡| 日韩av不卡免费在线播放| 精品亚洲成a人片在线观看| 日韩成人av中文字幕在线观看| 精品福利永久在线观看| 黑人巨大精品欧美一区二区蜜桃 | 久久久国产精品麻豆| 国产乱来视频区| 午夜激情久久久久久久| 人人妻人人澡人人爽人人夜夜| 建设人人有责人人尽责人人享有的| 亚洲av福利一区| 波野结衣二区三区在线| 久久99热6这里只有精品| 久久人人爽人人爽人人片va| 国产精品一区二区在线观看99| 51国产日韩欧美| 各种免费的搞黄视频| 五月开心婷婷网| 天堂俺去俺来也www色官网| 一级片'在线观看视频| 久久国产精品大桥未久av| 婷婷成人精品国产| 国国产精品蜜臀av免费| 中文字幕人妻丝袜制服| 久久久久久久国产电影| 久久热在线av|