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

    兩小型淺水湖泊水體葉綠素a濃度的季節(jié)變化及與理化因子和生物因子的關(guān)系*

    2017-11-08 05:18:47溫新利馮若楠薛穎昊席貽龍
    湖泊科學(xué) 2017年6期
    關(guān)鍵詞:鏡湖輪蟲淺水

    溫新利,馮若楠,張 根,薛穎昊,3,席貽龍

    (1:安徽師范大學(xué)生命科學(xué)學(xué)院, 安徽省高校生物環(huán)境與生態(tài)安全省級重點實驗室, 蕪湖 241000)(2:皖江城市帶退化生態(tài)系統(tǒng)的恢復(fù)與重建安徽省協(xié)同創(chuàng)新中心, 蕪湖 241000)(3:農(nóng)業(yè)部農(nóng)業(yè)生態(tài)與資源保護(hù)總站, 北京 100125)

    兩小型淺水湖泊水體葉綠素a濃度的季節(jié)變化及與理化因子和生物因子的關(guān)系*

    溫新利1,2,馮若楠1,2,張 根1,薛穎昊1,3,席貽龍1,2**

    (1:安徽師范大學(xué)生命科學(xué)學(xué)院, 安徽省高校生物環(huán)境與生態(tài)安全省級重點實驗室, 蕪湖 241000)(2:皖江城市帶退化生態(tài)系統(tǒng)的恢復(fù)與重建安徽省協(xié)同創(chuàng)新中心, 蕪湖 241000)(3:農(nóng)業(yè)部農(nóng)業(yè)生態(tài)與資源保護(hù)總站, 北京 100125)

    葉綠素a濃度;季節(jié)變動;非生物因子;生物因子;相對重要性;汀棠湖;鏡湖

    葉綠素a(Chl.a)是浮游植物生物體的重要組成成分之一,其濃度高低是衡量浮游植物現(xiàn)存量的重要指標(biāo)[1],也是評價水體富營養(yǎng)化程度的重要因素[2]. 鑒于此,人們非常關(guān)注水環(huán)境因子(如水溫和營養(yǎng)鹽)對Chl.a濃度變動的影響[3-9].

    由于大型淺水湖泊水域面積大,水體常用于水產(chǎn)養(yǎng)殖且與人們的生產(chǎn)生活關(guān)系密切,其Chl.a濃度的變動受到的關(guān)注度很高[3-5,8-9]. 研究表明,由于具有受風(fēng)浪攪動強烈、與外部水資源交換程度較高等特點,大型淺水湖泊中的Chl.a濃度常表現(xiàn)為較明顯的空間異質(zhì)性. 比如: 浮游植物在太湖的深水中央?yún)^(qū)受光照強度的影響大,在較淺的沿岸帶受氮的影響程度更高[10]. 與此相比,在探討小型淺水湖泊Chl.a濃度的時空變動與水環(huán)境間關(guān)系的相關(guān)研究中,為了闡述水體富營養(yǎng)化產(chǎn)生的原因、后果或者治理途徑,水溫和營養(yǎng)鹽[3,5-7,11]、水體異味物質(zhì)[12]及植物化感物質(zhì)[13]等因子備受重視. 自然界中,由于水體中通常都分布有輪蟲和枝角類等浮游動物[14],它們分別以粒徑為1~20和1~50 μm的浮游藻類、微型浮游生物和細(xì)菌等作為食物[15-17],營養(yǎng)生態(tài)位間存在著重疊,且它們在富營養(yǎng)化湖泊中往往占有優(yōu)勢[18],因而水體Chl.a濃度可能會受到輪蟲和枝角類牧食藻類的影響(即下行效應(yīng)),這個推測在南京玄武湖中的原位短期(15 d)實驗中已經(jīng)獲得了證實[19]. 但是,隨著季節(jié)的演替,在更長的時間內(nèi),特別是在藻類現(xiàn)存量較低的季節(jié),浮游動物對藻類的牧食會不會對淺水湖泊Chl.a濃度產(chǎn)生負(fù)面影響?至今,關(guān)于這個問題的研究報道并不多見.

    城市景觀湖泊由于位于城市內(nèi)部,多為小型封閉性淺水水體,易受稠密的人口等外界因素的干擾,因此城市淺水湖泊的富營養(yǎng)化問題也受到了人們的關(guān)注[20]. 為探明城市內(nèi)小型淺水湖泊中Chl.a濃度的季節(jié)變動及其影響因素,本文以蕪湖市內(nèi)的鏡湖和汀棠湖為研究水體,重點分析了非生物因子(如水溫、營養(yǎng)鹽)和生物因子(浮游動物)對2個湖泊Chl.a濃度季節(jié)動態(tài)的影響,為篩選和識別影響淺水湖泊水體Chl.a濃度變動的主要環(huán)境因子提供科學(xué)資料,同時也為Chl.a濃度的預(yù)測及富營養(yǎng)湖泊的生態(tài)治理提供一定的理論依據(jù).

    1 材料與方法

    1.1 研究水體

    鏡湖(31.334°N,118.380°E)和汀棠湖(31.367°N,118.392°E)分別位于蕪湖市(31.20°N,118.21°E)的中心和近郊區(qū),為城市公園景觀水域,其平均水深分別約為1.3和1.5~2.0 m,水域面積約為7.9和13.47 hm2[21],皆為小型淺水湖泊. 2個湖泊水體的營養(yǎng)水平接近,均處于富營養(yǎng)化初期階段[21-22]. 鏡湖水體中常年缺乏浮游甲殼動物(如枝角類和橈足類)[23].

    1.2 水樣采集

    水樣采集時,分別在鏡湖東、西2個湖區(qū),汀棠湖南、北水體的中心各設(shè)置1個采樣點(圖1).

    圖1 鏡湖和汀棠湖采樣站的設(shè)置Fig.1 Sampling sites in Lake Jinghu and Lake Tingtang

    鏡湖水樣的采集在2008年5月-2009年4月期間完成,除2008年11月份采集2次外,其余時間段均采取每4天或7天1次的樣品采集頻率. 汀棠湖水樣采集在2009年2月-2010年1月期間完成,除2009年10月份因天氣原因只進(jìn)行1次采集外,其他時間段均為每半月采集1次樣品. 用于水體總氮(TN)和總磷(TP)濃度測定的水樣每月采集1次,均在每月上旬進(jìn)行. 所有的樣品均在采樣日上午8:00-10:00之間進(jìn)行.

    水樣采集時,同步用水銀溫度計測定水溫,賽奇氏圓盤(Secchi disc)測定水體透明度(SD);分別用哈納HI-9143溶氧儀(意大利)和哈納HI-8424酸度計(意大利)測定溶解氧(DO)和水體pH值. TN和TP濃度的測定委托蕪湖市水文水資源局完成,其方法分別參照國標(biāo)GB 11894-1989和GB 11893-1989進(jìn)行. 為反映水體中粒徑大小不同的藻類生物量的多少,在Chl.a濃度測定時分為兩個程序進(jìn)行: 一是將來源于各水層中的等量水樣混合后,經(jīng)孔徑為25 μm 的篩絹過濾,測定過濾后水樣中的Chl.a濃度(Chl.a<25 μm),以間接反映水體中小型藻類的生物量; 二是將來源于各水層的混合水樣未經(jīng)過濾而直接帶回實驗室用于Chl.a的測定,以間接反映水體中藻類總體生物量. 兩類型Chl.a濃度的測定均采用丙酮萃取—分光光度法[24].

    包括輪蟲在內(nèi)的浮游動物的定性標(biāo)本采用25#浮游生物網(wǎng)大范圍拖撈獲得; 定量標(biāo)本采用2.5 L采水器,于3個不同深度各采集5 L湖水,混和后經(jīng)孔徑為25 μm的篩絹過濾收集,即刻用5%的福爾馬林溶液固定,實驗室內(nèi)靜置沉淀48 h后,收集沉淀物并濃縮至30 ml. 浮游動物計數(shù)采用1 ml計數(shù)框,在OLYMPUS顯微鏡下全片計數(shù),取3片計數(shù)的平均值作為其密度(橈足類分幼體和成體分別計數(shù)). 輪蟲的種類鑒定按Koste分類系統(tǒng)[25]; 枝角類的鑒定參照《中國動物志——淡水枝角類》[26],枝角類的分類地位鑒定至屬.

    1.3 數(shù)據(jù)分析

    季節(jié)的劃分參照氣象學(xué)上的氣溫法進(jìn)行. 具體是先按照水溫與氣溫間的關(guān)系[3],將測得的水溫?fù)Q算成氣溫,將連續(xù)5天的均溫高于22℃的時期視為夏季,而低于10℃的視為冬季,二者之間為春季和秋季. 根據(jù)輪蟲對不同粒徑藻類和其他食物資源的攝食情況[27-28],將2個湖泊中的輪蟲劃分成三種類型,其一為微型濾食性輪蟲,包括臂尾輪蟲(Brachionus)、龜甲輪蟲(Keratella)、龜紋輪蟲(Anuraeopsis)、葉輪蟲(Notholca)和三肢輪蟲(Filinia);其二為大型藻類掠食性輪蟲,包括異尾輪蟲(Trichocerca)、多肢輪蟲(Polyarthra)和疣毛輪蟲(Synchaeta)等; 第三為捕食性輪蟲,如肉食性的晶囊輪蟲(Asplanchna)等.

    應(yīng)用配對t檢驗法檢測每個湖泊非生物因子和生物因子在2個樣點間的差異性. 先對數(shù)據(jù)正態(tài)性進(jìn)行分析(對不符合正態(tài)分布的數(shù)據(jù)進(jìn)行對數(shù)轉(zhuǎn)換),然后在SPSS 16.0軟件中利用Pearson相關(guān)分析法分析2個湖泊中Chl.a和Chl.a<25 μm與各水環(huán)境因子及浮游動物密度間的相關(guān)性. 在R 2.13.0軟件中利用relweights函數(shù)來獲得具有相關(guān)性的被選環(huán)境因子對Chl.a濃度影響的相對權(quán)重[29],將Chl.a和Chl.a<25 μm分別與影響其季節(jié)變動相對權(quán)重最大的水環(huán)境因子進(jìn)行回歸分析. 另外,考慮到氮、磷濃度對于浮游植物生長的重要性,分別將各湖泊中兩類Chl.a濃度與水體氮、磷濃度進(jìn)行回歸分析.

    2 結(jié)果與分析

    2.1 鏡湖水體Chl.a濃度及理化因子的周年動態(tài)

    鏡湖2個樣點中Chl.a和Chl.a<25 μm的年均值分別為27.8和24.0 μg/L、25.2和21.4 μg/L. Chl.a和Chl.a<25 μm在兩樣點間均無顯著差異(P>0.05),其濃度的季節(jié)變動趨勢均表現(xiàn)為夏季>秋季>春季>冬季(圖2).

    鏡湖2個樣點的水溫、DO、pH、SD均不存在顯著性差異(P>0.05). 冬季水溫最低,夏季最高; 而DO濃度的周年變動表現(xiàn)為冬季高、夏季低. 水體pH值在7.0~8.5之間波動; SD的周年動態(tài)與Chl.a和水溫的周年動態(tài)總體變化趨勢均相反(圖3).

    圖2 鏡湖中兩種類型的Chl.a濃度的周年動態(tài)Fig.2 Annual dynamics of two fractions of Chl.a (Chl.a and Chl.a<25 μm) in Lake Jinghu

    圖3 鏡湖水溫、水體pH值、透明度和溶解氧濃度的周年動態(tài)Fig.3 Annual dynamics of water temperature, pH value and secchi depth and DO concentration in Lake Jinghu

    圖4 鏡湖和汀棠湖TN、TP濃度及TN/TP比的周年變動Fig.4 Annual dynamics of TN, TP concentrations and TN/TP ratio in Lake Jinghu and Lake Tingtang

    鏡湖TP濃度表現(xiàn)為秋季>夏季, 而春季與冬季的差別并不大. TN濃度表現(xiàn)為冬季最高, 而其余3個季節(jié)的差別也不大. TN/TP比的年均值為13.68, 整體上表現(xiàn)為冬季>春季>夏季>秋季(圖4).

    2.2 汀棠湖水體Chl.a濃度及理化因子的周年動態(tài)

    汀棠湖2個樣點的Chl.a和Chl.a<25 μm的年均值分別為25.7和25.1 μg/L、25.9和21.0 μg/L, 2個樣點間均沒有顯著差異(P>0.05), 均表現(xiàn)為夏季>春季>秋季>冬季(圖5);水溫、DO、pH、SD、N和P濃度在2個樣點也均不存在顯著差異(P>0.05), 且它們的周年變動趨勢與其在鏡湖中的變動均非常類似(圖4,圖6).

    圖5 汀棠湖兩種類型Chl.a濃度的周年動態(tài)Fig.5 Annual dynamics of two fractions of Chl.a (Chl.a and Chl.a<25 μm) in Lake Tingtang

    圖6 汀棠湖水溫、水體pH值、透明度和溶解氧濃度的周年動態(tài)Fig.6 Annual dynamics of water temperature, pH value, secchi depth and DO concentration in Lake Tingtang

    2.3 鏡湖和汀棠湖輪蟲和枝角類密度的周年動態(tài)

    汀棠湖中存在較多的枝角類和橈足類,其中枝角類主要有象鼻溞(Bosminaspp.)和秀體溞(Diaphanosomaspp.),均為小型枝角類. 象鼻溞的密度高峰期出現(xiàn)在2009年2月末至4月初; 秀體溞的密度高峰期出現(xiàn)在夏季(2009年6、7和8月). 與枝角類和橈足類密度相比,汀棠湖中輪蟲密度相對較高. 汀棠湖2個樣點中,微型濾食輪蟲密度占輪蟲總密度比例均較高,兩者均呈極顯著正相關(guān)(r=0.98,P<0.001). 輪蟲總密度和微型濾食輪蟲在冬季(2009年2月和2010年1月)密度較高,而在夏季(2009年6、7、8月)較低,春末和秋末密度最低.

    鏡湖微型濾食性輪蟲的類型與汀棠湖相同. 微型濾食性輪蟲密度在2個樣點中的變動趨勢與輪蟲總密度大體相似,且兩者呈極顯著相關(guān)(r=0.82,P<0.001). 在2008年9月末-10月末和2009年3月末-4月末,輪蟲總密度和微型濾食性輪蟲密度間均呈現(xiàn)逐步下降的趨勢.

    2.4 兩種類型的Chl.a濃度與非生物因子和生物因子間的關(guān)系

    從全年整體看,2個湖泊中的Chl.a和Chl.a<25 μm間均呈極顯著正相關(guān). 兩類Chl.a濃度均與水溫呈顯著正相關(guān),而與SD均呈顯著負(fù)相關(guān). 兩類Chl.a濃度均與鏡湖DO呈顯著負(fù)相關(guān),均與鏡湖輪蟲總密度和濾食性輪蟲密度呈顯著正相關(guān),與汀棠湖秀體溞密度呈顯著正相關(guān),但與汀棠湖輪蟲的密度均不存在顯著的相關(guān)性(表1).

    圖7 鏡湖和汀棠湖兩類Chl.a濃度與TP濃度,TN/TP比的關(guān)系Fig.7 Regression relationships between two fractions of Chl.a (Chl.a and Chl.a<25 μm) and TP as well as TN/TP ratio in Lake Jinghu and Lake Tingtang

    環(huán)境因子縮寫鏡湖汀棠湖Chl.aChl.a<25μmChl.aChl.a<25μmChl.a/(μg/L)Chl.a- 0.97**- 0.97**水溫/℃WT0.74**0.75**0.54**0.56**透明度/mSD-0.77**-0.78**-0.69*-0.70**pHpH0.010.02-0.22-0.18DO/(mg/L)DO-0.62**-0.63**-0.21-0.21TN/(mg/L)TN-0.15-0.05-0.35-0.36TP/(mg/L)TP0.50*0.63**0.050.11TN/TP比TN/TP-0.63**-0.66**-0.34-0.40象鼻溞密度/(ind./L)Bos- - -0.04 0 秀體溞密度/(ind./L)Dia- - 0.43**0.41**微型濾食輪蟲密度/(ind./L)Filter0.49**0.51**-0.14-0.16輪蟲總密度/(ind./L)Total0.44**0.46**-0.15-0.16

    *表示P<0.05,**表示P<0.01; 汀棠湖中浮游動物密度經(jīng)過lg轉(zhuǎn)化; 下同.

    鏡湖中,兩類Chl.a濃度與TP濃度均呈顯著或極顯著正相關(guān),而與TN/TP比均呈極顯著負(fù)相關(guān)(表1). lg Chl.a(Y1)和lg Chl.a<25 μm(Y2)均隨著lg TP(X1)的增加而顯著增加,回歸關(guān)系分別為Y1=1.17X1+2.78(R2=0.48,P<0.01),Y2=1.38X1+2.92(R2=0.59,P<0.01).Y1和Y2均隨著lg(TN/TP)(X2)的增加而顯著減小(表1),回歸關(guān)系分別為Y1=-1.01X2+2.53(R2=0.48,P<0.01),Y2=-1.14X2+2.52(R2=0.50,P<0.01)(圖7).

    汀棠湖中,兩類Chl.a濃度與TN、TP和TN/TP比間均沒有顯著相關(guān)性(表1). 然而,lg Chl.a(Y3)、lg Chl.a<25 μm(Y4)均隨著lg(TN/TP)(X3)的增加而減小的趨勢,回歸關(guān)系分別為Y3=-0.42X3+1.82(R2=0.13,P=0.09),Y4=-0.59X3+1.89(R2=0.16,P=0.056)(圖7).

    從季節(jié)來看,2個湖泊中兩類型的Chl.a濃度在每個樣點間也呈現(xiàn)極顯著正相關(guān). 鏡湖中,冬季與非冬季的水體中兩種類型的Chl.a濃度均與水溫呈極顯著正相關(guān),而與SD呈極顯著負(fù)相關(guān). 冬季鏡湖中2種類型的Chl.a濃度均與TN濃度呈顯著正相關(guān),而與濾食性輪蟲密度和輪蟲總密度呈極顯著負(fù)相關(guān); 其他3個季節(jié)兩種類型的Chl.a濃度均與DO濃度呈極顯著負(fù)相關(guān),而與濾食性輪蟲密度和輪蟲總密度呈極顯著正相關(guān)(表2).

    汀堂湖中,冬季與非冬季的水體中兩種類型的Chl.a濃度均與SD呈極顯著負(fù)相關(guān). 冬季汀堂湖中兩種類型的Chl.a濃度均與DO濃度呈顯著正相關(guān). 非冬季水體中兩類Chl.a濃度均與TN濃度呈負(fù)相關(guān)趨勢(rChl.a=-0.41,P=0.08;rChl.a<25 μm=-0.44,P=0.07). 盡管冬季與非冬季的水體中兩種類型的Chl.a濃度均與輪蟲總密度無顯著相關(guān)性,但冬季水體中Chl.a<25 μm與濾食性輪蟲密度間呈負(fù)相關(guān)(r=-0.49,P=0.057)(表2).

    各水環(huán)境因子對2個湖泊Chl.a濃度影響的相對權(quán)重見表3. 從全年整體和非冬季的情況看,鏡湖中兩種類型的Chl.a濃度均與SD的關(guān)系最為密切(表3). 但從冬季的情況看,輪蟲對藻類的牧食效應(yīng)(下行效應(yīng))對鏡湖中兩種類型的Chl.a濃度,尤其是Chl.a<25 μm(小型藻類的Chl.a濃度)的影響最大(表3). 汀棠湖中,無論是從全年整體看,還是分季節(jié)看,兩種類型的Chl.a濃度均與SD的關(guān)系最為密切(表2,表3). 但從冬季情況看,輪蟲對藻類的牧食效應(yīng)還是對Chl.a<25 μm產(chǎn)生了一定的影響(表3).

    表2 冬季與非冬季鏡湖和汀棠湖水體中兩種類型的Chl.a濃度與各環(huán)境因子間的相關(guān)系數(shù)

    各環(huán)境因子的縮寫見表1. 下同.

    鑒于透明度與2個湖泊Chl.a濃度間的關(guān)系均最為密切,故對它們之間的關(guān)系進(jìn)行回歸分析. 2個湖泊兩種類型的Chl.a濃度(Chl.a:Y1; Chl.a<25 μm:Y2)均與透明度(X)呈極顯著相關(guān)(圖8).

    表3 各環(huán)境因子對鏡湖和汀棠湖兩種類型Chl.a濃度變動影響的相對權(quán)重(%)

    圖8 鏡湖和汀棠湖兩種類型的Chl.a濃度與透明度間的回歸關(guān)系Fig.8 Regression relationships between two fractions of Chl.a (Chl.a and Chl.a<25 μm) and secchi depth in Lake Jinghu and Lake Tingtang

    3 討論

    3.1 水環(huán)境因子對鏡湖和汀棠湖Chl.a濃度影響的探討

    3.1.1 水體理化因子對Chl.a濃度的影響 眾多研究表明,水溫與Chl.a濃度呈正相關(guān)[6,11,23,30]. 與上述結(jié)論一致,本研究中,鏡湖和汀棠湖Chl.a濃度均與水溫呈顯著正相關(guān),說明水溫對浮游植物的生長和發(fā)育具有促進(jìn)作用. 在當(dāng)前全球氣候變暖的背景下,這種現(xiàn)象可引起環(huán)保工作者對水體富營養(yǎng)化成因的重視,因為水溫的升高既促進(jìn)了Chl.a濃度的上升,也有利于藍(lán)藻形成優(yōu)勢[31].

    有研究表明,水體Chl.a濃度與其pH值間呈顯著正相關(guān)[32]. 在浮游植物生長茂盛的水域,水體pH值較高,主要與浮游植物的光合作用吸收水體中的CO2、放出O2有關(guān). 本研究中,2個湖泊水體pH值均在7.0~8.5 之間波動,整體變化幅度很小,且無論從全年整體還是分季節(jié)看,Chl.a濃度與pH值間均沒有顯著相關(guān)性,這與阮曉紅等的結(jié)論[32]并不一致,而與范廷玉等[33]的研究一致. 究其原因,可能是阮曉紅等[32]的研究只在冬、夏季各采集1次樣品,且60個湖泊間營養(yǎng)水平相差較大有關(guān). 實際上,就單個水體pH的季節(jié)變動不大的湖泊來說,pH與水體藻類繁殖之間的關(guān)系可能不太大.

    Chl.a濃度與DO濃度可能存在著正相關(guān)性,因為Chl.a濃度越高,浮游植物數(shù)量越多,其通過光合作用釋放O2,會使DO濃度增加[30,33]. 但也有研究表明,Chl.a濃度與DO濃度卻呈負(fù)相關(guān)[4]. 本研究中,Chl.a濃度與DO濃度間的關(guān)系在鏡湖呈負(fù)相關(guān),而在冬季汀棠湖呈正相關(guān). 由于DO除了與浮游植物的光合作用放氧有關(guān)外,與水溫、氧分壓和有機體的呼吸消耗等密切相關(guān);一定范圍內(nèi),溫度越高,DO濃度越低[34]. 從本研究來看,汀棠湖和鏡湖均屬小型淺水湖泊,水體中的DO濃度受溫度的影響從而成為影響Chl.a濃度變動的被動因子的可能性更大.

    Chl.a濃度與營養(yǎng)鹽間存在著復(fù)雜的關(guān)系[35]. OECD認(rèn)為,氮、磷濃度均有可能影響Chl.a濃度[2]. 從國內(nèi)學(xué)者對我國湖泊和水庫等的研究發(fā)現(xiàn),Chl.a濃度與氮、磷濃度的關(guān)系也比較復(fù)雜,如對洞庭湖[4]、太湖[8]、蘇州平原河網(wǎng)區(qū)淺水湖泊[32]和珠江流域東江干流[36]的研究中均認(rèn)為P是影響Chl.a濃度變化的主導(dǎo)因子,而對三峽水庫香溪河庫灣的研究中卻發(fā)現(xiàn)Chl.a濃度受N的影響更大[5]. 有學(xué)者認(rèn)為,TN/TP比與藻類的生長有更直接的關(guān)系. 當(dāng)TN/TP比>7或<7時,P和N分別是可能的限制性營養(yǎng)鹽[37]. 本研究中,鏡湖中兩類Chl.a濃度與TP濃度均呈顯著正相關(guān),而與TN/TP比均呈顯著負(fù)相關(guān),且TN/TP比年均值為13.68,這表明鏡湖屬于P限制性水體,P對于Chl.a濃度的影響較大. 而在汀棠湖中,兩類Chl.a濃度均分別與TN/TP比和TN濃度呈負(fù)相關(guān),且TN/TP比的年均值為12.09,這表明汀棠湖可能偏向P限制性水體. 但汀棠湖中Chl.a濃度與TP濃度并沒有呈現(xiàn)正相關(guān)性,可能與汀棠湖在采樣期地處蕪湖市郊,周邊居民在水體放養(yǎng)較多的鰱鳙魚有關(guān),因為這些濾食性魚類對于藻類的濾食可能改變了Chl.a濃度與TP濃度的關(guān)系[38]. 當(dāng)然,這種推測還需要更進(jìn)一步的研究予以證實.

    3.1.2 浮游動物對Chl.a濃度的影響 輪蟲、枝角類等浮游動物均攝食浮游藻類,從而有可能會影響到水體中的Chl.a濃度. 本研究中,藻類的Chl.a總濃度與小型藻類的Chl.a濃度間呈極顯著相關(guān),且它們的均值差異都不大,這說明2個湖泊中大型藻類所占比例較小. 從全年整體來看,鏡湖的兩類Chl.a濃度與輪蟲密度均呈極顯著正相關(guān),一方面說明了浮游動物對于鏡湖Chl.a濃度的影響不大,另一方面也暗示了水體中可供輪蟲利用的食物資源是充足的. 汀棠湖中,兩類Chl.a濃度與輪蟲密度并未呈顯著相關(guān),這說明汀棠湖浮游動物的牧食對于水體Chl.a濃度的影響很小.

    但是從季節(jié)變化的特點看,在水溫最低、光照時間最短的冬季,水體透明度均最高,2個湖泊中微型濾食性輪蟲總密度和輪蟲總密度均較高,而Chl.a濃度卻最低. 此階段,輪蟲密度與Chl.a濃度呈顯著負(fù)相關(guān)性(鏡湖)或負(fù)相關(guān)趨勢(汀棠湖),且相比其他水環(huán)境因子而言,輪蟲密度對于鏡湖中藻類Chl.a總濃度,尤其對小型藻類的Chl.a濃度影響最大; 而這些現(xiàn)象在溫度較高、光照時間較長的非冬季中并沒有顯現(xiàn),這說明在低溫、短光照下的冬季水體中,藻類光合作用較弱,浮游輪蟲攝食藻類降低了Chl.a濃度. 但是,輪蟲密度對于汀棠湖Chl.a濃度的影響卻有限,一方面可能與汀棠湖樣品采集的頻率相對較低有關(guān),因為較長的采樣周期難以反映出對水環(huán)境變化敏感的輪蟲等浮游動物對藻類攝食的影響力; 另一方面也說明了浮游動物攝食藻類在影響小型淺水湖泊中Chl.a濃度變動的重要性方面需要更多的調(diào)查予以確認(rèn). 結(jié)合本實驗和張鎮(zhèn)等[19]的實驗結(jié)果,有理由相信,通過一定的措施增加浮游動物密度,將有可能有效地降低浮游植物的生物量和水華產(chǎn)生的可能性.

    3.2 小型封閉性水體與大型淺水湖泊、水庫和河流Chl.a濃度及其變動因素的比較

    大型淺水湖泊由于受風(fēng)浪的攪動強烈,其Chl.a濃度的變動除了與水溫及營養(yǎng)鹽等相關(guān)外,還與沉積物懸浮導(dǎo)致的濁度上升有關(guān),湖水渾濁使得水下光照成為影響水體Chl.a濃度變動的限制性因子[10]. 有些大型淺水湖泊(如鄱陽湖)并非封閉性水體,與外界的交換程度高,高水位的豐水期會稀釋營養(yǎng)鹽濃度,從而改變了其對水體Chl.a濃度的影響[9]. 對水庫的研究發(fā)現(xiàn),其中的Chl.a濃度除了受營養(yǎng)鹽的影響程度較大外,還與水庫的水文條件相關(guān),調(diào)水會影響Chl.a濃度與營養(yǎng)鹽間的關(guān)系[39]. 就河流而言,不同河段的Chl.a濃度往往不同,與水體的流動性也具有緊密的關(guān)系[36,40]. 這些現(xiàn)象說明沉積物的懸浮、水文條件和水體流動性是影響大型淺水湖泊、水庫或者河流中Chl.a濃度變動的不容忽視的因素,同時水體的不同區(qū)域受這些因素影響程度可能不同,導(dǎo)致了Chl.a濃度的分布具有空間異質(zhì)性. 至于小型淺水水體,其Chl.a濃度的分布是否具有空間異質(zhì)性與水體周圍的環(huán)境條件相關(guān). 就本研究而言,鏡湖和汀棠湖兩個樣點的Chl.a濃度間均沒有顯著差異,這可能是因為它們均屬城市內(nèi)湖,且與周邊的水交換程度均不高. 但是,小型淺水水體即使是封閉的,如果不同水域受風(fēng)的影響程度不同,也會影響到Chl.a濃度的空間分布特征[41].

    3.3 透明度與Chl.a濃度間關(guān)系的探討

    一般而言,透明度與該水體中懸浮物質(zhì)和Chl.a濃度均呈負(fù)相關(guān)[5,6,30,42]. 且有學(xué)者認(rèn)為,Chl.a濃度與透明度間呈現(xiàn)曲線關(guān)系[5,34]. 本研究中,無論從全年整體上看還是分季節(jié)看,除了冬季的鏡湖外,2個湖泊的Chl.a濃度均與透明度間的關(guān)系最為密切(圖8). 然而,有不少學(xué)者認(rèn)為,透明度可能是引起Chl.a濃度變動的被動因子[33],但作者認(rèn)為也有可能是誘導(dǎo)因素. 本研究的2個湖泊均在城市內(nèi)部,易受周邊人群的影響; 且它們均屬于小型淺水湖泊,其透明度易受地表水沖洗、雨水淋刷和人為活動等外源性因素的影響,這些因素可能會主動影響到水體中的浮游藻類的光合作用. 但從本研究并結(jié)合相關(guān)文獻(xiàn)看[5,34],不管透明度是作為Chl.a濃度變動的誘導(dǎo)因子還是被動因子,由于測定過程非常簡便,其可能是預(yù)測封閉性較好、與外界水交換程度低的水體Chl.a濃度變動的較為簡便的指標(biāo). 不過,由于透明度測定的過程受人為影響程度較大,因此就需要盡可能地提高測定的準(zhǔn)確性,比如: 選用黑白對照清晰的透明度盤、盡量選擇同一個人完成全部樣品的測定工作、保持透明度盤在水中的垂直狀態(tài)等減小測定誤差的措施. 當(dāng)然,在大型湖泊、水庫、河流及受風(fēng)影響較大的小型水體中,由于水體的流動性或者水位的變化較大,是否能夠利用透明度來預(yù)測Chl.a濃度的變動需要更多、更細(xì)致的研究加以證實.

    [1] Reynolds CS ed. The ecology of freshwater phytoplankton. London: Cambridge University Press, 1984.

    [2] Organization for Economic Co-operation and Development (OECD). Eutrophication of waters: Monitoring, assessment and control. Paris: OECD, 1982.

    [3] Xun Shangpei, Yang Yuanjian, He Binfangetal. Analysis of relationship between spring water temperature and chlorophyll-a in Lake Chaohu.JLakeSci, 2011,23(5): 767-772. DOI:10.18307/2011.0515. [荀尚培, 楊元建, 何彬方等. 春季巢湖水溫和水體葉綠素a濃度的變化關(guān)系. 湖泊科學(xué), 2011,23(5): 767-772.]

    [4] Zhang Guanggui. Spatial-temporal distribution of chlorophyll-a and its correlation with environment factors in Dongting Lake.EnvironmentalMonitoringinChina, 2016,32(4): 84-90. [張光貴. 洞庭湖水體葉綠素a時空分布及與環(huán)境因子的相關(guān)性. 中國環(huán)境監(jiān)測, 2016,32(4): 84-90.]

    [5] Han Xinqin, Ye Lin, Xu Yaoyangetal. Analysis of the spatial and temporal changes of chlorophyll-a concentration in Xiangxi Bay in spring and its impact factors.ActaHydrobiologicaSinica, 2006,30(1): 89-94. [韓新芹, 葉麟, 徐耀陽等. 香溪河庫灣春季Chl.a濃度動態(tài)及其影響因子分析. 水生生物學(xué)報, 2006,30(1): 89-94.]

    [6] Lu Huanchun, Wang Feier, Chen Yingxuetal. Multianalysis between chlorophyll-a and environmental factors in Qiandao Lake water.ChineseJournalofAppliedEcology, 2003,14(8): 1347-1350. [呂喚春, 王飛兒, 陳英旭等. 千島湖水體葉綠素a濃度與相關(guān)環(huán)境因子的多元分析. 應(yīng)用生態(tài)學(xué)報, 2003,14(8): 1347-1350.]

    [7] Liang Shuxuan, Wang Yunxiao, Qin Zhe. Chlorophyll-a and analysis of water quality of Lake Baiyangdian.TransactionsofOceanologyandLimnology, 2012,3: 66-73. [梁淑軒, 王云曉, 秦哲. 白洋淀葉綠素a及其水質(zhì)因子分析. 海洋湖沼通報, 2012,3: 66-73.]

    [8] Wang Zhen, Zou Hua, Yang Guijunetal. Spatial- temporal characteristics of chlorophyll-a and its relationship with environmental factors in Lake Taihu.JLakeSci, 2014,26(4): 567-575. DOI:10.18307/2014.0411. [王震, 鄒華, 楊桂軍等. 太湖葉綠素a的時空分布特征及其與環(huán)境因子的相關(guān)關(guān)系. 湖泊科學(xué), 2014,26(4): 567-575.]

    [9] Wu Zhaoshi, Zhang Lu, Liu Baoguietal. Spatial distribution of chlorophyll-a in Poyang Lake during wet season and its relationship with environmental factors.WetlandScience, 2014,12(3): 286-292. [吳召仕, 張路, 劉寶貴等. 鄱陽湖豐水期水體中葉綠素a濃度空間分布及其與環(huán)境因子的關(guān)系. 濕地科學(xué), 2014,12(3): 286-292.]

    [10] Thomas JR, Havens K, Zhu GWetal. Comparative analysis of nutrients, chlorophyll and transparency in two large shallow lakes (Lake Taihu, P.R. China and Lake Okeechobee, USA).Hydrobiologia, 2009,627: 211-231.

    [11] Zhou Hong. Quantity of chlorophyll-a and water quality in the West Lake, Hangzhou.JournalofZhejiangUniversity:ScienceEdition, 2001,28(4): 439-442. [周宏. 杭州西湖水體中葉綠素a含量與水質(zhì)的關(guān)系. 浙江大學(xué)學(xué)報: 理學(xué)版, 2001,28(4): 439-442.]

    [12] Hu Yufei, Yu Dezhao, Guo Longgenetal. Relationships of water taste and odor compounds and their related environmental factors in Lake Donghu, Wuhan.JLakeSci, 2017,29(1): 87-94. DOI:10.18307/2014.0411. [胡宇飛, 余得昭, 過龍根等. 武漢東湖水體異味物質(zhì)及其與水環(huán)境因子相互關(guān)系. 湖泊科學(xué), 2017,29(1): 87-94.]

    [13] Yao Yuan, He Feng, Hu Shenghuaetal. Effects of allelopathy of submerged macrophytes on the phytoplankton community collected from the west part of the West Lake wetland in Hangzhou, China.ActaEcologicaSinica, 2016,36(4): 971-978. [姚遠(yuǎn), 賀鋒, 胡勝華等. 沉水植物化感作用對西湖濕地浮游植物群落的影響. 生態(tài)學(xué)報, 2016,36(4): 971-978.]

    [14] Zhang Zongshe, Huang Xiangfei eds. Research methods of freshwater plankton. Beijing:Science Press, 1991. [章宗涉, 黃祥飛. 淡水浮游生物研究方法. 北京: 科學(xué)出版社, 1991.]

    [15] Pourriot R. Food and feeding habits of Rotifera.ArchivfürHydrobiologieBeiheftErgebnissederLimnologie, 1977,8: 243-260.

    [16] De Bernardi R, Giussani G, Manca M. Cladocera: predators and prey.Hydrobiologia, 1987,145: 225-243.

    [17] Yang Yufeng, Huang Xiangfei. Advances in ecological studies on zooplankton.JLakeSci, 2000,12(1): 82-89. DOI:10.18307/2000.0113. [楊宇峰, 黃祥飛. 浮游動物生態(tài)學(xué)研究進(jìn)展. 湖泊科學(xué), 2000,12(1): 82-89.]

    [18] Orcutt JD, Pace ML. Seasonal dynamics of rotifer and crustacean zooplankton populations in a eutrophic, monomictic lake with a note on rotifer sampling techniques.Hydrobiologia, 1984,119: 73-80.

    [19] Zhang Zhen, Chen Feizhou, Zhou Wanpingetal. Effects ofDaphniacarinatagrazing on the phytoplankton of Lake Xuanwu, Nanjing.JLakeSci, 2009,21(3): 415-419. DOI:10.18307/2009.0316. [張鎮(zhèn), 陳非洲, 周萬平等. 南京玄武湖隆腺溞(Daphniacarinata)牧食對浮游植物的影響. 湖泊科學(xué), 2009,21(3): 415-419.]

    [20] Zhou Kunpeng, Liu Yangchun, Wang Chongchen. Analysis on temporal-spatial variation of eutrophication in Shichahai area, Beijing.EnvironmentalChemistry,35(4): 703-712. [周坤朋, 劉陽春, 王崇臣. 北京什剎海區(qū)域水體富營養(yǎng)化時空演變特征分析. 環(huán)境化學(xué), 2016,35(4): 703-712.]

    [21] Wen XL, Xi YL, Qian FPetal. Comparative analysis of rotifer community structure in five subtropical shallow lakes in East China: role of physical and chemical conditions.Hydrobiologia, 2011,661: 303-316.

    [22] Qian Fangping, Xi Yilong, Wen Xinlietal. Eutrophication impact on community structure and species diversity of rotifers in five lakes of Anhui.BiodiversityScience, 2007,15(4): 344-355. [錢方平, 席貽龍, 溫新利等. 湖泊富營養(yǎng)化對輪蟲群落結(jié)構(gòu)及物種多樣性的影響. 生物多樣性, 2007,15(4): 344-355.]

    [23] Wen Xinli, Xi Yilong, Zhang Leietal. Analysis of community structure of rotifera and ecological assessment of water quality in Lake Jinghu, Wuhu city.ActaHydrobiologicaSinica, 2006,30(2): 152-158. [溫新利, 席貽龍, 張雷等. 蕪湖市鏡湖輪蟲群落結(jié)構(gòu)分析及水質(zhì)生態(tài)學(xué)評價. 水生生物學(xué)報, 2006,30(2): 152-158.]

    [24] Huang Xiangfei ed. Survey, observation and analysis of Lake Ecology. Beijing: Chinese Standard Press, 1999. [黃祥飛. 湖泊生態(tài)調(diào)查觀測與分析. 北京: 科學(xué)出版社, 1999.]

    [25] Koste W ed. Rotatoria: Die R?dertiere Mitteleuropas. Berlin: Gebrüder Borntraeger, 1978.

    [26] Jiang Xiezhi, Du Nanshan eds. Fauna of China animals—Freshwater Cladocer. Beijing: Science Press, 1979. [蔣燮治, 堵南山. 中國動物志——淡水枝角類. 北京: 科學(xué)出版社, 1979.]

    [28] Obertegger U, Smith HA, Flaim Getal. Using the guild ratio to characterize pelagic rotifer communities.Hydrobiologia, 2011,662: 157-162.

    [29] R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria: ISBN 3-900051-07-0. Available at: http://www. R-project. org., 2011.

    [30] Ge Dabing, Wu Xiaoling, Zhu Weilinetal. Chlorophyll-a and its relationship with water quality in Southlake, Yueyang city.EnvironmentalMonitoringinChina, 2005,21(4): 69-71. [葛大兵, 吳小玲, 朱偉林等. 岳陽南湖Chl.a濃度及其水質(zhì)關(guān)系分析. 中國環(huán)境監(jiān)測,21(4): 69-71.]

    [31] Lü H, Yang J, Liu LMetal. Temperature and nutrients are significant drivers of seasonal shift in phytoplankton community from a drinking water reservoir, subtropical China.EnvironmentalScienceandPollutionResearch, 2014,21:5917-5928.

    [32] Ruan Xiaohong, Shi Xiaodan, Zhao Zhenhuaetal. Correlation between chlorophyll-a concentration and environmental factors in shallow lakes in plain river network areas of Suzhou.JLakeSci, 2008,20(5): 556-562. DOI:10.18307/2008.0502. [阮曉紅, 石曉丹, 趙振華等. 蘇州平原河網(wǎng)區(qū)淺水湖泊葉綠素a濃度與環(huán)境因子的相關(guān)關(guān)系. 湖泊科學(xué), 2008,20(5): 556-562.]

    [33] Fan Tingyu, Cheng Fangkui, Yan Jiapingetal. Distribution and correlation factors of chlorophyll-a in Linhuan coal mining subsidence water in Huaibei.EnvironmentalChemistry, 2015,34(6): 1168-1176. [范廷玉, 程方奎, 嚴(yán)家平等. 淮北臨渙煤礦塌陷水域葉綠素a濃度與相關(guān)因子分析. 環(huán)境化學(xué), 2015,34(6): 1168-1176.]

    [34] Liu Jiankang ed. Advanced aquatic biology. Beijing: Science Press, 1999. [劉建康. 高級水生生物學(xué). 北京: 科學(xué)出版社, 1999.]

    [35] Hu Ren, Lin Qiuqi, Duan Shunshanetal. Distribution of chlorophyll-a and phosphorus in subtropical reservoirs, South China.EcologicalScience, 2002,21(4): 310-315. [胡韌, 林秋奇, 段舜山等. 熱帶亞熱帶水庫浮游植物Chl.a濃度與磷分布的特征. 生態(tài)科學(xué), 2002,21(4): 310-315.]

    [36] Tian Shimi, Yang Yang, Qiao Yongminetal. Temporal and spatial distribution of phytoplankton chlorophyll-a and its relationships with environmental factors in Dongjiang River, Pearl River basin.JLakeSci, 2015,27(1): 31-37. DOI:10.18307/2015.0104. [田時彌, 楊揚, 喬永民等. 珠江流域東江干流浮游植物葉綠素a濃度時空分布及與環(huán)境因子的關(guān)系. 湖泊科學(xué), 2015,27(1): 31-37.]

    [37] Liu X, Lu XH, Chen YW. The effects of temperature and nutrient ratios onMicrocystisblooms in Lake Taihu, China: An 11-year investigation.HarmfulAlgae, 2011,10(3): 337-343.

    [38] Li Peipei, Shi Wen, Liu Qigenetal. Spatial and temporal distribution patterns of chlorophyll-a and the correlation analysis with environmental factors in Lake Qiandao.JLakeSci, 2011,23(4): 568-574. DOI:10.18307/2011.0412. [李培培, 史文, 劉其根等. 千島湖葉綠素a的時空分布及其與影響因子的相關(guān)分析. 湖泊科學(xué), 2011,23(4): 568-574.]

    [39] Wang Wei, Gu Jiguang, Han Boping. Analysis of factors affecting chlorophyll-a concentration in small reservoirs in South China.ChineseJournalofAppliedEnvironmentalBiology, 2009,15(1): 64-71. [王偉, 顧繼光, 韓博平. 華南沿海地區(qū)小型水庫葉綠素a濃度的影響因子分析. 應(yīng)用與環(huán)境生物學(xué)報, 2009,15(1): 64-71.]

    [40] Zhou Beibei, Wang Guoxiang, Xu Yaoetal. Spatial distribution of chlorophyll-a and its relationships with environmental factors in Qin-huai River, Nanjing.JLakeSci, 2012,24(2): 267-272. DOI:10.18307/2012.0215. [周貝貝, 王國祥, 徐瑤等. 南京秦淮河葉綠素a濃度空間分布及其與環(huán)境因子的關(guān)系. 湖泊科學(xué), 2012,24(2): 267-272.]

    [41] Li Feipeng, Zhang Haiping, Chen Ling. Temporal and spatial distribution of environmental factors and chlorophyll-a and their correlation analysis in a small enclosed lake.EnvironmentalScience, 2013,34(10): 3854-3861. [李飛鵬, 張海平, 陳玲. 小型封閉水體環(huán)境因子與葉綠素的時空分布及相關(guān)性研究. 環(huán)境科學(xué), 2013,34(10): 3854-3861.]

    [42] Chen Y, Qin B, Teubner K. Long-term dynamics of phytoplankton assemblages:Miyrocystis-domination in Lake Taihu, a large shallow lake in China.JournalofPlanktonResearch, 2003,25: 445-453.

    Temporalvariationsofchlorophyll-aconcentrationandtheirrelationshipswithabioticandbioticfactorsintwosmallshallowlakes

    WEN Xinli1,2, FENG Ruonan1,2, ZHANG Gen1, XUE Yinghao1,3& XI Yilong1,2**

    (1:CollegeofLifeSciences,AnhuiNormalUniversity,KeyLaboratoryofBioticEnvironmentandEcologicalSafetyinAnhuiProvince,Wuhu241000,P.R.China)(2:CollaborativeInnovationCenterofRecoveryandReconstructionofDegradedEcosysteminWanjiangCityBelt,Wuhu241000,P.R.China)(3:RuralEnergy&EnvironmentAgency,MinistryofAgriculture,Beijing100125,P.R.China)

    Identifying the factors affecting the variations of chlorophyll-a concentration is very important in assessment and management of water quality in lakes. Most of previous studies focused on the relationships between chlorophyll-a concentration and abiotic factors, while paid little attention on the effects of biotic ones. To examine the relative importance of the abiotic and biotic factors on seasonal variation of chlorophyll-a in shallow lakes, water samples were collected semimonthly from February 2009 to January 2010 in Lake Tingtang, and four-day or one-week intervals from May 2008 to April 2009 in Lake Jinghu, respectively. Simultaneously, the abiotic factors including water temperature, pH, dissoloved oxygen, Secchi depth, total nitrogen(TN) and total phosphorus(TP) as well as the densities of rotifers and cladocerans were analyzed. The chlorophyll-a concentration in each water body was measured in two fractions: one filtered through a 25 μm plankton net (Chl.a<25 μm) and another unfiltered (Chl.a). During the monitoring period, the two fractions of chlorophyll-a in two lakes were both positively correlated with water temperature and negatively associated with secchi depth, and they did not show reverse relationship with the densities of rotifers and cladocerans in both lakes. Chl.a and Chl.a<25 μm increased significantly with TP but decreased with TN/TP ratio in Lake Jinghu, while they showed reducing trends with TN and TN/TP ratio in Lake Tingtang. During winter, the rotifer density was negatively associated with Chl.a and Chl.a<25 μm in Lake Jinghu, and the similar relationships between them were found in Lake Tingtang, indicating that the low chlorophyll-a concentration in winter might be partly due to the prey effect of rotifers on phytoplankton. Based on the relative weight of the effect intensity of each correlated factor on the variation of chlorophyll-a concentration, Chl.a and Chl.a<25 μm were both affected mostly by secchi depth. Additionally, the two fractions of chlorophyll-a decreased with an increase in secchi depth. Taken together, the results in this paper implied that Secchi depth might be a simple and convenient index for predicting chlorophyll-a concentration in the small enclosed lakes.

    Chlorophyll-a concentration; seasonal variation; biotic factor; abiotic factor; relative importance; Lake Tingtang; Lake Jinghu

    *國家自然科學(xué)基金項目(31470015)、安徽省高校自然科學(xué)研究重點項目(KJ2017A320)和重要生物資源的保護(hù)與利用安徽省重點實驗室專項基金聯(lián)合資助. 2016-11-08收稿; 2017-03-12收修改稿. 溫新利(1977~), 男, 博士, 副教授; E-mail: wenxinli1977@126.com.

    **通信作者; E-mail: ylxi1965@126.com.

    猜你喜歡
    鏡湖輪蟲淺水
    陳鏡湖烈士紀(jì)念館
    新型淺水浮托導(dǎo)管架的應(yīng)用介紹
    云南化工(2021年10期)2021-12-21 07:33:40
    溫度對不同種類臂尾輪蟲生長與繁殖的影響
    “逆天神偷”蛭形輪蟲
    悄悄地
    “逆天神偷”蛭形輪蟲
    知識窗(2019年4期)2019-04-26 03:16:02
    夜游天鏡湖
    帶阻尼的隨機淺水波方程的隨機吸引子
    (2+1)維廣義淺水波方程的Backlund變換和新精確解的構(gòu)建
    鏡湖:紹興歷史文化的一面鏡子
    国产欧美日韩精品亚洲av| 久久久国产成人免费| 国产精品美女特级片免费视频播放器 | 欧美日本中文国产一区发布| www.www免费av| 女人被狂操c到高潮| 中文字幕av电影在线播放| 亚洲欧美精品综合一区二区三区| 国产真人三级小视频在线观看| 在线天堂中文资源库| 麻豆av在线久日| 狂野欧美激情性xxxx| 真人做人爱边吃奶动态| 99久久久亚洲精品蜜臀av| 亚洲视频免费观看视频| 啪啪无遮挡十八禁网站| 少妇的丰满在线观看| 麻豆一二三区av精品| 亚洲中文字幕日韩| 国产主播在线观看一区二区| 精品免费久久久久久久清纯| 日本撒尿小便嘘嘘汇集6| 一级毛片女人18水好多| 12—13女人毛片做爰片一| 亚洲精品一二三| 亚洲国产欧美一区二区综合| 高清av免费在线| 成人永久免费在线观看视频| 成人三级黄色视频| 亚洲av美国av| tocl精华| 啦啦啦在线免费观看视频4| 国产亚洲精品久久久久久毛片| 男人舔女人下体高潮全视频| 久久亚洲真实| 中文字幕最新亚洲高清| 变态另类成人亚洲欧美熟女 | 日本免费a在线| 国产精品秋霞免费鲁丝片| xxxhd国产人妻xxx| 中文亚洲av片在线观看爽| 国产在线观看jvid| 日韩av在线大香蕉| 亚洲avbb在线观看| 变态另类成人亚洲欧美熟女 | 中亚洲国语对白在线视频| 91麻豆av在线| 一边摸一边抽搐一进一出视频| 欧美av亚洲av综合av国产av| 亚洲aⅴ乱码一区二区在线播放 | 在线免费观看的www视频| av超薄肉色丝袜交足视频| 亚洲第一欧美日韩一区二区三区| 亚洲av五月六月丁香网| 国产一区在线观看成人免费| 老鸭窝网址在线观看| 成熟少妇高潮喷水视频| 变态另类成人亚洲欧美熟女 | 国产黄色免费在线视频| 国产成人av激情在线播放| 欧美成人午夜精品| 一二三四社区在线视频社区8| 桃红色精品国产亚洲av| 欧美最黄视频在线播放免费 | 视频在线观看一区二区三区| 国产人伦9x9x在线观看| 超碰97精品在线观看| 天堂动漫精品| 国产色视频综合| 91大片在线观看| 女同久久另类99精品国产91| 一二三四社区在线视频社区8| 欧美日韩视频精品一区| 亚洲中文日韩欧美视频| 精品久久久久久,| 久久中文字幕人妻熟女| 免费不卡黄色视频| 免费日韩欧美在线观看| 另类亚洲欧美激情| 中国美女看黄片| 成人影院久久| 国产精品综合久久久久久久免费 | 久久人妻福利社区极品人妻图片| 大码成人一级视频| 欧美性长视频在线观看| 日韩 欧美 亚洲 中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 大陆偷拍与自拍| 美女 人体艺术 gogo| 91精品国产国语对白视频| 欧美一区二区精品小视频在线| 天堂√8在线中文| 国产熟女午夜一区二区三区| 丁香六月欧美| 女人高潮潮喷娇喘18禁视频| 99国产精品一区二区蜜桃av| 国产aⅴ精品一区二区三区波| 欧美成人午夜精品| 亚洲成av片中文字幕在线观看| 后天国语完整版免费观看| 国产99久久九九免费精品| 女人被躁到高潮嗷嗷叫费观| 国产成人精品在线电影| 91麻豆av在线| 天堂俺去俺来也www色官网| 免费在线观看完整版高清| 国产精品成人在线| 热re99久久国产66热| 亚洲少妇的诱惑av| cao死你这个sao货| 人人妻人人添人人爽欧美一区卜| 神马国产精品三级电影在线观看 | 97人妻天天添夜夜摸| 亚洲片人在线观看| 午夜免费激情av| xxx96com| 人人妻,人人澡人人爽秒播| 丝袜美腿诱惑在线| 国产高清videossex| 亚洲一区二区三区色噜噜 | 制服人妻中文乱码| ponron亚洲| 成人亚洲精品一区在线观看| av网站在线播放免费| 亚洲精品成人av观看孕妇| 日韩欧美三级三区| 日韩一卡2卡3卡4卡2021年| 黑人操中国人逼视频| 在线视频色国产色| a在线观看视频网站| 50天的宝宝边吃奶边哭怎么回事| 纯流量卡能插随身wifi吗| 亚洲av电影在线进入| 免费不卡黄色视频| 亚洲专区国产一区二区| 久久久水蜜桃国产精品网| 国产精品秋霞免费鲁丝片| 久久精品亚洲熟妇少妇任你| 日韩欧美一区视频在线观看| 中文字幕人妻熟女乱码| 欧美国产精品va在线观看不卡| 国产不卡一卡二| 亚洲一区二区三区欧美精品| 亚洲 欧美一区二区三区| 99国产极品粉嫩在线观看| 性色av乱码一区二区三区2| www.精华液| 黄网站色视频无遮挡免费观看| 一二三四社区在线视频社区8| 精品无人区乱码1区二区| 在线观看舔阴道视频| 成人av一区二区三区在线看| 天堂√8在线中文| 国产区一区二久久| 激情在线观看视频在线高清| 精品国产美女av久久久久小说| 日本五十路高清| 天堂√8在线中文| 国产亚洲欧美98| 久久草成人影院| 中文字幕高清在线视频| 欧美+亚洲+日韩+国产| 久久久久久久午夜电影 | 人人妻人人添人人爽欧美一区卜| 老汉色∧v一级毛片| 90打野战视频偷拍视频| 久9热在线精品视频| 夜夜躁狠狠躁天天躁| 欧美日韩黄片免| 国产成年人精品一区二区 | 激情视频va一区二区三区| 精品国产美女av久久久久小说| 最近最新免费中文字幕在线| 色播在线永久视频| 国产在线观看jvid| 最近最新免费中文字幕在线| 亚洲色图综合在线观看| 亚洲 国产 在线| 久久久久久久精品吃奶| 免费一级毛片在线播放高清视频 | 中文字幕另类日韩欧美亚洲嫩草| 9色porny在线观看| 亚洲专区中文字幕在线| 一边摸一边抽搐一进一小说| av国产精品久久久久影院| 午夜影院日韩av| 免费观看人在逋| 国产成人欧美在线观看| 一边摸一边做爽爽视频免费| 欧美激情 高清一区二区三区| 国产精品久久视频播放| 国产成人精品久久二区二区免费| 村上凉子中文字幕在线| 精品福利永久在线观看| 国产精品电影一区二区三区| 男人舔女人的私密视频| 国产伦人伦偷精品视频| 黄色丝袜av网址大全| 国产视频一区二区在线看| 成年人免费黄色播放视频| 国产麻豆69| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人免费av在线播放| 黄频高清免费视频| 午夜成年电影在线免费观看| 在线十欧美十亚洲十日本专区| 搡老岳熟女国产| 亚洲av熟女| 欧美中文综合在线视频| 亚洲精品国产区一区二| 成在线人永久免费视频| 国内久久婷婷六月综合欲色啪| 久热爱精品视频在线9| 午夜免费观看网址| 国产精品美女特级片免费视频播放器 | 久久中文字幕一级| 精品久久久久久成人av| 不卡一级毛片| 久久久精品欧美日韩精品| svipshipincom国产片| 国产亚洲精品久久久久久毛片| 十八禁网站免费在线| 亚洲五月天丁香| 国产精品 欧美亚洲| 在线观看舔阴道视频| 男女下面进入的视频免费午夜 | 中文字幕最新亚洲高清| 欧美日韩黄片免| 欧美激情极品国产一区二区三区| 日本撒尿小便嘘嘘汇集6| 精品少妇一区二区三区视频日本电影| 人人妻人人爽人人添夜夜欢视频| 一级作爱视频免费观看| 高清欧美精品videossex| 99re在线观看精品视频| 亚洲av美国av| 精品一区二区三区四区五区乱码| 久热爱精品视频在线9| 在线免费观看的www视频| 欧美日韩福利视频一区二区| 老司机深夜福利视频在线观看| 男人舔女人下体高潮全视频| 一区二区三区精品91| 久久 成人 亚洲| ponron亚洲| 国产精品自产拍在线观看55亚洲| 午夜成年电影在线免费观看| 999久久久国产精品视频| 午夜福利在线观看吧| 国产精品一区二区在线不卡| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧美一区二区三区黑人| 日韩欧美在线二视频| 亚洲精品国产区一区二| 久久中文字幕人妻熟女| 久久久国产成人精品二区 | 狠狠狠狠99中文字幕| 国产精品久久电影中文字幕| 欧美+亚洲+日韩+国产| 中文字幕人妻丝袜制服| 成人影院久久| 亚洲国产毛片av蜜桃av| 国产av又大| 在线观看免费高清a一片| 亚洲成人免费av在线播放| 久久香蕉激情| 欧美日韩精品网址| 亚洲av第一区精品v没综合| 高清在线国产一区| 黑人巨大精品欧美一区二区蜜桃| 亚洲人成电影免费在线| 丁香六月欧美| 成人三级黄色视频| 波多野结衣av一区二区av| 午夜日韩欧美国产| 97人妻天天添夜夜摸| 日韩成人在线观看一区二区三区| 村上凉子中文字幕在线| 午夜福利一区二区在线看| 精品人妻1区二区| 国产成人精品在线电影| 亚洲国产精品一区二区三区在线| 无人区码免费观看不卡| 国产伦一二天堂av在线观看| 国产不卡一卡二| 国产黄色免费在线视频| 法律面前人人平等表现在哪些方面| 国产熟女xx| 久久久久国产一级毛片高清牌| 午夜精品在线福利| 亚洲欧美日韩高清在线视频| 老司机深夜福利视频在线观看| 日韩三级视频一区二区三区| 国产在线精品亚洲第一网站| 国产精品国产av在线观看| 久久精品国产99精品国产亚洲性色 | tocl精华| 欧美日韩黄片免| 亚洲一码二码三码区别大吗| 老司机深夜福利视频在线观看| 久久精品人人爽人人爽视色| 美女福利国产在线| 女性生殖器流出的白浆| 午夜成年电影在线免费观看| 日韩精品免费视频一区二区三区| 大码成人一级视频| 啦啦啦免费观看视频1| 免费不卡黄色视频| 最好的美女福利视频网| 亚洲成a人片在线一区二区| 日韩免费av在线播放| 欧美精品一区二区免费开放| 99国产精品99久久久久| 亚洲视频免费观看视频| 美女扒开内裤让男人捅视频| 色在线成人网| 国产国语露脸激情在线看| 天堂俺去俺来也www色官网| 电影成人av| 岛国在线观看网站| av有码第一页| 高潮久久久久久久久久久不卡| 亚洲激情在线av| 久久久国产一区二区| 免费观看精品视频网站| 欧美在线黄色| 精品日产1卡2卡| 天天躁狠狠躁夜夜躁狠狠躁| tocl精华| 午夜视频精品福利| 一级,二级,三级黄色视频| 咕卡用的链子| 国产精品日韩av在线免费观看 | 亚洲欧美一区二区三区久久| 亚洲激情在线av| 国产精品 欧美亚洲| 丝袜在线中文字幕| 欧美日韩福利视频一区二区| 亚洲五月色婷婷综合| 日本免费a在线| av电影中文网址| 欧美激情 高清一区二区三区| www.熟女人妻精品国产| 日韩成人在线观看一区二区三区| 操美女的视频在线观看| 人人妻,人人澡人人爽秒播| a在线观看视频网站| 在线观看66精品国产| 露出奶头的视频| 成人精品一区二区免费| 国产精品免费一区二区三区在线| 久久伊人香网站| 女人高潮潮喷娇喘18禁视频| 国产97色在线日韩免费| 欧美亚洲日本最大视频资源| av福利片在线| 激情在线观看视频在线高清| 免费女性裸体啪啪无遮挡网站| 丰满迷人的少妇在线观看| 大码成人一级视频| 成人亚洲精品av一区二区 | 亚洲午夜理论影院| 久久中文看片网| 国产精华一区二区三区| 美女福利国产在线| 精品电影一区二区在线| 亚洲久久久国产精品| 久久久久国产一级毛片高清牌| 一边摸一边抽搐一进一出视频| 午夜福利在线免费观看网站| 国产亚洲欧美精品永久| 久久婷婷成人综合色麻豆| 国产精品一区二区三区四区久久 | 国产精品久久久久成人av| 女人高潮潮喷娇喘18禁视频| 一级,二级,三级黄色视频| 欧美另类亚洲清纯唯美| 97人妻天天添夜夜摸| 一进一出好大好爽视频| 精品免费久久久久久久清纯| 欧美一级毛片孕妇| 欧美日韩精品网址| 国产91精品成人一区二区三区| 一区二区三区精品91| 美女扒开内裤让男人捅视频| 侵犯人妻中文字幕一二三四区| 在线永久观看黄色视频| 久久午夜亚洲精品久久| 无人区码免费观看不卡| 9191精品国产免费久久| 好男人电影高清在线观看| 亚洲精品国产精品久久久不卡| 免费在线观看黄色视频的| 久久久久久久久免费视频了| 女性被躁到高潮视频| 热re99久久国产66热| 亚洲视频免费观看视频| 精品国产乱子伦一区二区三区| 99在线人妻在线中文字幕| 国产成人免费无遮挡视频| 午夜福利,免费看| 两人在一起打扑克的视频| 99久久久亚洲精品蜜臀av| 免费人成视频x8x8入口观看| 丝袜美足系列| 国产欧美日韩一区二区三| 久久久国产成人精品二区 | 亚洲一区二区三区欧美精品| 正在播放国产对白刺激| 亚洲av美国av| 亚洲欧洲精品一区二区精品久久久| 国产三级在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 天天躁夜夜躁狠狠躁躁| 久久香蕉激情| 神马国产精品三级电影在线观看 | 午夜免费鲁丝| 亚洲视频免费观看视频| 中文欧美无线码| 午夜视频精品福利| 中国美女看黄片| 无限看片的www在线观看| 久久久久久久精品吃奶| 久久精品亚洲熟妇少妇任你| 18禁美女被吸乳视频| 757午夜福利合集在线观看| 国产精品av久久久久免费| 高潮久久久久久久久久久不卡| 在线观看免费视频网站a站| 成人国语在线视频| 嫩草影视91久久| 一区二区三区精品91| 老汉色∧v一级毛片| 在线观看午夜福利视频| 亚洲成人久久性| 精品无人区乱码1区二区| 在线视频色国产色| www国产在线视频色| 成人av一区二区三区在线看| 在线观看免费日韩欧美大片| 久久香蕉精品热| 亚洲七黄色美女视频| 手机成人av网站| 91av网站免费观看| 久久精品国产综合久久久| 久久久久久久午夜电影 | 亚洲 欧美一区二区三区| 日韩欧美一区视频在线观看| 很黄的视频免费| 亚洲欧美激情在线| 精品乱码久久久久久99久播| 欧美黄色淫秽网站| 又黄又爽又免费观看的视频| 国产精品一区二区三区四区久久 | 热99re8久久精品国产| 亚洲中文字幕日韩| 国产三级黄色录像| 国产精品永久免费网站| 老司机亚洲免费影院| 天天躁狠狠躁夜夜躁狠狠躁| 成人手机av| 色婷婷久久久亚洲欧美| 亚洲中文av在线| 日本一区二区免费在线视频| 国产精品综合久久久久久久免费 | 欧美乱妇无乱码| 国产精品98久久久久久宅男小说| 精品国产乱码久久久久久男人| 国产成人精品久久二区二区免费| 淫秽高清视频在线观看| 免费观看精品视频网站| 俄罗斯特黄特色一大片| 亚洲午夜理论影院| 亚洲狠狠婷婷综合久久图片| 亚洲国产精品合色在线| 一夜夜www| 女人爽到高潮嗷嗷叫在线视频| av网站免费在线观看视频| 成人特级黄色片久久久久久久| 人人妻人人澡人人看| 国产又爽黄色视频| 国产精品一区二区免费欧美| 国产精品日韩av在线免费观看 | 亚洲自拍偷在线| 欧美中文综合在线视频| 成年人免费黄色播放视频| 亚洲专区字幕在线| 视频在线观看一区二区三区| 好男人电影高清在线观看| 一区在线观看完整版| 妹子高潮喷水视频| 天堂中文最新版在线下载| bbb黄色大片| 国产亚洲精品一区二区www| 亚洲成人免费av在线播放| 国产精品 国内视频| 亚洲精品国产色婷婷电影| 亚洲色图av天堂| 日韩欧美三级三区| 三级毛片av免费| 欧美精品啪啪一区二区三区| 99热只有精品国产| 亚洲人成伊人成综合网2020| 亚洲av五月六月丁香网| 香蕉丝袜av| 国产精品 欧美亚洲| av欧美777| 欧美日韩一级在线毛片| 少妇 在线观看| 亚洲欧美一区二区三区黑人| 欧美久久黑人一区二区| 国产主播在线观看一区二区| 三级毛片av免费| 91字幕亚洲| 女警被强在线播放| 又黄又爽又免费观看的视频| 日韩av在线大香蕉| 丁香六月欧美| 国产成年人精品一区二区 | 自线自在国产av| 一级片'在线观看视频| 天堂动漫精品| 丁香六月欧美| 啪啪无遮挡十八禁网站| 一边摸一边做爽爽视频免费| 操美女的视频在线观看| 99久久久亚洲精品蜜臀av| 精品午夜福利视频在线观看一区| 国产成年人精品一区二区 | 日韩人妻精品一区2区三区| 日本免费一区二区三区高清不卡 | 国产区一区二久久| 少妇粗大呻吟视频| 女警被强在线播放| 国产亚洲精品第一综合不卡| 99久久综合精品五月天人人| 身体一侧抽搐| 成人亚洲精品av一区二区 | 搡老熟女国产l中国老女人| 成在线人永久免费视频| 天堂俺去俺来也www色官网| e午夜精品久久久久久久| 亚洲欧美日韩无卡精品| 国产精品一区二区免费欧美| 欧美日韩乱码在线| 日韩精品中文字幕看吧| 99在线视频只有这里精品首页| 丝袜在线中文字幕| 亚洲avbb在线观看| 国产精品国产av在线观看| 久久精品国产综合久久久| 午夜免费激情av| 99久久久亚洲精品蜜臀av| 国产精品久久久av美女十八| 美女扒开内裤让男人捅视频| 久久精品人人爽人人爽视色| 黑人巨大精品欧美一区二区mp4| 97超级碰碰碰精品色视频在线观看| a级片在线免费高清观看视频| 9191精品国产免费久久| 国产精品偷伦视频观看了| 成人三级黄色视频| 国产精品永久免费网站| 99精国产麻豆久久婷婷| 老司机亚洲免费影院| 在线观看免费午夜福利视频| 国产精品香港三级国产av潘金莲| 久久影院123| 男女之事视频高清在线观看| 欧美日韩瑟瑟在线播放| 日本一区二区免费在线视频| 国产免费av片在线观看野外av| 在线国产一区二区在线| 黄片大片在线免费观看| 久久人人精品亚洲av| 美女扒开内裤让男人捅视频| 一个人免费在线观看的高清视频| 超碰成人久久| 国产欧美日韩一区二区精品| 好男人电影高清在线观看| 天堂影院成人在线观看| 久久国产精品男人的天堂亚洲| x7x7x7水蜜桃| 99精品欧美一区二区三区四区| 国产激情欧美一区二区| 亚洲av电影在线进入| 久久精品91蜜桃| 午夜老司机福利片| 精品国产美女av久久久久小说| 黄色a级毛片大全视频| 久久天堂一区二区三区四区| 亚洲一区二区三区不卡视频| 亚洲aⅴ乱码一区二区在线播放 | 国产精品香港三级国产av潘金莲| 国产av一区在线观看免费| 亚洲性夜色夜夜综合| 亚洲自偷自拍图片 自拍| 操出白浆在线播放| 美女大奶头视频| 人人澡人人妻人| 天天添夜夜摸| 超碰97精品在线观看| 亚洲avbb在线观看| 日韩视频一区二区在线观看| 美女国产高潮福利片在线看| 老司机午夜十八禁免费视频| 午夜福利,免费看| 久久青草综合色| 中国美女看黄片| 女人被躁到高潮嗷嗷叫费观| 国产熟女午夜一区二区三区| 女人高潮潮喷娇喘18禁视频| svipshipincom国产片| 亚洲精品美女久久久久99蜜臀| 麻豆一二三区av精品|