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

    莠去津?qū)σ吧参锶郝湎峦寥牢⑸锕δ芏鄻有缘挠绊?/h1>
    2017-08-09 00:42:45肖能文李俊生
    環(huán)境科學研究 2017年8期
    關(guān)鍵詞:除草劑利用差異

    閆 冰, 齊 月, 付 剛, 賀 婧, 肖能文, 李俊生*

    1.中國環(huán)境科學研究院, 北京 100012 2.北京師范大學水科學研究院, 北京 100875 3.中國人民大學環(huán)境學院, 北京 100872

    莠去津?qū)σ吧参锶郝湎峦寥牢⑸锕δ芏鄻有缘挠绊?/p>

    閆 冰1,2, 齊 月1,2, 付 剛1,2, 賀 婧1,3, 肖能文1, 李俊生1*

    1.中國環(huán)境科學研究院, 北京 100012 2.北京師范大學水科學研究院, 北京 100875 3.中國人民大學環(huán)境學院, 北京 100872

    為探明亞致死劑量農(nóng)藥對土壤微生物多樣性影響,基于大田控制試驗,采用Biolog-ECO微平板技術(shù)分析噴施亞致死劑量的莠去津?qū)σ吧参锶郝湎峦寥牢⑸锶郝涔δ芏鄻有缘挠绊? 結(jié)果表明:①莠去津的噴施濃度對土壤微生物群落利用全部碳源的能力無顯著影響(P>0.05),但與噴施處理后第30天相比,處理后第60天的土壤微生物群落利用全部碳源的能力有所降低,其中處理Ⅰ〔噴施濃度(以質(zhì)量濃度計)為1 200 ghm2〕、處理Ⅲ(噴施濃度為300 ghm2)顯著下降(P<0.05),而CK(對照組)、處理Ⅱ(噴施濃度為600 ghm2)降低程度不顯著(P>0.05),并且其各處理間土壤微生物群落對全部碳源的利用能力大小為CK>處理Ⅲ>處理Ⅱ>處理Ⅰ. ②噴施處理后第30天土壤微生物群落對碳水化合物類碳源的利用能力最強,對羧酸類碳源的利用能力相對較弱;噴施處理后第60天土壤微生物群落對碳水化合物類碳源的利用能力降幅最大,處理Ⅰ、處理Ⅱ、處理Ⅲ分別下降了0.63、0.78、0.85、0.76. ③噴施除草劑后第30天處理Ⅱ的Shannon-Wiener指數(shù)和Simpson指數(shù)顯著低于其他處理(P<0.05),而第60天各處理間多樣性指數(shù)無顯著差異;與噴施處理后第30天相比,第60天各處理多樣性指數(shù)都有所下降,但沒有顯著差異(P>0.05). ④對土壤微生物群落碳源利用能力進行主成分分析結(jié)果顯示,提取的與土壤微生物碳源利用相關(guān)的主成分累積貢獻率為86.6%,主成分1能夠區(qū)分噴施處理后兩個時期的土壤微生物群落特征. ⑤冗余分析(RDA)表明,pH是土壤微生物利用碳源能力差異的重要影響因素之一. 研究顯示,除草劑莠去津的施用可在一定程度上降低野生植物群落下土壤微生物群落功能多樣性.

    莠去津; 土壤微生物; 微生物功能多樣性; Biolog-ECO微平板技術(shù)

    土壤微生物是土壤及整個生態(tài)系統(tǒng)的重要組成部分,是生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動中的重要驅(qū)動者,對土壤肥力的形成、污染物降解和土壤結(jié)構(gòu)保持等起著積極作用[1],并且在各種生物地球化學循環(huán)中扮演著不可替代的角色[2]. 土壤中含有數(shù)量驚人的土壤微生物,并且具有極其豐富的微生物多樣性[3]. 土壤微生物多樣性對維持土壤生態(tài)系統(tǒng)健康穩(wěn)定及可持續(xù)性至關(guān)重要[4- 5],是評價土壤質(zhì)量的重要生物學指標[6].

    除草劑是當今廣泛使用的一類農(nóng)藥,由于具有高效和適合機械化作業(yè)等特點,其產(chǎn)量和使用量都呈現(xiàn)出逐年增加的趨勢[7- 10]. 除草劑的大量施用在給農(nóng)業(yè)生產(chǎn)帶來效益的同時,也給生態(tài)環(huán)境帶來一系列的影響. 農(nóng)田中噴施的除草劑通過直接和間接的方式進入土壤,對土壤微生物數(shù)量、多樣性、生態(tài)特征等產(chǎn)生影響[11- 14].有研究[15]顯示,土壤微生物多樣性對植物生長及產(chǎn)量都有強烈的影響. 并且,有研究[16]認為土壤微生物在土壤中農(nóng)藥的降解上扮演著重要角色. 目前,化學農(nóng)藥的施用對土壤微生物及其活性的影響,已成為評價農(nóng)藥生態(tài)安全的重要指標之一[17- 18].

    農(nóng)田周邊野生植物群落的生物多樣性是農(nóng)田生態(tài)系統(tǒng)生物多樣性中的重要組成部分,其為農(nóng)業(yè)生產(chǎn)提供了重要的服務功能[19]. 野生植物多樣性促進作物害蟲的天敵和傳粉動物群落多樣性和密度的發(fā)展,還可以作為作物害蟲防治的植物誘餌,并促進農(nóng)田及周邊土壤的養(yǎng)分循環(huán)[20- 21]. 一直以來人們關(guān)注的主要是除草劑的除草效果[22- 23]及其對田內(nèi)土壤微生物的影響[12- 13],而對于農(nóng)田除草劑施用過程中飄散到周邊的亞致死劑量的除草劑對生物多樣性的影響研究較少. 亞致死劑量的除草劑可能會影響農(nóng)田周邊野生植物群落的生物多樣性,也可能對土壤微生物產(chǎn)生影響,因此,研究亞致死劑量的除草劑對土壤微生物影響至關(guān)重要.

    該研究基于室外大田控制試驗,采用Biolog-ECO微平板技術(shù)對不同亞致死劑量除草劑施用下土壤微生物群落功能多樣性變化進行研究,以期為評價亞致死劑量除草劑在土壤中的生態(tài)安全提供參考,為探索更為經(jīng)濟合理地使用除草劑提供基礎數(shù)據(jù),為減少環(huán)境污染提供參考.

    1 材料與方法

    1.1 研究區(qū)概況

    試驗地設置在中國環(huán)境科學研究院北京順義試驗基地內(nèi),該試驗基地位于北京市(39.4°N~41.6°N、115.7°E~117.4°E)順義區(qū)趙全營鎮(zhèn). 該地屬典型的北溫帶半濕潤大陸性季風氣候,夏季高溫多雨,冬季寒冷干燥,春、秋短促,全年無霜期180~200 d,年均氣溫13.4 ℃,年均降水量595 mm. 降水季節(jié)分配不均勻,全年80%的降水集中在6—8月.

    1.2 試驗設計

    試驗所用地塊于2011—2013年種植過玉米和油菜,并且未施用過任何農(nóng)藥,于2014年3月進行人工翻地,充分均勻混合使整個試驗地塊土壤處于勻質(zhì)狀態(tài). 試驗采用隨機區(qū)組設計,設置四組處理:CK(對照組,噴施等量的清水)、處理Ⅰ〔大田推薦劑量,噴施濃度(以質(zhì)量濃度計)為 1 200 ghm2,以ai(有效成分)計,下同〕、處理Ⅱ(大田推薦劑量的12,噴施濃度為600 ghm2)、處理Ⅲ(大田推薦劑量的14,噴施濃度為300 ghm2),每個處理三次重復,試驗共設置12個,小區(qū)大小為3 m×5 m,各小區(qū)間設置寬度為2 m 的隔離帶. 試驗于2014年6月20日進行除草劑噴施處理,噴施前用塑料薄膜圍建起高度為1.5 m的隔離網(wǎng),防止除草劑飄散導致各試驗處理間產(chǎn)生相互影響.

    試驗所用除草劑為莠去津(山東勝邦綠野化學有限公司生產(chǎn)),具有土壤處理作用兼有莖葉處理作用,能被植物根吸收,華北地區(qū)大田推薦劑量為 1 200 ghm2. 試驗所用噴霧器為上海生產(chǎn)的worth牌NS- 5型手動農(nóng)藥噴霧器,噴頭壓力設定為2×105Pa.

    各試驗小區(qū)內(nèi)的植物群落類型都為野生植物群落,主要物種組成為馬唐(Digitariasanguinalis)、反枝莧(Amaranthusretroflexus)、苘麻(AbutilontheophrastiMedicus)、藜(Chenopodiumalbum)等.

    1.3 樣品采集

    試驗在除草劑噴施處理后第30天、第60天進行樣品采集,于2014年7月20日進行第一次樣品采集,于2014年8月20日進行第二次樣品采集. 取樣時除去地表的凋落物層,垂直取0~15 cm的土壤,每個小區(qū)內(nèi)用五點取樣法取樣并進行充分混勻后記做一個樣品,除去植物殘體、根系、石塊等雜物,用無菌自封袋封裝再用冰盒帶回實驗室放入4 ℃冰箱中保存待用.

    1.4 研究方法

    采用電位法測定土壤pH,稱取10 g過2 mm篩的風干土,加入25 mL去CO2的超純水,土水比為1∶2.5,振蕩2 min后靜置30 min進行土壤pH測定,用元素分析儀(Costech ECS 4012,Italy)測量土壤w(TC)、w(TN).

    利用Biolog-ECO微平板法研究施用除草劑對土壤微生物功能多樣性的影響. 稱取相當于10 g干質(zhì)量的鮮土置于無菌錐形瓶中,加入90 mL濃度為0.85%的無菌生理鹽水,將混合液放置在恒溫搖床上振蕩30 min(200 rmin),靜止沉淀15 min,吸取上清液,采用十倍稀釋法將其用無菌生理鹽水稀釋至濃度為10-3. 將制備好的微生物菌懸浮液倒入無菌移液槽中,用八通道移液器將菌懸液接種于Biolog-ECO微平板中,每孔150 μL,全部接種工作在超凈工作臺上完成. 最后,將Biolog-ECO微平板置于28 ℃恒溫培養(yǎng)箱中避光培養(yǎng)7 d,每隔24 h用酶標儀讀取590 nm處的吸光值.

    根據(jù)讀取的吸光值數(shù)據(jù)計算Biolog-ECO微平板板孔的AWCD(average well color development,平均顏色變化率)〔見式(1)〕,用以表示土壤微生物群落利用全部碳源的能力[24]. 采用培養(yǎng)第72小時時的數(shù)據(jù)計算Shannon-Wiener指數(shù)(H′)〔見式(2)〕、Simpson指數(shù)(D)〔見式(3)〕、McIntos指數(shù)(U)〔見式(4)〕,豐富度指數(shù)(S)用有顏色變化的微孔數(shù)(微孔光密度值>0.15)表示[25- 27],用于反映土壤微生物群落的功能多樣性.

    (1)

    (2)

    (3)

    (4)

    式中:Ci為第i孔在590 nm下的光密度值;R為對照孔的光密度,若Ci-R≤0,記為0;n為Biolog-ECO微平板的碳源種類數(shù)目,n=31;Pi為第i孔相對光密度值(C-R)與整個平板相對光密度總和的比率;ni為第i孔的相對吸光度.

    利用SPSS 22.0對土壤微生物功能多樣性指數(shù)進行單因素方差分析(one-way ANOVA),采用最小顯著差異法(LSD)進行多重比較. 根據(jù)培養(yǎng)第72小時數(shù)據(jù),利用R軟件vegan包進行主成分分析(PCA)及冗余分析(RDA),并用ggplot2包進行繪圖.

    2 結(jié)果與分析

    2.1 土壤微生物利用全部碳源的代謝活性變化

    噴施處理后第30天、第60天土壤微生物群落AWCD隨培養(yǎng)時間的變化情況如圖1所示. 由圖1可見,隨著培養(yǎng)時間的延長,土壤微生物對全部碳源的代謝能力呈現(xiàn)出增加的趨勢. 噴施處理后第30天、培養(yǎng)第168小時各處理組之間無顯著差異(P>0.05),但處理Ⅲ的AWCD值相對較高,而其他三組處理的AWCD值非常接近. 噴施處理后第60天,各處理組間AWCD值大小為CK>處理Ⅲ>處理Ⅱ>處理Ⅰ,但單因素方差分析結(jié)果表明,培養(yǎng)第168小時各處理組間AWCD值無顯著差異(P>0.05). 表明噴施處理后第60天,不同處理對土壤微生物活性沒有顯著影響. 對比分析噴施處理后第30天和第60天兩個階段土壤微生物群落AWCD值(見表1)發(fā)現(xiàn),處理后第60天各處理土壤微生物群落AWCD值都有所下降,其中處理Ⅰ和處理Ⅲ顯著降低(P<0.05),而CK和處理Ⅱ沒有顯著降低(P>0.05).

    2.2 土壤微生物利用不同種類碳源的動力學特征

    Biolog-ECO微平板中的全部碳源可劃分為六大類,即多聚物類、糖類、氨基酸類、胺類、羧酸類、酚酸類. 培養(yǎng)第72小時時,不同處理下的土壤微生物對不同碳源利用情況如圖2所示. 由圖2可知,噴施處理后第30天土壤微生物對碳水化合物類碳源的利用

    圖1 不同處理土壤微生物群落平均顏色變化率(AWCD)隨培養(yǎng)時間的變化Fig.1 Variation of soil microbial communities average well color development(AWCD) with incubation time in different treatment

    Table 1 ANOVA results of AWCD for different treatments incubation 168 hours under 30 d and 60 d after spraying

    能力最強,對羧酸類碳源的利用能力最弱,并且對碳水化合物類碳源和羧酸類碳源的利用能力上存在顯著差異(P<0.05),說明利用碳水化合物類碳源的土壤微生物群落顯著多于利用羧酸類碳源的土壤微生物群落. 噴施處理后第30天,各處理間土壤微生物利用羧酸類和氨基酸類碳源能力都存在顯著差異(P<0.05),對其他類型碳源利用能力無顯著差異(P>0.05),說明不同處理對利用羧酸類和氨基酸類碳源的微生物影響較大. 噴施處理后第60天土壤微生物對多聚物、碳水化合物、氨基酸和胺類碳源的利用能力相對較強,對酚酸類和羧酸類碳源的利用能力相對較弱. 噴施處理后第60天,各處理間土壤微生物利用各類碳源均無顯著差異,說明噴施處理后第60天各處理造成的差異減弱. 對比噴施處理后第30天和第60天土壤微生物利用各類碳源的能力發(fā)現(xiàn),碳水化合物類碳源和氨基酸類碳源在兩個時期內(nèi)的AWCD值降幅最大.

    注:同一碳源類型中不同字母表示差異顯著(P≤0.05).圖2 不同處理土壤微生物對各類碳源利用情況Fig.2 Utilization of different carbon sources of soil microorganisms under the different treatments

    2.3 土壤微生物對碳源利用的多樣性指數(shù)分析

    由表2可知,噴施處理后后第30天土壤微生物群落Shannon-Wiener指數(shù)和Simpson指數(shù)在處理Ⅱ中顯著低于其他處理(P<0.05),不同處理對McIntosh指數(shù)和豐富度指數(shù)沒有顯著影響(P>0.05). 噴施除草劑處理后第60天,各處理對土壤微生物群落多樣性指數(shù)都無顯著影響(P>0.05). 同時,對噴施處理后后第30天和第60天各處理的多樣性指數(shù)進行t檢驗(數(shù)據(jù)未列出)發(fā)現(xiàn),與第30天相比,第60天的土壤微生物多樣性指數(shù)降低,但所有處理的第30天和第60天土壤微生物多樣性指數(shù)都無顯著差異(P>0.05). 說明噴施除草劑處理后,隨著時間的變化土壤微生物群落功能多樣性并不存在顯著變化的情況. 噴施處理后第30天,各處理對土壤微生物群落Shannon-Wiener指數(shù)和Simpson指數(shù)有顯著影響(P<0.05),而噴施處理后第60天這種影響不顯著,說明噴施除草劑對土壤微生物群落功能多樣性有影響,但60 d后這種影響開始減退或各處理間差異消失.

    2.4 土壤微生物群落利用碳源的主成分分析

    表2 不同處理土壤微生物群落功能多樣性指數(shù)

    注:數(shù)據(jù)為平均值±標準偏差. 同一行中不同字母表示差異顯著(P≤0.05).

    基于培養(yǎng)第72 小時時的吸光值數(shù)據(jù),對噴施處理后第30天和第60天的土壤微生物代謝碳源能力進行主成分分析(PCA),按照主成分提取原則[28]共提取了7個主成分(principal component,PC). 結(jié)果顯示,7個主成分累積方差貢獻率為86.6%,其中PC1(主成分1)~PC7(主成分7)的方差貢獻率分別為60.36%、7.48%、5.95%、3.79%、3.52%、3.10%和2.44%.

    不同處理土壤微生物碳源代謝的PCA排序圖如圖3所示. 由圖3可知,PC1和PC2共解釋了67.84%的變異. 在PC1上噴藥后各處理明顯分為兩簇,第30天的所有樣地聚為一簇,第60天的所有樣地聚為另一簇. PCA分析發(fā)現(xiàn),土壤微生物群落結(jié)構(gòu)在噴施處理后第30天和第60天間呈現(xiàn)出顯著差異(P<0.05). 在PC2上,噴施處理后第30天處理Ⅱ的土壤微生物群落結(jié)構(gòu)和其他處理存在顯著差異(P<0.05),噴施處理后第60天各處理間的土壤微生物群落結(jié)構(gòu)沒有顯示出顯著差異(P>0.05).

    注:30 d和60 d分別表示噴施處理后第30天和第60天,下同.圖3 不同處理土壤微生物群落利用碳源能力的主成分分析Fig.3 Principal component analysis of carbon sources utilization of soil microbial communities under different treatments

    通過對主成分相關(guān)性矩陣進行分析(見表3)可知,對PC1貢獻較大的碳源共12種,其中碳水化合物類碳源占66.67%,氨基酸類碳源占16.67,多聚物類碳源和酚酸類碳源均占8.33%,說明影響PC1的主要是碳水化合物類碳源,其次是氨基酸類碳源. 對PC2貢獻較大的碳源只有一種,為碳水化合物類碳源.

    表3 與PC1、PC2顯著相關(guān)的主要碳源種類

    2.4 土壤微生物群落功能冗余分析

    利用培養(yǎng)第72 小時的吸光值數(shù)據(jù)進行降趨對應分析(detrended correspondence analysis,DCA),結(jié)果發(fā)現(xiàn),排序軸中的最大值為1.1,因此選擇冗余分析(RDA)方法分析環(huán)境因子與土壤微生物碳源代謝能力之間的關(guān)系(見圖7). 結(jié)果顯示,軸1和軸2分別解釋碳源代謝能力總變異的38.32%和2.51%,pH、w(TC)和w(TN)與軸1和軸2都呈顯著正相關(guān)(P<0.05),其中pH與土壤微生物碳源代謝能力相關(guān)性最強(箭頭連線長度最長). 蒙特卡羅置換檢驗(Monte Carlo permutation test)顯示顯著性為0.02,則認為排序結(jié)果可以接受環(huán)境因子對碳源代謝能力變異的解釋量.

    注:TC—w(TC);TN—w(TN).圖4 土壤微生物群落功能特征與環(huán)境因子的冗余分析Fig.4 RDA for environmental factors and soil microbial function diversity

    3 討論

    土壤微生物群落利用全部碳源的能力隨著處理后時間的延長而降低,與噴施處理后第30天相比,噴施處理后第60天土壤微生物利用碳源的能力降低,處理Ⅰ和處理Ⅲ的土壤微生物群落AWCD都顯著下降,但同一時期內(nèi)各處理間土壤微生物群落利用碳源的能力沒有顯著差異. 說明莠去津在不同濃度梯度下對土壤微生物活性沒有產(chǎn)生顯著影響,而隨著時間增加除草劑對土壤微生物活性的影響呈現(xiàn)出增加的趨勢. 噴施處理后第30天各處理土壤微生物群落對碳源的利用能力表現(xiàn)為處理Ⅰ相對較高. 有研究認為,在土壤中有農(nóng)藥開始介入時,土壤微生物群落對進入到土壤中的農(nóng)藥具有一定的耐受性,并且一些微生物會利用農(nóng)藥為碳源進行分解作用,這會導致微生物進行繁殖[29]. 該研究發(fā)現(xiàn),莠去津?qū)ν寥牢⑸锶郝涞挠绊懺趪娛┨幚砗蟮?0天開始顯現(xiàn),表現(xiàn)為CK處理的土壤微生物群落利用碳源能力相對較大,其次是處理Ⅲ、處理Ⅱ、處理Ⅰ. 因此,噴施處理后第60天各噴施農(nóng)藥處理組相比,低濃度莠去津處理的土壤微生物群落活性相對較高.

    土壤微生物對不同碳源利用能力的差異在一定程度上指示土壤微生物群落的結(jié)構(gòu)組成[30]. 噴施處理后兩個時期(第30天、第60天)土壤微生物群落對羧酸類碳源的利用能力都較低,并且各處理下土壤微生物群落對碳水化合物類碳源的利用能力降幅最大,隨著時間的增加對利用碳水化合物類碳源的微生物的影響加劇. 研究[31- 32]發(fā)現(xiàn),干擾后不同時期偏好利用不同類型碳源的土壤微生物可能增加或減少. 不同處理間土壤微生物對某些類型碳源利用存在顯著差異,同時也對一些類型碳源的利用沒有顯著差異. 姚斌等[33]認為,土壤微生物對不同類型碳源的利用不同,可能是除草劑進入土壤后土壤微生物能利用的碳源種類發(fā)生了改變. 孫約兵等[34]研究發(fā)現(xiàn),施用一定濃度的硝磺草酮后,土壤微生物群落對Biolog-ECO微平板中全部碳源的代謝能力整體上呈現(xiàn)出升高的狀態(tài),但從土壤微生物對各類型的碳源利用能力上看,不同濃度的硝磺草酮處理間存在一定的差異.

    Shannon-Wiener指數(shù)、Simpson指數(shù)、McIntosh指數(shù)和豐富度指數(shù)分布用來反映土壤微生物的豐富度、最常見的物種優(yōu)勢度和均勻度[30,35]. 噴施除草劑后,隨著時間的增加各多樣性指數(shù)都有所降低,但是差異并未達到顯著水平. 楊永華等[36]研究認為,農(nóng)藥污染的土壤中土壤微生物群落Shannon指數(shù)、Simpson指數(shù)和McIntosh指數(shù)顯著降低,主要原因是由于受農(nóng)藥污染土壤微生物群落中利用碳源的微生物物種的均一程度受到顯著影響.

    噴施處理后第30天,土壤微生物群落Shannon-Wiener指數(shù)和Simpson指數(shù)在處理Ⅱ中顯著低于其他處理,而噴施處理后第60天,各處理間土壤微生物群落多樣性指數(shù)表現(xiàn)出無顯著差異. 分析認為,在莠去津施用初期,土壤微生物群落對不同濃度的處理響應不同,碳源利用的多樣性存在明顯差異,而隨著時間的推移,這種影響被微生物群落自身的修復調(diào)劑能力抵消. 除草劑進入土壤后可能造成土壤微生物群落功能發(fā)生變化,適量的農(nóng)藥處理將使得土壤中原有的微生物群落內(nèi)部種群間的競爭關(guān)系發(fā)生變化,原有的平衡遭到破壞,優(yōu)勢類群改變[37]. Ros等[38]研究發(fā)現(xiàn),高濃度莠去津刺激土壤細菌生長,這些都可能使微生物多樣性降低. Moreno等[39]通過研究不同濃度莠去津處理下的土壤微生物量碳和呼吸發(fā)現(xiàn),高濃度條件下土壤微生物量碳減少,同時土壤微生物呼吸增強. 但是筆者并沒有發(fā)現(xiàn)高濃度的莠去津能夠顯著減少土壤微生物群落的功能多樣性指數(shù). 鄭麗萍等[40]認為,低濃度有機氯會造成逆境刺激微生物代謝加強,導致土壤微生物對碳源的利用率較高,而高濃度有機氯污染可能對微生物的活性起極大的抑制作用,減少微生物多樣性. 農(nóng)藥作為土壤微生物群落的外來干擾因素,濃度較小時,其干擾被微生物群落自身強大的功能性冗余所抵消,而當其濃度增大時,可能使得某些微生物菌種的生長受到抑制,最終導致土壤微生物群落均勻性下降[41]. 對土壤微生物代謝特性進行RDA分析,表明pH是土壤微生物碳源利用能力差異的重要因子,推測噴施莠去津?qū)ν寥览砘再|(zhì)的影響是造成土壤微生物碳源利用能力差異的重要原因. 但筆者發(fā)現(xiàn),包括pH在內(nèi)的環(huán)境因子只能解釋變量的42.21%,因此除此之外還存在其他影響因素,需要進一步研究.

    該研究對亞致死劑量除草劑莠去津處理后的土壤微生物群落功能多樣性進行了研究,得到了一些初步結(jié)果. 雖然Biolog-ECO微平板技術(shù)能夠簡單、方便、快捷靈敏地反映出土壤微生物碳源代謝能力,進而評估土壤微生物群落碳源代謝功能多樣性. 但由于其自身存在的局限性,即微平板中只有31種碳源,導致其不能反映自然條件下土壤中的復雜環(huán)境,致使所得的研究結(jié)果也并不能夠完整地反映土壤中微生物群落的功能多樣性和復雜性[42]. 因此,在后續(xù)的研究中,還應結(jié)合高通量測序技術(shù),深入研究亞致死劑量除草劑對非靶標植物群落下的土壤微生物群落的影響,以期為闡釋除草劑對土壤微生物多樣性的影響提供理論依據(jù).

    4 結(jié)論

    a) 噴施不同濃度莠去津后對土壤微生物代謝活性存在一定程度的影響,但是差異不明顯. 全劑量的莠去津(處理Ⅰ)噴施處理后,對土壤微生物代謝活性的影響程度最大,亞致死劑量的影響程度相對較弱.

    b) 噴施除草劑莠去津后,隨著時間的增加對利用不同類型碳源的土壤微生物群落都存在一定的影響,其中對利用碳水化合物類碳源的微生物的影響最為強烈.

    c) 噴施莠去津后,土壤微生物群落的Shannon-Wiener指數(shù)、Simpson指數(shù)和McIntosh指數(shù)都呈下降趨勢,但沒有達到顯著水平(P>0.05). 除草劑的施用在一定程度上降低土壤微生物群落功能多樣性.

    d) 土壤理化因子是土壤微生物代謝能力差異的重要影響因素,該研究中pH、TN、TC都對土壤微生物群落代謝碳源能力有顯著影響,而pH的影響最強烈.

    [1] ZAK D R,HOLMES W E,WHITE D C,etal.Plant diversity,soil microbial communities,and ecosystem function:are there any links?[J].Ecology,2003,84(8):2042- 2050.

    [2] FALKOWSKI P G,FENCHEL T,DELONG E F.The microbial engines that drive Earth′s biogeochemical cycles[J].Science,2008,320(5879):1034- 1039.

    [3] GANS J,WOLINSKY M,DUNBAR J.Computational improvements reveal great bacterial diversity and high metal toxicity in soil[J].Science,2005,309(5739):1387- 1390.

    [4] BAUHUS J,KHANNA P.Carbon and nitrogen turnover in two acid forest soils of southeast Australia as affected by phosphorus addition and drying and rewetting cycles [J].Biology and Fertility of Soils,1994,17(3):212- 218.

    [5] WARDLE D.A comparative assessment of factors which influence microbial biomass carbon and nitrogen levels in soil[J].Biological reviews,1992,67(3):321- 358.

    [6] SUN Ruibo,ZHANG Xuexian,GUO Xisheng,etal.Bacterial diversity in soils subjected to long-term chemical fertilization can be more stably maintained with the addition of livestock manure than wheat straw[J].Soil Biology and Biochemistry,2015,88(9):9- 18.

    [7] GIANESSI L P,REIGNER N P.The value of herbicides in US crop production[J].Weed Technology,2007,21(2):559- 566.

    [8] DAVIS A S,HILL J D,CHASE C A,etal.Increasing cropping system diversity balances productivity,profitability and environmental health[J].PloS One,2012,7(10):e47149.

    [9] YOUNG B G,GIBSON D J,GAGE K L,etal.Agricultural weeds in glyphosate-resistant cropping systems in the United States[J].Weed Science,2013,61(1):85- 97.

    [10] 強勝,宋小玲,戴偉民.抗除草劑轉(zhuǎn)基因作物面臨的機遇與挑戰(zhàn)及其發(fā)展策略[J].農(nóng)業(yè)生物技術(shù)學報,2010,18(1):114- 125. QIANG Sheng,SONG Xiaoling,DAI Weimin The opportunity and challenge faced by transgenic herbicide-resistant crops and their development strategy[J].Journal of Agricultural Biotechnology,2010,18(1):114- 125.

    [11] 姚斌,徐建民,尚鶴,等.阿特拉津除草劑對土壤微生物生態(tài)特征的影響[J].水土保持學報,2005,19(3):46- 49. YAO Bin,XU Jianmin,SHANG He,etal.Study on ecological effect of atrazine on soil microbial activity[J].Journal of Soil and Water Conservation,2005,19(3):46- 49.

    [12] 趙蘭,黎華壽.四種除草劑對稻田土壤微生物類群的影響[J].農(nóng)業(yè)環(huán)境科學學報,2008,27(2):508- 514. ZHAO Lan,LI Huashou.The impacts of four typical herbicides on microbial population in paddy soil[J].Journal of Agro-Environm ent Science,2008,27(2):508- 514.

    [13] 王正貴,封超年,郭文喜,等.除草劑異丙隆對麥田土壤微生物數(shù)量及酶活性的影響[J].應用與環(huán)境生物學報,2010,16(5):688- 691. WANG Zhenggui,FENG Chaonian,GUO Wenxi,etal.Effects of isoproturon on soil microbial populations and enzyme activities[J].Chinese Journal of Applied and Environmental Biology,2010,16(5):688- 691.

    [14] 李鑫,張景云,張萌萌,等.化學除草劑不同施用方法對紫花苜蓿根際土壤微生物群落碳源利用的影響[J].草地學報,2014,22(1):57- 64. LI Xin,ZHANG Jingyun,ZHANG Mengmeng,etal.Effects of chemical herbicide at different spraying ways on the carbon source metabolism of soil microbial community in alfalfa rhizosphere[J].Acta Agrestia Sinica,2014,22(1):57- 64.

    [15] BONKOWSKI M,ROY J.Soil microbial diversity and soil functioning affect competition among grasses in experimental microcosms[J].Oecologia,2005,143(2):232- 240.

    [16] YU Y L,CHEN Y X,LUO Y M,etal.Rapid degradation of butachlor in wheat rhizosphere soil[J].Chemosphere,2003,50(6):771- 774.

    [17] MARX MC,WOOD M,JARVIS S.A microplate fluorimetric assay for the study of enzyme diversity in soils[J].Soil Biology and Biochemistry,2001,33(12):1633- 1640.

    [18] TRASAR C C,LEIROS M,SEOANE S,etal.Limitations of soil enzymes as indicators of soil pollution[J].Soil Biology and Biochemistry,2000,32(13):1867- 1875.

    [19] ISAACS R,TUELL J,FIEDLER A,etal.Maximizing arthropod-mediated ecosystem services in agricultural landscapes:the role of native plants [J].Frontiers in Ecology and the Environment,2009,7(4):196- 203.

    [20] MARSHALL E,BROWN V,BOATMAN N,etal.The role of weeds in supporting biological diversity within crop fields [J].Weed Research,2003,43(2):77- 89.

    [21] PETITE S,BOURSAULT A,LE Guilloux M,etal.Weeds in agricultural landscapes:a review[J].Agronomy for Sustainable Development,2010,31(2):309- 317.

    [22] 張盼盼,王君杰,陳凌,等.不同除草劑對糜子田雜草的防除效果[J].西北農(nóng)業(yè)學報,2013,22(10):208- 212. ZHANG Panpan,WANG Junjie,CHEN Ling,etal.Control effect of different herbicides on weed in broomcorn millet field[J].Acta Agriculturae Boreali-occidentalis Sinica,2013,22(10):208- 212.

    [23] DUKE S,DAYAN F,ROMAGNI J,etal.Natural products as sources of herbicides:current status and future trends[J].Weed Research,2000,40(1):99- 112.

    [24] O′DONNELL A G,SEASMAN M,MACRAE A,etal.Plants and fertilisers as drivers of change in microbial community structure and function in soils [J].Plant and Soil,2001,232(1/2):135- 145.

    [25] ZAK J C,WILLING M R,MOORHEAD D L,etal.Functional diversity of microbial communities:a quantitative approach[J].Soil Biology and Biochemistry,1994,26(9):1101- 1108.

    [26] SCHIPPER L,DEGENS B,SPARLING G,etal.Changes in microbial heterotrophic diversity along five plant successional sequences[J].Soil Biology and Biochemistry,2001,33(15):2093- 2103.

    [27] 林先貴.土壤微生物研究原理與方法[M].北京:高等教育出版社,2010.

    [28] 郝黎仁,樊元,郝哲歐.SPSS 使用統(tǒng)計分析[M].北京:中國水利水電出版社,2003.

    [29] CROUZET O,WISZNIOWSKI J,DONNADIEU F,etal.Dose-dependent effects of the herbicide mesotrione on soil cyanobacterial communities[J].Archives of Environmental Contamination and Toxicology,2013,64(1):23- 31.

    [30] 李春艷,劉標,韓正敏,等.轉(zhuǎn)Bt水稻土壤微生物多樣性對O3濃度升高的響應[J].中國環(huán)境科學,2014,34(11):2922- 2930. LI Chunyan,LI Biao,Han Zhengmin,etal.Responses to elevated O3of soil microbial diversity from transgenic Bt rice[J].China Environmental Science,2014,34(11):2922- 2930.

    [31] 曹成有,姚金冬,韓曉姝,等.科爾沁沙地小葉錦雞兒固沙群落土壤微生物功能多樣性[J].應用生態(tài)學報,2011,22(9):2309- 2315. CAO Chengyou,YAO Jindong,HAN Xiaoshu,etal.Soil microbes functional diversity in sand-fixingCaraganamicrophyllacommunities in HorqinSandy Land[J].Chinese Journal of Applied Ecology,2011,22(9):2309- 2315.

    [32] 李小容,韋金玉,陳云,等.海南島不同林齡的木麻黃林地土壤微生物的功能多樣性[J].植物生態(tài)學報,2014,38(6):608- 618. LI Xiaorong,WEI Jinyu,CHEN Yun,etal.Functional diversity of soil microorganisms inCasuarinaequisetifoliawoodlands of different stand ages in Hainan Island[J].Chinese Journal of Plant Ecology,2014,38(6):608- 618.

    [33] 姚斌,徐建民,尚鶴,等.甲磺隆污染土壤的微生物生態(tài)效應[J].農(nóng)業(yè)環(huán)境科學學報,2005,24(3):557- 561. YAO Bin,XUE Jianmin,SHANG He,etal.Ecological effect of metsulfuron-methyl on soil microbe[J].Journal of Agro-Envionment Science,2005,24(3):557- 561.

    [34] 孫約兵,王潤瓏,徐應明,等.除草劑硝磺草酮對土壤微生物生態(tài)效應研究[J].中國環(huán)境科學,2016,36(1):190- 196. SUN Yuebing,WANG Runlong,XU Yingming,etal.Ecological effects of herbicide mesotrione on soil microbial communities[J].China Environmental Science,2016,36(1):190- 196.

    [35] HU Junli,LIN Xiangui,WANG Junhua,etal.Microbial functional diversity,metabolic quotient,and invertase activity of a sandy loam soil as affected by long-term application of organic amendment and mineral fertilizer[J].Journal of Soils and Sediments,2011,11(2):271- 280.

    [36] 楊永華,姚健,華曉梅.農(nóng)藥污染對土壤微生物群落功能多樣性的影響[J].微生物學雜志,2000,20(2):23- 25. YANG Yonghua,YAO Jian,HUA Xiaomei.Effect of pesticide pollution against functional microbial diversity in soil[J].Journal of Microbiology,2000,20(2):23- 25.

    [37] JOHNSEN K,JACOBSEN C S,TORSVIK V,etal.Pesticide effects on bacterial diversity in agricultural soils:a review[J].Biology and Fertility of Soils,2001,33(6):443- 453.

    [38] ROS M,GOBERNA M,MORENO J,etal.Molecular and physiological bacterial diversity of a semi-arid soil contaminated with different levels of formulated atrazine[J].Applied Soil Ecology,2006,34(2):93- 102.

    [39] MORENO J L,ALIAGA A,NAVAARRO S,etal.Effects of atrazine on microbial activity in semiarid soil[J].Applied Soil Ecology,2007,35(1):120- 127.

    [40] 鄭麗萍,龍濤,林玉鎖,等.Biolog-ECO解析有機氯農(nóng)藥污染場地土壤微生物群落功能多樣性特征[J].應用與環(huán)境生物學報,2013,19(5):759- 765. ZHENG Liping,LONG Tao,LIN Yusuo,etal.Biolog-ECO analysis of microbial community functional diversity in organochlorine contaminated soil[J].Chinese Journal of Applied and Environmental Biology,2013,19(5):759- 765.

    [41] 張紅,呂永龍,辛曉云,等.殺蟲劑類POPs對土壤中微生物群落多樣性的影響[J].生態(tài)學報,2005,25(4):937- 942. ZHANG Hong,LV Yonglong,XIN Xiaoyun,etal.Effects of organochlorine pesticides on soil microbial community functional diversity[J].Acta Ecologica Sinica,2005,25(4):937- 942.

    [42] KONG W D,ZHU Y G,FU B J,etal.The veterinary antibiotic oxytetracycline and Cu influence functional diversity of the soil microbial community[J].Environmental Pollution,2006,143(1):129- 137.

    Effects of Atrazine on Soil Microbial Functional Diversity of Wild Weed Communities

    YAN Bing1,2, QI Yue1,2, FU Gang1,2, HE Jing1,3, XIAO Nengwen1, LI Junsheng1*

    1.Chinese Research Academy of Environmental Sciences, Beijing 100012, China 2.College of Water Sciences, Beijing Normal University, Beijing 100875, China 3.School of Environment and Natural Resources, Renmin University of China, Beijing 100872, China

    To evaluate the impacts of atrazine on soil microbial communities under wild weed communities, the functional diversity of soil microbial communities in different treatments was investigated with Biolog-ECO micro-plates. The results showed that spraying concentrations of atrazine had no significant influence on soil microbial carbon utilization (P>0.05). Compared with the 30thday after spraying, the soil microbial carbon utilization ability was reduced in the 60thday. The I and III treatments were reduced significantly (P<0.05), while the CK and II treatments had no significant reduction (P>0.05). After spraying herbicide, the order of the ability of soil microbial carbon utilization after 60 days was CK>III>II>I. The utilization abilities of different types of carbon sources were different after spraying herbicide. 30 days after spraying, the soil microbial communities mostly used carbohydrates as carbon source, and used Carboxylic acid least. 60 days after spraying, the ability to use carbohydrates as carbon source was largely declined, with treatments declined by 0.63, 0.78, 0.85, 0.76 respectively. 30 days after spraying herbicide, the Shannon-Wiener and Simpson indices of soil microbial community in treatment II were significant less than other treatments (P<0.05), but the diversity indices showed no significant differences among treatments 30 days after spraying herbicide. Compared with the 30thday after spraying, the diversity index was decreased in the 60thday after spraying, but the differences were not significant (P>0.05). Principal component analysis demonstrated that the principal components extracted from 31 carbon sources were related to soil microbial carbon source utilization; the principal components accounted for 86.6% of the total variables variance. Principal component 1 could separate characteristics of soil microbial community after spraying. Redundancy Analysis showed that pH was the main factor that was most highly correlated with variance of soil microbial carbon utilization. The functional diversity of soil microbial communities will be changed by spraying the herbicide atrazine.

    atrazine; soil microorganisms; functional diversity; Biolog-ECO micro-plate

    2017-02-14

    2017-04-18

    科技部重點研發(fā)計劃項目(2016YFC1201100);北京市財政資金項目(PXM2014- 178218- 000005)

    閆冰(1986-),男,黑龍江綏化人,yanbing718@126.com.

    *責任作者,李俊生(1968-),男,安徽巢湖人,研究員,博士,博導,主要從事生物多樣性保護、自然保護區(qū)管理、生物安全評估以及氣候變化影響評價研究,lijsh@craes.org.cn

    X172

    1001- 6929(2017)08- 1246- 09

    A

    10.13198j.issn.1001- 6929.2017.02.54

    閆冰,齊月,付剛,等.莠去津?qū)σ吧参锶郝湎峦寥牢⑸锕δ芏鄻有缘挠绊慬J].環(huán)境科學研究,2017,30(8):1246- 1254.

    YAN Bing,QI Yue,FU Gang,etal.Effects of atrazine on soil microbial functional diversity of wild weed communities [J].Research of Environmental Sciences,2017,30(8):1246- 1254.

    猜你喜歡
    除草劑利用差異
    利用min{a,b}的積分表示解決一類絕對值不等式
    相似與差異
    音樂探索(2022年2期)2022-05-30 21:01:37
    利用一半進行移多補少
    找句子差異
    封閉式除草劑什么時間噴最合適
    利用數(shù)的分解來思考
    Roommate is necessary when far away from home
    生物為什么會有差異?
    如何正確選擇使用農(nóng)藥及除草劑
    玉米田除草劑的那些事
    營銷界(2015年23期)2015-02-28 22:06:18

    精品久久久久久电影网| 别揉我奶头 嗯啊视频| 国产精品久久久久久av不卡| 久久久久久国产a免费观看| 国产午夜精品久久久久久一区二区三区| 欧美最新免费一区二区三区| 久久久久久久久久成人| 中文字幕久久专区| 内地一区二区视频在线| 精品99又大又爽又粗少妇毛片| 麻豆乱淫一区二区| 97精品久久久久久久久久精品| 国产精品久久久久久av不卡| 简卡轻食公司| 国产人妻一区二区三区在| 亚洲人与动物交配视频| a级毛色黄片| 最近的中文字幕免费完整| 亚洲欧美清纯卡通| 国产免费视频播放在线视频 | 夜夜爽夜夜爽视频| 亚洲av成人精品一区久久| 亚洲图色成人| 伊人久久国产一区二区| 久久久久久久国产电影| www.色视频.com| 免费观看a级毛片全部| 国产精品无大码| 国产极品天堂在线| 在线a可以看的网站| 黄色配什么色好看| .国产精品久久| 国产有黄有色有爽视频| 精品国产一区二区三区久久久樱花 | 欧美激情久久久久久爽电影| 国产毛片a区久久久久| 日本熟妇午夜| 欧美日韩综合久久久久久| 大陆偷拍与自拍| 亚洲色图av天堂| 亚洲电影在线观看av| 亚洲性久久影院| 国产成人freesex在线| 岛国毛片在线播放| 精品国产露脸久久av麻豆 | 亚洲av免费在线观看| 一级a做视频免费观看| 国产高清有码在线观看视频| 日韩av在线大香蕉| 亚洲高清免费不卡视频| 人妻少妇偷人精品九色| 五月天丁香电影| 国产69精品久久久久777片| 国产精品蜜桃在线观看| 亚洲在线观看片| 亚洲无线观看免费| 久久久久久久久久成人| 国产伦理片在线播放av一区| 内射极品少妇av片p| 国产成人精品福利久久| 精品不卡国产一区二区三区| 人人妻人人澡人人爽人人夜夜 | 国内揄拍国产精品人妻在线| 97超视频在线观看视频| 精品一区在线观看国产| 超碰av人人做人人爽久久| 六月丁香七月| 午夜精品在线福利| 成年人午夜在线观看视频 | 黄色配什么色好看| 人妻制服诱惑在线中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 国产成人a∨麻豆精品| 一区二区三区乱码不卡18| 午夜福利在线观看吧| 国产黄色视频一区二区在线观看| 亚洲精品456在线播放app| 午夜久久久久精精品| 日韩一区二区视频免费看| 国产午夜精品一二区理论片| 日韩,欧美,国产一区二区三区| 亚洲欧美一区二区三区黑人 | 少妇的逼水好多| 日日摸夜夜添夜夜添av毛片| 一个人看的www免费观看视频| 国产永久视频网站| 少妇丰满av| 日韩三级伦理在线观看| 听说在线观看完整版免费高清| 日韩欧美精品免费久久| 少妇熟女欧美另类| 亚洲精品一区蜜桃| 国产欧美日韩精品一区二区| 少妇熟女aⅴ在线视频| 人人妻人人澡人人爽人人夜夜 | 97在线视频观看| 日韩av在线大香蕉| 能在线免费观看的黄片| 亚洲不卡免费看| 秋霞伦理黄片| 少妇熟女aⅴ在线视频| 91午夜精品亚洲一区二区三区| 亚洲美女视频黄频| 欧美另类一区| av国产久精品久网站免费入址| 成年人午夜在线观看视频 | 成人高潮视频无遮挡免费网站| 免费观看av网站的网址| 97人妻精品一区二区三区麻豆| 久久人人爽人人爽人人片va| 亚洲精品久久午夜乱码| 亚洲四区av| 国产精品久久久久久久久免| 国产一区有黄有色的免费视频 | 国产精品综合久久久久久久免费| 国产精品福利在线免费观看| eeuss影院久久| 亚洲内射少妇av| videossex国产| 国产精品国产三级国产专区5o| 精品久久久噜噜| 欧美不卡视频在线免费观看| 欧美另类一区| 日韩强制内射视频| 婷婷色综合www| 干丝袜人妻中文字幕| 国产免费又黄又爽又色| 欧美激情国产日韩精品一区| 身体一侧抽搐| 国产一区二区在线观看日韩| 夫妻午夜视频| 在线观看免费高清a一片| 国产成人aa在线观看| 欧美高清性xxxxhd video| 黄色欧美视频在线观看| 亚洲精品一二三| 寂寞人妻少妇视频99o| 亚洲成人中文字幕在线播放| 亚洲国产欧美人成| 免费观看无遮挡的男女| 国产精品久久久久久精品电影小说 | 亚洲欧美日韩东京热| 久久久久久九九精品二区国产| 精品国产三级普通话版| 视频中文字幕在线观看| 国产人妻一区二区三区在| 国产综合懂色| 亚洲精品一二三| 边亲边吃奶的免费视频| 七月丁香在线播放| 日日干狠狠操夜夜爽| 少妇人妻精品综合一区二区| 国产乱人偷精品视频| 免费看光身美女| freevideosex欧美| 中文精品一卡2卡3卡4更新| 好男人视频免费观看在线| 精品人妻熟女av久视频| 国产中年淑女户外野战色| 天堂影院成人在线观看| 伊人久久精品亚洲午夜| 男女啪啪激烈高潮av片| 亚洲一级一片aⅴ在线观看| 69av精品久久久久久| 三级国产精品片| 国产成人精品一,二区| 亚洲欧美日韩东京热| 男女啪啪激烈高潮av片| 亚洲一级一片aⅴ在线观看| 一夜夜www| 天堂俺去俺来也www色官网 | 国产精品综合久久久久久久免费| 久久久久久久久大av| 免费观看性生交大片5| 国产中年淑女户外野战色| 天堂网av新在线| 夫妻性生交免费视频一级片| 黄片无遮挡物在线观看| 久久久久久久久久黄片| 国产黄片视频在线免费观看| 国产麻豆成人av免费视频| 一级爰片在线观看| 国内精品一区二区在线观看| 91午夜精品亚洲一区二区三区| 国产一区亚洲一区在线观看| 久久韩国三级中文字幕| 啦啦啦啦在线视频资源| 久久这里有精品视频免费| 国产熟女欧美一区二区| 久久99蜜桃精品久久| 男女边吃奶边做爰视频| 2021少妇久久久久久久久久久| 精品久久久久久成人av| 国产精品人妻久久久久久| 国产免费福利视频在线观看| 久久久久久伊人网av| 一夜夜www| 人妻系列 视频| 联通29元200g的流量卡| 国产黄片美女视频| av国产免费在线观看| 99九九线精品视频在线观看视频| 国产老妇伦熟女老妇高清| 天堂影院成人在线观看| 成人高潮视频无遮挡免费网站| 日产精品乱码卡一卡2卡三| 国产国拍精品亚洲av在线观看| 卡戴珊不雅视频在线播放| 国产精品一二三区在线看| 国产探花在线观看一区二区| 插逼视频在线观看| 亚洲欧洲国产日韩| 日韩成人av中文字幕在线观看| 久久精品夜色国产| 免费少妇av软件| 免费观看av网站的网址| 精品国产三级普通话版| 久久97久久精品| 麻豆av噜噜一区二区三区| 尾随美女入室| 神马国产精品三级电影在线观看| 久久人人爽人人爽人人片va| 男插女下体视频免费在线播放| 中文字幕制服av| 免费av观看视频| 神马国产精品三级电影在线观看| 久久人人爽人人爽人人片va| 亚洲欧美清纯卡通| 国产在线一区二区三区精| 国内精品一区二区在线观看| 黄片无遮挡物在线观看| 久久久精品欧美日韩精品| 乱系列少妇在线播放| 精品人妻一区二区三区麻豆| 女人久久www免费人成看片| 美女主播在线视频| 啦啦啦中文免费视频观看日本| 国产视频内射| 免费av观看视频| 男女边摸边吃奶| 能在线免费看毛片的网站| 亚洲人与动物交配视频| 亚洲丝袜综合中文字幕| 精品一区二区三区人妻视频| 亚洲怡红院男人天堂| 日韩av在线大香蕉| 一级av片app| 日韩成人伦理影院| 婷婷色麻豆天堂久久| 一级片'在线观看视频| 亚洲精品国产av成人精品| 成人亚洲精品av一区二区| 国产精品99久久久久久久久| 久久99蜜桃精品久久| 久久久精品欧美日韩精品| 亚洲精品第二区| 午夜精品一区二区三区免费看| 亚洲av二区三区四区| 精华霜和精华液先用哪个| 欧美日韩视频高清一区二区三区二| 韩国高清视频一区二区三区| 97超碰精品成人国产| 婷婷色麻豆天堂久久| 日韩强制内射视频| 99热全是精品| 日日摸夜夜添夜夜添av毛片| 国产亚洲5aaaaa淫片| 亚洲最大成人av| 国产视频内射| 熟妇人妻久久中文字幕3abv| 哪个播放器可以免费观看大片| 亚洲欧美日韩卡通动漫| 直男gayav资源| 一级毛片黄色毛片免费观看视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 人人妻人人澡欧美一区二区| 欧美性猛交╳xxx乱大交人| 青春草国产在线视频| 亚洲,欧美,日韩| av免费观看日本| 欧美成人a在线观看| 直男gayav资源| 伊人久久精品亚洲午夜| 91午夜精品亚洲一区二区三区| 成人特级av手机在线观看| 青春草国产在线视频| 免费看光身美女| 在线观看一区二区三区| 97热精品久久久久久| 高清在线视频一区二区三区| 最近中文字幕2019免费版| 丰满人妻一区二区三区视频av| 久久国内精品自在自线图片| 国产成人a区在线观看| 成人毛片a级毛片在线播放| 国产av不卡久久| 狂野欧美白嫩少妇大欣赏| 国产亚洲精品av在线| 2021天堂中文幕一二区在线观| 视频中文字幕在线观看| 亚洲精品456在线播放app| 大香蕉久久网| 亚洲国产精品专区欧美| 日本一二三区视频观看| 免费人成在线观看视频色| 久久久亚洲精品成人影院| 久久99热这里只频精品6学生| 99久久人妻综合| 成年女人看的毛片在线观看| 少妇的逼好多水| 成人性生交大片免费视频hd| 国产伦理片在线播放av一区| 五月伊人婷婷丁香| 色播亚洲综合网| 久久精品国产亚洲av天美| 国产精品一区二区在线观看99 | 一级毛片我不卡| 国产色爽女视频免费观看| 纵有疾风起免费观看全集完整版 | 国产成人aa在线观看| 国产精品不卡视频一区二区| 欧美日韩在线观看h| 69av精品久久久久久| 日日啪夜夜爽| av网站免费在线观看视频 | 乱系列少妇在线播放| 尾随美女入室| 菩萨蛮人人尽说江南好唐韦庄| 国产在视频线精品| 久久这里有精品视频免费| 亚洲国产精品sss在线观看| videossex国产| 成年女人在线观看亚洲视频 | 人人妻人人澡欧美一区二区| 亚洲aⅴ乱码一区二区在线播放| 真实男女啪啪啪动态图| 精品人妻偷拍中文字幕| 国产综合精华液| 极品少妇高潮喷水抽搐| 午夜福利网站1000一区二区三区| 综合色av麻豆| 精品99又大又爽又粗少妇毛片| 日本-黄色视频高清免费观看| 久久精品综合一区二区三区| 中文字幕亚洲精品专区| 高清午夜精品一区二区三区| 亚洲欧洲日产国产| 国产单亲对白刺激| 你懂的网址亚洲精品在线观看| 国产精品三级大全| 久热久热在线精品观看| 欧美区成人在线视频| av一本久久久久| av线在线观看网站| 国内精品美女久久久久久| 欧美一区二区亚洲| 99久久人妻综合| 国产乱人视频| 国产黄色小视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 久久热精品热| 日日干狠狠操夜夜爽| 又黄又爽又刺激的免费视频.| 亚洲精品日韩在线中文字幕| 久久6这里有精品| 成人综合一区亚洲| 亚洲精品久久久久久婷婷小说| 一区二区三区免费毛片| 成人午夜精彩视频在线观看| 91aial.com中文字幕在线观看| 日韩欧美精品免费久久| 18+在线观看网站| 日韩制服骚丝袜av| 亚洲人与动物交配视频| 国产老妇女一区| 嫩草影院精品99| 日本三级黄在线观看| 国产男女超爽视频在线观看| 亚洲国产欧美在线一区| 欧美3d第一页| 欧美区成人在线视频| av在线天堂中文字幕| 极品少妇高潮喷水抽搐| 91午夜精品亚洲一区二区三区| 少妇猛男粗大的猛烈进出视频 | 18+在线观看网站| 色吧在线观看| 色视频www国产| 国产免费一级a男人的天堂| 久久综合国产亚洲精品| 中文字幕人妻熟人妻熟丝袜美| 只有这里有精品99| 亚洲精品国产av成人精品| 亚洲国产日韩欧美精品在线观看| www.色视频.com| 美女高潮的动态| 在线天堂最新版资源| 免费大片18禁| 亚洲av成人精品一区久久| 亚洲真实伦在线观看| 丰满少妇做爰视频| 日日摸夜夜添夜夜添av毛片| 久久久久久久久久黄片| 啦啦啦啦在线视频资源| 插逼视频在线观看| 久久久精品欧美日韩精品| 国产久久久一区二区三区| 黄片无遮挡物在线观看| 禁无遮挡网站| 久久精品国产亚洲av涩爱| 亚洲熟妇中文字幕五十中出| 国产乱来视频区| 高清午夜精品一区二区三区| 综合色丁香网| 久久久精品免费免费高清| 深夜a级毛片| 在线 av 中文字幕| 国内精品一区二区在线观看| 97在线视频观看| 观看免费一级毛片| 日韩欧美精品v在线| 大香蕉97超碰在线| 看免费成人av毛片| 国产精品日韩av在线免费观看| 精品久久久久久久久久久久久| eeuss影院久久| videossex国产| 午夜激情福利司机影院| 99视频精品全部免费 在线| 91精品伊人久久大香线蕉| 久久精品国产亚洲av涩爱| 久久久久久久午夜电影| 国产大屁股一区二区在线视频| 天堂av国产一区二区熟女人妻| 免费av观看视频| 亚洲精品亚洲一区二区| 在线观看免费高清a一片| 国产色婷婷99| 美女大奶头视频| 成人av在线播放网站| 黄色配什么色好看| 久久久国产一区二区| 欧美成人精品欧美一级黄| 日韩欧美 国产精品| 成人无遮挡网站| 国产伦在线观看视频一区| av又黄又爽大尺度在线免费看| 欧美日韩亚洲高清精品| 天天躁夜夜躁狠狠久久av| 国产在视频线精品| 亚洲精品一二三| 免费电影在线观看免费观看| 国产精品伦人一区二区| 成人毛片60女人毛片免费| 国产精品久久久久久av不卡| 日日摸夜夜添夜夜添av毛片| 亚洲精品乱久久久久久| 中文字幕亚洲精品专区| 亚洲av中文av极速乱| 久久国内精品自在自线图片| 白带黄色成豆腐渣| 国产精品精品国产色婷婷| 国产色婷婷99| 搡女人真爽免费视频火全软件| 国产亚洲一区二区精品| av天堂中文字幕网| 男人舔女人下体高潮全视频| 水蜜桃什么品种好| 精品不卡国产一区二区三区| 久久韩国三级中文字幕| 91精品一卡2卡3卡4卡| 高清av免费在线| 精品久久久久久久久久久久久| 精品人妻视频免费看| 99热6这里只有精品| 免费观看av网站的网址| 国产成人a∨麻豆精品| 国产午夜精品论理片| 国产一区二区亚洲精品在线观看| 亚洲在线自拍视频| 99久久精品热视频| 成人亚洲精品av一区二区| 亚洲乱码一区二区免费版| 国产精品久久久久久精品电影小说 | 中文字幕制服av| 天天躁夜夜躁狠狠久久av| 女人十人毛片免费观看3o分钟| 婷婷色av中文字幕| 午夜福利高清视频| 99久国产av精品国产电影| 国产精品1区2区在线观看.| 69av精品久久久久久| 91精品一卡2卡3卡4卡| 又粗又硬又长又爽又黄的视频| 亚洲国产精品国产精品| 亚洲欧美精品专区久久| 日韩 亚洲 欧美在线| 成人毛片60女人毛片免费| 99久国产av精品国产电影| 国产老妇女一区| av在线天堂中文字幕| 免费看不卡的av| 丰满乱子伦码专区| 丝袜喷水一区| 久久这里只有精品中国| 美女高潮的动态| 久久精品人妻少妇| 久久久久精品性色| 国产免费福利视频在线观看| 日韩欧美精品免费久久| 国产黄片美女视频| 午夜福利视频精品| 人妻制服诱惑在线中文字幕| 久久久精品欧美日韩精品| 成年女人在线观看亚洲视频 | 床上黄色一级片| 秋霞伦理黄片| 国产国拍精品亚洲av在线观看| 在线观看人妻少妇| 搡女人真爽免费视频火全软件| 国产精品久久久久久久久免| 国产午夜精品久久久久久一区二区三区| 尾随美女入室| 一级a做视频免费观看| 久久久精品免费免费高清| 一级片'在线观看视频| 黄色欧美视频在线观看| 亚洲av成人av| 国产日韩欧美在线精品| 亚洲av福利一区| 成人美女网站在线观看视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产综合精华液| 国产伦在线观看视频一区| 一级毛片久久久久久久久女| 尾随美女入室| 99九九线精品视频在线观看视频| 麻豆成人av视频| 成人漫画全彩无遮挡| 国产又色又爽无遮挡免| 国国产精品蜜臀av免费| 人人妻人人看人人澡| av一本久久久久| 亚洲精品一二三| 听说在线观看完整版免费高清| 国产伦一二天堂av在线观看| 一级毛片久久久久久久久女| 国产伦一二天堂av在线观看| 欧美性感艳星| 亚洲国产日韩欧美精品在线观看| 成人漫画全彩无遮挡| 高清视频免费观看一区二区 | 日日摸夜夜添夜夜添av毛片| 乱人视频在线观看| 国产伦理片在线播放av一区| 18禁在线播放成人免费| 麻豆成人av视频| 久久97久久精品| 国产精品一区二区三区四区免费观看| 一个人观看的视频www高清免费观看| 免费无遮挡裸体视频| 天天躁日日操中文字幕| 女人久久www免费人成看片| 天天躁日日操中文字幕| 国产男女超爽视频在线观看| 欧美日韩国产mv在线观看视频 | 国产麻豆成人av免费视频| 亚洲av中文av极速乱| 少妇人妻精品综合一区二区| 国产爱豆传媒在线观看| 成人午夜高清在线视频| 亚洲丝袜综合中文字幕| 日韩不卡一区二区三区视频在线| 久久韩国三级中文字幕| 国产免费又黄又爽又色| 五月玫瑰六月丁香| 午夜精品国产一区二区电影 | 色5月婷婷丁香| 少妇裸体淫交视频免费看高清| 欧美激情国产日韩精品一区| 日韩强制内射视频| 国产成人精品一,二区| 日韩欧美精品v在线| 五月天丁香电影| 国内精品一区二区在线观看| 搡老乐熟女国产| 久久久久久久大尺度免费视频| 国产精品无大码| 日本黄色片子视频| 2022亚洲国产成人精品| 51国产日韩欧美| 美女大奶头视频| 欧美成人午夜免费资源| 99热网站在线观看| 日韩电影二区| 激情 狠狠 欧美| 黄色配什么色好看| 国产精品久久久久久av不卡| 永久网站在线| 少妇熟女欧美另类| 国产一区二区在线观看日韩| 中文字幕免费在线视频6| 色吧在线观看| 久久久久久久久中文| 亚洲精品久久午夜乱码| 大陆偷拍与自拍| 国产成人aa在线观看| 少妇人妻一区二区三区视频| 别揉我奶头 嗯啊视频| 黑人高潮一二区| 免费看不卡的av| 国产高清有码在线观看视频| 在线天堂最新版资源| 日韩 亚洲 欧美在线| 深夜a级毛片|