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

    南黃海春季大型底棲動物優(yōu)勢種生態(tài)位

    2015-03-11 06:41:44彭松耀李新正王洪法張寶琳
    生態(tài)學(xué)報(bào) 2015年6期
    關(guān)鍵詞:沙蠶黃海生境

    彭松耀,李新正,王洪法,張寶琳

    1 中國科學(xué)院海洋研究所, 青島 266071 2 中國科學(xué)院大學(xué), 北京 100039

    南黃海春季大型底棲動物優(yōu)勢種生態(tài)位

    彭松耀1,2,李新正1,*,王洪法1,張寶琳1

    1 中國科學(xué)院海洋研究所, 青島 266071 2 中國科學(xué)院大學(xué), 北京 100039

    2011年4月對南黃海大型底棲動物進(jìn)行定量采樣調(diào)查,獲得大型底棲動物105種,篩選獲得優(yōu)勢種20種,隸屬3門3綱17科。南黃海大型底棲動物優(yōu)勢種以運(yùn)動、有頜、表層沉積物取食者相對豐度較高,占27.6%。對選取的20個(gè)優(yōu)勢種以Shannon-Wiener指數(shù)為基礎(chǔ)進(jìn)行了生態(tài)位寬度測定,以Pianka 重疊指數(shù)為基礎(chǔ)進(jìn)行了生態(tài)位重疊值分析,以平均邊緣指數(shù)(OMI)和耐受指數(shù)(TOI)為基礎(chǔ)分析了物種生態(tài)位與環(huán)境之間的相互關(guān)系。結(jié)果表明:南黃海大型底棲動物群落優(yōu)勢種生態(tài)位寬度變化范圍為1.24—2.15,生態(tài)位寬度值較高的有薄索足蛤、蜈蚣?xì)W努菲蟲、掌鰓索沙蠶、淺水薩氏真蛇尾、黃海刺梳鱗蟲和寡節(jié)甘吻沙蠶;優(yōu)勢種OMI指數(shù)變化范圍為0.23—4.95,OMI指數(shù)值較高的有擬特須蟲和細(xì)弱吻沙蠶;優(yōu)勢種耐受指數(shù)值變化范圍為0.13—3.85,值較高的有深鉤毛蟲、寡鰓齒吻沙蠶和黃海刺梳鱗蟲;優(yōu)勢種生態(tài)位之間的重疊值不均一,在0—0.95之間,長葉尖索沙蠶和掌鰓索沙蠶生態(tài)位重疊值較高,為0.95;平均邊緣指數(shù)分析隨機(jī)置換顯著性檢驗(yàn)表明:水深、水溫、鹽度、沉積物粒徑、沉積物總有機(jī)碳和總氮與長吻沙蠶、細(xì)弱吻沙蠶、背蚓蟲、角海蛹和淺水薩氏真蛇尾生態(tài)位之間的相互關(guān)系較顯著(P<0.05)。對群落優(yōu)勢種進(jìn)行聚類和多維度排序分析,結(jié)果表明20個(gè)物種可分為廣布種、典型生境種和特化種,主坐標(biāo)分析反映了優(yōu)勢種在不同資源位點(diǎn)的分布狀況,支持聚類和排序的結(jié)果。研究表明,優(yōu)勢種生態(tài)位的寬度、優(yōu)勢種之間的重疊值與物種攝食功能類群、生活型、資源位上的分布、物種數(shù)量及生境狀況密切相關(guān),反映了大型底棲動物群落中各物種對生境資源的不同利用能力。

    大型底棲動物; 優(yōu)勢種; 生態(tài)位; 生態(tài)位寬度; 生態(tài)位重疊值; 南黃海

    生態(tài)位是生態(tài)學(xué)中重要概念和核心思想,生態(tài)位研究在理解群落結(jié)構(gòu)功能、群落內(nèi)物種種間關(guān)系、生物多樣性、群落動態(tài)演替和種群進(jìn)化等方面具有重要的作用[1]。生態(tài)位定義由Hutchinson于1957年提出[2],即某種生物的生態(tài)位是由所有生物和非生物環(huán)境條件所組成的n維空間中該生物得以繼續(xù)生存的范圍。生態(tài)位的定義促進(jìn)了生態(tài)位分化、生態(tài)位寬度和生態(tài)位重疊的研究[1,3]。例如Colwell和Futuyma[4]通過分析不同生境狀況下物種個(gè)體分布的均勻性來估算生態(tài)位的寬度。Feinsinger等[5]使用生態(tài)相似性比例指數(shù)來計(jì)算生態(tài)位的寬度。Pianka通過計(jì)算物種間生態(tài)位的重疊值[3],分析不同生物類群資源利用的交叉狀況。近年來,隨著多元統(tǒng)計(jì)分析的廣泛運(yùn)用,OMI指數(shù)(Outlying Mean Index,邊緣平均指數(shù))被提出[6],用于分析物種生態(tài)位與環(huán)境之間的相關(guān)性。例如Korsu等[7]對芬蘭北部溪流土著魚類和入侵種魚類的生態(tài)位進(jìn)行研究,發(fā)現(xiàn)生態(tài)位寬度及生態(tài)位重疊值反映了物種對生態(tài)環(huán)境資源的利用能力;Gascón等[8]對葉足動物生態(tài)位寬度與棲息地的退化進(jìn)行了分析;Olivier等[9]研究發(fā)現(xiàn)地下水暗層生物生態(tài)位可以反映生物對環(huán)境污染的敏感性。目前,國內(nèi)學(xué)者對大型底棲動物生態(tài)位的研究報(bào)道較多[10-12],但研究區(qū)域多集中在潮間帶。本研究選取南黃海為研究對象,對該海域大型底棲動物優(yōu)勢種的生態(tài)位進(jìn)行分析,以期了解這些物種種在群落中的關(guān)系和相對地位,以及對環(huán)境的適應(yīng)性和對資源的利用特征。

    1 材料與方法

    1.1 調(diào)查區(qū)域

    圖1 南黃海站位圖 Fig.1 Sampling stations of macrozoobenthos in the southern Yellow Sea

    南黃海為我國東部重要海域之一。它北起山東半島東端成山角至朝鮮半島西部白翎島一線(約北緯37.5°),南至長江口的啟東角,東鄰朝鮮半島,西接魯、蘇海岸。本研究在南黃海設(shè)置31個(gè)采樣站(圖1),分布于E121°—125°,N32°—37°之間水域。于2011年4月搭乘“科學(xué)三號”科學(xué)考察船對大型底棲動物進(jìn)行調(diào)查采樣。

    1.2 底棲動物的采集及鑒定

    每個(gè)采樣站用0.1m2Gray-O′Hara箱式采泥器重復(fù)采樣2次,計(jì)算為一個(gè)樣品,用孔徑為0.5 mm的套篩篩選,所獲樣品用75%的酒精固定后帶回實(shí)驗(yàn)室,隨后進(jìn)行種類鑒定、個(gè)體計(jì)數(shù)、稱重(使用0.001g精度天平),并對所獲數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。

    1.3 理化指標(biāo)的測定及方法

    底層水溫、水深和鹽度用多功能溫鹽深儀(CTD)現(xiàn)場同步測定。于每個(gè)樣點(diǎn)采集表層5 cm沉積物樣,低溫保存。帶回實(shí)驗(yàn)室后測定其物理和化學(xué)指標(biāo),沉積物中值粒徑分析采用Cilas(940L)型激光粒度儀進(jìn)行測定,沉積物總氮(TN)和總有機(jī)碳(TOC)的測定按照《海洋監(jiān)測規(guī)范(GB17378.5—2007)》[14]的測定方法。

    1.4 數(shù)據(jù)分析處理

    大型底棲動物優(yōu)勢種分析采用物種豐度等級曲線,分析時(shí)樣點(diǎn)物種豐度數(shù)據(jù)進(jìn)行缺失轉(zhuǎn)化[8];生態(tài)位寬度采用以Shannon-Wiener多樣性指數(shù)為基礎(chǔ)的生態(tài)位寬度指數(shù)[15]:

    生態(tài)位重疊指數(shù)采用Pianka指數(shù)[3]:

    式中,Pij和Pkj為種i和種k的個(gè)體數(shù)占j站位個(gè)體數(shù)的比例,R為環(huán)境資源位點(diǎn)的總數(shù)(本研究中R為31,每個(gè)取樣位點(diǎn)相當(dāng)于一個(gè)環(huán)境位點(diǎn))。Bi為生態(tài)位寬度,取值范圍為0—lnR,Bi值越大表示該種生態(tài)位寬度越大;Qik為生態(tài)位重疊值,取值范圍為0—1,Qik值越大表示生態(tài)位重疊值越大。為了分析物種生態(tài)位寬度和生境之間關(guān)系,同時(shí)運(yùn)用平均邊緣指數(shù)(OMI)和耐受指數(shù)(TOL)[6,16]分別對研究區(qū)內(nèi)的優(yōu)勢種的生態(tài)位和生態(tài)位寬度進(jìn)行計(jì)算,物種的生態(tài)位平均邊緣指數(shù)是通過計(jì)算該物種與假定“物種”之間的距離,該假定“物種”要求適合的生境是所調(diào)查區(qū)域的平均生境,物種的耐受指數(shù)是通過計(jì)算物種調(diào)查區(qū)域物種沿環(huán)境梯度分布范圍的寬窄,同時(shí)檢驗(yàn)各物種的顯著性(本文計(jì)算各優(yōu)勢種的P值),分析P值小于0.05并且生態(tài)位寬度和生態(tài)位對該物種變異解釋度達(dá)50%以上的物種的生態(tài)位、生態(tài)位寬度,計(jì)算時(shí)所有環(huán)境因子進(jìn)行l(wèi)og(x+1)轉(zhuǎn)換。

    以上所有分析都在R語言環(huán)境中完成。采用英國普利茅斯海洋研究所開發(fā)的PRIMER6.0軟件進(jìn)行物種等級聚類、非度量多維標(biāo)度排序(NMDS)和主坐標(biāo)分析[17]。分析時(shí)優(yōu)勢種豐度數(shù)據(jù)均進(jìn)行平方根轉(zhuǎn)換,聚類和排序采用種間歐式距離,主坐標(biāo)分析采用Bray-Curtis距離。

    1.5 底棲生物群落位置及優(yōu)勢種功能群類型的劃分

    南黃海底棲生物群落位置劃分根據(jù)Zhang等的研究[13],劃分為4個(gè)群落(圖1),分別為冷水團(tuán)群落、混合群落、廣溫群落和長江口群落。底棲生物優(yōu)勢種功能群類型(FFG, HBG)的劃分根據(jù)Fauchald 和Jumars[18]以及Bonsdorff[19]的研究,劃分為9種功能群類型(表2),即:(1)運(yùn)動、有頜、表層沉積物取食者(SMJ),(2)運(yùn)動、有頜、肉食者(CMJ),(3)半運(yùn)動、有頜、肉食者(CDJ),(4)運(yùn)動、無頜、表層沉積物取食者(SMX),(5)運(yùn)動、無頜、穴居者(BMX),(6)半運(yùn)動、有觸手、表層沉積物取食者(SDT),(7)固著生活、有觸手、表層沉積物取食者(SST),(8)運(yùn)動、無頜、肉食者(CMX),(9)半運(yùn)動、泵吸、濾食者(FDP)。

    2 結(jié)果

    2.1 種類組成

    調(diào)查共發(fā)現(xiàn)大型底棲動物105種,20個(gè)優(yōu)勢種隸屬3門3綱17科(表1),鑒于物種數(shù)較多,本文通過物種豐度等級曲線選取20個(gè)優(yōu)勢種(圖2)進(jìn)行生態(tài)位、生態(tài)位重疊度、OMI指數(shù)、等級聚類、非度量多維度排序和主坐標(biāo)分析。

    表1 大型底棲動物優(yōu)勢種Table 1 The dominant species of macrozoobenthos

    圖2 優(yōu)勢種豐度等級曲線Fig.2 Rank abundance curve of the dominant species

    2.2 優(yōu)勢種功能群類型相對豐度空間分布

    對南黃海大型底棲生物優(yōu)勢種功能群類型的相對豐度空間分布進(jìn)行研究,結(jié)果表明,南黃海運(yùn)動、有頜、表層沉積物取食者相對豐度較高,為27.6%,固著生活、有觸手、表層沉積物取食者相對豐度較低,為1.0%。黃海冷水團(tuán)群落區(qū)域運(yùn)動、無頜、肉食者和運(yùn)動、有頜、肉食者相對豐度較高,分別為32.3%和29.7%?;旌先郝鋮^(qū)域運(yùn)動、有頜、表層沉積物取食底泥者和運(yùn)動、無頜、穴居者相對豐度較高,分別為46%和30%。廣溫群落區(qū)域運(yùn)動有頜、表層沉積物取食者相對豐度較高,分別為32.8%。長江口群落區(qū)域運(yùn)動、無頜、穴居者和運(yùn)動、無頜、表層沉積物取食者相對豐度較高,分別為43.8%和31.3%。

    2.3 優(yōu)勢種的生態(tài)位

    對南黃海大型底棲動物優(yōu)勢種生態(tài)位進(jìn)行研究,結(jié)果表明(表2)大型底棲動物生態(tài)位寬度(Shannon-Wiener指數(shù))變化范圍為1.24—2.15(生態(tài)位寬度較高且均大于1),薄索足蛤、蜈蚣?xì)W努菲蟲、掌鰓索沙蠶、淺水薩氏真蛇尾、黃海刺梳鱗蟲和寡節(jié)甘吻沙蠶生態(tài)位較寬,分別為2.15、2.11、2.07、1.92、1.91和1.91;擬特須蟲、細(xì)弱吻沙蠶和背蚓蟲生態(tài)位較窄,分別為,1.24、1.33和1.35。大型底棲動物優(yōu)勢種的平均邊緣指數(shù)變化范圍為0.23—4.95,擬特須蟲和細(xì)弱吻沙蠶平均邊緣指數(shù)值較高,分別為4.95和4.78;掌鰓索沙蠶、異足科索沙蠶、長葉尖索沙蠶和副櫛蟲平均邊緣指數(shù)較低,分別為0.23、0.38、0.41和0.49。耐受指數(shù)變化范圍0.13—3.85,深鉤毛蟲、寡鰓齒吻沙蠶和黃海刺梳鱗蟲耐受指數(shù)值較高,分別為3.85、2.46和2.25;細(xì)弱吻沙蠶和副櫛蟲耐受指數(shù)較低,分別為0.13和0.41。OMI分析表明,沿給定環(huán)境梯度長吻沙蠶、細(xì)弱吻沙蠶、背蚓蟲、角海蛹和淺水薩氏真蛇尾生態(tài)位分化較為顯著(P<0.05)。OMI分析表明(圖3),第一和第二環(huán)境排序軸分別解釋了優(yōu)勢種生態(tài)位與環(huán)境因子間相互關(guān)系的變異程度,解釋的變異量分別為75.8%和14.7%。

    表2 南黃海大型底棲動物優(yōu)勢種功能群、生活型、生態(tài)幅、OMI指數(shù)和耐受指數(shù)

    Table 2 The functional feeding groups, habit/behavior groups, niche breadth, OMI index, Tolerance index of dominant species from southern Yellow Sea

    物種Species編號No.功能群類型FFGandHBG生態(tài)位寬度NichebreadthOMI指數(shù)OMI耐受指數(shù)TOLP異足科索沙蠶S1SMJ1.370.381.940.807長葉尖索沙蠶S2SMJ1.580.411.920.65掌鰓索沙蠶S3SMJ2.070.231.480.47蜈蚣?xì)W努菲蟲S4CMJ2.110.751.870.18黃海刺梳鱗蟲S5CMJ1.910.782.250.297長吻沙蠶S6CDJ1.892.810.550.006寡節(jié)甘吻沙蠶S7CDJ1.911.390.850.112寡鰓齒吻沙蠶S8CMJ1.522.522.460.07細(xì)弱吻沙蠶S9CDJ1.354.780.130.01擬特須蟲S10BMX1.244.950.600.054深鉤毛蟲S11CMJ1.391.473.850.362獨(dú)指蟲S12SMX1.561.031.500.528背蚓蟲S13SMX1.333.350.690.042角海蛹S14BMX1.683.391.460.005尖葉長手沙蠶S15SDT1.662.180.710.055副櫛蟲S16SST1.390.490.410.81淺水薩氏真蛇尾S17CMX1.923.111.920.001日本梯形蛤S18FDP1.591.151.760.269日本胡桃蛤S19FDP1.681.520.700.175薄索足蛤S20FDP2.150.690.870.199

    FFG:Functional feeding group;HBG:生境/行為類群Habit/Behavior Groups;SMJ: 運(yùn)動、有頜、表層沉積物取食者 Surface deposit-feeding, Motile, Jawed;CMJ: 運(yùn)動、有頜、肉食者Carnivore, Motile, Jawed;CDJ: 半運(yùn)動、有頜、肉食者Carnivore, Discretely Motile, Jawed;SMX: 運(yùn)動、無頜、表層沉積物取食者Surface deposit-feeding, Motile, Non-jawed;BMX: 運(yùn)動、無頜、穴居者Burrowing, Motile, Non jawed;SDT: 半運(yùn)動、有觸手、表層沉積物取食者Surface deposit-feeding, Discretely motile, Tentaculate;SST:固著生活、有觸手、表層沉積物取食者 Surface deposit-feeding, Sessile, Tentaculate;CMX: 運(yùn)動、無頜、肉食者Carnivore, Motile, Non jawed; FDP: 半運(yùn)動、泵吸、濾食者Filter-feeding, Discretely motile, Pumping; 字體加粗部分P值小于0.05

    表3為環(huán)境因子與OMI第一和第二排序軸的相關(guān)系數(shù),表4為各優(yōu)勢種在OMI分析第一和第二排序軸上的物種分?jǐn)?shù),圖4為各優(yōu)勢種在OMI分析第一和第二排序軸上的分布。

    圖3 優(yōu)勢種沿OMI分析環(huán)境梯度的分布Fig.3 The distribution of dominant species on environmental gradients using the OMI analysis

    表3 環(huán)境因子與OMI排序軸的相關(guān)關(guān)系Table 3 Correlation coefficients for the first two axes and 6 environmental variables of OMI analysis

    2.4 優(yōu)勢種的重疊值

    優(yōu)勢種種群生態(tài)位重疊值的范圍在0—0.95之間(表5),結(jié)果表明,長葉尖索沙蠶和掌鰓索沙蠶的生態(tài)位重疊值最大,為0.95。異足科索沙蠶與寡鰓齒吻沙蠶的重疊值為0.70,長葉尖索沙蠶、掌鰓索沙蠶分布與寡節(jié)甘吻沙蠶的重疊值分別為0.82和0.89,長吻沙蠶與背蚓蟲、細(xì)弱吻沙蠶和擬特須蟲的重疊值分別為0.87、0.77和0.67,細(xì)弱吻沙蠶和尖葉長手沙蠶的重疊值為0.69,擬特須蟲和深鉤毛蟲的重疊值為0.64,日本胡桃蛤和薄索足蛤的重疊值為0.71。

    表4 優(yōu)勢種在OMI分析第一和第二軸的物種分?jǐn)?shù)Table 4 11—5 Species scores from both the first and the second axes performed by OMI analysis

    表5 20個(gè)優(yōu)勢種群間優(yōu)勢種重疊值Table 5 Niche overlaps among dominant species

    黑體為生態(tài)位重疊值大于0.6者

    2.5 物種等級聚類和非度量多維度排序

    聚類、非度量多維度排序和主坐標(biāo)分析(圖5)的結(jié)果表明,20個(gè)優(yōu)勢種分為4組,第1組包括長葉尖索沙蠶和掌鰓索沙蠶,第2組包括淺水薩氏真蛇尾,第3組包括角海蛹,第4組包括其他優(yōu)勢種。非度量多維度排序分析的脅強(qiáng)系數(shù)(Stress)為0.12。聚類分析、非度量多維度排序和主坐標(biāo)分析結(jié)果與OMI分析優(yōu)勢種空間分布結(jié)果較一致(圖3),較好的反應(yīng)了南黃海大型底棲動物優(yōu)勢種分布、數(shù)量狀況和生境相似性程度。

    圖4 優(yōu)勢種生態(tài)位和生態(tài)位寬度在平均邊緣指數(shù)分析第一和第二軸上的分布Fig.4 Niche position and niche breadth distributed on the first two axes of OMI analysis

    3 討論

    3.1 大型底棲動物優(yōu)勢種的空間分布

    圖5 優(yōu)勢種的聚類和排序Fig.5 Cluster, NMDS and Pco of the dominant species

    優(yōu)勢類群的分布主導(dǎo)底棲動物分布格局[20],優(yōu)勢種分布是物種適應(yīng)生境條件的表現(xiàn)[21]。本研究結(jié)果與劉錄三等[22]和Zhang等[13]對南黃海大型底棲動物優(yōu)勢種的分布研究較一致,相比東海大型底棲動物優(yōu)勢種研究結(jié)果[23]則差異較大,僅有長吻沙蠶、黃海刺梳鱗蟲、背蚓蟲和角海蛹等優(yōu)勢種較一致。生境在決定南黃海大型底棲動物優(yōu)勢種的空間分布上具有重要作用:黃海冷水團(tuán)區(qū)域大型底棲動物優(yōu)勢種主要為冷水種,淺水薩氏真蛇尾、日本梯形蛤和蜈蚣?xì)W努菲蟲分布較廣;混合群落位于冷水團(tuán)邊緣區(qū)域,薄索足蛤、黃海刺梳鱗蟲、掌鰓索沙蠶、長葉尖索沙蠶和角海蛹分布較廣,大型底棲動物優(yōu)勢種呈現(xiàn)冷水種和暖水種共同出現(xiàn)的特點(diǎn);廣溫性群落位于南黃海近岸海域以及蘇北淺灘海域,長葉尖索沙蠶、掌鰓索沙蠶、長吻沙蠶、寡節(jié)甘吻沙蠶和尖手長葉沙蠶分布較廣;長江口群落區(qū)域擬特須蟲和背蚓蟲分布較廣,分布主要受長江徑流及底質(zhì)有機(jī)質(zhì)含量豐富的影響。此外,混合群落區(qū)域位于生態(tài)交錯(cuò)區(qū),大型底棲動物優(yōu)勢種功能群類型多樣,物種豐富,體現(xiàn)了邊緣效應(yīng),反映了生物對環(huán)境的適應(yīng)性。例如,薄索足蛤的分布與黃海冷水團(tuán)范圍變動密切相關(guān)[24]。除了棲息地的分化,食性分化也是生態(tài)位分化的一種表現(xiàn)形式[25]。功能攝食類群是物種食物關(guān)系的結(jié)合[26],能較好反映物種集群食性的分化。功能攝食類群分布與餌料豐富程度以及棲息地環(huán)境狀況密切相關(guān)[27]。南黃海大型底棲動物優(yōu)勢種功能群以肉食者和食底泥者為主,黃海冷水團(tuán)區(qū)域優(yōu)勢種功能群以肉食者為主,廣溫群落區(qū)域、混合群落區(qū)域和長江口區(qū)域底棲動物優(yōu)勢種功能群以食底泥者為主,攝食器官和運(yùn)動方式多樣,優(yōu)勢種食性和攝食行為空間分化明顯。

    3.2 優(yōu)勢種生態(tài)位、生態(tài)位寬度與生態(tài)位重疊值

    生態(tài)位假說[28]認(rèn)為一定區(qū)域內(nèi),可以利用大多數(shù)該區(qū)域內(nèi)環(huán)境和資源的物種比只能利用邊緣資源和生境的物種分布的更廣。南黃海大型底棲動物優(yōu)勢種擬特須蟲和細(xì)弱吻沙蠶OMI值較高(分別為4.95和4.78),為邊緣物種,各出現(xiàn)在4個(gè)環(huán)境位點(diǎn)中;掌鰓索沙蠶、異足科索沙蠶和長葉尖索沙蠶OMI值較低(分別為0.23、0.38和0.41)為非邊緣物種,分別出現(xiàn)在17、7和9個(gè)環(huán)境位點(diǎn)中。這一結(jié)果表明非邊緣物種比邊緣種分布要廣。此外,典型生境物種也比邊緣種分布更廣,豐度更高。例如,受黃海冷水團(tuán)區(qū)域低溫高影響,典型生境優(yōu)勢種淺水薩氏真蛇尾比邊緣種分布更廣,出現(xiàn)在黃海冷水團(tuán)區(qū)域8個(gè)環(huán)境位點(diǎn)中,豐度高值達(dá)100 個(gè)/m2,能夠形成以該種為主的冷水性物種群落;OMI分析也顯示資源可利用性對該種豐度影響顯著(P=0.001)。而黃海冷水團(tuán)區(qū)域其他優(yōu)勢種,例如蜈蚣?xì)W努菲蟲和日本梯形蛤生態(tài)位較小(OMI指數(shù)值較低)可利用的生境較多,但是局部豐度不高,這可能與這兩個(gè)物種與周圍物種的競爭有關(guān)。

    生態(tài)位寬度是生物利用資源多樣性的一個(gè)指標(biāo)。在現(xiàn)有資源譜中,僅能利用其一小部分的生物,稱為狹生態(tài)位;而利用其大部分的,則稱為廣生態(tài)位[29]。生態(tài)位寬度大小反映物種的分布和數(shù)量[30]。南黃海大型底棲動物優(yōu)勢種薄索足蛤、蜈蚣?xì)W努菲蟲、掌鰓索沙蠶、淺水薩氏真蛇尾、黃海刺梳鱗蟲和寡節(jié)甘吻沙蠶生態(tài)位較寬(Shannon-Wiener >1.9),其中黃海刺梳鱗蟲、蜈蚣?xì)W努菲蟲和寡節(jié)甘吻沙蠶總豐度(40 個(gè)/m2、70 個(gè)/m2和60 個(gè)/m2)較其他物種如尖手長葉沙蠶(110 個(gè)/m2)和角海蛹(225 個(gè)/m2)低,甚至比生態(tài)位較窄的背蚓蟲(80 個(gè)/m2)低,這一現(xiàn)象與物種分布位點(diǎn)及其豐度配置格局有關(guān),因?yàn)閿?shù)量較大的物種如尖手長葉沙蠶分布位點(diǎn)較少,為5個(gè),最終生態(tài)位寬度也較??;角海蛹分布位點(diǎn)較多,為10個(gè),生態(tài)位寬度卻低于黃海刺梳鱗蟲等優(yōu)勢種,這與除C3站(110 個(gè)/m2)以外其他位點(diǎn)角海蛹的豐度低于黃海刺梳鱗蟲等物種有關(guān),優(yōu)勢種生態(tài)位寬度與分布位點(diǎn)密切相關(guān),該結(jié)果與胡知淵等[11]和焦海峰等[12]的研究結(jié)果較一致。Shannon-Wiener指數(shù)表征的物種生態(tài)位寬度側(cè)重于物種和資源位點(diǎn)的數(shù)量,而耐受指數(shù)表征的物種生態(tài)位寬度側(cè)重于環(huán)境資源量對物種的影響[7],能夠較好描述物種生態(tài)位特征與分布格局之間關(guān)系[31],反映物種對環(huán)境的適應(yīng)狀況。本研究結(jié)果表明,南黃海大型底棲動物優(yōu)勢種深鉤毛蟲、寡鰓齒吻沙蠶和黃海刺梳鱗蟲的耐受指數(shù)值較高,生態(tài)位較寬,特化程度較小,對生態(tài)因子耐受范圍較大。且上述物種均為運(yùn)動、有頜、肉食者,運(yùn)動和攝食能力強(qiáng)于耐受指數(shù)較低的半運(yùn)動、有頜、肉食者細(xì)弱吻沙蠶和著生活、有觸手、表層沉積物取食者副櫛蟲,進(jìn)一步印證了功能群類型影響種類的生態(tài)位寬度。例如,Person報(bào)道[32]白毛蟲科Litocorsastremma(運(yùn)動、有頜、肉食者)有較強(qiáng)運(yùn)動攝食能力,生態(tài)位較寬,無需棲息地的分化消除種間競爭影響。

    物種生態(tài)位寬度與生境退化狀況密切相關(guān)[33],耐受性較弱的物種易受生境退化及喪失的影響,而耐受性較強(qiáng)的物種通過尋找合適棲息地來應(yīng)對生境退化的影響[34]。本研究表明耐受性指數(shù)較高的深鉤毛蟲、寡鰓齒吻沙蠶和黃海刺梳鱗蟲生態(tài)位較寬,為廣生性物種。深鉤毛蟲和黃海刺梳鱗蟲在南黃海出現(xiàn)頻率較高[22],寡鰓齒吻沙蠶在乳山近海分布廣泛[35]。上述物種多出現(xiàn)在擾動較強(qiáng)的水域,例如,渤海灣近岸海域[36]和連云港近岸海域[37],受生境退化影響較小。細(xì)弱吻沙蠶和副櫛蟲耐受指數(shù)較低,生態(tài)位較窄,為狹生性物種,生境退化可能對其影響較大。例如,多毛綱雙櫛蟲科Hobsoniaflorida具狹窄的生境生態(tài)位[38],該種可能通過低分散和棲息地特化來減弱生境強(qiáng)烈變化對其種群的影響[39]。

    OMI分析中的隨機(jī)置換顯著性檢驗(yàn)?zāi)軌蜉^好的反映環(huán)境梯度物種生態(tài)位分化狀況[16]。本研究結(jié)果表明,南黃海大型底棲動物優(yōu)勢種長吻沙蠶、細(xì)弱吻沙蠶、背蚓蟲、角海蛹和淺水薩氏真蛇尾生態(tài)位(P<0.05)與環(huán)境因子間的相互關(guān)系較顯著,底質(zhì)狀況、沉積物粒徑、沉積物有機(jī)質(zhì)含量、沉積物總氮含量、水溫和水深是上述物種分布的重要影響因素[40-43]。

    生態(tài)位重疊值與物種分布的環(huán)境位點(diǎn)的重合情況密切相關(guān)[10]。本研究發(fā)現(xiàn)南黃海優(yōu)勢種之間的生態(tài)重疊值差異較大,各優(yōu)勢種的生態(tài)位之間存在著不同程度的重疊,生態(tài)位較寬(耐受指數(shù)較高)的長葉尖索沙蠶、掌鰓索沙蠶和寡節(jié)甘吻沙蠶重疊值較高(重疊值均大于0.8),Wathne等[44]認(rèn)為,生態(tài)位重疊值大于0.6時(shí),種間生態(tài)位重疊較顯著,且這些種出現(xiàn)的環(huán)境位點(diǎn)重合度較高;而這些種與生態(tài)位寬度較小物種的重疊值則較低,最低為0,說明物種占據(jù)的資源位點(diǎn)完全不同,這一結(jié)果與葛寶明[10]和焦海峰等[12]的研究結(jié)果較一致。此外,研究發(fā)現(xiàn)生態(tài)位重疊值較高的長葉尖索沙蠶和掌鰓索沙蠶、日本胡桃蛤和薄索足蛤的功能群類型相同,反映這些物種對資源利用的方式相似,但物種間是否存在資源利用性競爭還要取決于共用資源的供應(yīng)狀況,此外海洋中物種間的共棲、互利和寄生等關(guān)系也會導(dǎo)致生態(tài)位重疊值較高。例如,漁山島潮間帶太平洋側(cè)花??鸄nthopleuranigrescens與疣荔枝螺Thaisluteostoma之間存在共棲關(guān)系,生態(tài)位重疊值較高(0.94)[12]。生態(tài)位重疊值較大的物種可能在食性、生活型等方面存在一定的差異,在群落中有不同的生態(tài)功能和位置,這樣的物種也可以很好的共存[11]。

    3.3 物種聚類和排序

    物種聚類和排序的結(jié)果反映物種分布和數(shù)量,并且聚類排序結(jié)果與生態(tài)位寬度和重疊值有關(guān),如排序(NMDS)得到的第1組(Group 1)的物種生態(tài)位較寬,生態(tài)位重疊值較大。第2組(Group 2)和第3組(Group 3)分別為1個(gè)種,與其他種類距離較遠(yuǎn),表明其在資源位點(diǎn)上分布的差異較大,緩解了物種間對資源的競爭。第4組(Group 4)的物種,由于其分布的環(huán)境位點(diǎn)差異較大,導(dǎo)致物種生態(tài)位寬度和重疊值變化很大,但是功能群和生活型的不同緩解了這些物種間的競爭。非度量多維度排序分析給出的脅強(qiáng)系數(shù)較小,為0.12,表明優(yōu)勢種排序圖與給定資料差異較小,較好的反映優(yōu)勢種分布狀況。主坐標(biāo)分析則直觀反映南黃海不同群落區(qū)域大型底棲動物優(yōu)勢種的在不同資源位點(diǎn)的數(shù)量分布狀況,支持OMI分析和NMDS排序結(jié)果,較好地反映物種對棲息地環(huán)境的適應(yīng)性和物種對環(huán)境資源的利用能力。

    4 展望

    本研究中生態(tài)位寬度和生態(tài)位重疊值是物種水平分布上的一個(gè)測度,OMI分析也揭示了一定環(huán)境梯度下物種生態(tài)位分布。但是隨時(shí)空尺度變化,物種的生態(tài)習(xí)性以及生境狀況會發(fā)生變化,研究結(jié)果也不同,因此準(zhǔn)確反映大型底棲生物種生態(tài)位在環(huán)境梯度上的分布,還需對時(shí)空尺度上物種生態(tài)習(xí)性、環(huán)境因素和一些重要的生態(tài)脅迫因素(如,水母暴發(fā)、赤潮和低氧區(qū)等)等變量與生物現(xiàn)存量數(shù)據(jù)進(jìn)行耦合。

    致謝:本課題組王金寶、董棟、馬林、孫德海、寇琦、隋吉星、甘志彬、黃慧等參加了野外調(diào)查和實(shí)驗(yàn)室樣品分析工作,特此致謝。

    [1] 張金屯. 數(shù)量生態(tài)學(xué). 北京: 科學(xué)出版社, 2004: 120-240.

    [2] Hutchinson G E. Concluding remarks. Cold Spring Harbor Symposia on Quantitative Biology, 1957, 22(2): 415-427.

    [3] Pianka E R. The structure of lizard communities. Annual Review of Ecology and Systematics, 1973, 4(1): 53-74.

    [4] Colwell R K, Futuyma D J. On the measurement of niche breadth and overlap. Ecology, 1971, 52(4): 567-576.

    [5] Feinsinger P, Spears E E, Poole R W. A simple measure of niche breadth. Ecology, 1981, 62(1): 27-32.

    [6] Dolédec S, Chessel D, Gimaret-Carpentier C. Niche separation in community analysis: a new method. Ecology 2000, 81(10): 2914-2927.

    [7] Korsu K, Heino J, Huusko A, Muotka T. Specific niche characteristics facilitate the invasion of an alien fish invader in boreal streams. International Journal of Ecology, 2012: 1-10.

    [8] Gascón S, Machado M, Sala J, da Fonseca L C, Cristo M, Boix D. Spatial characteristics and species niche attributes modulate the response by aquatic passive dispersers to habitat degradation. Marine and Freshwater Research, 2012, 63(3): 232-245.

    [9] Dole-Olivier M J, Malard F, Martin D, Lefébure T, Gibert J. Relationships between environmental variables and groundwater biodiversity at the regional scale. Freshwater Biology, 2009, 54(4): 797-813.

    [10] 葛寶明,鮑毅新,鄭祥,程宏毅. 靈昆島潮間帶大型底棲動物群落結(jié)構(gòu)與生態(tài)位分析. 生態(tài)學(xué)報(bào), 2005, 5(11): 3037-3043.

    [11] 胡知淵, 鮑毅新, 葛寶明, 李歡歡. 圍墾灘涂潮溝秋季大型底棲動物群落和生態(tài)位分析. 動物學(xué)報(bào), 2006, 52(4): 800-809.

    [12] 焦海峰, 施慧雄, 尤仲杰, 黃濱. 漁山島巖礁基質(zhì)潮間帶大型底棲動物優(yōu)勢種生態(tài)位. 生態(tài)學(xué)報(bào), 2011, 31(14): 3928-3936.

    [13] Zhang J L, Xu F S, Liu R Y. Community structure changes of macrobenthos in the South Yellow Sea. Chinese Journal of Oceanology and Limnology, 2012, 30(2): 248-255.

    [14] 國家海洋局. 海洋監(jiān)測規(guī)范第5部分-沉積物分析. 北京: 中國標(biāo)準(zhǔn)出版社, 2007: 50-84.

    [15] Shannon C E, Wiener W. The Mathematical Theory of Communication. Chicago: University of Illinois Press, 1963.

    [16] Reichard M, Jurajda P,imková A, Matějusová I. Size-related habitat use by bitterling (Rhodeussericeus) in a regulated lowland river. Ecology of Freshwater Fish, 2002, 11(2): 112-122.

    [17] Clarke K R, Gerler R N. PRIMER v6: User Manual/Tutorial. Plymouth: PRIMER-E, 2006.

    [18] Fauchald K, Jumars P A. The diet of worms: a study of polychaete feeding guilds. Oceanography and Marine Biology: An Annual Review, 1979, 17: 193-284.

    [19] Bonsdof E, Norlko A, Sandberg E. Structuring zoobenthos: the importance of predation, siphon cropping and physical disturbance. Journal of Experimental Marine Biology and Ecology, 1995, 192(1): 125-144.

    [20] 謝志才, 馬凱, 葉麟, 陳靜, 蔡慶華. 保安湖大型底棲動物結(jié)構(gòu)與分布格局研究. 水生生物學(xué)報(bào), 2007, 31(2): 174-183.

    [21] 周洪華,陳亞寧,李衛(wèi)紅. 塔里木河下游綠洲-荒漠過渡帶植物多樣性特征及優(yōu)勢種群分布格局. 中國沙漠, 2009, 29(4): 688-696.

    [22] 劉錄三, 李新正. 南黃海春秋季大型底棲動物分布現(xiàn)狀. 海洋與湖沼, 2003, 34(1): 26-32.

    [23] 劉錄三, 李新正. 東海春秋季大型底棲動物分布現(xiàn)狀. 生物多樣性, 2002, 10(4): 351-358.

    [24] 曲方圓, 于子山, 劉衛(wèi)霞, 隋吉星, 張志南. 北黃海春季大型底棲生物群落結(jié)構(gòu). 中國海洋大學(xué)學(xué): 自然科學(xué)版, 2009, 39(Suppl.): 109-114.

    [25] 沈國英, 施并章. 海洋生態(tài)學(xué). 北京: 科學(xué)出版社, 2008: 64-183.

    [26] 何志輝. 淡水生態(tài)學(xué). 北京:中國農(nóng)業(yè)出版社出版, 2000.

    [27] Morin J G., Kastendiek J E, Harrington A, Davis N. Organization and patterns of interactions in a subtidal sand community on an exposed coast. Marine Ecology-Progress Series, 1985, 27: 163-185.

    [28] Venier L A, Fahrig L. Habitat availability causes the species abundance-distribution relationship. Oikos, 1996, 76(3): 564-570.

    [29] 趙惠勛, 李俊清, 王鳳有. 群體生態(tài)學(xué). 哈爾濱:東北林業(yè)大學(xué)出版社, 1990: 13-28.

    [30] 葛寶明, 鄭祥, 程宏毅, 鮑毅新, 陳良. 靈昆島圍墾灘涂潮溝大型底棲動物群落和物種生態(tài)位分析.水生生物學(xué)報(bào), 2007, 31(5): 675-681.

    [31] Heino J, Mykr? H, Muotka T. Temporal variability of nestedness and idiosyncratic species in stream insect assemblages. Diversity and Distributions, 2009, 15(2): 198-206.

    [32] Pearson T H.Ltocorsastremmaa new genus and species of pilargid (Polychaeta: Annelida) from the west coast of Scotland, with notes on two other pilargid species. Journal of Natural History. 1970, 4(1): 69-77.

    [33] Bennett J R, Cumming B F, Ginn B K, Smol J P. Broad-scale environmental response and niche conservatism in lacustrine diatom communities. Global Ecology and Biogeography, 2010, 19(5): 724-732.

    [34] Fischer J, Lindenmayer D B. Landscape modification and habitat fragmentation: a synthesis. Global Ecology and Biogeography, 2007, 16(3): 265-280.

    [35] Li X Z, Li B Q, Wang H F, Wang J B, Wang X C, Zhou J, Han Q X, Ma L, Dong C, Zhang B L. Community structure of macrobenthos in coastal water off Rushan, southern Shandong Peninsula,and the relationships with environmental factors. Acta Oceanologica Sinica, 2009, 28(5): 81-93.

    [36] 王瑜, 劉錄三, 劉存歧, 朱延忠, 徐海明. 渤海灣近岸海域春季大型底棲動物群落特征. 環(huán)境科學(xué)研究, 2010, 23(4): 430-436.

    [37] 陳斌林, 方 濤, 李道季. 連云港近岸海域底棲動物群落組成及多樣性特征.華東師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 2007, (2): 1-10.

    [38] Kelly D W, MacIsaac H J, Heath D D. Vicariance and dispersal effects on phylogeographic structure and speciation in a widespread estuarine invertebrate. Evolution, 2006, 60(2): 257—267.

    [39] Olson M A, Zajac R N, Russello M A. Estuarine-Scale Genetic Variation in the PolychaeteHobsoniaflorida(Ampharetidae; Annelida) in Long Island Sound and Relationships to Pleistocene Glaciations. Biological Bulletin, 2009, 217(1): 86-94.

    [40] 王宗興, 韋欽勝, 劉軍, 冉祥斌. 乳山灣外海夏季大型底棲動物分布與環(huán)境因子的典范對應(yīng)分析, 應(yīng)用與環(huán)境生物學(xué)報(bào), 2012, 18(4): 599-604.

    [41] Ferrando A, Méndez N. Effects of organic pollution in the distribution of annelid communities in the Estero de Urías coastal lagoon, Mexico. Scientia Marina, 2011, 75(2): 351-358.

    [42] Hayward P J, Ryland J S. The Marine Fauna of the British Isles and North-West Europe. Volume I: Introduction and Protozoans to Arthropods. Oxford: Clarendon Press, 1990: 627-627.

    [43] 廖玉麟, 肖寧. 中國海棘皮動物的種類組成及區(qū)系特點(diǎn). 生物多樣性, 2011, 19(6): 729-736.

    [44] Wathne J A, Haug T, Lydersen C. Prey preference and niche overlap of ringed sealsPhocahispidaand harp sealsP.groenlandicain the Barents Sea. Marine Ecology Progress Series, 2000, 194: 233-239.

    Niche analysis of dominant species of macrozoobenthic community in the southern Yellow Sea in spring

    PENG Songyao1,2, LI Xinzheng1,*, WANG Hongfa1, ZHANG Baolin1

    1InstituteofOceanology,ChineseAcademyofSciences,Qingdao266071,China2GraduateUniversity,ChineseAcademyofSciences,Beijing100039,China

    A quantitative survey to determine the composition of the macrobenthic communities in the southern Yellow Sea was conducted in April 2011. Over 100 species of macrobenthos were identified. Twenty dominant species, belonging to 17 families, 3 classes and 3 phyla, were determined via rank abundance curve analysis. Surface deposit feeding, motile, jawed feeders (27.6%) were the feeding group with the highest relative abundance values, while surface deposit feeding, sessile, tentaculate feeders were the least abundant (1%). The Yellow Sea Cold Water Mass benthic community was dominated by both carnivorous, motile, non-jawed species and carnivorous, motile, jawed feeders, making up 32.3% and 29.7% of the community composition. In the Mixed community the trophic groups were dominated by both surface deposit feeding, motile, jawed and burrowing, motile, non-jawed feeders (46% and 30%, respectively). Surface deposit feeding, motile, jawed feeders (32.8%) was the most abundant macrobenthos category in the Eurythermal community, while in the Yangtze River Estuary community burrowing, motile, non-jawed and surface deposit-feeding, motile, non-jawed feeders represented 43.8% and 31.3% of the total abundance, respectively. The Shannon-Wiener Niche Breadth index, the Pianka Niche Overlap index, the outlying mean index (OMI) and the tolerance index (TOL) were used to analyze the relationships of the dominant species and environmental factors. Four indices were significantly different among species.Thyasiratokunagai,Onuphisgeophiliformis,Nin?epalmata,Ophiurasarsiivadicola,EhlersileanirahwanghaiensisandGlycindegurjanovaehad large niche breadths. The OMI varied from 0.23 to 4.95.ParalacydoniaparadoxaandGlyceratenuishad higher values of OMI (4.95 and 4.78) than that of any other species and the most marginal niche positions of all the dominant species examined.N.palmata(0.23) andKuwaitaheteropoda(0.38) had relatively non-marginal niche positions, and occurred in average conditions across the sampled sites. The tolerance index ranged from 0.13 to 3.85.Sigambrabassi(3.85),Nephtysoligobranchia(2.46) andE.hwanghaiensis(2.25) were the main contributors to the tolerance index. OMI analysis indicated that two species with lower tolerances (narrower niche breadths; more specialized species) wereG.tenuisandParamphicteisangustifolia. Niche overlaps varied from 0 to 0.95.ScoletomalongifoliaandN.palmatahad a niche overlap value close to 0.95. The significance of the OMI analysis random permutation test demonstrates that niche segregation ofGlycerachirori,G.tenuis,Notomastuslatericeus,OphelinaacuminateandO.sarsiivadicolais effective along a given environmental gradient (e.g., depth, water bottom temperature, salinity, median sediment size, total organic carbon and total nitrogen). Using hierarchical clustering (CLUSTER) with between-group linkage and non-metric multidimensional scaling (NMDS), the 20 dominant species were classified into three groups: generalist species, typical habitat species and specialist species. We also used Principal Coordinates Analysis to visualize the distribution patterns of dominant species in different habitats. Niche breadth and niche overlap of each species were strongly related to functional feeding groups, habit/behaviors of species, spatial distribution, species abundance and habitat conditions. This reveals the differing abilities of species to effectively utilize their environmental resources.

    macrozoobenthos; dominant species; niche; niche breadth; niche overlaps; the southern Yellow Sea

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(973)(2011CB403605); 國家自然科學(xué)基金面上項(xiàng)目(41176133)

    2013-05-31;

    日期:2014-04-25

    10.5846/stxb201305311254

    *通訊作者Corresponding author.E-mail: lixzh@qdio.ac.cn

    彭松耀,李新正,王洪法,張寶琳.南黃海春季大型底棲動物優(yōu)勢種生態(tài)位.生態(tài)學(xué)報(bào),2015,35(6):1917-1928.

    Peng S Y, Li X Z, Wang H F, Zhang B L.Niche analysis of dominant species of macrozoobenthic community in the southern Yellow Sea in spring.Acta Ecologica Sinica,2015,35(6):1917-1928.

    猜你喜歡
    沙蠶黃海生境
    核電廠取水口海域底棲生物起動研究
    水道港口(2023年3期)2023-09-06 09:22:44
    棗樹適應(yīng)干旱生境研究進(jìn)展
    落葉果樹(2021年6期)2021-02-12 01:29:10
    東方濕地 黃海明珠
    解讀沙蠶的養(yǎng)殖方式及其應(yīng)用開發(fā)現(xiàn)狀
    黃海簡介
    沙蠶與水稻共養(yǎng)模式,仿生態(tài)養(yǎng)殖促增產(chǎn)
    三角恒等變換去哪兒了
    金秀瑤族長鼓舞傳承的文化生境
    體育科技(2016年2期)2016-02-28 17:06:07
    不同生境中入侵植物喜旱蓮子草與本地種接骨草的光合特性
    天然和養(yǎng)殖沙蠶營養(yǎng)成分分析
    亚洲黑人精品在线| videosex国产| 久久久久久免费高清国产稀缺| 欧美丝袜亚洲另类 | 国产精品九九99| 国产高清视频在线播放一区| 在线观看66精品国产| 国产成人啪精品午夜网站| 精品免费久久久久久久清纯 | 亚洲精品自拍成人| 国产真人三级小视频在线观看| 国产精品98久久久久久宅男小说| 国产高清国产精品国产三级| 国产精品久久电影中文字幕 | 日韩人妻精品一区2区三区| 深夜精品福利| 伊人久久大香线蕉亚洲五| 757午夜福利合集在线观看| 丰满迷人的少妇在线观看| 视频区图区小说| 欧美精品人与动牲交sv欧美| 19禁男女啪啪无遮挡网站| 在线观看免费视频网站a站| 一级黄色大片毛片| 又黄又爽又免费观看的视频| 欧美日韩国产mv在线观看视频| svipshipincom国产片| 日本黄色日本黄色录像| 精品少妇一区二区三区视频日本电影| 国产成人av激情在线播放| 亚洲av欧美aⅴ国产| 日本精品一区二区三区蜜桃| 女人高潮潮喷娇喘18禁视频| 女性被躁到高潮视频| 国产高清国产精品国产三级| 午夜老司机福利片| 母亲3免费完整高清在线观看| 高清视频免费观看一区二区| 国产男靠女视频免费网站| 乱人伦中国视频| 欧美乱码精品一区二区三区| 两个人看的免费小视频| 日本wwww免费看| 亚洲人成电影免费在线| cao死你这个sao货| 99在线人妻在线中文字幕 | 51午夜福利影视在线观看| 久久中文看片网| videos熟女内射| 国产精品亚洲一级av第二区| 亚洲国产欧美网| 在线播放国产精品三级| 精品一区二区三区四区五区乱码| 午夜福利在线观看吧| 极品教师在线免费播放| 国产不卡av网站在线观看| 精品久久久久久,| 激情视频va一区二区三区| 欧美激情久久久久久爽电影 | 亚洲欧美色中文字幕在线| 免费在线观看影片大全网站| 女性被躁到高潮视频| 啦啦啦视频在线资源免费观看| 高清av免费在线| 国产精品美女特级片免费视频播放器 | www.熟女人妻精品国产| 免费看十八禁软件| 黄色片一级片一级黄色片| 午夜福利在线免费观看网站| 国产91精品成人一区二区三区| 亚洲少妇的诱惑av| 热re99久久精品国产66热6| 午夜福利影视在线免费观看| 人妻丰满熟妇av一区二区三区 | 欧美最黄视频在线播放免费 | 亚洲av日韩精品久久久久久密| 80岁老熟妇乱子伦牲交| 亚洲av成人不卡在线观看播放网| 国产精品免费视频内射| 好看av亚洲va欧美ⅴa在| 中文字幕最新亚洲高清| 中文字幕av电影在线播放| 最近最新中文字幕大全电影3 | 久久国产精品男人的天堂亚洲| 少妇粗大呻吟视频| 变态另类成人亚洲欧美熟女 | 国产又色又爽无遮挡免费看| 久久久久久人人人人人| 久久精品亚洲熟妇少妇任你| 国产精品一区二区免费欧美| 中文字幕色久视频| 成年人黄色毛片网站| 亚洲av日韩在线播放| 国产三级黄色录像| 国产不卡av网站在线观看| 国产一区有黄有色的免费视频| 91字幕亚洲| 亚洲精品在线美女| 亚洲成人免费电影在线观看| 老司机深夜福利视频在线观看| 欧美大码av| 久久影院123| www.精华液| 99久久综合精品五月天人人| 国产成人一区二区三区免费视频网站| 高清黄色对白视频在线免费看| 成年人免费黄色播放视频| 亚洲精品美女久久久久99蜜臀| 午夜福利乱码中文字幕| ponron亚洲| 在线观看www视频免费| 一边摸一边抽搐一进一出视频| 精品国产国语对白av| 一级毛片女人18水好多| 色尼玛亚洲综合影院| 中出人妻视频一区二区| 很黄的视频免费| 欧美日韩一级在线毛片| 国产精品成人在线| 亚洲黑人精品在线| 性色av乱码一区二区三区2| 成人18禁高潮啪啪吃奶动态图| а√天堂www在线а√下载 | 久久草成人影院| 大片电影免费在线观看免费| 99久久国产精品久久久| 操美女的视频在线观看| 久99久视频精品免费| 在线观看日韩欧美| 成年人黄色毛片网站| 国产精品偷伦视频观看了| 中国美女看黄片| 亚洲七黄色美女视频| 国产精品自产拍在线观看55亚洲 | 精品欧美一区二区三区在线| 久久久精品国产亚洲av高清涩受| 亚洲色图 男人天堂 中文字幕| 老司机在亚洲福利影院| 亚洲中文日韩欧美视频| 黄色怎么调成土黄色| 国产aⅴ精品一区二区三区波| 中文字幕最新亚洲高清| av一本久久久久| 精品欧美一区二区三区在线| 操美女的视频在线观看| 国产精品秋霞免费鲁丝片| 国产精品一区二区免费欧美| 欧美另类亚洲清纯唯美| 久久草成人影院| 人妻久久中文字幕网| 777米奇影视久久| 女人被狂操c到高潮| 国产亚洲欧美精品永久| 精品久久久久久久久久免费视频 | 99精品欧美一区二区三区四区| 欧美黑人精品巨大| 极品少妇高潮喷水抽搐| 乱人伦中国视频| av线在线观看网站| 国产欧美日韩一区二区三区在线| 可以免费在线观看a视频的电影网站| 日韩大码丰满熟妇| 后天国语完整版免费观看| 国产区一区二久久| 国产成人免费观看mmmm| 久久久国产成人精品二区 | 国产精华一区二区三区| 99riav亚洲国产免费| 亚洲五月婷婷丁香| 最新的欧美精品一区二区| 国产黄色免费在线视频| 日韩人妻精品一区2区三区| 电影成人av| 久久天躁狠狠躁夜夜2o2o| 大香蕉久久成人网| 亚洲色图 男人天堂 中文字幕| 99香蕉大伊视频| 国产成人免费观看mmmm| 亚洲欧洲精品一区二区精品久久久| 欧美国产精品一级二级三级| 久久国产精品人妻蜜桃| 王馨瑶露胸无遮挡在线观看| 欧美乱码精品一区二区三区| www日本在线高清视频| 97人妻天天添夜夜摸| 777久久人妻少妇嫩草av网站| 国产精品1区2区在线观看. | 亚洲国产看品久久| 国产av一区二区精品久久| 国产高清激情床上av| 久久亚洲精品不卡| 久久天堂一区二区三区四区| 久久中文看片网| 亚洲av成人av| 一夜夜www| 国产精品久久视频播放| 午夜精品久久久久久毛片777| 亚洲人成伊人成综合网2020| 777米奇影视久久| 成熟少妇高潮喷水视频| 女人被狂操c到高潮| 亚洲国产中文字幕在线视频| √禁漫天堂资源中文www| 欧美国产精品一级二级三级| 亚洲欧洲精品一区二区精品久久久| 久久影院123| 脱女人内裤的视频| 丝袜美足系列| 一本综合久久免费| 亚洲国产精品合色在线| 国产精品99久久99久久久不卡| 麻豆成人av在线观看| 身体一侧抽搐| 超碰成人久久| 免费日韩欧美在线观看| 久久ye,这里只有精品| 欧美最黄视频在线播放免费 | 国产伦人伦偷精品视频| 男女下面插进去视频免费观看| 亚洲一码二码三码区别大吗| 欧美黑人精品巨大| 高清视频免费观看一区二区| 久久精品人人爽人人爽视色| 国产一区二区三区在线臀色熟女 | 国产成人精品在线电影| 午夜亚洲福利在线播放| 国产精品 国内视频| 成人手机av| 很黄的视频免费| 亚洲综合色网址| 黄色丝袜av网址大全| 看免费av毛片| 咕卡用的链子| 亚洲午夜精品一区,二区,三区| 国产精品98久久久久久宅男小说| 麻豆成人av在线观看| 12—13女人毛片做爰片一| 亚洲专区字幕在线| 国产精品国产高清国产av | 老司机影院毛片| 高清毛片免费观看视频网站 | 黄色女人牲交| 91麻豆精品激情在线观看国产 | 十八禁高潮呻吟视频| 波多野结衣av一区二区av| 久久久精品免费免费高清| 中文字幕人妻丝袜一区二区| 成人亚洲精品一区在线观看| 人成视频在线观看免费观看| 欧美成狂野欧美在线观看| 国产精品秋霞免费鲁丝片| 亚洲精品国产一区二区精华液| 国产在线精品亚洲第一网站| 91九色精品人成在线观看| 变态另类成人亚洲欧美熟女 | 最近最新免费中文字幕在线| 久久久精品国产亚洲av高清涩受| 99香蕉大伊视频| 国产成人啪精品午夜网站| 免费一级毛片在线播放高清视频 | 男女床上黄色一级片免费看| 色94色欧美一区二区| а√天堂www在线а√下载 | 欧美激情久久久久久爽电影 | 国产精品偷伦视频观看了| 老司机午夜十八禁免费视频| 大片电影免费在线观看免费| 国产成人av激情在线播放| 9191精品国产免费久久| 麻豆成人av在线观看| 亚洲av成人不卡在线观看播放网| av欧美777| 亚洲性夜色夜夜综合| 18在线观看网站| 美女 人体艺术 gogo| 国产精品一区二区免费欧美| 手机成人av网站| 村上凉子中文字幕在线| 亚洲男人天堂网一区| 视频区欧美日本亚洲| 一区二区三区精品91| 性少妇av在线| 久久国产精品人妻蜜桃| 午夜激情av网站| 欧美av亚洲av综合av国产av| 久久久水蜜桃国产精品网| 男女免费视频国产| 91国产中文字幕| 99久久综合精品五月天人人| 欧美最黄视频在线播放免费 | 精品一区二区三区四区五区乱码| 高清在线国产一区| 色播在线永久视频| 成人特级黄色片久久久久久久| 男女午夜视频在线观看| 欧美黑人欧美精品刺激| 一夜夜www| 免费在线观看黄色视频的| 免费久久久久久久精品成人欧美视频| 日日爽夜夜爽网站| 侵犯人妻中文字幕一二三四区| av天堂在线播放| 免费少妇av软件| 天堂√8在线中文| 久久精品91无色码中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 王馨瑶露胸无遮挡在线观看| 在线观看免费视频日本深夜| 丝瓜视频免费看黄片| 丰满的人妻完整版| 日韩熟女老妇一区二区性免费视频| cao死你这个sao货| 99热国产这里只有精品6| 夫妻午夜视频| 免费一级毛片在线播放高清视频 | 国产欧美日韩精品亚洲av| 人成视频在线观看免费观看| 精品国产一区二区久久| 人妻丰满熟妇av一区二区三区 | 国产一区在线观看成人免费| 老汉色av国产亚洲站长工具| av超薄肉色丝袜交足视频| 无人区码免费观看不卡| 免费在线观看日本一区| 一边摸一边抽搐一进一出视频| 日韩免费av在线播放| 欧美老熟妇乱子伦牲交| 大型av网站在线播放| aaaaa片日本免费| 日日爽夜夜爽网站| 日韩中文字幕欧美一区二区| 村上凉子中文字幕在线| 国产av精品麻豆| 桃红色精品国产亚洲av| 国产精品免费大片| 国产精品电影一区二区三区 | 亚洲av日韩在线播放| 免费在线观看亚洲国产| 看免费av毛片| 午夜两性在线视频| 在线观看免费视频网站a站| 欧美激情高清一区二区三区| 在线观看免费午夜福利视频| 国产精品一区二区在线不卡| 欧美激情久久久久久爽电影 | 三级毛片av免费| 黄色a级毛片大全视频| 一级黄色大片毛片| 久久狼人影院| 精品免费久久久久久久清纯 | 欧美 日韩 精品 国产| 欧美老熟妇乱子伦牲交| 老鸭窝网址在线观看| 多毛熟女@视频| 女人被躁到高潮嗷嗷叫费观| 90打野战视频偷拍视频| 超色免费av| 亚洲情色 制服丝袜| 一边摸一边抽搐一进一出视频| 亚洲第一av免费看| 18禁美女被吸乳视频| 欧美黄色淫秽网站| 两性午夜刺激爽爽歪歪视频在线观看 | 精品国产一区二区三区四区第35| 成在线人永久免费视频| 国产欧美日韩综合在线一区二区| 黄色a级毛片大全视频| 亚洲欧美激情在线| 在线永久观看黄色视频| 亚洲aⅴ乱码一区二区在线播放 | 巨乳人妻的诱惑在线观看| 中文字幕另类日韩欧美亚洲嫩草| 女人被狂操c到高潮| 午夜福利影视在线免费观看| 精品久久久久久,| 精品国产乱码久久久久久男人| 夜夜夜夜夜久久久久| 在线观看免费高清a一片| 欧美日韩黄片免| 免费在线观看完整版高清| 欧美乱妇无乱码| 高潮久久久久久久久久久不卡| 热re99久久精品国产66热6| 欧美黄色片欧美黄色片| 男人操女人黄网站| 我的亚洲天堂| 欧美日韩亚洲高清精品| 大陆偷拍与自拍| 欧美一级毛片孕妇| 少妇 在线观看| 欧美黑人欧美精品刺激| 亚洲五月婷婷丁香| 国产在线精品亚洲第一网站| 久久久久久免费高清国产稀缺| 宅男免费午夜| 亚洲人成电影观看| 久久国产精品人妻蜜桃| 日韩中文字幕欧美一区二区| 99热只有精品国产| 如日韩欧美国产精品一区二区三区| 黄色a级毛片大全视频| 国产成人欧美| 国产三级黄色录像| 欧美黄色淫秽网站| 高清欧美精品videossex| 日本欧美视频一区| 18在线观看网站| 啦啦啦在线免费观看视频4| 波多野结衣av一区二区av| 国产精品久久视频播放| 国产精品亚洲av一区麻豆| 丝袜美足系列| 午夜免费成人在线视频| 亚洲av电影在线进入| 午夜精品久久久久久毛片777| www.熟女人妻精品国产| 国产精品一区二区精品视频观看| 99re在线观看精品视频| 99国产综合亚洲精品| 性少妇av在线| 天天躁日日躁夜夜躁夜夜| 午夜免费成人在线视频| 亚洲人成77777在线视频| 欧美色视频一区免费| 国产在线精品亚洲第一网站| 成人特级黄色片久久久久久久| 日韩欧美三级三区| 制服诱惑二区| 丁香六月欧美| 久久久久久人人人人人| 人人妻人人澡人人爽人人夜夜| 99香蕉大伊视频| 免费在线观看日本一区| 侵犯人妻中文字幕一二三四区| 欧美激情 高清一区二区三区| 欧美成狂野欧美在线观看| 亚洲情色 制服丝袜| 国产一区二区三区视频了| 亚洲一区二区三区欧美精品| av一本久久久久| 国产亚洲av高清不卡| 高清视频免费观看一区二区| 亚洲熟女精品中文字幕| 中文字幕av电影在线播放| 亚洲 国产 在线| 免费高清在线观看日韩| 可以免费在线观看a视频的电影网站| 女性被躁到高潮视频| 两个人免费观看高清视频| 亚洲人成电影观看| 热re99久久国产66热| 精品人妻1区二区| 交换朋友夫妻互换小说| 精品一区二区三区av网在线观看| 热re99久久精品国产66热6| 99re6热这里在线精品视频| 久久这里只有精品19| 两性午夜刺激爽爽歪歪视频在线观看 | 搡老岳熟女国产| 18在线观看网站| 777久久人妻少妇嫩草av网站| 久久中文看片网| 99久久精品国产亚洲精品| 12—13女人毛片做爰片一| a级毛片在线看网站| 一级毛片高清免费大全| 亚洲国产毛片av蜜桃av| 色精品久久人妻99蜜桃| 国精品久久久久久国模美| 久久国产乱子伦精品免费另类| 午夜福利在线观看吧| 久久 成人 亚洲| 国产在线观看jvid| 91国产中文字幕| 在线观看舔阴道视频| 少妇被粗大的猛进出69影院| 免费在线观看黄色视频的| 757午夜福利合集在线观看| 99精品欧美一区二区三区四区| 免费不卡黄色视频| 久久国产精品影院| 国产午夜精品久久久久久| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品在线观看二区| 99riav亚洲国产免费| 欧美人与性动交α欧美软件| 亚洲成av片中文字幕在线观看| 亚洲成国产人片在线观看| 亚洲精品国产区一区二| 日韩人妻精品一区2区三区| 亚洲av片天天在线观看| 国产xxxxx性猛交| 国产熟女午夜一区二区三区| 老司机在亚洲福利影院| videosex国产| 国产人伦9x9x在线观看| 国产日韩欧美亚洲二区| 国产色视频综合| 欧美精品啪啪一区二区三区| 啦啦啦 在线观看视频| www.精华液| 午夜免费成人在线视频| 大香蕉久久成人网| 色婷婷久久久亚洲欧美| 亚洲精品国产精品久久久不卡| 多毛熟女@视频| 国产男靠女视频免费网站| 中文欧美无线码| 国产精品国产高清国产av | 国内久久婷婷六月综合欲色啪| 涩涩av久久男人的天堂| 亚洲九九香蕉| 亚洲在线自拍视频| av欧美777| 91麻豆av在线| 韩国精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 91老司机精品| 80岁老熟妇乱子伦牲交| 中文字幕色久视频| 99久久人妻综合| 一级作爱视频免费观看| 一级黄色大片毛片| 亚洲国产精品一区二区三区在线| 日日摸夜夜添夜夜添小说| 在线永久观看黄色视频| 精品欧美一区二区三区在线| 嫩草影视91久久| 一区在线观看完整版| 国产精品亚洲一级av第二区| 69av精品久久久久久| 亚洲av成人一区二区三| 捣出白浆h1v1| 久久ye,这里只有精品| 欧美在线黄色| 操美女的视频在线观看| 欧美日韩亚洲高清精品| 18禁裸乳无遮挡免费网站照片 | 欧美成狂野欧美在线观看| 午夜激情av网站| 性色av乱码一区二区三区2| 男男h啪啪无遮挡| 亚洲av成人av| 国产aⅴ精品一区二区三区波| 一a级毛片在线观看| 首页视频小说图片口味搜索| 亚洲全国av大片| 后天国语完整版免费观看| 又紧又爽又黄一区二区| 亚洲av熟女| 高清欧美精品videossex| 国产成人精品在线电影| 免费在线观看日本一区| 国产精品免费一区二区三区在线 | 搡老熟女国产l中国老女人| 久久香蕉国产精品| 夜夜躁狠狠躁天天躁| 1024香蕉在线观看| 亚洲国产中文字幕在线视频| 涩涩av久久男人的天堂| 人人妻人人澡人人看| 丝袜在线中文字幕| 亚洲国产精品一区二区三区在线| 91九色精品人成在线观看| 国产精品 欧美亚洲| a级毛片在线看网站| 又黄又爽又免费观看的视频| 香蕉丝袜av| 多毛熟女@视频| 国产真人三级小视频在线观看| 亚洲一区二区三区欧美精品| 国产乱人伦免费视频| а√天堂www在线а√下载 | 日韩欧美国产一区二区入口| 亚洲人成电影免费在线| 国产99白浆流出| 亚洲av第一区精品v没综合| 欧美乱色亚洲激情| 极品教师在线免费播放| 91大片在线观看| 大片电影免费在线观看免费| 日本黄色视频三级网站网址 | 欧美日韩亚洲国产一区二区在线观看 | 日韩大码丰满熟妇| 一级毛片精品| 精品人妻在线不人妻| 午夜福利免费观看在线| 亚洲国产中文字幕在线视频| 巨乳人妻的诱惑在线观看| 在线观看免费午夜福利视频| 纯流量卡能插随身wifi吗| 亚洲熟女精品中文字幕| 欧美激情久久久久久爽电影 | 欧美成狂野欧美在线观看| 欧美 亚洲 国产 日韩一| 香蕉丝袜av| 国产精品欧美亚洲77777| 国产av一区二区精品久久| 老司机福利观看| 欧美+亚洲+日韩+国产| 美女 人体艺术 gogo| 女同久久另类99精品国产91| 国产成人影院久久av| 老汉色∧v一级毛片| 高清毛片免费观看视频网站 | 国产淫语在线视频| 欧美精品av麻豆av| 1024香蕉在线观看| 欧美另类亚洲清纯唯美| 亚洲片人在线观看| √禁漫天堂资源中文www| 亚洲美女黄片视频| 亚洲第一欧美日韩一区二区三区| 一区二区三区激情视频| 大陆偷拍与自拍|