• <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久久| 丁香六月天网| 亚洲国产欧美日韩在线播放| √禁漫天堂资源中文www| 伦理电影免费视频| 卡戴珊不雅视频在线播放| 欧美精品亚洲一区二区| 亚洲图色成人| 在现免费观看毛片| 久久久a久久爽久久v久久| 美女午夜性视频免费| 婷婷成人精品国产| 涩涩av久久男人的天堂| 欧美97在线视频| 黄色一级大片看看| 黑丝袜美女国产一区| 99精国产麻豆久久婷婷| 日韩精品有码人妻一区| 高清av免费在线| 午夜福利影视在线免费观看| 少妇 在线观看| 国产精品久久久久久精品古装| 少妇的逼水好多| 黄色一级大片看看| 午夜福利在线观看免费完整高清在| 国产一区二区在线观看av| 国产成人精品婷婷| 欧美日韩av久久| 狠狠精品人妻久久久久久综合| 尾随美女入室| 高清黄色对白视频在线免费看| 日产精品乱码卡一卡2卡三| 欧美国产精品一级二级三级| 丰满少妇做爰视频| tube8黄色片| 18禁国产床啪视频网站| 寂寞人妻少妇视频99o| 久久久久久伊人网av| 2018国产大陆天天弄谢| 最近的中文字幕免费完整| 国产深夜福利视频在线观看| 国产成人午夜福利电影在线观看| 大香蕉久久成人网| 久久久国产一区二区| 99久国产av精品国产电影| 亚洲精品在线美女| 国产伦理片在线播放av一区| 国产成人精品无人区| 国产一区二区三区av在线| 亚洲国产精品一区三区| 国产片特级美女逼逼视频| 国产日韩欧美在线精品| 90打野战视频偷拍视频| av片东京热男人的天堂| 天美传媒精品一区二区| 王馨瑶露胸无遮挡在线观看| 热re99久久精品国产66热6| 一区二区日韩欧美中文字幕| 99精国产麻豆久久婷婷| 亚洲欧洲国产日韩| 综合色丁香网| 亚洲四区av| 国产日韩欧美亚洲二区| 国产成人一区二区在线| 青草久久国产| 97在线视频观看| a级毛片黄视频| 中文字幕亚洲精品专区| 欧美变态另类bdsm刘玥| 视频在线观看一区二区三区| 久久久久国产一级毛片高清牌| 我要看黄色一级片免费的| 黄色 视频免费看| 亚洲在久久综合| 国产成人精品久久二区二区91 | 国产成人精品久久久久久| 街头女战士在线观看网站| 久久影院123| 中文字幕精品免费在线观看视频| 亚洲四区av| 日韩电影二区| 你懂的网址亚洲精品在线观看| 一级毛片 在线播放| 两个人看的免费小视频| av电影中文网址| 亚洲精品国产av蜜桃| 日韩欧美一区视频在线观看| 一级片免费观看大全| 看非洲黑人一级黄片| 青草久久国产| 国产成人av激情在线播放| 五月天丁香电影| 一区福利在线观看| 亚洲成色77777| 男女边摸边吃奶| 欧美精品av麻豆av| 夜夜骑夜夜射夜夜干| 欧美av亚洲av综合av国产av | 免费看不卡的av| 视频在线观看一区二区三区| 少妇的丰满在线观看| 超碰成人久久| 久久久久久人人人人人| 亚洲,欧美,日韩| 成人毛片a级毛片在线播放| 少妇人妻 视频| 国产男女内射视频| 狂野欧美激情性bbbbbb| 免费黄频网站在线观看国产| 国产在线视频一区二区| 亚洲精品aⅴ在线观看| 久久久久视频综合| 黄色毛片三级朝国网站| 久久久久国产精品人妻一区二区| 波多野结衣av一区二区av| 亚洲精品美女久久av网站| 中文字幕人妻熟女乱码| 成人影院久久| www日本在线高清视频| 国产精品麻豆人妻色哟哟久久| 国产xxxxx性猛交| 国产在线免费精品| 夫妻性生交免费视频一级片| 黄片无遮挡物在线观看| 91在线精品国自产拍蜜月| 一二三四在线观看免费中文在| 午夜福利乱码中文字幕| 在线看a的网站| 国产成人91sexporn| 欧美人与性动交α欧美软件| 最新的欧美精品一区二区| 日韩av免费高清视频| 日韩免费高清中文字幕av| 三上悠亚av全集在线观看| 欧美最新免费一区二区三区| 哪个播放器可以免费观看大片| 男女国产视频网站| 超碰成人久久| 免费看不卡的av| 纵有疾风起免费观看全集完整版| 97在线视频观看| 色婷婷久久久亚洲欧美| 999久久久国产精品视频| 欧美日韩综合久久久久久| 深夜精品福利| 国产成人一区二区在线| 国产欧美日韩综合在线一区二区| 国产精品 欧美亚洲| 最近手机中文字幕大全| 久久婷婷青草| 亚洲成色77777| 久久国产精品大桥未久av| 欧美bdsm另类| 夫妻午夜视频| 久热久热在线精品观看| 1024香蕉在线观看| av免费观看日本| 黄片播放在线免费| 亚洲四区av| 国产日韩欧美在线精品| 久久久久久久久久久免费av| 亚洲久久久国产精品| 日韩熟女老妇一区二区性免费视频| 色婷婷av一区二区三区视频| 日本欧美国产在线视频| 青春草视频在线免费观看| 欧美亚洲 丝袜 人妻 在线| 久久久久国产一级毛片高清牌| 亚洲精品成人av观看孕妇| 一区福利在线观看| av免费观看日本| 黄色怎么调成土黄色| 乱人伦中国视频| 在线天堂中文资源库| 超色免费av| 在线观看免费高清a一片| 亚洲精品中文字幕在线视频| 亚洲av电影在线进入| 欧美人与性动交α欧美精品济南到 | 亚洲av综合色区一区| 久久这里有精品视频免费| 国产女主播在线喷水免费视频网站| 如日韩欧美国产精品一区二区三区| 久久久久久久久免费视频了| 天堂俺去俺来也www色官网| 女人久久www免费人成看片| 哪个播放器可以免费观看大片| 一区二区三区乱码不卡18| 久久久a久久爽久久v久久| 精品国产一区二区三区久久久樱花| 大片电影免费在线观看免费| 亚洲精品国产av成人精品| 另类亚洲欧美激情| 精品人妻偷拍中文字幕| 成人国产av品久久久| 搡老乐熟女国产| 一本久久精品| 在线观看人妻少妇| 日本免费在线观看一区| 亚洲精品乱久久久久久| 国产成人精品婷婷| 国产精品99久久99久久久不卡 | 日韩,欧美,国产一区二区三区| 国产精品成人在线| 七月丁香在线播放| 亚洲欧美精品综合一区二区三区 | 成人国语在线视频| 久久午夜福利片| 国产av一区二区精品久久| 亚洲国产精品一区二区三区在线| 大码成人一级视频| 人人妻人人添人人爽欧美一区卜| 一二三四在线观看免费中文在| 亚洲伊人色综图| 波野结衣二区三区在线| 亚洲综合色惰| 亚洲五月色婷婷综合| 国产精品不卡视频一区二区| 婷婷色综合大香蕉| 卡戴珊不雅视频在线播放| 久久久久国产精品人妻一区二区| 人妻 亚洲 视频| 日韩一本色道免费dvd| 高清不卡的av网站| 国产爽快片一区二区三区| 人妻一区二区av| 欧美bdsm另类| 亚洲 欧美一区二区三区| 天天影视国产精品| 美女大奶头黄色视频| 亚洲精品av麻豆狂野| 天天躁夜夜躁狠狠躁躁| 亚洲精品aⅴ在线观看| 国产精品三级大全| av女优亚洲男人天堂| 久久鲁丝午夜福利片| 免费在线观看黄色视频的| 国产成人aa在线观看| 亚洲精品一二三| 国产精品亚洲av一区麻豆 | 欧美少妇被猛烈插入视频| 亚洲精品成人av观看孕妇| 国产毛片在线视频| 久久久久视频综合| av国产精品久久久久影院| 午夜福利视频在线观看免费| 成人手机av| 观看美女的网站| 巨乳人妻的诱惑在线观看| a级毛片在线看网站| videossex国产| 涩涩av久久男人的天堂| 久久国产精品大桥未久av| 亚洲av福利一区| 我要看黄色一级片免费的| 亚洲国产日韩一区二区| 视频区图区小说| 亚洲三区欧美一区| 不卡视频在线观看欧美| √禁漫天堂资源中文www| 国产无遮挡羞羞视频在线观看| 免费女性裸体啪啪无遮挡网站| 自拍欧美九色日韩亚洲蝌蚪91| 性色avwww在线观看| 深夜精品福利| 少妇人妻 视频| 建设人人有责人人尽责人人享有的| 欧美+日韩+精品| 日本wwww免费看| 一二三四在线观看免费中文在| 亚洲精品美女久久久久99蜜臀 | 99国产精品免费福利视频| 欧美日韩视频精品一区| 亚洲欧美一区二区三区国产| 宅男免费午夜| 日产精品乱码卡一卡2卡三| 免费日韩欧美在线观看| 久久久精品免费免费高清| 性少妇av在线| 亚洲国产色片| 国产成人精品在线电影| 亚洲精品乱久久久久久| 我要看黄色一级片免费的| 免费看不卡的av| 久久精品亚洲av国产电影网| 丝袜脚勾引网站| 国产精品不卡视频一区二区| 久久久国产一区二区| 丰满饥渴人妻一区二区三| 最近最新中文字幕大全免费视频 | freevideosex欧美| av一本久久久久| 欧美 日韩 精品 国产| 精品国产露脸久久av麻豆| 人人澡人人妻人| 国产在线一区二区三区精| 成年美女黄网站色视频大全免费| 黄色毛片三级朝国网站| 青春草视频在线免费观看| 在线观看免费日韩欧美大片| av电影中文网址| 尾随美女入室| 美女主播在线视频| 黄网站色视频无遮挡免费观看| 亚洲欧美成人精品一区二区| 日韩伦理黄色片| av天堂久久9| 亚洲精品国产av成人精品| 成人二区视频| 免费少妇av软件| 久久鲁丝午夜福利片| 免费人妻精品一区二区三区视频| 亚洲欧美一区二区三区久久| 美女中出高潮动态图| 伊人亚洲综合成人网| 一本久久精品| 熟女av电影| 亚洲内射少妇av| 黄色 视频免费看| 欧美亚洲 丝袜 人妻 在线| 午夜福利视频精品| 国产精品二区激情视频| 国产97色在线日韩免费| 国产深夜福利视频在线观看| 男女边摸边吃奶| 丝袜脚勾引网站| 国产欧美亚洲国产| 亚洲国产精品一区三区| 高清在线视频一区二区三区| 99香蕉大伊视频| www日本在线高清视频| 老熟女久久久| 欧美日韩精品成人综合77777| 亚洲婷婷狠狠爱综合网| 亚洲国产日韩一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 日本91视频免费播放| 亚洲精品视频女| 成年人午夜在线观看视频| 捣出白浆h1v1| 国产一区二区在线观看av| 欧美精品亚洲一区二区| 成人亚洲精品一区在线观看| 久久99热这里只频精品6学生| 1024香蕉在线观看| 精品一区在线观看国产| 久久久久久人妻| 免费观看a级毛片全部| 老司机亚洲免费影院| 精品国产一区二区三区久久久樱花| 国产视频首页在线观看| 日韩欧美一区视频在线观看| 性少妇av在线| 天天躁夜夜躁狠狠躁躁| 亚洲一区二区三区欧美精品| a 毛片基地| 久久女婷五月综合色啪小说| 啦啦啦视频在线资源免费观看| 青春草视频在线免费观看| 亚洲欧美一区二区三区黑人 | 久久精品aⅴ一区二区三区四区 | 国产免费一区二区三区四区乱码| 国产精品.久久久| 国产精品免费大片| 久久国产精品男人的天堂亚洲| 黄片无遮挡物在线观看| 国产无遮挡羞羞视频在线观看| 黄频高清免费视频| 超碰成人久久| 在线天堂中文资源库| 免费观看性生交大片5| 美女中出高潮动态图| 熟女少妇亚洲综合色aaa.| 男人舔女人的私密视频| 精品国产乱码久久久久久小说| 欧美日韩精品网址| 91精品三级在线观看| 有码 亚洲区| 天天躁狠狠躁夜夜躁狠狠躁| 哪个播放器可以免费观看大片| 久久99热这里只频精品6学生| 人妻一区二区av| 少妇人妻久久综合中文| 婷婷色av中文字幕| 免费女性裸体啪啪无遮挡网站| 国产成人av激情在线播放| 亚洲第一av免费看| 久久久久精品久久久久真实原创| 免费在线观看视频国产中文字幕亚洲 | 永久免费av网站大全| 春色校园在线视频观看| 国产一区二区激情短视频 | 久久久久国产精品人妻一区二区| 纵有疾风起免费观看全集完整版| 成人国产麻豆网| 欧美成人精品欧美一级黄| 国产福利在线免费观看视频| 亚洲欧美中文字幕日韩二区| 看免费成人av毛片| 国产av国产精品国产| 赤兔流量卡办理| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产av影院在线观看| 欧美日韩精品网址| av.在线天堂| 高清不卡的av网站| 亚洲,欧美精品.| 黄色怎么调成土黄色| 日韩欧美一区视频在线观看| 伦理电影免费视频| 国产日韩欧美在线精品| 少妇的丰满在线观看| 亚洲精品国产av蜜桃| 精品一区二区三区四区五区乱码 | 亚洲天堂av无毛| 国产精品久久久久久精品电影小说| 电影成人av| 9191精品国产免费久久| 国产精品av久久久久免费| 我的亚洲天堂| 看十八女毛片水多多多| 中国国产av一级| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 免费黄色在线免费观看| 成人亚洲欧美一区二区av| 久久97久久精品| 午夜福利,免费看| 久久午夜综合久久蜜桃| 国产午夜精品一二区理论片| 国产精品欧美亚洲77777| 老女人水多毛片| 少妇的丰满在线观看| 亚洲精品一区蜜桃| 99久久中文字幕三级久久日本| 免费黄网站久久成人精品| av免费在线看不卡| 久久久久久人妻| 成人国语在线视频| 九九爱精品视频在线观看| 日韩一卡2卡3卡4卡2021年| 黄色 视频免费看| 国产一区二区在线观看av| 丝袜在线中文字幕| 亚洲美女视频黄频| 女人被躁到高潮嗷嗷叫费观| 色婷婷av一区二区三区视频| 久久久久久免费高清国产稀缺| 久久精品久久久久久噜噜老黄| 一区二区三区激情视频| 欧美人与性动交α欧美精品济南到 | 亚洲国产毛片av蜜桃av| 日韩av免费高清视频| 女的被弄到高潮叫床怎么办| 日日撸夜夜添| 久久精品夜色国产| 天天躁夜夜躁狠狠久久av| 亚洲欧美日韩另类电影网站| 男女下面插进去视频免费观看| 久久人人爽人人片av| 91午夜精品亚洲一区二区三区| 青春草国产在线视频| 纵有疾风起免费观看全集完整版| 天堂俺去俺来也www色官网| 日韩av在线免费看完整版不卡| 成人二区视频| av在线观看视频网站免费| 男女啪啪激烈高潮av片| 亚洲国产精品国产精品| 国产探花极品一区二区| 国产一区二区在线观看av| 国产成人av激情在线播放| 久久久久国产精品人妻一区二区| 国产高清不卡午夜福利| 一本—道久久a久久精品蜜桃钙片| 亚洲成人一二三区av| 日本色播在线视频| 美女国产高潮福利片在线看| 亚洲av日韩在线播放| 午夜激情久久久久久久| 中文字幕制服av| 日本av免费视频播放| 咕卡用的链子| 亚洲人成网站在线观看播放| 日本wwww免费看| 精品少妇一区二区三区视频日本电影 | 国产视频首页在线观看| 女人被躁到高潮嗷嗷叫费观| 免费高清在线观看日韩| 成人国语在线视频| 亚洲美女搞黄在线观看| 天天操日日干夜夜撸| 成人漫画全彩无遮挡| h视频一区二区三区| 国产熟女午夜一区二区三区| 亚洲四区av| 精品国产乱码久久久久久男人| 欧美精品av麻豆av| 久久精品国产亚洲av天美| 五月伊人婷婷丁香| 99久国产av精品国产电影| 青春草视频在线免费观看| 日韩中文字幕视频在线看片| 久久久久久免费高清国产稀缺| 亚洲精品久久久久久婷婷小说| 99九九在线精品视频| 秋霞伦理黄片| xxx大片免费视频| 在线观看www视频免费| 亚洲在久久综合| 久久久久久伊人网av| 老司机亚洲免费影院| 美女高潮到喷水免费观看| 一本色道久久久久久精品综合| 久久精品夜色国产| 2022亚洲国产成人精品| 麻豆精品久久久久久蜜桃| 人人妻人人澡人人看| 巨乳人妻的诱惑在线观看| 在线观看美女被高潮喷水网站| 黑丝袜美女国产一区| 最近的中文字幕免费完整| 人人妻人人澡人人看| 一区福利在线观看| 91精品国产国语对白视频| 亚洲av日韩在线播放| 国产精品久久久久久精品古装| 各种免费的搞黄视频| 高清视频免费观看一区二区| 国产精品久久久久成人av| 91aial.com中文字幕在线观看| 免费黄频网站在线观看国产| 人人澡人人妻人| 免费观看av网站的网址| 一级黄片播放器| 国产精品一国产av| 人妻少妇偷人精品九色| 国产一级毛片在线| av在线观看视频网站免费| 啦啦啦啦在线视频资源| 视频在线观看一区二区三区| 啦啦啦啦在线视频资源| 久久这里只有精品19| 亚洲三级黄色毛片| 青春草亚洲视频在线观看| 看免费成人av毛片| 精品第一国产精品| 精品一品国产午夜福利视频| 男女高潮啪啪啪动态图| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产视频首页在线观看| 99re6热这里在线精品视频| 色视频在线一区二区三区| 久久av网站| 丰满迷人的少妇在线观看| 亚洲国产精品国产精品| 亚洲av男天堂| 国产无遮挡羞羞视频在线观看| 精品国产乱码久久久久久男人| 国产精品久久久久成人av| 精品酒店卫生间| 精品卡一卡二卡四卡免费| 热99久久久久精品小说推荐| 国产成人免费观看mmmm| 999精品在线视频| 久久久国产精品麻豆| 热re99久久国产66热| 中文字幕人妻丝袜一区二区 | 免费女性裸体啪啪无遮挡网站| 亚洲av免费高清在线观看| 看十八女毛片水多多多| 天天躁日日躁夜夜躁夜夜| 伦精品一区二区三区| 国产精品免费大片| 亚洲精华国产精华液的使用体验| xxx大片免费视频| kizo精华| 精品少妇久久久久久888优播| 国产女主播在线喷水免费视频网站| 国产精品99久久99久久久不卡 | 精品人妻在线不人妻| 色哟哟·www| 1024香蕉在线观看| 国产精品女同一区二区软件| 午夜福利网站1000一区二区三区| av不卡在线播放| 久久精品久久久久久噜噜老黄| 久久久精品免费免费高清| 午夜免费鲁丝| 电影成人av| 久久精品久久久久久久性| 又粗又硬又长又爽又黄的视频| 精品亚洲成国产av| av在线播放精品| 交换朋友夫妻互换小说| 啦啦啦啦在线视频资源| 国产男女内射视频| 久久久欧美国产精品| 欧美 亚洲 国产 日韩一| 五月伊人婷婷丁香| 老女人水多毛片| 啦啦啦在线免费观看视频4| 国产精品.久久久| 建设人人有责人人尽责人人享有的| 中文字幕色久视频| 亚洲精华国产精华液的使用体验| 2021少妇久久久久久久久久久| 你懂的网址亚洲精品在线观看| 国产欧美日韩一区二区三区在线| 国产精品.久久久| 日本vs欧美在线观看视频| 欧美人与善性xxx| 在线观看免费高清a一片| 精品国产乱码久久久久久男人| 一级a爱视频在线免费观看| 汤姆久久久久久久影院中文字幕|