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

    土壤Mehlich-3可浸提態(tài)鎳對(duì)大麥根伸長的毒性

    2018-08-23 05:45:20朱廣云李菊梅李合蓮馬義兵
    中國環(huán)境科學(xué) 2018年8期
    關(guān)鍵詞:大麥外源毒性

    朱廣云,蔣 寶,李菊梅,李合蓮,馬義兵

    ?

    土壤Mehlich-3可浸提態(tài)鎳對(duì)大麥根伸長的毒性

    朱廣云1,2,蔣 寶2,李菊梅2,李合蓮1,馬義兵3*

    (1.濟(jì)南大學(xué)資源與環(huán)境學(xué)院,山東 濟(jì)南 250022;2.中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所 北京 100081;3.廣東省生態(tài)環(huán)境技術(shù)研究所,廣東 廣州 510650)

    本研究選擇我國具有代表性的11種外源添加鎳(水溶性鎳鹽)污染土壤樣品,在模擬田間人工降水(淋洗和非淋洗)處理后,研究了Mehlich-3單次或連續(xù)3次浸提態(tài)鎳對(duì)大麥根伸長的毒性.結(jié)果發(fā)現(xiàn): 結(jié)合重要影響因素土壤pH值,Mehlich-3可浸提態(tài)鎳和外源添加鎳具有顯著相關(guān)關(guān)系;基于Mehlich-3可浸提態(tài)鎳的植物毒性閾值也顯著受土壤性質(zhì)影響(半抑制濃度可相差38倍以上).通過回歸分析發(fā)現(xiàn)土壤性質(zhì)對(duì)Mehlich-3浸提鎳的生物毒性的影響與Mehlich-3浸提次數(shù)有關(guān),主要影響因子分別是鐵鋁錳氧化物和粘粒含量(單次浸提)或土壤pH值和有機(jī)質(zhì)(3次浸提).研究建立了基于可浸提態(tài)鎳對(duì)大麥根伸長的毒性閾值預(yù)測模型.研究結(jié)果證明了重金屬的化學(xué)浸提態(tài)代替全量作為風(fēng)險(xiǎn)評(píng)價(jià)標(biāo)準(zhǔn)的可行性,為提高土壤中鎳的生態(tài)風(fēng)險(xiǎn)評(píng)估的準(zhǔn)確性和科學(xué)性提供支撐.

    鎳;Mehlich-3;可浸提態(tài);土壤;生物有效性

    目前,大多數(shù)土壤生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)和土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)多以重金屬總量為指標(biāo)來判斷其生態(tài)毒性,基于土壤重金屬總量的生態(tài)風(fēng)險(xiǎn)評(píng)估和土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)在評(píng)估重金屬的生態(tài)毒害方面往往具有很大的不確定性[1],主要是由于總量不能很好地表示重金屬的移動(dòng)性和生物有效性[2].基于總量建立的生物毒害模型在生物有效性低的高背景土壤或不同污染源重金屬土壤中具有很大的局限性.重金屬的生態(tài)風(fēng)險(xiǎn)評(píng)估不僅應(yīng)該考慮重金屬含量還需要考慮其生物有效性.因此研究基于重金屬生物有效性或可浸提性的重金屬毒性及其預(yù)測模型對(duì)于土壤重金屬生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)是非常必要的[3].

    重金屬的生物有效部分被認(rèn)為是真正對(duì)有機(jī)體造成毒害的部分[4].重金屬的生物有效性可以反應(yīng)土壤重金屬污染程度以及對(duì)生態(tài)系統(tǒng)的潛在危害.明確重金屬的生物有效性是土壤重金屬生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)的關(guān)鍵.化學(xué)浸提法是目前應(yīng)用最廣泛的重金屬有效性的評(píng)價(jià)方法[5-6].其中Mehlich-3作為浸提劑[7],適用范圍廣泛,它不僅適用于酸性和中性土壤,也適用于堿性或石灰性土壤[8-10],是一種較好的聯(lián)合通用浸提劑.浸提劑中含有無機(jī)鹽組成的緩沖溶液和絡(luò)合劑,可以浸提土壤樣品中重金屬的可交換態(tài),部分碳酸鹽結(jié)合態(tài)和有機(jī)結(jié)合態(tài)[11].土壤浸提劑Mehlich-3在土壤營養(yǎng)元素的生物有效性評(píng)價(jià)已經(jīng)得到了廣泛應(yīng)用,但是在土壤重金屬毒性方面的研究還比較缺乏[9-10].

    本研究結(jié)合土壤理化性質(zhì),量化總外源添加鎳濃度和Mehlich-3可浸提態(tài)鎳濃度之間的關(guān)系;建立基于Mehlich-3可浸提態(tài)鎳對(duì)大麥根伸長的毒性閾值和土壤性質(zhì)的生物有效性模型.研究重金屬的化學(xué)浸提態(tài)代替全量作為風(fēng)險(xiǎn)評(píng)價(jià)標(biāo)準(zhǔn)的可行性,提高土壤中鎳的生態(tài)風(fēng)險(xiǎn)評(píng)估的準(zhǔn)確性和科學(xué)性.

    1 材料和方法

    1.1 土壤樣品的采集和處理

    土壤樣品采集于11個(gè)地點(diǎn)的表層(0~20cm)我國不同土壤性質(zhì)的,風(fēng)干備用,其中土壤樣品的pH值和有機(jī)碳含量分布與傳統(tǒng)農(nóng)業(yè)土壤基本一致,土壤主要理化性質(zhì)、土壤類型及種植農(nóng)作情況見表1.土壤pH值: 4.93~8.90,CEC: 7.47coml+/kg~ 33.6coml+/kg,OC: 0.6%~2.2%,黏粒含量(<2μm): 16%~66%.

    取過尼龍篩(<2mm)的風(fēng)干土,以NiCl2溶液的形式向土壤中均勻添加鎳.添加鎳的濃度根據(jù)土壤pH值不同設(shè)置7個(gè)濃度(不包括對(duì)照),土壤pH值<5時(shí),鎳的添加量分別是0,12.5,25,50,100,200,400, 800mg/kg;pH值為5~7時(shí),鎳的添加量分別是0,25,50,100,200,400,800,1600mg/kg;pH值>7時(shí),鎳的添加量分別是0,37.5,75,150,300,600,1200, 2400mg/ kg.所有處理保持100%最大持水量培養(yǎng)2d后,風(fēng)干,過篩(<2mm),保存.用ICP-OES測定土壤中鎳添加的準(zhǔn)確含量[19].將一半的土壤樣品采用模擬田間條件做淋洗處理減少重金屬添加時(shí)伴隨的鹽害同時(shí)縮小實(shí)驗(yàn)室處理和田間處理的差別.淋洗液的組份包括: CaCl25×10-4mol/L、Ca(NO3)25×10-4mol/L、MgCl25×10-4mol/L、Na2SO410-4mol/L和KCl 10-4mol/L,pH=5.9[20].淋洗處理后的樣品風(fēng)干、過篩(<2mm),用ICP-OES測定淋洗處理后土壤樣品中鎳含量[19].

    表1 土壤的主要理化性質(zhì)

    注:a去離子水測定(水土比為5:1)[12];b氯化銨方法測定[12];c有機(jī)碳含量測定方法[13-14];d檸檬酸鹽可提取態(tài)鋁、鐵和錳測定方法[15-16];e草酸可提取態(tài)鋁和鐵測定方法[17-18].

    1.2 化學(xué)浸提方法

    1.2.1 Mehlich-3浸提劑的配制 NH4F-EDTA貯備液:稱取27.78g NH4F和14.61g EDTA,加入去離子水溶解,定容至120mL去離子水中,加入14.61g EDTA去離子水定容至200mL,貯于塑料瓶中保存(在冰箱中冷藏可長期使用);Mehlich-3試劑:稱取20.0g NH4NO3溶于約500mL去離子水中,加入4mL上述NH4F-EDTA貯備液,再加入11.5mL CH3COOH和0.82mL HNO3,去離子水定容至1L.pH值應(yīng)為(2.5±0.1),貯于塑料瓶中備用.

    1.2.2 土壤中鎳的浸提 土壤樣品稱取5.00g左右于50mL離心管中,加入25mL Mehlich-3溶液,擰緊離心管瓶蓋,室溫下(25℃)于振蕩機(jī)上200r/min振蕩2h.振蕩結(jié)束后,離心管放置在離心機(jī)上4000r/min離心20min,取上清液經(jīng)孔徑0.45μm的濾膜或?yàn)V紙過濾后,收集濾液于塑料小方瓶中,待測,1次Mehlich-3浸提完成.離心后的土壤沉積物中再次加入25mL Mehlich-3溶液,按照上述方法再次振蕩,離心,過濾,收集濾液.連續(xù)浸提3次,3次Mehlich-3浸提完成.采用原子吸收光譜法(ZEEnit 700;Analytik Jena Germany)測定濾液中Mehlich-3浸提態(tài)鎳的濃度.所有處理做2個(gè)重復(fù).總樣本數(shù)為1056個(gè).

    1.3 基于Mehlich-3浸提態(tài)鎳含量的植物毒性計(jì)算

    1.3.1 大麥很伸長試驗(yàn) 大麥根伸長試驗(yàn)根據(jù)國際標(biāo)準(zhǔn)方法(ISO)11269-1[21]進(jìn)行.將大麥在氣候培養(yǎng)箱中進(jìn)行培養(yǎng),生長條件為白天14h(22℃),夜間10h(18℃),光照強(qiáng)度為24000lux/m2,濕度保持在70%,保持整個(gè)試驗(yàn)過程中土壤含水量為最大持水量的65%,測量大麥5d的根伸長,相關(guān)試驗(yàn)細(xì)節(jié)和大麥根伸長數(shù)據(jù)參考Li等[19].

    1.3.2 劑量-效應(yīng)曲線擬合 根據(jù)劑量-效應(yīng)log- logistic[22]方程擬合得到抑制大麥根伸長10%、20%和50%的Mehlich-3浸提態(tài)鎳的濃度.

    式中:為大麥相對(duì)根伸長;為Mehlich-3浸提態(tài)鎳的濃度,mg/kg;為EC(EC10、EC20、EC50)取對(duì)數(shù);0和是模型擬合參數(shù).

    1.3.3 低劑量毒物刺激作用曲線擬合 根據(jù)方程得到毒性刺激作用曲線和毒性閾值EC10、EC20和EC50以及95%置信區(qū)間.

    式中:為大麥相對(duì)根伸長;為Mehlich-3浸提態(tài)鎳的濃度,mg/kg;、、和是曲線擬合參數(shù),是與EC相關(guān)的變量,當(dāng)分別等于10、20或50時(shí),參數(shù)分別定義為EC10、C20和EC50.關(guān)于低劑量毒物刺激作用的詳細(xì)描述可參考文獻(xiàn)[19,24].

    1.4 數(shù)據(jù)處理

    采用Excel 2010做數(shù)據(jù)規(guī)劃求解處理,結(jié)合土壤性質(zhì),建立土壤中Mehlich-3浸提態(tài)鎳濃度和總外源添加鎳濃度之間的關(guān)系.

    采用SPSS 19.0進(jìn)行逐步回歸分析,建立基于Mehlich-3浸提態(tài)鎳對(duì)大麥根伸長的毒性閾值和土壤性質(zhì)之間的關(guān)系.

    2 結(jié)果與討論

    2.1 土壤中可浸提態(tài)鎳與外源添加鎳的關(guān)系

    采用Mehlich-3可浸提態(tài)鎳占外源添加鎳的百分?jǐn)?shù)來表示Mehlich-3對(duì)土壤中鎳的浸提率.對(duì)淋洗和非淋洗土壤而言,1次Mehlich-3浸提率平均分別為36%和40%;3次Mehlich-3浸提率平均分別為64%和75%.浸提率隨著pH值的升高而降低,在pH<7的土壤中浸提率平均最高為92%,而在pH>8的土壤中浸提率平均最高為67%.

    為了進(jìn)一步量化總外源添加鎳濃度和Mehlich-3可浸提態(tài)鎳濃度之間的關(guān)系,以及土壤性質(zhì)的影響,本研究結(jié)合土壤理化性質(zhì),建立了總外源添加鎳濃度和Mehlich-3可浸提態(tài)鎳濃度之間的線性回歸模型,其中,總外源添加鎳濃度分為0~300mg/ kg和0~2000mg/kg 2個(gè)濃度梯度(表2).1次Mehlich-3浸提量每提高1個(gè)單位,外源添加總量在非淋洗和淋洗土壤中分別提高了2.238和3.355個(gè)單位;3次Mehlich-3浸提量每提高1個(gè)單位,外源添加總量在非淋洗和淋洗土壤中分別提高了1.596和1.967個(gè)單位.當(dāng)總外源添加鎳濃度在0~300mg/kg時(shí),基于1次Mehlich-3可浸提態(tài)鎳濃度與外源添加鎳濃度的線性回歸方程的決定系數(shù)(2)在非淋洗和淋洗土壤中可以達(dá)到0.65和0.80,而基于3次Mehlich-3可浸提態(tài)鎳濃度與外源添加鎳濃度的線性回歸方程的2可以分別提高到0.86和0.90.當(dāng)總外源添加鎳的濃度在0~2000mg/kg時(shí),在非淋洗和淋洗土壤中,1次Mehlich-3可浸提態(tài)鎳濃度可以解釋線性回歸模型中總外源添加鎳的濃度變異的68%和76%,3次Mehlich-3可浸提態(tài)鎳濃度可以解釋線性回歸模型中總外源添加鎳濃度變異的92%和93%.在擬合的線性回歸方程中加入pH值作為影響因子,線性回歸方程的2有較小提高,而其他的土壤理化性質(zhì),如陽離子交換量(CEC)和有機(jī)質(zhì)(OC)等,則沒有起到提高線性回歸模型擬合效果的作用.基于1次或3次Mehlich-3可浸提態(tài)鎳和pH值建立線性回歸模型預(yù)測土壤中總外源添加鎳,結(jié)果見圖1.結(jié)果顯示在非淋洗和淋洗土壤中,基于3次Mehlich-3浸提建立的線性回歸模型2更高,模型擬合程度更好.

    表2 土壤總外源添加鎳的濃度和Mehlich-3浸提態(tài)鎳的濃度的線性回歸關(guān)系

    注: Total-Ni: 總外源添加鎳的濃度(mg/kg),M3-Ni: Mehlich-3浸提態(tài)鎳的濃度(mg/kg),r: 決定系數(shù).

    圖1 土壤中總外源添加鎳濃度的預(yù)測值和測量值

    (A)由表2的方程2和方程4計(jì)算得出;(B)由表2的方程6和方程8計(jì)算得出;(C)由表2的方程10和方程12計(jì)算得出;(D)由表2的方程14和方程16計(jì)算得出.UL和L分別代表非淋洗和淋洗土壤

    在大多數(shù)文獻(xiàn)中,土壤中重金屬的總量和植物體中的含量沒有線性關(guān)系[25],也就是說,不是所有的重金屬形態(tài)都能被植物吸收利用,植物對(duì)重金屬的吸收與重金屬的有效態(tài)部分顯著相關(guān)[26].Zhang等[2]研究發(fā)現(xiàn),植物可吸收利用的重金屬主要來自土壤中的水溶態(tài)和可交換態(tài).傳統(tǒng)的以總量作為標(biāo)準(zhǔn)判斷重金屬毒性是存在缺陷的,難以評(píng)估背景值高但有效性低或者不同來源重金屬的毒性,會(huì)出現(xiàn)高估或者低估重金屬毒性的情況[27].因此,在土壤重金屬生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)中應(yīng)該更多的關(guān)注重金屬的生物有效性.

    結(jié)合土壤性質(zhì)的,總外源添加鎳濃度和Mehlich-3可浸提態(tài)鎳濃度之間存在顯著的數(shù)量關(guān)系為土壤中鎳的可浸提態(tài)替代總外源添加鎳濃度作為重金屬風(fēng)險(xiǎn)評(píng)價(jià)的依據(jù)提供理論基礎(chǔ).為了保證線性回歸模型的穩(wěn)定性和準(zhǔn)確性,基于3次Mehlich-3可浸提態(tài)鎳濃度建立線性回歸模型更加合理.

    2.2 基于可浸提態(tài)鎳對(duì)大麥的毒性

    表3 基于1次Mehlich-3浸提態(tài)的鎳對(duì)大麥根伸長的毒性閾值

    注: EC(=10、20、50): 抑制大麥根伸長10%、20%和50%的Mehlich-3浸提態(tài)鎳的濃度;CI: ±95%的置信區(qū)間;—: CI無法確定;NC: 毒性閾值無法確定;a低劑量毒物刺激作用;*淋洗前后的毒性閾值有顯著性差異,<0.05.

    表4 基于3次Mehlich-3浸提態(tài)的鎳對(duì)大麥根伸長的毒性閾值

    注: EC(=10、20、50): 抑制大麥根伸長10%、20%和50%的Mehlich-3浸提態(tài)鎳的濃度;CI: ±95%的置信區(qū)間;—: CI無法確定;NC: 毒性閾值無法確定;a低劑量毒物刺激作用;*淋洗前后的毒性閾值有顯著性差異,<0.05.

    2.2.1 毒性閾值 1次和3次Mehlich-3可浸提態(tài)鎳濃度分別作為劑量,計(jì)算大麥根伸長作為效應(yīng)的毒性閾值及其95%的置信區(qū)間(表3和4).在非淋洗土壤中,基于1次Mehlich-3可浸提態(tài)鎳濃度的EC10、EC20、EC50的變化范圍分別是1~441mg/ kg,2~459mg/kg,9~490mg/kg,毒性閾值變化及其范圍倍數(shù)如下: EC10(441倍)>EC20(230倍)>EC50(54倍);在淋洗土壤中,基于1次Mehlich-3可浸提態(tài)鎳濃度的EC10從0.42mg/kg增加到384mg/kg,EC20從1mg/ kg增加到576mg/kg,EC50從4mg/kg增加到505mg/ kg,毒性閾值變化及其范圍倍數(shù)EC20(576倍)>EC10(384倍)> EC50(126倍).在非淋洗的重慶土壤和淋洗后的廣州土壤中出現(xiàn)低劑量毒物刺激作用.淋洗土壤中,鄭州土壤的EC20,烏魯木齊、楊凌、鄭州和德州土壤的EC50無法計(jì)算得到,為淋洗降低鎳毒性,盡管1次Mehlich-3可浸提態(tài)鎳濃度為16~409mg/kg.淋洗降低鎳毒性也出現(xiàn)在基于總外源添加鎳濃度的毒性閾值[19],后面將進(jìn)一步討論.基于3次Mehlich-3可浸提態(tài)鎳濃度的EC10、EC20、EC50相比基于1次Mehlich-3可浸提態(tài)鎳濃度的EC10、EC20、EC50在非淋洗土壤中分別增加101%、114%、147%;在淋洗土壤中分別增加117%、158%、87%,同樣存在閾值變化范圍大的情況.在Li等[19]以總量表示鎳毒性時(shí)出現(xiàn)相似的閾值變化范圍大的情況,其中pH值作為單一影響因子可以解釋EC50變異的70%以上.Rooney等[29]基于歐洲土壤(pH 3.3~8.7)建立鎳毒性模型時(shí)發(fā)現(xiàn)土壤pH值、有效陽離子交換量、粘粒含量和有機(jī)碳含量是解釋毒性閾值變化范圍大的主要因子.結(jié)果說明,無論是以重金屬的總量還是浸提態(tài)表示毒性都不能忽略土壤性質(zhì)的影響.

    表5 基于1次Mehlich-3浸提的大麥根伸長的毒性閾值與土壤性質(zhì)之間的線性回歸關(guān)系

    注:(1),(2),(3),(4),(5)分別代表方程中從左到右的影響因子的顯著性水平;2:決定系數(shù);adj2:調(diào)整決定系數(shù).

    表6 基于3次Mehlich-3浸提的大麥根伸長的毒性閾值與土壤性質(zhì)之間的線性回歸關(guān)系

    注:(1),(2),(3)分別代表方程中從左到右的影響因子的顯著性水平;2:決定系數(shù);adj2:調(diào)整決定系數(shù).

    2.2 淋洗對(duì)鎳毒性的影響 本研究中,基于1次或者3次Mehlich-3可浸提態(tài)鎳濃度表示鎳毒性時(shí),淋洗處理降低鎳毒性(表3和表4),9%的EC10值、36%的EC20值和55%的EC50值在淋洗處理后顯著提高.

    淋洗處理對(duì)鎳毒性的影響與土壤理化性質(zhì)有關(guān),淋洗處理降低鎳毒性多出現(xiàn)在堿性土壤中,在酸性和中性土壤中并不顯著.淋洗作用在堿性或石灰性土壤,和在酸性或中性土壤中不同的原因可能是可溶性鎳的釋放過程不同.Ponizovsky等[29]研究發(fā)現(xiàn),在非石灰性土壤中鎳的釋放可以用吸附方程表示,在石灰性土壤中鎳的釋放可能受土壤表面沉淀的影響.當(dāng)鎳毒性以總量濃度表達(dá)時(shí),Oorts等[30]發(fā)現(xiàn),相比在酸性和中性土壤中(pH 4.6~6.1),淋洗處理顯著降低了堿性土壤中鎳對(duì)微生物潛在硝化率和玉米殘?jiān)V化的毒害.Li等[19]在pH38.2的土壤中發(fā)現(xiàn)淋洗處理降低鎳對(duì)大麥根伸長的毒害.這些結(jié)論與本研究的結(jié)果相似.金屬的分配系數(shù)是指恒溫下固相溶質(zhì)濃度與液相溶質(zhì)濃度達(dá)到平衡時(shí)的比值[31],與土壤pH值密切相關(guān),隨著pH值的升高而增大[32].也就是說,在高pH值的土壤中,淋洗處理很大程度上降低了土壤溶液中鎳的含量.淋洗處理降低鎳毒性的原因可能還有淋洗降低了土壤的pH值和鹽度,增加金屬鹽含量導(dǎo)致鎳在固液界面重新分配[19].同時(shí),淋洗降低重金屬毒性與土壤性質(zhì),重金屬種類和生態(tài)毒性終點(diǎn)都有關(guān)系[33].

    2.3 基于可浸提態(tài)鎳的植物毒性模型

    基于1次Mehlich-3可浸提態(tài)鎳對(duì)大麥根伸長的毒性閾值和土壤性質(zhì)的簡單和多元線性回歸模型(表5)中Fe(Al和Mn)氧化物含量、黏粒含量(clay)是植物毒性模型的主要參數(shù).在非淋洗土壤中,當(dāng)檸檬酸鹽提取態(tài)鐵(CDFe)作為單一影響因子能夠解釋EC10變異的70.4%,EC20變異的75.3%和EC50變異的82.0%.當(dāng)線性回歸模型中加入其他影響因子時(shí)決定系數(shù)(2)有所提高.擬合效果最好的是CDFe、草酸提取態(tài)鐵(OXFe)、檸檬酸鹽提取態(tài)鋁(CDAl)和clay作為影響因子構(gòu)建的多元線性回歸模型,能夠解釋EC50變異的98.9%.在淋洗土壤中,以CDFe作為影響因子的線性回歸模型分別解釋了EC10變異的77.3%、EC20變異的86.2%和EC50變異的87.3%.

    基于3次Mehlich-3可浸提態(tài)鎳對(duì)大麥根伸長的毒性閾值和土壤性質(zhì)的簡單和多元線性回歸模型(表6)中pH值和OC是主要參數(shù).在非淋洗土壤中,pH值作為單一影響因子能夠解釋EC20變異的38.4%,EC50變異的64.0%.當(dāng)線性回歸模型中加入OC因子時(shí)2分別提高到0.749(EC20)和0.838(EC50).在淋洗土壤中,以pH值作為影響因子的線性回歸模型分別解釋了EC10變異的52.7%、EC20變異的82.4%和EC50變異的66.4%.

    大量的研究表明,pH值是影響鎳毒性和鎳在土壤中分配最主要的因子[33-36].Li等[19]基于土壤中總外源添加鎳對(duì)大麥根伸長的毒性閾值和土壤性質(zhì)構(gòu)建經(jīng)驗(yàn)預(yù)測模型時(shí)發(fā)現(xiàn),土壤pH值是預(yù)測鎳毒性最重要的影響因子.在銅鎳對(duì)西紅柿和小白菜生物量毒性的研究中發(fā)現(xiàn),pH值對(duì)鎳毒性的影響要比對(duì)銅的強(qiáng)烈[33].pH值影響鎳在固液界面的分配[31],鎳在土壤中的老化過程也受pH值影響[37].無論是從化學(xué)行為,還是鎳毒性來說,pH值都是影響鎳毒性最重要的因子.而基于1次Mehlich-3可浸提態(tài)鎳對(duì)大麥根伸長的毒性閾值和土壤性質(zhì)建立生物有效性模型時(shí),pH值不是最主要的影響因子,可能是基于1次Mehlich-3浸提態(tài)鎳建立有效性模型時(shí),相比pH來說,土壤組分的影響更大.

    3 結(jié)論

    3.1 土壤pH值是影響土壤Mehlich-3可浸提鎳與外源添加鎳二者關(guān)系最主要的影響因子;結(jié)合土壤pH值,其二者存在顯著的數(shù)量關(guān)系.

    3.2 土壤性質(zhì)對(duì)Mehlich-3浸提鎳的生物毒性的影響與Mehlich-3浸提次數(shù)有關(guān).1次Mehlich-3浸提態(tài)鎳時(shí),檸檬酸鹽提取態(tài)鐵、鋁和錳以及粘粒含量是預(yù)測鎳毒性的主要影響因子;3次Mehlich-3浸提時(shí),土壤pH值和有機(jī)質(zhì)是預(yù)測鎳毒性的主要影響因子.

    3.3 為保證預(yù)測模型的穩(wěn)定性和準(zhǔn)確性,3次Mehlich-3浸提比1次Mehlich-3浸提更為合理.

    [1] Janssen C R, Schamphelaere K D, Heijerick D, et al. Uncertainties in the Environmental Risk Assessment of Metals [J]. Human and Ecological Risk Assessment An International Journal, 2000,6(6):1003-1018.

    [2] Zhang M K, Liu Z Y, Wang H. Use of Single Extraction Methods to Predict Bioavailability of Heavy Metals in Polluted Soils to Rice [J]. Communications in Soil Science and Plant Analysis, 2010,41(7):820-831.

    [3] Weng L P, Wolthoorn A, Lexmond T M, et al. Understanding the effects of soil characteristics on phytotoxicity and bioavailability of nickel using speciation models [J]. Environmental Science and Technology, 2004,38(1):156-162.

    [4] Peijnenburg W J, Jager T. Monitoring approaches to assess bioaccessibility and bioavailability of metals: matrix issues [J]. Ecotoxicology and Environmental Safety, 2003,56(1):63-77.

    [5] Zhou Y, Ning X A, Liao X, et al. Characterization and environmental risk assessment of heavy metals found in fly ashes from waste filter bags obtained from a Chinese steel plant [J]. Ecotoxicology and Environmental Safety, 2013,95(1):130-136.

    [6] Nannoni F, Protano G. Chemical and biological methods to evaluate the availability of heavy metals in soils of the Siena urban area (Italy) [J]. Science of the Total Environment, 2016,568(15):1-10.

    [7] Mehlich A. Mehlich 3soil test extractant: A modification of Mehlich 2extractant [J]. Communications in Soil Science and Plant Analysis, 1984,15(12):1409-1416.

    [8] 劉 肅,李酉開. Mehlich 3通用浸提劑的研究[J]. 土壤學(xué)報(bào), 1995,32(2):132-141.

    [9] Sims J T. Comparison of mehlich 1and mehlich 3extractants for P, K, Ca, Mg, Mn, Cu and Zn in atlantic coastal plain soils1 [J]. Communications in Soil Science and Plant Analysis, 1989,20 (17/18):1707-1726.

    [10] Monterroso C, Alvarez E, Fernández Marcos M L. Evaluation of Mehlich 3reagent as a multielement extractant in mine soils [J]. Land Degradation and Development, 2015,10(1):35-47.

    [11] 楊 潔,瞿 攀,王金生,等.土壤中重金屬的生物有效性分析方法及其影響因素綜述[J]. 環(huán)境污染與防治, 2017,39(2):217-223.

    [12] Rayment G E, Higginson F R. Australian laboratory handbook of soil and water chemical methods [M]. Melbourne, Australia: Inkata Press Pty Ltd, 1992:5-7.

    [13] Igor M. Determination of carbon and nitrogen in samples of various soils by the dry combustion [J]. Communications in Soil Science and Plant Analysis, 1997,28(17/18):1499-1511.

    [14] Sherrod L A, Dunn G, Peterson G A, et al. Inorganic Carbon Analysis by Modified Pressure-Calcimeter Method [J]. Soil Science Society of America Journal, 2002,66(1):299-305.

    [15] McKeague J A, Day J H. Differenzierung- and Oxalate-extractable Fe and Al as aids in differentiating various classes of soils [J]. Canadian Journal of Soil Science, 1966,46(1):13-22.

    [16] Schwertmann U. Dissolution Kinetics of Nickel Surface Precipitates on Clay Mineral and Oxide Surfaces [J]. Soil Science Society of America Journal, 2010,65(3):685-694.

    [17] Mehra O P, Jackson M L. Iron oxide removal from soils and clays by a dithionite-citrate system buffered with sodium bicarbonate [J]. Clays and Clay Minerals, 1960,7(1):317-327.

    [18] Jackson M L, Lim C H, Zelazny L W, et al. Oxides, hydroxides, and aluminosilicates. In Methods of Soil Analysis. Part 1. Physical and Mineralogical Methods [M]. USA: American Society of Agronomy, 1986:101-150.

    [19] Li B, Zhang H, Ma Y B, et al. Influences of soil properties and leaching on nickel toxicity to barley root elongation [J]. Ecotoxicology and environmental safety, 2011,74(3):459-466.

    [20] Oorts K, Ghesquiere U E. Leaching and aging decrease nickel toxicity to soil microbial processes in soils freshly spiked with nickel chloride [J]. Environmental Toxicology and Chemistry, 2010,26(6):1130-1138.

    [21] Soil-Quality - Determination of Effects of Pollutants on Soil Flora - Part 1: Method for the Measuredment of Inhibition of Root Growth [S]. International Organization for Standardization (ISO): Geneva, Switzerland, 1993.

    [22] Haanstra L, Doelman P, Voshaar J H O. The use of sigmoidal dose response curve in soil ecotoxicological research [J]. Plant Soil., 1985, 84(2):293-297.

    [23] Schabenberger O, Tharp B E, Kells J J, et al. Statistical tests for hormesis and effective dosages in herbicide dose response [J]. Agronomy Journal, 1999,91(1):713-721.

    [24] 郭雪雁,馬義兵,李 波.陸地生態(tài)系統(tǒng)中低劑量毒物刺激作用及擬合模型研究進(jìn)展[J]. 生態(tài)學(xué)報(bào), 2009,29(8):4408-4419.

    [25] Chojnacka K, Chojnacki A, Górecka H, et al. Bioavailability of heavy metals from polluted soils to plants [J]. Science of the Total Environment, 2005,337(1/3):175-182.

    [26] Xiao L, Guan D, Peart M R, et al. The influence of bioavailable heavy metals and microbial parameters of soil on the metal accumulation in rice grain [J]. Chemosphere, 2017,185:868-878.

    [27] Wang X Q, Wei D P, Ma Y B, et al. Soil ecological criteria for nickel as a function of soil properties [J]. Environmental Science and Pollution Research, 2018,25(3):2137-2146.

    [28] Rooney C P, Zhao F J, McGrath S P. Phytotoxicity of nickel in a range of European soils: Influence of soil properties, Ni solubility and speciation [J]. Environmental Pollution, 2007,145(2):596-605.

    [29] Ackerman A J. Effect of field aging on nickel concentration in soil solutions [J]. Communications in Soil Science and Plant Analysis, 2008,39(3/4):510-523.

    [30] Oorts K, Ghesquiere U E. Leaching and aging decrease nickel toxicity to soil microbial processes in soils freshly spiked with nickel chloride [J]. Environmental Toxicology and Chemistry, 2007,26(6):1130-1138.

    [31] Sébastien Sauvé, William Hendershot, Allen H E. Solid-Solution Partitioning of Metals in Contaminated Soils: Dependence on pH, Total Metal Burden, and Organic Matter [J]. Environmental Science and Technology, 2000,34(7):1125-1131.

    [32] Tipping E, Rieuwerts J, Pan G, et al. The solid-solution partitioning of heavy metals (Cu, Zn, Cd, Pb) in upland soils of England and Wales [J]. Environmental Pollution, 2003,125(2):213-225.

    [33] Li B, Zhang H T, Ma Y B, et al. Relationships between soil properties and toxicity of copper and nickel to bok choy and tomato in Chinese soils [J]. Environmental Toxicology and Chemistry, 2013,32(10):2372-2378.

    [34] Weng L, Lexmond T M, Wolthoorn A, et al. Phytotoxicity and bioavailability of nickel: chemical speciation and bioaccumulation [J]. Environmental Toxicology and Chemistry, 2003,22(9):2180-2187.

    [35] Echevarria G, Massoura S T, Sterckeman T, et al. Assessment and control of the bioavailability of nickel in soils [J]. Environmental Toxicology and Chemistry, 2006,25(6):643-651.

    [36] Everhart J L, Jr M N D, Peltier E, et al. Assessing nickel bioavailability in smelter-contaminated soils [J]. Science of the Total Environment, 2006,367(2/3):732-744.

    [37] Ma Y B, Lombi E, McLaughlin M J, et al. Aging of nickel added to soils as predicted by soil pH and time [J]. Chemosphere, 2013,92(8):962-968.

    Toxicity thresholds based on Mehlich-3 extractable nickel to barley root elongation.

    ZHU Guang-yun1,2, JIANG Bao2, LI Ju-mei2, LI He-lian1, MA Yi-bing3*

    (1.School of Resources and Environment, University of Jinan, Jinan 250022, China;2.Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing 100081, China;3.Guangdong Institute of Eco-environmental Science & Technology, Guangzhou 510650, China)., 2018,38(8):3143~3150

    After water-soluble nickel salt (NiCl2) was spiked to 11representative Chinese soil samples, the spiked soil samples were incubated and treated with or without artificial rainwater leaching in order to simulate field conditions. The toxicity of these soil samples to barley root elongation was studied by a bioassay method, and the extractability of the added Ni in the soils was evaluated using single or sequential three Mehlich-3extractions. The results showed that the extractability of Ni added to soil was affected by soil pH, and Mehlich-3extractable Ni is as a function of total added Ni in soils and soil pH significantly. The phytotoxicity thresholds based on Mehlich-3extractable Ni were significantly affected by soil properties. In other words, the effective concentration causing 50% inhibition based on Mehlich-3extractable Ni varied with338-fold difference among soils. Regression analysis indicated that the soil citrate dithionate extractable Fe, Al and Mn, and clay content in soils were the most important factors to affect the phytotoxicity based on single Mehlich-3extractable Ni, and that soil pH and organic content were the most important factors to affect the phytotoxicity based on sequential three Mehlich-3extractable Ni. The influence of soil properties on the phytotoxicity based on Mehlich-3extractable Ni was associated with the number of extraction times. Finally, the predictive models for phytotoxicity thresholds based on Mehlich-3extractable Ni to barley root elongation were developed. The results indicated that the extractable Ni in soil can be used as criteria in order to improve the accuracy and scientificity of Ni ecological risk assessment in soil.

    nickel;Mehlich-3;extractable;soil;bioavailability

    X171.5

    A

    1000-6923(2018)08-3143-08

    朱廣云(1992-),女,山東臨沂人,濟(jì)南大學(xué)碩士研究生,主要從事土壤修復(fù)技術(shù)方面的研究.發(fā)表論文3篇.

    2017-12-19

    廣東省科學(xué)院引進(jìn)高層次領(lǐng)軍人才專項(xiàng)資金項(xiàng)目(2016GDASRC-0210)和國家重點(diǎn)研發(fā)計(jì)劃《農(nóng)田系統(tǒng)重金屬遷移轉(zhuǎn)化和安全閾值研究》項(xiàng)目(2016YFD0800400)

    * 責(zé)任作者, 研究員, mayibing@caas.ac.cn

    猜你喜歡
    大麥外源毒性
    具有外源輸入的船舶橫搖運(yùn)動(dòng)NARX神經(jīng)網(wǎng)絡(luò)預(yù)測
    我的大麥哥哥
    大麥蟲對(duì)聚苯乙烯塑料的生物降解和礦化作用
    Global interest in Chinese baijiu
    外源鉛脅迫對(duì)青稞生長及鉛積累的影響
    動(dòng)物之最——毒性誰最強(qiáng)
    外源鈣對(duì)干旱脅迫下火棘種子萌發(fā)的影響
    外源添加皂苷對(duì)斑玉蕈生長發(fā)育的影響
    大麥若葉青汁
    RGD肽段連接的近紅外量子點(diǎn)對(duì)小鼠的毒性作用
    99热精品在线国产| av中文乱码字幕在线| h日本视频在线播放| 97热精品久久久久久| 在线观看一区二区三区| 淫妇啪啪啪对白视频| 久9热在线精品视频| 日本爱情动作片www.在线观看 | 国产欧美日韩精品一区二区| 国产精华一区二区三区| 国产一区二区激情短视频| 成人毛片a级毛片在线播放| 久久久久久久亚洲中文字幕| 51国产日韩欧美| 欧美成人性av电影在线观看| 日韩欧美国产一区二区入口| 亚洲性久久影院| 亚洲av电影不卡..在线观看| 欧美性猛交╳xxx乱大交人| 久久久午夜欧美精品| 麻豆一二三区av精品| 国产精品美女特级片免费视频播放器| 久久久久久大精品| 亚洲专区中文字幕在线| 老司机深夜福利视频在线观看| x7x7x7水蜜桃| 91久久精品国产一区二区三区| 国内精品宾馆在线| 国模一区二区三区四区视频| 偷拍熟女少妇极品色| 亚洲第一区二区三区不卡| 日本在线视频免费播放| 日本 av在线| 搡老岳熟女国产| 国产精品一区二区三区四区免费观看 | 成人永久免费在线观看视频| 夜夜夜夜夜久久久久| 亚洲黑人精品在线| 精品乱码久久久久久99久播| ponron亚洲| 国产午夜精品论理片| 蜜桃亚洲精品一区二区三区| av在线天堂中文字幕| 亚洲精品成人久久久久久| 人妻丰满熟妇av一区二区三区| 给我免费播放毛片高清在线观看| 亚洲自偷自拍三级| 精品人妻偷拍中文字幕| 亚洲精华国产精华液的使用体验 | 亚洲av中文字字幕乱码综合| 夜夜夜夜夜久久久久| 亚洲aⅴ乱码一区二区在线播放| 欧美bdsm另类| 熟女电影av网| 国产av麻豆久久久久久久| 国产成人av教育| 听说在线观看完整版免费高清| av中文乱码字幕在线| 日日摸夜夜添夜夜添av毛片 | 窝窝影院91人妻| 日本一二三区视频观看| 国产精品一及| 午夜福利在线观看免费完整高清在 | 亚洲av成人精品一区久久| 午夜免费激情av| 日本 av在线| 老熟妇仑乱视频hdxx| 不卡视频在线观看欧美| 亚洲图色成人| 韩国av在线不卡| 精品一区二区三区视频在线| 天美传媒精品一区二区| 综合色av麻豆| 九色成人免费人妻av| 黄色女人牲交| 午夜福利18| 国内少妇人妻偷人精品xxx网站| 看十八女毛片水多多多| 国产激情偷乱视频一区二区| 老师上课跳d突然被开到最大视频| 国产视频内射| 熟女人妻精品中文字幕| 日韩亚洲欧美综合| 在线观看午夜福利视频| 一个人免费在线观看电影| 琪琪午夜伦伦电影理论片6080| 国产蜜桃级精品一区二区三区| 精品久久久久久,| 一进一出抽搐gif免费好疼| 熟女电影av网| av福利片在线观看| 亚洲色图av天堂| 国产精品国产高清国产av| 亚洲国产日韩欧美精品在线观看| 欧美高清性xxxxhd video| 精品无人区乱码1区二区| 最近最新免费中文字幕在线| 国产一区二区三区在线臀色熟女| 他把我摸到了高潮在线观看| 日本免费一区二区三区高清不卡| 欧美人与善性xxx| 精品久久久噜噜| 精品午夜福利视频在线观看一区| 国产成人aa在线观看| 欧美成人a在线观看| 我要搜黄色片| 亚洲狠狠婷婷综合久久图片| 国产精品av视频在线免费观看| 国产三级中文精品| 全区人妻精品视频| 国产在视频线在精品| 99riav亚洲国产免费| 日韩欧美国产一区二区入口| 国产成年人精品一区二区| 欧美一区二区精品小视频在线| av在线蜜桃| 麻豆久久精品国产亚洲av| 搡老熟女国产l中国老女人| 午夜福利成人在线免费观看| 国产精品嫩草影院av在线观看 | 亚洲五月天丁香| 91久久精品国产一区二区三区| 日韩欧美一区二区三区在线观看| 国产精品福利在线免费观看| 精品一区二区三区视频在线观看免费| 国产激情偷乱视频一区二区| 别揉我奶头 嗯啊视频| 日本撒尿小便嘘嘘汇集6| 国产亚洲av嫩草精品影院| 一级毛片久久久久久久久女| 亚洲熟妇中文字幕五十中出| 亚洲精品在线观看二区| 免费av不卡在线播放| 麻豆一二三区av精品| 校园春色视频在线观看| 国产成人av教育| 国模一区二区三区四区视频| 99九九线精品视频在线观看视频| 高清日韩中文字幕在线| 色在线成人网| 免费观看在线日韩| netflix在线观看网站| 一卡2卡三卡四卡精品乱码亚洲| 日本黄色片子视频| 99在线视频只有这里精品首页| 亚洲av日韩精品久久久久久密| av中文乱码字幕在线| 免费av不卡在线播放| 18+在线观看网站| 真实男女啪啪啪动态图| 久久久久国内视频| av在线天堂中文字幕| 日本免费一区二区三区高清不卡| 国产老妇女一区| 久久久久国产精品人妻aⅴ院| 两性午夜刺激爽爽歪歪视频在线观看| 成人一区二区视频在线观看| 九九爱精品视频在线观看| 欧美高清性xxxxhd video| aaaaa片日本免费| 日本欧美国产在线视频| 日本 欧美在线| 欧美一区二区亚洲| 精品久久国产蜜桃| 一进一出好大好爽视频| 少妇被粗大猛烈的视频| 亚洲美女黄片视频| 午夜福利在线在线| 亚洲精品一卡2卡三卡4卡5卡| 国产精品99久久久久久久久| 日韩欧美国产在线观看| 免费看av在线观看网站| 夜夜爽天天搞| 国产精品自产拍在线观看55亚洲| 国产成年人精品一区二区| 色噜噜av男人的天堂激情| 日韩精品青青久久久久久| 国产精品一区二区三区四区免费观看 | 久久99热这里只有精品18| 在线观看av片永久免费下载| 婷婷精品国产亚洲av| 小蜜桃在线观看免费完整版高清| 干丝袜人妻中文字幕| xxxwww97欧美| 一级黄片播放器| 麻豆国产av国片精品| 黄色女人牲交| 国产av在哪里看| 国产爱豆传媒在线观看| 亚洲成人久久性| 97碰自拍视频| 国产av在哪里看| 国产 一区精品| 一a级毛片在线观看| 久久精品国产自在天天线| 婷婷色综合大香蕉| 一进一出好大好爽视频| 99热只有精品国产| 99热网站在线观看| 国产成年人精品一区二区| 亚洲内射少妇av| 亚洲av日韩精品久久久久久密| a级毛片a级免费在线| 啦啦啦啦在线视频资源| 麻豆精品久久久久久蜜桃| 淫妇啪啪啪对白视频| 国产毛片a区久久久久| h日本视频在线播放| 伊人久久精品亚洲午夜| 国产精品一区二区三区四区免费观看 | 亚洲最大成人av| 听说在线观看完整版免费高清| 超碰av人人做人人爽久久| 深爱激情五月婷婷| 国产精品女同一区二区软件 | 黄色欧美视频在线观看| 男女那种视频在线观看| 亚洲精品久久国产高清桃花| 看十八女毛片水多多多| 亚洲欧美日韩无卡精品| 3wmmmm亚洲av在线观看| 久久香蕉精品热| 三级男女做爰猛烈吃奶摸视频| 国产精品1区2区在线观看.| 亚洲专区国产一区二区| 国产视频一区二区在线看| 成人特级黄色片久久久久久久| 老熟妇仑乱视频hdxx| 精品人妻1区二区| 日本色播在线视频| 久久精品国产清高在天天线| 国产欧美日韩精品亚洲av| 中文字幕熟女人妻在线| 国产一区二区三区在线臀色熟女| 亚洲av一区综合| 免费人成视频x8x8入口观看| 成人无遮挡网站| 禁无遮挡网站| 99精品久久久久人妻精品| 精品无人区乱码1区二区| 日本黄色片子视频| 99久久成人亚洲精品观看| 成人美女网站在线观看视频| 国产伦精品一区二区三区四那| 亚洲中文日韩欧美视频| 日韩国内少妇激情av| 国产精品一区二区三区四区免费观看 | av女优亚洲男人天堂| 午夜激情欧美在线| 色av中文字幕| 大又大粗又爽又黄少妇毛片口| 伦精品一区二区三区| av天堂中文字幕网| 国产精品综合久久久久久久免费| 精品久久久久久久久亚洲 | 日本三级黄在线观看| 亚洲三级黄色毛片| 欧美精品啪啪一区二区三区| 永久网站在线| 少妇人妻一区二区三区视频| 国产私拍福利视频在线观看| 淫秽高清视频在线观看| 国产精品久久视频播放| 少妇人妻精品综合一区二区 | 成人国产综合亚洲| 色哟哟哟哟哟哟| 国内揄拍国产精品人妻在线| 国产精品自产拍在线观看55亚洲| 亚洲精品久久国产高清桃花| 一进一出抽搐动态| 精品人妻一区二区三区麻豆 | 偷拍熟女少妇极品色| 偷拍熟女少妇极品色| av女优亚洲男人天堂| 女人被狂操c到高潮| 国产aⅴ精品一区二区三区波| 亚洲国产日韩欧美精品在线观看| 国产老妇女一区| 自拍偷自拍亚洲精品老妇| 在线免费观看的www视频| 国产精品日韩av在线免费观看| 一区二区三区四区激情视频 | 成人永久免费在线观看视频| 亚洲电影在线观看av| 国产一区二区三区av在线 | 成年女人毛片免费观看观看9| 亚洲在线观看片| 91在线观看av| 一个人看的www免费观看视频| 久久久久久九九精品二区国产| 日本色播在线视频| 久久精品影院6| 免费搜索国产男女视频| 国产成人aa在线观看| 亚洲第一电影网av| 欧美高清成人免费视频www| av.在线天堂| 美女大奶头视频| 一级av片app| 色5月婷婷丁香| 中文字幕av成人在线电影| 男女边吃奶边做爰视频| 午夜激情欧美在线| 亚洲成av人片在线播放无| 自拍偷自拍亚洲精品老妇| 欧美日韩乱码在线| 成人二区视频| 亚洲,欧美,日韩| 一区福利在线观看| 久久亚洲真实| 天堂影院成人在线观看| 亚洲成人免费电影在线观看| 九色成人免费人妻av| 少妇人妻精品综合一区二区 | 搡女人真爽免费视频火全软件 | 琪琪午夜伦伦电影理论片6080| 国产中年淑女户外野战色| 久久久久久久久久成人| 久久久久久久久久黄片| 久久精品国产亚洲网站| 国产亚洲精品久久久久久毛片| av在线天堂中文字幕| 99热这里只有精品一区| 亚洲色图av天堂| 亚洲精品色激情综合| 国产精品美女特级片免费视频播放器| 日本熟妇午夜| 国国产精品蜜臀av免费| 精品午夜福利视频在线观看一区| 悠悠久久av| 成人午夜高清在线视频| 国产高清三级在线| 大型黄色视频在线免费观看| 最近中文字幕高清免费大全6 | 看十八女毛片水多多多| 国产精品一及| 亚洲成人久久爱视频| 欧美日韩瑟瑟在线播放| 日本精品一区二区三区蜜桃| 女人被狂操c到高潮| 国产高清视频在线播放一区| 欧美成人a在线观看| 久久久久免费精品人妻一区二区| 在线免费观看的www视频| 男人狂女人下面高潮的视频| 我的老师免费观看完整版| 亚洲精华国产精华精| 欧美激情国产日韩精品一区| 国产欧美日韩精品亚洲av| 特大巨黑吊av在线直播| 男人舔奶头视频| 18禁裸乳无遮挡免费网站照片| 日本熟妇午夜| 久久久久久久午夜电影| 别揉我奶头 嗯啊视频| 午夜免费成人在线视频| 国内精品宾馆在线| 国产精品98久久久久久宅男小说| 一级毛片久久久久久久久女| 成人美女网站在线观看视频| av国产免费在线观看| 亚洲成人久久爱视频| 亚洲第一区二区三区不卡| 久久亚洲真实| 欧美激情国产日韩精品一区| 亚洲av免费在线观看| 97超级碰碰碰精品色视频在线观看| 午夜福利在线观看免费完整高清在 | 99精品在免费线老司机午夜| 国产精品,欧美在线| 欧美日韩亚洲国产一区二区在线观看| 国产男靠女视频免费网站| 1024手机看黄色片| 国产精品久久久久久久电影| 亚洲aⅴ乱码一区二区在线播放| 国产精品99久久久久久久久| 亚洲最大成人av| 真人做人爱边吃奶动态| 九九久久精品国产亚洲av麻豆| 国产中年淑女户外野战色| 亚洲人成网站在线播放欧美日韩| 91久久精品国产一区二区三区| 国内久久婷婷六月综合欲色啪| 亚洲精品在线观看二区| 桃色一区二区三区在线观看| 搡老岳熟女国产| 免费看日本二区| 成人三级黄色视频| 特大巨黑吊av在线直播| 久久亚洲精品不卡| 欧美高清性xxxxhd video| 男女下面进入的视频免费午夜| 内地一区二区视频在线| 久久久久久伊人网av| 级片在线观看| 久久精品影院6| 99热精品在线国产| 亚洲最大成人av| 在线观看舔阴道视频| 免费观看人在逋| 欧美绝顶高潮抽搐喷水| 亚洲aⅴ乱码一区二区在线播放| 国产乱人伦免费视频| 亚洲精品影视一区二区三区av| av天堂中文字幕网| 春色校园在线视频观看| 麻豆久久精品国产亚洲av| 乱系列少妇在线播放| 国内少妇人妻偷人精品xxx网站| 可以在线观看的亚洲视频| 成人毛片a级毛片在线播放| 精品欧美国产一区二区三| 国内精品美女久久久久久| 麻豆国产av国片精品| h日本视频在线播放| 欧美一区二区国产精品久久精品| 国产又黄又爽又无遮挡在线| 高清在线国产一区| 色在线成人网| 综合色av麻豆| 能在线免费观看的黄片| 伦理电影大哥的女人| 少妇裸体淫交视频免费看高清| 久久久久久久久久成人| 亚洲av日韩精品久久久久久密| 美女 人体艺术 gogo| 亚洲真实伦在线观看| 亚洲最大成人手机在线| 1000部很黄的大片| 最近最新免费中文字幕在线| 成年女人看的毛片在线观看| 欧美成人a在线观看| 又紧又爽又黄一区二区| 99久久无色码亚洲精品果冻| 成人特级黄色片久久久久久久| 男人和女人高潮做爰伦理| 国产精品一区二区三区四区免费观看 | 动漫黄色视频在线观看| 99精品久久久久人妻精品| 91在线观看av| 国产精品自产拍在线观看55亚洲| 不卡视频在线观看欧美| 一边摸一边抽搐一进一小说| 欧美不卡视频在线免费观看| 国产成人影院久久av| 久久久久九九精品影院| 亚洲avbb在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 露出奶头的视频| 又黄又爽又刺激的免费视频.| 长腿黑丝高跟| 99热这里只有是精品50| 久久久久久久午夜电影| 久久久精品大字幕| 精品久久国产蜜桃| 精品一区二区三区av网在线观看| 男女做爰动态图高潮gif福利片| 欧美丝袜亚洲另类 | av视频在线观看入口| h日本视频在线播放| 国内毛片毛片毛片毛片毛片| 国产在线男女| 九九热线精品视视频播放| 久久国产精品人妻蜜桃| 精品久久国产蜜桃| 真实男女啪啪啪动态图| 亚洲乱码一区二区免费版| 一级av片app| 国产三级中文精品| 欧美成人免费av一区二区三区| 大型黄色视频在线免费观看| 蜜桃亚洲精品一区二区三区| 免费看a级黄色片| 麻豆国产av国片精品| 日韩精品青青久久久久久| 成人av在线播放网站| 亚洲欧美日韩无卡精品| aaaaa片日本免费| 黄片wwwwww| 在线看三级毛片| 搡老岳熟女国产| 欧美3d第一页| 亚洲av中文av极速乱 | 国产精品一区二区性色av| 久99久视频精品免费| 亚洲第一区二区三区不卡| 亚洲成人久久爱视频| av.在线天堂| 嫁个100分男人电影在线观看| av天堂中文字幕网| 国产高清有码在线观看视频| 熟女电影av网| 国模一区二区三区四区视频| 成人国产一区最新在线观看| 亚洲国产色片| 国产 一区精品| 麻豆精品久久久久久蜜桃| 色噜噜av男人的天堂激情| 国产69精品久久久久777片| 丰满的人妻完整版| 久久久久久久久久久丰满 | 国产av不卡久久| 国产在线男女| 精品一区二区免费观看| 国产真实乱freesex| 国内精品一区二区在线观看| 亚洲av第一区精品v没综合| 国产成人a区在线观看| 综合色av麻豆| 久久久久久久久中文| 91在线精品国自产拍蜜月| 欧美+日韩+精品| 久久久久九九精品影院| h日本视频在线播放| 精品欧美国产一区二区三| 免费观看在线日韩| 99在线视频只有这里精品首页| 身体一侧抽搐| 亚洲欧美清纯卡通| 欧美日本亚洲视频在线播放| 亚洲成人久久爱视频| 3wmmmm亚洲av在线观看| 又爽又黄无遮挡网站| 一进一出好大好爽视频| 亚洲图色成人| 亚洲无线在线观看| 亚洲av成人精品一区久久| 亚洲最大成人av| 蜜桃久久精品国产亚洲av| 男女啪啪激烈高潮av片| 99久国产av精品| 国产精品福利在线免费观看| 尤物成人国产欧美一区二区三区| 精品一区二区免费观看| 国产精品野战在线观看| 国产一区二区激情短视频| 亚洲国产精品合色在线| 观看美女的网站| 国产午夜福利久久久久久| 男插女下体视频免费在线播放| 成人国产综合亚洲| 精品欧美国产一区二区三| 国产高清视频在线播放一区| 极品教师在线视频| av在线观看视频网站免费| 一本精品99久久精品77| 国产一区二区三区在线臀色熟女| 变态另类成人亚洲欧美熟女| 毛片一级片免费看久久久久 | 一夜夜www| 国产精品亚洲美女久久久| 午夜老司机福利剧场| 神马国产精品三级电影在线观看| 色视频www国产| 久久精品国产清高在天天线| 国产在视频线在精品| 九九热线精品视视频播放| 看十八女毛片水多多多| 国产精品久久视频播放| 嫩草影视91久久| 午夜福利成人在线免费观看| 久久99热这里只有精品18| 久久久久久久久中文| 亚洲精华国产精华精| 中文亚洲av片在线观看爽| 热99re8久久精品国产| 2021天堂中文幕一二区在线观| 日本一本二区三区精品| 久久亚洲精品不卡| 国产av一区在线观看免费| 俄罗斯特黄特色一大片| 亚洲国产精品合色在线| 国产精品久久视频播放| 午夜影院日韩av| 国产精华一区二区三区| 99国产极品粉嫩在线观看| 日本一本二区三区精品| 变态另类丝袜制服| 老熟妇乱子伦视频在线观看| 国内精品久久久久久久电影| 久久久午夜欧美精品| 久久久久精品国产欧美久久久| 综合色av麻豆| 熟妇人妻久久中文字幕3abv| 精品欧美国产一区二区三| 一级a爱片免费观看的视频| 欧洲精品卡2卡3卡4卡5卡区| 淫秽高清视频在线观看| 国产老妇女一区| 成人性生交大片免费视频hd| 日韩精品有码人妻一区| 日韩在线高清观看一区二区三区 | 啦啦啦观看免费观看视频高清| 国产精品,欧美在线| 日韩国内少妇激情av| 麻豆国产av国片精品| 神马国产精品三级电影在线观看| 久久国产精品人妻蜜桃| 午夜亚洲福利在线播放| 成年版毛片免费区| 国内精品久久久久久久电影| 免费不卡的大黄色大毛片视频在线观看 | 制服丝袜大香蕉在线| 国内精品美女久久久久久| 1024手机看黄色片| 日韩精品有码人妻一区| 18+在线观看网站| 亚洲av第一区精品v没综合| 日本黄大片高清| 免费黄网站久久成人精品| 91麻豆精品激情在线观看国产| 欧美激情国产日韩精品一区| 亚洲欧美日韩卡通动漫| 日韩人妻高清精品专区| av黄色大香蕉|