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

    雙臺子河口砂質(zhì)潮間帶小型底棲動物的群落結(jié)構(gòu)及其對環(huán)境的響應(yīng)*

    2020-05-25 01:05:26楊純懿張龍輝杜永芬吳亞平張安國袁秀堂
    海洋與湖沼 2020年3期
    關(guān)鍵詞:潮間帶線蟲沉積物

    楊純懿 張龍輝 杜永芬 吳亞平 張安國 袁秀堂

    雙臺子河口砂質(zhì)潮間帶小型底棲動物的群落結(jié)構(gòu)及其對環(huán)境的響應(yīng)*

    楊純懿1張龍輝2杜永芬3①吳亞平1張安國4袁秀堂4

    (1. 河海大學(xué)海洋學(xué)院 南京 210098; 2. 南京大學(xué)地理與海洋科學(xué)學(xué)院 南京 210023; 3. 南京師范大學(xué)海洋科學(xué)與工程學(xué)院 南京 210023; 4. 國家海洋環(huán)境監(jiān)測中心 大連 116023)

    潮間帶是典型的陸海交互區(qū), 其生物群落及其環(huán)境具有典型性特征。海洋小型底棲動物具有極高的群落多樣性和環(huán)境適應(yīng)性, 但是其時空分布特征以及與環(huán)境因子的關(guān)系研究具有區(qū)域局限性。2017年春季(5月)、夏季(8月)和秋季(10月)于我國緯度最高的雙臺子河口潮間帶, 布設(shè)7個斷面, 21個站位, 進行了小型底棲動物和沉積環(huán)境的定量采樣與分析。共檢獲7個小型底棲動物主要類群, 總平均豐度和生物量分別為(3530±1474)ind./10cm2和(1836.5±976.3)μg dwt/10cm2, 其中海洋線蟲為最優(yōu)勢類群, 平均豐度和生物量分別占96.6%和72.1%, 其次是雙殼類和橈足類。豐度/生物量的季節(jié)波動小, 自陸向海略呈降低趨勢(高位≥中位>低位), 統(tǒng)計分析表明其與中值粒徑和和葉綠素含量顯著相關(guān)。小型底棲動物群落受季節(jié)和潮位影響顯著, Bio-env分析顯示相關(guān)性最高的環(huán)境因子組合為沉積物中值粒徑、葉綠素和脫鎂葉綠酸含量。海洋線蟲與橈足的比值(N/C)評估顯示該區(qū)域存在有機污染, 且研究區(qū)西部最嚴重, 與同區(qū)有機污染物的環(huán)境評價結(jié)果較為一致, 說明其在生物監(jiān)測中具有重要意義。

    雙臺子河口; 砂質(zhì)潮間帶; 小型底棲動物; 群落結(jié)構(gòu)

    海洋小型底棲動物(簡稱小型動物)指分選時通過0.5mm孔徑網(wǎng)篩但被0.042mm(深海研究用0.031mm)孔徑網(wǎng)篩截留的后生動物, 主要包括自由生活海洋線蟲、橈足類、渦蟲類、動吻類和多毛類等(Higgins,1988; Giere, 2009)。在底棲生態(tài)系統(tǒng)中, 小型動物與底棲細菌、微藻和原生動物等微型生物組成底棲小/微食物網(wǎng)(杜永芬等, 2011), 其本身又是許多經(jīng)濟魚、蝦和貝類幼體的優(yōu)質(zhì)餌料, 因此, 小型底棲動物在底棲生態(tài)系統(tǒng)的物質(zhì)循環(huán)和能量流動中發(fā)揮重要作用(張志南等, 2017)。另外, 小型動物組成多樣, 在普遍接受的34個動物門類中, 至少有22個門類屬于小型動物類群; 小型動物的優(yōu)勢類群海洋線蟲具有獨特的生殖對策, 生命周期短, 繁殖速度快, 終生底棲, 便于取樣和實驗室培養(yǎng)等特征, 使其成為海洋環(huán)境監(jiān)測、生態(tài)系統(tǒng)演變與健康評估的理想指示類群(Hua, 2016; Du, 2018)。因此海洋小型底棲動物的群落分布特征, 及其在生態(tài)環(huán)境評估的應(yīng)用在基礎(chǔ)研究和應(yīng)用領(lǐng)域一直備受重視(Alves, 2013, 2015)。

    潮間帶是指平均最高潮位和平均最低潮位之間的海岸區(qū)域, 由于地處海陸交界, 水體呈現(xiàn)周期性地漲落, 沉積物規(guī)律性的被淹沒或暴露在空氣中, 環(huán)境(溫度, 鹽度, 沉積物粒度分布等)梯度變化明顯。潮間帶海洋沉積物能有效聚集、儲存有機質(zhì), 為生物的生長繁殖提供重要的物質(zhì)和能量來源, 是碳循環(huán)中重要的源/匯區(qū)(Ridgwell, 2005)。該區(qū)域不僅為生物提供了多樣化的棲居場所, 也是凈化污水, 削弱陸源污染物入海通量的一道屏障(劉海濱等, 2007)。因此潮間帶生物資源特征、生態(tài)系統(tǒng)過程與機理、生態(tài)環(huán)境監(jiān)測與評估等研究, 被列為海岸帶研究的核心內(nèi)容, 其中小型底棲動物的資源動態(tài)和群落結(jié)構(gòu)是重要研究之一(張艷等, 2009; Giere, 2009)。

    雙臺子河口(40°45′—40°51′N, 121°28′—121°44′E)位于渤海遼東灣, 主要由三條河流(雙臺子河, 大遼河和大凌河)沖積作用形成, 面積約6.7×104hm2, 占灘涂總面積的56% (Yang, 2015)。作為我國最北端的河口潮間帶, 每年11月份中下旬進入冰期, 至來年3月結(jié)束。潮流性質(zhì)屬于正規(guī)半日潮流, 潮流基本屬往復(fù)流型。潮間帶較為平緩, 屬典型光灘(吳冠等, 1991), 文蛤等埋棲性貝類資源豐富(Zhang, 2016)。潮上帶被開墾修建了大量的蝦蟹田和海參養(yǎng)殖池, 資源開發(fā)和濫捕使灘涂生態(tài)系統(tǒng)遭到一定程度的破壞。本研究對雙臺子河口潮間帶小型動物群落結(jié)構(gòu)和時空分布, 進行了描述和定量分析, 是我國緯度最高并具有冰期的潮間帶數(shù)據(jù)的首次報道; 同時應(yīng)用小型動物特點對其環(huán)境狀況進行了環(huán)境評估。研究結(jié)果不僅豐富了我國海岸帶研究的基礎(chǔ)資料, 也為環(huán)境監(jiān)測與管理提供了基礎(chǔ)理論數(shù)據(jù)支持。

    1 材料與方法

    1.1 樣品采集

    圖1 研究區(qū)域與站位分布

    分別于2017年5月10—13日、8月21—25日和10月20—24日大潮期間, 在采樣區(qū)設(shè)定的站點, 選取未經(jīng)擾動的沉積物, 用內(nèi)徑為2.9cm的注射器改造的采樣管, 采集0—5cm芯樣8個, 以減少生物斑塊分布的影響, 其中每2個芯樣合并為1個平行樣, 現(xiàn)場用福爾馬林溶液(/, 10%)固定后用于小型底棲動物分析(Somerfield, 1996)。同時, 采集深度為5cm沉積物樣品, 用黑色密封袋包裝后, –20°C冷凍保存, 用于葉綠素(chl)、脫鎂葉綠素(Phae)、粒度、有機質(zhì)及含水量分析。將溫度計插入沉積物下約5cm深處, 現(xiàn)場測定沉積物溫度。以現(xiàn)場挖穴(25cm× 25cm×30cm)滲出的海水采集水樣, 低溫靜置24h后, 常溫下使用手持鹽度計(MASTER-S/Mill α, ATAGO)測定鹽度, 重復(fù)3次, 取平均值作為沉積物鹽度。

    1.2 實驗室分析

    小型底棲動物的提取采用Ludox密度梯度離心法(Giere, 2009; Liu, 2015)。據(jù)Somerfield等(1996)采樣原則和分析要求, 室內(nèi)分樣以不低于200條海洋線蟲為原則(低于200條取全樣)。對研究區(qū)單個野外平行樣抽樣, 分選獲得的海洋線蟲的數(shù)量為3000—6000條, 因此統(tǒng)一取分樣10%用于分析。樣品用福爾馬林溶液濃度(5%)定容到200mL, 顛倒混合均勻后取 20mL分樣, 經(jīng)由500μm和31μm的孔徑的套篩沖洗過濾后, 用LudoxTM硅膠溶液(1.15g/mL)將31μm孔徑網(wǎng)篩上的滯留物全部轉(zhuǎn)移至離心管中(硅膠液與滯留物的體積比例大于3︰1), 低速離心(1800r/min, 10min); 離心后的上清液移出至燒杯, 再加入3倍于沉積物的Ludox液, 與沉積物充分混合后離心, 取上清液; 重復(fù)離心提取3次。將3次離心獲取的上清液合并, 31μm孔徑網(wǎng)篩過濾, 用自來水反復(fù)沖洗, 最后將滯留物全部轉(zhuǎn)移至培養(yǎng)皿中, 在解剖鏡下按類群分選、鑒定和計數(shù)(Du, 2009)。

    粒度分析采用馬爾文激光粒度儀(Mastersizer 3000)進行, 取2—3g沉積物樣品, 添加分散劑分散完全后上機測定。chl和Phae用分光光度法測定, 取3—5g冷凍的沉積物樣品至50mL離心管, 加入10mL濃度90%的丙酮溶液, 低溫避光萃取24h后離心(4000r/min, 15min)后, 取上清液并用90%丙酮溶液定容后, 用分光光度儀(UV1750, SHIMADEZU)測定計算chl和Phea含量(Environmental Sciences Section Inorganic Chemistry Unit, 1991), 其中Phae采用王榮(1986)修正的方法。含水量(WC)以沉積物烘干(60°C, 72h)前后質(zhì)量的減少來計算。有機質(zhì)含量(TOM)采用馬弗爐燒失量法。沉積物樣品先置于60°C烘箱中烘干至恒重(約72h), 然后以烘干樣品在馬弗爐灼燒(550°C)2h后的燒失量, 確定樣品的有機質(zhì)含量(Heiri, 2001; 錢寶等, 2011)。

    1.3 數(shù)據(jù)統(tǒng)計與分析

    小型動物豐度換算為每10cm2的個體數(shù)(ind./10cm2); 生物量以各個類群的豐度, 乘以對應(yīng)類群的平均個體干重獲得(郭玉清等, 2002; 張志南等, 2004)。

    統(tǒng)計分析包括單變量和多變量分析。單變量包括各環(huán)境因子和豐度/生物量。豐度/生物量與環(huán)境因子之間的相關(guān)性采用Spearman相關(guān)分析。季節(jié)、斷面及位置對單變量的影響采用3因子ANOVA差異性分析。以上分析采用SPSS23.0軟件包, 對不滿足方差齊性的數(shù)據(jù)進行雙平方根變換后進行。

    季節(jié)、斷面及位置對小型動物群落的影響, 采用3因子多變量PERMANOVA分析, 對數(shù)據(jù)進行4999次重復(fù)取樣(<0.05為顯著;<0.001為極顯著); 采用主成分分析(PCA)區(qū)分沉積環(huán)境差異及其主要成分; 非度量多維標度(MDS)進行群落相似性分析(生物矩陣經(jīng)雙平方根轉(zhuǎn)化); 生物-環(huán)境分析(Bio-env), 分析影響生物群落分布的最佳環(huán)境要素組合(周紅等, 2003)。以上分析采用PRIMER 7.0軟件進行。

    2 結(jié)果

    2.1 環(huán)境因子

    采樣期間氣溫波動范圍分別為(15.2±1.0)°C (5月)、(24.2±1.3)°C (8月)和(8.9±3.2)°C (10月), 分別代表了春季、夏季和秋季。間隙水鹽度波動為26.0—49.0最高和最低值均出現(xiàn)在A斷面; 春季最高(A1), 秋季最低(A1, A3), 季節(jié)波動較大(=0.008)。鹽度的季節(jié)波動發(fā)生在A斷面, 其顯著差異導(dǎo)致整個區(qū)域的季節(jié)波動; 春季(36.48±2.35)>夏季(34.67±2.85)>秋季(32.05±4.59)。每個季節(jié)高、中和低位站位間無顯著差異(表1)。

    沉積物粉砂和砂的含量較高, 黏土含量極低, 屬于砂質(zhì)粉砂(S-St)、粉砂質(zhì)砂(St-S)和砂(S)(表1)。采樣區(qū)中值粒徑范圍為55.3—149.9μm, 季節(jié)波動小。自高位向低位平均中值粒徑和砂含量呈增大趨勢(中值粒徑, 高位<中位≤低位; 砂含量, 高≤中≤低(春季), 高<中≤低(夏季); 高≤中<低(秋季),<0.05)。沉積物含水量和有機質(zhì)含量均比較穩(wěn)定, 分別為20.80%±1.07%和1.55%±0.23%, 季節(jié)和空間變化均較小。其中有機質(zhì)含量從高位到低位略呈降低趨勢(高≥中≥低)(圖2)。

    表1 不同潮位間表層沉積物特征

    Tab.1 Characteristics of surface sediment in different tidal levels

    注: 字母(a, b, c)表示不同潮位間環(huán)境因子的ANOVA分析結(jié)果, 標注字母相同(a vs a)表示彼此間差異不顯著, 不同字母表示差異顯著(a vs b vs c)表示彼此差異顯著

    Chl含量為0.02—0.93mg/kg。最大值出現(xiàn)在春季高站位(F3), 最小值出現(xiàn)在秋季低站位(B5)。季節(jié)波動不明顯; 空間分布上, 中位的季節(jié)波動較大, 高位和低位季節(jié)波動小, 呈現(xiàn)高位≥中位≥低位趨勢。脫鎂葉綠酸含量為0.06—1.56mg/kg, 高值出現(xiàn)在A—B斷面低站位1.45—1.46mg/kg, 低值0.06— 0.08mg/kg出現(xiàn)在D—G斷面的低站位, 平均含量為(0.39± 0.31)mg/kg, 季節(jié)間及空間變異不顯著(圖2)。

    圖2 雙臺子河口砂質(zhì)潮間帶表層沉積物中有機質(zhì)、葉綠素a (chl a)和脫鎂葉綠酸(Phae)分布

    注: 柱狀圖上的字母(a, b)表示不同潮位間環(huán)境因子的ANOVA分析結(jié)果, 標注字母相同(a vs a)表示彼此間差異不顯著, 不同字母表示差異顯著(a vs b)表示彼此差異顯著

    綜合各環(huán)境因子進行PCA分析結(jié)果見圖3。PC1軸解釋環(huán)境變量的35.2%, PC1和PC2累計解釋環(huán)境變量的61.4%。PC1相關(guān)因子較高的因子是中值粒徑(0.487)、黏土含量(–0.480)、脫鎂葉綠素酸(–0.468)和有機質(zhì)0 (–0.437); PC2相關(guān)因子較高的因子是砂含量(0.627)和溫度(–0.543)。沉積物環(huán)境主要受沉積物粒度組成和季節(jié)影響較大。

    圖3 雙臺子河口潮間帶環(huán)境因子PCA分析結(jié)果

    注: Mgs. 中值粒徑; Sand. 砂含量; Silt. 粉砂含量; WC. 含水量; TOM. 有機質(zhì)含量; chl. 葉綠素含量; Phae. 脫鎂葉綠酸含量; Salinity. 鹽度; T. 溫度

    2.2 小型底棲動物的群落特征

    共鑒定7個小型底棲動物類群(表2), 其中自由生活海洋線蟲(以下簡稱海洋線蟲)為最優(yōu)勢類群, 占小型底棲動物總豐度的96.57%, 其次是雙殼類(1.47%)和橈足類(1.28%)。小型底棲動物的總平均豐度為(3530±1474)ind./10cm2, 最高值出現(xiàn)在春季的高站位(C1), 高達(8293±357)ind./10cm2, 最低值在春季的低站位(E9), 豐度為(1368±419)ind./10cm2(圖4)。三因子分析表明, 豐度分布只有受潮位影響, 呈高位[(4816±1425)ind./10cm2]>中位[(3384±1130)ind./10cm2] >低位[(2391±553)ind./10cm2]; 季節(jié)和采樣斷面對豐度分布無顯著影響, 季節(jié)波動不明顯[夏季(3871± 1480)ind./10cm2≥秋季(3532±1261)ind./10cm2≥春季(3189±1647)ind./10cm2]。海洋線蟲作為豐度上的優(yōu)勢類群, 與小型底棲動物的分布模式基本一致。

    小型底棲動物的平均生物量為(1836.5±976.3)μg dwt/10cm2, 最高值出現(xiàn)在夏季F3站位(7157.6μg dwt/10cm2), 最低值出現(xiàn)在春季E9站位(547.3μg dwt/10cm2)(圖5)。海洋線蟲在生物量上仍占優(yōu)勢, 約占總生物量的72.14%, 其次是雙殼類(10.45%)。時空分布模式與豐度相似, 季節(jié)波動同樣不明顯, 夏季(2523.8±1516.9)μg dwt/10cm2≥(2004.9±948.2)μg dwt/10cm2≥春季(1541.3±947.7)μg dwt/10cm2; 斷面間無差異; 潮位間差異顯著, 高位(2714.8±1506.6)μg dwt/10cm2≥中位(2120.8±1051.6)μg dwt/10cm2≥低位(1326.2±662.4)μg dwt/10cm2。

    表2 雙臺子河口砂質(zhì)潮間帶小型底棲動物主要類群的豐度和生物量

    Tab.2 Abundance and biomass of main meiofaunal taxa in the sandy intertidal zone of the Shuangtaizi estuary

    小型底棲動物主要類群的豐度矩陣進行非度量多維標度(MDS)發(fā)現(xiàn), 各站位的小型底棲動物群落在季節(jié)上可以區(qū)分: 春季和夏季, 高潮位與中低潮位的群落可以區(qū)分, 秋季, 高中潮位與低潮位可以區(qū)分(圖6)。

    2.3 生物與環(huán)境相關(guān)分析

    Spearman相關(guān)性分析表明, 小型底棲動物的豐度和生物量與沉積物中值粒徑負相關(guān), 與黏土含量、粉砂含量和chl含量正相關(guān)(表3)。此外, 豐度與砂含量、有機質(zhì)及含水量顯著相關(guān)。其中沉積物組分在3個季節(jié)均顯著相關(guān), 而chl對豐度和生物量的影響主要發(fā)生在夏季。

    對小型底棲動物群落的PERMANOVA分析結(jié)果顯示(表4), 季節(jié)和潮位對小型底棲動物結(jié)構(gòu)具有顯著性影響(=0.001), 采樣斷面對群落結(jié)構(gòu)無顯著性影響(=0.707); 季節(jié)與潮區(qū), 潮區(qū)與斷面, 季節(jié)×斷面×潮區(qū)的交互作用顯著, 而季節(jié)×斷面交互作用不顯著。不同季節(jié)間和潮區(qū)間的群落組成均存在極顯著差異(=0.001)。

    圖4 雙臺子河口砂質(zhì)潮間帶小型底棲動物的豐度分布

    3 討論

    砂質(zhì)潮灘小型底棲動物的緯度分布研究表明, 其豐度和多樣性在溫帶區(qū)域最高, 極地區(qū)域最低; 其多樣性從極地到熱帶逐漸降低(Lee and Riveros, 2012)。我國海岸帶跨越溫帶、亞熱帶和熱帶, 其中對潮間帶小型底棲動物的研究最北端為僅見于北緯39°大連夏家河子(Hua, 2016)。因此, 研究資料的區(qū)域局限性制約了大尺度的生物多樣性分布格局分析。本研究是我國高緯度潮間帶小型底棲動物分布資料的重要補充。

    圖5 采樣區(qū)域小型動物生物量分布

    3.1 小型底棲動物在我國高緯度砂質(zhì)潮間帶的分布特點

    研究區(qū)小型底棲動物的豐度與我國濱海潮灘處于同一數(shù)量級, 與大連夏家河子的數(shù)值較為接近(Hua, 2016), 但呈現(xiàn)了更高的豐度值。在季節(jié)分布上, 與溫帶青島灣砂質(zhì)潮間帶分布相似(杜永芬等, 2011), 呈現(xiàn)夏季高春季和秋季低的趨勢; 不同于青島灣顯著的季節(jié)差異, 雙臺子河口區(qū)季節(jié)間無顯著差異。這表明高緯度影響削弱了豐度分布的季節(jié)影響。此外, 本次研究沒有涉及冰期冬季的樣品, 希望在后續(xù)研究中加以補充討論。

    圖6 雙臺子河口砂質(zhì)潮間帶小型底棲動物的非度量多維標度(MDS)

    海洋線蟲的優(yōu)勢特征, 與國內(nèi)大部分砂質(zhì)海灘的研究相同(范士亮等, 2006; 杜永芬等, 2011; Sun, 2014)。本研究海洋線蟲的豐度優(yōu)勢高達96.8%, 較大連夏家河高(81.56%), 支持Hua等(2016)研究, 認為海洋線蟲的優(yōu)勢度在我國砂質(zhì)潮間帶隨緯度增高而增高, 在高緯區(qū)具有更高的豐度。不同于已有研究的橈足類和渦蟲居次地位(Sun, 2014; Hua, 2016), 雙殼類在本研究區(qū)居次。這與雙臺子河口潮灘環(huán)境有關(guān), 該區(qū)域貝類自然資源豐富, 也被選為人工蛤類育苗場; 文蛤、四角蛤蜊、光滑河藍蛤等雙殼類是大型底棲動物的優(yōu)勢種(Zhang, 2016, 貝類幼體構(gòu)成了小型動物的重要組分。橈足類喜好高鹽和較為溫暖的環(huán)境(Ingole, 1998; Priyalakshmi, 2014), 河口的淡水輸入和高緯度冰期(11月到次年3月)會削弱了其生存環(huán)境適宜度。渦蟲喜好生活在中度至粗砂和礫石為特征的海灘中, 在泥質(zhì)環(huán)境下渦蟲豐度較低(Delgado, 2009; Kotwicki, 2014)。研究區(qū)域沉積物類型以粉砂質(zhì)砂和砂為主, 不適合渦蟲生長。

    表3 小型動物變量與環(huán)境因子相關(guān)性分析結(jié)果(=63)

    Tab.3 Results of the correlation analysis between meiofaunal variables and environmental factors

    注: **在0.01水平上(雙側(cè))顯著相關(guān); *在0.05水平上(雙側(cè))顯著相關(guān)

    表4 小型動物群落的多變量PERMANOVA分析結(jié)果

    Tab.4 PERMANOVA results for the multivariate descriptors of meiofaunal communities

    注: **在0.01水平上(雙側(cè))顯著相關(guān); *在0.05水平上(雙側(cè))顯著相關(guān)

    3.2 環(huán)境因子對小型底棲動物分布的影響

    生物與環(huán)境協(xié)同進化的結(jié)果使得特定的沉積微環(huán)境生存特定的生物類群。潮灘生物受潮水沖擊擾動的周期影響, 具有一定的微環(huán)境梯度, 沉積物粒徑分布反映了沉積物所處的水文動力狀況, chl和有機質(zhì)含量則代表了生物的食物來源。其中, 沉積物粒徑對潮間帶沉積物微環(huán)境具有主導(dǎo)作用, 尤其是在大尺度的空間分布上, 被認為是重要的物理因子(Giere, 2009)。通常認為125—250μm粒徑的沉積物最有助于微/小型底棲動物的生長(Fenchel, 1969, 1978)。研究區(qū)域的沉積物粒徑為 60—150μm, 與大連河口夏季的粒徑250—500μm比較, 更具有粒度優(yōu)勢, 因此呈現(xiàn)更高的豐度值。

    在小尺度的潮帶分布上, 砂質(zhì)潮灘生物在水平尺度上受沉積物穩(wěn)定性、暴露時間造成的捕食壓力等影響, 在垂直分布上受化學(xué)因子(如氧氣濃度梯度)等影響(Armonies, 2000)。因此, 高潮區(qū)環(huán)境惡劣, 豐度較低(Knox, 2001); 中潮區(qū)是最佳生存環(huán)境(Gheskiere, 2004), 豐度最高, 這里不僅物化環(huán)境相對優(yōu)化, 并且微生物活動豐富, 提供更多的食物來源(Takada, 2015)。同樣地, Hua等(2016)在青島和大連研究結(jié)果顯示, 小型底棲動物在溫帶沙灘中以中潮區(qū)豐度較高。受灘涂圍墾的影響, 本研究中的高位、中位和低位代表了相對海岸線的距離, 從海水淹沒的角度屬于中低潮區(qū), 因此豐度的高位≥中位>低位分布趨勢仍然支持中潮區(qū)豐度較高的分布趨勢。豐度變化趨勢與中值粒徑由高位到低位增大趨勢相反, 其粒徑較小的沉積結(jié)構(gòu)維持更高的豐度與Hua等(2016)結(jié)果一致。

    沉積物粒度組成、有機質(zhì)和Phae含量是影響研究區(qū)沉積環(huán)境區(qū)劃的重要因子(PCA分析); 相關(guān)分析也表明了粒度組成、有機質(zhì)和chl是小型動物豐度和群落構(gòu)成的重要影響因子。有機質(zhì)和chl含量略呈降低趨勢, 與豐度分布趨勢一致。綜合表明溫帶砂質(zhì)潮間帶小型動物的分布受非生物和生物因素的共同影響, 其中粒度構(gòu)成是主要因素。值得注意的是, 潮灘頻繁的人類活動, 如潮上帶圍墾建設(shè)的養(yǎng)殖池, 不但會改變潮流動力過程, 伴隨著污染物的排放, 進而改變沉積環(huán)境, 生態(tài)系統(tǒng)趨于變得簡單化, 生態(tài)問題趨于復(fù)雜化。

    3.3 小型底棲動物對環(huán)境生態(tài)的響應(yīng)與評估

    線蟲與橈足的比值(N/C)被認為海洋富營養(yǎng)化和有機污染程度的指標; N/C大于100意味著環(huán)境狀況為有機污染區(qū)(Raffaelli, 1981)。隨著相關(guān)研究數(shù)據(jù)的增多, 雖然該指標不斷引起爭議, 仍然被廣泛采用(Alves, 2013)。研究區(qū)52.2%—76.2%的采樣站位比值大于100, 即這些站位存在有機污染, 且研究區(qū)最西面的F和G斷面的污染程度最高(表5)。該結(jié)果與Yuan等(2017)關(guān)于16種多環(huán)芳烴的研究結(jié)果較為一致, 潮灘的西側(cè)及近岸區(qū)域的沉積物中有機污染物濃度高, 生態(tài)風險高。表明N/C值作為一重要的生態(tài)指數(shù), 在一定程度上可以為生物監(jiān)測提供依據(jù)。但是, 正如上文討論的結(jié)果, 隨著緯度升高, 海洋線蟲的優(yōu)勢度會增高, 橈足類的生存適宜度降低, 也會影響N/C比值。因此, 在應(yīng)用中需要綜合考慮環(huán)境、指示生物和群落結(jié)構(gòu)等因素進行綜合評估。

    表5 雙臺子河口采樣區(qū)域站位的海洋線蟲/橈足豐度比值(N/C)

    Tab.5 Ratio of nematodes/copepods (N/C) in the Shuangtaizi estuary

    4 結(jié)論

    雙臺子河口潮間帶代表了我國高緯度砂質(zhì)潮間帶, 除溫度具有季節(jié)變化外, 其他環(huán)境因子主要與潮位有關(guān), 中值粒徑由陸向海呈增大趨勢, 葉綠素、有機質(zhì)略呈增高趨勢。沉積物組成, 有機質(zhì)和脫鎂葉綠素含量是沉積微環(huán)境形成的重要因素。

    共檢獲7個小型底棲動物主要類群, 其中海洋線蟲為最優(yōu)勢類群, 雙殼類次之; 小型底棲動物的總平均豐度和生物量分別為3530±1474ind./10cm2和1836.5±976.3μg dwt/10cm2, 豐度和線蟲的優(yōu)勢度均高于我國其他砂質(zhì)潮間帶報道。

    小型動物豐度季節(jié)波動小, 空間上自陸向海豐度值降低(高≥中>低); 小型底棲動物群落受季節(jié)和潮位影響顯著。小型動物的分布受非生物和生物因素的共同影響, 其中粒度構(gòu)成是主要因素。

    海洋線蟲與橈足的比值(N/C)評估顯示該區(qū)域存在有機污染, 且研究區(qū)西部最嚴重, 與同區(qū)有機污染物的環(huán)境評價結(jié)果較為一致, 可為生物監(jiān)測提供重要參考。

    王 榮, 1986. 熒光法測定浮游植物色素計算公式的修正. 海洋科學(xué), 10(3): 1—5

    劉海濱, 張志南, 范士亮等, 2007. 潮間帶小型底棲生物生態(tài)學(xué)研究的某些進展. 中國海洋大學(xué)學(xué)報, 37(5): 767—774

    杜永芬, 徐奎棟, 類彥立等, 2011. 青島灣小型底棲生物周年數(shù)量分布與沉積環(huán)境. 生態(tài)學(xué)報, 31(2): 431—440

    吳 冠, 王錫侯, 劉恒魁, 1991. 遼東灣頂淺海區(qū)海流分布特征. 海洋通報, 10(5): 8—13

    張 艷, 張志南, 華 爾, 2009. 南黃海小型底棲動物分布及其與環(huán)境因子的關(guān)系. 中國農(nóng)學(xué)通報, 25(19): 323—329

    張志南, 林巋旋, 周 紅等, 2004. 東、黃海春秋季小型底棲生物豐度和生物量研究. 生態(tài)學(xué)報, 24(5): 997—1005

    張志南, 周 紅, 華 爾等, 2017. 中國小型底棲生物研究的40年——進展與展望. 海洋與湖沼, 48(4): 657—671

    范士亮, 劉海濱, 張志南等, 2006. 青島太平灣砂質(zhì)潮間帶小型底棲生物豐度和生物量的研究. 中國海洋大學(xué)學(xué)報(自然科學(xué)版), 36(S): 98—104

    周 紅, 張志南, 2003. 大型多元統(tǒng)計軟件PRIMER的方法原理及其在底棲群落生態(tài)學(xué)中的應(yīng)用. 青島海洋大學(xué)學(xué)報, 33(1): 58—64

    錢 寶, 劉 凌, 肖 瀟, 2011. 土壤有機質(zhì)測定方法對比分析. 河海大學(xué)學(xué)報(自然科學(xué)版), 39(1): 34—38

    郭玉清, 張志南, 慕芳紅, 2002. 渤海小型底棲動物生物量的初步研究. 海洋學(xué)報, 24(6): 76—83

    Alves A S, Ad?o H, Ferrero T J, 2013. Benthic meiofauna as indicator of ecological changes in estuarine ecosystems: the use of nematodes in ecological quality assessment. Ecological Indicators, 24: 462—475

    Alves A S, Caetano A, Costa J L, 2015. Estuarine intertidal meiofauna and nematode communities as indicator of ecosystem's recovery following mitigation measures. Ecological Indicators, 54: 184—196

    Armonies W, Reise K, 2000. Faunal diversity across a sandy shore. Marine Ecology Progress Series, 196: 49—57

    Delgado J D, Riera R, Monterroso ó, 2009. Distribution and abundance of meiofauna in intertidal sand substrata around Iceland. Aquatic Ecology, 43(2): 221—233

    Du Y F, Gao S, Liu X S, 2018. Meiofauna and nematode community characteristics indicate ecological changes induced by geomorphic evolution: A case study on tidal creek systems. Ecological Indicators, 87: 97—106

    Du Y F, Xu K D, Lei Y L, 2009. Simultaneous enumeration of diatom, protozoa and meiobenthos from marine sediments using Ludox-QPS method. Chinese Journal of Oceanology and Limnology, 27(4): 775—783

    Environmental Sciences Section Inorganic Chemistry Unit, Wisconsin State Lab of Hygiene, 1991. ESS Method 150.1 Chlorophyll-spectrophotometric. Environmental Sciences Section Inorganic Chemistry Unit, Wisconsin State Lab of Hygiene, 3359—3363

    Fenchel T M, 1978. The ecology of micro-and meiobenthos. Annual Review of Ecology and Systematics, 9: 99—121

    Fenchel T, 1969. The ecology of marine microbenthos IV. Structure and function of the benthic ecosystem, its chemical and physical factors and the microfauna commuities with special reference to the ciliated protozoa. Ophelia, 6(1): 1—182

    Gheskiere T, Hoste E, Vanaverbeke J, 2004. Horizontal zonation patterns and feeding structure of marine nematode assemblages on a macrotidal, ultra-dissipative sandy beach (De Panne, Belgium). Journal of Sea Research, 52(3): 211—226

    Giere O, 2009. Meiobenthology: the Microscopic Motile Fauna of Aquatic Sediments. Berlin: Springer

    Heiri O, Lotter A F, Lemcke G, 2001. Loss on ignition as a method for estimating organic and carbonate content in sediments: reproducibility and comparability of results. Journal of Paleolimnology, 25(1): 101—110

    Higgins R P, Thiel H, 1988. Introduction to the Study of Meiofauna. Washington: Smithsonian

    Hua E, Zhang Z N, Zhou H, 2016. Meiofauna distribution in intertidal sandy beaches along China shoreline (18°–40° N). Journal of Ocean University of China, 15(1): 19—27

    Ingole B S, Parulekar A H, 1998. Role of salinity in structuring the intertidal meiofauna of a tropical estuarine beach: field evidence. Indian Journal of Marine Sciences, 27(3—4): 356—361

    Knox G A, 2001. The Ecology of Seashores. Boca Raton, USA: CRC Press, 557

    Kotwicki L, Deidun A, Grzelak K, 2014. A preliminary comparative assessment of the meiofaunal communities of Maltese pocket sandy beaches. Estuarine, Coastal and Shelf Science, 150: 111—119

    Lee M R, Riveros M, 2012. Latitudinal trends in the species richness of free-living marine nematode assemblages from exposed sandy beaches along the coast of Chile (18–42°S). Marine Ecology, 33(3): 317—325

    Liu X S, Huang D M, Zhu Y M, 2015. Bioassessment of marine sediment quality using meiofaunal assemblages in a semi-enclosed bay. Marine Pollution Bulletin, 100(1): 92—101

    Priyalakshmi G, Menon N R, 2014. Ecology of Interstitial faunal assemblage from the Beaches along the coast of Kerala, India. International Journal of Oceanography, 2014: 284979

    Raffaelli D G, Mason C F, 1981. Pollution monitoring with meiofauna, using the ratio of nematodes to copepods. Marine Pollution Bulletin, 12(5): 158—163

    Ridgwell A, Zeebe R E, 2005. The role of the global carbonate cycle in the regulation and evolution of the Earth system. Earth and Planetary Science Letters, 234(3—4): 299—315

    Somerfield P J, Warwick R M, 1996. Meiofauna in Marine Pollution Monitoring Programmes. A Laboratory Manual. Lowestoft: Ministry of Agriculture, Fisheries and Farming

    Sun X Y, Zhou H, Hua E, 2014. Meiofauna and its sedimentary environment as an integrated indication of anthropogenic disturbance to sandy beach ecosystems. Marine Pollution Bulletin, 88(1—2): 260—267

    Takada Y, Kajihara N, Mochizuki S et al, 2015. Effects of environmental factors on the density of three species of peracarid crustaceans in micro-tidal sandy shores in Japan. Ecological Research, 30(1): 101—109

    Yang X L, Yuan X T, Zhang A G, 2015. Spatial distribution and sources of heavy metals and petroleum hydrocarbon in the sand flats of Shuangtaizi Estuary, Bohai Sea of China. Marine Pollution Bulletin, 95(1): 503—512

    Yuan X T, Yang X L, Zhang A G, 2017. Distribution, potential sources and ecological risks of two persistent organic pollutants in the intertidal sediment at the Shuangtaizi Estuary, Bohai Sea of China. Marine Pollution Bulletin, 114(1): 419—427

    Zhang A G, Yuan X T, Yang X L, 2016. Temporal and spatial distributions of intertidal macrobenthos in the sand flats of the Shuangtaizi Estuary, Bohai Sea in China. Acta Ecologica Sinica, 36(3): 172—179

    MEIOFAUNAL COMMUNITY IN SANDY INTERTIDAL ZONE AND ITS ENVIRONMENTAL RESPONSE IN SHUANGTAIZI ESTUARY, NE CHINA

    YANG Chun-Yi1, ZHANG Long-Hui2, DU Yong-Fen3, WU Ya-Ping1, ZHANG An-Guo4, YUAN Xiu-Tang4

    (1. College of Oceanography, Hohai University, Nanjing 210098, China; 2. School of Geography and Ocean Science, Nanjing University, Nanjing 210023, China; 3. College of Marine Science and Engineering, Nanjing Normal University, Nanjing 210023, China; 4. National Marine Environmental Monitoring Center, Dalian 116023, China)

    Intertidal area is a typical land-sea interaction zone, where the characteristics in biological composition and their micro-environments are of regional representative. Among them, meiofauna has extremely high community diversity and environmental adaptability, while the data of their spatial-temporal distribution patterns and the related regulators are of regional limit. Cores sediments were quantitatively taken in the northernmost-intertidal zone of the Shuangtaizi estuary in the Bohai Sea, NE China in May (spring), August (summer), and October (autumn) in 2017. Twenty-one stations were designed in seven sections, along which the meiofauna and sedimentary environments were analyzed. Seven meiofaunal taxa were found. Free-living marine nematodes were the most dominant group, accounting for 96.6% and 72.1% of the total in abundance and biomass composition, respectively, followed by bivalves and copepods. The average abundance and biomass of meiofauna were (3530±1474)ind./10cm2and (1836.5±976.3)μg dwt/10cm2, respectively, with no seasonal variation and downward trends from land to sea. The medium grain size and chlwere correlated significantly to the abundance/biomass. The ratio of nematodes/copepod (N/C) indicated serious organic pollution, which agreed with the results of the macrofaunal analysis and heavy metal samplings in the similar areas, and could be used as a bio-monitoring indicator.

    Shuangtaizi estuary; sandy intertidal zone; meiofauna; community structure

    * 國家自然科學(xué)基金項目, 41576154號; 海洋公益性行業(yè)科研專項經(jīng)費項目, 201305043號。楊純懿, 碩士研究生, E-mail: ycyhohai@163.com

    杜永芬, 博士, 副教授, E-mail: duyongfen@nju.edu.cn

    2019-10-31,

    2020-03-09

    Q178.4

    10.11693/hyhz20191000202

    猜你喜歡
    潮間帶線蟲沉積物
    走進潮間帶
    夏季蔬菜換茬期線蟲防治要注意
    晚更新世以來南黃海陸架沉積物源分析
    海洋通報(2022年2期)2022-06-30 06:07:04
    渤海油田某FPSO污水艙沉積物的分散處理
    海洋石油(2021年3期)2021-11-05 07:43:12
    水體表層沉積物對磷的吸收及釋放研究進展
    基于海陸融合DEM的大尺度潮間帶提取方法
    福建閩江口潮間帶大型底棲動物次級生產(chǎn)力時空特征
    海洋通報(2020年3期)2020-10-20 09:26:36
    潮間帶4.0MW風機吊裝技術(shù)研究與應(yīng)用
    地黃花對秀麗線蟲壽命的影響
    中成藥(2018年2期)2018-05-09 07:20:04
    朝鮮孢囊線蟲——浙江省孢囊線蟲新記錄種
    一个人看视频在线观看www免费| 日本猛色少妇xxxxx猛交久久| av.在线天堂| 久久久久久久久中文| 熟女电影av网| 欧美成人精品欧美一级黄| 婷婷六月久久综合丁香| 26uuu在线亚洲综合色| 插逼视频在线观看| 成年免费大片在线观看| 国产精品熟女久久久久浪| 成人漫画全彩无遮挡| 色综合亚洲欧美另类图片| 在线天堂最新版资源| 欧美潮喷喷水| 97超视频在线观看视频| 美女大奶头视频| 人人妻人人澡人人爽人人夜夜 | 国产视频内射| 国产69精品久久久久777片| 春色校园在线视频观看| 在线观看美女被高潮喷水网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | .国产精品久久| 全区人妻精品视频| 国产一级毛片七仙女欲春2| 大香蕉久久网| 亚洲精华国产精华液的使用体验| 国产久久久一区二区三区| 免费看av在线观看网站| 色尼玛亚洲综合影院| 亚洲精品中文字幕在线视频 | 一边亲一边摸免费视频| 免费播放大片免费观看视频在线观看| 亚洲av日韩在线播放| 欧美激情久久久久久爽电影| 国产精品99久久久久久久久| 国产亚洲精品av在线| 成年版毛片免费区| 国产成人一区二区在线| 尤物成人国产欧美一区二区三区| 麻豆成人午夜福利视频| 午夜激情欧美在线| 777米奇影视久久| 精品亚洲乱码少妇综合久久| 成人二区视频| 99热全是精品| 一级黄片播放器| 久久精品夜色国产| 国产精品美女特级片免费视频播放器| 天美传媒精品一区二区| 国产av码专区亚洲av| 午夜福利视频精品| av一本久久久久| 国产成人91sexporn| 免费看a级黄色片| av在线观看视频网站免费| av福利片在线观看| 成人欧美大片| 久久久久性生活片| .国产精品久久| 日韩在线高清观看一区二区三区| 国产成人a区在线观看| 91精品国产九色| 三级毛片av免费| 国产成人a∨麻豆精品| 免费播放大片免费观看视频在线观看| 国产成人a∨麻豆精品| 欧美性感艳星| 精品久久久久久久久av| 欧美+日韩+精品| 免费在线观看成人毛片| 国产人妻一区二区三区在| 搡老妇女老女人老熟妇| 国产精品精品国产色婷婷| 成人毛片60女人毛片免费| 特级一级黄色大片| 午夜精品一区二区三区免费看| 中文在线观看免费www的网站| 丰满少妇做爰视频| 毛片女人毛片| 日本欧美国产在线视频| 国产片特级美女逼逼视频| 97热精品久久久久久| 免费黄色在线免费观看| 女人被狂操c到高潮| av在线老鸭窝| 国产淫语在线视频| 亚洲av成人精品一区久久| 成人综合一区亚洲| 日本熟妇午夜| 免费人成在线观看视频色| 国产大屁股一区二区在线视频| 日本一本二区三区精品| 日本三级黄在线观看| 中文字幕久久专区| 国产探花极品一区二区| 国产精品综合久久久久久久免费| 精品久久久久久久久亚洲| 亚洲精品乱久久久久久| 国产黄色小视频在线观看| 插阴视频在线观看视频| 国产中年淑女户外野战色| 熟妇人妻不卡中文字幕| 夜夜爽夜夜爽视频| 国产激情偷乱视频一区二区| 亚洲成色77777| 国产三级在线视频| av又黄又爽大尺度在线免费看| 日本爱情动作片www.在线观看| 午夜免费激情av| 久久久久性生活片| 天堂影院成人在线观看| 午夜免费男女啪啪视频观看| 搞女人的毛片| 欧美3d第一页| 国产精品麻豆人妻色哟哟久久 | av免费观看日本| 久久精品夜夜夜夜夜久久蜜豆| 国产一区二区在线观看日韩| 夫妻午夜视频| 熟女人妻精品中文字幕| 日本三级黄在线观看| 一个人看视频在线观看www免费| 男的添女的下面高潮视频| 日韩欧美精品v在线| 床上黄色一级片| 视频中文字幕在线观看| 亚洲国产最新在线播放| 久久久久精品性色| 精品一区二区三卡| 高清午夜精品一区二区三区| 日本与韩国留学比较| 亚洲怡红院男人天堂| 国国产精品蜜臀av免费| 美女内射精品一级片tv| 国产美女午夜福利| 成年免费大片在线观看| 日本一本二区三区精品| 欧美最新免费一区二区三区| 国产成人精品福利久久| 国产午夜福利久久久久久| 国产精品无大码| 亚洲精品久久午夜乱码| 99热这里只有是精品在线观看| 亚洲三级黄色毛片| 熟女电影av网| 久久这里只有精品中国| 一个人观看的视频www高清免费观看| 熟女人妻精品中文字幕| 国产欧美日韩精品一区二区| 亚洲成人av在线免费| 91av网一区二区| 一级毛片久久久久久久久女| 国产亚洲一区二区精品| 欧美高清成人免费视频www| 久久精品人妻少妇| 国产成人福利小说| 在线 av 中文字幕| 久久久久久久久久久丰满| 免费看光身美女| 亚洲精品日韩在线中文字幕| 黑人高潮一二区| 国产黄色小视频在线观看| 三级毛片av免费| 成人美女网站在线观看视频| 波野结衣二区三区在线| 亚洲四区av| 简卡轻食公司| 又黄又爽又刺激的免费视频.| 黄色日韩在线| 日韩欧美三级三区| 18禁裸乳无遮挡免费网站照片| 精品亚洲乱码少妇综合久久| 欧美 日韩 精品 国产| av黄色大香蕉| 亚洲色图av天堂| 亚洲av成人精品一二三区| 国产成人a∨麻豆精品| 99热这里只有精品一区| 亚洲最大成人中文| 日本色播在线视频| 精品久久久久久久久av| 国产黄色小视频在线观看| 久久久久久九九精品二区国产| 亚洲精品一二三| av在线蜜桃| 久久精品国产自在天天线| 亚洲av福利一区| 噜噜噜噜噜久久久久久91| 舔av片在线| 在线 av 中文字幕| 午夜久久久久精精品| 最近中文字幕2019免费版| 久久久久久伊人网av| 国产在视频线精品| 高清在线视频一区二区三区| kizo精华| av在线观看视频网站免费| 久久综合国产亚洲精品| 在线免费观看不下载黄p国产| 欧美潮喷喷水| 午夜日本视频在线| 狂野欧美激情性xxxx在线观看| 禁无遮挡网站| 777米奇影视久久| 男女下面进入的视频免费午夜| 精品一区在线观看国产| 高清欧美精品videossex| 亚洲国产色片| 高清在线视频一区二区三区| 日本-黄色视频高清免费观看| 国产精品人妻久久久影院| 国产麻豆成人av免费视频| 欧美变态另类bdsm刘玥| 亚洲乱码一区二区免费版| 免费观看性生交大片5| 久久99热6这里只有精品| 日韩成人伦理影院| 亚洲激情五月婷婷啪啪| 99热网站在线观看| 亚洲综合精品二区| 日本一本二区三区精品| 日韩人妻高清精品专区| a级毛片免费高清观看在线播放| 国产精品美女特级片免费视频播放器| 久久久欧美国产精品| 日产精品乱码卡一卡2卡三| 18禁在线播放成人免费| 久久精品久久久久久久性| 免费看不卡的av| 国产午夜精品久久久久久一区二区三区| 免费av毛片视频| 中文字幕av成人在线电影| 欧美最新免费一区二区三区| 蜜臀久久99精品久久宅男| 97在线视频观看| 国产视频内射| 国产伦一二天堂av在线观看| 久久国产乱子免费精品| 久久综合国产亚洲精品| 色视频www国产| 最近最新中文字幕大全电影3| 国产永久视频网站| 成人午夜高清在线视频| 亚洲国产精品国产精品| 亚洲av二区三区四区| 欧美bdsm另类| 亚洲aⅴ乱码一区二区在线播放| 国产精品久久久久久久久免| 国产在视频线在精品| 三级国产精品欧美在线观看| 舔av片在线| 国产久久久一区二区三区| 国产女主播在线喷水免费视频网站 | 人妻制服诱惑在线中文字幕| 国产亚洲5aaaaa淫片| 九色成人免费人妻av| 国产乱来视频区| 白带黄色成豆腐渣| 亚州av有码| 日韩一区二区三区影片| 伦理电影大哥的女人| 一级片'在线观看视频| 韩国av在线不卡| 色5月婷婷丁香| 亚洲欧美成人综合另类久久久| 日韩欧美精品免费久久| 中文字幕免费在线视频6| 亚洲电影在线观看av| 成人特级av手机在线观看| 干丝袜人妻中文字幕| 欧美激情久久久久久爽电影| 日本午夜av视频| 欧美zozozo另类| 亚洲欧美一区二区三区国产| 听说在线观看完整版免费高清| 校园人妻丝袜中文字幕| 青春草国产在线视频| 一个人看视频在线观看www免费| 亚洲成人中文字幕在线播放| 啦啦啦中文免费视频观看日本| av在线老鸭窝| 深夜a级毛片| 最后的刺客免费高清国语| 亚洲在线观看片| 国产精品.久久久| 色哟哟·www| 国产一级毛片在线| 一个人看的www免费观看视频| 国产成人91sexporn| 在现免费观看毛片| 成年女人在线观看亚洲视频 | 精品久久久久久久久av| 麻豆国产97在线/欧美| 少妇人妻精品综合一区二区| 白带黄色成豆腐渣| 夫妻午夜视频| 国产一级毛片在线| 麻豆久久精品国产亚洲av| 久久久久久久午夜电影| 亚洲精品456在线播放app| 亚洲成色77777| xxx大片免费视频| 国产成人aa在线观看| 亚洲av中文字字幕乱码综合| 国产精品福利在线免费观看| 97在线视频观看| 国产老妇伦熟女老妇高清| 精品人妻一区二区三区麻豆| 亚洲三级黄色毛片| 大香蕉久久网| 街头女战士在线观看网站| 少妇的逼好多水| 亚洲va在线va天堂va国产| 麻豆av噜噜一区二区三区| 欧美极品一区二区三区四区| 国产在线男女| 日韩欧美 国产精品| 热99在线观看视频| 欧美性感艳星| 男人和女人高潮做爰伦理| 午夜激情福利司机影院| 美女高潮的动态| 午夜福利视频精品| 舔av片在线| av.在线天堂| 午夜福利在线观看免费完整高清在| 日产精品乱码卡一卡2卡三| 一级毛片 在线播放| 欧美精品国产亚洲| 一本久久精品| 国产高清三级在线| av专区在线播放| 国产精品麻豆人妻色哟哟久久 | 一级毛片 在线播放| 免费在线观看成人毛片| 最近中文字幕2019免费版| 床上黄色一级片| 男人狂女人下面高潮的视频| .国产精品久久| 有码 亚洲区| 99热这里只有是精品在线观看| 黄色配什么色好看| 菩萨蛮人人尽说江南好唐韦庄| 在现免费观看毛片| 精品少妇黑人巨大在线播放| 久久这里只有精品中国| 婷婷色综合www| 亚洲欧美成人精品一区二区| 亚洲怡红院男人天堂| av卡一久久| 大又大粗又爽又黄少妇毛片口| 日韩av在线免费看完整版不卡| 综合色丁香网| 99久久人妻综合| 亚洲经典国产精华液单| 亚洲综合色惰| 日韩欧美精品免费久久| 欧美一级a爱片免费观看看| 久久久久久九九精品二区国产| 超碰av人人做人人爽久久| 中文字幕制服av| 国产有黄有色有爽视频| 亚洲18禁久久av| av女优亚洲男人天堂| 可以在线观看毛片的网站| ponron亚洲| 久久精品久久精品一区二区三区| 观看免费一级毛片| 国产精品无大码| 最近手机中文字幕大全| 中文资源天堂在线| 国产精品av视频在线免费观看| 国产淫片久久久久久久久| 内射极品少妇av片p| 777米奇影视久久| 久久精品国产亚洲av天美| 神马国产精品三级电影在线观看| 大香蕉97超碰在线| 欧美+日韩+精品| 欧美精品国产亚洲| 色尼玛亚洲综合影院| av免费观看日本| 国产中年淑女户外野战色| 在线a可以看的网站| 乱人视频在线观看| 国产淫语在线视频| 亚洲精品乱久久久久久| 亚洲美女搞黄在线观看| 亚洲精品视频女| 久久国内精品自在自线图片| 欧美变态另类bdsm刘玥| 七月丁香在线播放| 婷婷色av中文字幕| 三级国产精品欧美在线观看| 六月丁香七月| 久久久久久久久久成人| 一级毛片 在线播放| 精品久久久久久成人av| 国产一区二区三区av在线| 久久精品夜色国产| av天堂中文字幕网| 99九九线精品视频在线观看视频| 男女那种视频在线观看| 汤姆久久久久久久影院中文字幕 | 国产成人a∨麻豆精品| 2018国产大陆天天弄谢| 建设人人有责人人尽责人人享有的 | 国产日韩欧美在线精品| 中文天堂在线官网| 毛片一级片免费看久久久久| 免费观看精品视频网站| 精品久久久久久电影网| 国产一级毛片在线| 久久久久久九九精品二区国产| 毛片一级片免费看久久久久| 久久久久久久久久久免费av| av又黄又爽大尺度在线免费看| a级毛色黄片| 日本三级黄在线观看| 久久久午夜欧美精品| 男人和女人高潮做爰伦理| 乱码一卡2卡4卡精品| 亚洲欧美日韩无卡精品| 美女cb高潮喷水在线观看| 欧美人与善性xxx| 搡老乐熟女国产| 亚洲丝袜综合中文字幕| 晚上一个人看的免费电影| av网站免费在线观看视频 | 欧美另类一区| 日日啪夜夜爽| 亚洲精品乱码久久久久久按摩| 精品久久久精品久久久| 日本熟妇午夜| 日日撸夜夜添| 又粗又硬又长又爽又黄的视频| 久久久久久久大尺度免费视频| av在线蜜桃| 最近中文字幕高清免费大全6| 亚洲三级黄色毛片| 又大又黄又爽视频免费| 人妻系列 视频| 国产一级毛片七仙女欲春2| 嫩草影院入口| 日韩成人伦理影院| 国产成人福利小说| 天堂俺去俺来也www色官网 | 超碰av人人做人人爽久久| 青春草视频在线免费观看| 少妇被粗大猛烈的视频| 熟女人妻精品中文字幕| 国产精品1区2区在线观看.| 国产高清国产精品国产三级 | 久久热精品热| 亚洲欧美日韩卡通动漫| 成人欧美大片| 亚洲怡红院男人天堂| 国产91av在线免费观看| 国产精品福利在线免费观看| 国产精品一二三区在线看| 久久热精品热| 久久久久网色| 久久久精品免费免费高清| 亚洲国产欧美在线一区| 国产综合懂色| 爱豆传媒免费全集在线观看| 日韩欧美国产在线观看| 九九久久精品国产亚洲av麻豆| 99久国产av精品| 国产精品福利在线免费观看| 非洲黑人性xxxx精品又粗又长| 国产成人精品久久久久久| 欧美日韩视频高清一区二区三区二| 国产又色又爽无遮挡免| 精品久久久久久久久av| 天堂网av新在线| 建设人人有责人人尽责人人享有的 | 在线 av 中文字幕| 日韩国内少妇激情av| 一级片'在线观看视频| 午夜久久久久精精品| 久久精品国产自在天天线| 日本黄色片子视频| 男女啪啪激烈高潮av片| 成人无遮挡网站| 高清视频免费观看一区二区 | 亚洲最大成人中文| 久久人人爽人人爽人人片va| 日韩欧美国产在线观看| 五月伊人婷婷丁香| 又爽又黄无遮挡网站| 男人舔奶头视频| 国产人妻一区二区三区在| 国产免费一级a男人的天堂| 欧美成人一区二区免费高清观看| 人妻制服诱惑在线中文字幕| 一区二区三区四区激情视频| 欧美激情国产日韩精品一区| 久久久久久久午夜电影| 亚洲av在线观看美女高潮| 只有这里有精品99| 亚洲欧美清纯卡通| 毛片一级片免费看久久久久| 亚洲精品456在线播放app| 久久久久久久久久人人人人人人| 少妇的逼好多水| 日韩欧美国产在线观看| 国产伦精品一区二区三区视频9| 内地一区二区视频在线| 日韩中字成人| 99热全是精品| 亚洲真实伦在线观看| 日韩一区二区三区影片| 国产av在哪里看| 人妻系列 视频| 亚洲精品,欧美精品| 亚洲成人一二三区av| 只有这里有精品99| 韩国高清视频一区二区三区| 亚洲真实伦在线观看| 国产乱来视频区| 国产爱豆传媒在线观看| 久久精品夜色国产| 91午夜精品亚洲一区二区三区| 简卡轻食公司| 国产精品人妻久久久久久| 欧美+日韩+精品| 精品久久国产蜜桃| 欧美变态另类bdsm刘玥| 国产精品精品国产色婷婷| 欧美高清成人免费视频www| 男女国产视频网站| 成人午夜精彩视频在线观看| 久久久精品94久久精品| 国产精品人妻久久久久久| 狂野欧美激情性xxxx在线观看| 亚洲欧美成人精品一区二区| 亚洲欧美日韩卡通动漫| 国产男女超爽视频在线观看| 男人和女人高潮做爰伦理| 国产精品综合久久久久久久免费| 国产亚洲精品久久久com| 精品一区二区三区人妻视频| 乱系列少妇在线播放| 亚洲精品自拍成人| 麻豆av噜噜一区二区三区| 亚洲成人精品中文字幕电影| 男人狂女人下面高潮的视频| 最新中文字幕久久久久| 精品99又大又爽又粗少妇毛片| 自拍偷自拍亚洲精品老妇| 成人鲁丝片一二三区免费| 久久国产乱子免费精品| 女的被弄到高潮叫床怎么办| 国产精品无大码| 如何舔出高潮| 精品国产三级普通话版| 丝瓜视频免费看黄片| 菩萨蛮人人尽说江南好唐韦庄| 国产综合懂色| 国产亚洲91精品色在线| 夜夜爽夜夜爽视频| 国产一区二区亚洲精品在线观看| 青春草视频在线免费观看| 久久久久久久午夜电影| 精品久久久久久久人妻蜜臀av| 狂野欧美激情性xxxx在线观看| 天堂av国产一区二区熟女人妻| 国产精品日韩av在线免费观看| 又爽又黄a免费视频| 免费av毛片视频| 国产精品麻豆人妻色哟哟久久 | 蜜桃亚洲精品一区二区三区| 搡女人真爽免费视频火全软件| 亚洲图色成人| 爱豆传媒免费全集在线观看| 精品久久久久久久人妻蜜臀av| 亚洲av电影在线观看一区二区三区 | 国产成人aa在线观看| 国产精品国产三级专区第一集| 一级av片app| 免费播放大片免费观看视频在线观看| 99久久人妻综合| 一级毛片 在线播放| 免费观看精品视频网站| 国产精品女同一区二区软件| 国产精品一区二区三区四区免费观看| 美女内射精品一级片tv| 国产av码专区亚洲av| 简卡轻食公司| 美女高潮的动态| 久久精品国产亚洲av涩爱| 国产黄片视频在线免费观看| 成人亚洲精品一区在线观看 | 亚洲久久久久久中文字幕| 日韩大片免费观看网站| 日韩视频在线欧美| 久久精品综合一区二区三区| 免费看不卡的av| 在线观看美女被高潮喷水网站| 少妇熟女aⅴ在线视频| 女人十人毛片免费观看3o分钟| 少妇熟女aⅴ在线视频| 18禁裸乳无遮挡免费网站照片| 成人毛片a级毛片在线播放| 国产成人aa在线观看| 亚洲精品久久久久久婷婷小说| 国产精品蜜桃在线观看| 久久这里只有精品中国| 中文字幕久久专区| 99九九线精品视频在线观看视频| 午夜免费男女啪啪视频观看| 熟妇人妻久久中文字幕3abv|