• <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
    朝鮮孢囊線蟲——浙江省孢囊線蟲新記錄種
    直男gayav资源| 综合色丁香网| 伦理电影免费视频| 91精品国产九色| 亚洲av成人精品一区久久| 肉色欧美久久久久久久蜜桃| 人妻制服诱惑在线中文字幕| 亚洲丝袜综合中文字幕| 国产一级毛片在线| 亚洲av在线观看美女高潮| 国产成人午夜福利电影在线观看| 国产精品一区二区在线观看99| 国内少妇人妻偷人精品xxx网站| 爱豆传媒免费全集在线观看| 激情 狠狠 欧美| 色视频在线一区二区三区| 91精品伊人久久大香线蕉| 国产片特级美女逼逼视频| 99久久综合免费| 黄色视频在线播放观看不卡| 亚洲精品国产av蜜桃| 成人一区二区视频在线观看| 国产免费福利视频在线观看| 成人毛片60女人毛片免费| 婷婷色综合大香蕉| 久久久久精品久久久久真实原创| 美女主播在线视频| 日本vs欧美在线观看视频 | 国产免费一级a男人的天堂| 国产成人免费无遮挡视频| 在线观看一区二区三区激情| a级毛色黄片| a 毛片基地| 最近中文字幕高清免费大全6| 日本免费在线观看一区| 久久久久久久精品精品| 日本黄色日本黄色录像| 热99国产精品久久久久久7| 一级毛片黄色毛片免费观看视频| 秋霞在线观看毛片| 在线观看一区二区三区| 三级经典国产精品| 国产亚洲一区二区精品| 激情五月婷婷亚洲| 男女免费视频国产| 国产一区有黄有色的免费视频| 啦啦啦中文免费视频观看日本| 国产精品秋霞免费鲁丝片| 少妇高潮的动态图| 精品一区二区免费观看| 久久久久久久久久成人| av在线app专区| 91久久精品国产一区二区三区| 插逼视频在线观看| 九九爱精品视频在线观看| 啦啦啦中文免费视频观看日本| 熟女电影av网| 99久久人妻综合| 97超视频在线观看视频| 国产高清不卡午夜福利| 亚洲人与动物交配视频| 99久久精品热视频| 中文天堂在线官网| av专区在线播放| 欧美丝袜亚洲另类| 欧美另类一区| 欧美日韩视频精品一区| 亚洲无线观看免费| 妹子高潮喷水视频| 午夜老司机福利剧场| 欧美三级亚洲精品| 免费黄网站久久成人精品| 啦啦啦中文免费视频观看日本| 国产午夜精品久久久久久一区二区三区| 亚洲无线观看免费| 男人和女人高潮做爰伦理| 日日撸夜夜添| 女人久久www免费人成看片| videos熟女内射| 国产亚洲精品久久久com| 国产男人的电影天堂91| 午夜福利网站1000一区二区三区| 久久久久久久亚洲中文字幕| videossex国产| 深爱激情五月婷婷| 啦啦啦在线观看免费高清www| 插逼视频在线观看| 免费看不卡的av| 亚洲成人一二三区av| 亚州av有码| 久久人妻熟女aⅴ| 在线精品无人区一区二区三 | av女优亚洲男人天堂| 男女边吃奶边做爰视频| 国产女主播在线喷水免费视频网站| 久久久久国产网址| 国产精品一二三区在线看| 国产精品99久久99久久久不卡 | 国产淫语在线视频| 99热这里只有精品一区| 亚洲av免费高清在线观看| 青春草亚洲视频在线观看| 日本wwww免费看| 欧美亚洲 丝袜 人妻 在线| 日韩欧美 国产精品| 亚洲精华国产精华液的使用体验| a级一级毛片免费在线观看| 亚洲人成网站在线观看播放| 亚洲,欧美,日韩| 波野结衣二区三区在线| 两个人的视频大全免费| 国产精品国产三级专区第一集| 亚洲美女黄色视频免费看| 视频区图区小说| 色吧在线观看| 久久人人爽人人爽人人片va| 日本猛色少妇xxxxx猛交久久| 一级毛片电影观看| 久久久久久久精品精品| 一区二区三区乱码不卡18| 亚洲精品第二区| 在现免费观看毛片| 欧美xxxx性猛交bbbb| 日本黄色片子视频| 涩涩av久久男人的天堂| 成人综合一区亚洲| 欧美+日韩+精品| 天堂中文最新版在线下载| 精品久久久久久久久亚洲| 日韩制服骚丝袜av| 一级a做视频免费观看| 97超视频在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩欧美 国产精品| 身体一侧抽搐| 这个男人来自地球电影免费观看 | 久久鲁丝午夜福利片| 亚洲成人av在线免费| 国产极品天堂在线| 日韩制服骚丝袜av| 一级毛片电影观看| 美女xxoo啪啪120秒动态图| 精品久久久精品久久久| 在线免费观看不下载黄p国产| 国产精品蜜桃在线观看| 最近最新中文字幕大全电影3| 成人午夜精彩视频在线观看| 老司机影院毛片| 久久热精品热| 亚洲欧美成人综合另类久久久| 欧美少妇被猛烈插入视频| 日韩欧美精品免费久久| 亚洲人与动物交配视频| 青春草国产在线视频| 日本欧美国产在线视频| 欧美丝袜亚洲另类| 观看免费一级毛片| 亚洲性久久影院| 人妻夜夜爽99麻豆av| 精品人妻视频免费看| 欧美精品国产亚洲| 精品少妇久久久久久888优播| 2021少妇久久久久久久久久久| 午夜激情福利司机影院| 中文在线观看免费www的网站| 在线精品无人区一区二区三 | 欧美另类一区| 黄色怎么调成土黄色| 最近的中文字幕免费完整| 国产大屁股一区二区在线视频| 亚洲无线观看免费| 久久国产精品大桥未久av | 三级国产精品欧美在线观看| 精品人妻视频免费看| 亚洲内射少妇av| 国精品久久久久久国模美| 在线免费十八禁| 国产大屁股一区二区在线视频| 亚洲无线观看免费| 午夜福利在线观看免费完整高清在| 街头女战士在线观看网站| 中文精品一卡2卡3卡4更新| 国产精品久久久久久av不卡| 51国产日韩欧美| 高清毛片免费看| 国产又色又爽无遮挡免| 欧美xxⅹ黑人| 国产 一区精品| 97在线视频观看| 一本一本综合久久| 中文字幕免费在线视频6| 18+在线观看网站| 亚洲欧美精品自产自拍| 搡女人真爽免费视频火全软件| 成人18禁高潮啪啪吃奶动态图 | 久久久久久久国产电影| 欧美老熟妇乱子伦牲交| 寂寞人妻少妇视频99o| 国产精品秋霞免费鲁丝片| 亚洲国产精品一区三区| 一本一本综合久久| 女的被弄到高潮叫床怎么办| 日韩成人伦理影院| 国产精品一区二区性色av| 少妇人妻 视频| 国产精品精品国产色婷婷| 亚洲国产精品成人久久小说| 国产乱人偷精品视频| 国产av一区二区精品久久 | 国产精品秋霞免费鲁丝片| 成人影院久久| 午夜老司机福利剧场| 国产精品国产三级国产专区5o| 男女无遮挡免费网站观看| 日韩,欧美,国产一区二区三区| 我的老师免费观看完整版| 免费观看a级毛片全部| 成人特级av手机在线观看| 高清不卡的av网站| 国产精品一二三区在线看| h视频一区二区三区| 人体艺术视频欧美日本| 国产黄色免费在线视频| 国产欧美另类精品又又久久亚洲欧美| 91久久精品国产一区二区成人| a级毛色黄片| 赤兔流量卡办理| 最近中文字幕2019免费版| 国内揄拍国产精品人妻在线| 亚洲色图综合在线观看| 18禁动态无遮挡网站| 国产在线视频一区二区| 日本av免费视频播放| 久久影院123| 中文字幕精品免费在线观看视频 | 亚洲美女搞黄在线观看| 国产精品久久久久久精品古装| 一级爰片在线观看| 波野结衣二区三区在线| 久久久久性生活片| 青春草视频在线免费观看| 亚洲成人av在线免费| 搡女人真爽免费视频火全软件| 卡戴珊不雅视频在线播放| 免费看av在线观看网站| 亚洲人成网站在线观看播放| 成人毛片a级毛片在线播放| 国产成人免费观看mmmm| 精品久久久久久久末码| 老熟女久久久| 亚洲成色77777| 国产成人91sexporn| 王馨瑶露胸无遮挡在线观看| 九九爱精品视频在线观看| 少妇高潮的动态图| 哪个播放器可以免费观看大片| 日韩成人av中文字幕在线观看| 久热久热在线精品观看| 日本黄色日本黄色录像| 日本色播在线视频| 日本av手机在线免费观看| 国产无遮挡羞羞视频在线观看| 日韩在线高清观看一区二区三区| 午夜日本视频在线| 一级毛片电影观看| 一级黄片播放器| 欧美日本视频| 国产精品福利在线免费观看| 国产精品精品国产色婷婷| 少妇熟女欧美另类| 嘟嘟电影网在线观看| 赤兔流量卡办理| 男的添女的下面高潮视频| 亚洲精品乱码久久久久久按摩| 91狼人影院| 欧美精品一区二区免费开放| 中文字幕精品免费在线观看视频 | 我的老师免费观看完整版| 欧美少妇被猛烈插入视频| 国产一级毛片在线| 亚洲成色77777| 国产精品三级大全| 亚洲精品456在线播放app| 亚洲欧美精品自产自拍| 天堂8中文在线网| 亚洲国产精品国产精品| 亚洲av中文av极速乱| 卡戴珊不雅视频在线播放| 日韩欧美一区视频在线观看 | 一边亲一边摸免费视频| 狠狠精品人妻久久久久久综合| 国产高潮美女av| 男男h啪啪无遮挡| 国产大屁股一区二区在线视频| 久久99蜜桃精品久久| 国产亚洲精品久久久com| 伦精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 老司机影院成人| 成年女人在线观看亚洲视频| 精品国产三级普通话版| 少妇人妻精品综合一区二区| 一个人看的www免费观看视频| 麻豆国产97在线/欧美| 成人二区视频| 一级毛片黄色毛片免费观看视频| 97在线视频观看| 如何舔出高潮| 色5月婷婷丁香| 成人二区视频| 在线观看免费视频网站a站| 国产乱人视频| 日韩强制内射视频| 成年免费大片在线观看| 国产精品人妻久久久久久| 婷婷色麻豆天堂久久| 国产精品精品国产色婷婷| 香蕉精品网在线| 精品国产一区二区三区久久久樱花 | 嘟嘟电影网在线观看| 女人久久www免费人成看片| 亚洲一级一片aⅴ在线观看| 国产成人aa在线观看| 人妻一区二区av| 亚洲欧美成人综合另类久久久| 中文字幕精品免费在线观看视频 | 尾随美女入室| 亚洲美女视频黄频| 午夜福利影视在线免费观看| 欧美97在线视频| 欧美亚洲 丝袜 人妻 在线| 少妇猛男粗大的猛烈进出视频| 亚洲国产精品国产精品| 午夜福利在线观看免费完整高清在| 日韩av免费高清视频| 久久99热6这里只有精品| kizo精华| 免费观看a级毛片全部| 性色avwww在线观看| 欧美一级a爱片免费观看看| 亚洲一区二区三区欧美精品| 久久人人爽av亚洲精品天堂 | 在线观看一区二区三区| 中文字幕精品免费在线观看视频 | 美女视频免费永久观看网站| 国产黄片美女视频| 一区二区三区精品91| 久久综合国产亚洲精品| 美女高潮的动态| 美女中出高潮动态图| 国产一区二区三区av在线| 亚洲va在线va天堂va国产| 免费人成在线观看视频色| 人体艺术视频欧美日本| 久久久久久久久大av| 精品人妻熟女av久视频| 亚洲久久久国产精品| 少妇人妻 视频| 水蜜桃什么品种好| 女性生殖器流出的白浆| 国产亚洲91精品色在线| 亚洲国产精品成人久久小说| 欧美日本视频| 久久99蜜桃精品久久| 大片免费播放器 马上看| 成人综合一区亚洲| av播播在线观看一区| 精品久久久精品久久久| av视频免费观看在线观看| 欧美日韩在线观看h| 插阴视频在线观看视频| 蜜桃在线观看..| 草草在线视频免费看| 亚洲欧美中文字幕日韩二区| 国产成人freesex在线| 九色成人免费人妻av| 嘟嘟电影网在线观看| 欧美xxⅹ黑人| 欧美一区二区亚洲| tube8黄色片| 91aial.com中文字幕在线观看| 丰满少妇做爰视频| 中国美白少妇内射xxxbb| 免费不卡的大黄色大毛片视频在线观看| 免费观看无遮挡的男女| 一级毛片电影观看| 国产av精品麻豆| 最近的中文字幕免费完整| 国产伦精品一区二区三区视频9| 在线亚洲精品国产二区图片欧美 | 丰满人妻一区二区三区视频av| 亚洲av免费高清在线观看| 国产伦理片在线播放av一区| av免费在线看不卡| 在线免费观看不下载黄p国产| 国产黄片美女视频| 欧美另类一区| 精品99又大又爽又粗少妇毛片| 交换朋友夫妻互换小说| 亚洲aⅴ乱码一区二区在线播放| 国产成人精品一,二区| 精华霜和精华液先用哪个| 在线 av 中文字幕| 久久鲁丝午夜福利片| 国产亚洲91精品色在线| 欧美极品一区二区三区四区| 纵有疾风起免费观看全集完整版| av国产精品久久久久影院| 青春草视频在线免费观看| 亚洲天堂av无毛| 亚洲欧美中文字幕日韩二区| 国产一区二区三区综合在线观看 | 欧美zozozo另类| 久热这里只有精品99| 国产伦精品一区二区三区视频9| 亚洲精品第二区| 久久久久久久久久人人人人人人| 成人特级av手机在线观看| 中文在线观看免费www的网站| 免费高清在线观看视频在线观看| 好男人视频免费观看在线| 亚洲丝袜综合中文字幕| 婷婷色综合大香蕉| 一本久久精品| 国产探花极品一区二区| 亚洲欧洲日产国产| 尤物成人国产欧美一区二区三区| 国产精品成人在线| 国产精品久久久久成人av| 国产成人精品福利久久| 久久人人爽人人爽人人片va| 天堂俺去俺来也www色官网| 全区人妻精品视频| 亚洲精品久久午夜乱码| 久热这里只有精品99| 大码成人一级视频| 国产欧美日韩一区二区三区在线 | 香蕉精品网在线| 久久久亚洲精品成人影院| 一本一本综合久久| 国产精品.久久久| 欧美国产精品一级二级三级 | 精品亚洲乱码少妇综合久久| 在线观看免费日韩欧美大片 | 日韩成人伦理影院| 亚洲av中文av极速乱| 美女xxoo啪啪120秒动态图| 嘟嘟电影网在线观看| 国产精品偷伦视频观看了| 久久久久网色| 99精国产麻豆久久婷婷| 一区二区三区四区激情视频| 国产精品国产av在线观看| 亚州av有码| 欧美极品一区二区三区四区| 三级国产精品欧美在线观看| 亚洲aⅴ乱码一区二区在线播放| 中文资源天堂在线| 国产亚洲91精品色在线| 亚洲不卡免费看| 亚洲成人中文字幕在线播放| 久久精品夜色国产| 高清午夜精品一区二区三区| 一二三四中文在线观看免费高清| 18+在线观看网站| 国产亚洲91精品色在线| 内地一区二区视频在线| 国产精品偷伦视频观看了| 日韩 亚洲 欧美在线| 亚洲欧美一区二区三区黑人 | 亚洲欧美日韩无卡精品| 亚洲欧美一区二区三区国产| 国产精品久久久久久久电影| 一区二区三区四区激情视频| 国产伦精品一区二区三区四那| 免费播放大片免费观看视频在线观看| 熟女人妻精品中文字幕| 天堂中文最新版在线下载| 大片免费播放器 马上看| 精品酒店卫生间| 亚洲色图综合在线观看| 天天躁日日操中文字幕| 大又大粗又爽又黄少妇毛片口| 亚洲国产av新网站| 久久久久久久国产电影| 亚洲欧美日韩无卡精品| 中文乱码字字幕精品一区二区三区| 成人亚洲精品一区在线观看 | 中文字幕制服av| av免费观看日本| 色视频www国产| 亚洲精品日韩av片在线观看| 国产精品一二三区在线看| freevideosex欧美| 99久久精品国产国产毛片| 色网站视频免费| 能在线免费看毛片的网站| 黄色一级大片看看| 午夜免费鲁丝| 18+在线观看网站| 国产精品三级大全| 亚洲国产毛片av蜜桃av| 草草在线视频免费看| 一级av片app| 亚洲欧美日韩无卡精品| 啦啦啦啦在线视频资源| 精品人妻偷拍中文字幕| 在现免费观看毛片| 国产成人精品久久久久久| 一区二区三区乱码不卡18| 夫妻性生交免费视频一级片| 欧美zozozo另类| 久久久久久久久久成人| 狂野欧美激情性xxxx在线观看| 联通29元200g的流量卡| 午夜激情久久久久久久| 精品久久久噜噜| 久久久成人免费电影| 久久精品国产a三级三级三级| 少妇被粗大猛烈的视频| 久热久热在线精品观看| 一区在线观看完整版| 人妻系列 视频| 日韩av不卡免费在线播放| 老司机影院成人| 一级毛片黄色毛片免费观看视频| 97热精品久久久久久| 色综合色国产| 一二三四中文在线观看免费高清| 黄色一级大片看看| 女人十人毛片免费观看3o分钟| 久久鲁丝午夜福利片| 国产高清有码在线观看视频| 国产在线视频一区二区| 免费观看性生交大片5| 人人妻人人爽人人添夜夜欢视频 | 精品国产乱码久久久久久小说| 在线亚洲精品国产二区图片欧美 | 亚洲图色成人| 精品国产乱码久久久久久小说| 韩国高清视频一区二区三区| 日产精品乱码卡一卡2卡三| 国产精品一二三区在线看| 国产 一区精品| 婷婷色av中文字幕| 亚洲成人中文字幕在线播放| 国产精品蜜桃在线观看| 五月伊人婷婷丁香| 熟妇人妻不卡中文字幕| 国产精品av视频在线免费观看| 九九爱精品视频在线观看| 日本vs欧美在线观看视频 | 国产色爽女视频免费观看| 一本—道久久a久久精品蜜桃钙片| 涩涩av久久男人的天堂| 国产精品一及| 欧美一级a爱片免费观看看| 久久人人爽人人爽人人片va| 午夜视频国产福利| 久久影院123| 麻豆成人av视频| 国产av国产精品国产| 网址你懂的国产日韩在线| 王馨瑶露胸无遮挡在线观看| 插逼视频在线观看| 精品久久久精品久久久| 九九在线视频观看精品| 色婷婷久久久亚洲欧美| 国产爱豆传媒在线观看| 一级黄片播放器| 日本av免费视频播放| 蜜桃久久精品国产亚洲av| 色婷婷av一区二区三区视频| 国产久久久一区二区三区| 99久久精品热视频| 国产成人a区在线观看| 国国产精品蜜臀av免费| 国产精品不卡视频一区二区| 嫩草影院新地址| 久久久久久久久久久免费av| 亚洲综合精品二区| 交换朋友夫妻互换小说| 国产精品国产三级国产专区5o| 哪个播放器可以免费观看大片| 狠狠精品人妻久久久久久综合| 国产成人91sexporn| 国产色婷婷99| 国产视频首页在线观看| 人妻一区二区av| 亚洲aⅴ乱码一区二区在线播放| 国产av码专区亚洲av| 国产精品久久久久久精品电影小说 | 26uuu在线亚洲综合色| 极品少妇高潮喷水抽搐| 五月玫瑰六月丁香| 另类亚洲欧美激情| 丰满迷人的少妇在线观看| www.av在线官网国产| 最后的刺客免费高清国语| 中文天堂在线官网| 美女视频免费永久观看网站| 久久精品国产自在天天线| 新久久久久国产一级毛片| 久久久久久久久大av| 国产色婷婷99| 久久ye,这里只有精品| 国产成人freesex在线| 国产成人a∨麻豆精品| av.在线天堂| 亚洲国产精品成人久久小说| 亚洲av国产av综合av卡| 免费黄色在线免费观看| 亚洲成人中文字幕在线播放| 极品少妇高潮喷水抽搐| 2018国产大陆天天弄谢|