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

    石油烴污染對隴東黃土高原土壤生物學(xué)及非生物學(xué)特性的影響

    2017-03-27 10:14:19王金成井明博段春燕王國鋒石國璽周天林
    水土保持通報 2017年1期
    關(guān)鍵詞:隴東土壤環(huán)境樣地

    王金成, 井明博, 段春燕, 王國鋒, 石國璽,4, 周天林

    (1.甘肅省高校隴東生物資源保護與利用省級重點實驗室, 甘肅 慶陽 745000; 2.隴東學(xué)院 生命科學(xué)與技術(shù)學(xué)院,甘肅 慶陽 745000; 3.陜西師范大學(xué) 生命科學(xué)學(xué)院, 陜西 西安 710062; 4.中國科學(xué)院 西北高原生物研究所, 青海 西寧810008)

    石油烴污染對隴東黃土高原土壤生物學(xué)及非生物學(xué)特性的影響

    王金成1,2, 井明博1,2, 段春燕2, 王國鋒3, 石國璽1,2,4, 周天林1,2

    (1.甘肅省高校隴東生物資源保護與利用省級重點實驗室, 甘肅 慶陽 745000; 2.隴東學(xué)院 生命科學(xué)與技術(shù)學(xué)院,甘肅 慶陽 745000; 3.陜西師范大學(xué) 生命科學(xué)學(xué)院, 陜西 西安 710062; 4.中國科學(xué)院 西北高原生物研究所, 青海 西寧810008)

    [目的] 探索石油污染與土壤生物學(xué)及非生物學(xué)特性間相互關(guān)系,為隴東黃土高原地區(qū)石油污染土壤的綜合治理提供基礎(chǔ)數(shù)據(jù)和背景資料。[方法] 采用常規(guī)方法測定甘肅省慶陽市5縣,1區(qū),8個采油區(qū)油污土壤化性質(zhì)、酶活性和土壤微生物遺傳多樣性。[結(jié)果] 土壤總石油烴含量(TPHs)與堿解氮、速效磷、pH值和脲酶呈負(fù)相關(guān)(p<0.01),與土壤多酚氧化酶、脫氫酶和過氧化氫酶呈正相關(guān)(p<0.01);樣地土壤TPHs含量低的土壤微生物Shannon-Wiener指數(shù)明顯高于TPHs含量高的樣地;NMDS排序結(jié)果顯示,石油污染對該地區(qū)土壤影響最大的因子為土壤堿解氮和速效磷,其次為TPHs,pH值、脲酶、Shannon-Wiener指數(shù)及過氧化氫酶。[結(jié)論] 隴東地區(qū)石油烴污染對土壤環(huán)境因子的影響主要表現(xiàn)在包括氮、磷等營養(yǎng)元素的有效循環(huán)、土壤脲酶活性和微生物群落多樣性,而對土壤微生物群落豐富度、均勻度和速效鉀的影響未達到顯著水平。

    隴東黃土高原; 石油烴污染; 土壤理化因子; 微生物群落

    文獻參數(shù): 王金成, 井明博, 段春燕, 等.石油烴污染對隴東黃土高原土壤生物學(xué)及非生物學(xué)特性的影響[J].水土保持通報,2017,37(1):009-016.DOI:10.13961/j.cnki.stbctb.2017.01.002; Wang Jincheng, Jing Mingbo, Duan Chunyan, et al. Impact of crude oil pollution on soil biological and abiological properties in Eastern Gansu Province[J]. Bulletin of Soil and Water Conservation, 2017,37(1):009-016.DOI:10.13961/j.cnki.stbctb.2017.01.002

    自20世紀(jì)70年代長慶油田公司從隴東地區(qū)誕生并逐漸發(fā)展壯大,目前成為全國第一大油田,有“西部大慶”之稱。隴東油區(qū)是長慶油田中心產(chǎn)油區(qū),境內(nèi)分布14個油田,63個開發(fā)區(qū),2 000多口油井,正在建設(shè)成為國家級能源化工基地。然而在原油開采過程中已造成隴東地區(qū)石油污染土壤面積達1.02×104hm2,其北部地區(qū)土壤總石油烴(TPHs)含量高達55 137±3 076.84 mg/kg[1]。

    原油污染土壤生態(tài)修復(fù)已成為環(huán)保領(lǐng)域研究熱點,而土壤理化環(huán)境因子所組成的非生物學(xué)特性和土壤微生物群落構(gòu)成的生物學(xué)特性是制約石油生物降解速率和修復(fù)效果的關(guān)鍵因素,因而石油污染土壤環(huán)境因子會直接或間接影響生物修復(fù)技術(shù)的選擇與修復(fù)效應(yīng)的準(zhǔn)確定位[2]。但目前隴東地區(qū)油污土壤相關(guān)研究多集中在修復(fù)效率、修復(fù)方式、降解效果和影響機制以及修復(fù)過程中的微生物群落變化等方面[3],而在干旱少雨、晝夜及四季溫差極大的隴東黃土高原地區(qū)土壤環(huán)境因子對石油污染如何響應(yīng)則缺乏系統(tǒng)分析和研究,因而在該地區(qū)開展生物修復(fù)前對其土壤理化性質(zhì)及微生物生態(tài)進行前期分析顯得尤為重要。

    為此,本研究擬以隴東黃土高原地區(qū)甘肅省慶陽市5縣1區(qū)8個采油區(qū)油污土壤為研究對象,分析上述區(qū)域土壤TPHs,pH值、堿解氮、速效鉀、速效磷、土壤脲酶、脫氫酶、多酚氧化酶、過氧化氫酶活性以及土壤微生物群落遺傳多樣性等13項環(huán)境因子,基于Bray-Curtis相異指數(shù),將所測定的13項土壤環(huán)境因子擬合到8個石油污染樣地的非度量多維尺度(NMDS)排序圖上,旨在揭示石油污染與土壤生物學(xué)及非生物學(xué)特性間相互關(guān)系,對石油污染土壤的綜合治理提供基礎(chǔ)數(shù)據(jù)和背景資料,亦為隴東黃土高原特定區(qū)域下開展石油污染土壤生物修復(fù)打下基礎(chǔ)。

    1 研究區(qū)概況

    于2012年7—8月采樣,采樣點分布在隴東黃土高原慶陽市5縣1區(qū),由南至北,從西向東共計8個鄉(xiāng)鎮(zhèn)的采油區(qū)。采樣點土壤類型為黃綿土、黑爐土和沙土,地理坐標(biāo)地處東經(jīng)107°47′—108°54′,北緯35°65′—36°63′,海拔1 098~1 581 m(表1)。

    表1 研究區(qū)域土樣采集地點概況

    上述區(qū)域?qū)冱S土高原溝壑地貌,是典型水土流失區(qū)之一,侵蝕面積25 050.5 km2,占全區(qū)總面積的92.4%,侵蝕模數(shù)6 383~9 000 t/hm2。屬暖溫帶半干旱季風(fēng)氣候區(qū),光、熱等氣候資源豐富,作物種植制度為旱地一年一熟。由于季風(fēng)的強烈影響,冬季寒冷干燥,春季多風(fēng)且干旱少雨,夏季溫暖多雨,水熱同期,年降水量400~600 mm,雨量分布不均,且主要集中在7—9月,多年平均蒸發(fā)量1 000 mm以上,地面蒸發(fā)量為350~560 mm,地下水位深,土壤瘠薄,水土流失嚴(yán)重[3]。

    2 研究方法

    2.1 土壤樣品采集

    為保證研究區(qū)域內(nèi)每個鄉(xiāng)鎮(zhèn)采樣點代表性,每個樣地設(shè)3個采樣點,采樣點均位于各縣境內(nèi)隨機分布的采油機附近。利用取土鉆,按對角線五點取樣法,在每一樣地中采集5—25 cm深度的土壤,將5個點采集的土壤混勻成1個土樣,所采土樣放在冷藏箱中,及時帶回實驗室根土分離,除利用鮮土提取土壤微生物細(xì)菌基因組總DNA外,其余土樣在室溫背陰處風(fēng)干,過2 mm篩備用[3]。

    2.2 土壤酶活性、理化性質(zhì)及土壤微生物群落遺傳多樣性測定

    過氧化氫酶采用高錳酸鉀滴定法,脲酶和多酚氧化酶采用7230G分光光度計比色法[4];脫氫酶采用唐景春等的分析方法[5]。土壤堿解氮采用堿解擴散法;土壤速效磷、速效鉀分別采用0.5 mol/LNaHCO3浸提—鉬銻抗比色法和火焰光度計法[6];電位法(pHS-3酸度計)測定土樣pH值;采用超聲—索氏萃取—重量法測定土壤中總石油烴含量[7];采用變性凝膠梯度電泳(DGGE)技術(shù)分析微生物群落遺傳多樣性[3]。

    2.3 數(shù)據(jù)處理與計算

    利用SPSS 16.0和R2.15.2(http:∥www.r-project.org/)對數(shù)據(jù)進行處理分析,用Origin 8.0和R 2.15.2作圖。

    為了闡述石油污染樣地與土壤環(huán)境因子間關(guān)系,基于Bray-Curtis相異指數(shù),使用R語言“vegan”程序包中“envfit”程序,將所測定的13項土壤環(huán)境因子擬合到8個石油污染樣地的非度量多維尺度(NMDS)排序圖上,為防止樣點在各象限的分布發(fā)生移位,在擬合時對NMDS的排序軸進行了旋轉(zhuǎn),以確保NMDS第一軸能夠最大程度地代表群落相異性的變異。8個石油污染樣地矩陣及環(huán)境因子數(shù)據(jù)集分別用Bray-Curtis與Euclidean距離表示[8]。

    3 結(jié)果與分析

    3.1 不同污染樣地土壤理化性質(zhì)變化情況

    表2為不同污染樣地土壤理化性質(zhì)、土壤TPHs含量、速效鉀、堿解氮、速效磷、pH值和海拔在各樣地間差異均達極顯著水平(FTPHs=112.398>F0.01,F(xiàn)速效鉀=24.956>F0.01,F(xiàn)堿解氮=50.159>F0.01,F(xiàn)速效磷=49.493>F0.01,F(xiàn)pH=38.132>F0.01,F(xiàn)海拔=54.02>F0.01,p<0.01)。

    不同樣地間土壤TPHs含量的高低所引起的理化環(huán)境因子的變化不盡相同(p<0.05),進一步對不同污染土樣的土壤TPHs含量與其他環(huán)境因子進行相關(guān)性分析顯示,土壤TPHs含量與土壤pH、堿解氮和速效磷含量均呈極顯著負(fù)相關(guān)(表2),而與速效鉀和樣地海拔間無顯著相關(guān),說明土壤TPHs的增加可能是制約隴東地區(qū)油污土壤氮、磷元素有效化的主要因素。

    表2 研究區(qū)不同污染樣地土壤化學(xué)性質(zhì)的變化

    注:TPHs為土壤總石油烴含量; 數(shù)字為平均值±標(biāo)準(zhǔn)誤; 每組同列的不同字母表示差異顯著性;小寫字母為5%顯著水平。下同。

    3.2 不同污染樣地土壤酶活性變化情況

    表3為不同污染樣地土壤酶活性變化情況。不同樣地間土壤脫氫酶、多酚氧化酶、過氧化氫酶和脲酶差異極顯著(F脫氫酶=173.512>F0.01,F(xiàn)多酚氧化酶=89.104>F0.01,F(xiàn)過氧化氫酶=111.041>F0.01,F(xiàn)脲酶=50.023>F0.01,p<0.01)。與土壤TPHs含量相關(guān)性分析后發(fā)現(xiàn),除與土壤脲酶呈負(fù)相關(guān)外,土壤TPHs與土壤多酚氧化酶、脫氫酶和過氧化氫酶均呈正相關(guān)關(guān)系(表3),說明土壤TPHs可增加土壤中脫氫酶、過氧化氫酶、多酚氧化酶的活性,參與TPHs降解。而參與土壤養(yǎng)分有效化的脲酶活性與土壤TPHs呈負(fù)相關(guān)說明隴東地區(qū)石油污染土壤養(yǎng)分有效化水平低下,這一結(jié)果亦可解釋土壤理化性質(zhì)中堿解氮和速效磷含量低下的原因。

    表3 研究區(qū)不同污染樣地土壤酶活性的變化

    3.3 不同污染樣地土壤微生物群落遺傳多樣性變化情況

    本文采用變性凝膠梯度電泳(DGGE)技術(shù)分析微生物群落遺傳多樣性,利用3種多樣性指標(biāo)進行土壤微生物群落遺傳多樣性分析,即Patrick 豐富度指數(shù)(S):S=物種數(shù),即條帶數(shù);Shannon-Wiener多樣性指數(shù)(H):H=-ΣPilnPi,Pi代表土壤樣品中第i個條帶的吸光度占所有條帶吸光度總和的比例;Pielou 均勻度指數(shù)(J):J=H/Hmax,Hmax代表最大多樣性值(Hmax=lnS)[3]。圖1為不同樣地油污土壤PCR-DGGE指紋圖譜及其圖譜聚類結(jié)果。圖譜上條帶信息反映了土壤微生物群落多樣性變化情況,其中電泳條帶的多少反映出土壤中細(xì)菌群落多樣性;條帶的粗細(xì)則反映了種群密度的差異。

    圖2反映了不同樣地油污土壤樣品微生物群落遺傳多樣性指數(shù)變化情況。F檢驗結(jié)果顯示,3種多樣性指數(shù)在各樣地間差異均達極顯著水平(FH=25.018>F0.01,F(xiàn)J=16.121>F0.01,F(xiàn)S=16.121>F0.01,p<0.01),上述結(jié)果體現(xiàn)了TPHs對土壤微生物群落結(jié)構(gòu)的影響。多重比較(Duncan,α=0.05)結(jié)果顯示,不同樣地間土壤TPHs含量的高低所引起的土壤微生物群落遺傳多樣性的變化不盡相同(p<0.05)。

    由圖2可見,土壤TPHs含量介于7 995.61~8 992.74 mg/kg樣地Shannon-Wiener指數(shù)(HYCZ=2.79±0.01,HCHX=2.77±0.02,HGCX=2.88±0.01,HTZZ=2.75±0.03)高于土壤TPHs含量介于17 139.95至22 348.03 mg/kg的樣地(HFJCX=2.62±0.01,HDZZ=2.59±0.02,HTCX=2.52±0.02,HQZZ=2.55±0.01)(p<0.05)。Patrick指數(shù)方面,除樣地DZZ以外,其余樣地均表現(xiàn)出與Shannon-Wiener指數(shù)相似的變化趨勢,即Shannon-Wiener指數(shù)越高,物種條帶數(shù)越多。此外,各樣地間Pielou 均勻度指數(shù)偏低,說明TPHs污染造成土壤微生物群落分布不均。

    圖1 研究區(qū)不同樣地油污土壤PCR-DGGE 指紋圖譜及其圖譜聚類

    3.4 不同污染樣地與其環(huán)境因子的非度量多維尺度(NMDS)排序

    為了闡述石油污染樣地與土壤環(huán)境因子間關(guān)系,以及環(huán)境因子在驅(qū)動不同樣地分異過程中的相對重要性,基于Bray-Curtis相異指數(shù),將所測定的13項土壤環(huán)境因子擬合到8個石油污染樣地的非度量多維尺度(NMDS)排序圖上,脅強系數(shù)Stress=0.134 3,說明該二維NMDS分析具有一定的擬合效果,可解釋樣地環(huán)境因子信息量的86.57%。

    圖2 研究區(qū)不同樣地油污土壤樣品微生物群落遺傳多樣性指數(shù)變化

    由圖3及表4可見,基于Bray-Curtis相異性指數(shù),本研究所調(diào)查的8個樣地在NMDS排序圖上較為分散,說明8個樣地相似性較低,即樣地間存在顯著差異,表現(xiàn)為8個樣地在NMDS排序圖分異為8個點集,相對而言樣地YCZ,CHX,GCX和TZZ之間相似性較高,位于慶陽市中部的樣地DZZ和TCX與其他樣地差異較大,而位于北部環(huán)縣境內(nèi)的樣地FJCX和QZZ呈現(xiàn)較高的相似性(圖3,表4)。13個環(huán)境因子變量擬合至NMDS排序圖的結(jié)果顯示,其中有9個環(huán)境變量與樣地NMDS排序結(jié)果顯著相關(guān)。

    注:排序圖上顯示了與各樣地分異有關(guān)的土壤環(huán)境因子,橢圓代表平行樣品間的標(biāo)準(zhǔn)差。圖3 研究區(qū)不同污染樣地的非度量多維尺度(NMDS)排序

    環(huán)境因子r2p值TPHs0.54480.002**速效鉀0.24740.058堿解氮0.65150.001***速效磷0.60230.001***pH值0.52210.002**脫氫酶0.31480.029*多酚氧化酶0.34070.026*過氧化氫酶0.40250.003**脲酶0.51130.005**海拔高度0.20950.085Shannon-Wiener多樣性指數(shù)0.41540.002**物種數(shù)0.20830.064Pielou均勻度指數(shù)0.19690.091

    注:*代表與各樣地分異相關(guān)的土壤環(huán)境因子顯著水平:***表示p<0.001,**表示p<0.01,*表示p<0.05時相關(guān)性顯著。

    由圖3可見樣地在NMDS排序圖上的分異主要與土壤理化因子中的堿解氮含量(r2=0.651 5,p=0.001),速效磷(r2=0.602 3,p=0.001),土壤TPHs含量(r2=0.544 8,p=0.002)及土壤pH(r2=0.522 1,p=0.002)顯著相關(guān),其中堿解氮和速效磷含量是導(dǎo)致樣地分異的主要環(huán)境因子;土壤酶活性方面,與過氧化氫酶(r2=0.402 5,p=0.003),脲酶(r2=5 113,p=0.005),多酚氧化酶(r2=0.340 7,p=0.003)以及脫氫酶(r2=0.314 8,p=0.003)顯著相關(guān),其中過氧化氫酶和脲酶在樣地分異過程中起主要作用;就土壤生物學(xué)特性而言,3個土壤微生物群落遺傳多樣性指數(shù)中只有Shannon-Wiener多樣性指數(shù)(H)與樣地分異有關(guān)(r2=0.415 4,p=0.002)。上述結(jié)果說明,隴東黃土高原地區(qū)受石油污染影響最大的環(huán)境因子為土壤堿解氮和速效磷,其次為土壤TPHs,pH值、脲酶、微生物群落多樣性和過氧化氫酶,而速效鉀、海拔、微生物群落豐富度(S)與均勻度(J)所受影響較小。

    4 討論與結(jié)論

    4.1 石油烴污染對隴東黃土高原土壤非生物學(xué)特性的影響

    由土壤理化性質(zhì)及酶活性所構(gòu)成的非生物學(xué)環(huán)境因子是生物地球化學(xué)循環(huán)的主要參與者,可反映土壤環(huán)境狀況[2],是評價土壤質(zhì)量及其污染程度的重要指示。本文對隴東黃土高原地區(qū)慶陽市境內(nèi)5縣1區(qū)石油污染土壤非生物學(xué)環(huán)境因子的分析結(jié)果顯示不同濃度TPHs所引起的土壤理化及酶活性變化不盡相同。相關(guān)性分析結(jié)果顯示土壤TPHs含量與土壤pH值、堿解氮和速效磷含量均呈極顯著負(fù)相關(guān)(表3),與土壤速效鉀及樣地海拔間無顯著性相關(guān)。前人研究指出石油污染土壤的氮素匱乏為微生物提供了選擇優(yōu)勢,進而引起石油污染土壤中堿解氮含量升高,而王傳遠(yuǎn)等[2]研究認(rèn)為石油污染可降低土壤pH值,但與堿解氮含量無顯著相關(guān)。劉五星等[9]對南京某煉油廠的油污土壤進行理化性質(zhì)分析后指出石油污染對土壤的速效氮、磷、鉀含量無顯著影響。上述結(jié)果與本研究結(jié)果不盡相同,說明土壤TPHs對隴東黃土高原地區(qū)土壤理化性質(zhì)的影響具有明顯的獨特性和地域性。究其原因可能是由于長慶油田隴東產(chǎn)油區(qū)原油微量元素相對較低,而該地區(qū)土壤本身缺磷、少氮、多鉀,因而石油烴污染土壤后導(dǎo)致外源碳源增加[10],打破了隴東黃土高原土壤原有營養(yǎng)水平比例,進而改變土壤微生物群落結(jié)構(gòu),引起土著石油降解菌的大量生長,而這些土著石油降解菌在生長的同時大量消耗了土壤中原有的氮、磷等營養(yǎng)物質(zhì)進而造成土壤氮、磷含量嚴(yán)重不足;其次,本研究結(jié)果顯示土壤TPHs與土壤脲酶活性與呈負(fù)相關(guān)(表3),而脲酶可加速土壤養(yǎng)分有效化[11],因此土壤脲酶活性低下可能是隴東黃土高原地區(qū)堿解氮和速效磷含量低下的又一原因,而該地區(qū)油污土壤速效鉀含量在不同樣地間高低不一的原因有待進一步分析研究;第三,前人研究指出石油烴污染濃度也是影響土壤營養(yǎng)水平有效化的又一因素,張曉陽等[12]對陜北石油污染對土壤理化性質(zhì)影響時發(fā)現(xiàn),當(dāng)土壤TPHs超過5 000 mg/kg時,土壤營養(yǎng)水平有效化顯著降低,而本研究8各樣地土壤TPHs含量介于7 795.61~22 348.03 mg/kg(表1),加之該地區(qū)土壤本身營養(yǎng)元素水平特征,說明重度污染時土壤速效養(yǎng)分水平低下。

    在土壤TPHs降解過程中氧化還原酶系占有重要地位,對土壤氧化還原酶的研究多集中在土壤脫氫酶、過氧化氫酶和多酚氧化酶。李慧等[13]通過對沈撫灌區(qū)含油污水灌溉對土壤酶活性的影響中發(fā)現(xiàn),灌溉含石油烴的污水能增加上述3種酶活性,且隨著污染程度的增加而增加。朱凡等[14]在利用植物修復(fù)PAHs污染土壤過程發(fā)現(xiàn),土壤過氧化氫酶、多酚氧化酶活性隨著修復(fù)時間增加呈先升后降的趨勢。本研究結(jié)果發(fā)現(xiàn),除與土壤脲酶呈負(fù)相關(guān)外,土壤TPHs與土壤多酚氧化酶、脫氫酶和過氧化氫酶均呈正相關(guān)關(guān)系(表4),說明土壤TPHs可提高上述3種酶活性進而加速TPHs降解,這與眾多學(xué)者研究結(jié)果較為一致,土壤TPHs對隴東黃土高原地區(qū)土壤酶活性的影響沒有明顯的地域性。本研究中,相對于正在采油作業(yè)的樣地FJCX,DZZ,TCX和QZZ,樣地YCZ,CHX,GCX和TZZ均已停止采油廢棄2~4 a,且樣地YCZ,CHX,GCX和TZZ的TPHs含量明顯低于在采油作業(yè)的樣地,生態(tài)毒性相對較弱,使得TPHs的自然降解使得土壤微生物群落多樣性和數(shù)量趨于穩(wěn)定,而隴東黃土高原地區(qū)降雨量極低,加之TPHs在土壤遷移和擴散程度相對緩慢[15],造成石油降解菌富集的局域微環(huán)境中土壤TPHs含量明顯降低,因而樣地YCZ,CHX,GCX和TZZ土壤脫氫酶活性相對偏低。過氧化氫酶活性可以反映土壤除去呼吸過程中產(chǎn)生的過氧化氫的能力,若過氧化氫持續(xù)積累且無法有效分解將會對土壤微生物產(chǎn)生不利影響。由表4可見,土壤TPHs含量的增加有助于當(dāng)?shù)剡m生土著石油降解菌生長繁殖,進而引起其群落結(jié)構(gòu)和組成發(fā)生改變,加速了土壤TPHs降解,而在生物降解過程中過氧化氫的積累使得微生物及土壤生境擴大,解除了由于呼吸作用而產(chǎn)生的過氧化氫,進而使得過氧化氫酶活性得到提升[16]。因此,隴東黃土高原地區(qū)開展石油烴污染生態(tài)修復(fù)首先應(yīng)外援添加相應(yīng)的N,P營養(yǎng)元素加快對石油烴的分解,其次如何提高土壤脲酶活性,提升營養(yǎng)元素有效化程度是決定該地區(qū)修復(fù)效果的關(guān)鍵因素。

    4.2 石油烴污染對隴東黃土高原土壤生物學(xué)特性的影響

    微生物種群結(jié)構(gòu)與多樣性都是表征生態(tài)系統(tǒng)群落結(jié)構(gòu)的重要參數(shù),其對環(huán)境污染物的反應(yīng)表現(xiàn)為多種形式。TPHs進入土壤后對生態(tài)環(huán)境的影響首先表現(xiàn)為土壤微生物群落結(jié)構(gòu)組成、數(shù)量及多樣性的改變。通常認(rèn)為石油污染導(dǎo)致土壤中微生物多樣性降低,而劉五星等[9]對南京某煉油廠石油污染土壤微生物群落遺傳多樣性分析后指出,土壤TPHs與微生物群落多樣性及均勻度指數(shù)多樣性呈正相關(guān)關(guān)系;姜睿玲等[15]研究了土壤PAHs對桑園土壤微生物群落遺傳多樣性的影響,結(jié)果顯示低濃度PAHs對土壤微生物多樣性及豐富度指數(shù)具有刺激作用,而高濃度污染時可明顯抑制土壤微生物多樣性及豐富度指數(shù)。Li等[17]研究指出水稻土中石油污染濃度與細(xì)菌多樣性正相關(guān)。本研究中土壤TPHs含量介于7 995.61~8 992.74 mg/kg樣地(YCZ,CHX,GCX和TZZ)土壤微生物群落Shannon-Wiener指數(shù)明顯高于土壤TPHs含量介于17 139.95~22 348.03 mg/kg的樣地(FJCX,DZZ,TCX和QZZ)。Patrick指數(shù)方面,除樣地DZZ以外,其余樣地均表現(xiàn)出與Shannon-Wiener指數(shù)相似的變化趨勢,即Shannon-Wiener指數(shù)越高,物種條帶數(shù)越多(圖1)。此外,各樣地間Pielou 均勻度指數(shù)偏低,說明TPHs污染造成了土壤微生物群落分布不均。上述試驗結(jié)果與本研究結(jié)果不盡相同,究其原因除土壤質(zhì)地和污染程度不同外,亦可能是土壤TPHs污染對土壤微生物群落結(jié)構(gòu)產(chǎn)生了影響,且樣地YCZ,CHX,GCX和TZZ均已廢棄,因而TPHs在長期自然降解過程中生態(tài)毒性相對減弱,土壤微生物逐漸適應(yīng)該種生態(tài)條件,驅(qū)使其群落多樣性和數(shù)量趨于穩(wěn)定,因而微生物群落種類和數(shù)量相對較高,表現(xiàn)為Patrick豐富度指數(shù)高于樣地FJCX,DZZ,TCX和QZZ。其次,新開采油井中,土壤TPHs污染對土壤環(huán)境因子的改變強化了生態(tài)選擇,進而改變了土壤微生物群落結(jié)構(gòu),引起多數(shù)微生物種類及數(shù)量上的下降和分布不均,但這種變化可能會隨著石油烴自然降解過程的進行而逐漸緩解[18];第三,隴東黃土高原地區(qū)生態(tài)脆弱,南低北高,氣候干旱少雨,寒暑溫差顯著,土壤瘠薄造就了該地區(qū)土壤環(huán)境因子地域性變化,進而也造成了TPHs在土壤遷移和擴散程度不同,各樣地間Pielou 均勻度指數(shù)總體偏低。因此,隴東地區(qū)石油烴污染生物修復(fù)過程中應(yīng)外援投加“適生”的土著石油降解菌群,強化土著降解菌在該地區(qū)油污土壤生態(tài)位亦為影響修復(fù)效果的關(guān)鍵因素。

    4.3 石油烴污染對隴東黃土高原土壤因子的影響

    基于上述土壤環(huán)境因子測定結(jié)果,以Bray-Curtis相異性指數(shù)為基礎(chǔ)構(gòu)建隴東黃土高原地區(qū)石油污染樣地間相似性矩陣,即土壤環(huán)境因子相似程度高,其樣地排序距離就小,采用NMDS方法分析該地區(qū)石油污染后土壤環(huán)境因子驅(qū)動樣地分異的影響機制。結(jié)果表明(圖3),樣地YCZ,CHX,GCX和TZZ之間相似性較高,位于慶陽市中部的樣地DZZ和TCX與其他樣地差異較大,而位于北部環(huán)縣境內(nèi)的樣地FJCX和QZZ呈現(xiàn)較高的相似性。不難看出樣地YCZ,CHX,GCX和TZZ均為易廢棄2~4 a油井,而這4各樣地相似性高可能是由于廢棄后土著降解菌自然降解后土壤環(huán)境因子趨同所致;而樣地DZZ和TCX均位于慶陽市中部,且均為正在采油作業(yè)油井;環(huán)縣境內(nèi)的樣地FJCX和QZZ亦為采油作業(yè)1 a的新油井。根據(jù)NMDS排序結(jié)果看,雖然這4個樣地均為正在采油作業(yè)的新油井但卻劃分為2個板塊,說明造成這一結(jié)果的原因必然與土壤環(huán)境因子有關(guān)。為了闡述石油污染樣地與土壤環(huán)境因子間關(guān)系,使用R語言“vegan”程序包中“envfit”程序,將所測定的13項土壤環(huán)境因子擬合到8個石油污染樣地的非度量多維尺度(NMDS)排序圖上(圖3),環(huán)境因子數(shù)據(jù)用Euclidean距離表示[8]。首先樣地YCZ,CHX,GCX和TZZ之間相似性較高主要與土壤環(huán)境因子中的土壤堿解氮含量(r2=0.651 5,p=0.001),速效磷(r2=0.602 3,p=0.001),土壤脲酶活性(r2=0.511 3,p=0.005),脫氫酶活性(r2=0.314 8,p=0.029),土壤pH值(r2=0.522 1,p=0.002)及土壤微生物群落Shannon-Wiener指數(shù)(H)(r2=0.415 4,p=0.002)有關(guān),其中堿解氮和速效磷含量是造成上述樣地相似性高的最主要環(huán)境因子。然而樣地YCZ,CHX,GCX和TZZ之間相較之下,樣地CHX和GCX間與樣地YCZ和TZZ之間排序距離更小,究其原因,由圖3可見,決定樣地CHX和GCX相似性更高的原因與土壤堿解氮,速效磷和土壤脲酶活性有關(guān),說明樣地CHX和GCX的土壤脲酶提高了土壤營養(yǎng)元素的有效化水平,進而以較高的土壤堿解氮和速效磷含量使得樣地CHX和GCX不同于其他樣地。而樣地YCZ和TZZ方面,決定這2個相似性更高的原因與土壤微生物群落Shannon-Wiener指數(shù)(H),pH值和土壤脫氫酶活性有關(guān),其中土壤pH值最為相關(guān)(r2=0.522 1,p=0.002)。微生物總數(shù)與土壤脫氫酶活性成正相關(guān)關(guān)系,其活性高低可以反映微生物量對有機物的降解效果[16]。本研究結(jié)果與此一致,決定樣地YCZ和TZZ不同于其他樣地的環(huán)境因子主要包括土壤脫氫酶活性和微生物群落Shannon-Wiener指數(shù)(H)以及pH值,說明在自然降解過程中這2個樣地的土壤微生物群落遺傳多樣性較為相似,且與土壤脫氫酶活性呈正相關(guān)關(guān)系;而影響土壤微生物群落遺傳多樣性的環(huán)境因子中土壤pH值與之相關(guān),前人研究指出認(rèn)為石油污泥的pH值呈堿性,超過了微生物生長的最適范圍,進而對土壤微生物的生長有很強的抑制作用。土壤pH值可能會通過影響降解過程中關(guān)鍵酶的酶活性來影響降解率,說明土壤pH值的改變有助于自然降解過程中土著石油降解菌的生長繁殖,樣地YCZ和TZZ土壤TPHs自然降解進程優(yōu)于其他樣地。

    就正在采油作業(yè)油井樣地FJCX和QZZ與DZZ和TCX而言,決定慶陽市中部的樣地DZZ和TCX相似性更高的原因與土壤過氧化氫酶(r2=0.402 5,p=0.003)與多酚氧化酶活性(r2=0.340 7,p=0.026)有關(guān),相較之下過氧化氫酶活性最為相關(guān),而這兩種酶均參與石油烴及其代謝中間產(chǎn)物在土壤中的降解過程,特別是多酚氧化酶能把土壤中芳香族化合物和酚類氧化成分子量不等的有機質(zhì)和色素,完成土壤芳香族化合物的循環(huán)[19]。而該樣地開采期均為2 a,石油烴污染后強化了對該樣地土壤微生物種群的選擇,有助于石油降解菌群的生長繁殖,刺激石油降解菌產(chǎn)生多酚氧化酶。石油降解菌群的大量繁殖加強對該樣地石油烴TPHs的生物降解,而降解過程中產(chǎn)生的過氧化氫產(chǎn)生使得土壤石油降解菌生境逐漸擴大,因此慶陽市中部的樣地DZZ和TCX不同于其北部環(huán)縣境內(nèi)的樣地FJCX和QZZ。圖3中與環(huán)縣境內(nèi)的樣地FJCX和QZZ相似性較高有關(guān)的土壤環(huán)境因子只有土壤TPHs(r2=0.544 8,p=0.002),說明土壤石油烴含量是決定該樣地不同于其他樣地的環(huán)境因子??v觀上述結(jié)果,影響隴東黃土高原地區(qū)石油污染土壤的環(huán)境因子中,雖然土壤微生物群落豐富度(S)(r2=0.208 3,p=0.064)、均勻度(J)(r2=0.196 9,p=0.091)、速效鉀含量(r2=0.247 4,p=0.058)和海拔(r2=0.209 5,p=0.085)并未達到顯著水平,但慶陽市南低北高的地理特點和降雨量不均及晝夜溫差大等因素在本次研究中并未涉及,而上述環(huán)境因子可能是解釋環(huán)縣境內(nèi)樣地FJCX和QZZ不同于慶陽市中部的樣地DZZ和TCX的原因,而這些環(huán)境因素對隴東地區(qū)石油污染土壤的影響有待后續(xù)分析研究。

    [1] 王金成,周天林,井明博,等.隴東黃土高原地區(qū)石油污泥原位修復(fù)過程中土壤主要肥力指標(biāo)動態(tài)變化分析[J].環(huán)境科學(xué)學(xué)報,2015,35(1):280-287.

    [2] 王傳遠(yuǎn),楊翠云,孫志高,等.黃河三角洲生態(tài)區(qū)土壤石油污染及其與理化性質(zhì)的關(guān)系[J].水土保持學(xué)報,2010,24(2):214-217.

    [3] 王金成,井明博,段春艷,等.隴東黃土高原石油污染土壤環(huán)境因子對金盞菊(Calendulaofficinalis)—微生物聯(lián)合修復(fù)的響應(yīng)[J].環(huán)境科學(xué)學(xué)報,2015,35(9):2971-2981.

    [4] 嚴(yán)昶升.土壤肥力研究方法[M].北京:農(nóng)業(yè)出版社,1988:263-280.

    [5] 唐景春,王斐,褚洪蕊,等.玉米草(ZeaMexicana)與海藻寡糖聯(lián)合修復(fù)石油烴污染土壤的研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2010,29(11):2107-2113.

    [6] 鮑士旦.土壤農(nóng)化分析[M].3版.北京:中國農(nóng)業(yè)出版社,1999:245-353.

    [7] 王如剛,王敏,牛曉偉,等.超聲—索氏萃取—重量法測定土壤中總石油烴含量[J].分析化學(xué),2010,38(3):417-420.

    [8] Shi Guoxi, Liu Yongjun, Johnson N C, et al., Interactive influence of light intensity and soil fertility on root-associated arbuscular mycorrhizal fungi. [J]. Plant and Soil, 2014,378(1/2):173-188.

    [9] 劉五星,駱永明,滕應(yīng),等.石油污染土壤的生態(tài)風(fēng)險評價和生物修復(fù)(Ⅱ):石油污染土壤的理化性質(zhì)和微生物生態(tài)變化研究[J].土壤學(xué)報,2007,44(5):848-853.

    [10] 王金成,井明博,肖朝霞,等.隴東黃土高原地區(qū)石油污染土壤微生物群落及其與環(huán)境因子的關(guān)系[J].水土保持通報,2012,32(5):145-151.

    [11] Antonio D A, Francesca B, Massimo G,etal., Bioremediation of petroleum hydrocarbons in anoxicmarine sediments:Consequences on the speciation of heavy metals [J]. Marine Pollution Bulletin, 2009,58(12):1808-1814.

    [12] 張曉陽,李凱榮,張麟君.陜北石油污染對土壤理化性質(zhì)的影響[J].水土保持研究,2013,20(3):32-38.

    [13] 李慧,陳冠雄,楊濤,等.沈撫灌區(qū)含油污水灌溉對稻田土壤微生物種群及土壤酶活性的影響[J].應(yīng)用生態(tài)學(xué)報, 2005,16(7):1355-1359.

    [14] 朱凡,洪湘琦,閆文德,等. PAHs污染土壤植物修復(fù)對酶活性的影響[J].生態(tài)學(xué)報,2014,34(3):581-588.

    [15] 姜睿玲,楊統(tǒng)一,唐玉斌,等.多環(huán)芳烴污染對桑園土壤微生物結(jié)構(gòu)及種群多樣性的影響[J].中國環(huán)境科學(xué),2012,32(9):1655-1661.

    [16] Gnseppe T, Maria L C. Guido G. Oxidative polymerisation of phenols by a phenol oxidase from green olives [J]. Enzyme and Microbial technology, 2003,33(1):47-54.

    [17] Li Hui, Zhang Chenggang, Cheng Guanxiong, et al., Effect of petroleum-containing wastewater irrigation on bacterial diversities and enzymatic activities in a paddy soil irrigation area [J]. Journal of Environmental Quality, 2005,34(3):1073-1080.

    [18] 楊茜,吳蔓莉,曹碧霄,等.石油降解菌的篩選、降解特性及其與基因的相關(guān)性研究[J].安全與環(huán)境學(xué)報,2014,14(1):187-192.

    [19] Gnseppe T, Maria L C. Guido G. Oxidative polymerisation of phenols by a phenol oxidase from green olives[J]. Enzyme and Microbial technology, 2003,33(1):47-54.

    Impact of Crude Oil Pollution on Soil Biological and Abiological Properties in Eastern Gansu Province

    WANG Jincheng1,2, JING Mingbo1,2, DUAN Chunyan2, WANG Guofeng3, SHI Guoxi1,2,4, ZHOU Tianlin1,2

    (1.ProvincialKeyUniversityLaboratoryforProtectionandUtilizationofLongdongBioresourcesinGansuProvince,Qingyang,Gansu, 745000,China; 2.CollegeofLifeScienceandTechnology,LongdongUniversity,Qingyang,Gansu, 745000,China; 3.CollegeofLifeScience,ShaanxiNormalUniversity,Xi’an,Shaanxi710062,China; 4.ChineseAcademyofSciences,NorthwestInstituteofPlateauBiology,Xining,Qinghai810008,China)

    [Objective] The objective of this study is to explain relationship between soil biology and abiological characteristics on crude-oil contaminated soil and provide basic data and background information for its comprehensive remediation in eastern Gansu Province. [Methods] The changes of soil physical and chemical properties as well as enzyme activities in crude oil-contaminated soil in Qingyang City of Gansu Province were investigated using the conventional methods, and the genetic diversity of microbial community was analyzed based on PCR-DGGE. [Results] The contents of soil available N, available P, pH value and activity of soil urease decreased with the increase of soil TPHs concentration respectively(p<0.01), while the activities of soil polyphenoloxidase, dehydrogenase and catalase had a significant positive correlation with soil TPHs concentration respectively(p<0.01). As for the genetic diversity of soil microbial communities, Shannon-Wiener index in the studied sites with lower TPHs concentration was significantly higher than that of the sites with higher TPHs concentration(p<0.05). Moreover, nonmetric multidimensional scaling(NMDS) showed that petroleum-contaminated soil in Gansu Province were strongly affected by soil available N and available P, and following-affected by soil TPHs concentration, pH value, urease activity, Shannon-Wiener index of soil microbial communities and activity of soil catalase.[Conclusion] Impact of crude oil pollution on soil environmental factors in Eastern Gansu Province of Loess Plateau was mainly reflected by the nutrients effective circulation of nitrogen and phosphorus as well as soil microbial-communities diversity and soil urease activity, but there were no significant impact on available K, richness and evenness of soil microbial communities.

    Eastern Gansu Province of Loess Plateau; crude oil pollution; soil environmental factors; microbial community

    2016-06-21

    2016-09-20

    甘肅省科技支撐計劃項目“隴東黃土高原地區(qū)石油污染土壤原位修復(fù)技術(shù)研究與示范”(1204FKCM173); 甘肅省教育廳科技計劃項目“隴東黃土高原石油污染土壤微生物群落動態(tài)分析及石油降解菌的篩選”(2014B-091); 國家自然科學(xué)資助項目(31500427); 甘肅省慶陽市科技支撐計劃項目(GC2011-16); 隴東學(xué)院青年科技創(chuàng)新項目(XYZK1603)

    王金成(1985—),男(漢族),甘肅省兩當(dāng)縣人,碩士,講師,主要從事環(huán)境微生物學(xué)及分子微生物生態(tài)學(xué)領(lǐng)域研究工作。E-mail:ldxywjc@163.com。

    周天林(1961—),男(漢族),甘肅省寧縣人,大學(xué)本科,教授,主要從事土壤生態(tài)學(xué)研究工作。E-mail:wensent20002 002@163.com。

    A

    1000-288X(2017)01-0009-08

    X53, X172

    猜你喜歡
    隴東土壤環(huán)境樣地
    隴東能源大數(shù)據(jù)中心
    額爾古納市興安落葉松中齡林植被碳儲量研究
    不同種植區(qū)隴東苜蓿營養(yǎng)價值的比較研究
    中國飼料(2021年17期)2021-11-02 08:15:24
    昆明市主要綠化樹種閾值測定與分析
    基于角尺度模型的林業(yè)樣地空間結(jié)構(gòu)分析
    隴東黃土高原軟棗獼猴桃栽培技術(shù)
    有機氯農(nóng)藥對土壤環(huán)境的影響
    土壤環(huán)境安全及其污染防治對策
    土壤環(huán)境質(zhì)量監(jiān)測的現(xiàn)狀及發(fā)展趨勢
    一个人免费在线观看的高清视频| 可以在线观看毛片的网站| 免费观看人在逋| 少妇丰满av| 99久久久亚洲精品蜜臀av| 色5月婷婷丁香| 国产在线精品亚洲第一网站| 亚洲成a人片在线一区二区| 99精品久久久久人妻精品| 最近中文字幕高清免费大全6 | 十八禁国产超污无遮挡网站| av黄色大香蕉| 国产熟女xx| 欧美日韩黄片免| 亚洲国产精品999在线| 亚洲av成人av| 欧美成人性av电影在线观看| 真人一进一出gif抽搐免费| 麻豆av噜噜一区二区三区| 老女人水多毛片| 久久精品国产亚洲av涩爱 | 国产野战对白在线观看| 久久精品国产亚洲av香蕉五月| 99久久精品热视频| 亚州av有码| 国产精品自产拍在线观看55亚洲| 国产精品爽爽va在线观看网站| 一夜夜www| 久久午夜福利片| 国产精品美女特级片免费视频播放器| 精品午夜福利视频在线观看一区| 最近最新中文字幕大全电影3| 天堂av国产一区二区熟女人妻| 国产伦在线观看视频一区| 波野结衣二区三区在线| 我的老师免费观看完整版| 久久精品人妻少妇| eeuss影院久久| 亚洲激情在线av| a级毛片免费高清观看在线播放| 欧美性感艳星| 色播亚洲综合网| 久久午夜亚洲精品久久| 夜夜夜夜夜久久久久| 真实男女啪啪啪动态图| 可以在线观看毛片的网站| 丁香六月欧美| 日本熟妇午夜| 1000部很黄的大片| 久久人人精品亚洲av| 色吧在线观看| 国产精品精品国产色婷婷| 日本 欧美在线| 精品人妻熟女av久视频| 亚洲精品乱码久久久v下载方式| 91在线观看av| 两个人的视频大全免费| 国产男靠女视频免费网站| 国产欧美日韩一区二区三| 中文字幕熟女人妻在线| 夜夜爽天天搞| 欧美xxxx性猛交bbbb| 蜜桃久久精品国产亚洲av| 成人精品一区二区免费| 国产白丝娇喘喷水9色精品| 成人高潮视频无遮挡免费网站| 亚洲电影在线观看av| 午夜影院日韩av| 高清日韩中文字幕在线| 一个人免费在线观看的高清视频| 欧美日韩乱码在线| 一边摸一边抽搐一进一小说| 久9热在线精品视频| 首页视频小说图片口味搜索| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看日本一区| 99国产精品一区二区三区| 国产精品一区二区性色av| 国产三级在线视频| 99热只有精品国产| 51午夜福利影视在线观看| 麻豆成人午夜福利视频| 人妻久久中文字幕网| 九色国产91popny在线| av专区在线播放| 白带黄色成豆腐渣| 身体一侧抽搐| 午夜福利18| eeuss影院久久| 一二三四社区在线视频社区8| 在线观看美女被高潮喷水网站 | 亚洲第一欧美日韩一区二区三区| 少妇人妻精品综合一区二区 | 亚洲五月婷婷丁香| 久久久国产成人免费| 午夜福利成人在线免费观看| 午夜久久久久精精品| 免费av观看视频| 免费观看精品视频网站| 亚洲国产精品合色在线| 亚洲成人久久性| 美女大奶头视频| 国产av不卡久久| 免费搜索国产男女视频| 亚洲欧美日韩高清专用| 男人狂女人下面高潮的视频| 免费人成在线观看视频色| 成人国产综合亚洲| 99久久久亚洲精品蜜臀av| 最新中文字幕久久久久| 舔av片在线| 1024手机看黄色片| 真人做人爱边吃奶动态| 成人三级黄色视频| 日韩高清综合在线| 中文字幕av在线有码专区| 婷婷精品国产亚洲av在线| 听说在线观看完整版免费高清| 最新在线观看一区二区三区| 亚洲av二区三区四区| 高清日韩中文字幕在线| 亚洲欧美激情综合另类| 精品福利观看| 色吧在线观看| 日本撒尿小便嘘嘘汇集6| 欧美日本视频| 黄色女人牲交| 精品欧美国产一区二区三| 可以在线观看毛片的网站| 国产三级中文精品| 久久天躁狠狠躁夜夜2o2o| 久久国产精品人妻蜜桃| 色哟哟哟哟哟哟| 免费电影在线观看免费观看| 18禁黄网站禁片免费观看直播| 在线观看66精品国产| 一区二区三区高清视频在线| 久久精品国产清高在天天线| 舔av片在线| 久久久久九九精品影院| 国产69精品久久久久777片| 久久国产精品人妻蜜桃| 老女人水多毛片| 国产精品一及| 久久精品夜夜夜夜夜久久蜜豆| 久久这里只有精品中国| 在线免费观看的www视频| 淫秽高清视频在线观看| 欧美在线一区亚洲| 非洲黑人性xxxx精品又粗又长| 桃色一区二区三区在线观看| 男女那种视频在线观看| 99久久99久久久精品蜜桃| 国产 一区 欧美 日韩| 91狼人影院| 成人国产一区最新在线观看| АⅤ资源中文在线天堂| 亚洲五月天丁香| 能在线免费观看的黄片| 中文字幕人妻熟人妻熟丝袜美| 赤兔流量卡办理| 国产精品一区二区三区四区免费观看 | 欧美成人一区二区免费高清观看| 成年人黄色毛片网站| 99久久成人亚洲精品观看| 欧美最黄视频在线播放免费| 国产欧美日韩精品一区二区| 又紧又爽又黄一区二区| 日本与韩国留学比较| 18禁在线播放成人免费| 夜夜夜夜夜久久久久| 亚洲欧美日韩高清在线视频| 免费av不卡在线播放| 我的女老师完整版在线观看| 一区二区三区免费毛片| 欧美色欧美亚洲另类二区| 国产又黄又爽又无遮挡在线| 亚洲第一区二区三区不卡| 99在线人妻在线中文字幕| 性色avwww在线观看| 亚洲精品456在线播放app | 又黄又爽又刺激的免费视频.| 激情在线观看视频在线高清| av国产免费在线观看| 亚洲熟妇中文字幕五十中出| 在线观看午夜福利视频| 久久99热6这里只有精品| 国产又黄又爽又无遮挡在线| 九九在线视频观看精品| 五月伊人婷婷丁香| 国产国拍精品亚洲av在线观看| 久久天躁狠狠躁夜夜2o2o| 精品99又大又爽又粗少妇毛片 | 91狼人影院| 97碰自拍视频| 露出奶头的视频| 一级作爱视频免费观看| 在线播放无遮挡| 男女做爰动态图高潮gif福利片| 成人性生交大片免费视频hd| 欧美激情久久久久久爽电影| 噜噜噜噜噜久久久久久91| 欧美成人a在线观看| 一级毛片久久久久久久久女| 亚洲一区高清亚洲精品| 亚洲在线观看片| 变态另类成人亚洲欧美熟女| 成人国产一区最新在线观看| 99久久精品热视频| 哪里可以看免费的av片| 免费大片18禁| 成人鲁丝片一二三区免费| 91在线观看av| 日本五十路高清| 亚洲最大成人av| 日本撒尿小便嘘嘘汇集6| 午夜日韩欧美国产| 婷婷精品国产亚洲av在线| 国产精品久久视频播放| 欧美激情国产日韩精品一区| 国产高清视频在线观看网站| 亚洲精品乱码久久久v下载方式| 国内少妇人妻偷人精品xxx网站| 噜噜噜噜噜久久久久久91| 国产日本99.免费观看| 日日夜夜操网爽| av在线蜜桃| 国产欧美日韩一区二区三| 亚洲成av人片免费观看| 美女被艹到高潮喷水动态| 黄色日韩在线| eeuss影院久久| 国产成人影院久久av| 日韩大尺度精品在线看网址| 99国产极品粉嫩在线观看| 在线免费观看的www视频| 国产精品久久久久久人妻精品电影| 日本撒尿小便嘘嘘汇集6| 亚洲午夜理论影院| 最好的美女福利视频网| 国产精品女同一区二区软件 | 亚洲欧美精品综合久久99| 成人国产综合亚洲| 变态另类丝袜制服| 十八禁国产超污无遮挡网站| 一a级毛片在线观看| 亚洲人与动物交配视频| 成年女人看的毛片在线观看| 日韩av在线大香蕉| 看黄色毛片网站| 九九久久精品国产亚洲av麻豆| 热99re8久久精品国产| 一个人看的www免费观看视频| 亚洲av日韩精品久久久久久密| 99久久成人亚洲精品观看| 最好的美女福利视频网| 最近在线观看免费完整版| 白带黄色成豆腐渣| 亚洲av第一区精品v没综合| 欧美潮喷喷水| 一进一出好大好爽视频| 日本黄色片子视频| 97碰自拍视频| 欧美日韩福利视频一区二区| 欧美日韩黄片免| 久久6这里有精品| 亚洲中文日韩欧美视频| 国产精品一区二区三区四区久久| 一本精品99久久精品77| 少妇丰满av| 无人区码免费观看不卡| 国产精品综合久久久久久久免费| 欧美乱妇无乱码| 国产毛片a区久久久久| a级毛片免费高清观看在线播放| 在线播放国产精品三级| 国产主播在线观看一区二区| 色av中文字幕| 精品人妻视频免费看| 亚洲av免费高清在线观看| 伦理电影大哥的女人| 综合色av麻豆| 99久久精品一区二区三区| 国产91精品成人一区二区三区| 亚洲精品456在线播放app | 久久久久久久午夜电影| 精品久久久久久久久久免费视频| 丰满的人妻完整版| 日韩 亚洲 欧美在线| 日韩国内少妇激情av| 亚洲av美国av| 国产成人av教育| 少妇熟女aⅴ在线视频| av欧美777| 男女床上黄色一级片免费看| 一进一出抽搐动态| 成人永久免费在线观看视频| 一级a爱片免费观看的视频| 日本 欧美在线| 国产探花在线观看一区二区| 精品一区二区免费观看| 久久伊人香网站| 欧美日本亚洲视频在线播放| 很黄的视频免费| 三级国产精品欧美在线观看| 成年人黄色毛片网站| 国产精品1区2区在线观看.| 国产伦人伦偷精品视频| 国产熟女xx| 国产成人福利小说| 国产精品电影一区二区三区| АⅤ资源中文在线天堂| 免费观看精品视频网站| 国产av不卡久久| 国产色爽女视频免费观看| 国产主播在线观看一区二区| 中文字幕免费在线视频6| av在线观看视频网站免费| 亚洲人成电影免费在线| 我的老师免费观看完整版| 色精品久久人妻99蜜桃| 欧美最新免费一区二区三区 | 国产欧美日韩精品一区二区| 男女床上黄色一级片免费看| 亚洲在线自拍视频| 97超视频在线观看视频| 国产 一区 欧美 日韩| 草草在线视频免费看| 午夜精品久久久久久毛片777| 亚洲成av人片在线播放无| 丰满的人妻完整版| 亚洲国产色片| 亚洲va日本ⅴa欧美va伊人久久| 亚洲激情在线av| 美女高潮的动态| 日韩中文字幕欧美一区二区| 亚洲最大成人av| 美女 人体艺术 gogo| 91狼人影院| 高清日韩中文字幕在线| 午夜激情福利司机影院| 国产毛片a区久久久久| 十八禁网站免费在线| 性欧美人与动物交配| 国产午夜福利久久久久久| 国产精品久久久久久人妻精品电影| 免费看光身美女| 日韩欧美免费精品| 亚洲,欧美,日韩| or卡值多少钱| 国产熟女xx| 亚洲最大成人av| 国产在线男女| 欧美一区二区国产精品久久精品| 日日夜夜操网爽| 免费高清视频大片| 淫妇啪啪啪对白视频| 不卡一级毛片| 日本a在线网址| 99久久99久久久精品蜜桃| 有码 亚洲区| 如何舔出高潮| 午夜福利高清视频| 99热6这里只有精品| 好男人电影高清在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲国产高清在线一区二区三| 亚洲精品久久国产高清桃花| 日本成人三级电影网站| 18禁裸乳无遮挡免费网站照片| 精品福利观看| 91在线精品国自产拍蜜月| 欧美日韩国产亚洲二区| 性色av乱码一区二区三区2| 亚洲国产日韩欧美精品在线观看| av在线天堂中文字幕| 免费av观看视频| 亚洲最大成人av| 日韩欧美在线乱码| 国产精品精品国产色婷婷| 9191精品国产免费久久| 我要搜黄色片| 一个人看视频在线观看www免费| 亚洲性夜色夜夜综合| 91在线观看av| 性色av乱码一区二区三区2| 熟妇人妻久久中文字幕3abv| 性色av乱码一区二区三区2| 一个人看视频在线观看www免费| 我要搜黄色片| 婷婷亚洲欧美| 日韩高清综合在线| 麻豆成人av在线观看| 狠狠狠狠99中文字幕| a级一级毛片免费在线观看| 午夜福利欧美成人| 国产亚洲精品久久久com| 国内精品久久久久精免费| 午夜影院日韩av| av女优亚洲男人天堂| 成年版毛片免费区| 亚洲最大成人av| 亚洲av成人精品一区久久| 五月伊人婷婷丁香| 国内精品久久久久久久电影| 午夜两性在线视频| 嫩草影院新地址| 看免费av毛片| 亚洲美女视频黄频| 久99久视频精品免费| 国产精品1区2区在线观看.| 欧美性猛交╳xxx乱大交人| 久久久久久久亚洲中文字幕 | 一级黄片播放器| 91久久精品电影网| 性色avwww在线观看| 国产毛片a区久久久久| 欧美性感艳星| 亚洲午夜理论影院| 国产三级黄色录像| 久久国产精品人妻蜜桃| 黄色女人牲交| 成年女人毛片免费观看观看9| 好男人电影高清在线观看| 国模一区二区三区四区视频| 欧美+日韩+精品| or卡值多少钱| 久久久精品大字幕| 欧美黄色片欧美黄色片| 亚洲精品乱码久久久v下载方式| 美女大奶头视频| 久久99热6这里只有精品| 老熟妇仑乱视频hdxx| 一边摸一边抽搐一进一小说| 欧美激情国产日韩精品一区| 国产亚洲欧美在线一区二区| 亚洲国产精品久久男人天堂| 老熟妇仑乱视频hdxx| 午夜福利视频1000在线观看| 国产毛片a区久久久久| 99国产综合亚洲精品| 免费观看精品视频网站| 99久久精品热视频| 人妻久久中文字幕网| 2021天堂中文幕一二区在线观| 成人av在线播放网站| 动漫黄色视频在线观看| 在线观看美女被高潮喷水网站 | 国产精品综合久久久久久久免费| 97热精品久久久久久| 又黄又爽又免费观看的视频| 嫩草影院精品99| 国产精品野战在线观看| 亚洲av.av天堂| 精品99又大又爽又粗少妇毛片 | 亚洲精品成人久久久久久| 欧美乱妇无乱码| 免费无遮挡裸体视频| 狠狠狠狠99中文字幕| 亚洲一区二区三区不卡视频| 变态另类丝袜制服| 淫妇啪啪啪对白视频| 久久精品国产亚洲av香蕉五月| 成人毛片a级毛片在线播放| 搞女人的毛片| 国产一区二区亚洲精品在线观看| 黄色丝袜av网址大全| 精品人妻视频免费看| 热99re8久久精品国产| 国产高清视频在线播放一区| 在线观看午夜福利视频| 国产探花极品一区二区| 欧美性猛交黑人性爽| 午夜精品一区二区三区免费看| 免费电影在线观看免费观看| 国内久久婷婷六月综合欲色啪| 搡老岳熟女国产| 国内精品久久久久久久电影| 亚洲中文字幕日韩| 亚洲av第一区精品v没综合| 国产午夜精品论理片| 99在线人妻在线中文字幕| 国产在线精品亚洲第一网站| 九九久久精品国产亚洲av麻豆| 成熟少妇高潮喷水视频| 精品久久久久久,| 亚洲人与动物交配视频| 每晚都被弄得嗷嗷叫到高潮| 青草久久国产| 亚洲人成网站在线播放欧美日韩| 欧美黑人巨大hd| 99久久九九国产精品国产免费| 天美传媒精品一区二区| 好看av亚洲va欧美ⅴa在| 午夜福利免费观看在线| 日日摸夜夜添夜夜添小说| 天堂动漫精品| 日韩大尺度精品在线看网址| 十八禁网站免费在线| 熟女人妻精品中文字幕| 日韩精品中文字幕看吧| 国产高清有码在线观看视频| 熟妇人妻久久中文字幕3abv| 18禁在线播放成人免费| 免费av观看视频| 90打野战视频偷拍视频| 天堂av国产一区二区熟女人妻| 国产av不卡久久| 午夜精品一区二区三区免费看| 91麻豆精品激情在线观看国产| 久久伊人香网站| 婷婷丁香在线五月| 日韩人妻高清精品专区| www.www免费av| 成年版毛片免费区| 搡老妇女老女人老熟妇| 国产 一区 欧美 日韩| 亚洲片人在线观看| 国产69精品久久久久777片| 人人妻人人澡欧美一区二区| 午夜激情福利司机影院| 欧美又色又爽又黄视频| 非洲黑人性xxxx精品又粗又长| 欧美+亚洲+日韩+国产| 精品熟女少妇八av免费久了| 免费在线观看影片大全网站| 亚洲精品影视一区二区三区av| 欧美三级亚洲精品| 无人区码免费观看不卡| 男女做爰动态图高潮gif福利片| 国产av麻豆久久久久久久| 51午夜福利影视在线观看| 成人特级黄色片久久久久久久| 99热这里只有是精品在线观看 | 亚洲人成网站在线播| 国产成人福利小说| 欧美不卡视频在线免费观看| 欧美色视频一区免费| 午夜福利18| 亚洲人成网站高清观看| 制服丝袜大香蕉在线| 国产极品精品免费视频能看的| 成人鲁丝片一二三区免费| 亚洲av美国av| 免费看光身美女| 18禁裸乳无遮挡免费网站照片| 日本黄大片高清| 亚洲成人久久爱视频| 成人av一区二区三区在线看| 男人和女人高潮做爰伦理| 日韩中文字幕欧美一区二区| 一级作爱视频免费观看| 不卡一级毛片| 国产淫片久久久久久久久 | 亚洲国产欧美人成| 搡老熟女国产l中国老女人| 日韩欧美免费精品| 婷婷亚洲欧美| 国产精品久久视频播放| 99视频精品全部免费 在线| 久久九九热精品免费| 不卡一级毛片| 精品无人区乱码1区二区| 无遮挡黄片免费观看| 国产久久久一区二区三区| 男女床上黄色一级片免费看| 免费一级毛片在线播放高清视频| 免费看美女性在线毛片视频| 国产精品一区二区三区四区免费观看 | 国产中年淑女户外野战色| 欧美又色又爽又黄视频| 亚洲一区高清亚洲精品| eeuss影院久久| 又爽又黄a免费视频| 97超视频在线观看视频| 国产色婷婷99| 国产亚洲欧美在线一区二区| 免费av观看视频| 国内毛片毛片毛片毛片毛片| 熟妇人妻久久中文字幕3abv| 久久6这里有精品| 亚洲专区国产一区二区| 色精品久久人妻99蜜桃| 国产老妇女一区| 99在线视频只有这里精品首页| 亚洲专区中文字幕在线| 欧美黄色淫秽网站| 一a级毛片在线观看| 乱码一卡2卡4卡精品| 亚洲熟妇熟女久久| 人人妻人人看人人澡| 两个人视频免费观看高清| 三级男女做爰猛烈吃奶摸视频| 国产淫片久久久久久久久 | 亚洲欧美日韩卡通动漫| 久久国产精品影院| 久久久国产成人精品二区| 国产精品一区二区三区四区久久| 老熟妇乱子伦视频在线观看| 51国产日韩欧美| 国产私拍福利视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产精品美女特级片免费视频播放器| 免费av毛片视频| 99国产精品一区二区三区| 国产精品美女特级片免费视频播放器| 一进一出好大好爽视频| 久久精品国产亚洲av涩爱 | 五月伊人婷婷丁香| 亚洲中文日韩欧美视频| 国产精华一区二区三区| avwww免费| 一本综合久久免费| 一级作爱视频免费观看| 亚洲激情在线av|