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

    馬鈴薯種植區(qū)不同土壤類型鎘賦存形態(tài)及風險評估

    2020-03-24 23:50:49婁飛符東順付天嶺蔡娜楊世梅何騰兵
    南方農(nóng)業(yè)學報 2020年12期
    關(guān)鍵詞:風險評估馬鈴薯有效性

    婁飛 符東順 付天嶺 蔡娜 楊世梅 何騰兵

    摘要:【目的】研究喀斯特馬鈴薯種植區(qū)土壤鎘的賦存形態(tài)及農(nóng)產(chǎn)品安全問題,為馬鈴薯種植區(qū)風險評估提供理論依據(jù)。【方法】以貴州省馬鈴薯主栽區(qū)域某縣馬鈴薯及相應(yīng)3種土壤類型(黃棕壤、棕壤和石灰土)為研究對象,采用CaCl2/BCR提取法測定土壤鎘的賦存形態(tài),以相關(guān)性和多因素方差分析探討土壤pH、有機質(zhì)、陽離子交換量(CEC)、總鎘及其交互作用對鎘賦存形態(tài)的影響;以富集系數(shù)和有效性評估土壤鎘累積能力;采用地質(zhì)累積指數(shù)法和潛在生態(tài)危害指數(shù)法比較土壤總鎘和有效鎘的生態(tài)風險?!窘Y(jié)果】研究區(qū)域土壤總鎘含量均值為1.394 mg/kg,變異系數(shù)為0.25;土壤鎘分級提取的5種賦存形態(tài)中,弱酸態(tài)鎘占總鎘比例最大,達28%~58%,可還原態(tài)鎘僅占4%~13%,其余3種形態(tài)介于二者之間??傛k對有效鎘、弱酸態(tài)鎘、可氧化態(tài)鎘和殘渣態(tài)鎘等4種鎘賦存形態(tài)均有極顯著影響(P<0.01,下同),總鎘與pH的交互作用對有效鎘影響大于總鎘和pH對有效鎘的主效影響,而pH與有機質(zhì)、pH與CEC交互影響弱于pH對有效鎘的主效影響,土壤鎘賦存形態(tài)的差異主要受總鎘和土壤pH、有機質(zhì)、CEC共同影響。種植區(qū)馬鈴薯塊莖鎘含量均值為0.086 mg/kg,富集系數(shù)為6%,有效性為15%。種植區(qū)土壤總鎘地質(zhì)累積指數(shù)(Igeo)為2~3級,表明土壤處于中度污染至強污染,綜合潛在生態(tài)風險指數(shù)(RI)為100.571~281.430,表明土壤處于輕微污染至中度污染。在馬鈴薯種植區(qū)3種土壤類型中,黃棕壤和棕壤pH均值分別為5.90和5.63,相對較低,但黃棕壤中有機質(zhì)和CEC含量最低,黃棕壤馬鈴薯的鎘有效性均值為17%,主要以弱酸態(tài)鎘形態(tài)存在,有利于馬鈴薯吸收富集,其富集系數(shù)均值為7%,且其潛在風險均值為212.371。反之,石灰土中pH較高,為6.87,石灰土和棕壤中有機質(zhì)和CEC較高,馬鈴薯富集能力較低,所以其鎘活性較低,潛在風險較低?!窘Y(jié)論】馬鈴薯種植區(qū)3種土壤類型中黃棕壤總鎘含量最高,活性最強,對應(yīng)馬鈴薯中鎘富集能力最強,其潛在風險較大,應(yīng)注重對該類型土壤在馬鈴薯種植區(qū)域的監(jiān)管。

    關(guān)鍵詞: 馬鈴薯;鎘;賦存形態(tài);有效性;土壤類型;風險評估

    中圖分類號: S155;X825? ? ? ? ? ? ? ? ? ? ? ?文獻標志碼: A 文章編號:2095-1191(2020)12-2901-10

    Abstract:【Objective】To provide reasonable theoretical basis for risk assessment of potato growing areas,speciation of cadmium in soil of karst potato origin and Safety of agricultural products were studied. 【Method】Potato and three soil types(yellow brown soil, brown soil and calcareous soil) in a main planting area of Guizhou were selected as the research objects. The speciations of cadmium in soil were determined by CaCl2 /BCR extraction method. The effects of soil pH,organic matter,cation exchange capacity(CEC),total cadmium and their interaction on the occurrence forms of cadmium were discussed by correlation and multi-factor analysis of variance. The accumulation of cadmium in soil was evaluated by enrichment coefficient and availability. The ecological risks of total cadmium and available cadmium in soil were compared by using geological accumulation index and potential ecological hazard index. 【Result】The average content of total cadmium in the soil was 1.394 mg/kg,and the coefficient of variation was 0.25. Among the five forms of soil cadmium,weak acid cadmium accounted for the largest proportion of total cadmium,accounting for 28%-58%,reducible cadmium only accounted for 4%-13%, and the other three forms were between the two. Total cadmium had extremely effects on effective cadmium, weak acid cadmium, oxidizable cadmium and residual cadmium(P<0.01, the same below). The intera-ction between total cadmium and pH had a greater effect on available cadmium than major effect of total cadmium and pH on available cadmium,while the interaction between pH and organic matter,pH and CEC was weaker than major effect of pH on available cadmium. The difference of occurrence forms of cadmium in soil was mainly influenced by total cadmium, soil pH, organic matter and CEC. The average content of cadmium in potato tubers was 0.086 mg/kg,the enrichment factor was 6%,and the availability was 15%. The total cadmium geological accumulation index(Igeo) of the soil in the potato planting area was 2-3,indicating that the soil was in moderate to strong pollution, and the comprehensive potential ecological risk index(RI) was 100.571-281.430,indicating that the soil was in mild to moderate pollution.In the three soil types of potato planting area,the average pH values of yellow brown soil and brown soil were 5.90 and 5.63respectively,which were relatively low,but the contents of organic matter and CEC in yellow brown soil were the lowest. The mean availability of cadmium in yellow brown soil was 17%,and it was mainly in the form of weak acid cadmium,which was conducive to the absorption and enrichment of potato. The mean enrichment factor was 7%,and the mean potential risk was 212.371. On the contrary,the pH value of calcareous soil was high(6.87). The organic matter and CEC in calcareous soil and brown soil were higher,and the accumulation ability of potato was low, so their cadmium activity was lower,and the potential risk was low. 【Conclusion】Among the three types of soil in potato planting area,yellow brown soil has the highest total cadmium content and the strongest activity,corresponding to the highest cadmium enrichment capacity in potato, which has great potential risk. Therefore,attention should be paid to the supervision of this type of soil in potato planting area.

    Key words: potato; cadmium; speciation; effectiveness; soiltype; risk assessment

    Foundation item: National Natural Science Foundation of China(U1612442); Construction Project of Science and Technology Infrastructure Platform in Guizhou(QKHPTRC〔2019〕5701)

    0 引言

    【研究意義】鎘是我國土壤最主要的污染物,給受污染區(qū)域的農(nóng)作物安全生產(chǎn)帶來極大的超標風險(徐佳慧等,2020),在鎘—土壤—生物這一復雜的動態(tài)鏈式反應(yīng)中,僅有部分土壤鎘能被生物吸收利用。因此,研究土壤鎘中各賦存形態(tài)及運用鎘的形態(tài)組分評估生物有效性,有利于進一步定量評估重金屬的危害性和治理修復重金屬污染土壤。【前人研究進展】化學提取法是常用于研究土壤鎘賦存形態(tài)的手段(Colzato et al.,2017),能以微量樣品展現(xiàn)土壤鎘的存在形態(tài)和比例。目前已有大量關(guān)于采用化學提取法分析水稻田、大白菜種植區(qū)等土壤鎘賦存形態(tài)和生物有效性的研究報道(袁波等,2011;葉興銀等,2018;彭云霄等,2019),同時深入探討酸堿度(劉丹青等,2013)、氧化還原條件(王成文等,2016)、多種金屬離子(臧彩云等,2017)、有機碳(楊美玉等,2018)等對鎘賦存形態(tài)的影響,并從土壤鎘與生物的反應(yīng)環(huán)境(Colzato et al.,2017)、不同形態(tài)鎘的化學性質(zhì)(Castillo-Michel et al.,2016)及不同形態(tài)鎘與植物鎘含量的相關(guān)性(袁波等,2011)等方面,論證了有效鎘能更精確表征土壤中鎘的生態(tài)風險。衡量不同pH條件下農(nóng)用地總鎘的農(nóng)產(chǎn)品質(zhì)量安全、農(nóng)作物生長或土壤生態(tài)環(huán)境風險,主要根據(jù)《土壤環(huán)境質(zhì)量農(nóng)用地土壤污染風險管控標準》(GB 15618—2018)中篩選值和管控值;評估土壤重金屬環(huán)境質(zhì)量常用內(nèi)梅羅(Nemerow)綜合污染指數(shù)法和地質(zhì)累積指數(shù)法。何騰兵等(2018)采用內(nèi)梅羅指數(shù)法評估貴州喀斯特地區(qū)施用10年豬糞肥的黃壤剖面土壤環(huán)境質(zhì)量及潛在風險,結(jié)果表明施用10年豬糞肥的表層土壤中鎘風險高于未施用豬糞肥的表層土壤;張俊輝等(2019)采用內(nèi)梅羅指數(shù)法評估植被葉子的鎘污染程度,結(jié)果表明部分植被樣點鉻、鎳和鎘污染嚴重超標,需加以重視?!颈狙芯壳腥朦c】目前環(huán)境質(zhì)量的評估主要針對環(huán)境中重金屬全量的生態(tài)風險,與有效鎘比較,總鎘無法精確表征土壤中鎘的生態(tài)風險。馬鈴薯為貴州省糧食的主要來源之一,貴州省近一半的耕地面積種植馬鈴薯,而有關(guān)馬鈴薯種植區(qū)土壤鎘有效性的研究鮮見報道。【擬解決的關(guān)鍵問題】針對貴州省馬鈴薯主種植區(qū)不同土壤類型開展鎘賦存形態(tài)研究,比較主種植區(qū)土壤總鎘和有效鎘與馬鈴薯鎘含量的相關(guān)性,以及土壤總鎘和有效鎘風險評估結(jié)果,為馬鈴薯種植區(qū)中不同土壤類型鎘污染的風險評估提供參考依據(jù)。

    1 材料與方法

    1. 1 研究區(qū)概況

    研究區(qū)域位于貴州省西北部某縣(東經(jīng)103°36′~104°45′,北緯26°36′~27°26′),屬亞熱帶季風性濕潤氣候區(qū),年平均氣溫10 ℃左右,年降水量890 mm,平均海拔2200 m。該區(qū)域成土母巖以砂頁巖和玄武巖為主,主要土壤類型是黃棕壤、棕壤和石灰土??h內(nèi)耕地面積22.2×104 ha(航測),其馬鈴薯年均種植面積11.0×104 ha。

    1. 2 樣品采集與制備

    通過前期充分調(diào)研該縣馬鈴薯種植用地中重金屬基本情況,根據(jù)種植力度最大的4個鎮(zhèn)內(nèi)單位面積主栽品種位點數(shù),按照點位權(quán)重,結(jié)合樣點地形等因素挑選30塊田地作為采樣對象,由該縣西南向東北方向土壤類型依次為石灰土、棕壤和黃棕壤,為耕種熟化的旱地。每塊田地采用五點法采集0~20 cm 耕層土樣,混勻后采用四分法反復取舍,保留2 kg左右的土樣裝入布袋,共采集土樣30份(黃棕壤15份,采樣點編號為HZ-1~HZ-15;棕壤7份,采樣點編號為Z-1~Z-7;石灰土8份,采樣點編號為S-1~S-8,下同)。室內(nèi)風干后四分法取縮分研磨10目和100目尼龍篩制備土壤樣品,其中10目用于pH測定;100目用于有機質(zhì)、總鎘和各賦存形態(tài)鎘含量測定。供試土壤理化性質(zhì)見表1。同時在對應(yīng)土壤樣品采集位點上協(xié)同采集馬鈴薯樣品30份運回實驗室,經(jīng)自來水洗凈后,再用超純水多次潤洗,然后切片,65 ℃下烘干至恒重,研磨粉碎,過40目尼龍篩用于馬鈴薯中鎘含量測定。

    1. 3 樣品測定

    土壤pH采用電位法測定(水土比為2.5∶1),有機質(zhì)含量采用重鉻酸鉀—外加熱法測定(魯如坤,1999)。土壤和植物樣品中總鎘含量運用高壓密閉消解法消解,使用ICP-MS(美國,Thermo Fisher X2)測定。土壤不同形態(tài)鎘采用CaCl2/BCR法測定(高煥方等,2017),采用石墨爐原子吸收法測定樣品中鎘含量,土壤樣品采用國家標準物質(zhì)GBW07443(GSS-5)和GBW07443(GSF-3)進行質(zhì)量控制(何騰兵等,2018)。

    1. 4 評估方法

    1. 4. 1 地質(zhì)累積指數(shù) 計算公式如下:

    Igeo=log2[Cn/1.5BEn]

    式中,Igeo為地質(zhì)累積指數(shù);Cn為樣品中元素n的濃度;BEn為環(huán)境背景濃度值,貴州土壤鎘背景值為 0.21 mg/kg(何騰兵等,2018);1.5為修正系數(shù),通常表征沉積特征、巖石地質(zhì)及其他影響。地質(zhì)累積指數(shù)分為7個級別:Igeo<0,表示無污染,污染級別是0級;0≤Igeo<1,表示無污染~中度污染,污染級別是1級;1≤Igeo<2,表示中度污染,污染級別是2級;2≤Igeo<3,表示中度污染~強污染,污染級別是3級;3≤Igeo<4,表示強污染,污染級別是4級;4≤Igeo<5,表示強污染~極強污染,污染級別是5級;Igeo≥5,表示極強污染,污染級別是6級(Müller,1969)。

    1. 4. 2 Hakanson潛在生態(tài)危害指數(shù) 計算公式如下:

    式中,RI為綜合潛在生態(tài)風險指數(shù);E[ir]為某一重金屬的潛在生態(tài)風險系數(shù);T[ir]為重金屬i的毒性影響系數(shù);C[if]為某一重金屬的污染系數(shù);C[is]為不同土壤層次重金屬i的實測值;C[in]為重金屬參照背景值,貴州土壤鎘背景值為0.21 mg/kg。鎘的毒性影響系數(shù)為30(何騰兵等,2018)。RI<150,表明生態(tài)危害程度為輕微;150≤RI<300,表明生態(tài)危害程度為中等;300≤RI<600,表明生態(tài)危害程度為強;RI≥600,表明生態(tài)危害程度為極強。

    1. 4. 3 鎘富集系數(shù)/有效性計算 鎘富集系數(shù)/有效性是馬鈴薯中鎘總含量與土壤中鎘總含量/有效含量之比,一定程度上反映土壤—馬鈴薯系統(tǒng)中鎘遷移的難易程度(袁波等,2011)。

    鎘富集系數(shù)=馬鈴薯中鎘總含量/土壤中鎘總含量×100

    鎘有效性=馬鈴薯中鎘總含量/土壤中鎘有效態(tài)含量×100

    1. 5 統(tǒng)計分析

    所有試驗數(shù)據(jù)采用Excel 2007進行整理分析,運用SPSS 20.0進行相關(guān)分析,用R語言軟件作圖。

    2 結(jié)果與分析

    2. 1 區(qū)域土壤鎘賦存形態(tài)分析結(jié)果

    從圖1可看出,研究區(qū)域土壤鎘的賦存形態(tài)具有明顯差異;土壤總鎘平均含量為1.394 mg/kg,是研究區(qū)域背景值的4~9倍,變異系數(shù)為0.25。采用CaCl2/BCR(回收率為86.40%~107.98%)分級提取的有效鎘、弱酸態(tài)鎘、可氧化態(tài)鎘、可還原態(tài)鎘和殘渣態(tài)鎘含量中,有效鎘和弱酸態(tài)鎘占總鎘含量比例最大,分別為26%~56%和28%~58%,有效鎘與弱酸態(tài)鎘的相關(guān)系數(shù)為0.97,可見單一提取劑CaCl2提取有效鎘與BCR提取的含量幾乎一致;其次是可氧化態(tài)鎘和殘渣態(tài)鎘,可還原態(tài)鎘所占比例最小,為4%~13%。以鎘的穩(wěn)定性來看,固定態(tài)鎘(可氧化態(tài)鎘、可還原態(tài)鎘和殘渣態(tài)鎘)的比值更高,達42%~72%;土壤4種形態(tài)鎘變異系數(shù)大小順序為:可還原態(tài)鎘(0.41)>可氧化態(tài)鎘(0.40)>殘渣態(tài)鎘(0.30)>弱酸態(tài)鎘(0.26)。說明樣點間4種形態(tài)鎘均處于中等強度變異。從有效鎘的比例來看,區(qū)域土壤鎘以固定態(tài)為主,鎘在土壤中化學性質(zhì)穩(wěn)定、遷移能力弱。

    從圖2可看出,不同土壤類型中鎘賦存形態(tài)同樣存在差異。黃棕壤總鎘均值為1.487 mg/kg,其弱酸態(tài)鎘、可氧化態(tài)鎘、可還原態(tài)鎘和殘渣態(tài)鎘含量所占比例分別為45%、27%、7%和17%;棕壤總鎘均值為1.267 mg/kg,其弱酸態(tài)鎘、可氧化態(tài)鎘、可還原態(tài)鎘和殘渣態(tài)鎘含量所占比例分別為51%、21%、9%和18%;石灰土總鎘均值為1.331 mg/kg,其弱酸態(tài)鎘、可氧化態(tài)鎘、可還原態(tài)鎘和殘渣態(tài)鎘含量所占比例分別為42%、30%、8%和15%。比較不同土壤類型 中鎘賦存形態(tài)比例,其差異主要體現(xiàn)在弱酸態(tài)鎘和可氧化態(tài)鎘,兩者所占比例是決定不同土壤鎘穩(wěn)定性的關(guān)鍵;種植區(qū)中黃棕壤總鎘含量最大,棕壤總鎘含量最小,但土壤中鎘活躍性最高,石灰土中鎘最穩(wěn)定。

    2. 2 土壤鎘賦存形態(tài)與理化性質(zhì)相關(guān)分析結(jié)果

    將土壤pH、有機質(zhì)、陽離子交換量(CEC)和總鎘與5種形態(tài)鎘進行多因素方差分析,由表2可知,總鎘對有效鎘、弱酸態(tài)鎘、可氧化態(tài)鎘和殘渣態(tài)鎘的主效應(yīng)最明顯;pH僅對有效鎘和弱酸態(tài)鎘產(chǎn)生極顯著(P<0.01)或顯著(P<0.05,下同)影響;有機質(zhì)和CEC不能顯著影響鎘賦存形態(tài)(P>0.05)。兩因子交互作用對鎘賦存形態(tài)的影響程度取決于自身主效影響,當兩因子均對鎘賦存形態(tài)產(chǎn)生顯著影響時,其交互作用會得到加強,例如PpH×總鎘

    pH有利于土壤固態(tài)鎘的分解與有效鎘的形成,但不同pH條件下,其對鎘賦存形態(tài)影響效果存在差異。由表3可知,pH為4.5~5.5、5.5~6.5和6.5~7.5時,pH與有效鎘的相關(guān)系數(shù)分別為-0.09、-0.44和-0.61,表明H+對土壤有效鎘的影響程度隨pH的升高而加強;此外,在pH為6.5~7.5時,土壤樣品均為石灰土,說明pH對石灰土中有效鎘影響高于黃棕壤和棕壤。土壤受總鎘影響的程度取決于土壤類型,其中棕壤中有效鎘受總鎘影響最大(R2=0.861),其次為石灰土(R2=0.434),而總鎘對黃棕壤中有效鎘影響不顯著(R2=0.090)??傮w而言,有效鎘易在總鎘高、pH趨于6.5的石灰土中高表達,在總鎘高、pH 5.5~6.5的棕壤和黃棕壤區(qū)域高表達。

    2. 3 種植區(qū)馬鈴薯富集系數(shù)、有效性和風險評估結(jié)果

    如表4所示,馬鈴薯塊莖鎘均值為0.086 mg/kg;富集系數(shù)均值為6%;CaCl2/BCR提取的有效性分別為8%~22%和7%~19%,其有效性均值分別為15%和14%。種植區(qū)土壤有效鎘與馬鈴薯塊莖中鎘含量相關(guān)性R有效鎘=0.714,線性擬合效果為R2有效鎘=0.51;總鎘與馬鈴薯塊莖中鎘含量相關(guān)性R總鎘=0.631,線性擬合效果為R2總鎘=0.40。土壤有效鎘與馬鈴薯塊莖中鎘含量的相關(guān)性、線性擬合效果均優(yōu)于總鎘。與富集系數(shù)相比,生物有效性更能準確指示環(huán)境污染對土壤的沖擊。雖然人類活動引起種植區(qū)土壤鎘含量高于背景值4~9倍,但馬鈴薯塊莖鎘含量均值小于《食品安全國家標準食品中污染物限量》(GB 2762—2017)中莖塊類蔬菜限定值(0.1 mg/kg),說明區(qū)域總鎘雖高,但遷移到馬鈴薯塊莖中的鎘含量較小。此外,黃棕壤、棕壤和石灰土中馬鈴薯塊莖鎘含量均值分別為0.099、0.068和0.078 mg/kg,變異系數(shù)分別為0.27、0.48和0.26。從馬鈴薯塊莖鎘含量均值來看,黃棕壤中馬鈴薯生物有效性最大,其次為石灰土,最小為棕壤。但從不同類型土壤有效鎘均值來看,棕壤有效鎘均值(0.646 mg/kg)高于石灰土(0.559 mg/kg)??赡苁怯绊懽厝乐墟k遷移到馬鈴薯的因素較石灰土復雜(棕壤中馬鈴薯塊莖鎘含量變異系數(shù)0.48>石灰土0.26),從而導致棕壤有效鎘含量高而生物有效性低。

    從圖3和表5可看出,馬鈴薯種植區(qū)土壤總鎘含量(0.704~1.970 mg/kg)均超出篩選值(0.3 mg/kg)(GB 15618—2018),pH 4.5~5.5,有6.67%樣品(占全量比)超過管控值(1.5 mg/kg);pH 5.5~6.5和pH 6.5~7.5,樣品均小于管控值(2.0和3.0 mg/kg)??傛k變異系數(shù)(0.25)較低,說明樣品間總鎘連續(xù)性較好,在0.139 mg/kg左右波動;而管控值是隨著pH升高進行梯試劃分,pH越高,管控值越高,使得高pH區(qū)域的土壤潛在風險偏低,表明石灰土的潛在風險較棕壤和黃棕壤低。以土壤有效鎘作為評價指標,樣品均在0.5 mg/kg左右波動并有6.67%樣品低于篩選值。采用Igeo和RI評價種植區(qū)土壤鎘生態(tài)風險:以總鎘作為評價指標,土壤總鎘的Igeo處于2~3級,說明土壤處于中度污染至強污染,土壤總鎘的RI為100.571~281.430,表明土壤處于輕微污染至中度污染;以有效鎘作為評價指標,土壤有效鎘的Igeo處于0~2級,為無污染至中度污染;土壤有效鎘的RI為38.429~128.140,表明土壤處于輕微污染。相比土壤總鎘,有效鎘的評估結(jié)果明顯較低,更準確地反映種植區(qū)馬鈴薯塊莖鎘含量(均值0.086 mg/kg)和生物有效性低(均值15%)的特點。

    3 討論

    貴州省馬鈴薯種植區(qū)土壤鎘賦存形態(tài)特征:一是種植區(qū)土壤鎘賦存形態(tài)以固定態(tài)鎘為主,由于黃棕壤、棕壤和石灰土受淋溶影響程度的不同,使得表層土壤鐵錳含量存在差異,導致鎘賦存形態(tài)在不同土壤中差異表達;二是總體而言,土壤總鎘高、pH低的區(qū)域有利于有效鎘集聚,土壤pH和總鎘高的區(qū)域有利于固定態(tài)鎘集聚。由于黃棕壤、棕壤和石灰土的土壤特點不同,其鎘賦存形態(tài)受土壤理化性質(zhì)的影響有差別,導致研究區(qū)域整體情況更加復雜。從區(qū)域土壤鎘賦存形態(tài)比例(42%~72%固定態(tài)鎘,28%~58%可提取態(tài)鎘)來看,土壤中鎘主要與土壤膠體結(jié)合(Khaokaew et al.,2011)。通過原子結(jié)合水鈉錳礦、高嶺石、有機質(zhì)等,或與硫化物、氫氧化物、碳酸結(jié)合物等造成固定態(tài)鎘,導致區(qū)域土壤有效鎘含量低于Hashimoto和Yamaguchi(2013)從水稻土中提取的有效鎘;黃棕壤、棕壤和石灰土中鎘賦存形態(tài)的差異主要體現(xiàn)在有效鎘和可氧化態(tài)鎘上,由于淋溶淀積作用,黃棕壤和棕壤中鐵錳含量隨深度而增加(寧婧等,2009),表層土壤鐵錳缺乏致使黃棕壤和棕壤表層土壤可氧化態(tài)鎘含量下降,有效鎘含量上升。種植區(qū)土壤鎘賦存形態(tài)與理化性質(zhì)相關(guān)性和多因素方差分析結(jié)果表明,總鎘和pH對鎘賦存形態(tài)影響最大,與袁波等(2011)對菜地土壤有效鎘的集聚結(jié)果、Colzato等(2017)利用化學提取法研究水稻土不同形態(tài)鎘的集聚趨勢一致??傛k為不同形態(tài)鎘存在的基礎(chǔ)。pH可改變鎘的溶解性,鎘—膠體復合物或與硫化物、氫氧化物、碳酸等結(jié)合的沉淀物均在酸性條件下易于溶解,從而在酸性條件下可溶解性鎘含量增加,導致游離鎘(有效鎘/弱酸態(tài)鎘)含量隨pH降低而升高,同時結(jié)合態(tài)鎘(可氧化)含量因pH降低而減少。土壤pH為4.5~6.5時,土壤pH與Cd2+線性擬合結(jié)果為R2=0.93(Cheng et al.,2017);土壤pH為4.0~9.0 時,土壤pH與可交換態(tài)鎘、可還原態(tài)鎘和殘渣態(tài)鎘含量呈極顯著線性相關(guān)(陳楠等,2018)等研究結(jié)果也表明土壤pH對土壤鎘賦存形態(tài)具有極大影響。與此同時,pH能影響鎘各賦存形態(tài)的存在形式,在堿性土壤中主要形態(tài)是CdCO3和Cd-土壤膠體,但在酸性土壤中CdS和Cd-土壤膠體占主要優(yōu)勢。

    種植區(qū)黃棕壤處于經(jīng)強度淋溶初級階段、為鹽基不飽和的弱富鋁化土壤;棕壤有較強烈的生物累積作用、土壤處于硅鋁化階段、有明顯的淋溶淀積、鐵錳含量隨深度而增加、鹽基不飽和、弱富鋁化等現(xiàn)象;石灰土發(fā)育于碳酸鹽巖母質(zhì),CaCO3和MgCO3含量在90%以上,其他成分主要是SiO2、Al2O3和Fe2O3等。在成土過程中,土壤脫鈣和覆鈣作用反復進行,使得土壤鹽基飽和度較高,碳酸鹽巖上發(fā)育的石灰土向游離鐵含量較高且CEC較低的鐵鋁土轉(zhuǎn)化(寧婧等,2009)。從土壤成分來看,石灰土成分較單一,90%以上為碳酸鹽,其鎘賦存形態(tài)90%左右為CdCO3(Colzato et al.,2017),因而研究的3種土壤中,石灰土的有效鎘與pH的相關(guān)系數(shù)(R=-0.610)大于棕壤(R=-0.142)和黃棕壤(R=-0.136)。由于石灰土中碳酸根、pH競爭與鎘結(jié)合;棕壤、黃棕壤因淋溶、脫硅富鋁化作用產(chǎn)生大量的鐵、鋁與鎘競爭H+,使喀斯特地區(qū)各鎘之間轉(zhuǎn)化變得更加復雜。土壤CEC和有機質(zhì)含量被視為與土壤表面中陰離子數(shù)量有關(guān),土壤表面中陰離子量越多,與游離鎘結(jié)合點越大。而種植區(qū)土壤有機質(zhì)的來源眾多,部分鎘元素與有機質(zhì)的外球面結(jié)構(gòu)通過非特異性結(jié)合,其結(jié)合鍵易受到外界的影響而斷裂,導致在酸性條件下增加土壤有效鎘含量。而內(nèi)球面的特異性結(jié)合會降低游離鎘含量,這種拮抗作用使得有機質(zhì)、CEC與鎘相關(guān)性不顯著(宋波和曾煒銓,2015;Bradl,2004)??傮w而言,種植區(qū)土壤鎘賦存形態(tài)以固定態(tài)鎘為主,碳酸和pH對鎘影響可能存在拮抗作用,土壤有機質(zhì)種類繁雜等因素使得pH、有機質(zhì)和CEC對鎘賦存形態(tài)影響較其他研究結(jié)果(曾煒銓等,2015)偏低。

    綜合生物有效性、富集系數(shù)、有效鎘和總鎘的區(qū)域風險評估結(jié)果,明確有效鎘能更精準表示污染物對土壤的沖擊(曾煒銓等,2015)。研究區(qū)域土壤主要以黃棕壤、棕壤和石灰土交錯分布,其碳酸鹽巖分布廣泛,土壤中鎘的平均背景含量遠高于我國其他地區(qū),具有明顯鎘地球化學高背景現(xiàn)象。鎘是分散元素之一,主要呈類質(zhì)同象伴生于閃鋅礦中,鉛鋅礦床一直是鎘的最主要來源,以石炭系、二疊系和寒武系地層分布區(qū)發(fā)育的土壤和水系沉積物中的鎘含量較高,即海相沉積碳酸鹽巖的巖石風化產(chǎn)物含鎘量較高(韓存亮,2012;龍家寰等,2014)。鎘主要集聚在土壤表層,在剖面垂直方向上鎘含量隨著土層深度增加而明顯降低(何騰兵等,2018)。本研究區(qū)為多年農(nóng)用地,磷肥與復合肥的施用量較大,但種植區(qū)馬鈴薯鎘含量均值(0.086 mg/kg)小于食品安全準則的限定值,可能是土壤淋溶、脫硅富鋁化作用產(chǎn)生大量離子與鎘進行鰲合、拮抗等作用,或者是馬鈴薯對污染物抵御(符東順等,2018)使得種植區(qū)馬鈴薯有效性僅為15%,較西紅柿(10.5倍)(周貴宇等,2016)、生菜(25倍)(陳展祥等,2018)和大田小白菜(1.08倍)(韓熙等,2018)等蔬菜中鎘的生物有效性低,可見種植區(qū)土壤中鎘的遷移能力較弱。對比馬鈴薯種植區(qū)土壤中有效鎘、總鎘含量與馬鈴薯塊莖中鎘含量相關(guān)系數(shù)及線性擬合度,結(jié)果為R有效鎘>R總鎘和R2有效鎘>R2總鎘。Colzato等(2017)認為無論鎘與土壤,還是土壤與植物體的相互作用均在水界面完成,Cd2+、CdOH-和CdCl+等有效鎘易被植物根毛吸收運往植物體內(nèi);Naidu等(1997)、葉宏萌等(2016)研究表明土壤中鎘不同形態(tài)活躍性大小順序為:有效態(tài)鎘>氧化態(tài)鎘>可還原態(tài)鎘>殘渣態(tài)鎘,活性越高越容易從土壤中釋放,從而被生物吸收,固定態(tài)鎘則是傾向于土壤內(nèi)部穩(wěn)定化;根據(jù)土壤鎘賦存形態(tài)的化學性質(zhì),Devesa和Vélez(2016)、葉興銀等(2018)認為土壤有效鎘是評估農(nóng)產(chǎn)品質(zhì)量安全、農(nóng)作物生長或土壤生態(tài)環(huán)境風險的關(guān)鍵。從不同土壤有效鎘及植物鎘含量綜合來看,黃棕壤有效鎘及其馬鈴薯生物有效性均最大,石灰土有效鎘含量最低而馬鈴薯鎘含量高于棕壤。說明種植區(qū)黃棕壤中馬鈴薯塊莖對有效鎘吸收較高,石灰土中馬鈴薯塊莖對鎘吸收較低,在棕壤表現(xiàn)厭鎘現(xiàn)象。從不同土壤植物鎘含量變異系數(shù)來看,影響棕壤中馬鈴薯鎘含量的因素更復雜,可能存在某些作用機制導致棕壤有效鎘含量高而生物有效性低。從不同土壤生物有效性和富集系數(shù)的比較,3種土壤中均表現(xiàn)出有效鎘與馬鈴薯塊莖中鎘含量相關(guān)性均高于全鎘的特點??傮w而言,雖然不同土壤中鎘賦存形態(tài)、遷移能力各具特色,但土壤鎘與馬鈴薯的反應(yīng)環(huán)境、土壤鎘賦存形態(tài)的化學性質(zhì)導致種植區(qū)有效鎘與馬鈴薯塊莖中鎘含量相關(guān)性優(yōu)于總鎘,以有效鎘表征馬鈴薯種植區(qū)土壤中鎘的遷移能力更科學客觀。種植區(qū)土壤總鎘在1.394 mg/kg左右波動,以農(nóng)用地標準劃分,土壤總鎘含量高于篩選值,且pH低的區(qū)域(黃棕壤、棕壤)較pH高的區(qū)域(石灰土)易超出管制值。以RI劃分,區(qū)域處于輕微污染至中度污染。以Igeo劃分,區(qū)域處于中度污染至強污染。但無論從種植區(qū)馬鈴薯鎘含量還是不同類型土壤中馬鈴薯鎘含量來看,樣品整體處于食品安全準則限定值內(nèi)。以總鎘作為馬鈴薯種植區(qū)潛在風險指標評估結(jié)果與實際相違背。從種植區(qū)土壤有效鎘、總鎘與馬鈴薯塊莖中鎘含量相關(guān)系數(shù)及線性擬合比較結(jié)果均表明以有效鎘作為貴州馬鈴薯種植區(qū)土壤風險指標更科學客觀。

    4 結(jié)論

    馬鈴薯種植區(qū)3種土壤類型中黃棕壤總鎘含量最高,活性最強,對應(yīng)馬鈴薯中鎘富集能力最強,其潛在風險較大,應(yīng)注重對該類型土壤馬鈴薯種植區(qū)域的監(jiān)管。

    參考文獻:

    陳楠,張昊,楊慧敏,魏祥東. 2018. 土壤pH對土壤鎘形態(tài)及稻米鎘積累的影響[J]. 湖南農(nóng)業(yè)大學學報(自然科學版),44(2):176-182. [Chen N,Zhang H,Yang H M,Wei X D. 2018. Effects of soil pH on soil cadmium formations and its accumulation in rice[J]. Journal of Hunan Agricultural University (Natural Sciences),44(2):176-182.]

    陳展祥,陳傳勝,陳衛(wèi)平,焦文濤. 2018. 凹凸棒石及其改性材料對土壤鎘生物有效性的影響與機制[J]. 環(huán)境科學,39(10):4744-4751. [Chen Z X,Chen C S,Chen W P,Jiao W T. 2018. Effect and mechanism of attapulgite and its modified materials on bioavailability of cadmium in soil[J]. Environmental Sciences,39(10):4744-4751.]

    符東順,何冠諦,付天嶺,何騰兵. 2018. 植物重金屬抗性基因挖掘和作用機制研究進展[J]. 江蘇農(nóng)業(yè)科學,46(10):18-23. [Fu D S,He G D,F(xiàn)u T L,He T B. 2018. Advances in the research of gene mining and mechanism of plant heavy metal resistance[J]. Jiangsu Agricultural Sciences,46(10):18-23.]

    高煥方,曹園城,何爐杰,許桂英,嚴歡,李聰,白薇揚,陳志. 2017. Tessier法和BCR法對比磷酸二氫鈉處置含鉛污染土壤形態(tài)分析[J]. 環(huán)境工程學報,11(10):5751-5756. [Gao H F,Cao Y C,He L J,Xu G Y,Yan H,Li C,Bai W Y,Chen Z. 2017. Speciation analysis of lead-contaminated soil treated with sodium dihydrogen phosphate using Tessier and BCR[J]. Chinese Journal of Environmental Engineering,11(10):5751-5756.]

    韓存亮. 2012. 地球化學異常與豬糞施用條件下土壤中鎘的分布、有效性與風險控制[D]. 北京:中國科學院研究生院. [Han C L. 2012. Soil cadmium distribution,availa-bility and risk control under conditions of geochemical anomaly or pig manure application[D]. Beijing:University of Chinese Academy of Sciences.]

    韓熙,張錫洲,余海英. 2018. 鈍化材料復配對土壤Cd生物有效性的影響[J]. 環(huán)境工程學報,12(10):2884-2891. [Han X,Zhang X Z,Yu H Y. 2018. Effects of compounds of passivation materials on bioavailability of Cd in soils[J]. Chinese Journal of Environmental Engineering,12(10):2884-2891.]

    何騰兵,黃會前,付天嶺,成劍波,石振情,林昌虎. 2018. 施用10年豬糞肥的黃壤剖面重金屬分布及風險評價[J]. 安全與環(huán)境學報,18(2):789-794. [He T B,Huang H Q,F(xiàn)u T L,Cheng J B,Shi Z Q,Lin C H. 2018. Analysis of the heavy metal risk content rate through the vertical profile of the yellow soil a decade late due to the swine manure application[J]. Journal of Safety and Environment,18(2):789-794.]

    劉丹青,陳雪,楊亞洲,王淑,李玉姣,胡浩,張春華,葛瀅. 2013. pH值和Fe、Cd處理對水稻根際及根表Fe、Cd吸附行為的影響[J]. 生態(tài)學報,33(14):4306-4314. [Liu D Q,Chen X,Yang Y Z,Wang S,Li Y J,Hu H,Zhang C H,Ge Y. 2013. Effects of pH,F(xiàn)e and Cd concentrations on the Fe and Cd adsorption in the rhizosphere and on the root surfaces of rice[J]. Acta Ecologica Sinica, 33(14):4306-4314.]

    龍家寰,劉鴻雁,劉方,朱恒亮,趙志鵬. 2014. 貴州省典型污染區(qū)土壤中鎘的空間分布及影響機制[J]. 土壤通報,45(5):1252-1259. [Long J H,Liu H Y,Liu F,Zhu H L,Zhao Z P. 2014. Spatial distribution and effect mechanisms of cadmium in soils in typical contaminated areas,Guizhou Province[J]. Chinese Journal of Soil Science,45(5):1252-1259.]

    魯如坤. 1999. 土壤農(nóng)業(yè)化學分析[M]. 北京:中國農(nóng)業(yè)科技出版社. [Lu R K. 1999. Soil argrochemistry analysis [M]. Beijing:China Agriculture Science and Technology Press.]

    寧婧,何騰兵,林昌虎,顏德宏,舒英格,陳祖擁. 2009. 黔北碳酸鹽巖發(fā)育土壤的發(fā)生特征與診斷特性及系統(tǒng)分類[J]. 貴州農(nóng)業(yè)科學,37(3):76-81. [Ning J,He T B,Lin C H,Yan D H,Shu Y G,Chen? Z Y. 2009. The genetic character and diagnostic characteristic of the soil derived from carbonate rock and its taxonomic classification in Northern Guizhou[J]. Guizhou Agricultural Sciences,37(3):76-81.]

    彭云霄,彭煒東,余江,鄧思維,賀玉龍,常家華,謝世前. 2019. 大田與盆栽條件下重金屬鎘賦存形態(tài)差異[J]. 江蘇農(nóng)業(yè)學報,35(6):1368-1373. [Peng Y X,Peng W D,Yu J,Deng S W,He Y L,Chang J H,Xie S Q. 2019. Di-fferences of heavy metal cadmium fractions in field-pot planting[J]. Jiangsu Journal of Agricultural Sciences,35(6):1368-1373.]

    宋波,曾煒銓. 2015. 土壤有機質(zhì)對鎘污染土壤修復的影響[J]. 土壤通報,46(4):1018-1024. [Song B,Zeng W Q. 2015. Effects of organic matter on the remediation of cadmium-contaminated soil-A review[J]. Chinese Journal of Soil Science, 46(4):1018-1024.]

    王成文,許模,張俊杰,丁慶忠,凌睿雯,郭蕾蕾. 2016. 土壤pH和Eh對重金屬鉻(Ⅵ)縱向遷移及轉(zhuǎn)化的影響[J]. 環(huán)境工程學報,10(10):6035-6041. [Wang C W,Xu M,Zhang J J,Ding Q Z,Ling R W,Guo L L. 2016. Influence of soils pH and Eh on vertical migration and transformation of Cr(Ⅵ)[J]. Chinese Journal of Environmental Engineering,10(10):6035-6041.]

    徐佳慧,王萌,張潤,吳玲玲. 2020. 土壤鎘污染的生物毒性研究進展[J]. 生態(tài)毒理學報,15(5):82-91. [Xu J H,Wang M,Zhang R,Wu L L. 2020. Toxicity of cadmium pollution in soil to organisms:A review[J]. Asian Journal of Ecotoxicology,15(5):82-91.]

    楊美玉,羅力,陳衛(wèi)鋒,倪進治,魏然. 2018. 生物炭對污泥施用土壤中三氯生生物有效性和微生物群落結(jié)構(gòu)的影響[J]. 環(huán)境科學學報,38(12):4789-4795. [Yang M Y,Luo L,Chen W F,Ni J Z,Wei R. 2018. Effects of biochars on the bioavailability of triclosan and microbial community structure in the soil amended with sewage sludge[J]. Acta Scientiae Circumstantiae,38(12):4789-4795.]

    葉宏萌,李國平,鄭茂鐘,袁旭音,常雪花. 2016. 武夷山茶園土壤中五種重金屬的化學形態(tài)和生物有效性[J]. 環(huán)境化學,35(10):2071-2078. [Ye H M,Li G P,Zheng M Z,Yuan X Y,Chang X H. 2016. Speciation and bioavailabi-lity of five toxic heavy metals in the tea garden soils of Wuyishan[J]. Environmental Chemistry,35(10):2071-2078.]

    葉興銀,張衛(wèi),龍精華,常文靜,曾輝. 2018. 人工納米顆粒輸入對稻田土壤Cd形態(tài)轉(zhuǎn)化及生物有效性的影響[J]. 環(huán)境工程學報,12(12):3426-3432. [Ye X Y, Zhang W, Long J H, Chang W J, Zeng H. 2018. Impact of engineered nanoperticles on Cd chemical speciation and bioavailability in paddy soil[J]. Chinese Journal of Environmental Engineering,12(12):3426-3432.]

    袁波,傅瓦利,藍家程,張婷,彭景濤. 2011. 菜地土壤鉛、鎘有效態(tài)與生物有效性研究[J]. 水土保持學報,25(5):130-134. [Yuan B,F(xiàn)u W L,Lan J C,Zhang T,Peng J T. 2011. Study on the available and bioavailability of lead and cadmium in soil of vegetable plantation[J]. Journal of Soil and Water Conservation,25(5):130-134.]

    臧彩云,沈袁玲,顧敏京,黃顧林,左文剛,柏彥超,單玉華. 2017. 零價鐵對灘涂土壤重金屬鎘形態(tài)及有效性的影響[J]. 江蘇農(nóng)業(yè)科學,45(18):249-251. [Zang C Y,Shen Y L,Gu M J,Huang G L,Zuo W G,Bo Y C,Shan Y H. 2017. Effects of zero-valent iron on the morphology and availability of cadmium,a heavy metal,in intertidal soils[J]. Jiangsu Agricultural Sciences ,45(18):249-251.]

    曾煒銓,宋波,袁立竹,黃宇妃,伏鳳艷. 2015. 改良劑對廣西環(huán)江強酸鉛鋅污染土壤的修復作用[J]. 環(huán)境科學,36(6):2306-2313. [Zeng W Q,Song B,Yuan L Z,Huang Y F,F(xiàn)u F Y. 2015. Effects of remedies on the remediation of typical Pb and Zn-contaminated soil in Huan-jiang,Guangxi[J]. Environmental Science,36(6):2306-2313.]

    張俊輝,林青,姜珊,劉濱,李東興,王彥虎. 2019. 寶雞市植被葉子重金屬分布規(guī)律及生態(tài)風險評價[J]. 環(huán)境科學,40(2):978-986. [Zhang J H,Lin Q,Jiang S,Liu B,Li D X,Wang Y H. 2019.? Analysis of heavy metal pollution and ecological risk assessment on vegetation leaves in Baoji City[J]. Environmental Science,40(2):978-986.]

    周貴宇,姜慧敏,楊俊誠,張建峰,張水勤,梁雷. 2016. 幾種有機物料對設(shè)施菜田土壤Cd、Pb生物有效性的影響[J]. 環(huán)境科學,37(10):4011-4019. [Zhou G Y,Jiang H M,Yang J C,Zhang J F,Zhang S Q,Liang L. 2016. Effects of different organic materials on bio-availability of Cd,Pb in a contaminated greenhouse soil[J]. Environmental Sciences,37(10):4011-4019.]

    Bradl H B. 2004. Adsorption of heavy metal ions on soils and soils constituents[J]. Journal of Colloid & Interface Science,277(1):1-18.

    Castillo-Michel H A,Larue C,del Real A E P,Cotte M,Sarrret G. 2016. Practical review on the use of synchrotron based micro- and nano- X-ray fluorescence mapping and X-ray absorption spectroscopy to investigate the interactions between plants and engineered nanomaterials[J]. Plant Physiology and Biochemistry,110(10):13-32.

    Cheng M ,Wang A,Tang C. 2017. Ammonium-based fertili-zers enhance Cd accumulation in Carpobrotus rossii grown in two soils differing in pH[J]. Chemosphere,188(11):689-696.

    Colzato M,Kamogawa M Y,Carvalho H W P,Alleoni L R F,Hesterberg D. 2017. Temporal changes in cadmium speciation in Brazilian soils evaluated using Cd LIII-edge XANES and chemical fractionation[J]. Journal of Environmental Quality,46(6):1206-1214.

    Devesa V,Vélez D. 2016. Cadmium:Properties and determination[J]. Encyclopedia of Food & Health,46(6):543-549.

    Hashimoto Y,Yamaguchi N. 2013. Chemical speciation of cadmium and sulfur K-edge XANES spectroscopy in flooded paddy soils amended with zerovalent iron[J]. Soil Science Society of America Journal,77(4):1189-1198.

    Khaokaew S,Chaney R L,Landrot G,Ginder-Vogel M,Sparks D L. 2011. Speciation and release kinetics of cadmium in an alkaline paddy soil under various flooding periods and draining conditions[J]. Environmental Scien-ce & Technology,45(10):4249-4255.

    Müller G. 1969. Index of geoaccumulation in sediments of the Rhine River[J]. GeoJournal,2(3):109-118.

    Naidu R,Kookana R S,Sumner M E,Harter R D,Tiller K G. 1997. Cadmium sorption and transport in variable charge soils:A review[J]. Journal of Environmental Quality,26(3):602-617.

    (責任編輯 羅 麗)

    猜你喜歡
    風險評估馬鈴薯有效性
    馬鈴薯有功勞
    如何提高英語教學的有效性
    甘肅教育(2020年6期)2020-09-11 07:45:28
    制造業(yè)內(nèi)部控制有效性的實現(xiàn)
    提高家庭作業(yè)有效性的理論思考
    甘肅教育(2020年12期)2020-04-13 06:24:56
    定邊馬鈴薯
    陜西畫報(2016年1期)2016-12-01 05:35:30
    我國養(yǎng)老保險基金投資運營的風險評估
    時代金融(2016年23期)2016-10-31 13:25:28
    天然氣開采危險有害因素辨識及風險評估管理
    大眾理財顧問(2016年8期)2016-09-28 13:55:43
    中小企業(yè)財務(wù)管理問題研究
    胖胖的馬鈴薯
    色视频在线一区二区三区| 亚洲国产毛片av蜜桃av| 亚洲成人av在线免费| 亚洲欧美精品专区久久| 亚洲精品国产成人久久av| 王馨瑶露胸无遮挡在线观看| 亚洲精品国产av蜜桃| 色5月婷婷丁香| 国产成人精品久久久久久| 日韩制服骚丝袜av| av卡一久久| 免费看日本二区| 中文资源天堂在线| 婷婷色麻豆天堂久久| 精品少妇久久久久久888优播| 久久久久国产精品人妻一区二区| 免费人成在线观看视频色| 免费观看在线日韩| 亚洲av成人精品一二三区| 亚洲精品久久午夜乱码| 香蕉精品网在线| 国产毛片在线视频| 精华霜和精华液先用哪个| 偷拍熟女少妇极品色| 欧美xxxx性猛交bbbb| 99热这里只有是精品50| 亚洲第一区二区三区不卡| 色婷婷av一区二区三区视频| 欧美 亚洲 国产 日韩一| 高清在线视频一区二区三区| 肉色欧美久久久久久久蜜桃| 一本—道久久a久久精品蜜桃钙片| 六月丁香七月| 高清毛片免费看| 国产成人精品福利久久| 色婷婷av一区二区三区视频| 人妻人人澡人人爽人人| 国产免费又黄又爽又色| 伊人亚洲综合成人网| 国产成人精品婷婷| 下体分泌物呈黄色| 中文资源天堂在线| 少妇被粗大猛烈的视频| 国模一区二区三区四区视频| 亚洲av在线观看美女高潮| 国产精品.久久久| 婷婷色综合大香蕉| av国产久精品久网站免费入址| 久久久久久久国产电影| 成人亚洲欧美一区二区av| 亚洲欧洲国产日韩| 亚洲成人手机| 亚洲av成人精品一区久久| 夜夜看夜夜爽夜夜摸| 久久久久久伊人网av| 成人无遮挡网站| 又爽又黄a免费视频| 久久国产精品大桥未久av | 国产视频内射| 人妻制服诱惑在线中文字幕| 久久久久久久久大av| 亚洲国产最新在线播放| 少妇精品久久久久久久| 我的老师免费观看完整版| 国产亚洲91精品色在线| 亚洲欧美日韩另类电影网站| 99国产精品免费福利视频| 少妇人妻一区二区三区视频| 校园人妻丝袜中文字幕| 又爽又黄a免费视频| 黑丝袜美女国产一区| www.色视频.com| av国产久精品久网站免费入址| 涩涩av久久男人的天堂| 成人毛片a级毛片在线播放| 啦啦啦在线观看免费高清www| 国内精品宾馆在线| 国产精品人妻久久久久久| 婷婷色麻豆天堂久久| 国产日韩一区二区三区精品不卡 | 纯流量卡能插随身wifi吗| av有码第一页| 亚洲伊人久久精品综合| 丝袜脚勾引网站| 亚洲人成网站在线观看播放| 精品久久久久久久久亚洲| 啦啦啦中文免费视频观看日本| 特大巨黑吊av在线直播| 国产精品久久久久久久电影| 男女国产视频网站| 18禁动态无遮挡网站| 国产黄片美女视频| 国产男人的电影天堂91| 三级国产精品片| 人妻一区二区av| 成人黄色视频免费在线看| 亚洲,一卡二卡三卡| 免费黄网站久久成人精品| 老司机影院毛片| 18禁裸乳无遮挡动漫免费视频| 亚洲四区av| 男女无遮挡免费网站观看| 天天操日日干夜夜撸| 欧美国产精品一级二级三级 | 成人国产麻豆网| 国语对白做爰xxxⅹ性视频网站| av在线老鸭窝| 免费高清在线观看视频在线观看| 亚洲婷婷狠狠爱综合网| 人妻 亚洲 视频| 最新的欧美精品一区二区| 久久人人爽av亚洲精品天堂| av在线app专区| 免费观看在线日韩| 久久久精品94久久精品| 日本黄色片子视频| 国产黄片美女视频| 国产精品免费大片| 欧美日韩亚洲高清精品| 日本黄大片高清| 国产亚洲午夜精品一区二区久久| 中文字幕亚洲精品专区| 国产av精品麻豆| 亚州av有码| 黄色视频在线播放观看不卡| 免费看光身美女| 少妇丰满av| 精品卡一卡二卡四卡免费| 久久久a久久爽久久v久久| 国产女主播在线喷水免费视频网站| 男女国产视频网站| 久久精品国产鲁丝片午夜精品| 寂寞人妻少妇视频99o| 成人特级av手机在线观看| 丰满饥渴人妻一区二区三| 国产欧美亚洲国产| 久久久国产一区二区| 狂野欧美激情性bbbbbb| 久久久久国产网址| 亚洲av.av天堂| 曰老女人黄片| 少妇高潮的动态图| 有码 亚洲区| 热re99久久精品国产66热6| 黄色一级大片看看| 日韩大片免费观看网站| 国产精品久久久久久精品古装| 少妇高潮的动态图| 久久久精品94久久精品| 久久精品国产亚洲av天美| 男女免费视频国产| av在线播放精品| 老熟女久久久| 男女免费视频国产| 人人妻人人添人人爽欧美一区卜| 免费黄网站久久成人精品| 午夜福利在线观看免费完整高清在| 国产精品欧美亚洲77777| 最近的中文字幕免费完整| 亚洲在久久综合| 久久久久久伊人网av| 国产精品一区www在线观看| 中国美白少妇内射xxxbb| 又大又黄又爽视频免费| 黄色视频在线播放观看不卡| 黄色日韩在线| 噜噜噜噜噜久久久久久91| 亚洲国产成人一精品久久久| 伊人亚洲综合成人网| 日本爱情动作片www.在线观看| 又爽又黄a免费视频| 男的添女的下面高潮视频| 天美传媒精品一区二区| 久久国产精品男人的天堂亚洲 | 欧美精品一区二区大全| 黄色毛片三级朝国网站 | 欧美性感艳星| 亚洲第一av免费看| 美女主播在线视频| 亚洲人与动物交配视频| 欧美日韩综合久久久久久| 色哟哟·www| 久热久热在线精品观看| 亚洲精品aⅴ在线观看| 一级毛片久久久久久久久女| 亚洲三级黄色毛片| 日本午夜av视频| 高清欧美精品videossex| 18+在线观看网站| 99久久精品热视频| 色视频www国产| 女人久久www免费人成看片| 久久久国产一区二区| 国产精品久久久久久久久免| 国产免费福利视频在线观看| 18禁在线播放成人免费| 国产免费又黄又爽又色| 搡女人真爽免费视频火全软件| 精品人妻熟女av久视频| 久久久久久久久久人人人人人人| 国产日韩一区二区三区精品不卡 | 80岁老熟妇乱子伦牲交| 国产精品秋霞免费鲁丝片| 欧美精品亚洲一区二区| 国产av码专区亚洲av| 精品久久国产蜜桃| 99热这里只有精品一区| 亚洲精品,欧美精品| 能在线免费看毛片的网站| 蜜桃在线观看..| 又大又黄又爽视频免费| 麻豆乱淫一区二区| 两个人的视频大全免费| 国产精品久久久久久av不卡| av视频免费观看在线观看| 成人国产麻豆网| 曰老女人黄片| 欧美日韩综合久久久久久| 激情五月婷婷亚洲| 一个人免费看片子| 精品国产乱码久久久久久小说| 亚洲欧美中文字幕日韩二区| 免费av不卡在线播放| 男女国产视频网站| 成人综合一区亚洲| 少妇人妻久久综合中文| 欧美3d第一页| 午夜日本视频在线| 高清欧美精品videossex| 国产欧美日韩精品一区二区| 女人精品久久久久毛片| 久热久热在线精品观看| 亚洲精品,欧美精品| 一区二区三区乱码不卡18| 中文资源天堂在线| 精品午夜福利在线看| 在线观看www视频免费| 成年美女黄网站色视频大全免费 | 亚洲成人一二三区av| av.在线天堂| 美女xxoo啪啪120秒动态图| 免费黄频网站在线观看国产| 人人妻人人澡人人看| 黄色视频在线播放观看不卡| 高清午夜精品一区二区三区| 热re99久久国产66热| 日韩欧美一区视频在线观看 | 精品一品国产午夜福利视频| 亚洲av.av天堂| 国产精品一二三区在线看| 大陆偷拍与自拍| 亚洲av电影在线观看一区二区三区| 伊人久久国产一区二区| 哪个播放器可以免费观看大片| 成年人午夜在线观看视频| 人妻制服诱惑在线中文字幕| 欧美少妇被猛烈插入视频| 亚洲美女搞黄在线观看| 中文欧美无线码| 国产精品熟女久久久久浪| 免费大片黄手机在线观看| 精品久久久久久久久av| 日本91视频免费播放| 国产精品无大码| 热99国产精品久久久久久7| 免费看光身美女| 精品久久久久久久久av| 国产色婷婷99| 99久久精品热视频| 免费黄网站久久成人精品| 这个男人来自地球电影免费观看 | 丝袜喷水一区| 亚洲av成人精品一二三区| 麻豆乱淫一区二区| 久久久久久久久久久久大奶| 欧美 亚洲 国产 日韩一| 黄色配什么色好看| 99九九线精品视频在线观看视频| 国产精品国产av在线观看| 99热网站在线观看| 欧美成人午夜免费资源| 少妇人妻一区二区三区视频| 日韩一区二区三区影片| 久久久亚洲精品成人影院| 能在线免费看毛片的网站| 不卡视频在线观看欧美| 久久久精品94久久精品| 老女人水多毛片| 免费播放大片免费观看视频在线观看| 国产无遮挡羞羞视频在线观看| 国产免费福利视频在线观看| 高清欧美精品videossex| 九草在线视频观看| 久久综合国产亚洲精品| 欧美高清成人免费视频www| 爱豆传媒免费全集在线观看| 国产日韩欧美亚洲二区| 色婷婷av一区二区三区视频| 免费黄网站久久成人精品| 亚洲欧美精品专区久久| 99热全是精品| 在线观看av片永久免费下载| 在线 av 中文字幕| 久久人妻熟女aⅴ| 18禁在线无遮挡免费观看视频| 麻豆成人av视频| 国产男女内射视频| 美女中出高潮动态图| 美女国产视频在线观看| 亚洲美女搞黄在线观看| 搡女人真爽免费视频火全软件| 色视频www国产| 91午夜精品亚洲一区二区三区| 人妻人人澡人人爽人人| 国产男女超爽视频在线观看| 精品少妇内射三级| 亚洲精品自拍成人| 狂野欧美激情性bbbbbb| 日韩av在线免费看完整版不卡| 国产精品熟女久久久久浪| 久久久久久久久久成人| 视频中文字幕在线观看| 亚洲在久久综合| 伦理电影免费视频| 国产在线男女| 色哟哟·www| 99热这里只有是精品50| 亚洲国产欧美在线一区| 亚洲av在线观看美女高潮| 99热国产这里只有精品6| 国产日韩欧美视频二区| 最新的欧美精品一区二区| 日韩 亚洲 欧美在线| 久久人妻熟女aⅴ| 王馨瑶露胸无遮挡在线观看| 亚洲成人av在线免费| 久久6这里有精品| 国产在线男女| 成人影院久久| 成年女人在线观看亚洲视频| 色哟哟·www| 精品久久久噜噜| 一本—道久久a久久精品蜜桃钙片| 成年av动漫网址| 日本色播在线视频| 久久久久久久大尺度免费视频| 最新的欧美精品一区二区| 国产精品久久久久久精品电影小说| 亚洲电影在线观看av| av国产久精品久网站免费入址| 国产白丝娇喘喷水9色精品| 最近中文字幕高清免费大全6| 午夜免费男女啪啪视频观看| 国产高清国产精品国产三级| 久久久精品免费免费高清| 中文精品一卡2卡3卡4更新| 91成人精品电影| 亚洲欧美成人综合另类久久久| 久久久久久久大尺度免费视频| 中文欧美无线码| 韩国av在线不卡| 亚洲精品日韩在线中文字幕| 亚洲精品第二区| 天天躁夜夜躁狠狠久久av| 欧美日韩av久久| 亚洲av日韩在线播放| 国产精品嫩草影院av在线观看| 成人免费观看视频高清| 午夜福利视频精品| 精品国产露脸久久av麻豆| 精品卡一卡二卡四卡免费| 女人久久www免费人成看片| 黄色毛片三级朝国网站 | 日本欧美国产在线视频| 高清黄色对白视频在线免费看 | 久久99精品国语久久久| 亚洲成人一二三区av| 国产黄频视频在线观看| 欧美 日韩 精品 国产| 一区在线观看完整版| 嘟嘟电影网在线观看| 黄片无遮挡物在线观看| 中文字幕制服av| 免费观看a级毛片全部| 欧美3d第一页| 国产极品粉嫩免费观看在线 | 国产成人精品一,二区| 一本—道久久a久久精品蜜桃钙片| 精品人妻一区二区三区麻豆| 久久精品国产自在天天线| 午夜激情福利司机影院| 久久久国产精品麻豆| 大码成人一级视频| 人妻系列 视频| 免费观看av网站的网址| 国产精品女同一区二区软件| 99热这里只有精品一区| 人人妻人人澡人人看| 亚洲av在线观看美女高潮| 一区二区av电影网| 中文在线观看免费www的网站| 日韩大片免费观看网站| 夜夜爽夜夜爽视频| 亚洲欧洲精品一区二区精品久久久 | 黄色视频在线播放观看不卡| 午夜日本视频在线| 在线观看www视频免费| 日本色播在线视频| 久久精品国产鲁丝片午夜精品| 99久久精品热视频| 欧美一级a爱片免费观看看| 午夜激情久久久久久久| 日韩在线高清观看一区二区三区| 人人澡人人妻人| 亚洲国产日韩一区二区| 嫩草影院入口| 大话2 男鬼变身卡| 亚洲精品成人av观看孕妇| 在线观看免费日韩欧美大片 | 午夜福利,免费看| 久久99热这里只频精品6学生| 免费大片黄手机在线观看| 久久国产亚洲av麻豆专区| 色网站视频免费| 热99国产精品久久久久久7| 亚洲av日韩在线播放| 亚洲欧美一区二区三区国产| 久久99精品国语久久久| 99热这里只有精品一区| 久久久久视频综合| 亚洲丝袜综合中文字幕| 青青草视频在线视频观看| 人妻制服诱惑在线中文字幕| 国语对白做爰xxxⅹ性视频网站| 肉色欧美久久久久久久蜜桃| 狂野欧美白嫩少妇大欣赏| 久久久午夜欧美精品| 热re99久久精品国产66热6| 国产一区二区在线观看日韩| 精品人妻一区二区三区麻豆| 又粗又硬又长又爽又黄的视频| 国产午夜精品一二区理论片| 视频中文字幕在线观看| 成人亚洲欧美一区二区av| 99久久精品国产国产毛片| 精品人妻偷拍中文字幕| 三级国产精品片| 精品国产乱码久久久久久小说| 亚洲图色成人| 九色成人免费人妻av| 国产午夜精品一二区理论片| 国精品久久久久久国模美| 人人妻人人澡人人爽人人夜夜| 日韩欧美精品免费久久| 天天躁夜夜躁狠狠久久av| 久久午夜福利片| 欧美日韩综合久久久久久| 99久久精品国产国产毛片| 国内少妇人妻偷人精品xxx网站| 免费黄色在线免费观看| 69精品国产乱码久久久| 丁香六月天网| 国国产精品蜜臀av免费| 国产成人freesex在线| 美女内射精品一级片tv| 日本黄大片高清| 免费在线观看成人毛片| 69精品国产乱码久久久| 久久久久久久亚洲中文字幕| 男女国产视频网站| 国产在线免费精品| 国产精品久久久久成人av| 日韩av在线免费看完整版不卡| 亚洲国产欧美日韩在线播放 | 国产av一区二区精品久久| 一级二级三级毛片免费看| av有码第一页| 一本一本综合久久| 亚洲真实伦在线观看| 久久久久久久国产电影| 又粗又硬又长又爽又黄的视频| 国产成人午夜福利电影在线观看| 在线观看免费高清a一片| 天堂中文最新版在线下载| 精品人妻熟女av久视频| 少妇被粗大猛烈的视频| 国产高清三级在线| 偷拍熟女少妇极品色| 夫妻性生交免费视频一级片| 下体分泌物呈黄色| 日韩中字成人| 97在线人人人人妻| 有码 亚洲区| 欧美+日韩+精品| 欧美日韩国产mv在线观看视频| 天天操日日干夜夜撸| 亚洲精品中文字幕在线视频 | 国产精品蜜桃在线观看| 中国三级夫妇交换| 久久亚洲国产成人精品v| 亚洲高清免费不卡视频| 91久久精品电影网| 国产精品久久久久久久久免| 一本色道久久久久久精品综合| 国产精品人妻久久久久久| 久久6这里有精品| 精品一区二区三区视频在线| 街头女战士在线观看网站| av黄色大香蕉| 亚洲精品国产色婷婷电影| 黄色一级大片看看| 国产老妇伦熟女老妇高清| 在线观看免费日韩欧美大片 | av有码第一页| 在线 av 中文字幕| 午夜av观看不卡| 国产亚洲av片在线观看秒播厂| 久久久精品94久久精品| 亚洲成色77777| 中文乱码字字幕精品一区二区三区| 国产一区二区三区av在线| tube8黄色片| 精品一区在线观看国产| 波野结衣二区三区在线| 国产中年淑女户外野战色| 亚洲怡红院男人天堂| 精品一品国产午夜福利视频| 日日摸夜夜添夜夜爱| 午夜福利在线观看免费完整高清在| 国产亚洲一区二区精品| 亚洲精品乱久久久久久| 在线精品无人区一区二区三| 亚洲欧美清纯卡通| 国产一区二区三区av在线| 国产精品麻豆人妻色哟哟久久| 男女国产视频网站| 亚洲成人av在线免费| 国产精品国产三级国产专区5o| 一区二区av电影网| h视频一区二区三区| 精品少妇久久久久久888优播| 午夜福利在线观看免费完整高清在| 青春草国产在线视频| 色吧在线观看| 日韩成人伦理影院| 亚洲精品国产av成人精品| 2022亚洲国产成人精品| 高清在线视频一区二区三区| 国产亚洲最大av| 精品国产一区二区三区久久久樱花| 最近的中文字幕免费完整| 极品人妻少妇av视频| 国产午夜精品一二区理论片| 午夜老司机福利剧场| 亚洲精品日韩在线中文字幕| 精品亚洲成a人片在线观看| 精品一区在线观看国产| 欧美区成人在线视频| 亚洲精品成人av观看孕妇| 亚洲av男天堂| 在线播放无遮挡| 亚洲欧美清纯卡通| 国产午夜精品久久久久久一区二区三区| 91成人精品电影| 人人妻人人添人人爽欧美一区卜| 国产成人精品无人区| 久久久久久久大尺度免费视频| 亚洲欧洲日产国产| 亚洲经典国产精华液单| 麻豆精品久久久久久蜜桃| 亚洲一级一片aⅴ在线观看| 午夜福利,免费看| 亚洲av免费高清在线观看| 少妇人妻一区二区三区视频| 日韩精品免费视频一区二区三区 | 亚洲伊人久久精品综合| 免费看不卡的av| 中文字幕精品免费在线观看视频 | 久久久亚洲精品成人影院| 丝瓜视频免费看黄片| 久久午夜综合久久蜜桃| 成年美女黄网站色视频大全免费 | 国产在线免费精品| 水蜜桃什么品种好| 亚洲国产av新网站| 人人妻人人爽人人添夜夜欢视频 | 最近中文字幕2019免费版| 三级国产精品片| 亚洲国产欧美在线一区| 女人久久www免费人成看片| 免费在线观看成人毛片| 中国三级夫妇交换| 人妻制服诱惑在线中文字幕| 少妇人妻精品综合一区二区| 成人国产av品久久久| 国产精品一区二区性色av| 亚洲国产成人一精品久久久| 亚洲精品国产色婷婷电影| 免费观看在线日韩| 亚洲精品中文字幕在线视频 | 免费黄频网站在线观看国产| 啦啦啦中文免费视频观看日本| 欧美精品人与动牲交sv欧美| 亚洲国产精品国产精品| 午夜免费鲁丝| 久久久亚洲精品成人影院| 精品人妻一区二区三区麻豆| 嫩草影院新地址| 91午夜精品亚洲一区二区三区| 免费黄频网站在线观看国产| 免费观看性生交大片5| 亚洲精品视频女| 午夜影院在线不卡| 六月丁香七月| 国产亚洲5aaaaa淫片|