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

    土壤中銅生態(tài)閾值的影響因素及其預(yù)測模型

    2014-04-26 05:58:55王小慶李菊梅韋東普陳世寶馬義兵黃占斌洛陽理工學(xué)院環(huán)境工程與化學(xué)系河南洛陽4710中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所國家土壤肥力與肥料效益監(jiān)測站網(wǎng)北京100081中國礦業(yè)大學(xué)北京化學(xué)與環(huán)境工程學(xué)院北京10008
    中國環(huán)境科學(xué) 2014年2期
    關(guān)鍵詞:中銅毒理學(xué)老化

    王小慶,李菊梅,韋東普,陳世寶,馬義兵*,黃占斌(1.洛陽理工學(xué)院環(huán)境工程與化學(xué)系,河南 洛陽4710;.中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,國家土壤肥力與肥料效益監(jiān)測站網(wǎng),北京 100081;.中國礦業(yè)大學(xué)(北京)化學(xué)與環(huán)境工程學(xué)院,北京 10008)

    土壤中銅生態(tài)閾值的影響因素及其預(yù)測模型

    王小慶1,2,李菊梅2,韋東普2,陳世寶2,馬義兵2*,黃占斌3(1.洛陽理工學(xué)院環(huán)境工程與化學(xué)系,河南 洛陽471023;2.中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,國家土壤肥力與肥料效益監(jiān)測站網(wǎng),北京 100081;3.中國礦業(yè)大學(xué)(北京)化學(xué)與環(huán)境工程學(xué)院,北京 100083)

    利用中國土壤的銅毒理學(xué)數(shù)據(jù)通過物種敏感性分布法推導(dǎo)了土壤中銅的不同風險水平(P%)的毒害濃度值(P% Hazardous concentration, HCp),并利用淋洗-老化因子校正HCp獲得不同土地類型土壤中銅的生態(tài)閾值,探討了土壤性質(zhì)對銅生態(tài)閾值的影響并建立了兩者之間的量化關(guān)系.結(jié)果表明,土壤性質(zhì)對銅的生態(tài)閾值有顯著影響,土壤pH值和陽離子交換量(CEC)是影響土壤銅生態(tài)閾值的最主要因子,可達銅生態(tài)閾值變異的80%以上.基于土壤pH值和CEC的兩因子回歸模型能較好地預(yù)測銅生態(tài)閾值,其決定系數(shù)R2為0.820~0.913;增加土壤有機碳含量(OC)的三因子模型具有更高的準確性,其決定系數(shù)R2為0.852~0.988.研究結(jié)果可為科學(xué)合理地進行土壤中銅的生態(tài)風險評價和建立土壤環(huán)境質(zhì)量標準提供依據(jù).

    生態(tài)閾值;銅;土壤性質(zhì);淋洗-老化因子;預(yù)測模型

    土壤生態(tài)閾值是指為了對陸地生物及關(guān)鍵的土壤生態(tài)功能提供適當?shù)谋Wo而制定的土壤中污染物的濃度限值[1-2],其確定方法中近來最常用的是物種敏感性分布法(SSD)[3-5],該方法是利用累積概率分布函數(shù)擬合污染物的毒理學(xué)數(shù)據(jù)建立其物種敏感性分布曲線,依據(jù)不同的保護程度(風險水平) 獲取曲線上不同百分點所對應(yīng)的濃度值作為基準值[4,6],其風險水平的選取通常依據(jù)土地類型而定,如農(nóng)業(yè)用地選取 5%處所對應(yīng)的濃度HC5(5%毒害濃度),即保護95%生物物種的限量值;商業(yè)用地選取保護 60%物種的限量值HC40等.利用SSD法建立的土壤生態(tài)閾值不同于傳統(tǒng)的土壤環(huán)境質(zhì)量標準,其推導(dǎo)過程中考慮了物種敏感性、土壤性質(zhì)、生物有效性及污染來源等因素的差異,且可依據(jù)不同的風險水平選取相應(yīng)的限量值,具有科學(xué)性、基礎(chǔ)性和區(qū)域性的特點[7-8].

    銅(Cu)作為植物必需的微量營養(yǎng)元素在植物生理代謝上起著重要作用,但若土壤中Cu過量時,植物會吸收土壤中過量的 Cu2+并將其大部分淀積于根部表層,影響植物對其他養(yǎng)分的吸收,從而對植物產(chǎn)生嚴重毒害[9],Cu還可被植物吸收積累在可食部位,對人類健康形成潛在危害[10].目前工農(nóng)業(yè)生產(chǎn)活動導(dǎo)致的土壤 Cu污染已經(jīng)引起關(guān)注

    [11],準確科學(xué)地對土壤中的Cu進行風險識別和評價是土壤污染治理的重要任務(wù)之一.我國當前的土壤環(huán)境質(zhì)量標準主要基于對人類健康的保護,未有明確的生態(tài)環(huán)境質(zhì)量標準值[12],不利于農(nóng)業(yè)的可持續(xù)發(fā)展.該標準中的銅限量值僅依據(jù)土壤pH值進行簡單劃分[13],而銅對土壤生物的毒性大小與土壤性質(zhì)密切相關(guān)[14-16],進入土壤的時間長短對其毒性亦有著不容忽視的影響[17],即外源進入的Cu污染物與背景值部分的活性差異很大.

    基于不同物種對于污染物敏感性差異提出的物種敏感性分布法(SSD)[3],已被多個國家及機構(gòu)用于環(huán)境質(zhì)量基準的制定.我國關(guān)于SSD法的研究還處在起步階段,僅部分學(xué)者構(gòu)建了水體環(huán)境中的 DDT、林丹、多環(huán)芳烴等有機污染物的SSD曲線,且其毒理學(xué)數(shù)據(jù)多來自于國外的數(shù)據(jù)庫[18-21],本課題組利用此方法制定了中國土壤中鎳的生態(tài)閾值并對其影響因素進行了深入探討

    [22].本研究利用基于中國土壤的 21個物種的Cu毒理學(xué)數(shù)據(jù),結(jié)合Cu的生物毒害預(yù)測模型,應(yīng)用物種敏感性分布法推導(dǎo)出土壤中銅的 HCp(p為不同風險水平,如5、20等),并利用淋洗-老化因子對其進行校正獲得以外源添加量表示的銅生態(tài)閾值;探討了土壤性質(zhì)對銅生態(tài)閾值的影響并建立了兩者之間的量化關(guān)系,以期為中國土壤中銅的污染風險評價及其生態(tài)環(huán)境質(zhì)量標準的制定提供了科學(xué)依據(jù).

    1 材料與方法

    1.1 銅毒理學(xué)數(shù)據(jù) (EC10)的收集與篩選

    在中國知網(wǎng)(CNKI)數(shù)據(jù)庫中以摘要中有“土壤”和“銅”為搜索條件,收集所有基于中國土壤的銅毒理學(xué)數(shù)據(jù)并進行篩選,同時匯總有關(guān)研究項目中獲得的銅毒理學(xué)數(shù)據(jù)[23].數(shù)據(jù)首選 EC10(一定時間間隔內(nèi)引起 10%抑制效應(yīng)的濃度),利用文獻中報道的結(jié)果或通過有明顯劑量-效應(yīng)關(guān)系的原始數(shù)據(jù)計算.篩選后的數(shù)據(jù)應(yīng)滿足以下條件:試驗程序及數(shù)據(jù)處理規(guī)范,測試終點為植物、無脊椎動物、微生物或微生物過程的相應(yīng)評價終點,實驗結(jié)果以銅的總量(mg/kg以干重計)表示,同時附有必要的土壤性質(zhì)參數(shù)等.對于同一物種存在多個數(shù)據(jù)時,通過計算其幾何均值形成“物種平均”EC10或NOEC值(無效應(yīng)觀測濃度)作為最終值.篩選出的21個物種(19種植物、2種微生物)的生態(tài)毒理學(xué)數(shù)據(jù)均為基于外源添加量表示的毒性閾值EC10add,不包括試驗土壤本底值,除特別標明外,本文中的EC10均為EC10add.

    1.2 淋洗因子計算

    由于 EC50(一定時間間隔內(nèi)引起 50%抑制效應(yīng)的濃度)處于劑量-效應(yīng)“S”型曲線的中部,不受個別數(shù)據(jù)點所導(dǎo)致的斜率變化的影響,比EC10誤差更小[24],因此利用“亞洲金屬”項目中的數(shù)據(jù),分別計算大麥根伸長、西紅柿及小白菜生長3個評價終點在17種淋洗土壤[20]中的EC50與同類型未淋洗土壤中的EC50的比值,取3個比值的算術(shù)均值作為該土壤的淋洗因子(LF)值.并將LF與對應(yīng)土壤的理化性質(zhì)參數(shù)進行多元線性回歸分析建立數(shù)值關(guān)系.

    1.3 銅毒理學(xué)數(shù)據(jù)歸一化及HCp值推導(dǎo)

    由于本研究中所選用的銅生物毒性預(yù)測模型均基于淋洗土壤建立,來源于未淋洗土壤的毒理學(xué)數(shù)據(jù)應(yīng)利用淋洗因子(LF)對其進行校正:即取基于未淋洗土壤獲得的毒性閾值與試驗土壤理化性質(zhì)參數(shù)對應(yīng)的LF值的乘積作為對應(yīng)類型淋洗土壤的毒性閾值.利用李波等[25]和韋東普等

    [26]建立的銅生物毒性預(yù)測模型將21個物種的EC10(基于淋洗土壤或經(jīng)淋洗因子校正)歸一化到不同土壤條件下,獲得對應(yīng)不同土壤條件的銅毒理學(xué)數(shù)據(jù)組,具體方法見參考文獻[8].利用BurrIII分布函數(shù)擬合經(jīng)歸一化后的銅 EC10,建立不同土壤條件下的物種敏感性分布曲線,根據(jù)SSD曲線確定HCp,即在此濃度下生境中(1-p)%的生物是(相對)安全的[4,23].HCp依據(jù)土地類型分別為:國家自然保護區(qū) HC1、農(nóng)林用地(包括草地)HC5、住宅區(qū)和公園等綠化用地HC20、商業(yè)及工礦業(yè)用地HC40.

    土壤pH值、陽離子交換量(CEC)、有機碳含量(OC)是影響土壤中銅的生物有效性/毒性及老化的主要影響因素[28].因此,設(shè)定以下土壤條件:pH 值 4.5~9、陽離子交換量(CEC) 5~30cmol/kg、有機碳含量(OC)1%~4%,推導(dǎo)出對應(yīng)不同土壤性質(zhì)的銅HC5值以確定土壤中銅生態(tài)閾值的主要影響因子.其中土壤 pH值在水土比為1:5條件下振蕩1h,靜置30min后測定[29],OC含量為總碳與無機碳含量之差[30],CEC 基于NH4Cl緩沖法測定[31].

    1.4 銅HCp值老化校正

    實驗室毒理學(xué)試驗結(jié)果與野外田間存在的差異來自于淋洗和老化的雙重作用[24],實驗室人工污染法獲得的毒理學(xué)數(shù)據(jù)應(yīng)同時經(jīng)淋洗因子和老化因子校正后才能更接近野外田間實際污染情況.目前尚無基于中國土壤的 Cu老化模型,采用基于歐洲土壤推導(dǎo)的 Cu老化因子模型[32],在老化模型研究中利用同位素稀釋技術(shù)確定添加到土壤中外源金屬經(jīng)過長時間老化后的活性(E值)變化,建立了以老化時間(t, d)和土壤pH值(1:5CaCl2)為參數(shù)的Cu老化因子計算模型(1).

    利用老化因子模型可將實驗室內(nèi)的短期試驗結(jié)果校正到一定老化時間的值,由于大多數(shù)毒理學(xué)試驗基于實驗室老化兩周(14d)后的土壤進行,本研究中老化因子(AF)取值為 14d E值與360d E值的比值,將HCp校正至老化時間為1a (360d)的老化HCp值,以老化HCp值作為銅的生態(tài)閾值.本文中除特別標明外 HCp值均為老化HCp,其老化時間為1a.

    2 結(jié)果與分析

    2.1 淋洗及老化因子

    17種土壤[25]中銅的淋洗因子(LF)值在1.05~2.78之間,表明淋洗處理能降低銅的生物毒害,在高pH值土壤中表現(xiàn)更顯著,潮土中的LF值最大.將淋洗因子(LF)值與土壤性質(zhì)參數(shù)之間進行多元線性回歸分析,獲得 LF的預(yù)測公式.式(2)~式(4)分別適用于pH<7.0、pH 7.0~8.5和pH >8.5的土壤,式中 CEC為土壤陽離子交換量(cmol/kg),OC為土壤有機碳含量(%),Clay為土壤黏粒含量(%).

    上述的LF預(yù)測公式能較為準確地通過土壤pH值、CEC、Clay計算出LF值,但影響LF值的土壤性質(zhì)參數(shù)除pH值、CEC、OC、Clay外,還包括其他因素,如土壤電導(dǎo)率(EC)、生物種類等,且很多有關(guān)銅毒性閾值研究的報道中缺少黏土含量(Clay)參數(shù).為簡化起見,本研究中的LF依據(jù)pH分段取值為上述3個評價終點計算的LF值算術(shù)均值,即 pH≤7.0時 LF=1.4,pH=7~8.5時LF=1.7,pH≥8.5時LF=1.9.

    進入土壤中的重金屬的生物有效性/毒性會因重金屬在土壤中的微孔擴散、黏土礦物的層間固定、土壤中成核/沉淀、礦物表面氧化及與土壤中固相成分形成配合物等過程顯著降低[17,32],即老化效應(yīng).因此在利用實驗室獲得的毒理學(xué)數(shù)據(jù)推導(dǎo)銅的土壤生態(tài)閾值時必須要考慮老化效應(yīng).如前所述,老化因子(AF)值取14d的活性值與360d活性值比值.AF值在1.22~1.32之間,老化效應(yīng)隨著土壤 pH值的升高而略有增強,但總體變化不顯著.

    2.2 土壤性質(zhì)對HCp值的影響

    圖1 不同土壤條件下SSD曲線Fig.1 SSD curves in different soils (EC10 was corrected with L-A/F) EC10為經(jīng)淋洗-老化因子校正的老化EC10

    應(yīng)用SSD法推導(dǎo)出不同土壤條件下的p%毒害濃度值并利用 AF值對其進行校正,獲得老化HCp值.不同土壤條件下的SSD曲線(圖1)表明:隨著土壤 pH值或 CEC含量的增大,土壤中銅SSD曲線向X軸數(shù)值大的方向(右)偏移,HCp值隨著土壤pH值、CEC含量的增高而增大;OC含量對銅SSD曲線分布有一定的影響,隨著OC含量的增加SSD曲線亦逐步向X軸數(shù)值大的方向移動,但隨著OC含量的增加,其影響程度降低,且 OC對銅 SSD曲線的影響程度在低累積概率范圍內(nèi)(<10%)較小.在土壤 OC值一定的條件下, HCp值隨著pH值和CEC的增加顯著增大,如當OC=2%,pH=5.0,CEC=10cmol/kg的條件下,HC20為 34.72mg/kg;pH=8.0,CEC=10cmol/kg的條件下,HC20為65.37mg/kg;pH=5.0,CEC= 20cmol/kg的條件下,HC20則為57.66mg/kg.但當土壤pH和CEC一定時,隨著土壤OC含量的增加,HC5值變化不顯著,HC20和 HC40略有增加.我國土壤有機碳含量普遍較低,且土壤中 OC的含量變化不大,如李波等[25]采集的中國 17種典型農(nóng)田土壤樣品OC的變化范圍為0.6%~4.3%,除北京的棕壤外,其余16個土壤樣品的OC均小于3.0%[25,33],土壤中銅HCp值的主要影響因素為pH和CEC.

    2.3 HCp值的預(yù)測模型

    利用Excel以∑(預(yù)測HCp-推導(dǎo)HCp)2值最小為條件進行規(guī)劃求解,量化HCp值與對應(yīng)土壤性質(zhì)參數(shù)(pH值、CEC、OC)之間的關(guān)系,獲得HCp值的預(yù)測模型,如表1所示.

    表1中HCp為經(jīng)老化-淋洗因子校正后的老化HCp值,pH值為土水比為1:5時的pH[29],OC為土壤有機碳含量(%)[30],CEC為陽離子可交換量(cmol/kg)[31],從表中HCp的各預(yù)測模型中土壤性質(zhì)參數(shù)斜率的大小可看出:土壤pH值和CEC是影響土壤銅生態(tài)閾值的主要因子.三因子模型和兩因子模型均能較好地依據(jù)土壤性質(zhì)參數(shù)預(yù)測銅老化HCp值,其預(yù)測值與推導(dǎo)值之間的相關(guān)系數(shù)(R2)分別為0.820~0.913和0.852~0.988.兩因子模型在 CEC較高的堿性(pH≥7.5)土壤中預(yù)測效果較差,HCp預(yù)測值小于實際推導(dǎo)值;三因子模型只在土壤pH(≥8.5)和CEC(≥25cmol/kg)兩個土壤性質(zhì)參數(shù)值都大的極端情況下其預(yù)測效果較差,其預(yù)測值小于實際推導(dǎo)值.本研究建立的模型中CEC為1mol/L中性(pH=7)NH4Cl淋洗法測定值,其值大小與測定土壤的pH值和提取劑pH值(7.0)的差值大小密切相關(guān),中國17種代表性農(nóng)田土壤的性質(zhì)參數(shù)數(shù)據(jù)表明在土壤pH值非常高的情況下,其CEC值通常低于20cmol/kg.但將這兩種模型用于堿性土壤時仍需謹慎,預(yù)測公式計算的HCp值可能會導(dǎo)致過度保護,即當土壤中銅濃度達到預(yù)測公式計算的HCp值時,其受影響的物種少于(1-p)%.

    表1 不同土地類型的土壤銅生態(tài)閾值預(yù)測模型Table 1 Predictable models of ecological thresholds for copper in soils for different land types

    土壤pH值和CEC可控制HC1、HC5、HC20、HC40變量的80%以上,而logOC分別控制HC1、HC5、HC20、HC40變量的3.2%、6.7%、8.7%和7.5%,且土壤中OC的含量變化不大[25,34].有關(guān)污染土壤中銅賦存形態(tài)的研究表明由于吸附及螯合作用,土壤有機物會促進土壤中銅向穩(wěn)定形態(tài)轉(zhuǎn)變,從而降低銅的生物有效性[35],而有關(guān)銅的生物有效性研究則表明,銅的生物有效性與土壤有機碳含量成正相關(guān)[36].銅在土壤中的化學(xué)形態(tài)分配受土壤中總銅含量、pH、OC、黏粒礦物、養(yǎng)分和CEC等多因素的影響[37-39],我國大部分土壤黏粒礦物具有可變電荷表面且有機質(zhì)含量低,土壤CEC受土壤pH值、OC、土壤黏土類型及數(shù)量的共同影響,所以CEC對銅生態(tài)閾值的顯著影響可能是土壤 pH、OC等其他因素對其影響的體現(xiàn).在實際應(yīng)用中,若土壤性質(zhì)參數(shù)OC缺失,可直接利用兩因子模型(R2=0.820~0.913) 或?qū)⑼寥繭C設(shè)一個缺省值(如1%)來預(yù)測不同土壤中的銅生態(tài)閾值.

    此處建立的HCp預(yù)測模型體現(xiàn)了土壤性質(zhì)對銅生態(tài)閾值的影響,量化了土壤性質(zhì)參數(shù)與銅生態(tài)閾值的數(shù)值關(guān)系,也說明了建立土壤性質(zhì)對應(yīng)的污染物生態(tài)閾值的必要性.可基于上述兩個模型計算出的對應(yīng)土壤性質(zhì)的銅生態(tài)閾值,再利用前述老化因子的校正方法將其校正至實際污染時間,為田間土壤實際污染的生態(tài)風險評價提供科學(xué)依據(jù),同時也為制定適用性更強的土壤銅生態(tài)環(huán)境質(zhì)量標準提供科學(xué)依據(jù).

    2.4 銅生態(tài)閾值預(yù)測模型的田間驗證

    研究者們在“亞洲金屬”項目中為了驗證實驗室獲得的研究結(jié)果能否用于田間作物,選取湖南祁陽的酸性土(紅壤,26°45′N,111°52′E,pH5.3, CEC7.47,OC0.87%)、山東德州(潮土,37°20′N, 116°38′E,pH8.9,CEC 8.33,OC0.69%)的堿性土及浙江嘉興的中性土(水稻土,30°77'N,120o76'E,pH6.7, CEC19.33,OC1.42%)進行了連續(xù)兩年的田間試驗.將本研究建立的HC5預(yù)測模型的計算值與對應(yīng)土壤類型中不同田間作物的實際EC10值進行比較,其結(jié)果如圖2所示.除2008年的嘉興中性土壤中的田間油菜外,其他各物種的 EC10值均處于1:1線的上方,表明基于實驗室毒理學(xué)數(shù)據(jù)推導(dǎo)出的HC5值較田間作物的毒性閾值小.而本課題前期的有關(guān)鎳生態(tài)閾值的研究中也發(fā)現(xiàn),2008年的田間油菜鎳EC10小于其相應(yīng)土壤條件下的鎳HC5值[22],故推測可能是由于氣候條件、生產(chǎn)管理等其他因素導(dǎo)致其銅、鎳EC10異常偏小而導(dǎo)致的結(jié)果.

    圖2 HC5值與田間作物毒性閾值(EC10)的比較Fig.2 Comparison of HC5 derived from with the ecotoxicity data from different field crops in three field sites

    重金屬污染多以非可溶鹽的形式進入到土壤里,如污泥、礦渣及顆粒的大氣沉降等,如 Luo等[40]的研究表明中國農(nóng)田土壤中的銅污染輸入量中 69%來自于畜禽糞便,18%來自于含銅顆粒的大氣沉降,以這類途徑進入土壤中的重金屬的生物有效性較低,其對生物的毒害遠小于可溶性金屬鹽.且當金屬進入土壤后,與土壤發(fā)生吸附、解吸、絡(luò)合、沉淀和溶解等一系列反應(yīng),形成不同的化學(xué)形態(tài),有研究表明土壤中的銅多以穩(wěn)定的殘渣態(tài)存在[35-36],本研究中用于推導(dǎo)銅生態(tài)閾值的毒理學(xué)數(shù)據(jù)是利用實驗室外源添加可溶性銅鹽模擬銅污染的方法獲得的,且對于風險水平的選取也較為嚴格,如農(nóng)林用地中的園地、林地、草地等均選取與耕地同等風險水平的 HC5,所以將上述的兩類預(yù)測模型應(yīng)用于具有多種銅污染來源的實際土壤是較為保守的.

    3 結(jié)論

    3.1 淋洗與老化均能降低土壤中銅的生物毒性,除個別土壤外,17種代表性的中國土壤中銅的淋洗因子(LF)均大于 1,其老化因子(AF)隨著土壤pH增大略有增強,但變化不明顯.

    3.2 土壤中銅SSD曲線隨著pH值和CEC值的增大向濃度值(X軸)增大方向移動,OC含量大小對SSD曲線的位置有一定的影響,但在累積概率較低范圍內(nèi)其影響程度較小.

    3.3 土壤銅生態(tài)閾值受土壤pH值和陽離子交換量(CEC)、有機碳含量(OC)的共同影響,CEC、pH值是主要影響因子,可控制其值變異的 80%以上.

    [1] UKEA (United Kingdom Environment Agency). Soil screening values for use in UK ecological risk assessment [R]. Warrington: United Kingdom Environment Agency, 2004.

    [2] USEPA (U.S. Environmental Protection Agency). Guidance for developing ecological soil screening levels [R]. OSWER Directive 9285.7-55. Washington D C: U.S. Environmental Protection Agency, 2003.

    [3] Kooijman S A L M. A safety factor for LC50values allowing for differences in sensitivity among species [J]. Water Research, 1987,21:269-276.

    [4] Newman M C, Ownby D R, Mezin L C A, et al. Applying species-sensitivity distributions in ecological risk assessment: Assumptions of distribution type and sufficient numbers of species [J]. Environ. Toxicol. Chem., 2000,19: 508-515.

    [5] Hall L W, Scott M C, Killen W D. Ecological risk assessment of copper and cadmium in surface waters of Chesapeake Bay watershed [J]. Environ. Toxicol. Chem., 1998,17:1172-1189.

    [6] Hose G C, Vandenbrink P J. Confirming the species-sensitivity distribution concept for Endosulfan using laboratory, mescom, and field data [J]. Arch. Environ. Contam. Toxicol., 2004,47:511-520.

    [7] 顏增光,谷慶寶,周 娟,等.構(gòu)建土壤生態(tài)篩選基準的技術(shù)關(guān)鍵及方法學(xué)概述 [J]. 生態(tài)毒理學(xué)報, 2008,3(5):417-427.

    [8] 王小慶,韋東普,黃占斌,等.物種敏感性分布在土壤中鎳生態(tài)閾值建立中的應(yīng)用研究 [J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2012,(1):92-99. [9] 邵興華,張建忠,林國衛(wèi),等.銅脅迫對油麥菜生長和土壤酶活性的影響 [J]. 中國農(nóng)學(xué)通報, 2010,26(4):157-161.

    [10] 金 進,葉亞新,李 周,等.重金屬銅對玉米的影響 [J]. 玉米科學(xué), 2006,14(3):83-86.

    [11] 林義章,徐 磊.銅污染對高等植物的生理毒害作用研究 [J].中國生態(tài)農(nóng)業(yè)學(xué)報, 2007,15(1):201-204.

    [12] 周國華,秦緒文,董巖翔.土壤環(huán)境質(zhì)量標準的制定原則與方法[J]. 地質(zhì)通報, 2005,24(8):721-727.

    [13] GB15618-1995 土壤環(huán)境質(zhì)量標準 [S].

    [14] Rooney C P, Zhao F J, Mcgrath S P. Soil factors controlling the expression of copper toxicity to plants in a wide range of European soils [J]. Environ. Toxicol. Chem., 2006,25:726-732.

    [15] Warne M S, Hhhmsbergen D A, Stevens D, et al. Modeling the toxicity of copper and zinc salts to wheat in 14soils [J]. Environ. Toxicol. Chem., 2008,27:786-792.

    [16] Li X F, Sun J W, Qiao M, et al. Copper toxicity thresholds in Chinese soils based on substrate-induced nitrification assay [J]. Environ. Toxicol. Chem., 2010,29:294-300.

    [17] Oorts K, Bronckaers H, Smolders E. Discrepancy of the microbial response to elevated Cu between freshly spiked and long-term contaminated soils. [J]. Environ. Toxicol. Chem., 2006,25:845-853.

    [18] 王 印,王軍軍,秦 寧,等.應(yīng)用物種敏感性分布評估DDT和林丹對淡水生物的生態(tài)風險 [J]. 環(huán)境科學(xué)學(xué)報, 2009:29(11): 2407-2414.

    [19] 劉 良,顏小品,王 印,等.應(yīng)用物種敏感性分布評估多環(huán)芳烴對淡水生物的生態(tài)風險 [J]. 生態(tài)毒理學(xué)報, 2009,4(5):647-654.

    [20] 蔣丹烈,尹大強.應(yīng)用物種敏感性分布法對太湖沉積物中多環(huán)芳烴的生態(tài)風險分析 [J]. 生態(tài)毒理學(xué)報, 2011,6(1):60-66.

    [21] 陳波宇,鄭斯瑞,牛希成,等.水生生物對三唑磷的物種敏感度分布研究 [J]. 環(huán)境科學(xué), 2011,32(4):1101-1107.

    [22] 王小慶,馬義兵,黃占斌.土壤中鎳生態(tài)閾值的影響因素及預(yù)測模型 [J]. 農(nóng)業(yè)工程學(xué)報, 2012,28(5):220-225.

    [23] Ma Y B, Mclaughlin M J, Zhu Y G, et al. Final Report for Metals in Asia [R]. CSIRO Land and Water, Beijing, 2009.

    [24] Smolders E, Oorts K, Vansprang P, et al. Toxicity of trace metals in soil as affected by soil type and aging after contamination: calibrated bioavailability models to set ecological soil standards [J]. Environ. Toxicol. Chem., 2009,28(8):1633-1642.

    [25] 李 波.外源重金屬銅、鎳的植物毒害及預(yù)測模型研究 [D]. 北京:中國農(nóng)業(yè)科學(xué)院, 2010,25-42.

    [26] 韋東普.應(yīng)用發(fā)光細菌法測定中國土壤中銅、鎳毒性的研究[D]. 北京:中國農(nóng)業(yè)科學(xué)院, 2010,75-76.

    [27] 金小偉,雷炳莉,許宜平,等.水生態(tài)基準方法學(xué)概述及建立我國水生態(tài)基準的探討 [J]. 生態(tài)毒理學(xué)報, 2009,4(5):609-616.

    [28] Lock k, Janssen C. Influence of aging on copper bioavailability in soils [J]. Environ. Toxicol. Chem., 2003,22(5):1162-1166.

    [29] Rayment G E, Higginson F R. Ion-exchange properties [M]// Australian laboratory handbook of soil and water chemical methods. Inkata press, Melbourne, VIC, Australia, 1992,137-194.

    [30] 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.

    [31] Gillman G P, Skjemstad J O, Bruce R C. A comparison of methods used in Queensland for determining cation exchange properties [R]. CSIRO Australia, Division of Soils, 1982,44-45.

    [32] Ma Y B, Lombi E, Oliver I W, et al. Long-term aging of copper added to soils [J]. Environ. Sci. Technol., 2006,40(20):6310-6317.

    [33] 張洪濤,李 波,劉繼芳,等.西紅柿鎳毒害的土壤主控因子和預(yù)測模型研究 [J]. 生態(tài)毒理學(xué)報, 2009,4(4):569-576.

    [34] 吳樂知,蔡祖聰.中國土壤有機質(zhì)含量變異性與空間尺度的關(guān)系[J]. 地球科學(xué)進展, 2006,21(9):965-972.

    [35] 郭 平.長春市土壤重金屬污染機理與防治對策研究 [D]. 吉林大學(xué),長春, 2005,67-73.

    [36] 丁 園,宗良綱,舒紅英,等.土壤銅鎘的植物可利用性評估模型[J]. 環(huán)境科學(xué)與技術(shù), 2011,34(2):188-192.

    [37] Conesa H M, Robinson B H, Schulin R, et al. Growth of Lygeum spartum in acid mine tailings: response of plants developed from seedlings, rhizomes and at field conditions [J]. Environmental Pollution, 2007,145(3):700-707.

    [38] Deng H, Yea Z H, Wong M H. Accumulation of lead, zinc, copper and cadmium by 12 wetland plant species thriving in metalcontaminated sites in China [J]. Environmental Pollution, 2004, 132(1):29-40.

    [39] Ju X T, Kou C L, Christie P, et al. Changes in the soil environment from excessive application of fertilizers and manures to two contrasting intensive cropping systems on the North China Plain [J]. Environmental Pollution, 2007,145(2):497-506.

    [40] Luo L, Ma Y B, Zhang S Z, et al. An inventory of trace element inputs to agricultural soils in China [J]. Journal of Environmental Management, 2009,90(8):2524-2530.

    Major soil factors affecting ecological threshold for copper and the predictable models.

    WANG Xiao-qing1,2, LI Ju-mei2,

    WEI Dong-pu2, CHEN Shi-bao2, MA Yi-bing2*, HUANG Zhan-bin3(1.Department of Environmental Engineering and Chemistry, Luoyang Institute of Science and Technology, Luoyang 471023, China;2.National Soil Fertility and Fertilizer Effects Long-term Monitoring Network, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing 100081, China;3.School of Chemical and Environmental Engineering, University of Mining and Technology-Beijing, Beijing 100083, China). China Environmental Science, 2014,34(2):445~451

    The hazardous concentrations (HCp) are the concentrations of a substance in soils which can be toxic to p% of bio-species in ecosystems and, in other words, protect (1-p)% of bio-species in ecosystem. The values of HCp are the scientific basis for setting up soil ecological quality standards. In the present study, the values of HCp for copper in soil were derived using a method of species sensitivity distribution based on copper toxicity data from Chinese soils. Correcting the derived HCp with aging factor (AF) to obtain soil ecological thresholds of copper, i e., aged HCp. Predictable models for copper ecological thresholds based on soil properties were also developed. The results showed that soil pH and cation exchange capacity were the main soil factors affecting soil ecological thresholds of copper. Soil cation exchange capacity and pH could explain the >80% variation of it. The two-factor predictable models based on soil pH and cation exchange capacity could predict copper ecological thresholds with determination coefficients (R2) from 0.820 to 0.913. While the three-factor predictable models with taking into account the effect of organic carbon on copper ecological threshold were more accurate than the two-factor models with R2from 0.852 to 0.988. The predictable models developed in the present study proved the quantitative relationship between soil properties with copper ecological thresholds and the necessity of developing and implementing soil specific or regional soil environmental quality standards.

    ecological threshold;copper;soil properties;leaching-aging factor;predictable model

    X53

    :A

    :1000-6923(2014)02-0445-07

    王小慶(1979-),女,安徽懷寧人,講師,博士,主要從事土壤重金屬生態(tài)風險評價管理及形態(tài)分析發(fā)面的研究.發(fā)表論文10余篇.

    2013-06-20

    公益性行業(yè)(農(nóng)業(yè))科研專項(200903015);國際銅業(yè)協(xié)會資助項目

    * 責任作者, 研究員, ybma@caas.ac.cn

    猜你喜歡
    中銅毒理學(xué)老化
    延緩大腦老化,要怎樣吃
    豬飼料中銅和鋅的替代來源
    節(jié)能技術(shù)在開關(guān)電源老化測試中的應(yīng)用
    電子制作(2018年10期)2018-08-04 03:24:30
    PM2.5毒理學(xué)實驗染毒方法及毒理學(xué)效應(yīng)
    杜絕初春老化肌
    Coco薇(2016年2期)2016-03-22 02:40:06
    鹽酸改性蒙脫土對廢水中銅離子的吸附研究
    災(zāi)害毒理學(xué)理論研究初探
    火災(zāi)毒理學(xué)若干問題的探討
    X65管線鋼根焊工藝中銅襯墊滲銅研究
    焊接(2015年7期)2015-07-18 10:59:15
    微波消解-分光光度法測定茶葉中銅的含量
    久久 成人 亚洲| 一边摸一边抽搐一进一出视频| 亚洲五月色婷婷综合| 无遮挡黄片免费观看| 美女高潮喷水抽搐中文字幕| 国产精品久久久久久人妻精品电影 | 亚洲成人手机| 精品人妻在线不人妻| 中文字幕另类日韩欧美亚洲嫩草| 桃花免费在线播放| 男女床上黄色一级片免费看| 汤姆久久久久久久影院中文字幕| 欧美一级毛片孕妇| 国产老妇伦熟女老妇高清| 亚洲伊人色综图| 啦啦啦 在线观看视频| 免费人妻精品一区二区三区视频| 丰满人妻熟妇乱又伦精品不卡| av电影中文网址| 99久久人妻综合| 亚洲男人天堂网一区| 免费看十八禁软件| 法律面前人人平等表现在哪些方面| 精品久久久精品久久久| 欧美精品高潮呻吟av久久| 午夜福利在线免费观看网站| 久久久精品国产亚洲av高清涩受| 一本综合久久免费| 亚洲少妇的诱惑av| 亚洲av第一区精品v没综合| 亚洲成a人片在线一区二区| 国产免费视频播放在线视频| 好男人电影高清在线观看| 捣出白浆h1v1| 人人妻人人澡人人看| 一级毛片女人18水好多| 久久精品国产a三级三级三级| 男女高潮啪啪啪动态图| 久久久久网色| 午夜福利一区二区在线看| 69av精品久久久久久 | 少妇被粗大的猛进出69影院| 嫁个100分男人电影在线观看| 伦理电影免费视频| 亚洲精品国产区一区二| 美女高潮到喷水免费观看| 久久久国产成人免费| 色视频在线一区二区三区| 另类精品久久| 午夜福利免费观看在线| 少妇的丰满在线观看| 无限看片的www在线观看| 亚洲美女黄片视频| 亚洲av第一区精品v没综合| 久久久久国产一级毛片高清牌| 夜夜骑夜夜射夜夜干| 无遮挡黄片免费观看| 少妇裸体淫交视频免费看高清 | 日韩欧美免费精品| 亚洲精品国产一区二区精华液| 国产精品熟女久久久久浪| 国产精品麻豆人妻色哟哟久久| 91成年电影在线观看| 菩萨蛮人人尽说江南好唐韦庄| 精品少妇久久久久久888优播| 不卡av一区二区三区| 黑人操中国人逼视频| 免费观看a级毛片全部| 国产欧美日韩综合在线一区二区| 国产激情久久老熟女| 搡老乐熟女国产| www日本在线高清视频| 好男人电影高清在线观看| 欧美黄色片欧美黄色片| 首页视频小说图片口味搜索| 天堂8中文在线网| 老司机深夜福利视频在线观看| 久久久精品区二区三区| 在线观看一区二区三区激情| 日本黄色日本黄色录像| 久久久久久久精品吃奶| 一进一出抽搐动态| www.熟女人妻精品国产| 99九九在线精品视频| 中文亚洲av片在线观看爽 | 国产欧美日韩一区二区三区在线| 久久精品熟女亚洲av麻豆精品| 五月开心婷婷网| av网站免费在线观看视频| 丝袜人妻中文字幕| 夜夜爽天天搞| 香蕉国产在线看| 嫁个100分男人电影在线观看| 婷婷丁香在线五月| 看免费av毛片| 久久国产精品影院| 精品视频人人做人人爽| 岛国毛片在线播放| 亚洲七黄色美女视频| 亚洲一卡2卡3卡4卡5卡精品中文| 法律面前人人平等表现在哪些方面| 亚洲午夜理论影院| 精品国产亚洲在线| 精品一区二区三卡| 亚洲av电影在线进入| 中文欧美无线码| av又黄又爽大尺度在线免费看| 中文字幕av电影在线播放| 建设人人有责人人尽责人人享有的| 午夜免费成人在线视频| 欧美日本中文国产一区发布| 久久精品国产亚洲av高清一级| 国产黄色免费在线视频| 香蕉国产在线看| 国产亚洲午夜精品一区二区久久| 国产日韩一区二区三区精品不卡| 亚洲欧美一区二区三区久久| 999精品在线视频| 人妻 亚洲 视频| 99精国产麻豆久久婷婷| 麻豆国产av国片精品| 男女高潮啪啪啪动态图| 国产视频一区二区在线看| 天天躁夜夜躁狠狠躁躁| 欧美黄色片欧美黄色片| 久久亚洲真实| 成年女人毛片免费观看观看9 | 亚洲熟女精品中文字幕| 男女高潮啪啪啪动态图| 91九色精品人成在线观看| 国产精品99久久99久久久不卡| 久久精品国产a三级三级三级| 国产在线免费精品| 2018国产大陆天天弄谢| 老司机靠b影院| 91成年电影在线观看| 亚洲国产欧美在线一区| 亚洲成人免费电影在线观看| 国产精品自产拍在线观看55亚洲 | 女性被躁到高潮视频| 日韩欧美一区二区三区在线观看 | 色婷婷av一区二区三区视频| 亚洲色图av天堂| 一本一本久久a久久精品综合妖精| 欧美日韩福利视频一区二区| 99国产精品一区二区三区| 国产主播在线观看一区二区| 国产成人av教育| 欧美 亚洲 国产 日韩一| 久久精品aⅴ一区二区三区四区| 日韩大码丰满熟妇| 黑人操中国人逼视频| avwww免费| 日韩大码丰满熟妇| 一本大道久久a久久精品| 欧美黑人精品巨大| 99久久国产精品久久久| 黄色成人免费大全| 性高湖久久久久久久久免费观看| 国产精品免费大片| 精品国产超薄肉色丝袜足j| 久久久国产精品麻豆| 老司机午夜福利在线观看视频 | 51午夜福利影视在线观看| 久久ye,这里只有精品| 亚洲国产成人一精品久久久| 久久精品国产亚洲av高清一级| 最近最新中文字幕大全电影3 | 国产精品久久久人人做人人爽| 69精品国产乱码久久久| 激情在线观看视频在线高清 | 在线观看人妻少妇| 国产精品麻豆人妻色哟哟久久| 欧美激情极品国产一区二区三区| 亚洲七黄色美女视频| 日韩有码中文字幕| 夫妻午夜视频| 一区二区三区乱码不卡18| 999精品在线视频| 国产亚洲精品第一综合不卡| 国产1区2区3区精品| 亚洲精品成人av观看孕妇| 午夜福利乱码中文字幕| 国产精品av久久久久免费| 日本一区二区免费在线视频| 成人手机av| 18禁黄网站禁片午夜丰满| 成人三级做爰电影| 一边摸一边抽搐一进一出视频| 国产又爽黄色视频| 肉色欧美久久久久久久蜜桃| 日韩有码中文字幕| 久久国产亚洲av麻豆专区| 韩国精品一区二区三区| 欧美在线黄色| 国产精品av久久久久免费| 精品久久蜜臀av无| 女性被躁到高潮视频| kizo精华| 免费看十八禁软件| 精品第一国产精品| 国产亚洲av高清不卡| 高清视频免费观看一区二区| 中国美女看黄片| 熟女少妇亚洲综合色aaa.| 亚洲一区中文字幕在线| 亚洲第一青青草原| 一级,二级,三级黄色视频| 免费人妻精品一区二区三区视频| 黑人猛操日本美女一级片| 下体分泌物呈黄色| 精品亚洲成国产av| 久久中文字幕人妻熟女| 1024视频免费在线观看| 亚洲综合色网址| 国产在视频线精品| 老熟妇仑乱视频hdxx| 久久亚洲精品不卡| 欧美精品一区二区大全| 亚洲精品国产一区二区精华液| 国产老妇伦熟女老妇高清| 一区福利在线观看| 国产精品 国内视频| 男女高潮啪啪啪动态图| 国产日韩欧美视频二区| 久久久欧美国产精品| 中文欧美无线码| 精品少妇久久久久久888优播| 国产一区二区三区视频了| 国产一区二区三区综合在线观看| 欧美日韩福利视频一区二区| 国产成人精品久久二区二区免费| 精品亚洲乱码少妇综合久久| 亚洲人成77777在线视频| 如日韩欧美国产精品一区二区三区| 女警被强在线播放| 国产亚洲精品一区二区www | 9色porny在线观看| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利一区二区在线看| 久久天躁狠狠躁夜夜2o2o| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av成人一区二区三| 97在线人人人人妻| 亚洲七黄色美女视频| www.熟女人妻精品国产| 激情视频va一区二区三区| 亚洲成人国产一区在线观看| 欧美乱妇无乱码| 黄色毛片三级朝国网站| 丁香六月天网| 日韩视频一区二区在线观看| 免费不卡黄色视频| 丁香六月欧美| 少妇粗大呻吟视频| 久久狼人影院| 日韩人妻精品一区2区三区| 国产欧美日韩综合在线一区二区| 亚洲中文日韩欧美视频| 精品国产一区二区三区四区第35| av超薄肉色丝袜交足视频| av一本久久久久| 不卡一级毛片| 视频区欧美日本亚洲| 桃花免费在线播放| 国产精品国产高清国产av | 黄片小视频在线播放| 欧美日韩福利视频一区二区| 99九九在线精品视频| 国产在线观看jvid| 国产精品偷伦视频观看了| 99国产精品99久久久久| 久久亚洲精品不卡| 操出白浆在线播放| 国产成人精品久久二区二区免费| 亚洲,欧美精品.| 日韩大片免费观看网站| 99国产极品粉嫩在线观看| 国产av又大| 欧美日韩精品网址| 69av精品久久久久久 | 中国美女看黄片| 欧美 亚洲 国产 日韩一| 亚洲av日韩精品久久久久久密| 色尼玛亚洲综合影院| 丝袜喷水一区| 久久这里只有精品19| 久久婷婷成人综合色麻豆| 亚洲国产精品一区二区三区在线| 欧美黄色淫秽网站| 成人精品一区二区免费| 日本五十路高清| 一夜夜www| 99久久精品国产亚洲精品| 欧美日韩国产mv在线观看视频| 精品国产一区二区三区四区第35| 电影成人av| 最新美女视频免费是黄的| 亚洲熟女毛片儿| 午夜福利,免费看| 午夜精品久久久久久毛片777| 丰满少妇做爰视频| 人妻 亚洲 视频| 亚洲天堂av无毛| 美国免费a级毛片| 国产不卡av网站在线观看| 久久香蕉激情| 老司机深夜福利视频在线观看| 国产av一区二区精品久久| 久久久国产成人免费| 高清欧美精品videossex| 日韩欧美一区二区三区在线观看 | 精品国产乱码久久久久久男人| 伦理电影免费视频| 亚洲伊人色综图| 18在线观看网站| 亚洲午夜理论影院| 国产免费av片在线观看野外av| 精品熟女少妇八av免费久了| 中文字幕制服av| 欧美成狂野欧美在线观看| 黑人巨大精品欧美一区二区mp4| 亚洲国产精品一区二区三区在线| 欧美人与性动交α欧美软件| 黄频高清免费视频| 夜夜骑夜夜射夜夜干| 天天添夜夜摸| 99久久99久久久精品蜜桃| 人妻 亚洲 视频| 欧美另类亚洲清纯唯美| 午夜福利欧美成人| 亚洲综合色网址| 亚洲欧美色中文字幕在线| 最新在线观看一区二区三区| 亚洲美女黄片视频| 成年女人毛片免费观看观看9 | 少妇被粗大的猛进出69影院| 精品熟女少妇八av免费久了| 夜夜骑夜夜射夜夜干| 在线观看人妻少妇| 男女下面插进去视频免费观看| 老司机午夜福利在线观看视频 | 免费黄频网站在线观看国产| 久久精品国产亚洲av高清一级| av一本久久久久| 老熟妇仑乱视频hdxx| 亚洲精品中文字幕一二三四区 | 在线亚洲精品国产二区图片欧美| 国产免费av片在线观看野外av| 亚洲五月婷婷丁香| 咕卡用的链子| 在线永久观看黄色视频| 精品国产乱码久久久久久小说| 老司机亚洲免费影院| 亚洲精品美女久久av网站| 久久久久视频综合| 老汉色∧v一级毛片| 天堂中文最新版在线下载| 久久久久久久久免费视频了| 99国产极品粉嫩在线观看| 国产成人影院久久av| 国产亚洲精品一区二区www | 亚洲中文av在线| 一级片'在线观看视频| 午夜精品国产一区二区电影| 精品福利永久在线观看| 脱女人内裤的视频| 午夜福利乱码中文字幕| 美女高潮喷水抽搐中文字幕| 婷婷丁香在线五月| 国产精品久久久久久精品古装| av有码第一页| 亚洲综合色网址| 香蕉久久夜色| 欧美黑人欧美精品刺激| 国产97色在线日韩免费| 超碰97精品在线观看| 夜夜爽天天搞| 久久精品国产综合久久久| 高清毛片免费观看视频网站 | 亚洲成人国产一区在线观看| 国产单亲对白刺激| 在线亚洲精品国产二区图片欧美| 色精品久久人妻99蜜桃| 色播在线永久视频| 在线永久观看黄色视频| 久久久久久亚洲精品国产蜜桃av| 久久久久久免费高清国产稀缺| 国产成人精品久久二区二区91| 在线观看免费高清a一片| 久久av网站| 美女扒开内裤让男人捅视频| 十八禁高潮呻吟视频| 欧美乱妇无乱码| 精品福利观看| 国产有黄有色有爽视频| 精品国产乱码久久久久久小说| 一区在线观看完整版| 欧美激情久久久久久爽电影 | 国产亚洲欧美在线一区二区| 国产亚洲午夜精品一区二区久久| 人人妻,人人澡人人爽秒播| 中国美女看黄片| 亚洲欧美一区二区三区久久| 精品少妇黑人巨大在线播放| 欧美成人免费av一区二区三区 | 亚洲午夜精品一区,二区,三区| 亚洲中文av在线| 91大片在线观看| 99re在线观看精品视频| 三级毛片av免费| 999久久久精品免费观看国产| 亚洲国产精品一区二区三区在线| 一区二区三区国产精品乱码| 美女高潮喷水抽搐中文字幕| 美女扒开内裤让男人捅视频| 国产日韩欧美亚洲二区| 男人操女人黄网站| 久久久久网色| 精品国产乱码久久久久久男人| 美女福利国产在线| 日韩欧美一区视频在线观看| bbb黄色大片| 人成视频在线观看免费观看| 日韩一卡2卡3卡4卡2021年| 亚洲三区欧美一区| 国产一区二区激情短视频| 亚洲欧美一区二区三区黑人| 精品欧美一区二区三区在线| 欧美日韩亚洲高清精品| 亚洲精品自拍成人| 欧美人与性动交α欧美软件| 成人亚洲精品一区在线观看| 青青草视频在线视频观看| 国产不卡一卡二| 亚洲成a人片在线一区二区| 欧美日韩亚洲高清精品| 亚洲av国产av综合av卡| 久久精品亚洲熟妇少妇任你| 午夜福利免费观看在线| 午夜免费鲁丝| 他把我摸到了高潮在线观看 | 麻豆成人av在线观看| 男女之事视频高清在线观看| 大片免费播放器 马上看| 精品少妇黑人巨大在线播放| 成人影院久久| 人妻 亚洲 视频| 欧美一级毛片孕妇| 美女午夜性视频免费| 伦理电影免费视频| 日韩 欧美 亚洲 中文字幕| 黑人操中国人逼视频| 人人妻人人爽人人添夜夜欢视频| 夜夜夜夜夜久久久久| 亚洲国产成人一精品久久久| 男人操女人黄网站| 黑人巨大精品欧美一区二区mp4| 最近最新中文字幕大全免费视频| 欧美日韩国产mv在线观看视频| 亚洲精品粉嫩美女一区| 搡老乐熟女国产| 亚洲精品在线美女| 男女之事视频高清在线观看| 丝袜美足系列| 老汉色av国产亚洲站长工具| 国产高清国产精品国产三级| 精品国产一区二区三区四区第35| 日韩大片免费观看网站| 黑人巨大精品欧美一区二区蜜桃| 欧美日韩成人在线一区二区| 在线观看66精品国产| 交换朋友夫妻互换小说| 不卡av一区二区三区| 欧美激情高清一区二区三区| 午夜福利影视在线免费观看| 国产av精品麻豆| 黑人操中国人逼视频| 露出奶头的视频| 色综合婷婷激情| 91av网站免费观看| 一本久久精品| 精品国产一区二区三区久久久樱花| 丝袜美足系列| 亚洲av日韩精品久久久久久密| 淫妇啪啪啪对白视频| 热99久久久久精品小说推荐| 多毛熟女@视频| 日韩欧美一区视频在线观看| 久久精品国产99精品国产亚洲性色 | 18禁黄网站禁片午夜丰满| 亚洲国产欧美日韩在线播放| 亚洲色图 男人天堂 中文字幕| 变态另类成人亚洲欧美熟女 | 夫妻午夜视频| 国产淫语在线视频| 777米奇影视久久| 国产日韩欧美在线精品| av又黄又爽大尺度在线免费看| tube8黄色片| 两人在一起打扑克的视频| 亚洲av成人不卡在线观看播放网| 最黄视频免费看| 999久久久国产精品视频| 在线观看免费午夜福利视频| 久久久久久久国产电影| 国产欧美日韩综合在线一区二区| 日韩有码中文字幕| 99国产极品粉嫩在线观看| 少妇的丰满在线观看| 搡老乐熟女国产| 一区福利在线观看| av一本久久久久| 热99久久久久精品小说推荐| 久久国产精品人妻蜜桃| 热99久久久久精品小说推荐| 天天添夜夜摸| 久久人妻av系列| 男女床上黄色一级片免费看| 国产高清国产精品国产三级| 国产不卡av网站在线观看| 国产日韩一区二区三区精品不卡| 午夜福利视频在线观看免费| 一本色道久久久久久精品综合| 在线观看人妻少妇| 国产国语露脸激情在线看| 999久久久国产精品视频| 三级毛片av免费| 韩国精品一区二区三区| 国产一区二区在线观看av| videos熟女内射| 久久久久国产一级毛片高清牌| 国产精品久久久av美女十八| 亚洲成国产人片在线观看| 日本欧美视频一区| 精品欧美一区二区三区在线| avwww免费| av一本久久久久| 欧美激情 高清一区二区三区| 色视频在线一区二区三区| 国产一区二区三区综合在线观看| 午夜免费鲁丝| av有码第一页| 天天躁日日躁夜夜躁夜夜| 天堂中文最新版在线下载| 成年人午夜在线观看视频| 午夜福利在线观看吧| 麻豆乱淫一区二区| 99香蕉大伊视频| 成年动漫av网址| 91精品国产国语对白视频| 亚洲黑人精品在线| 亚洲全国av大片| 亚洲精品国产色婷婷电影| 99精国产麻豆久久婷婷| 999久久久精品免费观看国产| 亚洲国产欧美在线一区| 黄色视频不卡| 欧美日韩精品网址| 性少妇av在线| 久久中文看片网| 超碰成人久久| 国产欧美日韩综合在线一区二区| 欧美变态另类bdsm刘玥| 精品国产国语对白av| 国产精品成人在线| 欧美在线黄色| 精品高清国产在线一区| 国产黄色免费在线视频| 欧美av亚洲av综合av国产av| 99精品在免费线老司机午夜| 亚洲国产欧美一区二区综合| 老熟女久久久| 9色porny在线观看| 久久久久久久精品吃奶| 日韩制服丝袜自拍偷拍| 丰满迷人的少妇在线观看| 亚洲精品中文字幕一二三四区 | 欧美精品一区二区免费开放| 日韩 欧美 亚洲 中文字幕| 99国产精品99久久久久| 超色免费av| 国产xxxxx性猛交| 黄色毛片三级朝国网站| 18在线观看网站| 亚洲第一av免费看| 国精品久久久久久国模美| 色精品久久人妻99蜜桃| 一个人免费看片子| 一级毛片精品| 高清毛片免费观看视频网站 | a级毛片黄视频| 视频区图区小说| 午夜福利,免费看| 下体分泌物呈黄色| 两人在一起打扑克的视频| 捣出白浆h1v1| 又紧又爽又黄一区二区| 亚洲第一欧美日韩一区二区三区 | 国产在线精品亚洲第一网站| 日本a在线网址| 国产深夜福利视频在线观看| 青青草视频在线视频观看| 啪啪无遮挡十八禁网站| 精品亚洲成a人片在线观看| 99精品在免费线老司机午夜| 这个男人来自地球电影免费观看| 亚洲国产毛片av蜜桃av| 亚洲精品av麻豆狂野| 亚洲av国产av综合av卡| 巨乳人妻的诱惑在线观看| 午夜日韩欧美国产| 成人黄色视频免费在线看| 另类亚洲欧美激情| 又黄又粗又硬又大视频|