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

    我國土壤中重金屬銅的生物配體模型的建立與應(yīng)用

    2014-09-27 07:32:13馬義兵王學(xué)東
    生態(tài)毒理學(xué)報 2014年4期
    關(guān)鍵詞:土壤溶液活度大麥

    李 波,馬義兵,王學(xué)東

    1. 遼寧省農(nóng)業(yè)科學(xué)院植物營養(yǎng)與環(huán)境資源研究所,沈陽 1101612. 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 1000813. 首都師范大學(xué)資源環(huán)境與旅游學(xué)院,北京 100048

    我國土壤中重金屬銅的生物配體模型的建立與應(yīng)用

    李 波1,馬義兵2,*,王學(xué)東3

    1. 遼寧省農(nóng)業(yè)科學(xué)院植物營養(yǎng)與環(huán)境資源研究所,沈陽 1101612. 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 1000813. 首都師范大學(xué)資源環(huán)境與旅游學(xué)院,北京 100048

    關(guān)于重金屬污染的風(fēng)險評價及其預(yù)測模型日益成為環(huán)境領(lǐng)域的研究熱點。近年來,一種用于預(yù)測和評價環(huán)境中重金屬生物毒性的機(jī)理性模型—生物配體模型(BLM)被廣泛應(yīng)用于水體及陸地生態(tài)系統(tǒng)。本研究以我國土壤的陸地生物配體模型(TBLM)建構(gòu)為目標(biāo),以土壤溶液系統(tǒng)為媒介,通過17種土壤上重金屬銅離子與大麥根長的相互作用關(guān)系,發(fā)現(xiàn)了土壤中Cu-TBLM的主要影響因素為Cu2+、CuOH+、Mg2+以及銅離子與大麥根系表面的專性結(jié)合能力?;谀P痛篼湼L預(yù)測值與實測值之間的相關(guān)關(guān)系,通過數(shù)學(xué)擬合功能求得TBLM中各參數(shù)值為logKCuBL= 4.87、logKCuOH+= 7.62、logKMgBL= 1.91、f50%= 0.103、β= 1.09。本研究所得到的TBLM模型能很好地預(yù)測我國土壤中銅對大麥根長的毒害程度,預(yù)測值與實測值的相關(guān)性達(dá)到了90%。本研究結(jié)果不僅可以為我們降低重金屬離子生物有效性提供有力的理論借鑒,更對我國土壤環(huán)境質(zhì)量保護(hù)和長期良性可持續(xù)發(fā)展具有重要意義。

    銅;大麥根長;陸地生物配體模型;毒性預(yù)測

    眾所周知,土壤中重金屬離子的全量不能很好地評價其生物有效性,真正作用于土壤生物體的重金屬離子只占全量的很小部分[1-3]。土壤溶液是金屬離子與生物體作用的媒介,水溶性金屬離子被認(rèn)為是作用于植物根系最直接、最有效的形態(tài)[4-5]。以銅為例,土壤溶液中總銅含量與植物生長效應(yīng)具有很好的相關(guān)性[4,6],但這種經(jīng)驗的劑量-效應(yīng)關(guān)系很難從機(jī)理層面解釋銅離子與土壤中共存生物體的相互作用關(guān)系。

    近幾年來,應(yīng)用生物配體模型(BLM)理論解釋溶液中金屬離子對植物、動物、微生物的毒性作用的報道逐漸增多[7-10]。該模型將生物受體位點作為生物配體,假設(shè)當(dāng)結(jié)合在具有生理活性的生物受體位點的重金屬達(dá)到一定量時,毒性就可能發(fā)生。模型考慮了影響生物毒性的溶液組成性質(zhì)以及具有生物有效性的金屬離子比例,在水質(zhì)、土壤毒性預(yù)測上取得了一定的成績[9-11]。與水體BLM相似,土壤中的生物配體模型(陸地生物配體模型,TBLM)把土壤溶液相作為金屬離子與生物體作用的媒介,一切作用的發(fā)生都以這個媒介為基礎(chǔ),然而,不同的是TBLM相對于水體BLM來說更加復(fù)雜,土壤作為土壤溶液成分的供給源引入了更多不可控制的因素。目前關(guān)于TBLM的研究還只處于探索階段[12-13]。TBLM模型以有效態(tài)金屬離子為研究出發(fā)點,任何影響金屬離子有效性的因子都應(yīng)該被模型所考慮。以銅離子為例,土壤溶液中銅離子有效性與自由銅離子活度、生物配體專性作用位點、以及其他可能具有生物有效性的銅離子形態(tài)有關(guān),同時,自由銅離子活度還受土壤溶液中共存競爭陽離子數(shù)量(如Ca2+、Mg2+等)、可溶性有機(jī)碳(DOC)濃度以及pH的影響。本研究以我國17種土壤中銅離子的TBLM建構(gòu)為研究目標(biāo),從土壤中有效性銅離子形態(tài)與植物毒性反應(yīng)關(guān)系入手,針對不同土壤類型中銅的有效性及影響因素,從機(jī)理角度解釋土壤中銅離子與植物根系的相互作用關(guān)系,探討TBLM在我國土壤上的預(yù)測能力及應(yīng)用性。

    1 材料與方法(Materids and methods)

    1.1 土壤溶液參數(shù)測量

    根據(jù)ISO 11269-1[16]大麥根伸長毒性試驗的標(biāo)準(zhǔn)方法,取風(fēng)干的污染土壤樣品每個處理360 g,裝入自制根長生長筒內(nèi)(內(nèi)徑3.5 cm,高12 cm),3次重復(fù),保持最大持水量(WHC)的60%培養(yǎng)7 d后,種植預(yù)發(fā)芽的大麥種子4粒到根長生長筒內(nèi),置于人工氣候箱內(nèi)生長,生長條件為白天14 h(22 °C),夜間10 h(18 °C),光照強(qiáng)度為24 000 lux·m-2,水分濕度為70%,保持整個試驗過程中土壤含水量為最大持水量的60%,5 d后收集大麥根,清洗后測定每個重復(fù)內(nèi)每株大麥的最長的根長,并求均值。

    1.3 TBLM的建立

    1.3.1 土壤溶液中銅離子活度測量

    使用Orion 94-29銅離子選擇電極(Thermo Electron Corporation, USA)測量土壤溶液中銅離子活度[17]。由于土壤樣本較大以及污染處理較多,本研究僅選取17個土壤樣品的兩個污染濃度,即約為大麥根伸長10%(EC10)及50%(EC50)抑制時的濃度水平[18]。電極標(biāo)準(zhǔn)溶液使用各5 mL 1 mol·L-1Cu(NO3)2、0.01 mol·L-1IDA、0.025 mol·L-1KHC8H4O4、0.1 mol·L-1KNO3及不同體積的0.02 mol·L-1NaOH(1~10 mL)配置[13]。根據(jù)電極標(biāo)準(zhǔn)溶液測定結(jié)果繪制電極電位mV對pCu(自由離子活度的負(fù)對數(shù))標(biāo)準(zhǔn)曲線圖,標(biāo)準(zhǔn)溶液的pCu結(jié)果使用WHAM 6[19]模型計算。測量土壤溶液中銅離子的電極電位,根據(jù)標(biāo)準(zhǔn)曲線圖查出相應(yīng)的pCu值。

    1.3.2 土壤溶液中銅離子活度的預(yù)測

    表1 17 種土壤的基本理化性狀Table 1 Selected physicochemical properties of soil used in this experiment

    表2 17個土壤樣本的土壤溶液參數(shù)Table 2 Soil solution parameters measured for Cu added to 17 soils by extracting soil solution

    1.3.3 TBLM的建立

    根據(jù)Wang等[20]研究結(jié)果,即除二價銅離子(Cu2+)之外,羥基銅離子(CuOH+)在高pH條件下具有很強(qiáng)的毒性,考慮到我國土壤背景中大部分土壤的pH高于7,因此在TBLM建模時Cu2+與CuOH+被同時考慮為生物有效性的銅離子形態(tài)。此外,Mg2+可以與Cu2+競爭大麥根表面的結(jié)合位點,因此被作為Cu2+毒性消減作用的唯一離子[20]。最終輸入到TBLM模型的變量鎖定為{Cu2+}(離子活度形式)、{CuOH+}、{Mg2+}以及大麥根長相對長度RE,TBLM方程為[13]:

    2 結(jié)果與討論(Results and discussion)

    2.1 WHAM 6 模型校正結(jié)果

    根據(jù)測量的與模型預(yù)測的pCu的一致程度,采用最大相關(guān)系數(shù)及最小殘差平方和法校正WHAM 6化學(xué)模型中FA默認(rèn)參數(shù)。研究發(fā)現(xiàn)我國農(nóng)業(yè)土壤的DOC對銅的絡(luò)合能力與WHAM 6中默認(rèn)的絡(luò)合能力相差很大,即我國農(nóng)業(yè)土壤的DOC遠(yuǎn)小于模型默認(rèn)能力。通常來說,WHAM 6中65%的FA被認(rèn)為是活性的FA[21-24],而對于我國土壤來說,只有當(dāng)銅與FA的結(jié)合常數(shù)logKCuFA從默認(rèn)的2.1(圖1a)下降到1.5(圖1b)時,測量與預(yù)測的pCu才最吻合,此時兩者的相關(guān)系數(shù)R2從0.51增加到0.74,方程的斜率值從0.931增加到1.005(圖1)。因此,輸入到WHAM 6中DOC參數(shù)為FA = 2 × DOC × 0.65,并logKCuFA=1.5。

    2.2 TBLM模型建立

    輸入表2中136個土壤溶液變量組到校正的WHAM 6模型,得到的銅離子和其他陽離子的活度見表3。從表3發(fā)現(xiàn),在pH<7.3的土壤里,銅離子活度占總?cè)芙庑糟~的40%以內(nèi),且與土壤pH呈負(fù)相關(guān);pH>7.3的土壤里,銅離子活度僅占總?cè)芙庑糟~的15%以內(nèi),大部分銅離子與DOM結(jié)合起來,失去活性。此外,總體上{Cu2+}顯著高于{CuOH+},{Ca2+}離子活度高于{Mg2+}。輸入表3中變量到DPS 9.0數(shù)據(jù)表中,根據(jù)大麥根長的預(yù)測值與測量值之間的相關(guān)關(guān)系,擬合得到各參數(shù)值為logKCuBL= 4.87、logKCuOHBL=7.62、logKMgBL=1.91、f50%= 0.103、β= 1.09。預(yù)測和測量的大麥根長的相關(guān)關(guān)系如圖2a所示,其中殘差平方和RMSE = 11.4,相關(guān)系數(shù)R2= 0.90??傮w看來大麥根長的預(yù)測結(jié)果很好,基本所有的預(yù)測值都落在實際測量值2倍的范圍內(nèi)(圖2b),意味著TBLM能夠很好地預(yù)測我國土壤中銅對大麥根長生長的毒害程度。

    圖1 銅離子選擇電極測量的pCu與WHAM 6模型預(yù)測的pCu的相關(guān)性Fig. 1 Comparison between measured and predicted pCu. The predictions were made with WHAM 6 using default parameter values and 65% AFA (a), optimized parameter values and 65% AFA (b).

    表3 土壤溶液離子活度預(yù)測值及大麥根長值Table 3 Relative cation activity predicated using soil solution parameter for 17 Cu added soils

    2.3 TBLM驗證

    為了進(jìn)一步確定本研究結(jié)果的準(zhǔn)確性,我們應(yīng)用王學(xué)東等[12]大麥根長水培試驗得到的參數(shù)值即logKCuBL= 6.57、logKCuOHBL= 7.03、logKMgBL= 3.00代入到以上TBLM方程中,通過DPS軟件數(shù)據(jù)擬合得到f50%= 0.047、β= 0.625,預(yù)測值與實測值的相關(guān)系數(shù)R2= 0.51,RMSE = 22.47,說明模型預(yù)測能力很差(表4)。此外,我們也使用Thakali等[9]中得到的參數(shù)f50%= 0.05、β= 0.96以及王學(xué)東等(2008)得到的結(jié)合常數(shù)logKCuBL= 6.57、logKCuOHBL= 7.03、logKMgBL= 3.00代入到TBLM方程中,結(jié)果發(fā)現(xiàn)模型的預(yù)測能力也很差,即R2= 0.56,RMSE = 23.71。盡管本研究得到的TBLM模型不盡完美,仍有一些數(shù)值偏離實測值兩倍的范圍區(qū)間,但相對于Thakali等[9]基于歐洲非碳酸鈣土壤的TBLM模型(RMSE =15.4;R2= 0.85)卻有很大的提高(表4)。此外,Thakali等[9]選取的土壤樣本pH較低,有機(jī)質(zhì)含量較高,而我國碳酸鈣土壤比例較大,有機(jī)質(zhì)含量較低,基于我國土壤背景得到的TBLM更能針對我國土壤的分布特點,從機(jī)理上解釋土壤溶液系統(tǒng)中銅與植物根系作用的方式。

    2.4 TBLM建模的影響因素

    由于DOC對銅有很強(qiáng)的絡(luò)合能力,在污染較輕的土壤上,幾乎90%以上的銅以DOC絡(luò)合形式存在,因此,DOC的濃度及結(jié)構(gòu)差異對銅的生物有效性影響很大。通常來說,酸性土壤溶液中的DOC比中性與堿性土壤絡(luò)合能力低,主要因為酸性土壤中H+含量高,H+可以占據(jù)DOC與Cu2+的結(jié)合位點,降低DOC對銅的絡(luò)合能力[9]。本研究中沒有引入不同pH土壤的DOC差異,主要是因為在模型校正時沒有發(fā)現(xiàn)極異常數(shù)據(jù)值。其次,DOC因來源不同,結(jié)構(gòu)差異很大,一般來源于農(nóng)業(yè)土壤的DOC要比來源于有機(jī)肥、秸稈腐熟物等DOC芳香度低,也就是說來源于農(nóng)業(yè)土壤的DOC含有較多的羧基、較少的芳香基[25-26],從而導(dǎo)致前者的絡(luò)合能力相對較低,因此,追溯農(nóng)業(yè)土壤的外源有機(jī)物來源及添加量對于模型準(zhǔn)確建構(gòu)具有一定的作用。再次,不同土壤溶液的提取方法產(chǎn)生的效果不同。MacDonald等[27]提出使用滲透方法提取土壤溶液,能夠較好地模擬了田間土壤的化學(xué)性質(zhì),而采用離心、壓力、抽氣等方法提取的土壤溶解性離子的濃度要大于滲透方法測得的離子濃度。因此,正確的土壤溶液提取方法對TBLM的準(zhǔn)確建立起到非常重要的作用。

    圖2 大麥根長的測量值和預(yù)測值的相關(guān)性。其中a為百分?jǐn)?shù)、b為對數(shù)轉(zhuǎn)換值Fig. 2 The relationship between the measured and predicted barley relative root lengths, and a (%) and b (log-transformed) represent Cu toxicity results

    表4 不同研究得到的大麥根系與銅離子結(jié)合的生物配體模型參數(shù)Table 4 The BLM parameters based on barley root and Cu ion binding from different studies

    此外,不同的生物配體對TBLM的建構(gòu)影響很大。Thakali等[9-10]基于pH ≤ 7.0的歐洲土壤的植物生長和微生物活性抑制試驗建立了適合歐洲非堿性土壤的TBLM,研究發(fā)現(xiàn)對于土壤硝化細(xì)菌來說降低銅害的競爭離子為H+與Mg2+,而對于大麥根長與西紅柿莖葉生長來說競爭離子僅為H+;然而,Lock等[28]基于非堿性溶液培養(yǎng)試驗的研究發(fā)現(xiàn)降低大麥根長銅害的競爭離子為H+與Mg2+。王學(xué)東等[12,20]通過大麥根長水培試驗,進(jìn)一步驗證了Mg2+的競爭作用;并通過pH范圍試驗,發(fā)現(xiàn)了除Cu2+之外,CuOH+的毒害作用,進(jìn)一步解釋了低pH條件下H+對Cu2+的競爭作用[9,28]只是一種假象,而實際上是由于pH的升高,生成了CuOH+這種比Cu2+毒性更強(qiáng)的形態(tài)。由于本研究中大部分土壤pH超過7,因此,選用王學(xué)東的研究結(jié)果構(gòu)建我國TBLM。值得注意的是,不同生物配體與銅離子的結(jié)合能力(以logK表示)有所不同,對于西紅柿莖葉生長來說,雖然競爭離子與大麥根長相同,但絡(luò)合常數(shù)logK卻有顯著差異[9-10]。Cheng等[29]研究了銅對生菜的毒性,估算生菜根系與銅的結(jié)合常數(shù)(logK)變化在11.15~12.16之間;Antunes等[30]通過小麥培養(yǎng)試驗發(fā)現(xiàn)了銅濃度不同時小麥根和銅的結(jié)合常數(shù)(logK)不同,變化范圍在4.07~9.06之間。本研究得到的結(jié)合常數(shù)遠(yuǎn)小于生菜根與銅的結(jié)合常數(shù),而與小麥根結(jié)合常數(shù)變化區(qū)間較為一致,可能由于大麥與小麥植物科屬相近所致。

    總之,BLM是一個機(jī)理性的模型,能夠從機(jī)理角度摸清土壤生物體與金屬離子的相互作用關(guān)系,是一個能夠替代生物毒性試驗來預(yù)測金屬毒性的有用工具。盡管模型現(xiàn)在存在著局限性,但隨著對痕量金屬吸收過程中的化學(xué)、生理學(xué)和生物學(xué)過程的深入理解,這些問題將會逐步得到解決,未來BLM既面臨挑戰(zhàn)又有著廣闊的應(yīng)用和發(fā)展前景。

    [1] Sauvé S, McBride M B, Norvell W A, et al. Copper solubility and speciation of in situ contaminated soils: Effects of copper level, pH and organic matter [J]. Water, Air and Soil Pollution, 1997, 100(1-2): 133-149

    [2] McLaughlin M J. Bioavailability of metals to terrestrial plants. In Allen HE (Ed.). Bioavailability of Metals in Terrestrial Ecosystems: Influence of Partitioning for Bioavailability to Invertebrates, Microbes and Plants [M]. Pensacola, FL: SETAC Press, 2000: 39-67

    [3] Nolan A L, Lombi E, McLaughlin M J. Metal bioaccumulation and toxicity in soils: Why bother with speciation? [J]. Australian Journal of Chemistry, 2003, 56(3): 77-91

    [4] Zhao F J, Rooney C P, Zhang H, et al. Comparison of soil solution speciation and diffusive gradients in thin-films measurement as an indicator of copper bioavailability to plants [J]. Environmental Toxicology and Chemistry, 2006, 25(3): 733-742

    [5] Song J, Zhao F J, Luo Y M, et al. Copper uptake by Elsholtzia splendens and Silene vulgaris and assessment of copper phytoavailability in contaminated soils [J]. Environmental Pollution, 2004, 128(3): 307-315

    [6] Sauvé S, Hendershot W, 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

    [7] Niyogi S, Wood C M. Biotic ligand model, a flexible tool for developing site-specific water quality guidelines for metals [J]. Environmental Science and Technology, 2004, 38(23): 6177-6192

    [8] Slaveykova V I, Wilkinson K J. Predicting the bioavailability of metal complexes: Critical review of the biotic ligand model [J]. Environmental Chemistry, 2005, 2(1): 9-24

    [9] Thakali S, Allen H E, Di Toro D M, et al. A terrestrial biotic ligand model. 1. Development and application to Cu and Ni toxicity to barley root elongation in soils [J]. Environmental Science and Technology, 2006, 40(22): 7085-7093

    [10] Thakali S, Allen H E, Di Toro D M, et al. A terrestrial biotic ligand model. 2. Application to Ni and Cu toxicities to plants, invertebrates, and microbes in soil [J]. Environmental Science and Technology, 2006, 40(22): 7094-7100

    [11] 黃圣彪. 水環(huán)境中銅形態(tài)與其生物有效性/毒性關(guān)系及其預(yù)測模型研究[D]. 北京, 中國科學(xué)院生態(tài)環(huán)境研究中心, 2003

    Huang S B. Measurement and prediction in the toxic/bioavailable concentration of copper in natural water [D]. Beijing, Research Center for Eco-Environmental Sciences, 2003 (in Chinese)

    [12] 王學(xué)東, 馬義兵, 華珞, 張璇. 銅對大麥的急性毒性預(yù)測模型-生物配體模型研究[J]. 環(huán)境科學(xué)學(xué)報, 2008, 28(8): 1704-1712

    Wang X D, Ma Y B, Hua L, et al. Development of Biotic Ligand Mod el (BLM) predicting copper acute toxicity to barley (Hordeum vulgare) [J]. Acta Scientiae Circumstantiae, 28(8): 1704-1712 (in Chinese)

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

    Li B. The Phytotoxicity of added copper and nickel to soils and predictive models [D]. Beijing, Chinese Academy of Agricultural Sciences, 2010 (in Chinese)

    [14] Thibault D H, Sheppard M I. A disposable system for soil pore-water extraction by centrifugation [J]. Communications in Soil Science and Plant Analysis, 1992, 23(13): 1629-1641

    [15] Zarcinas B A, McLaughlin M J, Smart M K. The effect of acid digestion technique on the performance of nebulisation systems used in inductively coupled plasma spectrometry [J]. Communications in Soil Science and Plant Analysis, 1996, 27(5-8): 1331-1354

    [16] International Organization for Standardization. Soil quality-Determination of the effects of pollutants on soil flora-Part 1: Method for the measurement of inhibition of root growth [S]. Geneva, Switzerland. 1993

    [17] Sauvé S, McBride M B. Hendershot W H. Ion-selective electrode measurements of copper (II) activity in contaminated soils [J]. Archives Environment Contamination Toxicology, 1995, 29(3): 373-379

    [18] Li B, Ma Y B, McLaughlin M J, et al. Influences of soil properties and leaching on copper toxicity to barley root elongation [J]. Environmental Toxicology and Chemistry, 2010, 29(4): 835-842

    [19] Tipping E. Humic Ion Binding Model VI: An improved description of the interactions of protons and metal ions with humic substances [J]. Aquatic Geochemistry,1998, 4(1): 3-48

    [20] Wang X D, Ma Y B, Hua L, et al. Identification of CuOH+toxicity to barley root elongation in solution culture [J]. Environmental Toxicology and Chemistry 2009, 28(3): 662-667

    [21] Weng L P, Temminghoff E J M, Lofts S, et al. Complexation with dissolved organic matter and solubility control of heavy metal in a sandy soil [J]. Environmental Science and Technology, 2002, 36(22): 4804-4810

    [22] Van Laer L, Smolders E, Degryse F, et al. Speciation of nickel in surface waters measured with the Donnan membrane technique [J]. Analytica Chimica Acta, 2006, 578(2): 195-202

    [23] Oorts K, Bronckaers H, Smolders E. Discrepancy of the microbial response to elevated Cu between freshly spiked and long-term contaminated soils [J]. Environmental Toxicology and Chemistry, 2006, 25(3): 845-853

    [24] Oorts K, Ghesquiere U, Smolders 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.

    [25] Weishaar J L, Aiken G R, Bergamaschi B A, et al. Evaluation of specific ultraviolet absorbance as an indicator of the chemical composition and reactivity of dissolved organic carbon [J]. Environmental Science and Technology, 2003, 37(20): 4702-4708

    [26] Grybos M, Davranche M, Gruau G, et al. Is trace metal release in wetland soils controlled by organic matter mobility or Fe-oxyhydroxides reduction? [J]. Journal of Colloid and Interface Science, 2007, 314(2): 490-501

    [27] MacDonald J D, Belanger N, Hendershot W H. Column leaching using dry soil to estimate solid - solution partitioning observed in zero - tension lysimeters. 2 Trace metals [J]. Soil and Sediment Contamination, 2004, 13(4): 375-390

    [28] Lock K, Criel P, De Schamphelaere K A C, et al. Influence of calcium, magnesium, sodium, potassium and pH on copper toxicity to barley (Hordeum vulgare) [J]. Ecotoxicology and Environmental Safety, 2007, 68: 299-304

    [29] Cheng T, Allen H E. Prediction of uptake copper from solution by lettuce (Lactuca sativa romance) [J]. Environmental Toxicology and Chemistry, 2000, 20(11): 2544-2551

    [30] Antunes P M C. Beyond the free-ion activity model: The biotic ligand model for estimating copper uptake by durum wheat [D]. Guelph, University of Guelph, 2004

    DevelopmentandApplicationofBioticLigandModelforCopperinChineseSoils

    Li Bo1,Ma Yibing2,*,Wang Xuedong3

    1. Institute of Plant Nutrition and Environmental Resources, Liaoning Academy of Agricultural Sciences, Shenyang 110161, China2. Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing 100081, China3. The Key Lab of Resource Environment and GIS, College of Resource Environment & Tourism, Capital Normal University, Beijing 100048, China

    15 May 2014accepted30 July 2014

    The terrestrial biotic ligand model (TBLM) is a newly developed method for evaluating bioavailability of heavy metal ion in soils. Seventeen soil samples covering China were selected in this study by adding a series of Cu concentrations. The soil solutions were extracted according to the method of maximum water holding capacity. The cation and anion concentrations as well as dissolved organic carbon (DOC) in soil solutions were measured for the predictions of Cu (Cu2+and CuOH+especially at high-pH conditions) and Mg ion activities. The WHAM 6 model was used to predict free Cu2+, CuOH+and Mg2+activities using 65% DOM as active fluvic acid (FA) and adjusted binding constant logKCuFA= 1.5. The predicted ion activities and measured barley relative root elongation as variables were input into TBLM equation, and binding constants and parameters were fitted out as logKCuBL= 4.87, logKCuOH+= 7.62, logKβ= 1.09 by DPS statistic software according the best agreement between measured and predicted barley relative root lengths. The results showed that predicted and measured barley root elongation can match quite well with an R2= 0.90, RMSE = 11.4 and almost all predicted values lied within 2-fold of measured values. In order to ensure accuracy of TBLM in our study, some parameters from relative studies were compared with the values in this study, which showed that our TBLM was better than others and more suitable for Cu prediction of bioavailability in Chinese soils because a series of high-pH soils samples were included. The study will help a lot on formulating guidance of soil environment risk assessments not only for China but for the other areas around world.

    : copper; barley; terrestrial biotic ligand model; toxicity prediction

    2014-631-15錄用日期:2014-07-30

    1673-5897(2014)4-632-08

    : X171.5

    : A

    馬義兵(1957—),男,博士,研究員,主要研究方向為土壤重金屬環(huán)境化學(xué)、風(fēng)險評價與治理,發(fā)表學(xué)術(shù)論文160余篇。

    國家自然科學(xué)基金項目(No.20677077; No.40620120436)和國際銅業(yè)協(xié)會,澳大利亞力拓礦業(yè)集團(tuán)資助

    李波(1979-), 女, 博士, 主要從事土壤重金屬污染毒理與修復(fù)研究. E-mail: libocaas@163.com

    *通訊作者(Corresponding author),E-mail: ybma@caas.ac.cn

    10.7524/AJE.1673-5897.20140510006

    李 波,馬義兵,王學(xué)東. 我國土壤中重金屬銅的生物配體模型的建立與應(yīng)用[J]. 生態(tài)毒理學(xué)報, 2014, 9(4): 632-639

    Li B, Ma Y B, Wang X D. Development and application of biotic ligand model for copper in Chinese soils [J]. Asian Journal of Ecotoxicology, 2014, 9(4): 632-639 (in Chinese)

    猜你喜歡
    土壤溶液活度大麥
    我的大麥哥哥
    氮添加對亞熱帶常綠闊葉林土壤溶液化學(xué)特性的影響
    大麥蟲對聚苯乙烯塑料的生物降解和礦化作用
    Global interest in Chinese baijiu
    臨桂區(qū)土壤溶液的水化學(xué)特征及其溶蝕能力淺析
    CaO-SiO2-FeO-P2O5-Al2O3脫磷渣系中組元活度的計算
    大麥若葉青汁
    華人時刊(2016年19期)2016-04-05 07:56:10
    核電廠惰性氣體排放活度濃度的估算
    土壤增溫對杉木幼林不同深度土壤溶液N03-濃度的影響
    土壤增溫對杉木幼林不同深度土壤溶液DOM的影響
    日日摸夜夜添夜夜爱| 韩国精品一区二区三区| 欧美av亚洲av综合av国产av | 日日摸夜夜添夜夜爱| 成人国产av品久久久| 99久久人妻综合| videossex国产| 美女xxoo啪啪120秒动态图| 色吧在线观看| 国产成人免费观看mmmm| 制服诱惑二区| 最新的欧美精品一区二区| 国产成人精品无人区| 熟女电影av网| 99精国产麻豆久久婷婷| av有码第一页| 日韩熟女老妇一区二区性免费视频| 久久狼人影院| 青春草视频在线免费观看| 亚洲少妇的诱惑av| 精品亚洲成国产av| 在线观看一区二区三区激情| 999精品在线视频| 色视频在线一区二区三区| 极品人妻少妇av视频| 97人妻天天添夜夜摸| 国产精品久久久av美女十八| 国产一区亚洲一区在线观看| 涩涩av久久男人的天堂| 免费看av在线观看网站| 美女国产高潮福利片在线看| 黄色毛片三级朝国网站| 国产男人的电影天堂91| 狠狠婷婷综合久久久久久88av| 欧美最新免费一区二区三区| 精品国产露脸久久av麻豆| 大香蕉久久网| 飞空精品影院首页| 亚洲 欧美一区二区三区| 国产精品 欧美亚洲| 久久99精品国语久久久| 亚洲人成77777在线视频| 国产白丝娇喘喷水9色精品| 日日啪夜夜爽| www.自偷自拍.com| 精品亚洲乱码少妇综合久久| 国产精品香港三级国产av潘金莲 | 不卡视频在线观看欧美| 国产精品蜜桃在线观看| 在线观看美女被高潮喷水网站| 欧美少妇被猛烈插入视频| 欧美日本中文国产一区发布| 韩国av在线不卡| 亚洲欧洲国产日韩| 亚洲人成网站在线观看播放| 亚洲欧美成人精品一区二区| 久久精品人人爽人人爽视色| 交换朋友夫妻互换小说| 男的添女的下面高潮视频| 亚洲av日韩在线播放| 久久女婷五月综合色啪小说| 亚洲精品aⅴ在线观看| 一区在线观看完整版| 精品午夜福利在线看| 日韩不卡一区二区三区视频在线| 日日撸夜夜添| 在线观看一区二区三区激情| 青春草亚洲视频在线观看| 各种免费的搞黄视频| 欧美精品av麻豆av| 咕卡用的链子| 天堂俺去俺来也www色官网| 久久人人97超碰香蕉20202| 亚洲五月色婷婷综合| 国产在线一区二区三区精| 男男h啪啪无遮挡| 国产激情久久老熟女| 婷婷色综合大香蕉| 欧美xxⅹ黑人| 国产一区有黄有色的免费视频| 国产精品久久久久久精品电影小说| 国产综合精华液| 久久综合国产亚洲精品| 免费在线观看黄色视频的| 亚洲av成人精品一二三区| 亚洲欧洲精品一区二区精品久久久 | 不卡视频在线观看欧美| 亚洲av国产av综合av卡| 久久 成人 亚洲| 自线自在国产av| 国产免费现黄频在线看| 视频区图区小说| 99九九在线精品视频| 精品亚洲成a人片在线观看| 一级毛片黄色毛片免费观看视频| 亚洲国产av影院在线观看| 精品国产国语对白av| 久久精品国产亚洲av天美| 精品少妇黑人巨大在线播放| 久久久久久伊人网av| 伦理电影免费视频| 我要看黄色一级片免费的| 999久久久国产精品视频| 国产深夜福利视频在线观看| 亚洲第一青青草原| 边亲边吃奶的免费视频| 亚洲av中文av极速乱| 亚洲国产欧美在线一区| 亚洲国产欧美网| 国产一级毛片在线| 久久精品国产亚洲av天美| 欧美av亚洲av综合av国产av | 欧美黄色片欧美黄色片| 欧美最新免费一区二区三区| 母亲3免费完整高清在线观看 | 国产白丝娇喘喷水9色精品| 蜜桃国产av成人99| 中文字幕精品免费在线观看视频| 毛片一级片免费看久久久久| 国产精品免费大片| 久久这里只有精品19| 亚洲成国产人片在线观看| 久久99蜜桃精品久久| 韩国高清视频一区二区三区| 婷婷色麻豆天堂久久| 黄片无遮挡物在线观看| 亚洲色图 男人天堂 中文字幕| 国产黄频视频在线观看| 亚洲成国产人片在线观看| 精品第一国产精品| 最近的中文字幕免费完整| 天天躁狠狠躁夜夜躁狠狠躁| 成人手机av| 少妇猛男粗大的猛烈进出视频| 乱人伦中国视频| 中文字幕人妻丝袜一区二区 | 十八禁高潮呻吟视频| 看免费av毛片| 国产麻豆69| 国产有黄有色有爽视频| 欧美变态另类bdsm刘玥| 久久国产亚洲av麻豆专区| 午夜免费鲁丝| 国精品久久久久久国模美| 久久青草综合色| 成人国产av品久久久| 天天影视国产精品| 男女午夜视频在线观看| 高清视频免费观看一区二区| 久久人妻熟女aⅴ| 制服丝袜香蕉在线| 69精品国产乱码久久久| 中文字幕最新亚洲高清| 欧美精品亚洲一区二区| 两个人免费观看高清视频| 欧美精品一区二区大全| 精品国产一区二区久久| 人人妻人人添人人爽欧美一区卜| 午夜91福利影院| 99久久综合免费| 久久久久网色| 美国免费a级毛片| 嫩草影院入口| 9191精品国产免费久久| 亚洲国产最新在线播放| 午夜日韩欧美国产| 亚洲 欧美一区二区三区| 欧美人与性动交α欧美软件| 侵犯人妻中文字幕一二三四区| 女性被躁到高潮视频| 日韩欧美精品免费久久| 亚洲成人手机| videossex国产| 国产免费一区二区三区四区乱码| 国产精品免费视频内射| 欧美激情 高清一区二区三区| 只有这里有精品99| 国产精品一区二区在线不卡| 国产 精品1| 国产片特级美女逼逼视频| 国产男人的电影天堂91| 国产精品三级大全| 婷婷成人精品国产| 欧美精品一区二区大全| 黄色怎么调成土黄色| 精品亚洲乱码少妇综合久久| 91午夜精品亚洲一区二区三区| 丰满迷人的少妇在线观看| 自线自在国产av| 一级片免费观看大全| 天美传媒精品一区二区| 精品视频人人做人人爽| 亚洲av男天堂| 欧美成人午夜精品| 看十八女毛片水多多多| 久久久亚洲精品成人影院| 波野结衣二区三区在线| 人妻一区二区av| 性高湖久久久久久久久免费观看| 国产欧美亚洲国产| 日韩精品有码人妻一区| 成人漫画全彩无遮挡| xxx大片免费视频| 街头女战士在线观看网站| 久久99一区二区三区| 街头女战士在线观看网站| 国产精品成人在线| 91国产中文字幕| 欧美国产精品一级二级三级| 乱人伦中国视频| 久久99一区二区三区| 色视频在线一区二区三区| 久久久a久久爽久久v久久| 中文天堂在线官网| 午夜老司机福利剧场| 久久婷婷青草| 极品人妻少妇av视频| 亚洲精品成人av观看孕妇| 考比视频在线观看| 国产精品偷伦视频观看了| 中文字幕av电影在线播放| 男女国产视频网站| 亚洲av男天堂| av免费在线看不卡| 国产精品欧美亚洲77777| 天天躁日日躁夜夜躁夜夜| 亚洲精品av麻豆狂野| av在线播放精品| 久久久久久久久免费视频了| 免费观看性生交大片5| 国产老妇伦熟女老妇高清| 欧美精品高潮呻吟av久久| 久久精品亚洲av国产电影网| 欧美国产精品va在线观看不卡| 深夜精品福利| 国产人伦9x9x在线观看 | 十八禁网站网址无遮挡| 国产男人的电影天堂91| 超碰97精品在线观看| 你懂的网址亚洲精品在线观看| 国产高清国产精品国产三级| 男女无遮挡免费网站观看| 久久精品久久久久久噜噜老黄| 久久久久网色| 免费高清在线观看日韩| 日韩成人av中文字幕在线观看| 亚洲精品中文字幕在线视频| 亚洲少妇的诱惑av| 天天影视国产精品| 国产成人一区二区在线| 色94色欧美一区二区| 久久久久精品久久久久真实原创| 久久99热这里只频精品6学生| 狂野欧美激情性bbbbbb| 亚洲第一青青草原| 肉色欧美久久久久久久蜜桃| 成人黄色视频免费在线看| 久久久久人妻精品一区果冻| 亚洲成人一二三区av| 捣出白浆h1v1| 99精国产麻豆久久婷婷| 最新中文字幕久久久久| 国产精品免费大片| 精品一区在线观看国产| 亚洲欧美一区二区三区久久| 国产一级毛片在线| 美女脱内裤让男人舔精品视频| 青春草亚洲视频在线观看| 新久久久久国产一级毛片| 免费高清在线观看视频在线观看| 国产精品熟女久久久久浪| 岛国毛片在线播放| 欧美日韩亚洲国产一区二区在线观看 | 国产福利在线免费观看视频| 又黄又粗又硬又大视频| 国产高清不卡午夜福利| 有码 亚洲区| 久久久久网色| 不卡av一区二区三区| 2021少妇久久久久久久久久久| 天堂俺去俺来也www色官网| 亚洲视频免费观看视频| av在线观看视频网站免费| 国产av码专区亚洲av| 成人毛片60女人毛片免费| 成人亚洲欧美一区二区av| 亚洲一区二区三区欧美精品| 26uuu在线亚洲综合色| 在线观看国产h片| 一级毛片我不卡| 亚洲av免费高清在线观看| 18禁观看日本| 男女免费视频国产| 精品久久久精品久久久| 看十八女毛片水多多多| 久久综合国产亚洲精品| 美女中出高潮动态图| av免费在线看不卡| 熟妇人妻不卡中文字幕| 日韩精品有码人妻一区| 2021少妇久久久久久久久久久| 极品少妇高潮喷水抽搐| 80岁老熟妇乱子伦牲交| 久久久久久久久久久免费av| 国产精品国产三级专区第一集| 9191精品国产免费久久| 女性生殖器流出的白浆| 亚洲精品av麻豆狂野| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | www日本在线高清视频| 国产精品.久久久| 亚洲欧洲精品一区二区精品久久久 | 三级国产精品片| 男女下面插进去视频免费观看| 嫩草影院入口| 精品少妇黑人巨大在线播放| 日韩成人av中文字幕在线观看| 午夜免费观看性视频| 精品久久久久久电影网| 国产深夜福利视频在线观看| 久久精品久久久久久噜噜老黄| 人妻系列 视频| 国产精品av久久久久免费| 青青草视频在线视频观看| 日韩 亚洲 欧美在线| 丝袜脚勾引网站| 久久久久久久大尺度免费视频| 久久国产精品男人的天堂亚洲| 久久久久精品人妻al黑| 高清黄色对白视频在线免费看| 26uuu在线亚洲综合色| 中文精品一卡2卡3卡4更新| 国产精品女同一区二区软件| 中文字幕制服av| 亚洲一区二区三区欧美精品| 久久人人爽人人片av| 亚洲国产精品一区二区三区在线| 高清欧美精品videossex| 亚洲成人av在线免费| 国产日韩一区二区三区精品不卡| 精品人妻在线不人妻| 久久久国产一区二区| 国产欧美日韩一区二区三区在线| 看免费成人av毛片| 免费高清在线观看日韩| 男女下面插进去视频免费观看| 亚洲伊人久久精品综合| 美女福利国产在线| 亚洲精品日本国产第一区| 中文字幕亚洲精品专区| 只有这里有精品99| 国产精品免费视频内射| 亚洲人成电影观看| 一区二区三区激情视频| 国产精品成人在线| 亚洲国产欧美日韩在线播放| 久久免费观看电影| 久久国产亚洲av麻豆专区| 在线观看www视频免费| 日韩精品免费视频一区二区三区| 国产一区二区三区综合在线观看| 熟女少妇亚洲综合色aaa.| 欧美在线黄色| 欧美国产精品一级二级三级| 亚洲精华国产精华液的使用体验| 亚洲欧洲日产国产| 伊人久久国产一区二区| freevideosex欧美| 天堂8中文在线网| 热re99久久精品国产66热6| 极品少妇高潮喷水抽搐| 成人亚洲精品一区在线观看| av免费观看日本| 看免费av毛片| 免费大片黄手机在线观看| 天美传媒精品一区二区| 女性生殖器流出的白浆| 好男人视频免费观看在线| 欧美激情极品国产一区二区三区| 日韩成人av中文字幕在线观看| a级片在线免费高清观看视频| 伊人久久国产一区二区| 国产福利在线免费观看视频| 老司机亚洲免费影院| 亚洲av日韩在线播放| 男人爽女人下面视频在线观看| 丝袜美腿诱惑在线| 色婷婷久久久亚洲欧美| 天天躁夜夜躁狠狠久久av| 老司机亚洲免费影院| 国产成人精品久久二区二区91 | 国产精品久久久久成人av| 午夜福利在线观看免费完整高清在| 一区二区三区乱码不卡18| 黄色视频在线播放观看不卡| 国产精品 国内视频| 一二三四在线观看免费中文在| 成人毛片a级毛片在线播放| 在现免费观看毛片| 国产成人精品无人区| a 毛片基地| 亚洲激情五月婷婷啪啪| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av涩爱| 极品人妻少妇av视频| 中文字幕精品免费在线观看视频| 一区二区三区激情视频| 亚洲久久久国产精品| 日本猛色少妇xxxxx猛交久久| 飞空精品影院首页| 丝袜美腿诱惑在线| 91久久精品国产一区二区三区| 韩国av在线不卡| 成人亚洲精品一区在线观看| 国产人伦9x9x在线观看 | 男的添女的下面高潮视频| 热re99久久国产66热| 人妻一区二区av| 免费久久久久久久精品成人欧美视频| 丝袜在线中文字幕| 只有这里有精品99| 免费日韩欧美在线观看| 五月开心婷婷网| 成人影院久久| 久久人人97超碰香蕉20202| 美女大奶头黄色视频| 亚洲,一卡二卡三卡| 久久这里有精品视频免费| 欧美人与性动交α欧美精品济南到 | 不卡av一区二区三区| 久久精品人人爽人人爽视色| 丰满迷人的少妇在线观看| 天天影视国产精品| 亚洲欧美成人综合另类久久久| 水蜜桃什么品种好| 亚洲国产精品一区二区三区在线| 国产男女内射视频| 亚洲少妇的诱惑av| 亚洲熟女精品中文字幕| av女优亚洲男人天堂| 老汉色av国产亚洲站长工具| 亚洲精品乱久久久久久| 欧美日韩亚洲高清精品| 久久精品亚洲av国产电影网| 热99国产精品久久久久久7| 大陆偷拍与自拍| 99国产综合亚洲精品| 亚洲激情五月婷婷啪啪| 亚洲综合色网址| 国产乱来视频区| 美女高潮到喷水免费观看| 色哟哟·www| 亚洲av综合色区一区| 国产女主播在线喷水免费视频网站| 国产人伦9x9x在线观看 | 男人操女人黄网站| 老司机影院成人| 久久精品人人爽人人爽视色| 欧美国产精品一级二级三级| 热re99久久精品国产66热6| 一级毛片黄色毛片免费观看视频| 搡老乐熟女国产| 丝瓜视频免费看黄片| 成人毛片a级毛片在线播放| 在线观看美女被高潮喷水网站| 国产又爽黄色视频| 亚洲av.av天堂| 狠狠精品人妻久久久久久综合| 大香蕉久久成人网| 水蜜桃什么品种好| 精品酒店卫生间| 成人漫画全彩无遮挡| 人妻 亚洲 视频| 一边亲一边摸免费视频| 高清在线视频一区二区三区| 91国产中文字幕| 亚洲国产色片| 美女脱内裤让男人舔精品视频| 日韩成人av中文字幕在线观看| 国产又色又爽无遮挡免| 日韩一卡2卡3卡4卡2021年| 一边摸一边做爽爽视频免费| 久久人妻熟女aⅴ| 青青草视频在线视频观看| 亚洲国产欧美日韩在线播放| 亚洲精品视频女| 午夜老司机福利剧场| 国产一区二区激情短视频 | 国产精品三级大全| 亚洲第一av免费看| 色吧在线观看| 婷婷成人精品国产| 免费高清在线观看视频在线观看| 春色校园在线视频观看| 久久久欧美国产精品| av又黄又爽大尺度在线免费看| 97在线人人人人妻| 亚洲美女视频黄频| 香蕉精品网在线| 男的添女的下面高潮视频| 国产免费视频播放在线视频| 久久久精品免费免费高清| av网站在线播放免费| 在线观看免费视频网站a站| 国产精品欧美亚洲77777| 午夜福利一区二区在线看| 国产视频首页在线观看| 久久久久久久亚洲中文字幕| 欧美日韩视频精品一区| 最近的中文字幕免费完整| 亚洲成人av在线免费| 久久久久久久久免费视频了| 少妇精品久久久久久久| 中文欧美无线码| 成年人午夜在线观看视频| 久久久久久久久久人人人人人人| 欧美国产精品一级二级三级| 日韩一区二区视频免费看| 极品人妻少妇av视频| 黄色 视频免费看| 99热国产这里只有精品6| 三级国产精品片| 最新中文字幕久久久久| 精品人妻在线不人妻| 卡戴珊不雅视频在线播放| 免费观看在线日韩| 麻豆乱淫一区二区| 亚洲欧洲国产日韩| 精品亚洲成国产av| 啦啦啦在线观看免费高清www| 丁香六月天网| 菩萨蛮人人尽说江南好唐韦庄| 欧美成人午夜精品| 男的添女的下面高潮视频| 国产在视频线精品| 五月开心婷婷网| 啦啦啦视频在线资源免费观看| 精品少妇黑人巨大在线播放| 婷婷色综合www| 国产精品偷伦视频观看了| 午夜福利视频在线观看免费| 天堂8中文在线网| 午夜福利在线观看免费完整高清在| 波多野结衣一区麻豆| 国产又爽黄色视频| 免费观看性生交大片5| 高清黄色对白视频在线免费看| 最近的中文字幕免费完整| 欧美bdsm另类| 国产av一区二区精品久久| 国产免费视频播放在线视频| 男女边摸边吃奶| 日日啪夜夜爽| www.熟女人妻精品国产| 精品酒店卫生间| 国产精品久久久久成人av| 久久久久久久久久久久大奶| 亚洲欧美日韩另类电影网站| 天天躁夜夜躁狠狠久久av| 国产在视频线精品| 欧美xxⅹ黑人| 久久久久国产精品人妻一区二区| 人成视频在线观看免费观看| 99久久精品国产国产毛片| 肉色欧美久久久久久久蜜桃| 狂野欧美激情性bbbbbb| 搡女人真爽免费视频火全软件| 国产淫语在线视频| av在线播放精品| 欧美xxⅹ黑人| 天天躁夜夜躁狠狠躁躁| 欧美国产精品一级二级三级| 最近中文字幕2019免费版| 妹子高潮喷水视频| 亚洲内射少妇av| 免费高清在线观看日韩| 久久精品久久久久久久性| 精品久久久久久电影网| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲国产毛片av蜜桃av| www.熟女人妻精品国产| 久久久久精品人妻al黑| 久久 成人 亚洲| 成人亚洲精品一区在线观看| 日韩欧美精品免费久久| 老司机影院毛片| 欧美日韩国产mv在线观看视频| 考比视频在线观看| 色哟哟·www| 日韩一区二区视频免费看| 777米奇影视久久| 亚洲国产日韩一区二区| 777久久人妻少妇嫩草av网站| 人人妻人人澡人人看| 亚洲第一青青草原| 久久久久网色| 观看av在线不卡| 久久精品aⅴ一区二区三区四区 | 国产成人午夜福利电影在线观看| videossex国产| av在线观看视频网站免费| 中文字幕亚洲精品专区| 成人影院久久| 久久99一区二区三区| 亚洲av在线观看美女高潮| 国产一区亚洲一区在线观看| 久热这里只有精品99| 午夜免费鲁丝| 国产在线免费精品| 99九九在线精品视频| 日本免费在线观看一区| 成年人午夜在线观看视频| 国产精品偷伦视频观看了| 热re99久久精品国产66热6| 欧美日韩一级在线毛片|