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

    七星湖濕地磷素分布特征及擴(kuò)散通量估算

    2017-03-15 09:27:56鐘馨田永蘭鐘美芳張化永
    關(guān)鍵詞:明水磷素七星

    鐘馨,田永蘭,鐘美芳,張化永

    華北電力大學(xué)工程生態(tài)學(xué)與非線性科學(xué)研究中心,北京 102206

    七星湖濕地磷素分布特征及擴(kuò)散通量估算

    鐘馨,田永蘭,鐘美芳,張化永*

    華北電力大學(xué)工程生態(tài)學(xué)與非線性科學(xué)研究中心,北京 102206

    為研究山東省滕州市七星湖濕地磷素的分布特征及其遷移規(guī)律,植被在磷的遷移轉(zhuǎn)化中的作用,采集秋季七星湖表層沉積物及其上覆水體樣品,分析總磷(TP)、總?cè)芙庑粤?TDP)、溶解性活性反應(yīng)磷(SRP)、顆粒態(tài)磷(PP)的分布特征和相關(guān)性,計(jì)算沉積物-上覆水體界面的擴(kuò)散通量。結(jié)果表明:TP濃度在上覆水體明水區(qū)域沿途總體變化不大(0.05~0.06 mg/L),在蘆葦群落沿程降低(0.04~0.09 mg/L),香蒲群落中有波動(dòng)變化(0.02~0.09 mg/L)。上、中、下游3個(gè)斷面明水區(qū)域沉積物間隙水中各形態(tài)磷素濃度均比植物區(qū)域高,說明植物對沉積物中磷素的截留吸收轉(zhuǎn)化有重要作用。上、中、下游大型挺水植物各器官磷素所占比例有明顯差異,且地上部分高于地下部分。七星湖濕地經(jīng)過3年的磷素沉降,沉積物TP濃度增加42.2%。沉積物-上覆水體間TDP和SRP的擴(kuò)散通量均為正通量,表明磷素從沉積物向上覆水體擴(kuò)散,且明水區(qū)域擴(kuò)散通量均高于植物區(qū)域。表層沉積物間隙水和上覆水體各形態(tài)磷的相關(guān)性表明,間隙水中的磷主要以TDP和SRP的形態(tài)向上覆水體擴(kuò)散。

    間隙水;濕地;水生植物;磷形態(tài);擴(kuò)散通量

    濕地科學(xué)研究受到越來越多的關(guān)注,已成為重點(diǎn)學(xué)科和重要的研究領(lǐng)域。天然濕地常由于自身的演變,以及人類活動(dòng)的干預(yù),發(fā)生了結(jié)構(gòu)和功能的退化,逐漸失去作為濕地的特征或直接轉(zhuǎn)變成其他類型用地。人工濕地污水處理系統(tǒng)在長期運(yùn)行過程中存在著彈性小、抗沖擊能力差、水質(zhì)凈化效果不穩(wěn)定、壽命短等現(xiàn)象[1]。因此一些研究者提出了近自然濕地的工程手段,即根據(jù)生態(tài)學(xué)原理,結(jié)合具體地形地貌、水文和植被情況,模仿、接近自然濕地的一種可持續(xù)的、具有生物多樣性的生態(tài)工程技術(shù)。

    濕地中磷素的存在形態(tài)主要包括有機(jī)磷(生物態(tài)和非生物態(tài)的)、磷酸、可溶性磷酸鹽和不溶性磷酸鹽,磷素的主要去除機(jī)理包括化學(xué)作用(如沉淀作用和吸附作用)、生物作用(如植物吸收作用和微生物吸收與積累作用)和物理作用(如沉積作用)[2]。沉積物-上覆水體界面是水環(huán)境的一個(gè)特殊而重要的區(qū)域。磷素在沉積物-上覆水體界面的釋放和擴(kuò)散量是營養(yǎng)鹽在濕地中循環(huán)的重要研究內(nèi)容。沉積物中的磷素主要以間隙水作為介質(zhì),通過表面擴(kuò)散層向上覆水體擴(kuò)散遷移,從而影響上覆水體的磷濃度。在水生植物群落,植物通過吸收上覆水中的磷[3],或者通過根系吸收沉積物中的生物有效磷[4],或者通過莖葉截留吸收部分磷[5],或者通過自身生長變化改變水體理化性質(zhì)等方式影響磷的遷移轉(zhuǎn)化[6]。

    近年來對于海灣、河流和湖泊等區(qū)域沉積物間隙水和上覆水磷素分布特征的研究較多:如宋金明[7]對南沙群島沉積物間隙水中N、P、Si在沉積物-海水界面間的擴(kuò)散通量進(jìn)行研究;張彥等[8]對太湖沉積物及其間隙水中氮的時(shí)空分布以及垂向和水平方向的變化進(jìn)行了研究;楊楠等[9]通過研究白塔堡河干流平水期上覆水和沉積物間隙水N、P的分布特征,揭示其營養(yǎng)鹽的遷移規(guī)律。目前對濕地,尤其是近自然濕地,沉積物間隙水和上覆水磷素分布特征的研究較少,而近自然濕地作為具有巨大的環(huán)境功能與效益的生態(tài)系統(tǒng),探索其沉積物-上覆水磷素遷移轉(zhuǎn)化對于揭示濕地凈化機(jī)理方面意義重大。為全面了解近自然濕地磷素分布特征及其遷移規(guī)律,在山東省滕州市七星湖,在上、中、下游斷面分別設(shè)明水區(qū)域、香蒲群落、蘆葦群落,采集秋季表層沉積物及其上覆水體樣品,考察各形態(tài)磷的分布特征,計(jì)算沉積物-上覆水體界面的擴(kuò)散通量,及間隙水和上覆水體中各形態(tài)磷之間的相關(guān)性,分析近自然濕地磷素的遷移轉(zhuǎn)化規(guī)律,以期為近自然濕地的管理和機(jī)制研究提供依據(jù)。

    1 研究區(qū)域與方法

    1.1 試驗(yàn)區(qū)域

    七星湖濕地位于山東省滕州市城河南岸,地處暖溫帶半濕潤地區(qū)南部,季風(fēng)型大陸性氣候明顯,四季分明,雨量充沛,光照充足。濕地由煤礦塌陷區(qū)改造而成,包括5個(gè)子庫塘,通過涵管連通,總規(guī)模約80 hm2??傮w上,七星湖濕地地形多變、景觀結(jié)構(gòu)多樣、水流復(fù)雜。濕地是南四湖流域“截蓄導(dǎo)”生態(tài)工程的一部分,2009年后通過涵管將城河水引入濕地進(jìn)行水質(zhì)凈化。從2009年起,在城河上建造橡膠壩,抬升水位,開始向濕地引入城河河水,濕地中還建有生態(tài)過水堰、人工浮島等多項(xiàng)生態(tài)技術(shù),保證凈化城河入水的同時(shí)也能保障濕地的生態(tài)安全,強(qiáng)化濕地水體凈化效果。

    研究區(qū)域選取七星湖濕地北段,東西寬約1 200 m,南北長約680 m,面積為0.2 km2,水流方向自東向西,經(jīng)末端涵洞向南流出。該庫塘主要挺水植物有蘆葦、香蒲、荷花等;浮葉植物有水花生、菱角、槐葉萍、荇菜等;沉水植物有狐尾藻、苦草、菹草等。濕地進(jìn)水口水質(zhì)如表1所示。

    表1 濕地進(jìn)水口水質(zhì)

    1.2 試驗(yàn)布點(diǎn)及方法

    沿濕地內(nèi)水流方向?qū)竦胤譃樯嫌巍⒅杏?、下?個(gè)斷面,每個(gè)斷面分別布設(shè)香蒲、蘆葦、明水(無植物)3個(gè)采樣點(diǎn)(圖1)。上游斷面水域面積較大,且水較深,水體的流速最?。恢杏螖嗝嬗捎诘匦渭爸参锶郝涞纳L,使得河道較狹窄,采樣點(diǎn)離植物群落較近,水流速度相對較快;下游斷面周圍有大量的植物,水流較湍急。

    圖1 研究區(qū)域采樣點(diǎn)分布Fig.1 Sample sites in the near natural wetland

    于2015年10月底分別在明水區(qū)域、香蒲群落和蘆葦群落中采集水樣(水面下30 cm)和對應(yīng)位置的底泥,將樣品放于低溫箱中保存,并盡快帶至實(shí)驗(yàn)室分析測試。水樣pH、溫度、DO濃度等指標(biāo)利用多功能水質(zhì)儀(YSI professional plus)現(xiàn)場測定。水深用帶有重錘的皮尺測定。表層水的DO濃度在水面下30 cm處測定。水樣磷素指標(biāo)按照《水和廢水檢測分析方法(4版)》測定。水體中SRP為溶解性活性反應(yīng)磷(solubility reactive phosphorus),即為溶解性正磷酸鹽;顆粒態(tài)磷(particulate phosphrous,PP)濃度由水體中TP與總?cè)芙庑粤?total dissolved phosphorous,TDP)差值計(jì)算得來。采集的底泥樣品在實(shí)驗(yàn)室內(nèi)進(jìn)行恒溫離心(4 000 r/min,20 min)得到間隙水樣,測定方法與上覆水體相同。每個(gè)樣品分別取2個(gè)平行樣。底泥自然風(fēng)干后,過100目(0.15 mm)篩用密封袋封裝以供測定。底泥的測定采用SMT(標(biāo)準(zhǔn)測量和測試)分級提取NaOH-P、HCl-P、無機(jī)磷(IP)、有機(jī)磷(OP)和TP。TP主要由IP和OP組成,其中IP占主要部分;NaOH-P主要包括Fe、Mn、Al等金屬氧化物上吸附包裹的磷,pH、氧化還原電位等環(huán)境因素的改變?nèi)菀讓?dǎo)致這部分磷釋放到水體中,從而被生物利用;HCl-P主要包括被離子結(jié)合的磷,一般不易釋放到水體。每個(gè)底泥樣品分別取3個(gè)平行樣。

    在每個(gè)采樣點(diǎn)分別選取3株長勢相近的典型挺水植物,用水小心洗掉根部泥土,統(tǒng)計(jì)鮮重后用密封袋封裝。帶回實(shí)驗(yàn)室后將每種植物分為根、莖、葉3個(gè)部分,迅速置于105 ℃下殺青10~15 min,然后置于80 ℃烘箱中烘至恒重,測定干重。將烘干后的植物樣品各器官分別用碾磨機(jī)粉碎,過0.25 mm篩后保存。處理后的植物樣品經(jīng)H2SO4-H2O2消煮后,用鉬銻抗比色法(GB 9837—88)測全磷濃度。每個(gè)植物樣品的根、莖、葉分別測3個(gè)平行樣。

    2 結(jié)果與分析

    2.1 上覆水體磷素分布特征

    上覆水體理化性質(zhì)見表2。

    表2 上覆水體理化性質(zhì)

    由表2可以看出,明水區(qū)域的水深遠(yuǎn)大于植物區(qū)域,因此植物群落的水力因素以及擾動(dòng)對上覆水體和間隙水的磷素分布影響較大。明水區(qū)域DO濃度較挺水植物區(qū)域偏高,其可能的原因:1)明水區(qū)域浮水植物生物量比挺水植物區(qū)域大,光合作用明顯,導(dǎo)致明水區(qū)域DO濃度偏高;2)挺水植物區(qū)域枝葉較為茂密,擋住了部分陽光,影響了植物的光合作用。有研究表明[10],DO濃度與pH呈極顯著正相關(guān),明水區(qū)域的pH也明顯高于植物區(qū)域,光合作用強(qiáng)導(dǎo)致CO2消耗大,所以明水區(qū)域pH偏高。由表2可以看出,所有采樣點(diǎn)的pH均高于7,呈弱堿性。水體中H+濃度不宜過高,當(dāng)其超過一定濃度時(shí),便會(huì)對植物產(chǎn)生直接的毒害作用,一般對根的影響最大。但水體中OH-濃度也不宜過高,高pH不利于水生植物的生長,其是菹草在夏季死亡的原因之一[11]。

    上覆水體各形態(tài)磷素的濃度變化如圖2所示。由圖2可以看出,明水區(qū)域TP濃度沿途總體變化不大,而TDP和SRP濃度在中游斷面突然降低,在下游斷面又大幅增加。蘆葦群落各形態(tài)磷濃度的普遍變化趨勢是沿程降低,而香蒲群落各形態(tài)磷濃度變化較大,TDP和SRP濃度在入水時(shí)均不高,但在中游斷面大幅增加??赡苁侵杏嗡蛎娣e較小,人為擾動(dòng)較大,因此導(dǎo)致各形態(tài)磷濃度的大幅增加或降低。在上游斷面,植物區(qū)域的各形態(tài)磷濃度明顯高于明水區(qū)域;在下游斷面,植物區(qū)域的各形態(tài)磷濃度低于明水區(qū)域。其中TDP和SRP在沿程去除效果尤為明顯,這表明植物對磷的去除作用主要表現(xiàn)在TDP和SRP 2種形態(tài)上。

    圖2 上覆水體各形態(tài)磷素濃度 Fig.2 Concentration of different forms of phosphorus in the overlying water

    SRP主要來源于外源輸入和內(nèi)源釋放。水中SRP濃度很小的原因:1)磷化物不易溶于水;2)其易與鈣、鐵及鋁等離子發(fā)生反應(yīng),生成難溶性的沉積物?;钚粤资悄軌蛑苯颖簧锢玫娜芙鈶B(tài)磷,其易被生物吸收同化而進(jìn)入生物循環(huán),難溶性的磷酸鹽和有機(jī)磷必須經(jīng)微生物作用,轉(zhuǎn)化為溶解性的磷酸鹽,才能被生物體直接吸收利用[12]。水體中的磷主要通過基質(zhì)吸附、沉降、植物吸收、微生物轉(zhuǎn)換利用等方式去除,影響因素包括水體中的磷濃度、水力擾動(dòng)、濁度、DO濃度、pH、Eh、溫度、水生植物等[13]。在下游斷面明水區(qū)域TDP濃度的突然增加可能是由于采樣點(diǎn)的流域面積大、流速過慢引起的,由于采樣時(shí)正值秋季,大量浮水植物衰敗腐爛,向水體釋放部分可溶解態(tài)磷,而流速過慢不利于磷素的沿程遷移擴(kuò)散。在上游斷面,植物區(qū)域的各形態(tài)磷濃度明顯高于明水區(qū)域,其原因可能是大型挺水植物群落對上覆水體中的磷素不僅有一定截留作用,同時(shí)對流速也有一定的減緩作用,而明水區(qū)域流速較高,所以上游斷面植物群落各形態(tài)磷濃度較高。

    Zhang等[14]在2012年10月對七星湖濕地的采樣分析中發(fā)現(xiàn),水體TP濃度沿水流方向從入水口的1.16 mg/L逐漸降至出水口的0.45 mg/L,SRP為水體TP的主要成分。而在本次采樣分析中,各采樣點(diǎn)的TP濃度均低于0.10 mg/L,PP為水體TP的主要成分,占TP比例為38%~95%。在這3年人為影響和環(huán)境影響的共同作用下,七星湖濕地采樣區(qū)域的TP濃度得到了明顯的控制,其主要成分由SRP變?yōu)镻P,其中很重要的原因是七星湖濕地在無系統(tǒng)收割的情況下水生植物繁多,增加了對水體TP的吸收和截留,且水生植物在水體中截留的磷形態(tài)主要是SRP。隨著季節(jié)的變化,大多數(shù)植物在秋冬季逐漸衰敗腐爛而繼續(xù)留在濕地生態(tài)系統(tǒng)中,但植物的殘?bào)w會(huì)向水體排放顆粒態(tài)磷,因此濕地水體的主要成分由SRP變成PP。

    2.2 植物體內(nèi)磷素的分布特征

    植物在一定條件下可以吸收極少部分水中的溶解磷,而植物體內(nèi)的磷大多是通過根系微生物的分解作用,吸收沉積物中的部分形態(tài)磷。由圖3(a)可以看出,蘆葦各器官磷素所占比例沿程變化不大,始終是根和葉所占比例較大,而沿程各株蘆葦莖部的磷素所占比例為16.0%~24.3%,且下游斷面蘆葦莖部磷素所占比例明顯低于上游和中游斷面,而葉部比例較高。下游水域開闊,植物生長環(huán)境較好,因此蘆葦群落葉落期比上游和中游稍緩,這可能是影響下游段蘆葦根莖葉磷含量分布的因素之一。由圖3(b)可以看出,上游斷面香蒲根、莖、葉磷素所占比相差不大,但在中游和下游斷面根部所占比例明顯大于莖部和葉部,且中游和下游各器官磷素所占比例變化不大。在同一斷面上的蘆葦和香蒲各器官磷素所占比例也有一定的區(qū)別,蘆葦莖部磷素所占比例均比香蒲低。

    圖3 植物根、莖、葉中磷素所占比例 Fig.3 Concentration of phosphorus in all organs of plan

    圖4 間隙水各形態(tài)磷素的濃度Fig.4 Concentration of different forms of phosphorus in the pore water of sediment

    蘆葦和香蒲均屬于大型挺水植物,根系發(fā)達(dá),但蘆葦?shù)那o部中空,可能是導(dǎo)致其磷素所占比例較低的原因之一。此次采樣是在秋季,濕地大型挺水植物正處于由旺盛到衰敗的過渡期。目前植物地上部分的磷含量高于地下部分,因?yàn)樵谥参锼ネ銎跔I養(yǎng)鹽會(huì)從地上部分轉(zhuǎn)移到地下部分,因此在植物衰亡期到來之前收割比在生長期收割能增加因收割而帶出系統(tǒng)的磷量,在夏季初期收割植物可能會(huì)導(dǎo)致以后幾年里該植物的長期活力[4]。

    2.3 間隙水-沉積物的磷素分布特征

    2.3.1 間隙水磷素分布特征

    在內(nèi)源性營養(yǎng)鹽釋放吸附過程中,沉積物-間隙水-上覆水三者之間的關(guān)系異常密切,間隙水是水土界面交換的重要介質(zhì)[15]。由圖4可以看出,在上、中、下游3個(gè)斷面間隙水中明水區(qū)域各形態(tài)磷素濃度均比植物區(qū)域高,可能是大型挺水植物對磷素的截留吸收轉(zhuǎn)化起到了關(guān)鍵作用。而在水流沿途中,中游斷面營養(yǎng)鹽濃度稍微下降,但下降趨勢不大,這可能是中游水域面積小,水力擾動(dòng)較大而加速了沉積物-上覆水體界面的磷素?cái)U(kuò)散,導(dǎo)致間隙水中磷素濃度降低。而在中游斷面植物區(qū)域,香蒲群落各形態(tài)磷濃度均高于蘆葦群落。雖然蘆葦和香蒲都屬于根系發(fā)達(dá)的大型挺水植物,但香蒲的根系比蘆葦旁支更多,因此根系的吸收轉(zhuǎn)化作用更強(qiáng)。

    在與本試驗(yàn)同一季節(jié)開展的室內(nèi)試驗(yàn)中,對比沉水植物對水體磷素的影響,測得黑藻組間隙水SRP濃度明顯高于對照組(無植物組),SRP濃度最大值為0.61 mg/L[16]。七星湖濕地明水區(qū)域沉積物間隙水SRP濃度明顯高于植物區(qū)域,SRP濃度最大值為1.51 mg/L。由此可見,室內(nèi)試驗(yàn)和室外試驗(yàn)得出的結(jié)果具有明顯的差異性,甚至可能會(huì)是完全相反的結(jié)論。在同樣是室外試驗(yàn)的研究中發(fā)現(xiàn),白塔河沉積物間隙水中TP濃度在農(nóng)村帶河段最低,城鎮(zhèn)帶和城市帶河段較高[9]。農(nóng)村帶河段植被較豐富,而城鎮(zhèn)帶和城市帶河段水域較為開闊,這與本研究得出的結(jié)論相近。

    2.3.2 沉積物中磷素的組成

    七星湖濕地各采樣點(diǎn)的沉積物TP濃度為511.976~1 292.098 mg/kg,平均值為954.901 mg/kg。根據(jù)崔松波[17]在2012年10月對七星湖沉積物各形態(tài)磷的研究可知,當(dāng)時(shí)TP濃度平均值為671.500 mg/kg。七星湖近自然濕地經(jīng)過3年的磷素沉降,沉積物TP濃度增加了42.2%。目前國內(nèi)對濕地沉積物的研究較多,不同濕地沉積物的TP濃度差異也非常大。邵學(xué)新[20]以我國杭州灣南岸典型潮灘濕地為研究區(qū),得到其0~5 cm沉積物的TP濃度為560~680 mg/kg;謝凱[21]在對淮南潘謝采煤沉陷區(qū)的研究中指出,其表層沉積物的TP濃度為311.2~388.4 mg/kg;張巖[22]的研究表明,人工濕地基質(zhì)中TP濃度為139.73~242.03 mg/kg。與已有的研究相比,七星湖沉積物TP濃度偏高,有更大的內(nèi)源污染的危機(jī),對水體造成二次污染的風(fēng)險(xiǎn)較大,加劇了水體富營養(yǎng)化。而七星湖濕地沉積物TP濃度偏高的一個(gè)原因可能是在濕地建立初期曾人為地向濕地中增加燒結(jié)法赤泥以達(dá)到水質(zhì)凈化的目的。赤泥是制鋁工業(yè)提取氧化鋁時(shí)排出的污染性廢渣,燒結(jié)法赤泥具有高鈣、高硅而低鐵的特點(diǎn),高濃度鈣離子容易與磷素結(jié)合生成鈣結(jié)合態(tài)磷(Ca-P),沉降在濕地沉積物中,從而增加了七星湖濕地沉積物TP濃度。

    采用SMT分級方法分析得到濕地土壤磷組成,結(jié)果如圖5所示。由圖5可見,IP占沉積物TP的69%以上,其中,HCl-P占36%~68%,NaOH-P占11%~38%,OP占11%~31%??傮w而言,七星湖濕地沉積物中的磷以IP為主,而IP又多以HCl-P為主。上游和中游的明水區(qū)域沉積物IP濃度明顯低于植物區(qū)域,其中一個(gè)重要原因就是大型挺水植物的截留作用。生活在淺水區(qū)的挺水植物,一方面具有較大的生物量,相對較長的生長周期;另一方面與濕地水體有較多的直接接觸,因此大型挺水植物是濕地磷素遷移轉(zhuǎn)化的重要環(huán)節(jié)。而下游的植物群落OP濃度高于明水區(qū)域,OP主要來源于外源輸入,在試驗(yàn)區(qū)域的下游附近有家禽養(yǎng)殖,家禽平時(shí)活動(dòng)區(qū)域主要在淺水區(qū)的植物群落中,這可能是導(dǎo)致下游植物群落OP濃度偏高的主要原因。同時(shí),在植物衰亡期,不同的植物葉落速度不同,腐敗的葉富含大量的磷,葉落后經(jīng)過微生物的分解會(huì)增加沉積物IP濃度,這也是導(dǎo)致各采樣點(diǎn)不同形態(tài)磷濃度有明顯差距的原因之一。

    2.4 沉積物-上覆水體擴(kuò)散通量估算

    間隙水中的營養(yǎng)鹽向沉積物表層擴(kuò)散進(jìn)而向上覆水混合擴(kuò)散的過程,主要是由濃度差支配的。由于間隙水與上覆水中營養(yǎng)鹽存在濃度差異,必然存在由高濃度向低濃度進(jìn)行的分子擴(kuò)散作用,因此,研究沉積物-上覆水體界面間營養(yǎng)鹽的擴(kuò)散通量具有重要的環(huán)境意義[9]。根據(jù)二者之間的濃度梯度及其物理化學(xué)性質(zhì)可以估算沉積物-上覆水體界面處營養(yǎng)鹽的擴(kuò)散通量。根據(jù)Fick第一擴(kuò)散定律及文獻(xiàn)[21],其改進(jìn)公式為:

    F=φ·Ds·αc/αx

    (1)

    式中:F為分子擴(kuò)散通量,mg/(m2·d);φ為表層沉積物的孔隙度,%;Ds為表層沉積物中物質(zhì)的擴(kuò)散系數(shù),cm2/s;αc/αx為界面濃度梯度,mg/(L·cm)。

    當(dāng)φ≤70%時(shí),Ds=φ·D0;當(dāng)φ≥70%時(shí),Ds=φ2·D0[22]。D0為無限稀釋溶液中溶質(zhì)的擴(kuò)散系數(shù),6.12×10-6cm2/s[23-25]。根據(jù)文獻(xiàn)[26]估計(jì),φ平均值為65%。αc/αx用表層沉積物(0~5 cm)間隙水與上覆水營養(yǎng)鹽濃度差估算求得。目前,磷從沉積物向水體釋放可分為2個(gè)過程:微粒態(tài)磷的活化和沉積物中可溶態(tài)磷的擴(kuò)散作用,其中微粒態(tài)磷的活化主要是跟沉積物中鐵、鋁等結(jié)合,而這種擴(kuò)散作用主要發(fā)生在沉積物表層的間隙水與上覆水之間直接的交換作用。本文主要研究沉積物-上覆水體界面的TDP和SRP擴(kuò)散作用。

    圖5 各采樣點(diǎn)SMT分級所得不同結(jié)合態(tài)磷所占比例Fig.5 Proportion of various phosphorus fractions in soils by SMT protocol

    圖6 間隙水-上覆水體界面TDP和SRP的擴(kuò)散通量Fig.6 Diffusion fluxes of SRP and TDP across the sediment-water interface

    由圖6可以看出,在間隙水-上覆水界面TDP和SRP濃度的擴(kuò)散通量變化趨勢基本相同,且均為正值,說明3個(gè)區(qū)域各采樣點(diǎn)磷素都是由沉積物向上覆水體擴(kuò)散。在上、中、下游斷面,明水區(qū)域擴(kuò)散通量都高于植物區(qū)域,其原因:1)明水區(qū)域上覆水體pH均高于植物區(qū)域,pH在一定范圍內(nèi)越大磷釋放量也越大;2)植物群落會(huì)通過根系以及植物腐爛形成的有機(jī)質(zhì)吸收部分磷素,導(dǎo)致擴(kuò)散到上覆水體中的磷素減少。沉積物中磷的釋放與沉積物氧化還原條件密切相關(guān),好氧條件會(huì)促進(jìn)沉積物對磷的吸附,厭氧條件會(huì)促進(jìn)沉積物對磷的釋放。在中游斷面,香蒲群落TDP和SRP的擴(kuò)散通量陡然升高,由表2可以看出,該采樣點(diǎn)的DO濃度為6.06 mg/L,是所有采樣點(diǎn)中最低值,同時(shí)水力擾動(dòng)較大也會(huì)促進(jìn)磷素在沉積物-上覆水體中的擴(kuò)散。

    此外,沉積物-上覆水體水界面間營養(yǎng)鹽的遷移并不是單純由濃度梯度擴(kuò)散來控制的,同時(shí)還受到來水、水動(dòng)力等因素(如生物的擾動(dòng)作用、營養(yǎng)鹽在水體中的水平遷移擴(kuò)散、風(fēng)浪作用造成的紊流擴(kuò)散等)的影響。由表3可以看出:各形態(tài)磷在間隙水與上覆水體之間具有較好的顯著相關(guān)性;上覆水體與間隙水中的TDP相關(guān)系數(shù)為0.693(P<0.05),SRP相關(guān)系數(shù)為0.692(P<0.05);上覆水體中的TDP與間隙水中的SRP相關(guān)系數(shù)為0.677(P<0.05);上覆水體中的SRP與間隙水中的TDP相關(guān)系數(shù)為0.726(P<0.05)。由此也證明了在間隙水和上覆水體中擴(kuò)散的主要磷形態(tài)是TDP和SRP。

    表3 上覆水體和間隙水各形態(tài)磷之間的相關(guān)性

    注:*在0.05水平上顯著相關(guān),**在0.01水平上顯著相關(guān)(雙尾檢驗(yàn));樣本數(shù)為9。

    3 結(jié)論

    (1) 明水區(qū)域水深以及上覆水體的pH和DO濃度均高于植物區(qū)域,其可能原因:1)明水區(qū)域浮水植物多,光合作用強(qiáng);2)大型挺水植物枝葉茂密,影響了植物光合作用。

    (2)明水區(qū)域上覆水體中各形態(tài)磷沿程波動(dòng)較大,蘆葦群落沿程降低,而香蒲群落各形態(tài)磷濃度在中游變化較大,可能是中游水域面積較小,人為擾動(dòng)較大。明水區(qū)域與植物區(qū)域相比,TDP和SRP在沿程去除效果尤為明顯,表明植物吸收轉(zhuǎn)化的主要磷形態(tài)是TDP和SRP。

    (3)上、中、下游斷面大型挺水植物各器官磷素所占比例有明顯差異,且地上部分高于地下部分。這與采樣季節(jié)有直接關(guān)系,秋季是植物衰亡期,營養(yǎng)鹽會(huì)從地上部分轉(zhuǎn)移到地下部分,因此在植物衰亡期到來之前收割比在生長期收割能增加因收割而帶出系統(tǒng)的磷量。

    (4)上、中、下游3個(gè)斷面間隙水中明水區(qū)域各形態(tài)磷素濃度均比植物區(qū)域高,說明植物對沉積物中磷素的截留吸收轉(zhuǎn)化起到了重要作用。中游水域面積小,水力擾動(dòng)強(qiáng),導(dǎo)致間隙水中各形態(tài)磷濃度降低。七星湖近自然濕地經(jīng)過3年的磷素沉降,沉積物TP濃度增加了42.2%,七星湖沉積物TP濃度偏高,有更大的內(nèi)源污染危機(jī),對水體造成二次污染的風(fēng)險(xiǎn)較大,加劇了水體富營養(yǎng)化。

    (5)七星湖沉積物-上覆水體間TP、TDP、SRP、PP的擴(kuò)散通量均為正通量,表明磷素從沉積物向上覆水體擴(kuò)散,且明水區(qū)域擴(kuò)散通量都高于植物區(qū)域。由表層沉積物間隙水和上覆水體各形態(tài)磷的相關(guān)性表明,間隙水中的磷主要以TDP和SRP的形態(tài)向上覆水體擴(kuò)散。

    [1] 黃錦樓,陳琴,許連煌.人工濕地在應(yīng)用中存在的問題及解決措施[J].環(huán)境科學(xué),2013,34(1):401-408. HUANG J L,CHEN Q,XU L H.Problems and countermeasures in the application of constructed wetlands [J].Environmental Science,2013,34(1):401-408.

    [2] 盧少勇,金相燦,余剛.人工濕地的磷去除機(jī)理[J].生態(tài)環(huán)境,2006,15(2):391-396. LU S Y,JIN X C,YU G.Phosphorus removal mechanism of constructed wetland[J].Ecology and Environment,2006,15(2):391-396.

    [3] JIN X C,WANG S R,PANG Y.The adsorption of phosphate on different trophic lake sediments[J].Colloids & Surfaces A:Physicochemical & Engineering Aspects,2005,254:241-248.

    [4] MEULEMAN A F M,van LOGTESTIJN R,RIJS G B J,et al.Water and mass budgets of a vertical-flow constructed wetland used for wastewater treatment[J].Ecological Engineering,2003,20(1):31-44.

    [5] 歐維新,楊桂山,高建華.鹽城潮灘濕地對N、P營養(yǎng)物質(zhì)的截留效應(yīng)研究[J].濕地科學(xué),2006,4(3):179-186. OU W X,YANG G S,GAO J H.Retention effect of wetland for nitrogen and phosphorus nutrients in the coastal zone of the Yancheng [J].Wetland Science,2006,4(3):179-186.

    [6] PANT H K,REDDY K R,LEMON E.Phosphorus retention capacity of root bed media of sub-surface flow constructed wetlands[J].Ecological Engineering,2001,17:345-355.

    [7] 宋金明.中國近海沉積物-海水界面化學(xué)[M].北京:海洋出版社,1997.

    [8] 張彥,張遠(yuǎn),于濤,等.太湖沉積物及孔隙水中氮的時(shí)空分布特征[J].環(huán)境科學(xué)研究,2010,23(11):1333-1341. ZHANG Y,ZHANG Y,YU T,et al.Spatial and temporal distribution of nitrogen species in sediment and interstitial waters of Taihu Lake[J].Research of Environmental Sciences,2010,23(11):1333-1341.

    [9] 楊楠,于會(huì)彬,宋永會(huì),等.白塔堡河上覆水與沉積物間隙水N、P分布特征[J].環(huán)境科學(xué)研究,2013,26(7):728-735. YANG N,YU H B,SONG Y H,et al.Distribution of nitrogen and phosphorus in overlying water and pore water of sediment in Baitabuhe River[J].Reserach of Environmental Sciences,2013,26(7):728-735.

    [10] 張軍毅,黃君,嚴(yán)飛,等.梅梁湖水體溶解氧特征及其與pH的關(guān)系分析[J].復(fù)旦學(xué)報(bào)(自然科學(xué)版),2009,48(5):623-627. ZHANG J Y,HAUGN J,YAN F,et al.Preliminary study on characters of dissolved oxygen and the relationship with pH in Meiliang Lake[J].Journal of Fudan University(Natural Science),2009,48(5) 623-627.

    [11] 金送笛,李永函,王永利.幾種生態(tài)因子對菹草光合作用的影響[J].水生生物學(xué)報(bào),1991,15(4):295-301. JIN S D,LI Y H,WANG Y L.Effect of ecological factors on photosynthesis ofPotamogetoncrispus[J].Acta Hydrobiologica Sinica,1991,15(4):295-302.

    [12] 傅玲.水生植物群落對水體中氮磷凈化效果研究[D].南京:南京師范大學(xué),2014.

    [13] 吳曉磊.人工濕地廢水處理機(jī)理[J].環(huán)境科學(xué),1995,16(3):83-86. WU X L.Mechanism of wastewater treatment in constructed wetland[J].Environmental Science,1995,16(3):83-86.

    [14] ZHANG H Y,TIAN Y L,CUI S B,et al.Influence of macrophytes on phosphorus fractionation in surface sediments in a constructed wetland:insight from sediment compositions[J].Ecological Engineering,2016,97:400-409.

    [15] 王志齊,李寶,梁仁君,等.南四湖沉積物磷形態(tài)及其與間隙水磷的相關(guān)性分析[J].環(huán)境科學(xué)學(xué)報(bào),2013,33(1):139-146. WANG Z Q,LI B,LAING R J,et al.Correlation analyses of sediment phosphorus forms and soluble phosphorus concentrations in the interstitial water of Nansi Lake,China[J].Acta Scientiae Circumstantiae,2013,33(1) 139-146.

    [16] 王立志.沉水植物對湖泊沉積物-上覆水中磷遷移轉(zhuǎn)化的影響[D].南京:南京師范大學(xué),2011.

    [17] 崔松波.濕地表層底泥磷的分布及影響因素研究[D].北京:華北電力大學(xué),2014. CUI S B.Research on distribution and influencing factor of phosphorus in surficial sediments of wetland[D].Beijing:North China Electric Power University,2014.

    [18] 邵學(xué)新.潮灘濕地植物對磷素遷移轉(zhuǎn)化及截留的影響機(jī)制研究[D].杭州:浙江大學(xué),2014. SHAO X X.Effects of macrophytes on phosphorus transfer and speciation in intertidal march[D].Hangzhou:Zhejiang University,2014.

    [19] 謝凱.淮南潘謝采煤沉陷積水區(qū)沉積物磷遷移轉(zhuǎn)化過程與機(jī)制研究[D].北京:中國礦業(yè)大學(xué),2014. XIE K.Processes and mechanisms of phosphorus transform and transport in the sediments of aquatic zones around the Huainan Panxie coal mine subsidence areas[D].Beijing:China University of Mining Technology,2014.

    [20] 張巖.自由表面流-水平潛流復(fù)合人工濕地磷去除動(dòng)力學(xué)分析[D].北京:中國林業(yè)科學(xué)研究院,2014. ZHANG Y.Phosphorus dynamics in integrated constructed wetland with surface and subsurface flow[D].Beijing:Chinese Academy of Forestry,2014.

    [21] MCCOMB A J,QIU S,LUKATELICH R J,et al.Spatial and temporal heterogeneity of sediment phosphorus in the Pee Hlarvey estuarine system[J].Estuarine,Coastal and Shelf Science,1998,47(5):561-577.

    [22] UlLMAN W J,SANDSTRON M W.Dissolved nutrient fluxes from the nearshore sediments of Bowling Green Bay,Central Great Barrier Reef Lagoon(Australia)[J].Estuarine,Coastal and Shelf Science,1987,24(3):289-303.

    [23] FALCAO M,VALE C.Sediment-water exchanges of ammonium and phosphate in intertidal and subtidal areas of a Mesotidal Coastal Lagoon(Ria Formosa)[J].Hydrobiologia,1998,373/374:193-201.

    [24] KROM M D,BERNER R A.The diffusion coefficients of sulfate,ammonium,and phosphate ions in anoxic marine sediments[J].Limnology & Oceanography,1980,25(2):327-337.

    [25] WARNKEN K W,GILL G A,SANTSCHI P H,et al.Benthic exchange of nutrients in Galveston Bay,Texas[J].Estuaries and Coasts,2000,23(5):647-661.

    [26] 盧博,張福生,黃韶健,等.大亞灣表層沉積物性質(zhì)及其對海水養(yǎng)殖的影響[J].臺灣海峽,2002,21(4):489-495. LU B,ZHANG F S,HUANG S J,et al.Surface sediment properties and their influence on sea water culture in Daya Bay[J].Journal of Oceanography in Taiwan Strait,2002,21(4):489-495. ?

    Distribution of phosphorus and estimation of diffusion flux in the Qixinghu wetland

    ZHONG Xin, TIAN Yonglan, ZHONG Meifang, ZHANG Huayong

    Research Center for Engineering Ecology and Nonlinear Science, North China Electric Power University, Beijing 102206, China

    To study the distribution and migration of phosphorus in the Qixinghu wetland in Tengzhou City of Shangdong Province, the samples of overlying water and sediment of Qixinghu wetland were collected, the distribution characteristics of total phosphorus(TP), total dissolved phosphorus(TDP), solubility reactive phosphorus(SRP) and particulate phosphorus(PP) analyzed, and the relationship of P in overlying water and pore water of sediment studied. The results showed that, TP concentrations in overlying water of open water area were stable along with the flow direction (0.05-0.06 mg/L), declined in overlying water of reed community area (0.04-0.09 mg/L) and fluctuated in overlying water of cattail community area (0.02-0.09 mg/L). All measured phosphorus fractionations (TP, TDP, SRP and PP) in pore water of sediment of open water area were higher than those of plant communities, indicating that plants played an important role in the absorption and transformation of phosphorus in sediments. Phosphorus contents in aerial organs of macrophytes were significantly higher than the underground part. After settlement of phosphorus in three years, the content of TP in sediments increased by 42.2% in the wetland. The diffusion flux across sediment-water interface of TDP and SRP were positive, suggesting that phosphorus diffuses from sediment to overlying water. Moreover, diffusion fluxes of open water were higher than those of the plant communities. As indicated by the various forms of phosphorus correlation of pore water and overlying water, the phosphorus of pore water mainly diffuses into the overlying water by TDP and SRP.

    pore water of sediment; wetland; macrophytes; phosphorus fractionations; diffusion flux

    2016-09-09

    國家水體污染控制與治理科技重大專項(xiàng)(2009ZX07210-009,2015ZX07203-011,2015ZX07204-007);國家自然科學(xué)基金項(xiàng)目(39560023)

    鐘馨(1992—),女,碩士研究生,研究方向?yàn)樗廴究刂婆c生態(tài)修復(fù),zhongxin0531@126.com

    *責(zé)任作者:張化永(1963—),男,教授,博士,主要研究方向?yàn)榛謴?fù)生態(tài)學(xué),rceens@ncepu.edu.cn

    文章編號:1674-991X(2017)02-0152-09 doi:10.3969/j.issn.1674-991X.2017.02.023

    鐘馨,田永蘭,鐘美芳,等.七星湖濕地磷素分布特征及擴(kuò)散通量估算[J].環(huán)境工程技術(shù)學(xué)報(bào),2017,7(2):152-160.

    ZHONG X,TIAN Y L,ZHONG M F,et al.Distribution of phosphorus and estimation of diffusion flux in the Qixinghu wetland[J].Journal of Environmental Engineering Technology,2017,7(2):152-160.

    猜你喜歡
    明水磷素七星
    磷素添加對土壤水分一維垂直入滲特性的影響
    明水(甘新界)至哈密段公路高速預(yù)檢稱重系統(tǒng)設(shè)計(jì)
    定心冰糖
    七星湖
    草原歌聲(2020年2期)2020-09-25 08:38:00
    七星瓢蟲
    大灰狼(2018年5期)2018-06-20 14:49:32
    七星瓢蟲
    七星瓢蟲
    墻頭記
    作物高效利用土壤磷素的研究進(jìn)展
    土壤與作物(2013年4期)2013-03-11 18:47:35
    磷素營養(yǎng)對大豆磷素吸收及產(chǎn)量的影響
    亚洲一卡2卡3卡4卡5卡精品中文| 99国产精品免费福利视频| 成人漫画全彩无遮挡| 成人三级做爰电影| 一级毛片 在线播放| 极品人妻少妇av视频| 亚洲av日韩在线播放| 两个人免费观看高清视频| 在线观看www视频免费| 国产一区二区 视频在线| 精品午夜福利在线看| 熟妇人妻不卡中文字幕| 国产av国产精品国产| 精品一区在线观看国产| 中文字幕最新亚洲高清| 国产激情久久老熟女| 菩萨蛮人人尽说江南好唐韦庄| 国产免费现黄频在线看| 丝袜在线中文字幕| 亚洲欧美激情在线| 最近2019中文字幕mv第一页| 亚洲成人av在线免费| 久久精品aⅴ一区二区三区四区| 日本91视频免费播放| 韩国高清视频一区二区三区| 精品国产国语对白av| av卡一久久| 亚洲情色 制服丝袜| 一级毛片电影观看| 中文字幕色久视频| 日韩不卡一区二区三区视频在线| 国产黄色免费在线视频| 天堂中文最新版在线下载| 国产av精品麻豆| 看十八女毛片水多多多| 成人亚洲精品一区在线观看| 啦啦啦 在线观看视频| 少妇被粗大的猛进出69影院| 免费少妇av软件| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品国产区一区二| 老司机深夜福利视频在线观看 | 黑人欧美特级aaaaaa片| 18禁裸乳无遮挡动漫免费视频| 国产一区亚洲一区在线观看| 欧美另类一区| 人人澡人人妻人| 午夜日韩欧美国产| 咕卡用的链子| 一区二区三区乱码不卡18| 无限看片的www在线观看| 夜夜骑夜夜射夜夜干| 老司机深夜福利视频在线观看 | 午夜福利,免费看| 1024香蕉在线观看| 国产在线一区二区三区精| av在线老鸭窝| 午夜免费观看性视频| 久久毛片免费看一区二区三区| h视频一区二区三区| 亚洲第一av免费看| 80岁老熟妇乱子伦牲交| 亚洲精品美女久久久久99蜜臀 | 九草在线视频观看| 国产精品.久久久| 看十八女毛片水多多多| 国产av精品麻豆| 免费女性裸体啪啪无遮挡网站| 精品国产乱码久久久久久男人| 亚洲精品aⅴ在线观看| 18在线观看网站| 99久久精品国产亚洲精品| 男女高潮啪啪啪动态图| 亚洲av电影在线观看一区二区三区| 欧美日韩视频高清一区二区三区二| 91精品伊人久久大香线蕉| 一级毛片我不卡| 色婷婷久久久亚洲欧美| 国产 精品1| 啦啦啦啦在线视频资源| 亚洲av成人不卡在线观看播放网 | 久久影院123| 久久影院123| 国语对白做爰xxxⅹ性视频网站| 国产精品一区二区精品视频观看| 一二三四在线观看免费中文在| 另类精品久久| 超碰97精品在线观看| 亚洲成人免费av在线播放| 男女无遮挡免费网站观看| 在线观看免费视频网站a站| 99re6热这里在线精品视频| 90打野战视频偷拍视频| 视频在线观看一区二区三区| 免费在线观看黄色视频的| 狂野欧美激情性xxxx| 成人影院久久| 午夜福利视频精品| 中文乱码字字幕精品一区二区三区| 亚洲三区欧美一区| 97在线人人人人妻| 赤兔流量卡办理| 90打野战视频偷拍视频| 中文字幕人妻丝袜制服| 亚洲av欧美aⅴ国产| av福利片在线| 最新在线观看一区二区三区 | 午夜福利网站1000一区二区三区| 只有这里有精品99| 最黄视频免费看| 自拍欧美九色日韩亚洲蝌蚪91| 欧美另类一区| 日本av手机在线免费观看| 日本av手机在线免费观看| 成人免费观看视频高清| 日韩一区二区三区影片| 性色av一级| 亚洲欧美日韩另类电影网站| 国产亚洲av高清不卡| 亚洲精品中文字幕在线视频| 乱人伦中国视频| 在线观看免费视频网站a站| 国产精品国产av在线观看| 午夜福利一区二区在线看| 亚洲国产毛片av蜜桃av| 午夜免费观看性视频| 国产精品久久久人人做人人爽| av视频免费观看在线观看| 亚洲成av片中文字幕在线观看| 90打野战视频偷拍视频| 国产精品久久久久久人妻精品电影 | 国产男女内射视频| 国产男女内射视频| 亚洲色图综合在线观看| 色婷婷久久久亚洲欧美| 精品少妇一区二区三区视频日本电影 | 在线观看一区二区三区激情| 欧美av亚洲av综合av国产av | 99久久精品国产亚洲精品| 成人国产av品久久久| 国产国语露脸激情在线看| 最新在线观看一区二区三区 | 色精品久久人妻99蜜桃| 日本一区二区免费在线视频| 在线观看三级黄色| 看免费av毛片| 91精品国产国语对白视频| 国产99久久九九免费精品| 男女高潮啪啪啪动态图| 国产免费视频播放在线视频| 一二三四在线观看免费中文在| 国产爽快片一区二区三区| 一本一本久久a久久精品综合妖精| 叶爱在线成人免费视频播放| 亚洲国产欧美一区二区综合| 看非洲黑人一级黄片| 高清av免费在线| 久久青草综合色| 亚洲欧美日韩另类电影网站| 女人高潮潮喷娇喘18禁视频| a级片在线免费高清观看视频| 精品国产一区二区三区四区第35| 国产又爽黄色视频| 不卡视频在线观看欧美| 满18在线观看网站| 亚洲国产av新网站| 青青草视频在线视频观看| av国产久精品久网站免费入址| 日韩中文字幕欧美一区二区 | 悠悠久久av| 国产xxxxx性猛交| 国产又色又爽无遮挡免| 欧美日本中文国产一区发布| 欧美xxⅹ黑人| 国产国语露脸激情在线看| 欧美精品av麻豆av| 夜夜骑夜夜射夜夜干| 国产深夜福利视频在线观看| 久久久精品免费免费高清| 天天添夜夜摸| 十八禁人妻一区二区| 色综合欧美亚洲国产小说| 久久韩国三级中文字幕| 久久精品久久久久久久性| 成年动漫av网址| 国产一区二区三区综合在线观看| 永久免费av网站大全| 日日撸夜夜添| 国产av国产精品国产| 丰满饥渴人妻一区二区三| 精品亚洲成a人片在线观看| 亚洲,欧美,日韩| 一级片免费观看大全| 最近最新中文字幕免费大全7| xxxhd国产人妻xxx| 狂野欧美激情性bbbbbb| 国产日韩欧美亚洲二区| 欧美激情极品国产一区二区三区| 永久免费av网站大全| 亚洲国产av影院在线观看| 看免费av毛片| 国产精品免费视频内射| 国产亚洲精品第一综合不卡| 亚洲中文av在线| 交换朋友夫妻互换小说| 一区二区av电影网| 国产av国产精品国产| 丰满饥渴人妻一区二区三| 久久久欧美国产精品| 成年美女黄网站色视频大全免费| 高清视频免费观看一区二区| 日韩,欧美,国产一区二区三区| 国产成人午夜福利电影在线观看| 久久鲁丝午夜福利片| 18禁国产床啪视频网站| 最近最新中文字幕大全免费视频 | 久久精品亚洲av国产电影网| 成人国语在线视频| 国产99久久九九免费精品| av不卡在线播放| 亚洲精品一区蜜桃| 午夜影院在线不卡| 在线精品无人区一区二区三| 制服丝袜香蕉在线| 国产麻豆69| 久久久久久久精品精品| 亚洲五月色婷婷综合| 欧美日韩亚洲高清精品| kizo精华| 亚洲第一区二区三区不卡| 久久久久精品久久久久真实原创| 十分钟在线观看高清视频www| 亚洲精品av麻豆狂野| 亚洲精品在线美女| 欧美黄色片欧美黄色片| 日本av免费视频播放| 亚洲欧洲精品一区二区精品久久久 | 欧美中文综合在线视频| 亚洲国产欧美网| 韩国精品一区二区三区| 久久久国产精品麻豆| 日本色播在线视频| 亚洲 欧美一区二区三区| 国产爽快片一区二区三区| 欧美黑人欧美精品刺激| 亚洲激情五月婷婷啪啪| 欧美日韩一级在线毛片| 天堂8中文在线网| 一区二区三区激情视频| 在线观看免费视频网站a站| 国产成人精品久久久久久| 黄色视频不卡| 看免费成人av毛片| 色精品久久人妻99蜜桃| 日韩中文字幕视频在线看片| 亚洲一区中文字幕在线| 免费看不卡的av| 亚洲av中文av极速乱| 免费av中文字幕在线| 亚洲五月色婷婷综合| 欧美变态另类bdsm刘玥| 欧美日韩综合久久久久久| 18禁裸乳无遮挡动漫免费视频| 看免费av毛片| 久久精品久久久久久久性| 欧美老熟妇乱子伦牲交| 亚洲欧美成人精品一区二区| 超碰成人久久| 亚洲四区av| 久久精品亚洲熟妇少妇任你| 最近手机中文字幕大全| 2021少妇久久久久久久久久久| 亚洲欧美清纯卡通| 丝袜美足系列| 亚洲国产欧美网| 国产片特级美女逼逼视频| 97精品久久久久久久久久精品| 青春草国产在线视频| 亚洲成人一二三区av| 九九爱精品视频在线观看| 亚洲国产毛片av蜜桃av| 蜜桃在线观看..| 欧美日韩综合久久久久久| 黑丝袜美女国产一区| 日本一区二区免费在线视频| 91精品伊人久久大香线蕉| 伊人久久国产一区二区| 中文字幕人妻丝袜制服| 国产人伦9x9x在线观看| 午夜精品国产一区二区电影| 18禁国产床啪视频网站| 青草久久国产| 熟女少妇亚洲综合色aaa.| 丝袜脚勾引网站| 日韩电影二区| 国产精品久久久久久人妻精品电影 | 成人午夜精彩视频在线观看| 一级a爱视频在线免费观看| 久久精品国产a三级三级三级| 亚洲成av片中文字幕在线观看| 两性夫妻黄色片| 黄片播放在线免费| 亚洲av综合色区一区| www.av在线官网国产| 日本午夜av视频| 精品少妇内射三级| 97在线人人人人妻| 欧美最新免费一区二区三区| 中国三级夫妇交换| 国产成人免费无遮挡视频| 久久久久久人人人人人| 看免费成人av毛片| 如日韩欧美国产精品一区二区三区| 伦理电影免费视频| 精品一区在线观看国产| 国产成人免费无遮挡视频| 成人18禁高潮啪啪吃奶动态图| 久久青草综合色| 中文字幕av电影在线播放| 人人澡人人妻人| 一本久久精品| 欧美人与性动交α欧美软件| 一级黄片播放器| 美女视频免费永久观看网站| 熟女av电影| 免费人妻精品一区二区三区视频| 亚洲av福利一区| 超碰成人久久| 精品一区二区三卡| 伊人亚洲综合成人网| 久久ye,这里只有精品| 日韩中文字幕欧美一区二区 | 亚洲精品久久成人aⅴ小说| av在线app专区| 五月开心婷婷网| 精品久久蜜臀av无| 欧美 亚洲 国产 日韩一| 日本欧美视频一区| 交换朋友夫妻互换小说| 制服人妻中文乱码| 美女视频免费永久观看网站| 国产一区二区激情短视频 | 国产日韩欧美亚洲二区| 交换朋友夫妻互换小说| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久久精品古装| 亚洲av电影在线进入| 免费久久久久久久精品成人欧美视频| av不卡在线播放| 精品亚洲乱码少妇综合久久| 国产一区二区三区综合在线观看| 国产精品偷伦视频观看了| 国精品久久久久久国模美| 97在线人人人人妻| 久久99一区二区三区| 亚洲国产欧美在线一区| 国产精品一区二区精品视频观看| 中文字幕人妻丝袜一区二区 | 久久久久久久久免费视频了| 国产精品久久久久久久久免| 欧美激情高清一区二区三区 | 男人添女人高潮全过程视频| 亚洲av国产av综合av卡| 久久精品久久精品一区二区三区| www.自偷自拍.com| 免费观看av网站的网址| 亚洲欧洲日产国产| 超碰97精品在线观看| 人人妻人人澡人人爽人人夜夜| 国产亚洲欧美精品永久| 国产精品国产三级国产专区5o| 国产精品二区激情视频| 日本午夜av视频| 一本大道久久a久久精品| 蜜桃在线观看..| 精品一区二区三区四区五区乱码 | 在线观看人妻少妇| 一区福利在线观看| 亚洲,一卡二卡三卡| 国产在视频线精品| 建设人人有责人人尽责人人享有的| 欧美日韩亚洲高清精品| 欧美激情极品国产一区二区三区| 午夜日本视频在线| 视频区图区小说| 久久久久人妻精品一区果冻| 黄色视频在线播放观看不卡| 国产野战对白在线观看| 亚洲精品美女久久久久99蜜臀 | 男人舔女人的私密视频| 日韩精品免费视频一区二区三区| 久久久久久久久久久久大奶| 一区二区av电影网| 一区二区三区四区激情视频| 亚洲精品国产av蜜桃| 欧美日韩亚洲综合一区二区三区_| 国产 一区精品| 91老司机精品| 亚洲三区欧美一区| 久久精品aⅴ一区二区三区四区| 亚洲欧美一区二区三区黑人| 国产激情久久老熟女| 欧美日韩av久久| 涩涩av久久男人的天堂| 黄色视频不卡| 又大又黄又爽视频免费| 日日撸夜夜添| 国产成人免费无遮挡视频| 亚洲成人一二三区av| 亚洲国产日韩一区二区| 日韩精品有码人妻一区| 亚洲精品av麻豆狂野| 黑人巨大精品欧美一区二区蜜桃| 老司机影院毛片| 一本色道久久久久久精品综合| 老司机深夜福利视频在线观看 | 男女下面插进去视频免费观看| 国产精品偷伦视频观看了| 欧美日韩视频高清一区二区三区二| 啦啦啦啦在线视频资源| 久久久久精品人妻al黑| 久久久久精品久久久久真实原创| 99久久综合免费| 无遮挡黄片免费观看| 亚洲av日韩在线播放| 国产亚洲午夜精品一区二区久久| 少妇 在线观看| 一本一本久久a久久精品综合妖精| 91国产中文字幕| 天堂8中文在线网| 男人舔女人的私密视频| 黄频高清免费视频| 国产男人的电影天堂91| 亚洲精品第二区| 91精品三级在线观看| 午夜老司机福利片| 国产精品熟女久久久久浪| 亚洲欧美成人精品一区二区| 久久精品久久久久久噜噜老黄| 日本av手机在线免费观看| 亚洲欧美中文字幕日韩二区| 亚洲国产毛片av蜜桃av| 国产成人精品无人区| 国产人伦9x9x在线观看| 悠悠久久av| 美女脱内裤让男人舔精品视频| 在线观看三级黄色| 丰满饥渴人妻一区二区三| av有码第一页| 少妇的丰满在线观看| 少妇人妻精品综合一区二区| 久久精品国产亚洲av涩爱| 亚洲美女视频黄频| 亚洲四区av| 国产精品久久久人人做人人爽| 老鸭窝网址在线观看| 午夜免费鲁丝| 亚洲少妇的诱惑av| 高清欧美精品videossex| 精品国产超薄肉色丝袜足j| 女性生殖器流出的白浆| 少妇的丰满在线观看| 超碰97精品在线观看| 99国产综合亚洲精品| 国产男女内射视频| 狠狠精品人妻久久久久久综合| 高清黄色对白视频在线免费看| 亚洲,一卡二卡三卡| 日韩一区二区视频免费看| 日本午夜av视频| a 毛片基地| 午夜福利免费观看在线| 国产精品三级大全| 在线亚洲精品国产二区图片欧美| 香蕉丝袜av| 日韩av免费高清视频| 天天躁夜夜躁狠狠躁躁| 久久久久久久久久久久大奶| 丁香六月欧美| 欧美变态另类bdsm刘玥| 成人国产麻豆网| 日韩 亚洲 欧美在线| 女人久久www免费人成看片| 2018国产大陆天天弄谢| 国产乱来视频区| 你懂的网址亚洲精品在线观看| 久久久久久久精品精品| 亚洲精品一区蜜桃| 色精品久久人妻99蜜桃| 精品国产国语对白av| 黄片小视频在线播放| 欧美日韩亚洲国产一区二区在线观看 | 精品国产一区二区三区久久久樱花| 啦啦啦在线观看免费高清www| 如何舔出高潮| 日韩av不卡免费在线播放| 一级毛片我不卡| 欧美黑人精品巨大| 又大又黄又爽视频免费| 亚洲欧美精品综合一区二区三区| 久久精品人人爽人人爽视色| 夫妻午夜视频| 亚洲av成人不卡在线观看播放网 | 天天躁日日躁夜夜躁夜夜| 国产精品国产av在线观看| 精品一区二区免费观看| 侵犯人妻中文字幕一二三四区| 午夜av观看不卡| 成年人午夜在线观看视频| 黄色一级大片看看| 制服人妻中文乱码| 国产精品久久久av美女十八| 久久久久精品久久久久真实原创| 麻豆乱淫一区二区| 国产免费福利视频在线观看| 日韩成人av中文字幕在线观看| 国产精品久久久av美女十八| 少妇被粗大猛烈的视频| 97精品久久久久久久久久精品| 国产av码专区亚洲av| 久久久久久久久免费视频了| 男女无遮挡免费网站观看| 国产成人系列免费观看| 别揉我奶头~嗯~啊~动态视频 | 90打野战视频偷拍视频| 搡老乐熟女国产| av在线老鸭窝| 69精品国产乱码久久久| 免费黄色在线免费观看| 男女之事视频高清在线观看 | 超色免费av| 亚洲四区av| 国产男女超爽视频在线观看| 两个人看的免费小视频| tube8黄色片| 久久久久精品久久久久真实原创| 日本av手机在线免费观看| 亚洲美女搞黄在线观看| 国产精品 国内视频| 日韩欧美精品免费久久| 日韩电影二区| 视频区图区小说| 欧美成人午夜精品| 亚洲国产精品一区二区三区在线| 91精品伊人久久大香线蕉| 欧美黑人欧美精品刺激| 啦啦啦在线观看免费高清www| 涩涩av久久男人的天堂| 色婷婷av一区二区三区视频| 亚洲专区中文字幕在线 | 午夜福利影视在线免费观看| 大香蕉久久网| 欧美日韩亚洲国产一区二区在线观看 | 又黄又粗又硬又大视频| 亚洲成av片中文字幕在线观看| 国产精品 欧美亚洲| 涩涩av久久男人的天堂| 十分钟在线观看高清视频www| av在线老鸭窝| 亚洲一级一片aⅴ在线观看| 亚洲精品一区蜜桃| 最新的欧美精品一区二区| 亚洲一区中文字幕在线| 一级毛片黄色毛片免费观看视频| av国产精品久久久久影院| 久久久精品免费免费高清| 美女主播在线视频| 水蜜桃什么品种好| 51午夜福利影视在线观看| 日韩欧美精品免费久久| 91精品三级在线观看| 大片电影免费在线观看免费| 亚洲精品中文字幕在线视频| 亚洲精品久久久久久婷婷小说| 卡戴珊不雅视频在线播放| 在线观看一区二区三区激情| 一区在线观看完整版| a级毛片黄视频| 色94色欧美一区二区| 91精品国产国语对白视频| 国产成人一区二区在线| 色吧在线观看| 悠悠久久av| a级毛片在线看网站| 狂野欧美激情性bbbbbb| 亚洲精品aⅴ在线观看| 亚洲婷婷狠狠爱综合网| 精品少妇黑人巨大在线播放| 一本色道久久久久久精品综合| 女人被躁到高潮嗷嗷叫费观| 国产精品 国内视频| 欧美国产精品va在线观看不卡| 亚洲成色77777| 美女福利国产在线| 亚洲第一青青草原| 国产视频首页在线观看| 80岁老熟妇乱子伦牲交| 最近最新中文字幕免费大全7| 亚洲国产欧美日韩在线播放| 如日韩欧美国产精品一区二区三区| 日韩 亚洲 欧美在线| 天天躁夜夜躁狠狠久久av| 夫妻午夜视频| 久久久欧美国产精品| 国产日韩一区二区三区精品不卡| 亚洲精华国产精华液的使用体验| 男女国产视频网站| 欧美老熟妇乱子伦牲交| 黄片播放在线免费| 亚洲国产日韩一区二区| 啦啦啦中文免费视频观看日本| 亚洲精品久久成人aⅴ小说| 久久久久视频综合| 亚洲精品第二区| 男女下面插进去视频免费观看|