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

    紅壤不同粒徑團聚體對草甘膦農(nóng)藥降解動力學(xué)的影響*

    2021-05-08 03:18:20辛龍川吳文雪付格娟陳延華王學(xué)霞劉東生楊曉梅
    關(guān)鍵詞:紅壤半衰期草甘膦

    辛龍川, 吳文雪, 薛 萐, 付格娟, 陳延華, 王學(xué)霞, 劉東生, 楊曉梅,**

    紅壤不同粒徑團聚體對草甘膦農(nóng)藥降解動力學(xué)的影響*

    辛龍川1, 吳文雪1, 薛 萐2, 付格娟3, 陳延華4, 王學(xué)霞4, 劉東生4, 楊曉梅1,2**

    (1. 西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院 楊凌 712100; 2. 西北農(nóng)林科技大學(xué)水土保持研究所/黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室 楊凌 712100; 3. 西安市環(huán)境監(jiān)測站 西安 710100; 4. 北京市農(nóng)林科學(xué)院植物營養(yǎng)與資源研究所 北京 100097)

    草甘膦農(nóng)藥的大量噴施, 使其在環(huán)境特別是土壤中的殘留-累積風(fēng)險日益突出, 從團聚體粒徑角度研究紅壤不同粒徑團聚體中草甘膦的降解動力學(xué)及其相互作用特征仍鮮有報道?;诖? 本研究通過干篩篩分、室內(nèi)控制培養(yǎng)、液質(zhì)聯(lián)用定量分析相結(jié)合等探究草甘膦降解殘留, 并進一步分析團聚體理化性質(zhì)與草甘膦降解的關(guān)系。結(jié)果表明: 1)不同粒徑團聚體中, 草甘膦殘留量隨降解時間不斷減小, 且粒徑之間降解動力學(xué)差異不顯著。降解半衰期為15.8~20.6 d, 粒徑最小的團聚體(<0.25 mm)中草甘膦的降解半衰期最長, 為20.6 d。草甘膦在土壤中的主要降解產(chǎn)物氨甲基磷酸(AMPA)的含量隨著降解時間的增加而增加, 且在第5 d達(dá)到峰值, 而后不斷減小; 不同粒徑團聚體間AMPA含量差異顯著(<0.05)。2)相關(guān)分析及主成分分析發(fā)現(xiàn), 草甘膦殘留量與紅壤團聚體中速效磷含量呈顯著正相關(guān)(<0.05), 而其降解產(chǎn)物AMPA含量與團聚體中酸性磷酸酶活性及N-乙酰氨基-β-葡萄糖苷酶活性呈顯著正相關(guān)(<0.05)。團聚體粒徑與草甘膦殘留量間沒有顯著相關(guān)性, 但與AMPA含量顯著正相關(guān)(<0.05)。此外, 草甘膦降解過程中, 團聚體中有機質(zhì)含量及β-葡萄糖苷酶、N-乙酰氨基-β-葡萄糖苷酶、酸性磷酸酶活性與團聚體粒徑為顯著負(fù)相關(guān)關(guān)系(<0.05)。由此表明: 紅壤不同粒徑團聚體影響草甘膦降解速率, 粒徑最小的團聚體(<0.25 mm)中草甘膦農(nóng)藥的降解速率最慢, 但試驗結(jié)束時, 各粒徑紅壤團聚體中的草甘膦和AMPA含量均較高, 可能會影響土壤健康及生態(tài)環(huán)境安全; 此外, 草甘膦降解與土壤磷素密切相關(guān), 后續(xù)研究需探討磷虧缺或豐盈條件下, 草甘膦農(nóng)藥的土壤環(huán)境特征, 為后續(xù)農(nóng)田草甘膦環(huán)境風(fēng)險評估提供依據(jù)。

    土壤團聚體; 草甘膦; 降解動力學(xué); 土壤磷; 紅壤

    農(nóng)藥是當(dāng)前全球農(nóng)業(yè)生產(chǎn)中重要的生產(chǎn)資料。草甘膦(C3H8NO5P)因低毒性、非選擇性強、見效快等優(yōu)點被廣泛使用, 特別是隨著抗草甘膦作物的大面積種植, 草甘膦農(nóng)藥已成為全球使用量最大的除草劑[1]。據(jù)統(tǒng)計, 我國草甘膦農(nóng)藥2018年生產(chǎn)量為50萬t, 使用量為4萬t, 2019年生產(chǎn)量更是達(dá)55萬t[2]。然而, 隨著草甘膦的大量施用, 其在農(nóng)田[3-4]及周邊水域[5-6]中不斷累積, 特別是降解產(chǎn)物氨甲基磷酸(AMPA)的檢測率顯著增加, 并隨著雨滴擊濺、剝蝕搬運及徑流輸移進入下游生態(tài)系統(tǒng), 嚴(yán)重威脅著人類健康及生態(tài)環(huán)境安全。

    草甘膦進入土壤后, 因強烈的吸附力, 被固定在土壤顆粒中并參與土壤物理、化學(xué)、生物過程。土壤中草甘膦通過微生物降解, 降解主要有兩種途徑: 其一, 草甘膦在轉(zhuǎn)氨酶的作用下C—N鍵氧化斷裂生成中間產(chǎn)物AMPA; 其二, 草甘膦在C—P裂解酶的作用下生成肌氨酸[7-8], 其中絕大部分降解中間產(chǎn)物為AMPA。草甘膦在土壤中的持久性及降解在一定程度上受土壤理化性質(zhì)的影響, 與土壤類型、pH、有機質(zhì)含量、溫度、含水量等密切相關(guān), 從而草甘膦在土壤中的半衰期差異較大。Mamy等[9]研究發(fā)現(xiàn), 不同類型土壤中草甘膦的降解速率差異顯著(<0.05), 在pH小的砂質(zhì)黏壤土中的殘留量明顯高于其他土壤類型, 并且沙土含量較高的土壤中草甘膦的半衰期小于1 d, 黏土含量較高的土壤中草甘膦的半衰期為3.7 d; 而胡繼業(yè)等[10]研究發(fā)現(xiàn)草甘膦在沙質(zhì)土的半衰期為9.2 d; Yang等[11]研究表明, 草甘膦在中國楊凌黃土中的降解半衰期為3.5 d; 汪立高等[12]對比發(fā)現(xiàn)草甘膦在紅壤、沙壤和黃壤中的半衰期分別為15.71 d、12.86 d和13.30 d。微生物降解是進入土壤草甘膦衰減的重要途徑, 微生物活性及群落特征直接影響著其在土壤中的消散。謝怡等[13]研究發(fā)現(xiàn), 草甘膦在水稻土中的自然降解半衰期為38.72 d, 而滅菌除去原生微生物后, 半衰期可延長至96.27 d。長期施用草甘膦的土壤中, 存在種類繁多的能耐受或降解草甘膦的微生物類群, 如假單胞菌屬()、無色菌屬()、青霉屬()、嗜熱菌屬()等, 它們能夠以草甘膦為碳源、氮源、磷源等生長, 不僅對土壤的物質(zhì)轉(zhuǎn)化起到推動作用, 同時也可以加速草甘膦農(nóng)藥降解[14]。尚軼等[15]從草甘膦污染土壤中分離得到的菌株G1以草甘膦為唯一的碳源, 草甘膦降解率達(dá)71.76%; 石成春等[16]分離得到的黑曲霉B21對草甘膦的降解率達(dá)97%; 湯鳴強等[17]分離得到的以草甘膦為唯一碳源、氮源生長的酵母菌ZM-1對草甘膦降解率為85.38%。然而, 基于不同試驗研究的條件及研究手段, 草甘膦在土壤中的行為追蹤還存在較大的不足, 其衰減及消散仍需進一步研究與探討。

    目前, 針對草甘膦的研究主要集中在土壤類型、整個土體耕層以及環(huán)境因素對其降解衰減的影響。然而, 受田間環(huán)境的影響, 草甘膦及其降解產(chǎn)物會隨著土壤顆粒的剝蝕, 在土壤內(nèi)部或被地表徑流沖刷下搬運、轉(zhuǎn)移并在土壤顆粒上衰減消散, 因此土壤顆粒性質(zhì), 特別是土壤團聚體理化特征對其環(huán)境行為起著重要作用。潘吉秀等[18]研究表明, 六氯苯、多環(huán)芳烴等有機污染物在土壤不同粒徑團聚體中的含量隨粒徑減小先降低后升高, 滴滴涕含量隨粒徑減小先升高后降低; 黃玉芬等[19]研究了土壤不同粒徑團聚體對農(nóng)藥丁草胺的吸附, 發(fā)現(xiàn)其吸附能力由大到小依次為黏粒(<0.002 mm)、粉粒(0.002~0.02 mm)、砂粒(0.02~0.05 mm)。Bento等[20]通過風(fēng)洞模擬試驗研究了草甘膦在風(fēng)蝕過程中的分布及降解情況, 發(fā)現(xiàn)不同土壤粒徑上草甘膦和AMPA含量差異顯著, 且隨粒徑的增加而降低。然而, Bento等[20]的研究缺乏對不同顆粒團聚體性質(zhì)的分析, 且風(fēng)洞試驗?zāi)M時間較短, 無法進一步探明草甘膦在顆粒上的降解衰減動力學(xué)特征, 從而給后續(xù)污染評價帶來諸多的不確定性。鑒于此, 本研究采用室內(nèi)模擬培養(yǎng)試驗, 研究紅壤不同粒徑團聚體中草甘膦的降解動力學(xué)及不同粒徑團聚體理化性質(zhì)與草甘膦降解的關(guān)系, 為進一步探究草甘膦在土壤團聚體中的衰減規(guī)律以及風(fēng)險評估提供參考依據(jù)。

    1 材料與方法

    1.1 試驗土壤

    試驗選用中國農(nóng)業(yè)科學(xué)院祁陽紅壤實驗站周邊農(nóng)田土壤作為試驗土樣(0~30 cm, 田間持水量25%), 去除土壤中的枯落物、根系及礫石等雜質(zhì), 對大土塊進行破碎處理, 并在野外條件下全部過2 mm篩, 風(fēng)干后帶回實驗室, 在采集和運輸過程中盡量減少對土樣的擾動, 以免破壞土壤團聚體。在室內(nèi)采用干篩法, 將土壤進行篩分, 過篩孔徑分別為2 mm、1 mm、0.25 mm、0.05 mm, 其質(zhì)量比分別為26.01%、49.77%、23.80%、0.42%。篩分過程中, 由于<0.05 mm粒徑的土壤較少, 與上一級篩分粒徑0.25 mm的土壤混合使用, 從而得到待試土壤粒徑為: <2 mm (未篩分土樣)、1~2 mm、0.25~1 mm、<0.25 mm, 其基本化學(xué)性質(zhì)如表1所示。

    表1 試驗用紅壤不同粒徑團聚體的基本性質(zhì)

    1.2 研究方法

    將篩分好的土壤分別放入自封袋中, 調(diào)節(jié)土壤水分為田間持水量的60%, 然后放入培養(yǎng)箱中, 調(diào)節(jié)氣候箱溫度為(25±0.5) ℃, 濕度為(80±2)%, 光照16 h, 黑暗8 h, 培養(yǎng)7 d, 每天進行兩次水分調(diào)節(jié), 使其土壤含水量保持穩(wěn)定。預(yù)培養(yǎng)結(jié)束后, 分別稱取4種粒徑團聚體57.5 g放入白色塑料盒中(直徑為7 cm, 高10 cm), 每種粒徑團聚體稱取21盒, 噴施配制好的草甘膦溶液, 濃度為15 μg?g?1, 與土壤混合均勻, 同時設(shè)置空白對照試驗。在培養(yǎng)期間, 監(jiān)測土壤水分變化, 并進行每日兩次調(diào)節(jié), 校對土壤水分含量。在添加農(nóng)藥后的0 d、1 d、3 d、5 d、7 d、14 d和30 d采樣, 采樣時將培養(yǎng)盒中的土壤混合均勻后全部收集, 同種粒徑的團聚體每次采集3盒, 作為試驗重復(fù)樣品。迅速將采集的樣品分成兩部分, 一部分迅速冷凍保存(?20 ℃), 以測試不同粒徑團聚體草甘膦、AMPA含量以及酶活性; 另一部分風(fēng)干, 用來分析土壤團聚體的理化性質(zhì)。

    1.3 草甘膦與AMPA的測定

    1.3.1 試劑

    草甘膦和氨甲基磷酸(AMPA)(純度≥98%)購自美國Sigma-Aldrich公司, 100 μg·L?1的草甘膦同位素內(nèi)標(biāo)物(13C2-15N-GLY)和AMPA同位素內(nèi)標(biāo)物(13C-15N-AMPA)購自德國Augsburg公司, 氯甲酸-9-芴基甲酯(純度≥98%)購自瑞士Sigma-Aldrich公司, 十水合四硼酸鈉(純度≥99.5%)購自美國Sigma-Aldrich公司, 氫氧化鉀和鹽酸(分析純)購自德國Darmstadt公司, 甲醇和乙腈(色譜純)購自荷蘭All Chemicals公司, 甲酸(純度≥98%)購自荷蘭Gevaar公司, 乙酸銨(色譜純)購自美國Sigma公司, 氨溶液(純度25%)購自德國Merck KGaA公司, 試驗用水為超純水。

    1.3.2 草甘膦和AMPA的提取和衍生

    將2 g不同粒級土壤團聚體加入到50 mL塑料離心管中, 加入10 mL 0.6 mol?L?1的KOH對土壤樣本進行浸提, 經(jīng)震蕩機震蕩1 h后放入離心機在3500 r?min?1下離心15 min。吸取1 mL上清液至10 mL離心管中, 加入80 μL 6 mol?L?1的HCl調(diào)節(jié)pH至9左右, 然后加入40 μL的5 μg?mL?1同位素標(biāo)記的草甘膦/AMPA混合溶液, 充分混合后加入0.5 mL濃度為5%的硼酸鹽緩沖液及0.5 mL 6.5 mmol?L?1的FMOC-Cl, 震蕩試管使溶液充分混合均勻, 在室溫條件下衍生反應(yīng)30 min。衍生反應(yīng)結(jié)束后, 加入50 μL蟻酸, 再次混勻試管內(nèi)的溶液, 并將0.5 mL的衍生物轉(zhuǎn)移且過0.45 μm膜的測試瓶中待試。標(biāo)準(zhǔn)曲線采用0 μg?mL?1、0.01 μg?mL?1、0.02 μg?mL?1、0.04 μg?mL?1、0.1 μg?mL?1、0.2 μg?mL?1、0.4 μg?mL?1、1 μg?mL?1、2 μg?mL?1和4 μg?mL?1的標(biāo)準(zhǔn)溶液, 其衍生方式與樣品相同[21]。

    1.3.3 高效液相色譜-串聯(lián)質(zhì)譜法

    色譜柱: XBridge? Shield RP C18色譜柱(孔徑3.5 μm, 長度150 mm, 內(nèi)徑2.1 mm)。流動相為溶于超純水的5 mmol?L?1NH4Ac (溶劑A, pH=9)和比例為9∶1的MeOH∶H2O溶液(溶劑B, pH=9), 使用25% NH3將pH調(diào)整至約9。洗脫程序為: 0~1 min, 100% A∶0 B; 1~6 min, 溶劑B從0到100%線性增加; 6~8 min, 0 A∶100% B; 8~9 min, B從100%線性降低至0; 9~14 min, 100% A∶0 B。柱溫為35.0 ℃, 流速為0.4 mL?min?1, 在1 min內(nèi)重新建立初始條件, 總運行時間為15 min。使用以下離子源參數(shù)在負(fù)離子模式下進行電噴霧電離可獲得最佳響應(yīng): 毛細(xì)管電壓3500 V, 錐管電壓20 V, 離子源溫度120 ℃, 脫溶氣體溫度400 ℃, 氣流160~200 L?h?1(錐管)和580~600 L?h?1(脫溶), 具體參考Yang等[21]的研究。

    1.4 團聚體理化性質(zhì)測定

    團聚體的理化性質(zhì)測試方法如下: 采用0.5 mol?L?1NaHCO3溶液浸提, 鉬銻抗比色法測土壤樣品速效磷含量; 采用重鉻酸鉀容量法測土壤有機質(zhì)含量; 采用H2SO4-HClO4消解法測土壤全磷含量; 采用凱式法消解, 凱式定氮儀測定土壤全氮含量; 采用微孔板熒光法測量土壤中β-葡萄糖苷酶(BG)、N-乙酰氨基-β-葡萄糖苷酶(NAG)和酸性磷酸酶(ACP)活性。

    1.5 數(shù)據(jù)統(tǒng)計與分析

    基于草甘膦的化學(xué)特征, 其在土壤中衰減特征采用一級動力學(xué)降解模型擬合(SFO)。

    C=C0′e (1)

    式中:是時間的草甘膦濃度,0是初始草甘膦濃度(=0),是降解速率常數(shù)。

    化學(xué)污染物的半衰期是評估其是否對環(huán)境造成威脅的重要參數(shù), 因此用草甘膦的半衰期(DT50)表征草甘膦在土壤中的殘留狀況, 并分析其潛在威脅。草甘膦降解半衰期采用式(2)計算。

    DT50=ln2/(2)

    數(shù)據(jù)采用SPSS 26.0進行統(tǒng)計分析, 采用單因素方差分析(One-way ANOVA)比較紅壤不同粒徑團聚體中草甘膦、AMPA含量以及土壤理化性質(zhì)在相同培養(yǎng)時間的顯著性差異, 多重比較采用Duncan法(<0.05); 使用SPSS 26.0進行相關(guān)性分析和主成分分析; 運用OriginPro 9.1進行數(shù)據(jù)繪圖并進行數(shù)據(jù)的擬合。

    2 結(jié)果與分析

    2.1 不同粒徑團聚體中草甘膦及AMPA含量特征及降解動力學(xué)曲線

    紅壤不同粒徑團聚體中草甘膦與AMPA含量隨培養(yǎng)時間的變化差異顯著(圖1)。在培養(yǎng)初期, 4種粒徑團聚體中草甘膦殘留量快速下降, 到噴施后的第5 d, 4種粒徑團聚體中草甘膦殘留量分別降低30.0% (<2 mm)、33.4% (1~2 mm)、37.0% (0.25~ 1 mm)、34.8% (<0.25 mm); 然而隨著觀測時間的延長, 草甘膦殘留量逐漸降低, 各粒徑間草甘膦殘留量差異不顯著(圖1A)。相對應(yīng)的, 4種粒徑團聚體中AMPA含量在噴施初期呈快速增長趨勢, 各粒徑間差異不顯著, 且在第5 d達(dá)濃度峰值; 之后隨著觀測時間的推移, AMPA含量逐漸下降(圖1B), 各粒徑間的含量差異顯著(<0.05)。這是因為在加入草甘膦初期, 草甘膦降解生成AMPA, 土壤中AMPA含量開始增加, 但同時AMPA本身也會發(fā)生降解, 因此隨著觀測時間的推移, 其含量逐漸下降。4種粒徑團聚體中AMPA的含量峰值分別為2.21 μg·g–1(<2 mm)、2.43 μg·g–1(1~2 mm)、2.14 μg·g–1(0.25~1 mm)和2.47 μg·g–1(<0.25 mm); 但隨著時間的推移, AMPA含量隨團聚體粒徑的增加而增加, 以1~2 mm粒徑最高, <0.25 mm粒徑最低。

    對不同粒徑團聚體中草甘膦殘留量采用一級動力學(xué)模型擬和其降解過程(表2)。結(jié)果表明, 不同粒徑團聚體中草甘膦的降解動力學(xué)符合一級動力學(xué)模型,2分別為0.73 (<2 mm)、0.73 (1~2 mm)、0.62 (0.25~1 mm)、0.65 (<0.25 mm)。4種粒徑團聚體中, 草甘膦的半衰期分別為18.5 d (<2 mm)、15.8 d (1~2 mm)、17.7 d (0.25~1 mm)、20.6 d (<0.25 mm)。可以看出, 除未篩分土樣外, 草甘膦半衰期均隨著土壤團聚體粒徑的增加而降低, 而<2 mm的未篩分土樣的草甘膦半衰期介于0.25~1 mm和<0.25 mm團聚體之間。

    表2 紅壤不同粒徑團聚體中草甘膦降解動力學(xué)方程

    2.2 不同粒徑團聚體理化性質(zhì)及酶活性

    2.2.1 理化性質(zhì)

    4種粒徑團聚體土壤理化性質(zhì)如表3所示。由表3可知, 有機質(zhì)含量隨培養(yǎng)時間的增加波動明顯, 特別是噴施草甘膦1~3 d后, 有機質(zhì)含量下降(除未篩分土樣外), 此時其最小值分別為15.84 g?kg?1(1~2 mm)、16.37 g?kg?1(0.25~1 mm)和18.04 g?kg?1(<0.25 mm); 此后逐漸增加, 除未篩分土樣在噴施草甘膦7 d后達(dá)最大值17.30 g?kg?1外, 其余粒徑團聚體有機質(zhì)含量均在噴施后5 d達(dá)最大值, 分別為16.33 g?kg?1(1~2 mm)、17.34 g?kg?1(0.25~1 mm)和19.41 g?kg?1(<0.25 mm); 隨著培養(yǎng)時間的延長, 有機質(zhì)含量逐漸降低, 在噴施后30 d (未篩分土壤除外)各粒徑的含量低于噴施后0 d, 略有減少。<0.25 mm團聚體與0.25~1 mm團聚體、1~2 mm團聚體和未篩分土樣的有機質(zhì)含量差異顯著, 土壤粒徑小的有機質(zhì)含量高。

    添加草甘膦后, 4種粒徑團聚體的全氮含量整體變化相似, 隨培養(yǎng)時間增加先上升后下降, 但變化幅度有所不同(表3), <0.25 mm團聚體全氮含量在噴施后1 d達(dá)最大值(1.21 g?kg?1), 其余粒徑則在噴施后5~7 d達(dá)最大值, 分別為1.15 g?kg?1(<2 mm)、1.20 g?kg?1(1~2 mm)、1.10 g?kg?1(0.25~1 mm)。此后, 隨著培養(yǎng)時間的延長, 全氮含量逐漸減小, 噴施農(nóng)藥30 d后各粒徑中全氮含量高于噴施后0 d, 略有增加。4種粒徑團聚體全氮含量的變化可能是由于不是完全封閉試驗, 在草甘膦降解期間, 土壤中有氮的進入與釋放, 與土壤中利用氮源的微生物有關(guān)。<0.25 mm團聚體與0.25~1 mm、1~2 mm團聚體和未篩分土樣的全氮含量差異顯著, 土壤粒徑小的全氮含量高。

    由表3知, 紅壤各粒徑團聚體間全磷含量差異不顯著。同時, 在農(nóng)藥噴施0~1 d后, 速效磷含量達(dá)最大值, 分別為0.15 g?kg?1(<2 mm)、0.15 g?kg?1(1~2 mm)、0.14 g?kg?1(0.25~1 mm)和0.15 g?kg?1(<0.25 mm); 此后逐漸降低, 于噴施后5 d達(dá)最小值, 分別為0.12 g?kg?1(<2 mm)、0.12 g?kg?1(1~2 mm)、0.13 g?kg?1(0.25~1 mm)和0.12 g?kg?1(<0.25 mm); 然而, 隨著培養(yǎng)時間的延長, 速效磷含量逐漸增加, 噴施農(nóng)藥30 d后各粒徑中速效磷含量與噴施后5 d的含量接近, 且各種粒徑間速效磷含量差異并不顯著。

    表3 噴施草甘膦后不同時間紅壤不同粒徑團聚體的性質(zhì)

    不同小寫字母表示同一指標(biāo)同一時間不同粒徑間差異在<0.05水平顯著。Different lowercase letters indicate significant differences among different aggregate sizes for the same index at the same time at<0.05 level.

    2.2.2 酶活性

    4種粒徑團聚體噴施草甘膦后相對于對照的β-葡萄糖苷酶(BG)活性及變化率見圖2a與圖2b。結(jié)果表明: 各粒徑間BG活性及變化率變化趨勢相似, 隨培養(yǎng)時間波動明顯, 且存在一定的差異。噴施草甘膦后, 首先抑制BG活性, 之后隨培養(yǎng)時間延長, 可能由于微生物對環(huán)境刺激的適應(yīng)性加強, 再加上草甘膦提供了新的碳源, 草甘膦對BG活性由抑制變?yōu)榧せ? BG活性于噴施3 d后達(dá)最大值, 分別為38.18 μmol?kg?1?h?1(<2 mm)、33.46 μmol?kg?1?h?1(1~2 mm)、58.05 μmol?kg?1?h?1(0.25~1 mm)、77.35 μmol?kg?1?h?1(<0.25 mm)。然而噴施后7~14 d草甘膦對BG活性又表現(xiàn)為明顯的抑制作用, BG活性小幅度下降, 與剛噴施時接近。這可能是因為土壤中某些碳源消耗殆盡, 而利用這些碳源的微生物繁殖受到抑制, 但隨著培養(yǎng)時間的延長, 抑制作用逐漸減弱, BG活性逐漸恢復(fù)到對照水平。在相同培養(yǎng)時間, BG活性隨土壤粒徑減小而增大, 較小土壤團聚體中含有較高BG活性。

    由圖2c和圖2d可知, 4種粒徑團聚體N-乙酰氨基-β-葡萄糖苷酶(NAG)活性及其變化率隨時間變化波動明顯, 在噴施草甘膦當(dāng)天, 受草甘膦抑制, NAG活性逐漸下降, 活性變化率于噴施后1 d達(dá)最小值, 分別為?59.90% (<2 mm)、?90.64% (1~2 mm)、?73.55% (0.25~1 mm)、?24.48% (<0.25 mm)。之后, 可能部分微生物能夠利用草甘膦作為新的氮源, 草甘膦對NAG的抑制作用逐漸減弱并轉(zhuǎn)變?yōu)榧せ? NAG活性變化率于噴施后3 d達(dá)最大值, 分別為245.83% (<2 mm)、1047.17% (1~2 mm)、144.09% (0.25~1 mm)、233.82% (<0.25 mm); 此時, NAG活性達(dá)最大值, 分別為21.78 μmol?kg?1?h?1(<2 mm)、31.68 μmol?kg?1?h?1(1~2 mm)、20.48 μmol?kg?1?h?1(0.25~1 mm)、27.23 μmol?kg?1?h?1(<0.25 mm); 然而隨草甘膦不斷降解, NAG活性逐漸恢復(fù)到對照水平。在整個培養(yǎng)過程中, 各粒徑間NAG活性大部分時間隨土壤粒徑減小而增加; 未篩分土樣NAG活性大部分時間介于其他3種土壤粒徑之間。

    4種粒徑團聚體酸性磷酸酶(ACP)活性及變化率見圖2e和圖2f。ACP活性隨時間增加波動明顯, 噴施草甘膦后, ACP活性迅速增加, 于第1 d達(dá)到峰值, 分別為125.82 μmol?kg?1?h?1(<2 mm)、92.50 μmol?kg?1?h?1(1~2 mm)、87.34 μmol?kg?1?h?1(0.25~1 mm)、162.90 μmol?kg?1?h?1(<0.25 mm); 此后ACP活性逐漸下降。而噴施草甘膦后3~30 d與0~3 d變化趨勢相似, 于7 d達(dá)到峰值, 分別為157.90 μmol?kg?1?h?1(<2 mm)、86.14 μmol?kg?1?h?1(1~2 mm)、102.83 μmol?kg?1?h?1(0.25~1 mm)、193.58 μmol?kg?1?h?1(<0.25 mm), 與第1 d活性相比較高; 值得注意的是, 噴施農(nóng)藥30 d后各粒徑中ACP活性與噴施后第0 d的活性接近, 趨于增加。不同培養(yǎng)時期, 土壤ACP活性隨團聚體粒徑的減小而增加。由圖2f可知, 在培養(yǎng)過程中, 添加草甘膦對4種粒徑團聚體ACP活性大部分時間都表現(xiàn)出抑制作用, ACP活性在試驗結(jié)束后仍未恢復(fù)到對照水平。

    2.3 紅壤團聚體特征與草甘膦降解特征的相關(guān)性

    運用SPSS對草甘膦在團聚體土壤中的殘留量與團聚體粒徑、AMPA含量、速效磷含量、有機質(zhì)含量、相關(guān)酶活性等進行相關(guān)性分析(表4), 結(jié)果表明: 團聚體粒徑與土壤有機質(zhì)、BG活性、NAG活性、ACP活性顯著負(fù)相關(guān)(<0.05), 但團聚體粒徑與草甘膦殘留量間沒有顯著相關(guān)性, 團聚體粒徑與AMPA含量顯著正相關(guān)(<0.05)。微生物在降解草甘膦時, 會影響土壤性質(zhì)。從表4可知, 草甘膦在土壤中的殘留量與速效磷含量呈顯著正相關(guān)(<0.05), 與土壤有機質(zhì)含量和BG、NAG、ACP活性無顯著相關(guān)性。AMPA含量與NAG、ACP活性顯著正相關(guān)(<0.05)。

    對草甘膦殘留量、速效磷含量、有機質(zhì)含量、BG活性、NAG活性、ACP活性進行主成分分析, 并提取兩種主成分, 其貢獻率分別為41.56%和21.24%。不同粒徑團聚體草甘膦殘留量和速效磷含量、有機質(zhì)含量、BG活性、NAG活性、ACP活性的主成分得分在第2主成分出現(xiàn)明顯的分布差異: <0.25 mm團聚體分布于第2主成分正方向上, 1~2 mm團聚體、0.25~1 mm團聚體分布于第2主成分負(fù)方向上(圖3a)。而從圖3b看出, 對第1主成分起主要作用的指標(biāo)是BG酶活性、NAG酶活性等, 命名為土壤酶活性; 對第2主成分起主要作用的指標(biāo)是草甘膦殘留量與速效磷含量, 命名為土壤磷含量。

    表4 團聚體粒徑與草甘膦及氨甲基磷酸(AMPA)含量的相關(guān)性分析

    BG: β-葡萄糖苷酶; NAG: N-乙酰氨基-β-葡萄糖苷酶; ACP: 酸性磷酸酶。*和**分別表示顯著性水平為<0.05和<0.01。BG: β-glucosidase; NAG: N-acetylamino-β-glucosidase; ACP: acid phosphatase. * and ** indicate significant correlation at<0.05 and<0.01 levels, respectively.

    Gly: 草甘膦殘留量; OM: 有機質(zhì)含量; AP: 速效磷含量; BG: β-葡萄糖苷酶活性; NAG: N-乙酰氨基-β-葡萄糖苷酶活性; ACP: 酸性磷酸酶活性。Gly: glyphosate residue; OM: organic matter content; AP: available phosphorus content; BG: β-glucosidase activity; NAG: N-acetylamino-β-glucosidase activity; ACP: acid phosphatase activity.

    3 討論

    3.1 紅壤不同粒徑團聚體中草甘膦及AMPA的變化特征

    隨著草甘膦使用量的增加, 農(nóng)田土壤中草甘膦及其代謝產(chǎn)物AMPA的出現(xiàn)引發(fā)了人們對草甘膦環(huán)境風(fēng)險的高度關(guān)注與研究, 而從不同粒徑微觀尺度研究草甘膦在土壤中的衰減特征有助于更深一步了解草甘膦在土壤中的環(huán)境行為和后續(xù)評估農(nóng)田侵蝕性污染風(fēng)險。不同粒徑團聚體特性不同, 對有機污染物的降解、吸附等特征不同, 在養(yǎng)分的保持和供應(yīng)中的作用不同[22]。有研究發(fā)現(xiàn), 大部分重金屬如Cd、Cu、Pb和Zn等在土壤中的含量都隨著土壤粒徑的減小而增加[23-25], 部分農(nóng)藥如六六六(BHC)、滴滴涕(DDT)等的殘留量隨土壤粒徑的減小而增大[26], 而部分農(nóng)藥如氯丹等的殘留量則隨土壤粒徑的減小先減小后增大[27], 說明土壤粒徑對不同污染物的影響不同。本研究結(jié)果表明, 紅壤不同粒徑團聚體中的草甘膦殘留量差異不顯著, 說明粒徑對草甘膦在紅壤團聚體中的殘留量影響不明顯。但噴施后期, 在粒徑小的團聚體中, 草甘膦殘留量略高于較大粒級的土壤團聚體, 這一結(jié)果可能與土壤對草甘膦的吸附有關(guān)。研究表明, 草甘膦的吸附過程主要通過配位體交換, 與礦物化合物、有機質(zhì)等形成絡(luò)合物。一方面粒徑小的團聚體比表面積大, 吸附點位多, 對草甘膦的吸附量相比于粒徑大的團聚體不容易達(dá)到飽和; 另一方面, 有機質(zhì)含量與吸附量呈正相關(guān), 而粒徑小的團聚體有機質(zhì)含量高, 草甘膦更容易被吸附, 且吸附后的回收率也相對較低[28]。另外, 草甘膦進入土壤后易降解, 在施用草甘膦2 h內(nèi), 草甘膦會立即降解產(chǎn)生AMPA[29]; 同樣, Paudel等[30]發(fā)現(xiàn), AMPA會在反應(yīng)的1 min內(nèi)形成。本研究在噴施草甘膦當(dāng)天檢測到了AMPA則佐證了這一結(jié)論。

    Yang等[11]研究表明, 自然條件下草甘膦在黃土中的殘留量隨時間變化差異顯著, 初期快速下降, 隨后緩慢下降, 符合一級動力學(xué)降解模型, 而AMPA含量則隨時間增加先快速上升后逐漸下降, 這一結(jié)果與本試驗結(jié)果類似, 說明草甘膦衰減和AMPA形成與降解在不同土壤類型中的總體差異不大, 但擬合的降解模型的值卻與Yang等[11]相比相差1個數(shù)量級。本研究擬合結(jié)果可以看出, 不同粒徑團聚體草甘膦降解動力學(xué)擬合參數(shù)2均低于0.8; 從初始濃度0可以判斷, 模型預(yù)測值與其實測值有所差別, 因此在后續(xù)研究或模型應(yīng)用時, 應(yīng)該考慮預(yù)測效率。通過擬合方程, 本研究得到的草甘膦在紅壤4種粒徑團聚體中的半衰期分別為18.5 d (未篩分土樣)、15.8 d (1~2 mm)、17.7 d (0.25~1 mm)、20.6 d (<0.25 mm), 均比Yang等[11]計算的半衰期長, 介于沙土(16.9 d)和黏土(110 d)之間[31], 且半衰期隨著團聚體粒徑的減小而增大。草甘膦半衰期在不同研究中存在差異的原因可能有兩種: 一是在實驗室中估計的草甘膦半衰期不能真實反映田間草甘膦的降解, 特別是在環(huán)境系統(tǒng)的動態(tài)變化[32]。汪立高等[12]報道草甘膦在紅壤自然環(huán)境下的半衰期為15.71 d, 比本研究在實驗室中的模擬結(jié)果低; Litz等[33]報道在8 ℃的實驗室條件下草甘膦的半衰期為30 d, 其衰減率為0.023 d?1, 但是高溫野外條件下草甘膦的衰減更快, 半衰期為7.5~10.5 d。二是不同粒徑團聚體性質(zhì)存在差異(表1), 土壤性質(zhì)的復(fù)雜作用會影響草甘膦在土壤中的降解, 如土壤pH[34]等。

    3.2 不同粒徑團聚體中草甘膦降解與土壤性質(zhì)和酶活性的關(guān)系

    通常來說, 微生物的代謝會影響草甘膦在土壤中的吸附、遷移、降解等, 進而影響土壤的物質(zhì)轉(zhuǎn)化, 使土壤的一些理化性質(zhì)發(fā)生相應(yīng)的變化。研究表明, 短期使用93.7 μg?g?1的草甘膦, 使紅壤全氮含量增加, 速效磷、有機質(zhì)含量減少[34], 這與本研究結(jié)果類似。毛美紅等[35]研究發(fā)現(xiàn), 長期施用草甘膦的土壤, 有機質(zhì)有所減少, 但全氮、全磷、速效磷無顯著變化; 林德喜等[36]發(fā)現(xiàn)噴施1年草甘膦, 土壤性質(zhì)變化不明顯, 而噴施2年, 土壤有機質(zhì)、全氮等都有不同程度的提高。兩個研究同屬于長期施用草甘膦對土壤性質(zhì)的影響研究, 但結(jié)果卻有所不同,說明草甘膦在土壤中的不斷積累對土壤性質(zhì)的影響是復(fù)雜的, 進行風(fēng)險評估時要關(guān)注草甘膦的施用年限。另外, 草甘膦與磷酸鹽存在著吸附位點的競爭, 草甘膦的存在占據(jù)了磷在土壤上的吸附位點, 且草甘膦濃度越高競爭得到的吸附位點越多[37], 因此速效磷含量隨草甘膦的降解而減少, 兩者呈顯著正相關(guān)(<0.05, 表4)。劉合明等[38]研究發(fā)現(xiàn)在3種農(nóng)田土壤中不同粒徑團聚體有機碳含量差異顯著, 且隨粒徑減少而增加, 這與本研究結(jié)果一致。但馬瑞萍等[39]發(fā)現(xiàn)有機碳在森林帶黃土隨著團聚體粒徑的減小而先逐漸增加再減少, 在草原帶黃土中則隨粒徑的減小而增加; 同時黃榮珍等[40]研究發(fā)現(xiàn)不同作物紅壤團聚體有機質(zhì)隨粒徑減小先減少后增加再減少。造成這種差異的原因可能有兩種: 一是團聚體有機質(zhì)含量可能因植被方式的變化而不同; 二是篩分選擇的粒徑不同, 研究的團聚體粒徑大小不統(tǒng)一。

    土壤酶主要來自于土壤中微生物的生命活動, 參與了土壤中絕大多數(shù)生化過程, 在碳、氮、磷循環(huán)過程發(fā)揮著重要的作用, 與土壤微生物一起對土壤的物質(zhì)轉(zhuǎn)化起到推動作用[41]。通常認(rèn)為, β-葡萄糖苷酶(BG)、N-乙酰氨基-β-葡萄糖苷酶(NAG)和酸性磷酸酶(ACP)等是微生物分泌的碳、氮和磷獲取酶, 其活性可以作為微生物獲取碳、氮和磷的資源分配的替代指標(biāo)[42-43], 因為其他相關(guān)酶的活性通常低于這些酶并與之相關(guān)。草甘膦進入土壤后, 會與土壤酶發(fā)生反應(yīng), 酶活性可能被激活或被抑制, 其特征變化可反映草甘膦在土壤中的環(huán)境行為及土壤污染程度[44]。磷酸酶能促進分解各種有機磷化合物。呼蕾等[45]研究發(fā)現(xiàn)50~2000 μg?g?1的草甘膦、鄧曉等[46]研究發(fā)現(xiàn)50~500 μg?g?1的草甘膦, 都抑制ACP活性, 且隨著草甘膦濃度的增加抑制程度加深。而本研究噴施的草甘膦濃度為15 μg?g?1, 草甘膦對ACP也表現(xiàn)出抑制作用(圖2c), 說明草甘膦在低濃度也抑制ACP活性, 可能是因為草甘膦中存在磷酸基團, 與磷酸酶存在競爭關(guān)系, 但需要進一步研究。紅壤不同粒經(jīng)團聚體磷酸酶活性隨粒徑減小而增大, 這是因為粒徑小的團聚體表面吸附有大量的草甘膦, 降低了草甘膦對磷酸酶活性的抑制作用。另外, 草甘膦在降解過程中, BG和NAG活性波動變化大, 草甘膦對其既有激活作用, 也有抑制作用, 但最后兩種酶活性均恢復(fù)到未噴施草甘膦的水平。

    4 結(jié)論

    1)紅壤不同粒徑團聚體中草甘膦的衰減符合一級動力學(xué)降解模型, 其半衰期分別為18.5 d (未篩分土樣)、15.8 d (1~2 mm)、17.7 d (0.25~1 mm)、20.6 d (<0.25 mm), 半衰期隨團聚體粒徑的減小而增加(未篩分土樣外); 而AMPA含量隨培養(yǎng)時間增加先上升后下降。不同粒徑間草甘膦殘留量差異不顯著, 而不同粒徑間AMPA含量差異顯著(<0.05), 隨團聚體粒徑增大而增加。

    2)噴施草甘膦后, 不同粒徑團聚體有機質(zhì)、全氮及速效磷含量變化較大, 特別是速效磷含量隨著草甘膦農(nóng)藥的降解呈現(xiàn)強烈的波動, 這可能是因為草甘膦與速效磷在土壤顆粒表面吸附點位類似, 相互之間存在競爭導(dǎo)致的; 但后期隨著草甘膦的降解及礦化, 速效磷含量波動相對較小。而受噴施草甘膦影響, 土壤酸性磷酸酶活性明顯受到抑制, 而β-葡萄糖苷酶活性表現(xiàn)為先增加后減小的趨勢, N-乙酰氨基-β-葡萄糖苷酶活性波動較大。此外, 主成分及相關(guān)分析表明, 團聚體粒徑與有機質(zhì)含量及β-葡萄糖苷酶、N-乙酰氨基-β-葡萄糖苷酶、酸性磷酸酶活性呈顯著負(fù)相關(guān)(<0.05), 而團聚體粒徑與AMPA含量呈顯著正相關(guān)(<0.05), 草甘膦在土壤團聚體中的殘留量與速效磷含量呈顯著正相關(guān)(<0.05), 草甘膦的降解與土壤中的磷密切相關(guān)。因此, 后續(xù)還應(yīng)探討磷虧缺或豐盈條件下, 草甘膦的環(huán)境特征, 以期為后續(xù)農(nóng)田草甘膦環(huán)境風(fēng)險評估提供依據(jù)。

    [1] 顧安樂, 柏亞羅. 水稻用農(nóng)藥市場概況、產(chǎn)品開發(fā)及重點品種(Ⅱ)[J]. 現(xiàn)代農(nóng)藥, 2018, 17(6): 1–5 GU A L, BAI Y L. The market profile, product development and key products of the rice pesticide[J]. Modern Agrochemicals, 2018, 17(6): 1–5

    [2] 楊益軍, 張波, 吳江. 全球(中國)草甘膦行業(yè)(市場)發(fā)展?fàn)顩r及2020年展望[J]. 農(nóng)藥市場信息, 2020, (4): 31–32 YANG Y J, ZHANG B, WU J. Global (China) glyphosate industry (market) development status and 2020 outlook[J]. Pesticide Market News, 2020, (4): 31–32

    [3] ALONSO L L, DEMETRIO P M, ETCHEGOYEN M A, et al. Glyphosate and atrazine in rainfall and soils in agroproductive areas of the pampas region in Argentina[J]. Science of The Total Environment, 2018, 645: 89–96

    [4] PRIMOST J E, MARINO D J G, APARICIO V C, et al. Glyphosate and AMPA, “pseudo-persistent” pollutants under real-world agricultural management practices in the Mesopotamic Pampas agroecosystem, Argentina[J]. Environmental Pollution, 2017, 229: 771–779

    [5] OKADA E, ALLINSON M, BARRAL M P, et al. Glyphosate and aminomethylphosphonic acid (AMPA) are commonly found in urban streams and wetlands of Melbourne, Australia[J]. Water Research, 2020, 168: 115139

    [6] 魏琛, 宋麗婧, 楊衛(wèi)萍, 等. 貴陽市地表水中草甘膦農(nóng)藥殘留研究[J]. 環(huán)境科學(xué)與技術(shù), 2016, 39(3): 126–130 WEI C, SONG L J, YANG W P, et al. Research on glyphosate pesticide residue in surface water in Guiyang[J]. Environmental Science & Technology, 2016, 39(3): 126–130

    [7] CARLES L, GARDON H, JOSEPH L, et al. Meta-analysis of glyphosate contamination in surface waters and dissipation by biofilms[J]. Environment International, 2019, 124: 284–293

    [8] BORGGAARD O K, GIMSING A L. Fate of glyphosate in soil and the possibility of leaching to ground and surface waters: A review[J]. Pest Management Science, 2008, 64(4): 441–456

    [9] MAMY L, BARRIUSO E, GABRIELLE B. Environmental fate of herbicides trifluralin, metazachlor, metamitron and sulcotrione compared with that of glyphosate, a substitute broad spectrum herbicide for different glyphosate-resistant crops[J]. Pest Management Science, 2005, 61(9): 905–916

    [10] 胡繼業(yè), 趙殿英, 寧君, 等. 氣相色譜-氮磷檢測器測定草甘膦在土壤和蘋果中的殘留量[J]. 農(nóng)藥學(xué)學(xué)報, 2007, 9(3): 285–290 HU J Y, ZHAO D Y, NING J, et al. Determination of glyphosate residues in soil and apple by capillary gas chromatography with nitrogen-phosphorus detection[J]. Chinese Journal of Pesticide Science, 2007, 9(3): 285–290

    [11] YANG X M, WANG F, BENTO C P M, et al. Decay characteristics and erosion-related transport of glyphosate in Chinese loess soil under field conditions[J]. Science of The Total Environment, 2015, 530/531: 87–95

    [12] 汪立高, 楊仁斌, 魏鳳. 土壤中殘留草甘膦檢測方法及其消解動態(tài)研究[J]. 湖南農(nóng)業(yè)科學(xué), 2011, (23): 85–88 WANG L G, YANG R B, WEI F. Detection method and degradation dynamics of isopropylamine salt of glyphosate in soil[J]. Hunan Agricultural Sciences, 2011, (23): 85–88

    [13] 謝怡, 陳杰, 黃祥, 等. 降解菌ESG4對草甘膦污染土壤的生物修復(fù)研究[J]. 環(huán)境影響評價, 2018, 40(5): 84–87 XIE Y, CHEN J, HUANG X, et al. Research on bioremediation of glyphosate contaminated soil bysp. ESG4[J]. Environmental Impact Assessment, 2018, 40(5): 84–87

    [14] OBOJSKA A, TERNAN N G, LEJCZAK B, et al. Organophosphonate utilization by the thermophileT20[J]. Applied and Environmental Micro-biology, 2002, 68(4): 2081–2084

    [15] 尚軼, 沈慧敏, 楊順義, 等. 草甘膦降解菌的分離及其降解效能研究[J]. 湖北農(nóng)業(yè)科學(xué), 2011, 50(9): 1770–1772 SHANG Y, SHEN H M, YANG S Y, et al. Screening and isolation of bacteria degrading herbicide glyphosate and its degradation efficiency[J]. Hubei Agricultural Sciences, 2011, 50(9): 1770–1772

    [16] 石成春, 郭養(yǎng)浩, 王大奈, 等. 草甘膦曲霉生物降解的動力學(xué)研究[J]. 中國環(huán)境科學(xué), 2005, 25(3): 361–365 SHI C C, GUO Y H, WANG D N, et al. Studies on the kine-tics of the biodegradation of glyphosate by[J]. China Environmental Science, 2005, 25(3): 361–365

    [17] 湯鳴強, 尤民生. 抗草甘膦酵母菌ZM-1的分離鑒定及其生長降解特性[J]. 微生物學(xué)通報, 2010, 37(9): 1402–1409 TANG M Q, YOU M S. Isolation and identification of a glyphosate-resistant yeasty strain (ZM-1) and its optimal growth and degradation traits[J]. Microbiology China, 2010, 37(9): 1402–1409

    [18] 潘吉秀, 夏天翔, 姜林, 等. 工業(yè)污染土壤中不同粒徑下的三種典型有機污染物分布規(guī)律研究[J]. 土壤通報, 2014, 45(2): 462–468 PAN J X, XIA T X, JIANG L, et al. Distribution of organic pollutants in different size soil particles from industrial sites[J]. Chinese Journal of Soil Science, 2014, 45(2): 462–468

    [19] 黃玉芬, 劉忠珍, 魏嵐, 等. 土壤不同粒徑有機無機復(fù)合體對丁草胺的吸附特性[J]. 土壤學(xué)報, 2017, 54(2): 400–409 HUANG Y F, LIU Z Z, WEI L, et al. Effect of soil organo-inorganic compounds different in particle size on butachlor sorption[J]. Acta Pedologica Sinica, 2017, 54(2): 400–409

    [20] BENTO C P M, GOOSSENS D, REZAEI M, et al. Glyphosate and AMPA distribution in wind-eroded sediment derived from loess soil[J]. Environmental Pollution, 2017, 220: 1079–1089

    [21] YANG X M, WANG F, BENTO C P M, et al. Short-term transport of glyphosate with erosion in Chinese loess soil — A flume experiment[J]. Science of The Total Environment, 2015, 512/513: 406–414

    [22] 張菁菁, 溫蓓, 單孝全. 中國典型土壤不同粒徑級份對六氯苯的吸附行為研究[C]//持久性有機污染物論壇2008暨第三屆持久性有機污染物全國學(xué)術(shù)研討會論文集. 北京: 中國化學(xué)會, 中國環(huán)境科學(xué)學(xué)會, 2008: 83–85 ZHANG J J, WEN B, SHAN X Q. Adsorption behavior of hexachlorobenzene on different particle size fractions of typical soils in China[C]//Proceedings of the Persistent Organic Pollutants Forum 2008 and the Third National Symposium on Persistent Organic Pollutants. Beijing: Chinese Chemical Society, Chinese Society for Environmental Sciences, 2008: 83–85

    [23] ZONG Y T, XIAO Q, LU S G. Distribution, bioavailability, and leachability of heavy metals in soil particle size fractions of urban soils (northeastern China)[J]. Environmental Science and Pollution Research, 2016, 23(14): 14600–14607

    [24] 龔倉, 徐殿斗, 成杭新, 等. 典型熱帶林地土壤團聚體顆粒中重金屬的分布特征及其環(huán)境意義[J]. 環(huán)境科學(xué), 2013, 34(3): 1094–1100 GONG C, XU D D, CHENG H X, et al. Distribution characteristics and environmental significance of heavy metals in soil particle size fractions from tropical forests in China[J]. Environmental Science, 2013, 34(3): 1094–1100

    [25] FAN J L, DING W X, ZIADI N. Thirty-year manuring and fertilization effects on heavy metals in black soil and soil aggregates in northeastern China[J]. Communications in Soil Science and Plant Analysis, 2013, 44(7): 1224–1241

    [26] 李勛光. 農(nóng)藥在土壤中的分布——DDT、BHC量與土粒粒徑、有機質(zhì)量之間的關(guān)系[J]. 土壤學(xué)報, 1993, 30(1): 88–93 LI X G. Distribution of pesticides in soils — relationship among contents of DDT and BHC, size of soil particle and content of organic matter[J]. Acta Pedologica Sinica, 1993, 30(1): 88–93

    [27] 張勝田, 趙斌, 王風(fēng)賀, 等. 不同粒徑土壤對氯丹的吸附性能及其急性毒性[J]. 環(huán)境工程學(xué)報, 2017, 11(6): 3839–3845 ZHANG S T, ZHAO B, WANG F H, et al. Effects of soil particles on their adsorption performance of chloredane and their acute toxicity[J]. Chinese Journal of Environmental Engineering, 2017, 11(6): 3839–3845

    [28] 孫立思, 王娜, 孔德洋, 等. 土壤理化性質(zhì)對草甘膦殘留檢測的影響[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 2017, 33(9): 860–864 SUN L S, WANG N, KONG D Y, et al. Influence of soil physical and chemical properties on detection of glyphosate residue[J]. Journal of Ecology and Rural Environment, 2017, 33(9): 860–864

    [29] BENTO C P M, YANG X M, GORT G, et al. Persistence of glyphosate and aminomethylphosphonic acid in loess soil under different combinations of temperature, soil moisture and light/darkness[J]. Science of The Total Environment, 2016, 572: 301–311

    [30] PAUDEL P, NEGUSSE A, JAISI D P. Birnessite-catalyzed degradation of glyphosate: A mechanistic study aided by kinetics batch studies and NMR spectroscopy[J]. Soil Science Society of America Journal, 2015, 79(3): 815–825

    [31] BERGSTR?M L, B?RJESSON E, STENSTR?M J. Laboratory and lysimeter studies of glyphosate and aminomethylphosphonic acid in a sand and a clay soil[J]. Journal of Environmental Quality, 2011, 40(1): 98–108

    [32] LANCT?T C, ROBERTSON C, NAVARRO-MARTíN L, et al. Effects of the glyphosate-based herbicide Roundup WeatherMax?on metamorphosis of wood frogs () in natural wetlands[J]. Aquatic Toxicology, 2013, 140/141: 48–57

    [33] LITZ N T, WEIGERT A, KRAUSE B, et al. Comparative studies on the retardation and reduction of glyphosate during subsurface passage[J]. Water Research, 2011, 45(10): 3047–3054

    [34] 陳望舒. 草甘膦對土壤生態(tài)的影響和毒理研究[D]. 揚州: 揚州大學(xué), 2019: 87–88 CHEN W S. Research on the effects of glyphosate on soil ecology and toxicology[D]. Yangzhou: Yangzhou University, 2019: 87–88

    [35] 毛美紅, 俞婷婷, 傅柳方, 等. 草甘膦對毛竹筍用林土壤理化性質(zhì)的影響分析[J]. 竹子研究匯刊, 2011, 30(3): 29–32 MAO M H, YU T T, FU L F, et al. The influence of glyphosate on soil physicochemical properties of the shoot-oriented Moso bamboo forests[J]. Journal of Bamboo Research, 2011, 30(3): 29–32

    [36] 林德喜, 洪長福, 黃龍發(fā). 杉木幼林地噴施草甘膦后土壤的變化[J]. 土壤, 1998, (2): 100–102 LIN D X, HONG C F, HUANG L F. Soil changes after spraying glyphosate on young Chinese fir woodland[J]. Soils, 1998, (2): 100–102

    [37] 周垂帆, 林靜雯, 李瑩, 等. 磷與草甘膦在酸性土壤中吸附解吸交互作用機制[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2016, 35(12): 2367–2376 ZHOU C F, LIN J W, LI Y, et al. Competitive adsorption and desorption of glyphosate and phosphate on acid soil[J]. Journal of Agro-Environment Science, 2016, 35(12): 2367–2376

    [38] 劉合明, 楊志新, 劉樹慶. 不同粒徑土壤活性有機碳測定方法的探討[J]. 生態(tài)環(huán)境, 2008, 17(5): 2046–2049 LIU H M, YANG Z X, LIU S Q. Methods for determining labile orange matter in different sized soil particles of different soils[J]. Ecology and Environment, 2008, 17(5): 2046–2049

    [39] 馬瑞萍, 劉雷, 安韶山, 等. 黃土丘陵區(qū)不同植被群落土壤團聚體有機碳及其組分的分布[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2013, 21(3): 324–332 MA R P, LIU L, AN S S, et al. Soil organic carbon and its fractions in aggregates under different plant communities in the hill-gully region of the loess plateau[J]. Chinese Journal of Eco-Agriculture, 2013, 21(3): 324–332

    [40] 黃榮珍, 朱麗琴, 黃國敏, 等. 不同整治方式下紅壤坡耕地土壤團聚體有機碳的變化[J]. 南昌工程學(xué)院學(xué)報, 2018, 37(6): 57–62 HUANG R Z, ZHU L Q, HUANG G M, et al. Soil organic carbon change in soil aggregates after different farming practices on slope field in red soil area[J]. Journal of Nanchang Institute of Technology, 2018, 37(6): 57–62

    [41] 魯晶, 楊學(xué)春. 草甘膦對環(huán)境的影響研究進展[J]. 安徽農(nóng)學(xué)通報, 2017, 23(8): 71–75 LU J, YANG X C. Research progress on environmental impact of glyphosate[J]. Anhui Agricultural Science Bulletin, 2017, 23(8): 71–75

    [42] FANIN N, MOORHEAD D, BERTRAND I. Eco-enzymatic stoichiometry and enzymatic vectors reveal differential C, N, P dynamics in decaying litter along a land-use gradient[J]. Biogeochemistry, 2016, 129(1/2): 21–36

    [43] DENG L, PENG C H, HUANG C B, et al. Drivers of soil mic-robial metabolic limitation changes along a vegetation restoration gradient on the Loess Plateau, China[J]. Geoderma, 2019, 353: 188–200

    [44] 崔玉俠. 草甘膦與Cu的復(fù)合污染及其對土壤酶與微生物碳/氮的影響[D]. 重慶: 西南大學(xué), 2009: 35–38 CUI Y X. Study of combined pollution between glyphosate and Cu and its influence on the soil enzyme and microorga-nism carbon/nitrogen[D]. Chongqing: Southwest University, 2009: 35–38

    [45] 呼蕾, 和文祥, 王旭東, 等. 草甘膦的土壤酶效應(yīng)研究[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2009, 28(4): 680–685 HU L, HE W X, WANG X D, et al. Effect of glyphosate on soil enzyme[J]. Journal of Agro-Environment Science, 2009, 28(4): 680–685

    [46] 鄧曉, 武春媛, 李怡, 等. 土壤微生物種群與酶活性對草銨膦和草甘膦脅迫響應(yīng)的差異[J]. 農(nóng)藥, 2019, 58(8): 580–583 DENG X, WU C Y, LI Y, et al. Differences in the stress responses of soil microbial populations and enzyme activities to glufosinate-ammonium and glyphosate[J]. Agrochemicals, 2019, 58(8): 580–583

    Effects of aggregate size on kinetics of glyphosate degradation in red soil*

    XIN Longchuan1, WU Wenxue1, XUE Sha2, FU Gejuan3, CHEN Yanhua4, WANG Xuexia4, LIU Dongsheng4, YANG Xiaomei1,2**

    (1. College of Natural Resources and Environment, Northwest A&F University, Yangling 712100, China; 2. Institute of Soil and Water Conservation, Northwest A&F University / The State Key Laboratory of Soil Erosion and Dryland Farming on Loess Plateau, Yangling 712100, China; 3. Xi’an Environmental Monitoring Station, Xi’an 710100, China; 4. Institute of Plant Nutrition and Resources, Beijing Academy of Agricultural and Forestry Sciences, Beijing 100097, China)

    With intensive glyphosate application, its residues and consequent risks of soil health and ecological environment safety have received greater attention. The degradation kinetics of glyphosate in red soil aggregates with different sizes, as well as the interaction between physical and chemical properties of soil aggregates and the degradation of glyphosate, have rarely been studied. Thus, in this study, the degradation characteristics of glyphosate in red soil aggregates with different sizes were observed under laboratory conditions by particle pre-sieving, incubation in a controlled climatic chamber, and residue analysis via liquid chromatography-tandem mass spectrometry. The physical and chemical properties of the soil aggregates, such as contents of organic matter, total phosphorus, and available phosphorus, were tested according to the national approved methods and standards. The relationships between the physical and chemical properties of the aggregates and the degradation of glyphosate were further analyzed and compared in the same observation day. The results showed that 1) the glyphosate content decreased in the different aggregate particles during the observation period, following the single first-order kinetic degradation model. However, no significant differences were observed among different aggregate sizes. The half-life time of glyphosate in the different red soil aggregates ranged from 15.8 to 20.6 d, with a longer half-life time in the smallest aggregates (<0.25 mm, 20.6 d). The aminomethylphosphonic acid (AMPA) content, the main metabolite of glyphosate, increased immediately and peaked on the 5thday after glyphosate application, but no differences were found among different aggregates. However, the AMPA content changed and declined significantly in different aggregates after the 5thobservation day (<0.05). The contents of organic matter, total nitrogen, total phosphorus, and available phosphorus in different aggregates varied greatly, especially the available phosphorus content, which decreased with glyphosate degradation. 2) Correlation analysis and principal component analysis of glyphosate, aggregate size, and their properties showed that the residual glyphosate was significantly positively correlated with the content of available phosphorus (<0.05), and the AMPA content was significantly positively correlated with the activities of acid phosphatase and N-acetylamino-β-glucosidase (<0.05). There were no significant relationships between the aggregate size and the residuals of glyphosate, but a significant positive correlation was observed between the aggregate size and the AMPA content (<0.05). Furthermore, during the whole period of glyphosate degradation, the organic matter content, acid phosphatase, N-acetylamino-β-glucosidase, and β-glucosidase showed a significant negative relationship with the soil aggregate size (<0.05). In conclusion, the characteristics of the red soil aggregates affect the degradation kinetics of glyphosate, as well as the persistence of AMPA, especially the residuals in the smallest aggregates (<0.25 mm). The contents of glyphosate and AMPA in the red soil aggregates were still high after 30 days, which may affect soil health. Glyphosate degradation was also closely related to phosphorus in the soil. Therefore, the fate of glyphosate under conditions of phosphorus deficiency or abundant soil should be explored to provide detailed information on glyphosate risk assessment in red soil.

    Soil aggregates; Glyphosate; Degradation kinetics; Soil phosphorus; Red soil

    10.13930/j.cnki.cjea.200686

    辛龍川, 吳文雪, 薛萐, 付格娟, 陳延華, 王學(xué)霞, 劉東生, 楊曉梅. 紅壤不同粒徑團聚體對草甘膦農(nóng)藥降解動力學(xué)的影響[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(中英文), 2021, 29(5): 910-921

    XIN L C, WU W X, XUE S, FU G J, CHEN Y H, WANG X X, LIU D S, YANG X M. Effects of aggregate size on kinetics of glyphosate degradation in red soil[J]. Chinese Journal of Eco-Agriculture, 2021, 29(5): 910-921

    X53

    * 國家自然科學(xué)基金項目(41877072)和陜西省自然科學(xué)基金項目(2019JQ-639)資助

    楊曉梅, 主要研究方向為農(nóng)田污染物遷移機理與模型模擬、污染物環(huán)境風(fēng)險評估、土壤侵蝕與污染、土壤質(zhì)量評估與可持續(xù)發(fā)展。E-mail: xiaomei.yang@nwafu.edu.cn

    辛龍川, 主要研究方向為土壤污染物與微生物互饋作用。E-mail: xinlongchuan@163.com

    2020-08-21

    2021-01-01

    * The study was supported by the National Natural Science Foundation of China (41877072) and the Natural Science Foundation of Shaanxi Province (2019JQ-639).

    , E-mail: xiaomei.yang@nwafu.edu.cn

    Aug. 21, 2020;

    Jan. 1, 2021

    猜你喜歡
    紅壤半衰期草甘膦
    初探熱脫附技術(shù)在有機污染紅壤修復(fù)的應(yīng)用
    昆鋼科技(2021年6期)2021-03-09 06:10:26
    基于語言學(xué)中文學(xué)術(shù)圖書的半衰期分析研究*
    法國擬2020年廢除草甘膦的使用
    基于引用半衰期的我國五官學(xué)期刊文獻老化研究
    長江叢刊(2016年33期)2016-12-12 05:31:06
    長期不同施肥對赤紅壤稻田區(qū)肥力的影響
    基于JCR?的國外臨床醫(yī)學(xué)學(xué)科半衰期
    草甘膦生產(chǎn)廢水治理技術(shù)探討
    草甘膦 三季度有望吹響反攻號角
    營銷界(2015年23期)2015-02-28 22:06:12
    長期施肥下紅壤旱地土壤CO2排放及碳平衡特征
    長期施肥對紅壤pH、作物產(chǎn)量及氮、磷、鉀養(yǎng)分吸收的影響
    久热久热在线精品观看| 在线观看免费高清a一片| 国产成人91sexporn| 美女大奶头黄色视频| 汤姆久久久久久久影院中文字幕| 国产黄色免费在线视频| av在线观看视频网站免费| 波多野结衣一区麻豆| 七月丁香在线播放| 国产欧美日韩综合在线一区二区| 久久久亚洲精品成人影院| 日韩一本色道免费dvd| 成人黄色视频免费在线看| 国产国语露脸激情在线看| 国产麻豆69| 欧美日韩一级在线毛片| 亚洲综合色惰| 爱豆传媒免费全集在线观看| 亚洲国产av新网站| 国产成人免费观看mmmm| 高清av免费在线| 2018国产大陆天天弄谢| 夫妻午夜视频| 大片电影免费在线观看免费| 亚洲精品国产色婷婷电影| 久久ye,这里只有精品| 在线观看美女被高潮喷水网站| 国产免费一区二区三区四区乱码| 日韩免费高清中文字幕av| 一本大道久久a久久精品| 99久国产av精品国产电影| 欧美日韩精品网址| 91精品国产国语对白视频| 免费不卡的大黄色大毛片视频在线观看| 人人妻人人澡人人看| 丰满乱子伦码专区| 丝袜美足系列| 99久久综合免费| 亚洲欧美一区二区三区久久| 久久久久国产精品人妻一区二区| 丁香六月天网| 伊人久久国产一区二区| 国产成人精品一,二区| 精品人妻一区二区三区麻豆| 久久久国产欧美日韩av| 亚洲男人天堂网一区| av.在线天堂| 国产 一区精品| 中文字幕人妻熟女乱码| 18禁裸乳无遮挡动漫免费视频| 欧美亚洲 丝袜 人妻 在线| 美女大奶头黄色视频| 青草久久国产| 啦啦啦在线观看免费高清www| 午夜免费鲁丝| 欧美人与性动交α欧美软件| 黄色 视频免费看| 久久婷婷青草| 美女高潮到喷水免费观看| 色94色欧美一区二区| 国产午夜精品一二区理论片| 在线观看免费视频网站a站| 精品少妇黑人巨大在线播放| 电影成人av| 哪个播放器可以免费观看大片| 可以免费在线观看a视频的电影网站 | 精品99又大又爽又粗少妇毛片| 久久久久久久大尺度免费视频| 麻豆乱淫一区二区| 蜜桃国产av成人99| 丁香六月天网| 欧美黄色片欧美黄色片| 免费观看在线日韩| 一级,二级,三级黄色视频| 性少妇av在线| 国产精品女同一区二区软件| 黄色怎么调成土黄色| 中文字幕av电影在线播放| 在线观看免费高清a一片| 99热全是精品| 亚洲国产精品国产精品| 欧美+日韩+精品| 美女脱内裤让男人舔精品视频| 午夜福利网站1000一区二区三区| 18在线观看网站| 2022亚洲国产成人精品| 只有这里有精品99| 男人操女人黄网站| 精品午夜福利在线看| 少妇猛男粗大的猛烈进出视频| 韩国av在线不卡| 亚洲欧美色中文字幕在线| av有码第一页| 只有这里有精品99| 又粗又硬又长又爽又黄的视频| 卡戴珊不雅视频在线播放| 老司机亚洲免费影院| 熟妇人妻不卡中文字幕| 成年女人在线观看亚洲视频| 久久精品人人爽人人爽视色| 亚洲国产精品一区二区三区在线| 亚洲国产av新网站| 成人漫画全彩无遮挡| 好男人视频免费观看在线| 男男h啪啪无遮挡| 精品久久久精品久久久| 在线观看www视频免费| 久久精品国产亚洲av天美| 国产野战对白在线观看| 亚洲精品一二三| 中文字幕人妻丝袜制服| 日韩伦理黄色片| 99久久中文字幕三级久久日本| 久久久久国产一级毛片高清牌| 丝袜美足系列| 亚洲国产精品成人久久小说| 国产日韩欧美在线精品| 人人妻人人澡人人爽人人夜夜| 精品福利永久在线观看| 日本午夜av视频| 国产高清国产精品国产三级| 美女国产高潮福利片在线看| 如何舔出高潮| 国产亚洲精品第一综合不卡| 中文字幕人妻丝袜一区二区 | 免费高清在线观看日韩| 黑丝袜美女国产一区| av免费观看日本| 成人漫画全彩无遮挡| 最近2019中文字幕mv第一页| 国产日韩欧美亚洲二区| 考比视频在线观看| 五月开心婷婷网| 亚洲av电影在线进入| 亚洲av欧美aⅴ国产| av片东京热男人的天堂| 欧美精品国产亚洲| 91在线精品国自产拍蜜月| 9热在线视频观看99| 天堂俺去俺来也www色官网| 中文乱码字字幕精品一区二区三区| 女人被躁到高潮嗷嗷叫费观| 一个人免费看片子| 日韩,欧美,国产一区二区三区| 国产av精品麻豆| 国产女主播在线喷水免费视频网站| 搡女人真爽免费视频火全软件| 高清av免费在线| 午夜老司机福利剧场| 婷婷成人精品国产| 精品99又大又爽又粗少妇毛片| 18+在线观看网站| 免费高清在线观看日韩| 午夜日本视频在线| 极品人妻少妇av视频| 综合色丁香网| 亚洲精品日本国产第一区| 国产爽快片一区二区三区| 十八禁网站网址无遮挡| av在线观看视频网站免费| 桃花免费在线播放| 美女xxoo啪啪120秒动态图| 免费黄色在线免费观看| 久久久精品免费免费高清| 香蕉丝袜av| 看非洲黑人一级黄片| 亚洲精品国产色婷婷电影| 亚洲国产欧美网| 亚洲欧美精品综合一区二区三区 | 一个人免费看片子| 亚洲第一青青草原| 国产男女超爽视频在线观看| 男女边摸边吃奶| 26uuu在线亚洲综合色| 国产亚洲午夜精品一区二区久久| 黄色怎么调成土黄色| 亚洲国产色片| 99久久人妻综合| 国产精品国产av在线观看| 黄频高清免费视频| 男女边摸边吃奶| 成人国产av品久久久| 欧美97在线视频| 午夜福利在线免费观看网站| 欧美精品一区二区免费开放| 国产精品女同一区二区软件| 亚洲欧美一区二区三区国产| 国产精品久久久av美女十八| 国产毛片在线视频| 天天躁夜夜躁狠狠久久av| 女人精品久久久久毛片| 亚洲色图综合在线观看| 丝袜在线中文字幕| 99re6热这里在线精品视频| 亚洲色图 男人天堂 中文字幕| 国产综合精华液| 在线 av 中文字幕| 亚洲欧美精品综合一区二区三区 | 1024香蕉在线观看| 亚洲人成电影观看| 精品第一国产精品| 日韩中文字幕欧美一区二区 | 性少妇av在线| 亚洲av电影在线观看一区二区三区| 国产精品久久久久成人av| 日日啪夜夜爽| 国产男女超爽视频在线观看| 亚洲欧洲国产日韩| 毛片一级片免费看久久久久| 新久久久久国产一级毛片| 亚洲经典国产精华液单| 两个人看的免费小视频| 亚洲国产看品久久| 久久这里只有精品19| 国产黄色视频一区二区在线观看| 99国产综合亚洲精品| 男女国产视频网站| 亚洲av在线观看美女高潮| 美女大奶头黄色视频| 在线观看免费高清a一片| 久久这里有精品视频免费| 国产欧美日韩综合在线一区二区| 天堂俺去俺来也www色官网| 国产有黄有色有爽视频| av国产精品久久久久影院| 男人添女人高潮全过程视频| 日日撸夜夜添| 18禁动态无遮挡网站| av网站免费在线观看视频| 人人妻人人添人人爽欧美一区卜| av视频免费观看在线观看| 少妇人妻 视频| 少妇人妻精品综合一区二区| 最近最新中文字幕大全免费视频 | 18+在线观看网站| 一边摸一边做爽爽视频免费| 国产精品香港三级国产av潘金莲 | 一级a爱视频在线免费观看| 免费黄频网站在线观看国产| 久久人妻熟女aⅴ| 一级毛片电影观看| 国产成人aa在线观看| 午夜福利乱码中文字幕| 亚洲av欧美aⅴ国产| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品偷伦视频观看了| 性少妇av在线| 欧美亚洲日本最大视频资源| 亚洲av国产av综合av卡| av线在线观看网站| 亚洲av国产av综合av卡| 久久久久精品性色| 色94色欧美一区二区| 成人18禁高潮啪啪吃奶动态图| 午夜福利视频在线观看免费| 精品久久久久久电影网| 18禁观看日本| 午夜av观看不卡| 欧美日韩精品网址| 狠狠婷婷综合久久久久久88av| 久久青草综合色| 精品一品国产午夜福利视频| 一本色道久久久久久精品综合| 超碰97精品在线观看| 午夜福利网站1000一区二区三区| 我要看黄色一级片免费的| 久久鲁丝午夜福利片| 免费av中文字幕在线| 免费黄频网站在线观看国产| 亚洲综合色惰| 乱人伦中国视频| 国产黄频视频在线观看| 丝袜脚勾引网站| 老司机影院毛片| 观看美女的网站| 一区二区日韩欧美中文字幕| 熟女少妇亚洲综合色aaa.| 精品人妻一区二区三区麻豆| 一区在线观看完整版| 亚洲欧美日韩另类电影网站| 男的添女的下面高潮视频| 咕卡用的链子| 久久久久久久久久久免费av| 黄色毛片三级朝国网站| 久久鲁丝午夜福利片| 国产高清国产精品国产三级| 欧美变态另类bdsm刘玥| 97在线视频观看| 午夜影院在线不卡| 爱豆传媒免费全集在线观看| 老女人水多毛片| 国产欧美日韩综合在线一区二区| 18禁动态无遮挡网站| 国产精品国产三级专区第一集| 成年人午夜在线观看视频| 免费高清在线观看日韩| 午夜日韩欧美国产| 亚洲成av片中文字幕在线观看 | 春色校园在线视频观看| 国产xxxxx性猛交| 国产成人精品在线电影| 久久毛片免费看一区二区三区| 国产精品一区二区在线观看99| 亚洲av成人精品一二三区| 少妇猛男粗大的猛烈进出视频| 成年女人在线观看亚洲视频| 国产精品无大码| 日日摸夜夜添夜夜爱| 久久国产精品男人的天堂亚洲| 亚洲av欧美aⅴ国产| 中文字幕av电影在线播放| 国产精品麻豆人妻色哟哟久久| 寂寞人妻少妇视频99o| 免费av中文字幕在线| 亚洲精品第二区| 老司机影院毛片| 又粗又硬又长又爽又黄的视频| 日韩中文字幕欧美一区二区 | 免费日韩欧美在线观看| 日韩欧美一区视频在线观看| 欧美日韩精品成人综合77777| 国产老妇伦熟女老妇高清| 亚洲,欧美,日韩| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久精品电影小说| 人妻一区二区av| 日韩中文字幕视频在线看片| 亚洲欧美精品综合一区二区三区 | freevideosex欧美| 一本色道久久久久久精品综合| 国产精品99久久99久久久不卡 | 满18在线观看网站| 日韩中文字幕视频在线看片| 亚洲欧美日韩另类电影网站| 亚洲视频免费观看视频| 性少妇av在线| 在线 av 中文字幕| 人妻一区二区av| 亚洲综合色惰| 精品少妇一区二区三区视频日本电影 | 国产亚洲最大av| 久久久精品区二区三区| av网站免费在线观看视频| 久久久国产一区二区| 日韩制服丝袜自拍偷拍| 亚洲国产精品一区二区三区在线| 国产精品欧美亚洲77777| 亚洲欧美中文字幕日韩二区| 午夜福利在线免费观看网站| 69精品国产乱码久久久| 午夜福利影视在线免费观看| 日韩欧美一区视频在线观看| 最新中文字幕久久久久| 亚洲欧美一区二区三区国产| 亚洲精品美女久久av网站| 日日撸夜夜添| 欧美精品亚洲一区二区| 国产成人精品福利久久| 中国国产av一级| 老熟女久久久| 亚洲精品日韩在线中文字幕| 国精品久久久久久国模美| 精品亚洲成国产av| 只有这里有精品99| 国产在视频线精品| 少妇被粗大的猛进出69影院| 免费不卡的大黄色大毛片视频在线观看| 另类精品久久| 欧美最新免费一区二区三区| 久久久国产一区二区| 国产一级毛片在线| 国产精品国产三级专区第一集| 久久韩国三级中文字幕| 精品一品国产午夜福利视频| 日韩中字成人| 一级毛片我不卡| 大香蕉久久网| 久久久国产一区二区| 欧美 日韩 精品 国产| 男的添女的下面高潮视频| 免费观看性生交大片5| 欧美精品高潮呻吟av久久| 欧美日韩视频精品一区| 成年人午夜在线观看视频| 国产一区二区激情短视频 | 大片免费播放器 马上看| 18+在线观看网站| 亚洲成国产人片在线观看| 一边摸一边做爽爽视频免费| 精品国产乱码久久久久久男人| 天美传媒精品一区二区| 成年女人在线观看亚洲视频| 少妇人妻 视频| 韩国高清视频一区二区三区| 亚洲av日韩在线播放| 一二三四中文在线观看免费高清| 国产午夜精品一二区理论片| 十分钟在线观看高清视频www| 美女视频免费永久观看网站| 亚洲人成电影观看| 亚洲国产最新在线播放| 黄片小视频在线播放| 国产男人的电影天堂91| 丝瓜视频免费看黄片| 国产麻豆69| 国产亚洲精品第一综合不卡| 亚洲综合色网址| 午夜福利视频在线观看免费| 欧美bdsm另类| 亚洲三区欧美一区| 一级,二级,三级黄色视频| 99九九在线精品视频| 一区福利在线观看| 极品人妻少妇av视频| 欧美日韩精品网址| 满18在线观看网站| 色播在线永久视频| 熟女av电影| 国产成人免费观看mmmm| 在线观看www视频免费| 午夜福利在线观看免费完整高清在| 日韩欧美精品免费久久| 欧美+日韩+精品| 一级毛片电影观看| 天天躁狠狠躁夜夜躁狠狠躁| 一区二区av电影网| 国产亚洲一区二区精品| 国产成人精品久久二区二区91 | 高清av免费在线| 免费少妇av软件| 99热国产这里只有精品6| 亚洲一区中文字幕在线| 国产欧美亚洲国产| 青春草亚洲视频在线观看| 久久久欧美国产精品| av又黄又爽大尺度在线免费看| 欧美成人午夜免费资源| 中文欧美无线码| 午夜激情av网站| 最近最新中文字幕大全免费视频 | 一区二区三区精品91| 国产黄色视频一区二区在线观看| 熟女电影av网| 精品卡一卡二卡四卡免费| 国产成人精品福利久久| 久久99热这里只频精品6学生| 午夜福利一区二区在线看| 国产精品秋霞免费鲁丝片| 欧美精品亚洲一区二区| 中文字幕另类日韩欧美亚洲嫩草| 国语对白做爰xxxⅹ性视频网站| 亚洲,欧美,日韩| 久久影院123| 秋霞伦理黄片| 成年女人毛片免费观看观看9 | 欧美+日韩+精品| 国产日韩欧美在线精品| 亚洲av国产av综合av卡| 女性被躁到高潮视频| 久久人妻熟女aⅴ| av国产久精品久网站免费入址| 亚洲在久久综合| 亚洲国产精品一区三区| 国产午夜精品一二区理论片| 2018国产大陆天天弄谢| 日本欧美视频一区| 中国三级夫妇交换| 黑人巨大精品欧美一区二区蜜桃| 日韩av在线免费看完整版不卡| 超碰成人久久| 亚洲av日韩在线播放| 久久青草综合色| av有码第一页| 久久精品国产综合久久久| 国产精品av久久久久免费| 亚洲婷婷狠狠爱综合网| 国产熟女欧美一区二区| 丰满少妇做爰视频| 午夜福利在线免费观看网站| 嫩草影院入口| 大陆偷拍与自拍| 美国免费a级毛片| 免费观看在线日韩| 在现免费观看毛片| 精品人妻偷拍中文字幕| 韩国av在线不卡| 国产精品无大码| 纯流量卡能插随身wifi吗| 亚洲精品aⅴ在线观看| 国产亚洲欧美精品永久| 亚洲av在线观看美女高潮| 99久久中文字幕三级久久日本| 成年av动漫网址| 性少妇av在线| 一个人免费看片子| 777米奇影视久久| 精品人妻在线不人妻| 欧美另类一区| 亚洲精品久久成人aⅴ小说| 免费久久久久久久精品成人欧美视频| 久久久精品免费免费高清| 精品第一国产精品| 国产成人一区二区在线| av一本久久久久| 高清视频免费观看一区二区| 久久 成人 亚洲| 国产高清不卡午夜福利| 亚洲国产最新在线播放| av.在线天堂| 亚洲欧美一区二区三区国产| av在线观看视频网站免费| 日韩一区二区三区影片| 中文字幕另类日韩欧美亚洲嫩草| 观看美女的网站| 精品人妻熟女毛片av久久网站| 最近中文字幕高清免费大全6| 国产欧美亚洲国产| 热re99久久国产66热| 久久国产精品大桥未久av| 亚洲综合精品二区| 国产免费福利视频在线观看| 国产 精品1| 视频在线观看一区二区三区| 热re99久久精品国产66热6| 黄片无遮挡物在线观看| 啦啦啦啦在线视频资源| 国产精品秋霞免费鲁丝片| tube8黄色片| 在线观看三级黄色| 国产色婷婷99| 欧美日韩亚洲国产一区二区在线观看 | av电影中文网址| 久久久久久久久久久免费av| 国产精品 欧美亚洲| 亚洲国产精品999| 一二三四中文在线观看免费高清| 国产精品av久久久久免费| 国产老妇伦熟女老妇高清| 亚洲精品自拍成人| 不卡av一区二区三区| 国产精品一区二区在线观看99| 啦啦啦中文免费视频观看日本| 久久精品人人爽人人爽视色| 国语对白做爰xxxⅹ性视频网站| 天堂中文最新版在线下载| 热re99久久国产66热| 日韩在线高清观看一区二区三区| 男女边吃奶边做爰视频| 18+在线观看网站| 欧美精品国产亚洲| av线在线观看网站| 美女国产视频在线观看| 亚洲第一av免费看| 欧美国产精品一级二级三级| 18禁国产床啪视频网站| 啦啦啦中文免费视频观看日本| 90打野战视频偷拍视频| 最近手机中文字幕大全| 成人毛片a级毛片在线播放| 欧美成人午夜精品| 国产黄色视频一区二区在线观看| 午夜免费男女啪啪视频观看| 丰满迷人的少妇在线观看| 亚洲成av片中文字幕在线观看 | 99热全是精品| 69精品国产乱码久久久| 国产日韩欧美亚洲二区| 中文精品一卡2卡3卡4更新| 18+在线观看网站| 9色porny在线观看| www.av在线官网国产| 日韩免费高清中文字幕av| 国产女主播在线喷水免费视频网站| 日韩中文字幕欧美一区二区 | 一个人免费看片子| 自拍欧美九色日韩亚洲蝌蚪91| 久久毛片免费看一区二区三区| 国产精品三级大全| 另类亚洲欧美激情| 美女xxoo啪啪120秒动态图| 国产精品.久久久| 99久久精品国产国产毛片| 视频在线观看一区二区三区| 亚洲国产成人一精品久久久| 99热网站在线观看| 男女无遮挡免费网站观看| 在线观看美女被高潮喷水网站| 最近中文字幕2019免费版| 成年人午夜在线观看视频| 亚洲熟女精品中文字幕| videos熟女内射| 国产人伦9x9x在线观看 | 好男人视频免费观看在线| 伊人亚洲综合成人网| 女人被躁到高潮嗷嗷叫费观| 18禁裸乳无遮挡动漫免费视频| 1024香蕉在线观看| 9191精品国产免费久久| 日韩人妻精品一区2区三区| 美女视频免费永久观看网站| 亚洲成国产人片在线观看| 久久av网站| 亚洲精品国产av成人精品| 欧美日韩成人在线一区二区| 亚洲少妇的诱惑av| 国产激情久久老熟女| 捣出白浆h1v1| 免费日韩欧美在线观看| 91aial.com中文字幕在线观看| 精品亚洲成国产av| 丰满迷人的少妇在线观看| 精品国产露脸久久av麻豆| 国产乱来视频区| 亚洲综合色惰| 欧美老熟妇乱子伦牲交|