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

    上行與下行效應(yīng)對浮游動物的長期影響評價
    ——以滇池與撫仙湖沉積物象鼻溞(Bosmina)為例*

    2015-05-06 07:15:11盧慧斌陳光杰陳小林劉曉東顧兆炎張虎才趙帥營段立曾張自強
    湖泊科學(xué) 2015年1期
    關(guān)鍵詞:銀魚撫仙湖滇池

    盧慧斌,陳光杰,陳小林,劉曉東,顧兆炎,張虎才,趙帥營,陳 麗,段立曾,張自強

    (1:云南師范大學(xué)旅游與地理科學(xué)學(xué)院高原湖泊生態(tài)與全球變化重點實驗室,昆明 650500)(2:中國科技大學(xué)地球與空間學(xué)院極地環(huán)境研究室,合肥 230026)(3:中國科學(xué)院地質(zhì)與地球物理研究所,北京 100029)

    上行與下行效應(yīng)對浮游動物的長期影響評價
    ——以滇池與撫仙湖沉積物象鼻溞(Bosmina)為例*

    盧慧斌1,陳光杰1**,陳小林1,劉曉東2,顧兆炎3,張虎才1,趙帥營1,陳 麗1,段立曾1,張自強1

    (1:云南師范大學(xué)旅游與地理科學(xué)學(xué)院高原湖泊生態(tài)與全球變化重點實驗室,昆明 650500)(2:中國科技大學(xué)地球與空間學(xué)院極地環(huán)境研究室,合肥 230026)(3:中國科學(xué)院地質(zhì)與地球物理研究所,北京 100029)

    近年來云南高原湖泊面臨富營養(yǎng)化、漁業(yè)活動增強等多重環(huán)境壓力的疊加影響,對湖泊的有效治理與生態(tài)修復(fù)急需對多重壓力下生態(tài)系統(tǒng)的響應(yīng)模式進行系統(tǒng)了解.現(xiàn)有研究表明在系統(tǒng)生產(chǎn)力和捕食壓力的不同配置下,湖泊系統(tǒng)主要組成(如浮游動物)的響應(yīng)特征可能出現(xiàn)差異且捕食壓力可能隨營養(yǎng)水平的變化而改變,目前對云南湖泊生態(tài)系統(tǒng)的研究主要集中于單一環(huán)境壓力下的生態(tài)響應(yīng).本研究以目前分別處于重富營養(yǎng)和中-貧營養(yǎng)水平的滇池和撫仙湖為研究對象,應(yīng)用湖泊沉積物記錄進行多指標(biāo)分析,探討受外來魚類影響下兩個大型湖泊浮游動物長期響應(yīng)模式的異同.通過象鼻溞生物量與個體大小等指標(biāo),重建了近百年來滇池與撫仙湖典型浮游動物的變化歷史,結(jié)果表明隨著湖泊生產(chǎn)力水平(如沉積物色素生產(chǎn)量)的增加,2個湖泊中象鼻溞生物量顯著增加,同時物種相對組成出現(xiàn)明顯變化(如Bosminalongispina被B.longirostris取代),指示湖泊上行效應(yīng)對浮游動物的控制作用.同時象鼻溞的生物量、殼長與觸角長度的變化在1960s與1980s有明顯降低的趨勢,與同期外來魚類(如銀魚)引入與漁業(yè)產(chǎn)量增加的時間一致.進一步應(yīng)用多變量回歸分析與方差分解方法來定量評價上行與下行效應(yīng)對象鼻溞生物量變化的驅(qū)動強度,結(jié)果表明富營養(yǎng)化(沉積物色素)和捕食作用(象鼻溞殼長)對浮游動物長期變化的驅(qū)動強度比較相似(分別解釋了生物量變化的77.25%和83.59%),然而在滇池下行效應(yīng)對象鼻溞生物量的獨立影響比在貧營養(yǎng)的撫仙湖中更強(分別為15.46%和10.39%),上行效應(yīng)對象鼻溞生物量的影響在撫仙湖要明顯強于滇池(分別為69.74%和19.67%),而在滇池上行與下行效應(yīng)的相互作用強度明顯強于撫仙湖(分別為42.12%和3.46%).結(jié)果表明隨著湖泊營養(yǎng)水平的升高,浮游動物的生物量在2個湖泊中均顯著增加,而外來魚類(如太湖新銀魚)的引入和經(jīng)濟魚類數(shù)量的增加加劇了對浮游動物的捕食壓力,造成了浮游動物的生物量降低和個體減小.但在快速富營養(yǎng)化的滇池,對浮游動物的捕食壓力隨營養(yǎng)水平的變化出現(xiàn)較強的依賴性,而在總體處于中-貧營養(yǎng)水平的撫仙湖中此相互作用較弱.結(jié)果表明在不同營養(yǎng)水平的大型湖泊中,營養(yǎng)水平的變化幅度可以導(dǎo)致魚類捕食壓力的差異性變化,指示了對魚類捕食壓力的評價和浮游動物長期變化的特征分析需要考慮湖泊的營養(yǎng)水平與富營養(yǎng)化過程的差異.總之,富營養(yǎng)化和外來魚類的引入導(dǎo)致了高原湖泊生態(tài)系統(tǒng)的快速響應(yīng)與結(jié)構(gòu)變化,因此對高原湖泊的生態(tài)修復(fù)需要考慮湖泊營養(yǎng)水平對生態(tài)系統(tǒng)結(jié)構(gòu)與食物鏈作用的影響.

    古湖沼學(xué);象鼻溞;生物量;個體大?。簧闲行?yīng);下行效應(yīng);滇池;撫仙湖

    隨著人類活動的加劇和社會經(jīng)濟的持續(xù)發(fā)展,云南典型湖泊的生態(tài)系統(tǒng)結(jié)構(gòu)與功能在過去幾十年里發(fā)生了較大變化.例如滇池藍(lán)藻水華頻繁暴發(fā)[1],撫仙湖土著魚類大量減少[2],陽宗海水體砷污染等[3].各種環(huán)境壓力已經(jīng)對云南湖泊的生態(tài)系統(tǒng)產(chǎn)生了重要的影響,尤其是湖泊中的水生生物,而這些環(huán)境壓力的長期生態(tài)效應(yīng)可以被湖泊沉積物記錄下來[4].浮游動物處于湖泊食物鏈的中間環(huán)節(jié),既直接攝食浮游植物,其本身又是很多魚類的重要食物來源.同時浮游動物對湖泊環(huán)境的變化十分敏感,已成為評價湖泊生態(tài)環(huán)境變化的重要指標(biāo)之一.目前對云南湖泊生態(tài)系統(tǒng)中部分重要的水生生物指標(biāo)(如浮游動物)缺乏長期、連續(xù)的監(jiān)測數(shù)據(jù)和較為系統(tǒng)的研究,尤其是對多重環(huán)境壓力驅(qū)動湖泊生態(tài)系統(tǒng)的疊加和獨立效應(yīng)缺乏系統(tǒng)評價,因此急需開展對浮游動物長期變化的歷史重建與影響因子識別.

    枝角類是甲殼綱浮游動物的重要組成部分,是反映湖泊營養(yǎng)水平與捕食壓力變化的敏感生物指標(biāo)之一.象鼻溞(Bosmina)是云南湖泊中一種常見的浮游枝角類,其生物量與物種組成通常受湖泊初級生產(chǎn)力的控制,如B.longispina主要出現(xiàn)在較低的營養(yǎng)水平時期,而富營養(yǎng)時期則被B.longirostris所代替[5];同時由于枝角類作為一種重要的魚類食物對魚類捕食壓力的變化十分敏感,其個體大小已成為捕食壓力變化的良好指示器[6].大量研究表明大型食浮游動物魚類趨向于選擇性捕食個體更大、更可見的屬種(例如溞屬)和種群內(nèi)相對更大的個體,而一些個體相對較小的種類(例如象鼻溞屬)更容易受到小型食浮游動物魚類和大型無脊椎動物捕食者的攻擊和選擇性捕食[7].

    由于對云南湖泊水生生物的監(jiān)測存在時間尺度較短且缺乏連續(xù)、系統(tǒng)的調(diào)查等問題,使得在應(yīng)用現(xiàn)代監(jiān)測數(shù)據(jù)了解湖泊生態(tài)系統(tǒng)變化方面受到較大的限制.目前云南湖泊的污染修復(fù)與生態(tài)治理對獲取長期生態(tài)環(huán)境數(shù)據(jù)的需要變得急為迫切,古湖沼學(xué)的研究方法與技術(shù)手段可以為研究湖泊生態(tài)系統(tǒng)的長期變化提供重要信息[8].例如,湖泊食物網(wǎng)的一些重要組成部分(如硅質(zhì)藻類、幾丁質(zhì)無脊椎動物)可以較好地在湖泊沉積物中保存,使得它們成為古湖沼學(xué)研究中的重要生物指標(biāo),它們在沉積物中的生物量、群落組成等信息可以用來重建湖泊生物群落變化的歷史,并為評價不同環(huán)境壓力的長期生態(tài)效應(yīng)提供基礎(chǔ)數(shù)據(jù).在湖泊沉積物中象鼻溞的數(shù)量較大、保存良好,其響應(yīng)捕食壓力所表現(xiàn)的形態(tài)變化特征(殼長、刺長與第一觸角長度)在湖泊沉積物中可以較為完整地保存,同時根據(jù)其頭孔位置等特征可以鑒定到物種的水平,使得象鼻溞在古湖沼學(xué)的應(yīng)用中成為指示湖泊環(huán)境壓力變化的良好指示器[9].

    在云南湖泊的研究中,一些學(xué)者已經(jīng)利用湖泊沉積物中的枝角類提取出了湖泊環(huán)境變化的信息.例如,Liu等利用湖泊沉積物中枝角類群落結(jié)構(gòu)的變化探討了滇池中浮游動物對湖泊富營養(yǎng)化、外來魚類引入的響應(yīng)特征[10],也利用湖泊沉積物中的枝角類群落結(jié)構(gòu)變化以及象鼻溞的形態(tài)變化探討了撫仙湖富營養(yǎng)化和魚類引入對浮游動物的長期影響[11].然而,由于象鼻溞生物量的變化既受到湖泊富營養(yǎng)化又受到魚類捕食的影響,目前我們對上行效應(yīng)與下行效應(yīng)對不同湖泊的象鼻溞生物量變化是否存在差異性影響,以及這兩種效應(yīng)的相互作用強度如何還知之甚少.大量湖泊生態(tài)實驗與野外數(shù)據(jù)表明魚類對浮游動物的捕食壓力會隨著湖泊營養(yǎng)水平的變化而改變.如Jeppesen等發(fā)現(xiàn)在湖泊營養(yǎng)水平較低或較高時魚類的捕食壓力較大,而在中等的營養(yǎng)水平下魚類捕食壓力相對較弱.因此,在不同營養(yǎng)水平下湖泊浮游動物面臨的魚類捕食壓力隨著湖泊生產(chǎn)力的變化可能出現(xiàn)不同的變化模式,如在重富營養(yǎng)化的水體中魚類捕食壓力與營養(yǎng)水平的相互作用可能要強于營養(yǎng)水平總體較低且變化幅度較小的湖泊[12].

    本研究以目前處于重富營養(yǎng)水平的滇池和中-貧營養(yǎng)水平的撫仙湖為研究對象,通過分析沉積物色素、象鼻溞等生物記錄,識別象鼻溞生物量隨著湖泊營養(yǎng)水平變化的響應(yīng)模式,并對比2個湖泊上行與下行效應(yīng)對象鼻溞生物量的相對驅(qū)動強度及其相互作用強度.因此,通過對比不同營養(yǎng)水平及經(jīng)歷不同富營養(yǎng)化過程的湖泊沉積物記錄,可以幫助我們探討湖泊營養(yǎng)水平與捕食壓力對浮游動物變化的長期驅(qū)動過程與機制,為認(rèn)識多重環(huán)境壓力下湖泊生態(tài)系統(tǒng)的響應(yīng)模式并對污染湖泊開展生態(tài)修復(fù)提供重要的科學(xué)依據(jù).

    1 材料與方法

    1.1 樣品采集

    位于滇中地區(qū)的滇池為大型淺水湖泊,目前處于重富營養(yǎng)水平;而撫仙湖則是大型深水湖泊,目前處于中-貧營養(yǎng)水平(表1).利用重力采樣器于2012年6月采集滇池(24°46′7″N,102°39′51″E)樣品,水深5.4m,柱長33cm,以0.5cm間隔進行分樣;于2012年3月采集撫仙湖(24°28′10″N,102°52′12″E)樣品,水深58m,柱長21cm,按5cm以前間隔0.25cm,5cm以后間隔0.5cm進行分樣.樣品密封并避光保存,豎立帶回實驗室.

    表1 滇池與撫仙湖的湖泊學(xué)基本特征*

    *水域面積、容積、平均水深和最大水深數(shù)據(jù)來源于文獻(xiàn)[13];滇池水質(zhì)數(shù)據(jù)為外海資料,除標(biāo)注外其它數(shù)據(jù)均為本次調(diào)查結(jié)果.

    1.2 樣品處理與數(shù)據(jù)分析

    通過對樣品放射性核素210Pb和137Cs進行測試分析,利用美國Canberra公司生產(chǎn)的高純鍺探測器、數(shù)字化譜儀及多通道分析系統(tǒng),建立了兩個鉆孔的年代學(xué)序列.具體步驟包括:將凍干的樣品研磨,稱取3~5g左右干重樣品放入與標(biāo)準(zhǔn)源同一規(guī)格的塑料容器內(nèi)并蠟封2周左右,使226Ra與210Pb處于衰變平衡體系,然后測定不同深度樣品的210Pb和137Cs放射性強度,并應(yīng)用恒定放射性通量模式建立年代學(xué)系列[14].

    圖1 象鼻溞的長度測量Fig.1 Diagram showing the length measurement of Bosmina

    沉積物枝角類樣品的處理參照Korhola等的方法[15].首先將凍干的樣品充分混合均勻,稱取約0.1g(干重)樣品,倒入100ml燒杯中;再加入50ml的10%KOH,并在電熱板上加熱到70℃左右約30min,期間利用玻璃棒輕輕地攪拌樣品,以使樣品充分地分解絮凝;然后用78μm的篩網(wǎng)過濾,再用蒸餾水將篩網(wǎng)上的樣品清洗到帶有刻度的15ml離心管中,加入3滴95%的乙醇用于樣品的室溫保存,并定容.再將樣品充分搖勻,利用0.05ml的移液槍吸取樣品溶液,置于載玻片上,在電熱板上以50℃加熱,在樣品蒸發(fā)待干前,加入兩滴甘油-藏紅染料,用針輕輕攪拌使其混合并均勻分布,將蓋玻片置于其上,待其冷卻后用指甲油密封蓋玻片的邊緣以便長期保存.象鼻溞的鑒定與計數(shù)在萊卡顯微鏡(Leica DM2500)下進行,依據(jù)Goulden和Frey的鑒定標(biāo)準(zhǔn)[16],象鼻溞有Bosminalongispina與B.longirostris2個種.同時統(tǒng)計頭殼、觸角、后腹部的個數(shù)并使總個體數(shù)量不少于60個,最后換算成絕對豐度(每克干重沉積物中象鼻溞個體數(shù)量)作為生物量的指標(biāo).并測量玻片中所見象鼻溞的第一觸角長度(LA)、殼長(LC)與刺長(LM)(圖1)[17-18],每個樣品中象鼻溞不同部位的測量樣本個數(shù)不低于10.

    藻類植物死亡后產(chǎn)生的某些色素長期保存在湖泊沉積物中,由于不同植物具有不同的色素組合類型,因此根據(jù)湖泊沉積物中色素的含量和種類可恢復(fù)湖泊的初級生產(chǎn)力.沉積物色素分析方法主要參考文獻(xiàn)[19].具體步驟包括:將沉積物樣品攪拌均勻并凍干;稱重50~150mg,不同色素指標(biāo)按照標(biāo)準(zhǔn)方法用0.2mm孔隙先過濾,然后在最小光輻照下用氮氣干燥;類胡蘿卜素、葉綠素和衍生化合物含量用高效液相色譜儀(HPLC)測量,并用蘇丹染料做參照標(biāo)準(zhǔn).

    為探討上行效應(yīng)與下行效應(yīng)對象鼻溞生物量的影響強度,本文應(yīng)用沉積物Chl.a色素含量這一指標(biāo)指示湖泊初級生產(chǎn)力,用象鼻溞的殼長(LC)代表魚類的捕食壓力,對象鼻溞生物量的變化進行多元回歸分析,探討在滇池與撫仙湖中兩個自變量(Chl.a和LC)對因變量(象鼻溞生物量)變化的獨立解釋程度和相互作用強度,并應(yīng)用方差分解方法對結(jié)果進行圖形表達(dá).在統(tǒng)計分析中相關(guān)數(shù)據(jù)首先進行正態(tài)化分布的檢驗并進行相應(yīng)的數(shù)據(jù)轉(zhuǎn)換,本文所采用的統(tǒng)計分析在Rv 2.3.0軟件中完成,沉積物地層圖應(yīng)用C2 v1.4.3軟件進行處理分析.

    2 結(jié)果

    2.1 滇池和撫仙湖沉積物色素與象鼻溞變化的歷史重建

    近百年來沉積物色素記錄表明湖泊初級生產(chǎn)量在滇池和撫仙湖中均出現(xiàn)不同程度的增加,但其變化時間序列和幅度呈現(xiàn)明顯不同的特征. 在20世紀(jì)中期之前2個湖泊沉積物色素記錄持續(xù)處于較低的含量且變化幅度較??;1960s初期開始,滇池的Chl.a含量開始顯著增加并呈現(xiàn)持續(xù)快速增加的特征(圖2a),而撫仙湖從1950s開始出現(xiàn)小幅度的葉綠素含量上升,但總體增加幅度相對較小(圖2b).

    沉積物象鼻溞密度在2個湖泊中也出現(xiàn)較為顯著的變化.在湖泊生產(chǎn)力持續(xù)較低的時期,2個湖泊的象鼻溞密度總體維持在較低的水平,且以喜好貧營養(yǎng)生境的Bosminalongispina為主(如占全部象鼻溞密度的百分比大于60%~70%).隨著湖泊營養(yǎng)水平的升高,B.longispina在象鼻溞密度中所占的比例逐漸減少,而喜好富營養(yǎng)的B.longirostris占比則逐漸升高(圖2).同時2個湖泊中象鼻溞的密度總體上也隨湖泊生產(chǎn)力的增加而上升(圖3),但在撫仙湖其密度在過去40年里出現(xiàn)總體下降的趨勢(圖2b).

    與沉積物色素和象鼻溞密度的變化相比,沉積物象鼻溞的個體大小呈現(xiàn)較大的波動.滇池和撫仙湖中,象鼻溞總體大小的殼長呈現(xiàn)與較大個體物種B.longispina第一觸角長度相似的變化趨勢(滇池:r=0.578,P=0.01,n=18).與外來銀魚總捕撈量的歷史數(shù)據(jù)對比發(fā)現(xiàn),象鼻溞總體大小的殼長在銀魚引入時期出現(xiàn)明顯的減小.如在1980s初期開始2個湖泊中皆引入銀魚并很快成為建群種(圖2),同期象鼻溞的殼長和B.longispina第一觸角長度均出現(xiàn)明顯減小,指示了銀魚的引入導(dǎo)致對包括象鼻溞在內(nèi)的浮游動物捕食壓力的增加,從而對象鼻溞群落和個體大小產(chǎn)生抑制效應(yīng).

    圖2 滇池(a)和撫仙湖(b)沉積物指標(biāo)變化的地層圖Fig.2 Stratigraphic plots showing sedimentary pigment and bosminid changes at Lake Dianchi(a) and Lake Fuxian(b)

    2.2 上行與下行效應(yīng)對象鼻溞生物量長期變化的驅(qū)動分析

    為了評價上行與下行效應(yīng)對湖泊生態(tài)系統(tǒng)的驅(qū)動作用,應(yīng)用沉積物色素指標(biāo)指示湖泊初級生產(chǎn)力水平(即上行效應(yīng)),同時以沉積物象鼻溞殼長代表魚類捕食壓力的變化(即下行效應(yīng)),通過應(yīng)用多元線性回歸方法來定量探討兩種效應(yīng)對象鼻溞生物量變化的相對驅(qū)動強度及相互作用強度.方差分析結(jié)果表明,Chl.a、象鼻溞的LC共同解釋了滇池象鼻溞生物量變化的77.25%(P<0.001,n=18),Chl.a獨立解釋了象鼻溞生物量變化的19.67%,LC獨立解釋了象鼻溞生物量變化的15.46%,Chl.a與LC的相互作用可以解釋浮游動物生物量變化的42.12%.而在撫仙湖兩個自變量共同解釋了象鼻溞生物量變化的83.59%(P<0.001,n=15),Chl.a獨立解釋了象鼻溞生物量變化的69.74%,LC單一解釋了象鼻溞生物量變化的10.39%,Chl.a與LC的相互作用僅為3.46%.

    3 討論

    圖3 滇池和撫仙湖沉積物象鼻溞與Chl.a生產(chǎn)量的關(guān)系Fig.3 The relationship between the sedimentary bosminid and Chl.a productions at Lake Dianchi and Lake Fuxian

    多指標(biāo)沉積物記錄揭示了滇池和撫仙湖近百年來富營養(yǎng)化的歷史與浮游動物長期響應(yīng)的過程.20世紀(jì)中后期出現(xiàn)的富營養(yǎng)化導(dǎo)致了2個湖泊沉積物記錄的湖泊初級生產(chǎn)力水平上升(圖2),象鼻溞在湖泊沉積物中記錄的生物量也相應(yīng)上升(圖3),且都呈現(xiàn)喜貧營養(yǎng)水平的物種B.longispina被B.longirostris逐漸取代的趨勢,這些變化與富營養(yǎng)化之前持續(xù)較低的浮游植物與浮游動物生物量明顯不同.與現(xiàn)代記錄相一致的是沉積物記錄表明滇池經(jīng)歷了持續(xù)、快速富營養(yǎng)化過程,并與撫仙湖總體較弱的富營養(yǎng)化過程形成明顯差別.在對象鼻溞生物量變化的長期驅(qū)動強度上,上行效應(yīng)對滇池浮游動物的影響比撫仙湖弱,而下行效應(yīng)則比撫仙湖強.上行效應(yīng)與下行效應(yīng)的相互作用在富營養(yǎng)的滇池中較為顯著,表明在營養(yǎng)水平較高、水體較淺的滇池上行效應(yīng)影響了魚類(如引入的銀魚)對浮游動物的捕食壓力變化,而在中-貧營養(yǎng)的撫仙湖則表現(xiàn)為上行效應(yīng)較強的控制作用.

    3.1 滇池與撫仙湖近百年來浮游動物長期變化的歷史與特征

    20世紀(jì)中期以來滇池和撫仙湖沉積物記錄的Chl.a水平出現(xiàn)不同程度的增加,與現(xiàn)代監(jiān)測數(shù)據(jù)記錄的2個湖泊近30年的富營養(yǎng)化趨勢基本一致[20-21].隨著2個湖泊初級生產(chǎn)力水平的升高,B.longispina在2個湖泊中都呈現(xiàn)被B.longirostris逐漸取代的趨勢.我們對滇池沉積物中象鼻溞生物量分析的結(jié)果與Liu等的結(jié)果一致[10],主要是從1960s開始顯著增加的.我們的沉積物色素數(shù)據(jù)表明,撫仙湖湖泊初級生產(chǎn)力的提高,使得浮游動物的資源可獲得性增強,從而象鼻溞的生物量也開始明顯增加.例如Liu等在對撫仙湖的研究中發(fā)現(xiàn)營養(yǎng)鹽增加后,B.longirostris的絕對豐度明顯增加[11].

    在1960s與1980s期間,滇池與撫仙湖象鼻溞的生物量與殼長以及第一觸角長度出現(xiàn)顯著的波動,分別對應(yīng)了“四大家魚”、太湖新銀魚等外來魚類的引入.在撫仙湖與滇池,外來魚類的引入促進了湖泊中漁業(yè)捕撈量的增加,但同時對生態(tài)系統(tǒng)的影響之一是可能增加湖泊中的魚類捕食壓力.關(guān)于象鼻溞的形態(tài)變化與魚類捕食關(guān)系的研究在國外已取得大量成果,主要結(jié)果表明魚類的選擇性捕食會造成象鼻溞個體大小的減小[22-23],而象鼻溞刺長的變化比較復(fù)雜,可能更多的是響應(yīng)無脊椎動物捕食者的捕食壓力[24-27].例如在挪威北部湖泊中,一種高效率的食浮游動物魚類歐白鮭的引入,導(dǎo)致B.longispina和B.longirostris產(chǎn)卵母體的大小下降.由于B.longispina的個體大,使得它比B.longirostris更易受到攻擊,所以其個體大小呈現(xiàn)顯著的減小.雖然B.longirostris被認(rèn)為是不易受到選擇性捕食的魚類的攻擊,但是它平均個體大小的減小也被觀察到,表明它也受到歐白鮭的選擇性捕食[28].本研究發(fā)現(xiàn)象鼻溞生物量、殼長以及第一觸角長度的下降,時間上對應(yīng)于外來魚類的引入,尤其是象鼻溞的殼長的變化主要是由魚類捕食壓力的增加導(dǎo)致的,所以本文用其指示湖泊中魚類捕食壓力變化是合理的.Liu等在對滇池沉積物枝角類的研究中也發(fā)現(xiàn),在1980s中期到1990s末這一階段,幾乎所有枝角類屬種的絕對豐度都呈現(xiàn)下降的趨勢,而且也發(fā)現(xiàn)Daphinia休眠卵的長度變短[10].Liu等的這些研究結(jié)果反映了這一時期魚類捕食壓力的變化對浮游動物的明顯影響,但由于該研究沒有提取象鼻溞個體大小等直接反映魚類捕食壓力的信息,我們的結(jié)果是對該研究結(jié)論的一個重要補充.同時從Liu等的數(shù)據(jù)中可以看出,Bosminaspp.與B.longirostris的絕對豐度數(shù)據(jù)在1960s也出現(xiàn)顯著的下降,可能也間接反映了“四大家魚”引入這一事件.我們的分析結(jié)果與Liu等在對撫仙湖沉積物枝角類的研究中得到的1987年以后象鼻溞長度變化的信息一致[11].這些研究總體都表明象鼻溞在云南部分湖泊中是捕食壓力變化的良好指示者,特別是其個體大小對湖泊中魚類捕食壓力的變化十分敏感.

    表2 滇池枝角類與橈足類密度的現(xiàn)代調(diào)查數(shù)據(jù)[29,31-32]

    Tab.2 Modern survey data of cladoceran and copepod densities at Lake Dianchi

    年份橈足類/(ind./L)枝角類/(ind./L)19571355519805706021982165168198445.568199492420

    滇池等部分高原湖泊的現(xiàn)代調(diào)查數(shù)據(jù)也表明,湖泊浮游動物密度與魚類數(shù)量的變化呈現(xiàn)很好的一致性.1980s初期太湖新銀魚引入滇池后增殖較為成功,其產(chǎn)量迅速在湖泊中形成規(guī)模后對湖泊中的浮游動物密度有較大影響.尤其是在1984年太湖新銀魚的產(chǎn)量達(dá)到3600t,與此同時,水體中的橈足類、枝角類僅有45.5、68ind./L,比1980年分別下降了92%與88.7%(表2)[29].因此在太湖新銀魚捕食壓力增強的背景下,枝角類與橈足類生物量受到了捕食壓力的抑制作用而呈現(xiàn)協(xié)同下降的特征.關(guān)于太湖新銀魚食性的研究已有許多成果,主要研究在對浮游動物的選擇性捕食上都表明象鼻溞是太湖新銀魚的主要食物來源之一.高禮存等報道長刺溞(Daphnialongispina)在滇池銀魚(如平均體長為43~76mm)食物中數(shù)量最多,其次依次為象鼻溞、劍水蚤(Cyclops)、鏢水蚤(Diaptomus)、盤腸溞(Chydorus)和米蝦幼體[29].劉恩生等對太湖新銀魚的食物組成分析中發(fā)現(xiàn),太湖新銀魚的主要食物是橈足類,其次是枝角類,而輪蟲幾乎不被攝食[30].本文基于沉積物記錄的結(jié)果表明,象鼻溞在長時間尺度上表現(xiàn)為對太湖新銀魚引入十分敏感,與對該外來魚類的食性分析相一致.

    太湖新銀魚在撫仙湖中的增殖和種群數(shù)量與滇池表現(xiàn)出截然不同的模式,可能反映2個湖泊營養(yǎng)水平的差異導(dǎo)致魚類捕食壓力對浮游動物群落的差異性影響.滇池中的太湖新銀魚經(jīng)過1980s的增殖繁盛期后種群數(shù)量明顯下降,而在撫仙湖中卻逐漸成為優(yōu)勢種群(圖2).太湖新銀魚自1980s中期引入撫仙湖以來,成功增殖并形成穩(wěn)定產(chǎn)量,已成為撫仙湖漁業(yè)捕撈中的主要經(jīng)濟魚類.與此同時,撫仙湖中的太湖新銀魚在湖泊中增殖的成功對水體中的象鼻溞存在長期強烈的捕食壓力,使得象鼻溞的生物量、殼長與B.longispina的LA穩(wěn)定在低值水平(圖2b).而滇池中,可能由于營養(yǎng)水平的快速增加,使得湖泊中的魚類捕食壓力隨營養(yǎng)水平出現(xiàn)明顯的變化,同時滇池較強的漁業(yè)捕撈、魚類養(yǎng)殖等人類活動也可能影響浮游動物對魚類捕食的響應(yīng)模式.

    3.2 上行與下行效應(yīng)對浮游動物長期變化的驅(qū)動強度分析

    長期以來大量研究表明典型湖泊的食物網(wǎng)主要是由上行效應(yīng)控制,即資源可獲得性的驅(qū)動作用.例如浮游植物被營養(yǎng)鹽和光所控制,浮游動物被浮游植物控制,食浮游動物魚類被浮游動物所控制等[33].隨后的研究表明,食物網(wǎng)結(jié)構(gòu)可能同時受到下行效應(yīng)的控制,例如食浮游動物魚類控制浮游動物、浮游動物控制浮游植物等[34].然而,對食物網(wǎng)中上行效應(yīng)與下行效應(yīng)的相對影響程度研究還主要基于現(xiàn)代湖泊調(diào)查和中小型實驗的數(shù)據(jù).例如McQueen等的研究認(rèn)為,上行效應(yīng)在食物網(wǎng)底部是最強的,而且越往上影響越??;下行效應(yīng)在食物網(wǎng)的頂部是最強的,越往下影響越小[33,35].Jeppesen等在對從溫帶到北極地區(qū)的466個湖泊的下行效應(yīng)研究中發(fā)現(xiàn),在營養(yǎng)水平高的湖泊中,食浮游動物的魚類通過作用于浮游動物,進而對浮游植物產(chǎn)生顯著影響,而這一影響在營養(yǎng)水平低的湖泊中卻不明顯[12].本文方差分析的結(jié)果也發(fā)現(xiàn),在富營養(yǎng)的滇池中,代表魚類捕食壓力的殼長與代表湖泊初級生產(chǎn)力的Chl.a相互作用強度十分顯著,而這一相互作用在撫仙湖中較弱.所以,本研究結(jié)果拓展了Jeppesen等的研究區(qū)域,表明上行效應(yīng)與下行效應(yīng)的相互作用模式不僅在溫帶到北極地區(qū)的湖泊中存在,在低緯高原區(qū)的湖泊中也十分明顯.

    Jeppesen等還發(fā)現(xiàn)在空間尺度上對大個體浮游動物的捕食控制在淺水湖泊中比在深水湖泊中更高[12].我們的結(jié)果表明在大型重富營養(yǎng)的淺水湖泊滇池,下行效應(yīng)對象鼻溞生物量變化的影響比在大型貧營養(yǎng)的深水撫仙湖中大;上行效應(yīng)對象鼻溞生物量變化的影響在貧營養(yǎng)的深水撫仙湖中更大.Jeppesen等探討了淺水湖泊中下行效應(yīng)更強的原因,認(rèn)為可能是淺水湖泊中單位體積水柱中的食浮游動物魚類密度更大,或者底棲的無脊椎動物維持了食浮游動物魚類的生物量,進而對浮游動物造成高的捕食壓力,或者淺水湖泊對于浮游動物來講缺少避難所[36].

    Liu等在對撫仙湖的研究中提出,在貧營養(yǎng)湖泊中資源的可利用性可能對枝角類種群的發(fā)展更為重要[11].盡管太湖新銀魚的引入能通過糞便排放、尸體分解等過程,通過食物鏈的傳遞作用增加水體中營養(yǎng)鹽的釋放,增加藻類的初級生產(chǎn)力,然而對湖濱區(qū)面積狹小的深水撫仙湖中魚類的這種影響可能是非常微弱的.而在滇池這種湖濱區(qū)面積較大的大型淺水湖泊中,魚類通過此類作用促進湖濱區(qū)與底泥中營養(yǎng)鹽的釋放,進而促進藻類的生長,可能是比較顯著的.因此,在湖濱區(qū)面積不同與深度不同的湖泊中,魚類對湖泊中營養(yǎng)鹽循環(huán)所起的作用是不同的.

    4 結(jié)論

    通過對滇池與撫仙湖沉積物多指標(biāo)記錄的對比研究,本文在重建近百年來兩個湖泊初級生產(chǎn)力和象鼻溞長期變化歷史的基礎(chǔ)上,定量評價了上行與下行效應(yīng)對兩個大型湖泊浮游動物長期變化的驅(qū)動強度,并識別了營養(yǎng)水平等湖泊類型的差異可能對湖泊生態(tài)系統(tǒng)過程產(chǎn)生的差異性影響.具體來說,滇池與撫仙湖象鼻溞的生物量隨著湖泊營養(yǎng)水平和初級生產(chǎn)力的升高出現(xiàn)增加的趨勢,同時象鼻溞生物量以及個體大小變化特征與兩個湖泊外來魚類的引入時間基本一致.進一步的分析表明,下行效應(yīng)對象鼻溞生物量的影響在滇池中更大,而上行效應(yīng)則在撫仙湖中更強,上行效應(yīng)與下行效應(yīng)的相互作用在富營養(yǎng)的滇池中更顯著.總之,本結(jié)果表明沉積物象鼻溞記錄的生物量與個體大小變化能夠記錄湖泊富營養(yǎng)化與魚類捕食壓力的變化信息,是湖泊生態(tài)環(huán)境變化的良好指示者.因此,對不同類型的湖泊開展以浮游動物為代表的生態(tài)系統(tǒng)重要組分長期變化特征研究并進行對比分析,可以為我們認(rèn)識湖泊生態(tài)系統(tǒng)結(jié)構(gòu)與功能的長期變化模式與驅(qū)動機制提供重要的科學(xué)依據(jù),可為對污染湖泊開展生態(tài)修復(fù)并進行有效的湖泊管理提供科學(xué)依據(jù).

    致謝:感謝評審人對本文提出寶貴的意見,同時對在本研究中給予野外采樣、室內(nèi)分析等幫助的云南師范大學(xué)高原湖泊生態(tài)與全球變化重點實驗室的各位老師和同學(xué)表示感謝.

    [1] 郭懷成,王心宇,伊 璇.基于滇池水生態(tài)系統(tǒng)演替的富營養(yǎng)化控制策略.地理研究,2013,32(6):998-1006.

    [2] 熊 飛,李文朝,潘繼征等.云南撫仙湖魚類資源現(xiàn)狀與變化.湖泊科學(xué),2006,18(3):305-311.

    [3] 王振華,何 濱,潘學(xué)軍等.云南陽宗海砷污染水平、變化趨勢及風(fēng)險評估.中國科學(xué):化學(xué),2011,41(3):556-564.

    [4] Smol JP. Pollution of lakes and rivers: a paleoenvironmental perspective. New York:John Wiley & Sons, 2008.

    [5] Boucherle MM, Züllig H. Cladoceran remains as evidence of change in trophic state in three Swiss lakes.Hydrobiology, 1983,103:141-146.

    [6] Dodson SI. Adaptive change in plankton morphology in response to size-selective predation: a new hypothesis of cyclomorphosis.LimnologyandOceanography, 1974,19(5):721-729.

    [7] Brooks JL, Dodson SI. Predation, body size, and composition of plankton.Science, 1965,150(3692):28-35.

    [8] Smol JP. The power of the past: using sediments to track the effects of multiple stressors on lake ecosystems.FreshwaterBiology, 2010,55(s1):43-59.

    [9] Korosi JB, Kurek J, Smol JP. A review on utilizingBosminasize structure archived in lake sediments to infer historic shifts in predation regimes.JournalofPlanktonResearch, 2013,35(2):444-460.

    [10] Liu G, Liu Z, Chen Fetal. Response of the cladoceran community to eutrophication, fish introductions and degradation of the macrophyte vegetation in Lake Dianchi, a large, shallow plateau lake in southwestern China.Limnology, 2013,14(2):159-166.

    [11] Liu G, Liu Z, Li Yetal. Effects of fish introduction and eutrophication on the cladoceran community in Lake Fuxian, a deep oligotrophic lake in southwest China.JournalofPaleolimnology, 2009,42(3):427-435.

    [12] Jeppesen E, Jensen JP, Jensen Cetal. The impact of nutrient state and lake depth on top-down control in the pelagic zone of lakes: a study of 466 lakes from the temperate zone to the Arctic.Ecosystems, 2003,6(4):313-325.

    [13] 王蘇民,竇鴻身.中國湖泊志.北京:科學(xué)出版社,1998:368-377.

    [14] 王小雷,楊 浩,丁兆運等.云南撫仙湖近現(xiàn)代沉積速率變化研究.地理學(xué)報,2011,66(11):1551-1561.

    [15] Korhola A, Rautio M. Cladocera and other branchiopod crustaceans, tracking environmental change using lake sediments. Dordrecht:Springer Netherlands, 2001:5-41.

    [16] Goulden CE, Frey DG. The occurrence and significance of lateral head pores in the genusBosmina(Cladocera).Hydrobiology, 1963,48(3):513-522.

    [17] Sweetman JN, Finney BP. Differential responses of zooplankton populations (Bosminalongirostris) to fish predation and nutrient-loading in an introduced and a natural sockeye salmon nursery lake on Kodiak Island, Alaska, USA.JournalofPaleolimnology, 2003,30(2):183-193.

    [18] Ahlen E, Bystrom P, Korsman Tetal. Relationships between planktivore community capacity (PCC) and cladoceran microfossils in northern Swedish lakes.FundamentalandAppliedLimnology, 2011,178(4):315-324.

    [19] Leavitt PR, Hodgson DA. Sedimentary pigments, Tracking environmental change using lake sediments. Dordrecht:Springer Netherlands, 2001:295-325.

    [20] 李中杰,鄭一新,張大為等.滇池流域近20年社會經(jīng)濟發(fā)展對水環(huán)境的影響.湖泊科學(xué),2012,24(6):875-882.

    [21] 高 偉,陳 巖,徐 敏等.撫仙湖水質(zhì)變化(1980-2011年)趨勢與驅(qū)動力分析.湖泊科學(xué),2013,25(5):635-642.

    [22] Kerfoot WC. The divergence of adjacent populations.Ecology, 1975,56(6):1298-1313.

    [23] Stenson JAE. Significance of predator influence on composition ofBosminaspp. populations.LimnologyandOceanography, 1976,21(6):814-822.

    [24] Black RW. The nature and causes of cyclomorphosis in a species of theBosminalongirostriscomplex.Ecology, 1980,61(5):1122-1132.

    [25] Sprules WG, Carter JCH, Ramcharan CW. Phenotypic associations in the Bosminidae (Cladocera):zoogeographic patterns.LimnologyandOceanography, 1984,29(1):161-169.

    [26] de Stasio BT, Howard E, Hanson Retal. Cyclomorphosis inEubosminaspecies in New England lakes, USA.VerhandlungenderInternationalenVereinigungfürTheoretischeundAngewandteLimnologie,1990,24(1):481-486.

    [27] Post DM, Frost TM, Kitchell JF. Morphological responses byBosminalongirostrisandEubosminatubicento changes in copepod predator populations during a whole-lake acidification experiment.JournalofPlanktonResearch, 1995,17(8):1621-1632.

    [28] Amundsen PERA, Siwertsson A, Primicerio Retal. Long-term responses of zooplankton to invasion by a planktivorous fish in a subarctic watercourse.FreshwaterBiology, 2009,54(1):24-34.

    [29] 高禮存,莊大棟,遲金釗等.太湖短吻銀魚移殖滇池試驗研究.湖泊科學(xué),1989,1(1):79-87.

    [30] 劉恩生,鮑傳和,吳林坤等.太湖新銀魚的食物組成及與主要魚類間關(guān)系的初步分析.水利漁業(yè),2006,26(5):41-44.

    [31] 王忠澤.滇池浮游動物的初步調(diào)查.云南大學(xué)學(xué)報:自然科學(xué)版,1985,7:53-72.

    [32] 彭瓊英.滇池浮游動物的調(diào)查分析.水利漁業(yè),1995,(6):22-24.

    [33] McQueen DJ, Johannes MRS, Post JRetal. Bottom-up and top-down impacts on freshwater pelagic community structure.EcologicalMonographs, 1989,59(3):289-309.

    [34] Carpenter S, Kitchell J, Hodgson J. Cascading trophic interactions and lake productivity.BioScience, 1985,35(10):634-639.

    [35] McQueen DJ, Post JR, Mills EL. Trophic relationships in freshwater pelagic ecosystems.CanadianJournalofFisheriesandAquaticSciences, 1986,43(8):1571-1581.

    [36] Jeppesen E, Jensen JP, Sndergaard Metal. Top-down control in freshwater lakes: the role of nutrient state, submerged macrophytes and water depth.Hydrobiology, 1997,342/343:151-164.

    The long-term effects of bottom-up and top-down forcing on zooplankton: An example from sedimentary bosminid records of Lake Dianchi and Lake Fuxian

    LU Huibin1, CHEN Guangjie1, CHEN Xiaolin1, LIU Xiaodong2, GU Zhaoyan3, ZHANG Hucai1, ZHAO Shuaiying1, CHEN Li1, DUAN Lizeng1& ZHANG Ziqiang1

    (1:KeyLaboratoryofPlateauLakeEcologyandGlobalChange,SchoolofTourismandGeography,YunnanNormalUniversity,Kunming650500,P.R.China)(2:InstituteofPolarEnvironment,UniversityofScienceandTechnologyofChina,Hefei230026,P.R.China)(3:InstituteofGeologyandGeophysics,ChineseAcademyofSciences,Beijing100029,P.R.China)

    In recent years, lakes of Yunnan Plateau have increasingly experienced the pressure of multiple environmental stressors like eutrophication and fish invasion, as well as their synergic effects. To date, most studies have focused on the ecological impacts of a single stressor, in spite of the interaction of these multiple stressors. This paper selected Lake Dianchi and Lake Fuxian, which are currently hypo-eutrophic and oligo-mesotrophic respectively, and conducted multi-proxy analyses of lake sediments from both lakes. We aim to quantify the impacts of both eutrophication (bottom-up control) and predation pressure (top-down control) in driving the long-term changes of zooplankton production, and to compare the intensity of their single forcing and interaction between both lakes. Our results showed that both lakes experienced increased primary production (i.e. the production of sedimentary pigments) since the mid-20thcentury, which was consistent with the eutrophication history revealed by lake monitoring data. Concurrently, the production of zooplankton (i.e. sedimentary bosminid production) increases significantly with a general trend of the oligotrophic taxon ofBosminalongispinareplaced byB.longirostrisfound in both lakes. However, the bosminid production and body size (i.e. carapace and antennule length), displayed a decrease during the 1960s and 1980s, when exotic fish were introduced and/or established stable populations. With the methods of multiple regression and variance partitioning, our further analyses suggested that top-down control (as indicated by bosminid carapace length) exerted stronger impact on zooplankton production at Lake Dianchi than at Lake Fuxian (15.46% and 10.39%, respectively), while Lake Fuxian displayed a much stronger control of bottom-up (represented by sedimentary pigment concentration) forcing than Dianchi (69.74% and 19.67%, respectively). In addition, we found that the interaction between bottom-up and top-down is much more significant at Lake Dianchi than the Lake Fuxian (42.12 % and 3.46 %, respectively), despite that both forcing exerted a similar degree of control on zooplankton production at Lake Dianchi (77.25%) and Lake Fuxian (83.59%). These results indicate the significant control of lake eutrophication on zooplankton production at both lakes, while the introduction of exotic fish intensify the predation pressure on zooplankton, which also showed strong dependence on the rapid increase of trophic status at the hypo-eutrophic shallow lake of Lake Dianchi. In contrast, the strength of the interaction between bottom-up and top-down control was insignificant at the oligo-mesotrophic deep lake of Fuxian. Over all, our results infered that eutrophication and fish predation have led to the significant changes of ecosystem structure and trophic interaction at both large plateau lakes. Therefore, trophic status, the history of eutrophication and lake typology should be taken into account for sustainable management and ecological restoration of the plateau lakes in Yunnan.

    Paleolimnology;Bosmina; biomass; body size; bottom-up control; top-down control; Lake Dianchi; Lake Fuxian

    *國家自然科學(xué)基金項目(U1133601,41171048)、教育部科學(xué)技術(shù)研究項目(213034A)、云南省中青年學(xué)術(shù)技術(shù)帶頭人培養(yǎng)計劃項目(2012HB012)、云南省創(chuàng)新團隊項目(2011C1)和云南省高端人才項目(2010CI111)聯(lián)合資助. 2014-01-09收稿;2014-05-12收修改稿. 盧慧斌(1990~),男,碩士研究生;E-mail: luhuibin611@163.com.

    **通信作者;E-mail:guangjiechen@gmail.com.

    猜你喜歡
    銀魚撫仙湖滇池
    撫仙湖,萬頃平湖一鑒清
    金秋(2021年12期)2021-10-06 04:07:28
    浪漫抵達(dá)之處——讀《銀魚來》
    《銀魚來》創(chuàng)作談
    銀魚來
    在滇池東岸看西山
    文苑(2020年8期)2020-09-09 09:30:56
    撫仙湖戀歌
    青年歌聲(2019年12期)2019-12-17 06:32:34
    銀魚和螞蟻
    走進撫仙湖
    民族音樂(2018年1期)2018-04-18 03:24:28
    云南省撫仙湖流域水權(quán)制度建設(shè)研究
    滇池濕地,美回來了
    永久网站在线| 久久国产亚洲av麻豆专区| 女人精品久久久久毛片| 精品国产乱码久久久久久小说| 午夜av观看不卡| 99热国产这里只有精品6| av福利片在线| 少妇人妻久久综合中文| 亚洲成国产人片在线观看| 大码成人一级视频| 一本—道久久a久久精品蜜桃钙片| 自拍欧美九色日韩亚洲蝌蚪91| 午夜久久久在线观看| 天堂8中文在线网| 国产综合精华液| 卡戴珊不雅视频在线播放| 婷婷色av中文字幕| 岛国毛片在线播放| 最新中文字幕久久久久| 国精品久久久久久国模美| 亚洲欧美日韩另类电影网站| 熟女av电影| 日韩人妻精品一区2区三区| xxx大片免费视频| av电影中文网址| 久久久久久人妻| 免费看光身美女| 男女啪啪激烈高潮av片| 青春草亚洲视频在线观看| 一二三四在线观看免费中文在 | 免费女性裸体啪啪无遮挡网站| av福利片在线| 免费观看无遮挡的男女| 久久精品国产a三级三级三级| 女人久久www免费人成看片| 22中文网久久字幕| 国产国拍精品亚洲av在线观看| av不卡在线播放| 国产色爽女视频免费观看| 99精国产麻豆久久婷婷| 国产一区亚洲一区在线观看| 人妻少妇偷人精品九色| 国产精品不卡视频一区二区| 国产黄频视频在线观看| av在线播放精品| 婷婷成人精品国产| 精品酒店卫生间| 久久精品国产自在天天线| 看免费av毛片| 亚洲伊人色综图| 久久女婷五月综合色啪小说| 我的女老师完整版在线观看| 9热在线视频观看99| 美女中出高潮动态图| 国产成人精品久久久久久| 久久这里只有精品19| av国产精品久久久久影院| 婷婷色av中文字幕| 国产精品一区二区在线观看99| 又大又黄又爽视频免费| 人妻人人澡人人爽人人| 亚洲欧美色中文字幕在线| 你懂的网址亚洲精品在线观看| 亚洲精品乱久久久久久| 亚洲欧美日韩另类电影网站| 久久久久国产网址| 久久久久久伊人网av| 好男人视频免费观看在线| 老熟女久久久| 国产成人精品久久久久久| 午夜激情久久久久久久| 麻豆精品久久久久久蜜桃| 18禁国产床啪视频网站| 日韩伦理黄色片| 卡戴珊不雅视频在线播放| 高清毛片免费看| 亚洲人成77777在线视频| 七月丁香在线播放| 一本大道久久a久久精品| 亚洲第一av免费看| 婷婷色麻豆天堂久久| 中国国产av一级| 另类精品久久| 美女脱内裤让男人舔精品视频| 99久久中文字幕三级久久日本| 亚洲综合精品二区| 人妻 亚洲 视频| 日产精品乱码卡一卡2卡三| 国产日韩一区二区三区精品不卡| 18禁观看日本| 免费av不卡在线播放| av片东京热男人的天堂| av女优亚洲男人天堂| 日本欧美视频一区| 亚洲精品久久久久久婷婷小说| 精品久久国产蜜桃| 韩国av在线不卡| 午夜福利乱码中文字幕| 久久久久久久久久久久大奶| 黄片播放在线免费| 如何舔出高潮| 性高湖久久久久久久久免费观看| 亚洲国产av新网站| av不卡在线播放| 涩涩av久久男人的天堂| 亚洲成人av在线免费| 精品一区二区免费观看| 免费日韩欧美在线观看| 另类亚洲欧美激情| 日本色播在线视频| 好男人视频免费观看在线| av黄色大香蕉| 欧美日韩一区二区视频在线观看视频在线| 国产成人精品婷婷| 最近2019中文字幕mv第一页| 国产精品国产av在线观看| 亚洲激情五月婷婷啪啪| 综合色丁香网| 高清在线视频一区二区三区| 最近中文字幕2019免费版| 久久久久久伊人网av| 春色校园在线视频观看| 午夜福利乱码中文字幕| 国产黄频视频在线观看| 国产免费一区二区三区四区乱码| 久久精品国产a三级三级三级| av视频免费观看在线观看| 青春草视频在线免费观看| 熟女av电影| 老司机影院毛片| 人成视频在线观看免费观看| 蜜臀久久99精品久久宅男| 国产一区二区在线观看av| 久久久精品94久久精品| 人妻人人澡人人爽人人| 亚洲精品中文字幕在线视频| 美女中出高潮动态图| 精品国产一区二区三区久久久樱花| 超色免费av| 中文精品一卡2卡3卡4更新| av卡一久久| 久热久热在线精品观看| 成人午夜精彩视频在线观看| 国产免费一级a男人的天堂| 亚洲欧美成人综合另类久久久| 久久 成人 亚洲| 涩涩av久久男人的天堂| 黄片无遮挡物在线观看| 在现免费观看毛片| 看免费成人av毛片| 99热网站在线观看| 亚洲五月色婷婷综合| 亚洲精品美女久久av网站| 少妇人妻 视频| 日本av手机在线免费观看| 五月玫瑰六月丁香| 亚洲精品aⅴ在线观看| 日本与韩国留学比较| 国产精品国产三级国产av玫瑰| 久久精品国产综合久久久 | 国产成人午夜福利电影在线观看| 精品福利永久在线观看| 又大又黄又爽视频免费| 精品卡一卡二卡四卡免费| 久久久久久人人人人人| 高清视频免费观看一区二区| 制服丝袜香蕉在线| 亚洲综合色网址| 国产精品久久久久久久电影| 亚洲,欧美,日韩| 亚洲激情五月婷婷啪啪| 久久ye,这里只有精品| 99热这里只有是精品在线观看| 欧美+日韩+精品| 哪个播放器可以免费观看大片| 国产成人av激情在线播放| 一区二区三区乱码不卡18| 十八禁高潮呻吟视频| 一边亲一边摸免费视频| 夜夜骑夜夜射夜夜干| 乱人伦中国视频| 少妇高潮的动态图| 91午夜精品亚洲一区二区三区| 一级毛片 在线播放| 人妻人人澡人人爽人人| 一个人免费看片子| 男女下面插进去视频免费观看 | 午夜精品国产一区二区电影| 国产在线视频一区二区| 成人黄色视频免费在线看| 亚洲色图综合在线观看| 亚洲精品中文字幕在线视频| 97超碰精品成人国产| av在线观看视频网站免费| 天堂8中文在线网| 亚洲五月色婷婷综合| 亚洲,欧美,日韩| 一区二区av电影网| 国产免费一区二区三区四区乱码| 最近的中文字幕免费完整| 午夜福利乱码中文字幕| freevideosex欧美| 国产亚洲最大av| 人妻少妇偷人精品九色| 国产成人午夜福利电影在线观看| 久久人妻熟女aⅴ| 国产成人免费观看mmmm| av电影中文网址| 美女中出高潮动态图| 国产在线免费精品| 亚洲国产av新网站| 春色校园在线视频观看| 日本91视频免费播放| 女性生殖器流出的白浆| 少妇人妻 视频| 黄色怎么调成土黄色| 国产免费一区二区三区四区乱码| 如何舔出高潮| 在线精品无人区一区二区三| 91aial.com中文字幕在线观看| 一级a做视频免费观看| 亚洲国产色片| 成年动漫av网址| 亚洲av欧美aⅴ国产| 国产精品免费大片| 亚洲国产精品成人久久小说| 精品国产一区二区三区四区第35| 亚洲欧美精品自产自拍| 大片电影免费在线观看免费| 美女国产高潮福利片在线看| 精品一品国产午夜福利视频| 日韩人妻精品一区2区三区| av天堂久久9| 欧美亚洲日本最大视频资源| 国产在线免费精品| 桃花免费在线播放| 在现免费观看毛片| 成年女人在线观看亚洲视频| 亚洲国产精品专区欧美| 纯流量卡能插随身wifi吗| 亚洲av欧美aⅴ国产| 9191精品国产免费久久| 亚洲精品美女久久av网站| 九色成人免费人妻av| 久久久久久人妻| 18在线观看网站| 夜夜骑夜夜射夜夜干| 久久鲁丝午夜福利片| 一本色道久久久久久精品综合| 久久久久网色| av在线老鸭窝| 永久网站在线| 国产国语露脸激情在线看| 成人国语在线视频| 啦啦啦中文免费视频观看日本| 人妻系列 视频| 中国美白少妇内射xxxbb| 欧美日韩视频精品一区| 9191精品国产免费久久| 男女无遮挡免费网站观看| 性高湖久久久久久久久免费观看| 99久久中文字幕三级久久日本| 亚洲欧美中文字幕日韩二区| 精品人妻在线不人妻| 天天操日日干夜夜撸| 亚洲精品国产色婷婷电影| av女优亚洲男人天堂| 一区二区日韩欧美中文字幕 | 桃花免费在线播放| 一本久久精品| 亚洲丝袜综合中文字幕| 国产熟女午夜一区二区三区| 欧美国产精品va在线观看不卡| 精品一区二区三卡| 大香蕉久久成人网| 午夜免费鲁丝| 18禁观看日本| 90打野战视频偷拍视频| 精品久久久久久电影网| 99久国产av精品国产电影| 久久狼人影院| 欧美xxⅹ黑人| 伦理电影免费视频| 999精品在线视频| 亚洲欧美成人精品一区二区| 一级片免费观看大全| 一本色道久久久久久精品综合| 精品久久久精品久久久| 晚上一个人看的免费电影| 人人妻人人爽人人添夜夜欢视频| 精品国产乱码久久久久久小说| 久久韩国三级中文字幕| 欧美 日韩 精品 国产| 女人久久www免费人成看片| 美女国产视频在线观看| 免费日韩欧美在线观看| kizo精华| 免费在线观看完整版高清| 寂寞人妻少妇视频99o| 亚洲伊人久久精品综合| 99久久精品国产国产毛片| 成年女人在线观看亚洲视频| 免费观看在线日韩| 少妇人妻 视频| 日本欧美国产在线视频| 国产在线一区二区三区精| 亚洲人成77777在线视频| 国产国拍精品亚洲av在线观看| 国产乱来视频区| 国产在线视频一区二区| 美女国产视频在线观看| 下体分泌物呈黄色| 久久鲁丝午夜福利片| 久久影院123| 男人添女人高潮全过程视频| 久久人人爽人人爽人人片va| 国产xxxxx性猛交| 久久久久久伊人网av| 国产爽快片一区二区三区| 免费人妻精品一区二区三区视频| 亚洲欧美日韩另类电影网站| 视频在线观看一区二区三区| 最近的中文字幕免费完整| 满18在线观看网站| 男人爽女人下面视频在线观看| 99国产综合亚洲精品| 毛片一级片免费看久久久久| 性高湖久久久久久久久免费观看| 毛片一级片免费看久久久久| 免费大片黄手机在线观看| 在线观看美女被高潮喷水网站| av视频免费观看在线观看| 国产男女内射视频| 久久精品aⅴ一区二区三区四区 | 亚洲精品乱久久久久久| 在线观看一区二区三区激情| 中文精品一卡2卡3卡4更新| 午夜福利视频在线观看免费| 亚洲一级一片aⅴ在线观看| 日日啪夜夜爽| 成人国产麻豆网| 成人二区视频| 观看av在线不卡| 热re99久久精品国产66热6| 久久免费观看电影| 狂野欧美激情性bbbbbb| 男女啪啪激烈高潮av片| 最近的中文字幕免费完整| 22中文网久久字幕| 亚洲精品成人av观看孕妇| 男女啪啪激烈高潮av片| 麻豆乱淫一区二区| 久久99蜜桃精品久久| 少妇 在线观看| 国产福利在线免费观看视频| 日韩免费高清中文字幕av| 国产精品一区二区在线观看99| 亚洲欧美日韩卡通动漫| 国产一区二区激情短视频 | 国产欧美日韩一区二区三区在线| 成人黄色视频免费在线看| 另类精品久久| 国产一区二区在线观看日韩| 亚洲美女搞黄在线观看| 日韩制服骚丝袜av| 精品亚洲成a人片在线观看| 久久这里只有精品19| 亚洲美女搞黄在线观看| 亚洲av.av天堂| 亚洲av在线观看美女高潮| 亚洲欧洲国产日韩| 久久精品国产a三级三级三级| 亚洲国产av新网站| 久久久a久久爽久久v久久| 国产精品国产三级国产av玫瑰| 涩涩av久久男人的天堂| 最近最新中文字幕免费大全7| 亚洲国产精品国产精品| 最近最新中文字幕免费大全7| 老女人水多毛片| 成人国语在线视频| 精品视频人人做人人爽| 久久ye,这里只有精品| 国产精品国产av在线观看| 中文字幕最新亚洲高清| 国产探花极品一区二区| 卡戴珊不雅视频在线播放| 亚洲成色77777| 精品卡一卡二卡四卡免费| videosex国产| 一级毛片我不卡| 国产av码专区亚洲av| 亚洲综合精品二区| 天天操日日干夜夜撸| 国产精品99久久99久久久不卡 | 欧美日韩一区二区视频在线观看视频在线| www日本在线高清视频| 亚洲一码二码三码区别大吗| 亚洲五月色婷婷综合| 免费高清在线观看视频在线观看| 黄色 视频免费看| 免费观看a级毛片全部| 亚洲国产欧美日韩在线播放| 精品一区二区三区四区五区乱码 | av免费观看日本| 如何舔出高潮| 99热网站在线观看| 青春草视频在线免费观看| 男女边吃奶边做爰视频| 国产伦理片在线播放av一区| 两个人看的免费小视频| 在线天堂中文资源库| 毛片一级片免费看久久久久| 免费女性裸体啪啪无遮挡网站| 国产老妇伦熟女老妇高清| 男女下面插进去视频免费观看 | 看非洲黑人一级黄片| 99久久人妻综合| 亚洲国产精品专区欧美| 亚洲国产看品久久| 亚洲av电影在线观看一区二区三区| 高清视频免费观看一区二区| 国产av国产精品国产| av有码第一页| 99久久精品国产国产毛片| 人成视频在线观看免费观看| 青春草国产在线视频| 国产精品无大码| 午夜精品国产一区二区电影| 女人精品久久久久毛片| 青春草视频在线免费观看| 伦理电影免费视频| 黄色配什么色好看| 亚洲国产最新在线播放| 寂寞人妻少妇视频99o| 亚洲人与动物交配视频| 哪个播放器可以免费观看大片| 欧美成人午夜免费资源| 日韩三级伦理在线观看| 国产不卡av网站在线观看| 亚洲情色 制服丝袜| 免费大片黄手机在线观看| videossex国产| 日本黄大片高清| 国产黄色视频一区二区在线观看| 中文欧美无线码| 国产爽快片一区二区三区| 成人国产av品久久久| 王馨瑶露胸无遮挡在线观看| 国产成人精品无人区| 久久精品夜色国产| 中文字幕人妻熟女乱码| 国产亚洲精品第一综合不卡 | 91精品三级在线观看| 91精品伊人久久大香线蕉| 一级毛片黄色毛片免费观看视频| 女的被弄到高潮叫床怎么办| 五月开心婷婷网| av电影中文网址| 伊人亚洲综合成人网| 欧美丝袜亚洲另类| 免费播放大片免费观看视频在线观看| 国产色爽女视频免费观看| 国产亚洲午夜精品一区二区久久| 欧美激情国产日韩精品一区| av卡一久久| 寂寞人妻少妇视频99o| 亚洲精品成人av观看孕妇| 美女主播在线视频| 久久久a久久爽久久v久久| 黑人猛操日本美女一级片| 99re6热这里在线精品视频| 91精品伊人久久大香线蕉| 亚洲欧美精品自产自拍| 建设人人有责人人尽责人人享有的| 免费观看a级毛片全部| av在线app专区| 国产精品三级大全| 久久鲁丝午夜福利片| 夜夜骑夜夜射夜夜干| 1024视频免费在线观看| 男人舔女人的私密视频| 大码成人一级视频| 高清视频免费观看一区二区| 最后的刺客免费高清国语| 天天躁夜夜躁狠狠久久av| 免费看av在线观看网站| 成人国产av品久久久| 欧美 日韩 精品 国产| 一级毛片 在线播放| 国产av一区二区精品久久| 青春草亚洲视频在线观看| 成人亚洲精品一区在线观看| 成年人午夜在线观看视频| 黄片无遮挡物在线观看| 在线观看人妻少妇| 成人18禁高潮啪啪吃奶动态图| 久久久国产精品麻豆| 男人爽女人下面视频在线观看| 日本与韩国留学比较| 国产乱人偷精品视频| 夫妻午夜视频| 女性生殖器流出的白浆| 国产 精品1| av免费观看日本| 亚洲欧美色中文字幕在线| 精品第一国产精品| 在线天堂最新版资源| 少妇 在线观看| 亚洲欧洲国产日韩| 人妻 亚洲 视频| 国产成人免费无遮挡视频| 男女下面插进去视频免费观看 | 一本色道久久久久久精品综合| 亚洲欧美成人精品一区二区| 丰满少妇做爰视频| 91午夜精品亚洲一区二区三区| 天天躁夜夜躁狠狠躁躁| 国产精品三级大全| 久久久国产精品麻豆| 国产精品99久久99久久久不卡 | 香蕉国产在线看| 乱码一卡2卡4卡精品| 久久午夜综合久久蜜桃| 国产白丝娇喘喷水9色精品| 97人妻天天添夜夜摸| 久久精品国产自在天天线| 精品视频人人做人人爽| 国产1区2区3区精品| 男的添女的下面高潮视频| 中文字幕制服av| 国产亚洲欧美精品永久| 在线天堂中文资源库| 国产片特级美女逼逼视频| 日本午夜av视频| 天天影视国产精品| 午夜免费鲁丝| 国产福利在线免费观看视频| 18禁动态无遮挡网站| 高清av免费在线| 五月开心婷婷网| 九色亚洲精品在线播放| 久久这里只有精品19| 国产成人午夜福利电影在线观看| 日韩欧美精品免费久久| av播播在线观看一区| 丝袜喷水一区| 国产成人av激情在线播放| 制服丝袜香蕉在线| 久久久久久久大尺度免费视频| 国产国拍精品亚洲av在线观看| 国产在线免费精品| 51国产日韩欧美| 亚洲伊人久久精品综合| 亚洲,欧美精品.| 人妻人人澡人人爽人人| 中文欧美无线码| 巨乳人妻的诱惑在线观看| 国精品久久久久久国模美| 美女大奶头黄色视频| 高清毛片免费看| 成人影院久久| 精品国产露脸久久av麻豆| 超碰97精品在线观看| 精品午夜福利在线看| 在线天堂中文资源库| 黄色视频在线播放观看不卡| 极品人妻少妇av视频| 国产伦理片在线播放av一区| 99久国产av精品国产电影| 99香蕉大伊视频| 久久这里有精品视频免费| 18禁在线无遮挡免费观看视频| 这个男人来自地球电影免费观看 | 欧美成人午夜免费资源| 久久久久人妻精品一区果冻| 精品熟女少妇av免费看| 国产精品蜜桃在线观看| 97在线视频观看| 人体艺术视频欧美日本| 妹子高潮喷水视频| 亚洲内射少妇av| 老司机影院毛片| 欧美3d第一页| 午夜激情久久久久久久| 国产白丝娇喘喷水9色精品| 人妻系列 视频| 亚洲av福利一区| 色94色欧美一区二区| 少妇熟女欧美另类| xxx大片免费视频| 成人免费观看视频高清| 99热这里只有是精品在线观看| 777米奇影视久久| 久久久国产精品麻豆| 韩国高清视频一区二区三区| av又黄又爽大尺度在线免费看| 国产欧美亚洲国产| 国产 精品1| 十分钟在线观看高清视频www| 十八禁网站网址无遮挡| 久久久久国产精品人妻一区二区| 18禁裸乳无遮挡动漫免费视频| 黄色配什么色好看| 欧美成人午夜精品| 自拍欧美九色日韩亚洲蝌蚪91| 日韩熟女老妇一区二区性免费视频| a级毛色黄片| 大码成人一级视频| 国产精品一区二区在线观看99| 欧美3d第一页| √禁漫天堂资源中文www| 美女脱内裤让男人舔精品视频| 国产成人欧美| 亚洲欧洲精品一区二区精品久久久 |