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

    南沙群島海域長腹劍水蚤(Oithona spp.)的種類組成、數(shù)量分布及其與環(huán)境因子的關(guān)系*

    2016-01-15 11:22:37杜飛雁王亮根王雪輝寧加佳谷陽光李亞芳
    海洋與湖沼 2016年6期
    關(guān)鍵詞:組群水蚤南沙群島

    杜飛雁 王亮根 王雪輝 寧加佳 谷陽光 李亞芳

    (廣東省漁業(yè)生態(tài)環(huán)境重點(diǎn)開放實(shí)驗(yàn)室 農(nóng)業(yè)部南海漁業(yè)資源環(huán)境科學(xué)觀測(cè)實(shí)驗(yàn)站 中國水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所廣州 510300)

    橈足類(Copepods)是海洋中小型浮游動(dòng)物(0.2—20mm)中種類最豐富、數(shù)量最大的類群, 在海洋生態(tài)系統(tǒng)中具有重要的地位和作用。長腹劍水蚤屬(Oithona)屬于甲殼動(dòng)物亞門(Crustacea)顎足綱(Maxillopoda)橈足亞綱(Copepoda)劍水蚤目(Cyclopoidea)長腹劍水蚤科(Oithonidae), 是海洋中小型橈足類數(shù)量最大的類群之一, 廣泛分布于近岸和外海中(Paffenh?fer, 1993), 在海洋食物網(wǎng)、微食物環(huán)與碳循環(huán)中處于關(guān)鍵的位置(Nakamura et al, 1997;Saiz et al, 2003; Tumer, 2004)。長腹劍水蚤是海洋食物網(wǎng)中重要的中間捕食者, 以單細(xì)胞生物和橈足類幼蟲為食, 既為仔魚和其他食浮游動(dòng)物的動(dòng)物提供餌料, 又是微食物環(huán)中的食碎屑者, 在碳循環(huán)中細(xì)菌碳和藻類碳向高營養(yǎng)級(jí)轉(zhuǎn)化的過程里扮演著重要角色(Hwang et al, 2010)。我國海域目前已記錄的長腹劍水蚤有17種, 在南海均有分布(劉瑞玉, 2008)。長腹劍水蚤個(gè)體小且形態(tài)差異微小, 采用傳統(tǒng)形態(tài)學(xué)分類法對(duì)其準(zhǔn)確鑒定具有一定的困難, 存在大量的隱種, 且該屬種間具體的系統(tǒng)進(jìn)化關(guān)系尚不明確, 影響了海洋食物網(wǎng)結(jié)構(gòu)與功能和海洋生態(tài)系統(tǒng)相關(guān)研究的深入開展。

    南沙群島海域地處我國最南端的熱帶海域, 地理、地貌、氣候和生態(tài)環(huán)境獨(dú)特, 水文環(huán)境復(fù)雜(方文東, 2000), 島嶼眾多, 緊鄰大陸坡和深海盆地, 季風(fēng)盛行, 環(huán)流形勢(shì)復(fù)雜多變(楊海軍等, 1998), 是海洋生物多樣性最為豐富的區(qū)域之一。南沙群島海域目前已有記錄的長腹劍水蚤共有 12種, 但其數(shù)量分布和季節(jié)變化情況尚不清楚(陳清潮, 2003)。本文以2013年在南沙群島海域按季度月開展的 4個(gè)航次的生態(tài)調(diào)查數(shù)據(jù)為基礎(chǔ), 對(duì)該海域長腹劍水蚤的種類組成、組群結(jié)構(gòu)、數(shù)量的水平和垂直分布及其季節(jié)變化進(jìn)行研究, 為后續(xù)深入開展長腹劍水蚤的分子生物學(xué)研究提供基礎(chǔ)資料。

    此外, 浮游生物與影響因子之間關(guān)系復(fù)雜, 不是簡單的線性響應(yīng), 故引入非線性的廣義可加模型(generalized additive model, GAM)分析長腹劍水蚤與環(huán)境因子之間的關(guān)系。對(duì)于單物種與環(huán)境之間定量關(guān)系的研究, GAM 模型主要取決于原始數(shù)據(jù), 能更加深入地探討物種與環(huán)境變量間的關(guān)系, 提供更為豐富的信息, 是目前被認(rèn)為較先進(jìn)的方法(Hastie et al,1990; Bigelow et al, 1999;牛明香等, 2010)。GAM 模型雖然在生態(tài)學(xué)領(lǐng)域運(yùn)用廣泛(Christos et al, 1997;劉佳等, 2009; Chang et al, 2010; 唐浩等, 2013; Bai et al, 2015;), 但在浮游動(dòng)物生態(tài)學(xué)研究方面的應(yīng)用相對(duì)較少。

    1 材料與方法

    1.1 調(diào)查設(shè)置及取樣方法

    2013年按季度月分別于3—4月、6—7月、9—10月和 11—12月, 對(duì)南沙群島海域(111°30′E—112°30′E, 4°30′N—11°30′N)設(shè)置的 35 個(gè)站位, 進(jìn)行春、夏、秋、冬季4個(gè)航次生態(tài)環(huán)境的綜合調(diào)查。采樣過程中因惡劣海況和外界干擾, 各航次的作站情況有所差異, 其中春季完成23個(gè)站位(不含S1、S3、S5、S6、S8、S11、S17、S25、S26、S28、S29 和 S35站)、夏季完成21個(gè)站位(不含S7、S9、S10、S11、S12、S16、S18、S19、S23、S26、S27、S29、S33和S34站)、秋季完成24個(gè)站位(不含S1、S3、S8、S11、S17、S18、S25、S26、S29、S34 和 S35 站)、冬季完成18個(gè)站位(不含S1、S3、S6、S7、S8、S9、S13、S17、S19、S22、S23、S25、S7、S28、S33、S34和S35站), 站位地理位置見圖1和表1所示。每站分別使用中型浮游生物網(wǎng)(網(wǎng)長 280cm, 網(wǎng)口內(nèi)徑50cm, 網(wǎng)口面積 0.2m2, 網(wǎng)目大小 0.077mm)自 200m至表層垂直拖曳, 共采獲樣品86份。此外, 為了解南沙群島海域浮游動(dòng)物的垂直分布, 于2013年11月分別在 9°30′N 和 111°30′E 兩個(gè)斷面設(shè)置 3 個(gè)站位, 共 5個(gè)站位(S4、S14、S15、S16和S33站), 使用Multinet浮游生物連續(xù)采樣網(wǎng)(HYDRO-BIOS德國, 網(wǎng)口面積0.25m2, 網(wǎng)目大小 300μm), 分 0—30m、30—75m、75—150m、150—500m和500—750m等5層, 進(jìn)行浮游動(dòng)物垂直分層采集, 每站采集1次, 共采獲樣品22份。樣品的處理、保存、計(jì)數(shù)等均按《海洋調(diào)查規(guī)范》(GBT12763.6-2007)進(jìn)行。

    海水溫度(T)、鹽度(S)使用CTD(Sea-Bird美國)測(cè)定。采集5、10、25、50、75、100、150和 200m水樣, 使用10-AU-005-CE熒光儀(Turner Designs美國)測(cè)定葉綠素a(chl a)。

    圖1 調(diào)查站位設(shè)置示意圖Fig.1 Sketch map of sampling stations

    1.2 數(shù)據(jù)統(tǒng)計(jì)及分析

    1.2.1 出現(xiàn)頻率 通過出現(xiàn)頻率反映長腹劍水蚤單一種在調(diào)查期的出現(xiàn)情況, 計(jì)算式為:

    其中:a為該種出現(xiàn)站次,b為總采樣站次。

    1.2.2 種類更替 選用種類更替率(R)進(jìn)行種類組成的季節(jié)變化研究(連光山等, 1990), 其計(jì)算式為:

    其中:a與b分別為相鄰2個(gè)季節(jié)的種數(shù),c為相鄰2個(gè)季節(jié)共同的種數(shù)。

    1.2.3 組群分析 通過聚類分析研究長腹劍水蚤的組群結(jié)構(gòu)。首先將全年 0—200m垂直采集的在各采樣站次中出現(xiàn)次數(shù)大于或等于3的11種長腹劍水蚤, 按其出現(xiàn)與或缺情況進(jìn)行聚類分析和組群劃分。聚類分析采用距離系數(shù)(D=1-相關(guān)系數(shù)), 建立種類的二維矩陣, 進(jìn)而采用組間連接法(between-groups linkage)進(jìn)行聚類分析。數(shù)據(jù)處理采用PRIMER6.0統(tǒng)計(jì)軟件進(jìn)行。

    1.2.4 GAM 模型分析 GAM 模型的優(yōu)點(diǎn)是通過非參數(shù)平滑曲線建立預(yù)測(cè)變量和響應(yīng)變量之間的非線性關(guān)系(Maravelias, 2001; Wood, 2006), 且預(yù)測(cè)變量之間的關(guān)系是可加的。它能應(yīng)用非參數(shù)的方法檢測(cè)數(shù)據(jù)結(jié)構(gòu), 并找出其中的規(guī)律, 從而能得到更好的預(yù)測(cè)結(jié)果(Wood, 2006)。通過GAM模型分析長腹劍水蚤棲息密度與環(huán)境因子之間的關(guān)系, 表達(dá)式如下:

    其中,D為長腹劍水蚤調(diào)查期各站棲息密度,M為采樣月份, T為海水平均溫度,S為海水鹽度, chla為葉綠素a濃度(溫度、鹽度和葉綠素a濃度, 取長腹劍水蚤主要分布水層的平均值), lat為采樣緯度, lon為采樣經(jīng)度,s是光滑樣條函數(shù)(spline smoother),ε為隨機(jī)誤差項(xiàng), 利用t檢驗(yàn)評(píng)估預(yù)測(cè)變量的顯著性。以上模型的構(gòu)建過程均利用 R統(tǒng)計(jì)軟件(Version 3.3.0,https:// www.r-project.org/)中的mgcv軟件包實(shí)現(xiàn)。

    2 結(jié)果

    2.1 種類組成

    該海域在調(diào)查期間共出現(xiàn)長腹劍水蚤15種(包括1個(gè)未定種), 其中線長腹水蚤(Oithona linearis)、粗長腹劍水蚤(O. robusta)和長刺長腹劍水蚤(O.longispina)為南沙群島海域首次記錄。全部86個(gè)調(diào)查站位中, 除1個(gè)站沒有采獲外, 其余85個(gè)站均有長腹劍水蚤分布。

    各季中, 秋季出現(xiàn)種類數(shù)最多, 有12種; 春季次之, 出現(xiàn)11種; 冬季出現(xiàn)9種; 夏季最少, 僅有6種(表 1)。南沙群島四季均有出現(xiàn)的種類僅有羽長腹劍水蚤(O. plumifera)、瘦長腹劍水蚤(O. tenuis)、擬長腹劍水蚤(O. similis)、細(xì)長腹劍水蚤(O. attenuata)和堅(jiān)長腹劍水蚤(O. rigida)等 5種, 其年均出現(xiàn)頻率分別為98.9%、84.5%、90.6%、26.1%和26.8%, 其他種類均為季節(jié)性出現(xiàn)種。簡長腹劍水蚤(O. simplex)和長刺長腹劍水蚤除夏季外, 各季均有出現(xiàn), 年均出現(xiàn)頻率為7.7%和8.1%。偽長腹劍水蚤(O. fallax)、刺長腹劍水蚤(O. setigera)、敏長腹劍水蚤(O. vivida)和未定種等4種為2個(gè)季節(jié)出現(xiàn)的種類, 年均出現(xiàn)頻率在2.5%—18.9%。1個(gè)季節(jié)出現(xiàn)種包括線長腹水蚤(春季)、隱長腹劍水蚤(O. decipiens)(春季)、粗長腹劍水蚤(秋季)和短角長腹劍水蚤(O. brevicornis)(秋季)等4種, 其出現(xiàn)頻率也很低, 為3.3%和1.1%。全部15種長腹劍水蚤平均年出現(xiàn)頻率為25.8%。春夏季、夏秋季和秋冬季長腹劍水蚤的季節(jié)更替率分別為 45.5%、61.5%和25%, 平均更替率為44%。表明長腹劍水蚤屬雖然在南沙群島海域廣泛分布, 但單種的出現(xiàn)頻率較低, 其種類出現(xiàn)的季節(jié)性較強(qiáng)。

    表1 南沙群島海域長腹劍水蚤的出現(xiàn)頻率(%)Tab.1 The frequency of occurrence of Oithona spp. in Nansha Islands of South China Sea

    2.2 組群結(jié)構(gòu)

    群落中具有相似的分布特性, 在采樣中傾向于同時(shí)出現(xiàn)的一些物種可稱為物種組群, 鑒別物種組群并探討其形成的原因, 有助于判明群落結(jié)構(gòu)和物種生態(tài)特性, 以及物種之間生態(tài)上的相互聯(lián)系(邱永松, 1990), 對(duì)深入了解長腹劍水蚤的生態(tài)習(xí)性有重要作用。由聚類圖可以得出, 南沙群島海域的長腹劍水蚤可以劃分為2個(gè)組群(圖2)。組群Ⅰ組成較為簡單,由瘦長腹劍水蚤、擬長腹劍水蚤和羽長腹劍水蚤3種構(gòu)成, 但種間相似性較高(距離為 2.24—3.67); 組群Ⅱ由堅(jiān)長腹劍水蚤、細(xì)長腹劍水蚤、刺長腹劍水蚤、偽長腹劍水蚤、簡長腹劍水蚤、敏長腹劍水蚤、長刺長腹劍水蚤和線長腹水蚤等8種組成, 其組成較為復(fù)雜, 但種間相似性較低(距離為2.65—5.34)。

    圖2 南沙群島海域長腹劍水蚤種類聚類圖Fig.2 The clusters of concurrent Oithona species groups in Nansha Islands of South China Sea

    從各種出現(xiàn)頻率來看, 組群Ⅰ的羽長腹劍水蚤、擬長腹劍水蚤和瘦長腹劍水蚤是南沙群島海域中全年均有出現(xiàn), 且出現(xiàn)頻率位列前三位的長腹劍水蚤種類。可見, 組群Ⅰ雖然種類少, 但為常見廣布類群。組群Ⅱ的種類出現(xiàn)頻率均較低(F≤26.8%), 除堅(jiān)長腹劍水蚤和細(xì)長腹劍水蚤周年出現(xiàn)外, 其余6種均為季節(jié)性出現(xiàn)種類。因此, 組群Ⅱ?yàn)榈皖l類群。

    2.3 數(shù)量分布

    2.3.1 平面分布 南沙群島海域長腹劍水蚤年均棲息密度為30ind./m3, 僅占浮游動(dòng)物總數(shù)的11.6%,但長腹劍水蚤數(shù)量的季節(jié)變化較小(SD=15.7)。春季數(shù)量為22ind./m3;夏季增至28ind./m3, 秋季急劇上升至53ind./m3, 冬季降至全年最低, 為18ind./m3。

    從長腹劍水蚤總棲息密度的平面分布情況來看,長腹劍水蚤高數(shù)量區(qū)主要出現(xiàn)在近岸海域。春、夏和秋季的分布趨勢(shì)基本相同, 大致呈現(xiàn)西北和東南部高、東北和中部低的趨勢(shì)。除西北部數(shù)量較高外, 在東南部分別還有2個(gè)高數(shù)量區(qū)出現(xiàn)。東北和中部的低密度區(qū)除在夏季連為一片外, 春季和秋季為2個(gè)相對(duì)獨(dú)立的區(qū)域(圖 3)。冬季, 長腹劍水蚤棲息密度呈明顯的北低、南高的趨勢(shì), 海域內(nèi)除東南部的一個(gè)高數(shù)量區(qū)外, 其余區(qū)域數(shù)量均較低(圖3i)。

    春季、夏季和冬季, 南沙群島海域長腹劍水蚤的組群Ⅰ的棲息密度明顯高于組群Ⅱ, 棲息密度比值較為穩(wěn)定(2—2.4∶1)。而秋季, 2個(gè)組群的數(shù)量均明顯增加, 以組群Ⅱ的增幅較大, 其數(shù)量增至與組群Ⅰ相當(dāng)。從組群Ⅰ和組群Ⅱ各季的平面分布狀況來看, 2個(gè)組群各季的分布趨勢(shì)與總棲息密度基本保持一致,無明顯變化(圖3ii、3iii)。

    圖3 南沙群島海域長腹劍水蚤數(shù)量(ind./m3)的平面分布Fig.3 The density of Oithona spp. (ind./m3) in Nansha Islands of South China Sea

    2.3.2 垂直分布 垂直分層只采集到擬長腹劍水蚤、堅(jiān)長腹劍水蚤、羽長腹劍水蚤、瘦長腹劍水蚤、長刺長腹劍水蚤和敏長腹劍水蚤等6種。從種類的出現(xiàn)水層來看, 長腹劍水蚤屬種類組成沿水深變化的現(xiàn)象較為明顯。擬長腹劍水蚤僅出現(xiàn)于 0—30m層,長刺長腹劍水蚤僅出現(xiàn)于 30—75m 層, 堅(jiān)長腹劍水蚤僅出現(xiàn)于75—150m層, 敏長腹劍水蚤僅分布于30—500m 層; 瘦長腹劍水蚤雖然各水層均有出現(xiàn), 但主要分布于150m以淺的水層內(nèi)。而羽長腹劍水蚤則是在各水層均有較多數(shù)量分布的種類, 也是分布在50—500m層的主要種類。

    南沙群島750m以淺海域, 長腹劍水蚤平均棲息密度為 9.6ind./m3, 沿水深梯度的垂直變化明顯(SD=10.2)。0—30m層長腹劍水蚤棲息密度最高, 平均為 21.7ind./m3; 30—75m 降至 19.3ind./m3, 75—150m、150—500m 和 500—750m 的數(shù)量分別為6.1ind./m3、0.6ind./m3和 0.3ind./m3。111°30′E 斷面的平均數(shù)量為 11.7ind./m3, 高于 9°30′N 斷面的9.6ind./m3。

    如圖 4 所示, 9°30′N 和 111°30′E 兩個(gè)斷面長腹劍水蚤密度的垂直分布, 均呈現(xiàn)由淺至深數(shù)量明顯減少的趨勢(shì)。長腹劍水蚤主要集中分布在0—75m的水層中(D0—75m:D75—750m=8.9), 30—75m 其數(shù)量下降趨勢(shì)較為平緩, 75—500m 數(shù)量急劇降低至最低后, 在500—750m基本維持不變。

    2.4 GAM模型

    圖4 南沙群島海域長腹劍水蚤數(shù)量的垂直分布Fig.4 Vertical distribution of Oithona density in Nansha Islands of South China Sea

    表2 GAM模型分析結(jié)果Tab.2 Result of generalized additive models

    GAM模型的分析結(jié)果見表2和圖5。該模型對(duì)長腹劍水蚤密度的總偏差解釋率(deviance explained)為52.5%。t檢驗(yàn)表明, 緯度與長腹劍水蚤棲息密度呈顯著相關(guān)關(guān)系(P<0.01), 溫度和葉綠素a與長腹劍水蚤呈相關(guān)關(guān)系(P<0.05), 鹽度與棲息密度的相關(guān)性一般(P<0.1), 而經(jīng)度與棲息密度沒有相關(guān)性。

    緯度的變化對(duì)長腹劍水蚤數(shù)量變化的影響最為顯著。6°N以南海域長腹劍水蚤數(shù)量較低, 隨緯度的增加數(shù)量迅速上升, 緯度對(duì)數(shù)量的影響也最為顯著; 接近 6°N 時(shí), 緯度對(duì)長腹劍水蚤數(shù)量的影響開始下降; 在 6°N 以北區(qū)域, 緯度對(duì)數(shù)量的影響程度基本不變。6°N和11°N是長腹劍水蚤屬的主要分布區(qū)。海水平均溫度和鹽度分別低于29°C和33時(shí), 隨溫度和鹽度的升高, 其對(duì)長腹劍水蚤數(shù)量的影響程度逐漸增強(qiáng); 而海水平均溫度和鹽度分別高于29°C和33后, 隨溫度和鹽度的繼續(xù)升高, 其對(duì)數(shù)量的影響程度逐漸降低。從影響程度可以得出, 長腹劍水蚤屬最適溫度為 28.6—29.2°C、最適鹽度為 32.6—33.2、最適葉綠素a濃度為 510μg/L(圖 5)。

    圖5 各影響因子與南沙群島海域長腹劍水蚤數(shù)量關(guān)系的GAM圖Fig.5 Relationship between impact factors and density of Oithona in Nansha Islands of South China Sea based on generalized additive model

    3 討論

    3.1 長腹劍水蚤數(shù)量及其季節(jié)變化與環(huán)境因子之間的關(guān)系

    從GAM分析結(jié)果來看, 緯度、溫度、鹽度和葉綠素a濃度均與長腹劍水蚤數(shù)量有相關(guān)關(guān)系。雖然緯度的相關(guān)性最為顯著, 但在 6°N 以北的大部分區(qū)域,緯度對(duì)數(shù)量的影響程度基本不變。因此, 以溫度、鹽度和葉綠素a濃度為重點(diǎn), 討論南沙群島海域長腹劍水蚤數(shù)量及其季節(jié)變化與環(huán)境因子之間的關(guān)系。

    春季, 長腹劍水蚤主要棲息水層的溫度和葉綠素 a濃度在全年中最低, 而鹽度最高, 分別為28.21°C、0.13μg/L 和 33.19。長腹劍水蚤高數(shù)量區(qū)(≥30ind./m3)以相對(duì)低溫、低鹽和低生產(chǎn)力為主要特征。夏季, 溫度達(dá)到全年最高值(29.60°C), 全區(qū)海水溫度超過長腹劍水蚤的最適溫度(≥29.2°C)。但隨水溫的升高和鹽度的降低(33.14), 海域生產(chǎn)力水平明顯有所提高(Meanchla=12.54μg/L), 長腹劍水蚤數(shù)量也出現(xiàn)小幅上升。高數(shù)量區(qū)(≥30ind./m3)仍以相對(duì)低溫、低鹽和低生產(chǎn)力為主要特征。秋季是海域生產(chǎn)力水平最高的季節(jié)(Meanchla=16.99μg/L)。此時(shí), 海水溫鹽降至 28.84°C和 32.73, 進(jìn)入長腹劍水蚤的最適溫鹽區(qū)內(nèi), 而且全海域 91%的區(qū)域均在長腹劍水蚤的最適溫鹽范圍之內(nèi)。長腹劍水蚤的數(shù)量也急劇上升, 達(dá)到全年最高值。高數(shù)量區(qū)(≥60ind./m3)以相對(duì)低鹽和高生產(chǎn)力為主要特征。就長腹劍水蚤的 2個(gè)組群來看,春、夏和冬季, 組群Ⅱ的數(shù)量均為組群Ⅰ的 50%左右。但在環(huán)境最為適宜的秋季, 組群Ⅱ的數(shù)量急劇上升至與組群Ⅰ相當(dāng)??梢? 廣布的組群Ⅰ對(duì)環(huán)境的適應(yīng)性較強(qiáng), 而低頻的組群Ⅱ?qū)Νh(huán)境的敏感程度較高。冬季, 海水溫鹽和葉綠素a濃度均明顯降低, 分別為28.27°C、32.69和 15.34μg/L, 除鹽度外, 溫度和葉綠素 a濃度都遠(yuǎn)不在南沙群島海域長腹劍水蚤的最適范圍內(nèi)。各季中, 冬季全區(qū)溫鹽的變幅最大(SDT=0.76、SDS=0.33), 長腹劍水蚤數(shù)量急劇降至全年最低值, 海域內(nèi)絕大部分區(qū)域的數(shù)量低于30ind./m3。高數(shù)量區(qū)(≥30ind./m3)以相對(duì)高溫和低鹽為主要特征。

    南沙群島海域各季長腹劍水蚤最高數(shù)量區(qū)主要環(huán)境特征的季節(jié)差異明顯: 春季該區(qū)域的溫度是全海域最低(27.03°C), 夏季葉綠素a濃度是全海域次低(4.42μg/L), 秋、冬季海水鹽度均為全海域次低(32.31)??梢? 春、夏和秋冬季, 南沙群島海域制約長腹劍水蚤分布的重要因素分別為溫度、葉綠素a濃度和鹽度。同時(shí)也反映出南沙群島海域的長腹劍水蚤屬更適于低溫、低鹽和低生產(chǎn)力的環(huán)境, 與GAM分析圖所呈現(xiàn)的結(jié)果一致(圖5)。

    雖然在秋季長腹劍水蚤屬高數(shù)量區(qū)的葉綠素 a濃度也較高, 這與秋季的環(huán)境最為適宜有關(guān)。但從長腹劍水蚤的總體分布情況來看, 其更傾向于低葉綠素a的環(huán)境中。首先, 與長腹劍水蚤食性有關(guān)。作為小型橈足類, 長腹劍水蚤是海洋食物網(wǎng)中的中間捕食者, 以單細(xì)胞生物和橈足類幼蟲為食, 也是微食物環(huán)中的食屑者(Hwang et al, 2010), 不以浮游植物為食; 其次, 長腹劍水蚤是大型肉食性浮游動(dòng)物、仔魚和其他食浮游動(dòng)物生物的餌料。葉綠素a濃度高的高生產(chǎn)力區(qū)域, 也是小型橈足類攝食者的高數(shù)量區(qū),對(duì)長腹劍水蚤的捕食壓力大, 故長腹劍水蚤的數(shù)量也較低。

    3.2 長腹劍水蚤數(shù)量平面分布與季風(fēng)和海流的關(guān)系

    長腹劍水蚤與環(huán)境因子的GAM分析表明, 其更適于低溫和低鹽的環(huán)境, 反映出南沙群島海域的長腹劍水蚤與低溫低鹽的沿岸流有密切的關(guān)系。南沙海域上層(0—400m)主要有西部沿岸流、東部沿岸流、北南沙海流和南沙逆(風(fēng))流 4支流系(方文東等,1997)。而南沙群島海域長腹劍水蚤屬各季的高數(shù)量,均出現(xiàn)于受西部沿岸流和東部沿岸流影響的海域內(nèi)。

    浮游生物自主游動(dòng)能力弱, 受海流影響顯著。季風(fēng)是南海表層流場(chǎng)的主要驅(qū)動(dòng)力。南沙群島海域受季風(fēng)影響強(qiáng)烈, 該海域上層4支海流, 除東部沿岸流外,其他 3支主要流系均隨季風(fēng)逆轉(zhuǎn)而倒轉(zhuǎn)(方文東等,1997)。因此, 在受東部沿岸流影響的加里曼丹島東北近岸海域, 各季均會(huì)出現(xiàn)一個(gè)長腹劍水蚤的高數(shù)量區(qū)。春季是東北-西南季風(fēng)轉(zhuǎn)換期。此時(shí), 西部沿岸流受北南沙反氣旋影響, 其影響局限在海區(qū)西北部(劉巖松等, 2014), 該區(qū)域內(nèi)的長腹劍水蚤數(shù)量較高; 夏季盛行西南季風(fēng), 在其作用下西部沿岸流流幅變窄(薛惠潔等, 2001), 大部分海區(qū)受南沙上層反氣旋環(huán)流影響顯著(方文東等, 1997), 長腹劍水蚤的高數(shù)量區(qū)向西北方向有所回縮; 秋季是西南—東北季風(fēng)轉(zhuǎn)換期, 在南海中部氣旋式次海盆環(huán)流的作用下, 南沙西部沿岸流(越南沿岸流)為南向, 并形成東向的越南離岸流, 跨越至菲律賓沿岸(巴拉望島和呂宋島沿岸)(劉巖松等, 2014), 南沙大部分海域在其影響之下,長腹劍水蚤的數(shù)量急劇上升, 高數(shù)量區(qū)范圍也有明顯擴(kuò)大; 冬季, 在強(qiáng)東北風(fēng)作用下, 南?;旧嫌梢粴庑江h(huán)流控制, 整個(gè)南海海水表層溫度因受環(huán)流作用而明顯降低。此外, 強(qiáng)東北風(fēng)作用還使上層海水在南海海盆的西南角堆積。這種堆積在西南角產(chǎn)生一個(gè)正壓場(chǎng), 導(dǎo)致強(qiáng)的沿岸流, 再加上 15°N 以南風(fēng)應(yīng)力的作用下, 出現(xiàn)沿西邊界向南的高流速邊界流(劉秦玉等, 1997)。故在南沙群島海域東北部的來自南海中部的上層海水堆積區(qū)內(nèi), 長腹劍水蚤數(shù)量較低。由于正壓場(chǎng)和風(fēng)應(yīng)力的共同作用, 南沙群島海域沿岸流的勢(shì)力增強(qiáng), 控制了南沙西南部的大部分海域, 該區(qū)域內(nèi)長腹劍水蚤密度明顯高于上層海水堆積區(qū)。因此, 南沙群島海域長腹劍水蚤數(shù)量分布的季節(jié)變化與季風(fēng)驅(qū)動(dòng)的沿岸流勢(shì)力強(qiáng)弱密切相關(guān)。

    3.3 長腹劍水蚤數(shù)量的垂直變化

    1984年章淑珍和陳清潮指出, 南海浮游橈足類的垂直分布是不均勻的, 南海北部近海、北部陸架區(qū)、大陸坡和中部深海盆海域100m水深是橈足類垂直分布的主要轉(zhuǎn)折點(diǎn), 橈足類主要密集在100m以上水層, 而在100m至200m水層中, 其數(shù)量驟然下降;200m以深數(shù)量遞減較為緩慢(章淑珍等, 1984)。南沙西南大陸斜坡區(qū)橈足類主要密集出現(xiàn)在100m以上水深, 75—100m處其數(shù)量大幅下降, 75—100m水深是南海大部分海域橈足類垂直分布的主要轉(zhuǎn)折點(diǎn)(杜飛雁等, 2014)。從南沙群島長腹劍水蚤的垂直分布情況來看, 數(shù)量主要集中在 0—75m層, 75m以深數(shù)量急劇下降, 符合南海橈足類垂直分布規(guī)律。雖然幾次調(diào)查的分層有所差異, 但南沙群島海域長腹劍水蚤數(shù)量的垂直變化規(guī)律, 與南海西南大陸斜坡區(qū)、南沙群島東南部和南海中部海域浮游動(dòng)物一致, 均呈 0—100m層最高, 越往深層量值越低的分布趨勢(shì)(陳瑞祥等, 1988; 陳清潮等, 1991; 杜飛雁等, 2014)。

    南沙群島海域終年存在的溫躍層是影響浮游生物垂直分布的重要因素之一(邱章等, 1994; 杜飛雁等,2014), 溫躍層對(duì)浮游動(dòng)物的垂直分布具有阻隔作用(鄭重等, 1965; 李榮國等, 2004), 溫躍層的水溫變化梯度太大, 很少種類能穿過溫躍層上下移動(dòng)(Fragopoulo et al, 1990)。我國南黃海夏季和南海西南大陸斜坡區(qū)春季的浮游動(dòng)物, 也主要分布于溫躍層中(左濤等, 2004; 杜飛雁等, 2014)。南沙群島海域水深75m左右是溫躍層的下方(中國科學(xué)院南沙綜合科學(xué)考察隊(duì), 1989), 現(xiàn)場(chǎng)同步實(shí)測(cè)數(shù)據(jù)也表明海水溫度在 75m層降幅明顯增大, 而南沙群島長腹劍水蚤的數(shù)量, 也在 75m以深數(shù)量急劇下降。綜上所述表明,南沙群島長腹劍水蚤主要分布于溫躍層中, 受其阻隔作用, 長腹劍水蚤在溫躍層下方的漸變層內(nèi)數(shù)量較低, 溫躍層是影響其數(shù)量垂直分布的主要因素。

    致謝 樣品鑒定得到了中國科學(xué)院南海海洋研究所陳清潮研究員、廈門大學(xué)黃加祺教授和許振祖教授、國家海洋局第三海洋研究所陳瑞祥研究員和林茂研究員的幫助, 美國緬因大學(xué)陳勇教授、Kisei Tanaka和李增光在GAM模型應(yīng)用方面提供指導(dǎo)和幫助, 中國水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所黃洪輝研究員提供環(huán)境數(shù)據(jù), 課題組其他同志在樣品采集過程給予協(xié)作和幫助, 特致謝忱。

    中國科學(xué)院南沙綜合科學(xué)考察隊(duì), 1989. 南沙群島及其鄰近海區(qū)綜合調(diào)查研究報(bào)告(一), 下卷. 北京: 科學(xué)出版社, 392

    牛明香, 李顯森, 徐玉成, 2010. 基于廣義可加模型的時(shí)空和環(huán)境因子對(duì)東南太平洋智利竹 筴 魚漁場(chǎng)的影響. 應(yīng)用生態(tài)學(xué)報(bào), 21(4): 1049—1055

    方文東, 2000. 南海南部海洋環(huán)流的結(jié)構(gòu)與季節(jié)變化. 青島:中國科學(xué)院海洋研究所博士學(xué)位論文, 1—25

    方文東, 郭忠信, 黃羽庭, 1997. 南海南部海區(qū)的環(huán)流觀測(cè)研究. 科學(xué)通報(bào), 42(21): 2264—2271

    左 濤, 王 榮, 王 克等, 2004. 夏季南黃海浮游動(dòng)物的垂直分布與晝夜垂直移動(dòng). 生態(tài)學(xué)報(bào), 24(3): 524—530

    劉 佳, 黃清輝, 李建華, 2009. 崇明北湖葉綠素a濃度與環(huán)境因子的GAM回歸分析. 中國環(huán)境科學(xué), 29(12): 1291—1295

    劉巖松, 于 非, 刁新源等, 2014. 越南離岸流跨海盆特征初步分析. 海洋科學(xué), 38(7): 95—102.

    連光山, 林玉輝, 蔡秉及等, 1990. 大亞灣浮游動(dòng)物群落的特征.見: 莊一純主編. 大亞灣海洋生態(tài)文集(II). 北京: 海洋出版社, 180—186

    劉秦玉, 李 薇, 徐啟春, 1997. 東北季風(fēng)與南海海洋環(huán)流的相互作用. 海洋與湖沼, 28(5): 493—502

    劉瑞玉, 2008. 中國海洋生物名錄. 北京: 科學(xué)出版社, 624

    杜飛雁, 王雪輝, 谷陽光等, 2014. 南沙群島西南大陸斜坡海域浮游動(dòng)物的垂直分布. 海洋學(xué)報(bào), 36(6): 94—103

    李冠國, 范振剛, 2004. 海洋生態(tài)學(xué). 北京: 高等教育出版社,30—31

    楊海軍, 劉秦玉, 1998. 南海海洋環(huán)流研究綜述. 地球科學(xué)進(jìn)展, 13(4): 364—368

    邱 章, 黃企洲, 1994. 南沙群島海區(qū)溫躍層時(shí)空分布的分析.見: 南沙群島海區(qū)物理海洋學(xué)研究論文集 I. 北京: 海洋出版社, 64—80

    邱永松, 1990. 南海北部大陸架魚種組群初步分析. 水產(chǎn)學(xué)報(bào),14(4): 267—276

    陳清潮, 2003. 南沙群島海區(qū)生物多樣性名典. 北京: 科學(xué)出版社, 141

    陳清潮, 張谷賢, 尹健強(qiáng)等, 1991. 1988年夏季南沙群島海區(qū)東部和南部浮游動(dòng)物的分布. 見: 南沙群島及其鄰近海區(qū)海洋生物研究論文集(二). 北京: 海洋出版社, 186—193

    陳瑞祥, 蔡秉及, 林 茂等, 1988. 南海中部海域浮游動(dòng)物的垂直分布. 海洋學(xué)報(bào), 10(3): 337—341

    鄭 重, 鄭執(zhí)中, 王 榮等, 1965. 煙、威鮐魚漁場(chǎng)及鄰近水域浮游動(dòng)物生態(tài)的初步研究. 海洋與湖沼, 7(4): 329—354

    唐 浩, 許柳雄, 陳新軍等, 2013. 基于GAM模型研究時(shí)空及環(huán)境因子對(duì)中西太平洋鰹魚漁場(chǎng)的影響. 海洋環(huán)境科學(xué),32(4): 518—522

    章淑珍, 陳清潮, 1984. 南海浮游橈足類的生態(tài)研究I. 數(shù)量的分布和變化. 熱帶海洋, 3(1): 46—55

    薛惠潔, 柴 扉, 徐丹亞等, 2001. 南海沿岸流特征及其季節(jié)變化. 見: 中國海洋學(xué)文集-南海海流數(shù)值計(jì)算及中尺度特征研究. 北京: 海洋出版社, 64—74

    Bigelow K A, Boggs C H, He X, 1999. Environmental effects on swordfish and blue shark catch rates in the US North Pacific longline fishery. Fisheries Oceanography, 8(3): 178—198

    Chang J H, Chen Y, Holland D et al, 2010. Estimating spatial distribution of American lobster Homarus americanus using habitat variables. Marine Ecology Progress Series, 420: 145—156

    Fragopoulu N, Lykakis J J, 1990. Vertical distribution and nocturnal migration of zooplankton in relation to the development of the seasonal thermocline in Patraikos Gulf.Marine Biology, 104(3): 381—387

    Hastie T J, Tibshirani R J, 1990. Generalized Additive Models.London: Chapman and Hall

    Hwang J S, Kumar R, Dahms H U et al, 2010. Interannual,seasonal, and diurnal variations in vertical and horizontal distribution patterns of 6 Oithona spp. (Copepoda:Cyclopoida) in the South China Sea. Zoological Studies,49(2): 220—229

    Li B, Cao J, Chang J H et al, 2015. Evaluation of effectiveness of fixed-station sampling for monitoring American lobster settlement. North American Journal of Fisheries Management, 35(5): 942—957

    Maravelias C D, Reid D G, 1997. Identifying the effects of oceanographic features and zooplankton on prespawning herring abundance using generalized additive models.Marine Ecology Progress Series, 147(1—3): 1—9

    Maravelias C D, 2001. Habitat associations of Atlantic Herring in the Shetland area: influence of spatial scale and geographic segmentation. Fisheries Oceanography, 10(3): 259—267

    Nakamura Y, Turner J T, 1997. Predation and respiration by the small cyclopoid copepod Oithona similisr: how important is feeding on ciliates and heterotrophic flagellates? Journal of Plankton Research, 19(9): 1275—1288

    Paffenh?fer G A, 1993. On the ecology of marine cyclopoid copepods (Crustacea, Copepoda). Journal of Plankton Research, 15(1): 37—55

    Saiz E, Calbet A, Broglio E, 2003. Effects of small-scale turbulence on copepods: the case of Oithona davisae.Limnology and Oceanography, 48(3): 1304—1311

    Turner J T, 2004. The importance of small planktonic copepods and their roles in pelagic marine food webs. Zoological Studies, 43(2): 255—266

    Wood S N, 2006. Generalized Additive Models: An Introduction with R. Boca Raton, FL: Chapman and Hall, CRC Press

    猜你喜歡
    組群水蚤南沙群島
    73個(gè)傳統(tǒng)建筑組群組團(tuán)出道!帶你活進(jìn)從前的慢時(shí)光
    水蚤心率觀測(cè)的實(shí)驗(yàn)條件探究
    北方水庫水季節(jié)性浮游生物抑制方法與實(shí)踐
    遼寧化工(2021年5期)2021-06-03 05:15:50
    指向生命觀念培養(yǎng)的初中生物探究性實(shí)驗(yàn)設(shè)計(jì)
    ——以“不同咖啡因飲料對(duì)水蚤心率的影響”為例
    “組群”“妙比”“知人”:小學(xué)語文古詩群文閱讀的三個(gè)途徑
    我國美濟(jì)礁引種植物調(diào)查
    磁盤組群組及iSCSI Target設(shè)置
    論南沙群島的群島地位
    法治研究(2016年4期)2016-12-01 03:41:30
    水蚤和水蚯蚓的培育
    《美麗的南沙群島》(第二課時(shí))教學(xué)設(shè)計(jì)
    亚洲性夜色夜夜综合| 99久久综合免费| 国产成人精品在线电影| 飞空精品影院首页| 99精国产麻豆久久婷婷| 精品国产超薄肉色丝袜足j| 色老头精品视频在线观看| 中文字幕av电影在线播放| 国产精品秋霞免费鲁丝片| 午夜久久久在线观看| 91精品伊人久久大香线蕉| 九色亚洲精品在线播放| 另类精品久久| 国产成+人综合+亚洲专区| 欧美黄色淫秽网站| 一本色道久久久久久精品综合| 18禁裸乳无遮挡动漫免费视频| 日日夜夜操网爽| 欧美精品亚洲一区二区| 一区二区三区激情视频| av电影中文网址| 成年人午夜在线观看视频| 国产视频一区二区在线看| 又黄又粗又硬又大视频| 老司机靠b影院| 纯流量卡能插随身wifi吗| 亚洲伊人色综图| 少妇被粗大的猛进出69影院| 国产精品久久久人人做人人爽| 19禁男女啪啪无遮挡网站| 久久毛片免费看一区二区三区| 亚洲成av片中文字幕在线观看| 国产一卡二卡三卡精品| 久久久水蜜桃国产精品网| 午夜福利在线免费观看网站| 久久精品国产亚洲av高清一级| 久9热在线精品视频| 99精品久久久久人妻精品| 久久久久精品国产欧美久久久 | 桃红色精品国产亚洲av| 国产黄频视频在线观看| 国产亚洲午夜精品一区二区久久| 久久热在线av| 精品一区在线观看国产| 国产成人影院久久av| 亚洲精品日韩在线中文字幕| 热99re8久久精品国产| av片东京热男人的天堂| 久久久国产一区二区| 成人三级做爰电影| 青春草亚洲视频在线观看| av国产精品久久久久影院| 一级毛片精品| 国产免费一区二区三区四区乱码| 亚洲va日本ⅴa欧美va伊人久久 | 国产欧美日韩精品亚洲av| 亚洲免费av在线视频| 国产精品国产三级国产专区5o| 爱豆传媒免费全集在线观看| a在线观看视频网站| 男女之事视频高清在线观看| 亚洲成av片中文字幕在线观看| www.av在线官网国产| 久久女婷五月综合色啪小说| 亚洲人成77777在线视频| 美女高潮到喷水免费观看| 1024香蕉在线观看| 69av精品久久久久久 | 天天操日日干夜夜撸| 久久久久国内视频| av福利片在线| 爱豆传媒免费全集在线观看| 人人澡人人妻人| 国产一区二区 视频在线| 人妻 亚洲 视频| 老熟女久久久| 欧美激情极品国产一区二区三区| 国产免费福利视频在线观看| 又大又爽又粗| 精品免费久久久久久久清纯 | 亚洲男人天堂网一区| 18禁国产床啪视频网站| 在线天堂中文资源库| 美女高潮到喷水免费观看| 黄频高清免费视频| 亚洲国产欧美一区二区综合| 色精品久久人妻99蜜桃| 波多野结衣一区麻豆| 免费不卡黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 可以免费在线观看a视频的电影网站| 夫妻午夜视频| 一级毛片女人18水好多| 丰满少妇做爰视频| 国产在视频线精品| 又紧又爽又黄一区二区| 少妇粗大呻吟视频| 丁香六月天网| 欧美午夜高清在线| 国产一区二区三区在线臀色熟女 | 人成视频在线观看免费观看| 久久人妻福利社区极品人妻图片| 亚洲精品自拍成人| 在线观看免费午夜福利视频| 18禁裸乳无遮挡动漫免费视频| 久久国产精品人妻蜜桃| 91精品伊人久久大香线蕉| 侵犯人妻中文字幕一二三四区| 老熟妇仑乱视频hdxx| 精品一区二区三卡| 女人久久www免费人成看片| 美女大奶头黄色视频| 国产国语露脸激情在线看| 免费av中文字幕在线| 99国产精品一区二区三区| 天堂俺去俺来也www色官网| 国产日韩一区二区三区精品不卡| 免费看十八禁软件| 精品视频人人做人人爽| 欧美日韩亚洲高清精品| 午夜福利视频在线观看免费| 国产日韩欧美在线精品| 亚洲精品美女久久av网站| 亚洲九九香蕉| 香蕉国产在线看| 可以免费在线观看a视频的电影网站| 国产成人av激情在线播放| 一本久久精品| 国产日韩欧美视频二区| 不卡一级毛片| 伊人久久大香线蕉亚洲五| 亚洲成人免费av在线播放| 精品人妻在线不人妻| 性高湖久久久久久久久免费观看| 天堂中文最新版在线下载| 免费观看a级毛片全部| 亚洲综合色网址| 亚洲欧美日韩高清在线视频 | www.精华液| 一级毛片女人18水好多| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲精品一二三| 成人免费观看视频高清| 欧美xxⅹ黑人| 免费观看a级毛片全部| 国产人伦9x9x在线观看| 国产免费视频播放在线视频| 下体分泌物呈黄色| 精品福利观看| 在线观看www视频免费| 男女午夜视频在线观看| 国产精品久久久人人做人人爽| 国产一区二区激情短视频 | 一本综合久久免费| 韩国高清视频一区二区三区| 美女中出高潮动态图| 12—13女人毛片做爰片一| 97人妻天天添夜夜摸| 亚洲专区国产一区二区| 超碰成人久久| 国产亚洲av高清不卡| 免费日韩欧美在线观看| 热re99久久精品国产66热6| 在线观看人妻少妇| 久久亚洲国产成人精品v| 丝袜美腿诱惑在线| 日本av免费视频播放| 日韩中文字幕视频在线看片| 精品国产乱子伦一区二区三区 | 90打野战视频偷拍视频| 十八禁人妻一区二区| 日韩一卡2卡3卡4卡2021年| 男女国产视频网站| 91九色精品人成在线观看| 在线观看www视频免费| 日韩视频在线欧美| 日本精品一区二区三区蜜桃| 狠狠狠狠99中文字幕| 国产欧美日韩综合在线一区二区| 丰满人妻熟妇乱又伦精品不卡| 精品久久久精品久久久| 三上悠亚av全集在线观看| 人人妻,人人澡人人爽秒播| 欧美午夜高清在线| 欧美日本中文国产一区发布| 丝袜美腿诱惑在线| 男人操女人黄网站| 人妻久久中文字幕网| 夫妻午夜视频| av电影中文网址| av免费在线观看网站| 国产精品麻豆人妻色哟哟久久| 深夜精品福利| 亚洲av日韩在线播放| 国产成人一区二区三区免费视频网站| 伦理电影免费视频| 精品少妇内射三级| 国产成人欧美| 久久国产精品男人的天堂亚洲| 视频区欧美日本亚洲| 交换朋友夫妻互换小说| 亚洲伊人久久精品综合| 两个人看的免费小视频| 欧美人与性动交α欧美软件| 性色av一级| 精品人妻在线不人妻| 免费观看av网站的网址| 极品人妻少妇av视频| 亚洲av美国av| 考比视频在线观看| 俄罗斯特黄特色一大片| 中文精品一卡2卡3卡4更新| 成年女人毛片免费观看观看9 | 老司机福利观看| 看免费av毛片| 一级a爱视频在线免费观看| 精品亚洲成a人片在线观看| 国产一区二区激情短视频 | 精品福利观看| 飞空精品影院首页| 午夜福利一区二区在线看| 视频区欧美日本亚洲| 香蕉国产在线看| 免费久久久久久久精品成人欧美视频| 国产成+人综合+亚洲专区| 麻豆乱淫一区二区| 午夜影院在线不卡| 岛国毛片在线播放| 久久性视频一级片| 国产激情久久老熟女| 美国免费a级毛片| 激情视频va一区二区三区| 精品国产一区二区三区四区第35| 国产欧美亚洲国产| 最新的欧美精品一区二区| 国产精品免费大片| 在线十欧美十亚洲十日本专区| 国产精品影院久久| 成年人免费黄色播放视频| 视频在线观看一区二区三区| 一本—道久久a久久精品蜜桃钙片| 精品国产乱码久久久久久小说| 人人妻人人添人人爽欧美一区卜| 欧美精品亚洲一区二区| 一区福利在线观看| 成年人免费黄色播放视频| 国产精品影院久久| 国产一区二区三区综合在线观看| 久久九九热精品免费| 国产麻豆69| 无限看片的www在线观看| 国产又爽黄色视频| 中文字幕人妻丝袜制服| 大型av网站在线播放| 国产精品秋霞免费鲁丝片| 久久ye,这里只有精品| 999久久久精品免费观看国产| kizo精华| 99精品久久久久人妻精品| 啦啦啦视频在线资源免费观看| 97在线人人人人妻| 亚洲久久久国产精品| 制服诱惑二区| 美女福利国产在线| 韩国精品一区二区三区| 黄色视频不卡| 亚洲精华国产精华精| 欧美一级毛片孕妇| 日日夜夜操网爽| 侵犯人妻中文字幕一二三四区| 亚洲情色 制服丝袜| 欧美久久黑人一区二区| 免费在线观看视频国产中文字幕亚洲 | 大香蕉久久网| a在线观看视频网站| 欧美精品av麻豆av| 一二三四在线观看免费中文在| 午夜免费观看性视频| av网站免费在线观看视频| 最新在线观看一区二区三区| 亚洲av片天天在线观看| 美女视频免费永久观看网站| 中文字幕人妻丝袜一区二区| 欧美精品一区二区免费开放| 亚洲精品久久午夜乱码| 国产精品国产三级国产专区5o| 久久久欧美国产精品| 日韩人妻精品一区2区三区| 自拍欧美九色日韩亚洲蝌蚪91| 免费在线观看视频国产中文字幕亚洲 | 亚洲第一青青草原| 亚洲精品久久久久久婷婷小说| 久久精品国产综合久久久| 男女高潮啪啪啪动态图| 黄色a级毛片大全视频| 亚洲av欧美aⅴ国产| 精品国产超薄肉色丝袜足j| 手机成人av网站| 国产亚洲精品一区二区www | 国产精品欧美亚洲77777| 可以免费在线观看a视频的电影网站| 老司机午夜十八禁免费视频| 国产成人av教育| 欧美日韩一级在线毛片| 美女主播在线视频| 日本av免费视频播放| 亚洲国产欧美一区二区综合| 亚洲激情五月婷婷啪啪| 母亲3免费完整高清在线观看| 超色免费av| 亚洲精华国产精华精| 久久99一区二区三区| 最近中文字幕2019免费版| 国产免费av片在线观看野外av| 亚洲视频免费观看视频| 国产91精品成人一区二区三区 | 欧美日韩国产mv在线观看视频| 久久精品国产综合久久久| 下体分泌物呈黄色| 一区二区三区激情视频| 亚洲中文字幕日韩| 亚洲伊人色综图| 国产成人系列免费观看| 国产不卡av网站在线观看| 久久亚洲国产成人精品v| 十八禁网站网址无遮挡| 亚洲第一青青草原| 国内毛片毛片毛片毛片毛片| 午夜老司机福利片| 菩萨蛮人人尽说江南好唐韦庄| 国产成人精品久久二区二区91| 国产亚洲精品久久久久5区| 亚洲黑人精品在线| 巨乳人妻的诱惑在线观看| 久久精品国产a三级三级三级| 国产精品熟女久久久久浪| av不卡在线播放| 大片免费播放器 马上看| 中文字幕精品免费在线观看视频| 国产精品久久久久久精品电影小说| 午夜激情av网站| 韩国高清视频一区二区三区| 日韩制服丝袜自拍偷拍| 青春草视频在线免费观看| 久久天堂一区二区三区四区| 在线精品无人区一区二区三| 亚洲性夜色夜夜综合| 男女之事视频高清在线观看| 69精品国产乱码久久久| 18禁黄网站禁片午夜丰满| 可以免费在线观看a视频的电影网站| av在线播放精品| 一级黄色大片毛片| 欧美av亚洲av综合av国产av| 国产成人精品在线电影| 久久九九热精品免费| 自拍欧美九色日韩亚洲蝌蚪91| av国产精品久久久久影院| 久久久久久免费高清国产稀缺| 亚洲精品久久成人aⅴ小说| 久久国产精品大桥未久av| 纵有疾风起免费观看全集完整版| 国产欧美日韩综合在线一区二区| 国产精品一区二区在线不卡| www日本在线高清视频| 女人被躁到高潮嗷嗷叫费观| 在线永久观看黄色视频| 精品福利永久在线观看| 少妇的丰满在线观看| 亚洲成人手机| 欧美精品高潮呻吟av久久| 亚洲久久久国产精品| 18在线观看网站| 精品一区在线观看国产| 久久精品国产综合久久久| 十八禁高潮呻吟视频| 国产一区二区激情短视频 | 午夜福利影视在线免费观看| 日韩中文字幕欧美一区二区| 一区二区日韩欧美中文字幕| 中亚洲国语对白在线视频| 99久久精品国产亚洲精品| 精品欧美一区二区三区在线| 丝袜美足系列| 老熟妇仑乱视频hdxx| 1024香蕉在线观看| 51午夜福利影视在线观看| 人人妻人人澡人人看| 9色porny在线观看| 美女视频免费永久观看网站| 精品乱码久久久久久99久播| 久久久久久久久久久久大奶| 一本—道久久a久久精品蜜桃钙片| 高潮久久久久久久久久久不卡| 建设人人有责人人尽责人人享有的| 一区二区日韩欧美中文字幕| 最近最新免费中文字幕在线| 久久久久国产一级毛片高清牌| 国产免费一区二区三区四区乱码| 老汉色∧v一级毛片| kizo精华| 一级,二级,三级黄色视频| 午夜久久久在线观看| 亚洲精品第二区| 欧美黑人精品巨大| 亚洲av日韩在线播放| 菩萨蛮人人尽说江南好唐韦庄| 日日爽夜夜爽网站| 国产精品免费视频内射| h视频一区二区三区| 51午夜福利影视在线观看| 操美女的视频在线观看| 亚洲av片天天在线观看| 国产高清视频在线播放一区 | 免费在线观看黄色视频的| 9色porny在线观看| 99热全是精品| 亚洲av美国av| 99香蕉大伊视频| 久久99热这里只频精品6学生| 人妻人人澡人人爽人人| 12—13女人毛片做爰片一| 国产精品二区激情视频| 久久精品亚洲av国产电影网| 国产男人的电影天堂91| 啦啦啦视频在线资源免费观看| 丝袜在线中文字幕| 久热这里只有精品99| 最新的欧美精品一区二区| 一级片'在线观看视频| 国产在视频线精品| 亚洲欧美色中文字幕在线| 五月开心婷婷网| 国产一区二区三区av在线| 亚洲av电影在线观看一区二区三区| 黄色视频在线播放观看不卡| 新久久久久国产一级毛片| 色视频在线一区二区三区| 亚洲色图综合在线观看| 热re99久久精品国产66热6| 亚洲男人天堂网一区| 丝袜人妻中文字幕| 国产精品 欧美亚洲| 欧美国产精品va在线观看不卡| 亚洲成人国产一区在线观看| 成人国产av品久久久| 日韩中文字幕视频在线看片| 午夜激情av网站| 国产免费现黄频在线看| 国产精品 欧美亚洲| 乱人伦中国视频| 麻豆av在线久日| 99九九在线精品视频| 欧美黑人欧美精品刺激| 热99re8久久精品国产| 一边摸一边抽搐一进一出视频| 在线天堂中文资源库| 日韩中文字幕视频在线看片| 好男人电影高清在线观看| 大香蕉久久成人网| 午夜福利视频精品| 欧美在线黄色| 成人黄色视频免费在线看| 岛国在线观看网站| 日本黄色日本黄色录像| 亚洲精品国产色婷婷电影| 91成人精品电影| 中文欧美无线码| 亚洲精品国产av成人精品| 99精品欧美一区二区三区四区| 午夜久久久在线观看| 亚洲伊人色综图| 男人操女人黄网站| 五月天丁香电影| 黄色视频不卡| 熟女少妇亚洲综合色aaa.| 精品人妻一区二区三区麻豆| 欧美精品一区二区免费开放| 午夜福利影视在线免费观看| 大香蕉久久网| 国产片内射在线| 青春草视频在线免费观看| 性少妇av在线| 天堂俺去俺来也www色官网| 国产av精品麻豆| 亚洲成人手机| 国产亚洲欧美在线一区二区| 精品国产乱子伦一区二区三区 | 精品久久蜜臀av无| 国产在线一区二区三区精| 99香蕉大伊视频| 亚洲三区欧美一区| 中文字幕人妻丝袜制服| 亚洲成人手机| 久久久久精品人妻al黑| 精品第一国产精品| 国产av一区二区精品久久| 欧美黄色淫秽网站| 成人国产一区最新在线观看| 又大又爽又粗| 老司机影院成人| 男女之事视频高清在线观看| 在线观看舔阴道视频| 性高湖久久久久久久久免费观看| 国产高清国产精品国产三级| 日韩制服丝袜自拍偷拍| 国产熟女午夜一区二区三区| 亚洲午夜精品一区,二区,三区| 高清视频免费观看一区二区| 熟女少妇亚洲综合色aaa.| 国产在线观看jvid| 国产成人免费观看mmmm| 天堂俺去俺来也www色官网| 美女大奶头黄色视频| 两性午夜刺激爽爽歪歪视频在线观看 | 视频在线观看一区二区三区| 岛国在线观看网站| 美女脱内裤让男人舔精品视频| av在线播放精品| 欧美黑人欧美精品刺激| 伦理电影免费视频| 人妻 亚洲 视频| 高清av免费在线| 亚洲少妇的诱惑av| 亚洲国产av影院在线观看| 国产欧美日韩一区二区三区在线| 制服诱惑二区| 日韩中文字幕欧美一区二区| av超薄肉色丝袜交足视频| 我要看黄色一级片免费的| 女人久久www免费人成看片| 一级毛片精品| 高清视频免费观看一区二区| 成人18禁高潮啪啪吃奶动态图| 欧美变态另类bdsm刘玥| 欧美另类一区| 亚洲欧美色中文字幕在线| 18禁裸乳无遮挡动漫免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产av影院在线观看| 丝瓜视频免费看黄片| 国产一区二区三区av在线| 婷婷色av中文字幕| 欧美亚洲 丝袜 人妻 在线| 亚洲精品成人av观看孕妇| 1024视频免费在线观看| 国产精品一区二区精品视频观看| 天堂中文最新版在线下载| 高清欧美精品videossex| 精品国产乱码久久久久久男人| 一本综合久久免费| 中文字幕另类日韩欧美亚洲嫩草| 国产成+人综合+亚洲专区| 香蕉丝袜av| 日韩免费高清中文字幕av| 国产伦人伦偷精品视频| 午夜福利影视在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| av电影中文网址| 亚洲性夜色夜夜综合| 精品高清国产在线一区| 精品久久久久久久毛片微露脸 | 久久久久久久国产电影| 操出白浆在线播放| 国产一级毛片在线| 美女主播在线视频| 亚洲 国产 在线| 国产成人免费观看mmmm| 波多野结衣av一区二区av| av网站在线播放免费| 国产又爽黄色视频| 国产精品免费大片| 在线永久观看黄色视频| 91大片在线观看| 欧美老熟妇乱子伦牲交| 淫妇啪啪啪对白视频 | 视频在线观看一区二区三区| 国产亚洲av高清不卡| 日韩大码丰满熟妇| 日韩一区二区三区影片| tocl精华| www.自偷自拍.com| 日本vs欧美在线观看视频| tube8黄色片| 欧美一级毛片孕妇| av欧美777| 丰满人妻熟妇乱又伦精品不卡| 国产不卡av网站在线观看| 国产精品一区二区在线不卡| 男女床上黄色一级片免费看| av国产精品久久久久影院| 久久精品人人爽人人爽视色| 国产99久久九九免费精品| 亚洲伊人久久精品综合| 新久久久久国产一级毛片| 18禁裸乳无遮挡动漫免费视频| 国产淫语在线视频| 亚洲欧洲精品一区二区精品久久久| 国产91精品成人一区二区三区 | 亚洲精品久久午夜乱码| 91精品三级在线观看| 国产成人一区二区三区免费视频网站| 国产亚洲一区二区精品| 老鸭窝网址在线观看| 午夜激情av网站| 亚洲精品日韩在线中文字幕| 欧美亚洲日本最大视频资源| 午夜精品久久久久久毛片777| 不卡av一区二区三区| 国产精品自产拍在线观看55亚洲 | 亚洲精品一区蜜桃| 两人在一起打扑克的视频| 国产亚洲av片在线观看秒播厂| 夜夜夜夜夜久久久久| 欧美av亚洲av综合av国产av|