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

    集裝箱養(yǎng)殖生態(tài)凈化池塘中浮游植物群落結(jié)構(gòu)多樣性分析

    2024-01-01 00:00:00閔文武王龍燕陳飛雄周其椿
    水產(chǎn)科技情報 2024年5期
    關(guān)鍵詞:浮游植物物種多樣性群落結(jié)構(gòu)

    摘" 要:為了解集裝箱養(yǎng)殖排放水污染狀況及其生態(tài)凈化池塘中浮游植物群落結(jié)構(gòu)的變化動態(tài),于2020年對“集裝箱+生態(tài)池塘”循環(huán)水養(yǎng)殖系統(tǒng)的生態(tài)凈化池塘水體進(jìn)行了調(diào)查采樣,運(yùn)用物種多樣性指數(shù)分析浮游植物群落結(jié)構(gòu)特征,運(yùn)用密度和生物量比較曲線(ABC曲線)和聚類分析探究浮游植物群落的受干擾程度及其時間聚類差異,運(yùn)用冗余分析揭示浮游植物群落與環(huán)境因子的響應(yīng)關(guān)系。結(jié)果顯示,調(diào)查期間生態(tài)凈化池塘水體中共檢出浮游植物7門49屬80種,以綠藻門為主,占比為47.5%,其次是硅藻門,占17.5%;各月份浮游植物物種數(shù)的平均值為21種,其中7月份物種數(shù)最多(26種),10月份最少(13種)。賴氏澤絲藻(Limnothrix redekei)、湖泊假魚腥藻(Pseudanabaena limnetica)、尖針桿藻(Synedra acus)和叢球韋斯藻(Westella botryoides)為優(yōu)勢種。浮游植物密度、生物量、Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)和Margalef豐富度指數(shù)平均值分別為13.31×106 cells/L、7.16 mg/L、1.93、0.64和1.26。ABC曲線顯示,浮游植物群落結(jié)構(gòu)受中等干擾,群落組成主要以密度和生物量占優(yōu)勢的物種為主。聚類分析結(jié)果顯示,浮游植物群落呈現(xiàn)不顯著的季節(jié)分異特征。冗余分析結(jié)果顯示,水體氨氮含量與浮游植物群落結(jié)構(gòu)顯著相關(guān)。

    關(guān)鍵詞:水產(chǎn)養(yǎng)殖;生態(tài)池塘;浮游植物;群落結(jié)構(gòu);物種多樣性

    doi:10.16446/j.fsti.20230200122

    收稿日期:2023-02-22

    作者簡介:閔文武(1987—),男,助理研究員,主要從事水產(chǎn)養(yǎng)殖技術(shù)及水環(huán)境研究。E-mail:minwenwu@126.com

    項(xiàng)目資助:貴州省科技計劃項(xiàng)目[黔科合支撐(2020)1Y102;(2020)1Y105];貴州省農(nóng)業(yè)科學(xué)院青年科技基金項(xiàng)目[黔農(nóng)科院青年基金" (2019)25];貴州省科技計劃項(xiàng)目[黔科合重大專項(xiàng)字(2019)3010];貴陽市科技計劃項(xiàng)目[筑科合同(2022)3-29號]。

    浮游植物作為水生態(tài)系統(tǒng)中的主要初級生產(chǎn)者,對維持水生態(tài)系統(tǒng)健康穩(wěn)定及其完整性具有重要作用

    [1-2]。浮游植物生長周期短,對環(huán)境變化敏感,其群落結(jié)構(gòu)特征的變化對水環(huán)境健康狀況具有較好的指示作用[3-4]。通過研究浮游植物群落結(jié)構(gòu)特征來評價水生態(tài)環(huán)境的健康狀態(tài)及生物完整性,在河流、湖泊、養(yǎng)殖水域等水體已得到大量的運(yùn)用,為指導(dǎo)生態(tài)環(huán)境管理和提高水產(chǎn)養(yǎng)殖生態(tài)效益提供了大量參考依據(jù)[5-9]?!凹b箱+生態(tài)池塘”養(yǎng)殖技術(shù)具有養(yǎng)殖和水體污染物處理雙高效的特點(diǎn),是農(nóng)業(yè)農(nóng)村部2018—2020年主推的養(yǎng)殖模式之一[10]。有關(guān)該養(yǎng)殖技術(shù)及水體污染物處理工藝等方面的研究不斷涌現(xiàn),為推動該模式的日趨完善提供了大量參考依據(jù)[11-14]。在有關(guān)水產(chǎn)養(yǎng)殖污染物處理工藝的研究中,已有許多關(guān)于浮游植物群落的研究[15-17],但在“集裝箱+生態(tài)池塘”模式中相關(guān)研究尚不多見[18]。本試驗(yàn)對“集裝箱+生態(tài)池塘”養(yǎng)殖系統(tǒng)中生態(tài)凈化池塘內(nèi)的浮游植物進(jìn)行了逐月調(diào)查,分析浮游植物群落結(jié)構(gòu)與養(yǎng)殖區(qū)排放水體水質(zhì)指標(biāo)的響應(yīng)關(guān)系,了解浮游植物群落在生態(tài)凈化池塘中的生態(tài)效應(yīng),總結(jié)集裝箱養(yǎng)殖排放水污染物濃度的變化特點(diǎn),以期為評價養(yǎng)殖生態(tài)效益和提高養(yǎng)殖水體污染物凈化處理效率提供參考依據(jù)。

    1" 材料和方法

    1.1" 試驗(yàn)地點(diǎn)及采樣時間

    試驗(yàn)在位于貴州省遵義市某集裝箱水產(chǎn)養(yǎng)殖基地進(jìn)行。養(yǎng)殖集裝箱(長6.1 m、寬2.8 m、高2.4 m)共20個,每個集裝箱蓄水容積為25 m3,養(yǎng)殖魚類品種為鯉(Cyprinus carpio)、大口黑鱸(Micropterus salmoides)、黃顙魚(Pelteobagrus fulvidraco)、草魚(Ctenopharyngodon idella)和斑點(diǎn)叉尾魚回(

    Lctalurus punctatus)。異位處理生態(tài)池塘面積為3 400 m2,水深約2 m。于2020年1月—12月,每月中旬在生態(tài)池塘內(nèi)采樣1次(2月未采樣),共計11次采樣?!凹b箱+生態(tài)池塘”循環(huán)水養(yǎng)殖系統(tǒng)示意圖見圖1。

    1.2" 水樣采集與處理

    在集裝箱排水管的末端(圖1中A點(diǎn))檢測水體污染物濃度?,F(xiàn)場用水質(zhì)分析儀(Hachi Hq30d)測定水溫、溶解氧和pH,用2 L聚乙烯樣本瓶采集水樣送到實(shí)驗(yàn)室,在24 h內(nèi)測定總氮、總磷、氨氮和懸浮物等指標(biāo),具體測定方法參照文獻(xiàn)[19]進(jìn)行。

    在生態(tài)池塘中的對角(圖1中B、C點(diǎn))采集浮游植物平行樣本進(jìn)行混合。定性樣本用25#浮游生物網(wǎng)在水面呈“8”字型撈取并過濾3~5 min,將網(wǎng)底管水樣移入100 mL聚乙烯樣本瓶,加入體積分?jǐn)?shù)為10%的魯哥氏液固定。定量樣本用聚乙烯樣本瓶取水面下0.5 m處的水樣2 L,加入20 mL魯哥氏液,帶回實(shí)驗(yàn)室靜置48 h后,取100 mL底層濃縮液進(jìn)行鏡檢。浮游植物物種鑒定和分類參照文獻(xiàn)[20-23]進(jìn)行。

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

    按文獻(xiàn)[24]計算浮游植物的密度和生物量。優(yōu)勢種按照優(yōu)勢度指數(shù)(Y≥0.02)確定,群落結(jié)構(gòu)多樣性采用Shannon-Wiener多樣性指數(shù)(H’)、Pielou均勻度指數(shù)(J’)和Margalef豐富度指數(shù)(d)進(jìn)行分析。各指數(shù)的計算公式如下:

    Y=(ni/N)×fi(1)

    H’=- ∑(niN)ln∑(niN)(2)

    J’=H’lnS(3)

    d=S-1lnN (4)

    式(1)~(4)中,N為群落中物種的總密度,ni為第i個物種的密度,fi為該物種在采樣點(diǎn)i出現(xiàn)的頻率,S為總種類數(shù)。

    通過浮游植物密度和生物量比較曲線(ABC曲線)可對群落處于不同干擾狀態(tài)時的情況進(jìn)行分析。若密度優(yōu)勢度曲線高于生物量優(yōu)勢度曲線,則群落處于嚴(yán)重干擾狀態(tài);若兩條曲線相交,則群落處于中等干擾狀態(tài);若生物量優(yōu)勢度曲線高于密度優(yōu)勢度曲線,則群落處于未受干擾狀態(tài)[25-26]。

    計算不同月份之間的Bray-Curtis相似性系數(shù),得到相似矩陣。通過非度量多維尺度分析(non-metric multidimensional scaling,NMDS)和等級聚類(cluster)方法實(shí)現(xiàn)對浮游植物群落的排序和分類,對群落結(jié)構(gòu)格局進(jìn)行分析。用脅迫系數(shù)(stress)衡量NMDS二維點(diǎn)陣分布效果,一般情況下,當(dāng)stress<0.05時具有很好的代表性;0.05≤stress<0.1時結(jié)果基本可信;0.1≤stress<0.2時結(jié)果具有一定的解釋意義[27]。

    用浮游植物密度進(jìn)行去趨勢對應(yīng)分析(detrended correspondence analysis,DCA),其結(jié)果:第1軸梯度值大于4時,選擇典范對應(yīng)分析(canonical correspondence analysis,CCA);介于3~4時,選用CCA和線性模型的冗余分析(redundancy analysis,RDA)均可;小于3時選用RDA。在進(jìn)行分析前,對生物數(shù)據(jù)和環(huán)境數(shù)據(jù)(pH除外)進(jìn)行l(wèi)g(x+1)轉(zhuǎn)換,運(yùn)用前選法(forward selection)和999次蒙特卡洛置換檢驗(yàn)(Monte Carlo permutation test)識別各環(huán)境因子的貢獻(xiàn)率。

    1.4" 數(shù)據(jù)處理

    試驗(yàn)數(shù)據(jù)用平均值±標(biāo)準(zhǔn)差表示。采用EXCEL 2009計算浮游植物密度、生物量和優(yōu)勢度指數(shù)。采用PRIMER 5軟件計算H’、J’和d,繪制Bray-Curtis聚類圖、NMDS分析和ABC曲線。采用Canoco 4.5軟件進(jìn)行RDA、DCA和CCA分析。設(shè)差異顯著性水平為0.05,極顯著水平為0.01。

    2" 結(jié)果和分析

    2.1" 養(yǎng)殖區(qū)排放水體水質(zhì)指標(biāo)

    由表1可見,養(yǎng)殖區(qū)排放水體溫度呈現(xiàn)規(guī)律性的季節(jié)演替。水溫年平均值為18.41 ℃,各月份中,8月最高(27.67 ℃),12月最低(6.21 ℃)。全年水體pH呈弱堿性,平均值為7.62。

    養(yǎng)殖區(qū)排放水體溶解氧的平均值為6.19 mg/L,其變化規(guī)律表現(xiàn)為春冬季高于夏秋季。氨氮平均值為2.64 mg/L,變化范圍為0.06~12.01 mg/L,其中5月份為全年最高(P<0.05),且極顯著高于10月份(P<0.01)??偟骄禐?.03 mg/L,變化范圍為1.84~15.33 mg/L,5月份為全年最高(P<0.05),其他月份之間差異不顯著。全年總磷平均值為0.68 mg/L,變化范圍為0.09~2.81 mg/L,其中4月份和5月份總磷指標(biāo)顯著高于其他月份(P<0.05),而9月份和11月份顯著低于其他月份(P<0.05),4月份和5月份極顯著高于9月份和11月份(P<0.01)。懸浮物的平均值為29.57 mg/L,其變化范圍為9.67~79.67 mg/L,

    其中4月份最高(P<0.05),12月份最低(P<0.05),且4月份與12月份差異極顯著(P<0.01)。

    2.2" 浮游植物物種組成及群落結(jié)構(gòu)

    2.2.1" 物種組成

    調(diào)查期間共鑒定出浮游植物7門49屬80種,其中藍(lán)藻門6屬6種,占浮游植物物種組成的7.50%,硅藻門7屬14種,占17.50%,綠藻門25屬38種,占47.50%,隱藻門2屬3種,占3.75%,裸藻門6屬13種,占16.25%,甲藻門2屬5種,占6.25%,金藻門1屬1種,占1.25%。

    各月份浮游植物物種數(shù)平均值為21種(見圖2),其中7月份最多(26種),4月份其次(25種),10月份最少(13種)。各月份浮游植物物種組成均以綠藻門為主,其次是硅藻門,第3為藍(lán)藻門,其他門類物種數(shù)較少。5月份未檢出硅藻門物種,4月、10月和12月未檢出隱藻門物種,1月、3月和12月未檢出裸藻門物種,5月、10月、11月和12月未檢出甲藻門物種,金藻門物種僅1月和3月檢測到。

    2.2.2" 優(yōu)勢種

    表2顯示優(yōu)勢度Y>0.001的物種。其中Y>0.02的物種有4種,即藍(lán)藻門2種,為賴氏澤絲藻(Limnothrix redekei)和湖泊假魚腥藻(Pseudoanabaena limnetica),硅藻門1種,為尖針桿藻(Synedra acus),綠藻門1種,為叢球韋斯藻(Westella botryoides)。此外,史氏棒膠藻(Rhabdogloea smithii)、梅尼小環(huán)藻(Cyclotella meneghiniana)和四尾柵藻(Scenedesmus quadricauda)、沙角衣藻(Chlamydomonas sajao)、貴州衣藻(Chlamydomonas guizhounensis)、尖尾藍(lán)隱藻(Chroomonas

    acuta)和嚙蝕隱藻(Cryptomonas erosa)為Y>0.01的物種。

    2.2.3" 密度和生物量

    浮游植物密度的變化范圍為1.17×106~43.47×106 cells/L,平均值為13.31×106 cells/L,其中3月份最高,10月份最低。藍(lán)藻門、硅藻門、綠藻門和其他門類的平均密度分別為3.88×106、0.73×106、5.98×106、1.09×106 cells/L。浮游植物生物量變化范圍為0.72~21.39 mg/L,平均值為7.16 mg/L,其中3月份最高,10月份最低。藍(lán)藻門、硅藻門、綠藻門和其他門類的平均生物量分別為0.34、2.32、1.67、2.12 mg/L。

    浮游植物密度和生物量各月份的變化情況見圖3。藍(lán)藻門密度和生物量隨時間的變化趨勢不同,3、4、9月份藍(lán)藻門密度較高,1、6月份生物量較高;硅藻門密度和生物量隨時間的變化趨勢較一致,3、12月份硅藻門密度和生物量較高;綠藻門密度和生物量的最大值均出現(xiàn)在3月,最小值均出現(xiàn)在10月;除3月份外,其他門類的密度和生物量隨時間的變化趨勢較一致。

    2.2.4" ABC曲線

    由圖4可見,ABC曲線顯示浮游植物生物量和密度曲線接近并且出現(xiàn)交叉,且W(W=0.003)為正值,表明浮游植物群落結(jié)構(gòu)受中度干擾,群落穩(wěn)定性較差,群落失去大個體優(yōu)勢種,密度和生物量的不均等程度減弱,浮游植物組成主要以密度和生物量均占優(yōu)勢的大個體物種為主,如賴氏澤絲藻、湖泊假魚腥藻、尖針桿藻和叢球韋斯藻等。

    2.2.5" 浮游植物物種多樣性

    由表3可見,生態(tài)池塘浮游植物的Shannon-Wiener多樣性指數(shù)(H’)介于1.11~2.43,平均值為1.93,7月份最高,其次是6月份,11月份最低。Pielou均勻度指數(shù)(J’)介于0.38~0.83,平均值為0.64,10月最高,其次是6月份,11月份最低。Margalef豐富度指數(shù)(d)介于0.86~1.62,平均值為1.26,7月份最高,其次是4月份,10月份最低。

    2.3" 浮游植物群落與環(huán)境因子的關(guān)系

    2.3.1" 聚類分析

    由圖5可見,各月份浮游植物群落Bray-Curtis聚類分析結(jié)果表現(xiàn)出3個相似度較大的類群,其中10、11和12月為第Ⅰ類,8月為第Ⅱ類,其余月份為第Ⅲ類。第Ⅰ類群中表現(xiàn)出2個亞類群(11、12月;10月),第Ⅲ類群也表現(xiàn)出2個亞類群(1、3、4、5、6月;7、9月)。Bray-Curtis聚類和

    聚類圖和NMDS分析

    NMDS排序結(jié)果均顯示浮游植物群落呈現(xiàn)一定的季節(jié)變化趨勢,Ⅰ類主要為冬季,Ⅱ類為夏季,Ⅲ類為春夏季。NMDS排序的脅迫系數(shù)stress=0.1(0.1≤stress<0.2),表明NMDS的二維點(diǎn)陣排序結(jié)果置信度較低,說明浮游植物群落除了受季節(jié)影響外,還與其他影響因素有關(guān)。

    2.3.2" 冗余分析(RDA)

    選用Y>0.001的浮游植物群落DCA結(jié)果為第1排序軸的梯度長度值為3.56(<4),因此,選用RDA解析浮游植物群落與環(huán)境因子的響應(yīng)關(guān)系。由圖6可見,第1軸和第2軸特征值分別為0.268和0.131,兩軸物種與環(huán)境因子之間的相關(guān)系數(shù)分別為0.981和0.998,擬合度較高,物種與環(huán)境因子的累計百分比分別為34.0%和50.7%。各環(huán)境因子與浮游植物群落結(jié)構(gòu)的響應(yīng)關(guān)系分別為氨氮(F=2.95,P=0.000),溶解氧(F=1.28,P=0.207),總氮(F=1.14,P=0.390),總磷(F=1.13,P=0.397),pH(F=1.23,P=0.308),懸浮物(F=1.21,P=0.334),水溫(F=1.13,P=0.376)??梢姲钡秋@著影響生態(tài)池塘中浮游植物群落結(jié)構(gòu)的環(huán)境因子。

    注:代碼所表示物種見表2。

    3" 討論

    3.1" 養(yǎng)殖水體污染物濃度

    隨著魚體生長和氣候的變化,魚類攝食、代謝等活動都會發(fā)生相應(yīng)變化,養(yǎng)殖排放水體中的污染物濃度也會隨之改變。根據(jù)王龍燕等[18]的研究結(jié)果,集裝箱養(yǎng)殖系統(tǒng)排放水中主要污染物為固體懸浮物(SS)和總氮(TN),且5月份的總污染指數(shù)最高。本試驗(yàn)系統(tǒng)中,1—7月,隨著水溫的不斷升高,魚類代謝不斷增強(qiáng),污染物濃度逐漸增加;7月底—8月,養(yǎng)殖區(qū)氣溫和水溫均達(dá)到一年中的最高值;為了保證魚類的耗氧需求,將達(dá)到商品魚規(guī)格的魚類逐漸上市,養(yǎng)殖容量和密度在8月底明顯降低。岳冬冬等[28]研究表明,羅非魚(Oreochromis niloticus)養(yǎng)殖產(chǎn)量與排污量具有高度的同步變化趨勢,羅非魚養(yǎng)殖產(chǎn)量增長率快于同期排污量增長率0.45%,羅非魚養(yǎng)殖產(chǎn)量每增加1%,污染物排放量將增加0.077%。不同養(yǎng)殖模式排污量差異較大[29]。因此,9月—12月,排放水體污染物濃度顯著降低的原因主要與養(yǎng)殖容量和投飼量改變有關(guān)。

    本試驗(yàn)中,生態(tài)凈化池塘對排入污染物的凈化能力相對穩(wěn)定,原因主要有兩點(diǎn):一是集裝箱排放水體進(jìn)入生態(tài)池塘前,已通過2臺轉(zhuǎn)鼓微濾機(jī)(80 m3/h)進(jìn)行了物理過濾,去除了80%的懸浮顆粒物;二是在生態(tài)池塘中種植有粉綠狐尾藻(Myriophyllum aquaticum)和豆瓣菜(Nasturtium officinale)(約占池塘面積的25%),并放養(yǎng)部分鰱(Hypophthalmichthys molitrix)和鳙(Aristichthys nobilis)。試驗(yàn)證實(shí),組合運(yùn)用物理過濾和生態(tài)凈化等方法能有效降低水體污染物濃度。本試驗(yàn)發(fā)現(xiàn),生態(tài)凈化池塘的水質(zhì)狀況隨時間呈非規(guī)律性變化,未出現(xiàn)重度污染現(xiàn)象,表明生態(tài)凈化池塘水質(zhì)總體呈現(xiàn)較好的狀態(tài),但對污染物的凈化能力有限,還需要根據(jù)養(yǎng)殖區(qū)排放水體中污染物的濃度水平及變化特征因地制宜選取凈化方法,從而提高污染物去除效率,提高凈化能力的穩(wěn)定性。

    3.2" 浮游植物物種組成及多樣性

    浮游植物是水環(huán)境健康狀況的重要指示生物,其物種多樣性能反映群落結(jié)構(gòu)的復(fù)雜性與穩(wěn)定性。不同養(yǎng)殖模式水體中的浮游植物群落組成具有一定的差異,本研究在養(yǎng)殖排放水中檢出浮游植物7門80種,其物種數(shù)及物種組成與其他養(yǎng)殖模式存在一定差異:張曉蕾等[15]在養(yǎng)殖草魚的池塘循環(huán)流水系統(tǒng)中養(yǎng)殖槽的末端、凈化池中段和凈化池末段3個位點(diǎn)共檢出浮游植物7門99種,主要為綠藻門和藍(lán)藻門,凈化區(qū)藍(lán)藻門的密度和生物量高于養(yǎng)殖區(qū),不同區(qū)域的藍(lán)藻門相對密度均遠(yuǎn)高于綠藻門及硅藻門;楊卓等[16]在不同密度組合的鰱、鳙和綠狐尾藻(Myriophyllum elatinoides)處理主養(yǎng)吉富羅非魚(GIFT Oreochromis niloticus)養(yǎng)殖池塘中共檢出浮游植物7門168種。分析不同水產(chǎn)養(yǎng)殖模式中浮游植物的群落結(jié)構(gòu)特征對充分了解水環(huán)境中生態(tài)結(jié)構(gòu)動態(tài)和推動水產(chǎn)養(yǎng)殖技術(shù)高質(zhì)量發(fā)展具有重要意義。

    本研究中,浮游植物物種類組成和密度均以綠藻門為主,分別占總物種數(shù)和總密度的47.50%和44.92%,浮游植物生物量以硅藻門為主,占總生物量組成的32.40%,調(diào)查結(jié)果與趙秀俠等[30]的結(jié)果相似。賴氏澤絲藻、湖泊假魚腥藻、尖針桿藻和叢球韋斯藻為優(yōu)勢種,其中賴氏澤絲藻和湖泊假魚腥藻是藍(lán)藻門中的小型絲狀藻,在我國多數(shù)水域中均以優(yōu)勢種出現(xiàn),如果在富營養(yǎng)化的水體中大量出現(xiàn),容易形成較高的生物量,從而導(dǎo)致藍(lán)藻水華[31-32]。減少藍(lán)藻或者控制其大量繁殖,并評估其在水環(huán)境中的潛在風(fēng)險,是一項(xiàng)具有挑戰(zhàn)且十分重要的工作[33-34]。尖針桿藻是硅藻門中營養(yǎng)物質(zhì)豐富的藻類,在自然水體中廣泛分布,在冬春季低溫、高溶解氧和貧營養(yǎng)型的水體中大量出現(xiàn)[35]。本研究在3月和12月大量檢出針桿藻種類,與其適宜于低溫輕污染水體環(huán)境的研究結(jié)果相符。有研究表明,尖針桿藻在大量出現(xiàn)的時候會結(jié)成板狀群體,在過濾池等懸浮物豐富的水體中會板結(jié)成藻類團(tuán)狀泥糊,影響凈化水質(zhì)的混凝沉淀過程,造成堵塞[36],因此在排放水處理系統(tǒng)中應(yīng)注意避免針桿藻種類的暴發(fā)。叢球韋斯藻屬于綠藻門,通常為由4個細(xì)胞呈四方形排列在1個平面上的定形群體,在軟水湖泊中大量出現(xiàn)。沙角衣藻同屬于綠藻門,能主動合成有機(jī)物,具有活化和溶解難溶性磷源的作用,對改良土壤和調(diào)節(jié)水質(zhì)具有重要作用[37],也能較好地去除水體中的磷元素[38]。研究表明,隨著水環(huán)境條件的改善,嚙蝕隱藻的密度會逐漸增加,且能較好地反映水質(zhì)狀況的變化[39]。養(yǎng)殖區(qū)水體排入生態(tài)凈化池塘使水體趨于富營養(yǎng)化,主要表現(xiàn)為喜富營養(yǎng)物種大量繁殖(如賴氏澤絲藻、湖泊假魚腥藻、針桿藻和小環(huán)藻等)和脅迫具有活化和分解營養(yǎng)物質(zhì)的物種出現(xiàn)(如沙角衣藻等)。

    3.3" 浮游植物群落結(jié)構(gòu)與環(huán)境因子的關(guān)系

    水溫是影響浮游植物群落結(jié)構(gòu)的重要環(huán)境因子之一。有研究表明,水溫能夠影響藻類生命活動的活性和速率以及藻類對水中營養(yǎng)物的利用率和自身代謝率[40]。本研究中,浮游植物密度和生物量的變化總體表現(xiàn)為春冬季高于夏秋季,但差異不顯著,說明除了水溫外,集裝箱養(yǎng)殖區(qū)養(yǎng)殖容量和投飼量變化引起的排放水污染物濃度變化也可能是影響生態(tài)凈化池塘浮游植物群落結(jié)構(gòu)特征的主要因素。

    氮和磷是浮游植物生長必需的大量元素,排放水體中氮、磷營養(yǎng)鹽的含量直接影響浮游植物的群落結(jié)構(gòu)。RDA表明,氨氮是顯著影響浮游植物群落結(jié)構(gòu)的環(huán)境因子,氨和氮在水體中的形態(tài)受到pH、溫度和溶解氧等環(huán)境因素的影響,浮游植物能夠吸收溶解無機(jī)氮和一些小分子的有機(jī)氮。韓菲爾[41]的研究表明,藻類對氮素的吸收與其個體大小、生長周期等特性有關(guān),環(huán)境因子可以通過改變關(guān)鍵酶活和光合作用影響藻類的氮吸收速率,進(jìn)而改變浮游植物群落結(jié)構(gòu)。本研究中,5月份氨氮和總氮的濃度顯著高于其他月份,表明高濃度的氮素脅迫對浮游植物群落結(jié)構(gòu)影響顯著。

    相對于湖泊、河流等開放型水體,生態(tài)凈化池塘是一個相對封閉且穩(wěn)定的水體,其影響浮游植物群落的環(huán)境因素較為單一,主要為水環(huán)境理化因子和養(yǎng)殖區(qū)的排污量[42-43]。ABC曲線分析結(jié)果表明,浮游植物群落結(jié)構(gòu)受中度干擾,群落穩(wěn)定性較差,主要影響因素是魚類飼養(yǎng)活動造成的養(yǎng)殖排放水體中污染物濃度的變化。7月份之前,集裝箱養(yǎng)殖密度較大,魚類攝食量和排污量均處于最高值,污染物排放量相對較高;到7月中旬,隨著達(dá)到商品魚規(guī)格的魚類逐漸上市,養(yǎng)殖密度不斷下降,投飼量和排污量相應(yīng)降低,浮游植物群落表現(xiàn)出Bray-Curtis聚類和NMDS排序的結(jié)果。

    4" 結(jié)論

    本研究在“集裝箱+生態(tài)池塘”循環(huán)水養(yǎng)殖系統(tǒng)中的生態(tài)凈化池塘中檢出浮游植物7門49屬80種,物種組成以綠藻門為主,其中綠藻門種類隨著月份的遞增,其所占比例不斷下降,硅藻門種類則從5月份開始,隨著月份的增加,其所占比例不斷升高。浮游植物密度、生物量、Shannon-Wiener多樣性指數(shù)、Pielou均勻度指數(shù)和Margalef豐富度指數(shù)的平均值分別為13.31×106 cells/L、7.16 mg/L、1.93、0.64和1.26。浮游植物群落結(jié)構(gòu)受環(huán)境因素的影響中等,并表現(xiàn)出一定的季節(jié)分異差異,但置信度較低。氨氮是影響浮游植物群落結(jié)構(gòu)的主要環(huán)境因子之一。

    參考文獻(xiàn)

    [1]MESQUITA M C B,PRESTES A C C,GOMES A M A, et al. Direct effects of temperature on growth of different tropical phytoplankton species[J].Microbial Ecology,2020,79(1):1-11.

    [2]ZHANG M,QIN B L,YU Y, et al. Effects of temperature fluctuation on the development of cyanobacterial dominance in spring: implication of future climate change[J].Hydrobiologia, 2016,763(1):135-146.

    [3]RYABOV A,KERIMOGLU O,LITCHMAN E, et al. Shape matters: the relationship between cell geometry and diversity in phytoplankton[J].Ecology Letters,2021,24(4):847-861.

    [4]THRANE J E,HESSEN D O,ANDERSEN T. The impact of irradiance on optimal and cellular nitrogen to phosphorus ratios in phytoplankton[J].Ecology Letters,2016,19(8):880-888.

    [5]LISTMANN L,LEROCH M,SCHLTER L, et al. Swift thermal reaction norm evolution in a key marine phytoplankton species[J].Evolutionary Applications,2016,9(9):1156-1164.

    [6]WANG Z H,ZHAO J G,ZHANG Y J, et al. Phytoplankton community structure and environmental parameters in aquaculture areas of Daya Bay, South China Sea[J].Journal of Environmental Sciences,2009,21(9):1268-1275.

    [7]OLLI K,NYMAN E,TAMMINEN T. Half-century trends in alpha and beta diversity of phytoplankton summer communities in the Helsinki Archipelago, the Baltic Sea[J].Journal of Plankton Research,2022,45(1):146-162.

    [8]HOFMANN P,CLARK A,HOFFMANN P, et al. Beyond nitrogen: phosphorus-estimating the minimum niche dimensionality for resource competition between phytoplankton[J].Ecology Letters,2021,24(4):761-771.

    [9]陳思,李藝彤,蔡文貴,等.蝦蟹混養(yǎng)池塘浮游植物群落結(jié)構(gòu)的變化特征[J].南方水產(chǎn)科學(xué),2020,16(3):79-85.

    [10]陳學(xué)洲,舒銳,謝駿,等.“集裝箱+生態(tài)池塘” 集約養(yǎng)殖與尾水高效處理技術(shù)[J].中國水產(chǎn),2020(8):67-70.

    [11]李永鋒,姚東林,何偉華.受控式集裝箱高效節(jié)能循環(huán)水養(yǎng)殖紅羅非和寶石鱸試驗(yàn)[J].科學(xué)養(yǎng)魚,2019(6):39-40.

    [12]楊麗專,劉付永忠,李勇,等.羅非魚受控式高效循環(huán)水集裝箱養(yǎng)殖技術(shù)[J].海洋與漁業(yè),2017(6):60

    [13]王磊,賈松鵬,舒銳,等.陸基推水集裝箱式水產(chǎn)養(yǎng)殖模式解析及發(fā)展展望[J].中國水產(chǎn),2020(11):51-53.

    [14]彭磊,邢道超,王惠姍,等.3種填料對集裝箱循環(huán)水養(yǎng)殖廢水處理的初步研究[J].漁業(yè)現(xiàn)代化,2017,44(5):1-6.

    [15]張曉蕾,王強(qiáng),張國奇,等.池塘循環(huán)流水養(yǎng)殖模式中浮游植物群落結(jié)構(gòu)的空間變化研究[J].南方水產(chǎn)科學(xué),2021,17(3):36-45.

    [16]楊卓,陳日釗,李文紅,等.鰱鳙+綠狐尾藻處理的羅非魚養(yǎng)殖池塘浮游植物群落特征及水質(zhì)評價[J].西南農(nóng)業(yè)學(xué)報,2021,34(4):880-888.

    [17]衛(wèi)鵬,畢相東,戴偉,等.淡水養(yǎng)殖池塘微型和超微型浮游植物的群落結(jié)構(gòu)組成[J].大連海洋大學(xué)學(xué)報,2022,37(1):113-119.

    [18]王龍燕,閔文武,王金樂,等.集裝箱循環(huán)水養(yǎng)殖系統(tǒng)中浮游植物群落結(jié)構(gòu)特征及影響因素[J].漁業(yè)現(xiàn)代化,2022,49(4):26-36.

    [19]魏復(fù)盛.水和廢水監(jiān)測分析方法[M].4版.北京:中國環(huán)境科學(xué)出版社,2002:243-284.

    [20]胡鴻鈞,魏印心.中國淡水藻類:系統(tǒng)、分類及生態(tài)[M].北京:科學(xué)出版社,2006.

    [21]韓茂森,束蘊(yùn)芳.中國淡水生物圖譜[M].北京:海洋出版社,1995.

    [22]周鳳霞,陳劍虹.淡水微型生物與底棲動物圖譜:第3版[M]. 北京:化學(xué)工業(yè)出版社,2019.

    [23]章宗涉,黃祥飛.淡水浮游生物研究方法[M].北京:科學(xué)出版社,1991.

    [24]趙文.水生生物學(xué)[M].北京:中國農(nóng)業(yè)出版社,2005.

    [25]FIELD J G,CLARKE K R,WARWICK R M.A practical strategy for analysing multispecies distribution patterns[J].Marine Ecology Progress Series,1982,8:37-52.

    [26]顏文超,宋偉華,俞存根,等.甌江口海域春秋季魚類多樣性及群落結(jié)構(gòu)研究[J].海洋湖沼通報,2018(6):132-141.

    [27]KHALAF M A,KOCHZIUS M. Changes in trophic community structure of shore fishes at an industrial site in the Gulf of Aqaba, Red Sea[J].Marine Ecology Progress Series,2002,239:287-299.

    [28]岳冬冬,李旭君,耿瑞,等.中國羅非魚養(yǎng)殖業(yè)排污量估算探析[J].漁業(yè)信息與戰(zhàn)略,2018,33(3):157-164.

    [29]韓士群,周慶,姚東瑞,等.水產(chǎn)養(yǎng)殖模式對池塘水環(huán)境和環(huán)境負(fù)荷量的影響[J].江蘇農(nóng)業(yè)學(xué)報,2018,34(3):578-584.

    [30]趙秀俠,方婷,陳金良,等.陸基推水集裝箱循環(huán)水養(yǎng)殖系統(tǒng)中浮游植物群落結(jié)構(gòu)特征[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2023,42(4):869-878.

    [31]閆蘇蘇.長壽湖浮游植物群落與水環(huán)境因子相關(guān)性研究[D].重慶:重慶大學(xué),2017.

    [32]胡韌.珠海水庫富營養(yǎng)化現(xiàn)狀、浮游植物群落特征與藍(lán)藻水華風(fēng)險分析[D].廣州:暨南大學(xué),2009.

    [33]ROSE A K,F(xiàn)ABBRO L,KINNEAR S. Cyanobacteria breakthrough: Effects of Limnothrix redekei contamination in an artificial bank filtration on a regional water supply[J].Harmful Algae,2018,76:1-10.

    [34]ZHU M L,YU G L,LI X C, et al. Taxonomic and phylogenetic evaluation of Limnothrix strains (Oscillatoriales, Cyanobacteria) by adding Limnothrix planktonica strains isolated from central China[J].Hydrobiologia,2012,698(1):367-374.

    [35]張正德,吳世萍,許亞群,等.江蘇橫山水庫浮游植物群落結(jié)構(gòu)的動態(tài)變化[J].水生態(tài)學(xué)雜志,2016,37(1):53-61.

    [36]朱雷雷.尖針桿藻對砂濾池的堵塞過程分析及控制技術(shù)研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2013.

    [37]鄧益,蘭利瓊,卿人韋,等.沙角衣藻-雜交水稻實(shí)驗(yàn)體系在不同供磷條件下生長情況的研究[J].土壤通報,2003,34(2):133-135.

    [38]逯多.利用藻類除去污水中的磷[D].成都:四川大學(xué),2001.

    [39]楊姣.惠州西湖浮游植物群落對生態(tài)修復(fù)的響應(yīng)[D].廣州:暨南大學(xué),2009.

    [40]王志紅,崔福義,安全.水溫與營養(yǎng)值對水庫藻華態(tài)勢的影響[J].生態(tài)環(huán)境,2005,14(1):10-15.

    [41]韓菲爾.太湖水體浮游植物氮素吸收過程及其影響因素[D].蘇州:蘇州科技大學(xué),2019.

    [42]蒲煒佳,董世鵬,張東旭,等.三疣梭子蟹池塘綜合養(yǎng)殖系統(tǒng)浮游植物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系[J].中國水產(chǎn)科學(xué),2022,29(4):549-561.

    [43]ZHANG M,DONG J,GAO Y N, et al. Patterns of phytoplankton community structure and diversity in aquaculture ponds, Henan, China[J].Aquaculture,2021,544:737078.

    Diversity analysis of phytoplankton community structure in container

    aquaculture ecological purification pond

    MIN Wenwu1,2, WANG Longyan3, CHEN Feixiong1, ZHOU Qichun1,2

    (1. Fisheries Research Institute,Guizhou Academy of Agricultural Sciences,Guiyang" 550025,China;

    2. Guizhou Special Aquatic Products

    Engineering Technology Center,Guiyang" 550025,China;

    3. Guizhou Institute of Biology,Guiyang" 550009,China)

    Abstract: To understand the pollution status of water discharged from container aquaculture and the dynamic changes of phytoplankton community structure in the ecological purification pond,an investigation and sampling were conducted in the ecological purification pond of the “container+ecological pond” recirculating aquaculture system in 2020.The species diversity index was used to analyze the characteristics of the phytoplankton community structure.The abundance and biomass comparison curve(ABC curve) as well as cluster analysis were used to investigate the disturbance degree of the phytoplankton community and its temporal clustering difference,and the redundancy analysis(RDA) was used to explore the response relationship between the phytoplankton community and environmental factors.The results showed that a total of 80 species of phytoplankton community belonging to 49 genera and 7 phyla were detected in the ecological purification pond,with Chlorophyta being the predominant phytoplankton group(47.5%),followed by Bacillariophyta(17.5%).The number of phytoplankton species reached the highest level in July(26 species) and the lowest in October(13 species).Dominant species included Limnothrix redekei, Pseudoanabaena limnetica,Synedra acus and Westella botryoides.The average values of phytoplankton density,biomass,Shannon-Wiener diversity index,Pielou evenness index and Margalef richness index were 13.31×106cells/L,7.16 mg/L,1.93,0.64 and 1.26,respectively.The ABC curve indicated that the phytoplankton community structure was moderately disturbed,with community composition mainly dominated by species with high density and biomass.Cluster analysis showed that the phytoplankton community exhibited some insignificant seasonal differentiation characteristics.RDA revealed that ammonia nitrogen was significantly related to the phytoplankton community structure.

    Key words: aquaculture; ecological purification pond; phytoplankton; community structure; species diversity

    猜你喜歡
    浮游植物物種多樣性群落結(jié)構(gòu)
    海南省紅樹林區(qū)浮游動物多樣性的初步研究
    海南省紅樹林區(qū)大型底棲動物多樣性的初步研究
    江西大崗山木荷天然林群落物種多樣性分析
    綠色科技(2016年21期)2016-12-27 10:53:26
    納板河保護(hù)區(qū)不同林分群落喬木物種多樣性研究
    用固定樣地法監(jiān)測銀瓶山闊葉林物種多樣性
    塔里木河阿拉爾段及多浪水庫浮游生物的調(diào)查研究
    輝騰錫勒風(fēng)電場區(qū)域內(nèi)物種多樣性研究
    渭河陜西段浮游植物群落結(jié)構(gòu)及水質(zhì)評價
    濟(jì)南流域春季浮游植物多樣性的研究
    峽江樞紐庫區(qū)水生生物調(diào)查及水質(zhì)評價
    黄色日韩在线| 99热精品在线国产| 日韩,欧美,国产一区二区三区 | 国产高清激情床上av| 最近视频中文字幕2019在线8| 日韩亚洲欧美综合| 级片在线观看| 亚洲av二区三区四区| 亚洲国产精品sss在线观看| 国产蜜桃级精品一区二区三区| 国产精品一区二区性色av| 国产精华一区二区三区| 国产蜜桃级精品一区二区三区| av国产免费在线观看| 尾随美女入室| 99久久精品一区二区三区| 老师上课跳d突然被开到最大视频| 99热这里只有是精品在线观看| 少妇熟女欧美另类| 亚洲欧美中文字幕日韩二区| 久久国内精品自在自线图片| 91久久精品国产一区二区成人| 国产精品久久久久久av不卡| 岛国在线免费视频观看| 欧美zozozo另类| 亚洲婷婷狠狠爱综合网| 国产精品久久视频播放| 国产久久久一区二区三区| 成人高潮视频无遮挡免费网站| 国产精品福利在线免费观看| 亚洲七黄色美女视频| 美女高潮的动态| 免费av观看视频| 在线观看av片永久免费下载| 久久久精品94久久精品| 一进一出抽搐动态| 国产午夜精品一二区理论片| 日韩人妻高清精品专区| 亚洲欧洲日产国产| 久久99热6这里只有精品| 18+在线观看网站| 大型黄色视频在线免费观看| 麻豆国产97在线/欧美| 亚洲av免费高清在线观看| 真实男女啪啪啪动态图| 日本免费一区二区三区高清不卡| 天堂√8在线中文| 亚洲av免费高清在线观看| 国产成人a区在线观看| 高清午夜精品一区二区三区 | 国产精品久久电影中文字幕| 中文亚洲av片在线观看爽| 国产精品久久久久久精品电影| av卡一久久| 国产精品永久免费网站| 国产精品国产高清国产av| 久久人人精品亚洲av| 一边亲一边摸免费视频| 18禁在线播放成人免费| 99热6这里只有精品| 久久精品国产99精品国产亚洲性色| 亚洲精华国产精华液的使用体验 | 欧美极品一区二区三区四区| 夫妻性生交免费视频一级片| 国产成人91sexporn| videossex国产| 久久精品国产鲁丝片午夜精品| 国产精品一区二区在线观看99 | 成人亚洲欧美一区二区av| 日韩制服骚丝袜av| 日韩欧美精品免费久久| 午夜亚洲福利在线播放| 人妻久久中文字幕网| 少妇高潮的动态图| 99热只有精品国产| 99国产精品一区二区蜜桃av| 日韩欧美精品免费久久| 亚洲欧美中文字幕日韩二区| 欧美xxxx性猛交bbbb| 国产69精品久久久久777片| av黄色大香蕉| 亚洲三级黄色毛片| 亚洲三级黄色毛片| 国产成人精品一,二区 | 日韩欧美 国产精品| 在线观看av片永久免费下载| 少妇丰满av| 国产视频首页在线观看| av在线播放精品| 能在线免费看毛片的网站| 一边摸一边抽搐一进一小说| 97热精品久久久久久| 毛片女人毛片| 男人的好看免费观看在线视频| 国产精品一及| 精品一区二区免费观看| 亚洲av男天堂| 九色成人免费人妻av| 九色成人免费人妻av| 亚洲av第一区精品v没综合| 一级毛片我不卡| 一个人观看的视频www高清免费观看| 国产久久久一区二区三区| 免费无遮挡裸体视频| 2021天堂中文幕一二区在线观| 日日啪夜夜撸| 亚洲内射少妇av| 岛国在线免费视频观看| 91在线精品国自产拍蜜月| 国产老妇伦熟女老妇高清| 亚洲久久久久久中文字幕| 欧美人与善性xxx| 国产麻豆成人av免费视频| 亚洲成a人片在线一区二区| 一本久久精品| 晚上一个人看的免费电影| 婷婷色av中文字幕| 又爽又黄无遮挡网站| 久久精品国产99精品国产亚洲性色| 国产精品久久电影中文字幕| 国产av麻豆久久久久久久| 中国美白少妇内射xxxbb| 国产伦理片在线播放av一区 | 亚洲精品乱码久久久久久按摩| 黄片无遮挡物在线观看| 热99在线观看视频| 国产成年人精品一区二区| a级毛片a级免费在线| 岛国在线免费视频观看| 一级黄片播放器| 日日摸夜夜添夜夜添av毛片| av免费在线看不卡| 久久精品国产亚洲网站| 小说图片视频综合网站| 舔av片在线| 午夜精品国产一区二区电影 | 91久久精品国产一区二区成人| 中文字幕人妻熟人妻熟丝袜美| 日本-黄色视频高清免费观看| 国产黄a三级三级三级人| 青春草国产在线视频 | 激情 狠狠 欧美| 男人狂女人下面高潮的视频| 欧美三级亚洲精品| 欧美一区二区亚洲| 久久人妻av系列| 在线观看av片永久免费下载| 嫩草影院入口| 看黄色毛片网站| 婷婷色综合大香蕉| 久久久色成人| 亚洲国产精品合色在线| 国产亚洲91精品色在线| 久久精品国产99精品国产亚洲性色| 夫妻性生交免费视频一级片| 成人漫画全彩无遮挡| 成人二区视频| 赤兔流量卡办理| 精品人妻一区二区三区麻豆| 99热只有精品国产| 色哟哟·www| 舔av片在线| 非洲黑人性xxxx精品又粗又长| 天堂av国产一区二区熟女人妻| 永久网站在线| 亚洲18禁久久av| 最后的刺客免费高清国语| 女人被狂操c到高潮| 日韩视频在线欧美| 老熟妇乱子伦视频在线观看| eeuss影院久久| 免费人成视频x8x8入口观看| 成年av动漫网址| 亚洲av二区三区四区| 丰满的人妻完整版| 插阴视频在线观看视频| av黄色大香蕉| 国产熟女欧美一区二区| 99热只有精品国产| 久久久久网色| 精品久久久久久久人妻蜜臀av| 一区二区三区高清视频在线| 亚洲自偷自拍三级| 国产精品一区www在线观看| 久久这里有精品视频免费| 亚洲激情五月婷婷啪啪| 一级毛片电影观看 | 亚洲欧洲日产国产| 成人性生交大片免费视频hd| 国产激情偷乱视频一区二区| 国产高潮美女av| 亚洲自偷自拍三级| 欧美色视频一区免费| 久久久久久久亚洲中文字幕| 亚洲人成网站在线播放欧美日韩| 亚洲美女视频黄频| 日本爱情动作片www.在线观看| 国产亚洲欧美98| 久久久久久久久中文| 又粗又爽又猛毛片免费看| 久久这里只有精品中国| 精品免费久久久久久久清纯| 啦啦啦啦在线视频资源| 久久久a久久爽久久v久久| 99视频精品全部免费 在线| 美女大奶头视频| 女人十人毛片免费观看3o分钟| 99热这里只有精品一区| 国产精品.久久久| 国产一区二区在线观看日韩| 我要搜黄色片| a级毛色黄片| 99精品在免费线老司机午夜| 国产精品综合久久久久久久免费| 色吧在线观看| 好男人视频免费观看在线| 如何舔出高潮| 精品少妇黑人巨大在线播放 | 国产极品精品免费视频能看的| 狂野欧美激情性xxxx在线观看| 久久久色成人| 青春草视频在线免费观看| 两个人的视频大全免费| 亚洲国产欧美在线一区| 最后的刺客免费高清国语| 国内久久婷婷六月综合欲色啪| 国产精品.久久久| 18禁黄网站禁片免费观看直播| 激情 狠狠 欧美| 久久99热6这里只有精品| 欧美激情国产日韩精品一区| 国产午夜精品久久久久久一区二区三区| 69av精品久久久久久| 免费黄网站久久成人精品| 人体艺术视频欧美日本| АⅤ资源中文在线天堂| 久久人人爽人人片av| 一个人看视频在线观看www免费| 草草在线视频免费看| 久久6这里有精品| 校园春色视频在线观看| 看十八女毛片水多多多| 黄片wwwwww| 亚洲av第一区精品v没综合| 亚洲成a人片在线一区二区| 国产精品.久久久| 亚洲18禁久久av| 91精品国产九色| 99久久久亚洲精品蜜臀av| 91在线精品国自产拍蜜月| 天堂网av新在线| av国产免费在线观看| 欧美色视频一区免费| 欧美激情久久久久久爽电影| 乱人视频在线观看| 丝袜喷水一区| 2022亚洲国产成人精品| 一夜夜www| 人妻少妇偷人精品九色| 日韩欧美精品v在线| 国产精品精品国产色婷婷| 成人av在线播放网站| 成人美女网站在线观看视频| 在线观看av片永久免费下载| 18禁在线播放成人免费| 乱系列少妇在线播放| 99热这里只有精品一区| 中文亚洲av片在线观看爽| 99热网站在线观看| 一个人免费在线观看电影| 国产麻豆成人av免费视频| 日韩一区二区视频免费看| 麻豆久久精品国产亚洲av| 国产精品一及| 久久午夜亚洲精品久久| 精品人妻一区二区三区麻豆| 熟女电影av网| 亚洲一级一片aⅴ在线观看| 日本一二三区视频观看| 天堂影院成人在线观看| 黄色一级大片看看| 免费av毛片视频| 免费看av在线观看网站| 2021天堂中文幕一二区在线观| 久久午夜亚洲精品久久| 婷婷亚洲欧美| 综合色丁香网| 麻豆成人av视频| 99热只有精品国产| 久久人妻av系列| 成人国产麻豆网| 夜夜夜夜夜久久久久| 变态另类成人亚洲欧美熟女| 欧美在线一区亚洲| 免费看光身美女| 91久久精品国产一区二区三区| 国产久久久一区二区三区| 久久国产乱子免费精品| 国产单亲对白刺激| 日本与韩国留学比较| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产探花在线观看一区二区| 国产私拍福利视频在线观看| 亚洲第一电影网av| 丰满的人妻完整版| 精品久久久久久久久久久久久| 精品久久久久久久久亚洲| 久久精品国产亚洲网站| 日韩精品青青久久久久久| 久久午夜亚洲精品久久| 身体一侧抽搐| 99久国产av精品国产电影| 在线a可以看的网站| av在线播放精品| 亚洲无线观看免费| 欧美一区二区国产精品久久精品| 欧美潮喷喷水| 精品久久国产蜜桃| 国产探花在线观看一区二区| 在线免费观看的www视频| 国产久久久一区二区三区| www日本黄色视频网| 日本黄色片子视频| 欧美日韩在线观看h| 人妻系列 视频| 大型黄色视频在线免费观看| 51国产日韩欧美| 身体一侧抽搐| 成人性生交大片免费视频hd| 日韩一本色道免费dvd| 精品欧美国产一区二区三| 女同久久另类99精品国产91| 日韩高清综合在线| 中文亚洲av片在线观看爽| 中国美女看黄片| 国产精品.久久久| 两个人视频免费观看高清| 91狼人影院| 欧美性猛交黑人性爽| 啦啦啦韩国在线观看视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产三级在线视频| 日本撒尿小便嘘嘘汇集6| 亚洲精品亚洲一区二区| 少妇熟女欧美另类| 国产精品福利在线免费观看| 美女高潮的动态| 久久99精品国语久久久| a级毛色黄片| av天堂在线播放| 欧美xxxx性猛交bbbb| ponron亚洲| 在线免费十八禁| 永久网站在线| 国产人妻一区二区三区在| 两个人的视频大全免费| 国产麻豆成人av免费视频| 日韩一本色道免费dvd| 国产美女午夜福利| 男人舔奶头视频| 欧美xxxx性猛交bbbb| 中文字幕av成人在线电影| 亚洲成人久久爱视频| 麻豆久久精品国产亚洲av| 1024手机看黄色片| 日韩人妻高清精品专区| 国产视频内射| av国产免费在线观看| 极品教师在线视频| 久久精品国产自在天天线| 午夜福利成人在线免费观看| 精品一区二区免费观看| 国产在线精品亚洲第一网站| 久久人人爽人人爽人人片va| a级一级毛片免费在线观看| 男女边吃奶边做爰视频| 欧美日本视频| 一级毛片久久久久久久久女| 久久久久国产网址| 天天躁日日操中文字幕| 成年版毛片免费区| 一本久久中文字幕| 国产高潮美女av| 熟女人妻精品中文字幕| 99国产极品粉嫩在线观看| 99热这里只有是精品在线观看| 亚洲乱码一区二区免费版| 婷婷色av中文字幕| 一区福利在线观看| 尤物成人国产欧美一区二区三区| 97超视频在线观看视频| av专区在线播放| 成年女人永久免费观看视频| 国产精品不卡视频一区二区| 91在线精品国自产拍蜜月| 黄色一级大片看看| 国产成人影院久久av| 天堂中文最新版在线下载 | 网址你懂的国产日韩在线| 精品久久久久久久末码| 欧美三级亚洲精品| 国产精品美女特级片免费视频播放器| 神马国产精品三级电影在线观看| 国产一级毛片在线| 久久99热这里只有精品18| 中文在线观看免费www的网站| av免费在线看不卡| 成人一区二区视频在线观看| or卡值多少钱| 人妻制服诱惑在线中文字幕| 一级毛片久久久久久久久女| 99久久九九国产精品国产免费| 亚洲高清免费不卡视频| 亚洲中文字幕日韩| 国产片特级美女逼逼视频| 日韩一本色道免费dvd| 99热这里只有是精品在线观看| 精品人妻视频免费看| 欧美另类亚洲清纯唯美| 丰满人妻一区二区三区视频av| 国产在线精品亚洲第一网站| 成人特级av手机在线观看| 亚洲成av人片在线播放无| 免费观看人在逋| 久久久久国产网址| 午夜免费激情av| 不卡一级毛片| 欧美色欧美亚洲另类二区| 18禁黄网站禁片免费观看直播| 高清毛片免费看| 国产精品免费一区二区三区在线| 日韩大尺度精品在线看网址| 欧美性猛交╳xxx乱大交人| 天天一区二区日本电影三级| 日韩欧美精品免费久久| 男的添女的下面高潮视频| 美女 人体艺术 gogo| 寂寞人妻少妇视频99o| 欧美成人a在线观看| 国产精品不卡视频一区二区| 午夜免费男女啪啪视频观看| 99视频精品全部免费 在线| 欧美激情在线99| kizo精华| 高清日韩中文字幕在线| 亚洲天堂国产精品一区在线| 午夜a级毛片| 校园春色视频在线观看| av国产免费在线观看| 亚洲欧美精品综合久久99| 日韩中字成人| 日本黄色片子视频| 成人毛片60女人毛片免费| 国产亚洲av嫩草精品影院| av福利片在线观看| 91在线精品国自产拍蜜月| 国产亚洲精品av在线| 少妇的逼好多水| 又黄又爽又刺激的免费视频.| 最新中文字幕久久久久| 嘟嘟电影网在线观看| 国产三级中文精品| 国产成人aa在线观看| 色哟哟哟哟哟哟| 亚洲最大成人手机在线| 欧美日韩一区二区视频在线观看视频在线 | 波多野结衣高清作品| 国产精品伦人一区二区| 色哟哟哟哟哟哟| 亚洲欧美清纯卡通| 岛国毛片在线播放| 偷拍熟女少妇极品色| 国产成人午夜福利电影在线观看| 99在线人妻在线中文字幕| 亚洲人成网站高清观看| 久久99热这里只有精品18| 欧美高清性xxxxhd video| 在线观看一区二区三区| 成人鲁丝片一二三区免费| 此物有八面人人有两片| 亚洲va在线va天堂va国产| 九草在线视频观看| 日本av手机在线免费观看| h日本视频在线播放| 男人和女人高潮做爰伦理| 欧美又色又爽又黄视频| 亚洲国产精品久久男人天堂| 有码 亚洲区| 亚洲在久久综合| 亚洲精品国产av成人精品| 男人舔奶头视频| 午夜精品国产一区二区电影 | 日本三级黄在线观看| 天堂√8在线中文| 伦精品一区二区三区| 欧美成人a在线观看| 一进一出抽搐gif免费好疼| 色视频www国产| 精华霜和精华液先用哪个| 国产高清三级在线| 听说在线观看完整版免费高清| 看十八女毛片水多多多| 国产精品一区www在线观看| av免费观看日本| 男女那种视频在线观看| 国产精品精品国产色婷婷| 欧美最新免费一区二区三区| 简卡轻食公司| 男人舔女人下体高潮全视频| 啦啦啦韩国在线观看视频| 日本爱情动作片www.在线观看| 你懂的网址亚洲精品在线观看 | 日日摸夜夜添夜夜添av毛片| 欧美成人一区二区免费高清观看| 大型黄色视频在线免费观看| 高清毛片免费看| 国产一区二区亚洲精品在线观看| 白带黄色成豆腐渣| 麻豆国产av国片精品| 97人妻精品一区二区三区麻豆| 99久久成人亚洲精品观看| 久久久久久大精品| 免费av不卡在线播放| 最近2019中文字幕mv第一页| 亚洲欧洲国产日韩| 成人午夜精彩视频在线观看| 一个人观看的视频www高清免费观看| 日韩欧美 国产精品| 亚洲av不卡在线观看| 亚洲国产欧洲综合997久久,| 国产精品乱码一区二三区的特点| 在现免费观看毛片| 亚洲精品成人久久久久久| 国产日韩欧美在线精品| 日韩成人av中文字幕在线观看| 亚洲中文字幕日韩| 女人被狂操c到高潮| 亚洲av免费在线观看| 熟妇人妻久久中文字幕3abv| 精品欧美国产一区二区三| 日本色播在线视频| 欧美最黄视频在线播放免费| 丝袜美腿在线中文| 男女啪啪激烈高潮av片| 别揉我奶头 嗯啊视频| 精品久久国产蜜桃| 国产精品一区二区在线观看99 | 免费看美女性在线毛片视频| 日本-黄色视频高清免费观看| 国产视频内射| 亚洲精品色激情综合| 日本熟妇午夜| 成人性生交大片免费视频hd| 国产精品av视频在线免费观看| 国产在线男女| 日韩欧美精品免费久久| 中文欧美无线码| 亚洲av成人精品一区久久| а√天堂www在线а√下载| 日韩欧美 国产精品| 蜜桃久久精品国产亚洲av| 欧美性猛交黑人性爽| 国产老妇女一区| 国产69精品久久久久777片| 国产精品人妻久久久影院| 成人高潮视频无遮挡免费网站| 日韩亚洲欧美综合| 久久精品久久久久久噜噜老黄 | 亚洲真实伦在线观看| 国产高清三级在线| 日日干狠狠操夜夜爽| 22中文网久久字幕| 99久久成人亚洲精品观看| 偷拍熟女少妇极品色| 身体一侧抽搐| 蜜桃久久精品国产亚洲av| a级毛片免费高清观看在线播放| 欧美日韩综合久久久久久| 夜夜夜夜夜久久久久| 国产男人的电影天堂91| 国产白丝娇喘喷水9色精品| 一级黄片播放器| 97超视频在线观看视频| 亚洲国产精品合色在线| 少妇人妻精品综合一区二区 | 免费黄网站久久成人精品| or卡值多少钱| 午夜福利在线观看免费完整高清在 | www.av在线官网国产| 日韩高清综合在线| 成人毛片a级毛片在线播放| 国产精品一区www在线观看| 美女大奶头视频| 成熟少妇高潮喷水视频| 国产精品免费一区二区三区在线| 赤兔流量卡办理| 毛片一级片免费看久久久久| ponron亚洲| 一卡2卡三卡四卡精品乱码亚洲| 99热网站在线观看| 亚洲欧美日韩高清在线视频| 床上黄色一级片| 免费av观看视频| 亚洲七黄色美女视频| 在线观看免费视频日本深夜| av女优亚洲男人天堂| 人人妻人人澡欧美一区二区| 蜜桃久久精品国产亚洲av| 国产男人的电影天堂91| 亚洲无线观看免费| 美女内射精品一级片tv| 国产成人精品久久久久久| 99热6这里只有精品| 午夜a级毛片| 免费观看精品视频网站|