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

    南極夏季宇航員海浮游動(dòng)物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系*

    2021-07-21 11:58:28牟文秀徐志強(qiáng)李超倫
    海洋與湖沼 2021年4期
    關(guān)鍵詞:南大洋橈足類磷蝦

    牟文秀 楊 光 郝 鏘 徐志強(qiáng) 李超倫, 4

    南極夏季宇航員海浮游動(dòng)物群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系*

    牟文秀1, 2, 3楊 光1, 2, 4郝 鏘5徐志強(qiáng)1李超倫1, 2, 3, 4①

    (1. 中國科學(xué)院海洋研究所 青島 266071; 2. 青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室 青島 266237; 3. 中國科學(xué)院大學(xué) 北京 100049; 4. 中國科學(xué)院海洋大科學(xué)研究中心 青島 266071; 5. 自然資源部第二海洋研究所 杭州 310012)

    基于中國第36次南極科學(xué)考察期間(2019年12月至2020年1月)于南極宇航員海使用北太平洋網(wǎng)(0.5 m2, 330 μm)采集的浮游動(dòng)物樣品, 對(duì)浮游動(dòng)物的豐度、群落結(jié)構(gòu)、優(yōu)勢(shì)種種群結(jié)構(gòu)及其與環(huán)境因子之間的關(guān)系進(jìn)行了研究。結(jié)果表明: 2019/2020年南極夏季宇航員海浮游動(dòng)物總豐度在8.93—102.00 ind./m3之間, 最高豐度記錄在宇航員海中部靠北區(qū)域, 浮游動(dòng)物主要由橈足類、磷蝦、毛顎類、浮游被囊類等類群組成。調(diào)查海域橈足類豐度對(duì)浮游動(dòng)物總豐度的貢獻(xiàn)可達(dá)90%以上, 其中、大中型橈足類,,和分布廣泛, 而小型橈足類、、、等在豐度上占比較高。通過聚類分析可將宇航員海浮游動(dòng)物在區(qū)域上劃分為5個(gè)群落, 普里茲灣渦流、威德爾渦流、南極繞極流南界等物理海洋學(xué)過程是宇航員海浮游動(dòng)物群落分布模式的主要控制因素。另外, 生物環(huán)境匹配(Bio-env)分析顯示, 餌料(葉綠素)也對(duì)該海區(qū)浮游動(dòng)物的種群和群落結(jié)構(gòu)具有重要影響。

    群落結(jié)構(gòu); 浮游動(dòng)物; 宇航員海; 南大洋

    浮游動(dòng)物作為食物網(wǎng)中連接初級(jí)生產(chǎn)者和高級(jí)消費(fèi)者的中間環(huán)節(jié)(Mayzaud, 2014), 在南大洋生態(tài)系統(tǒng)的物質(zhì)循環(huán)和能量流動(dòng)過程中起到承上啟下的作用(Schmidt, 2016), 并且通過攝食、排糞、被動(dòng)沉降等過程調(diào)節(jié)顆粒有機(jī)碳從表層向深海輸送(Steinberg, 2017)。浮游動(dòng)物的空間分布受環(huán)境條件以及洋流等因素的影響, 也被認(rèn)為是全球變化的生物指示器(Beaugrand, 2002; Estrada, 2012)。因此, 研究南大洋浮游動(dòng)物的分布、群落結(jié)構(gòu)與環(huán)境因子的關(guān)系及其時(shí)空尺度的變化可以更好地探知極地海洋生態(tài)系統(tǒng)結(jié)構(gòu)及其對(duì)全球變化的響應(yīng)與反饋(Constable, 2014; Pinkerton, 2020)。

    氣候變化背景下的南大洋環(huán)境自從20世紀(jì)60年代以來發(fā)生了顯著變化: 西風(fēng)逐漸增強(qiáng)且其影響區(qū)域極向擴(kuò)張、南大洋環(huán)狀模正位相年份增加、冰架崩塌甚至消失、底層水及表層水淡化、部分海域酸化(Rintoul, 2018)。這些環(huán)境變化已經(jīng)對(duì)南大洋生態(tài)系統(tǒng)產(chǎn)生了深遠(yuǎn)影響(Constable, 2014), 南大洋西南極半島周邊海域自20世紀(jì)中期以來溫度升高、海冰減少(Saba, 2014), 浮游植物群落發(fā)生重組, 大粒級(jí)的硅藻逐漸被小粒級(jí)的隱藻所替代(Montes-Hugo, 2009)。海冰的減少、浮游植物粒級(jí)的變化、南大洋環(huán)狀模正相位年份的增加等導(dǎo)致浮游動(dòng)物優(yōu)勢(shì)類群由磷蝦和橈足類向紐鰓樽等浮游被囊動(dòng)物轉(zhuǎn)變(Atkinson, 2004), 而南極磷蝦分布區(qū)域向高緯度收縮(Atkinson, 2019)。受此影響, 高營養(yǎng)級(jí)捕食者如企鵝、海豹、鯨魚的聚居地和覓食區(qū)域也發(fā)生了變化(Henley, 2019)。這一系列的變化預(yù)示著整個(gè)西南極海域生態(tài)系統(tǒng)結(jié)構(gòu)和功能發(fā)生了顯著的變化(Constable, 2014)。

    南大洋環(huán)境的變化趨勢(shì)在不同的扇區(qū)之間存在很大的區(qū)域性差異(Constable, 2014; Convey, 2019)。溫度升高、海冰減少等現(xiàn)象主要發(fā)生在西南極半島海域, 南大洋其他扇區(qū)如東南極海域海冰變化相對(duì)穩(wěn)定, 羅斯海扇區(qū)在1979年以來甚至呈現(xiàn)了溫度降低和海冰增多的趨勢(shì)(Lenaerts, 2017)。盡管人們從20世紀(jì)初期的時(shí)代就開始進(jìn)行浮游動(dòng)物的采樣和研究, 但這些工作主要集中在南大洋的大西洋扇區(qū)、印度洋普里茲灣海區(qū)及太平洋的羅斯海海區(qū)(Atkinson, 2012; Moriarty, 2013)。南大洋宇航員海(Cosmonaut Sea)位于東南極恩德比地(Enderby Land)的西部, 介于30°—60°E之間(Hunt, 2007)。目前, 宇航員海海域的各方面研究均相對(duì)較少, 我們對(duì)其浮游動(dòng)物的認(rèn)識(shí)主要來自20世紀(jì)80年代末期前蘇聯(lián)針對(duì)大型浮游動(dòng)物的系列航次和2006年夏季澳大利亞的“BROKE-West”航次(Hunt, 2007; Swadling, 2010)。Hunt等(2007)發(fā)現(xiàn)宇航員海大型浮游動(dòng)物群落年際變化受到環(huán)南極流和近岸流等環(huán)流的影響, Swadling等(2010)用網(wǎng)目更精細(xì)的RMT1網(wǎng)(315 μm)對(duì)橈足類等中小型浮游動(dòng)物信息進(jìn)行了補(bǔ)充, 并發(fā)現(xiàn)葉綠素濃度、海冰消退時(shí)間等大尺度海洋過程決定了宇航員海浮游動(dòng)物群落的分布模式。

    在2019年12月至2020年1月中國第36次南極考察期間, 我們使用北太平洋網(wǎng)于宇航員海進(jìn)行了海洋上層(0—200 m)浮游動(dòng)物樣品的采集, 基于樣品鑒定數(shù)據(jù)和環(huán)境因子資料, 對(duì)宇航員海浮游動(dòng)物的群落結(jié)構(gòu)及其與環(huán)境因子的關(guān)系進(jìn)行了分析。本研究的目的在于: (1) 描述該海域浮游動(dòng)物種類組成、分布、豐度、群落結(jié)構(gòu); (2) 分析環(huán)境因子對(duì)浮游動(dòng)物群落及分布的影響; 本文對(duì)數(shù)據(jù)稀少的宇航員海進(jìn)行浮游動(dòng)物數(shù)據(jù)的補(bǔ)充, 為今后該海域生態(tài)系統(tǒng)研究工作的開展提供基礎(chǔ)數(shù)據(jù)。

    1 材料與方法

    1.1 現(xiàn)場(chǎng)取樣

    在中國第36次南極科學(xué)考察期間(CHINARE), 于2019年12月6日至2020年1月6日在宇航員海56個(gè)站位(圖1)利用北太平洋網(wǎng)(0.5 m2, 330 μm)對(duì)浮游動(dòng)物進(jìn)行200 m至表層的拖網(wǎng)。拖網(wǎng)結(jié)束后將樣品立即置于5%的中性福爾馬林溶液中保存。

    圖1 宇航員海采樣站位圖

    1.2 樣品和數(shù)據(jù)分析

    實(shí)驗(yàn)室內(nèi)對(duì)浮游動(dòng)物樣品進(jìn)行鑒定并計(jì)數(shù), 其中, 大型浮游動(dòng)物進(jìn)行全樣計(jì)數(shù), 其余個(gè)體根據(jù)樣品量進(jìn)行1/2至1/8分樣后在解剖鏡下計(jì)數(shù)。四種橈足類優(yōu)勢(shì)種、和按照成體、幼體期CⅠ—CⅢ、CⅣ—CⅤ進(jìn)行計(jì)數(shù), 對(duì)兩種磷蝦和則按照CⅠ—CⅢ、FⅠ—FⅢ和FⅣ—FⅥ進(jìn)行計(jì)數(shù)(Yang, 2011)。每個(gè)樣品至少保證300—500個(gè)個(gè)體被計(jì)數(shù), 每個(gè)站位浮游動(dòng)物各物種的豐度由物種總個(gè)數(shù)除以拖網(wǎng)時(shí)過濾水的體積得到, 單位為ind./m3。

    將通過鑒定計(jì)數(shù)得到的所有浮游動(dòng)物數(shù)據(jù)整合后對(duì)浮游動(dòng)物群落進(jìn)行劃分。豐度數(shù)據(jù)經(jīng)平方根轉(zhuǎn)化后基于Bray-Curtis相似性指數(shù)和組平均的方法進(jìn)行q-型聚類(Field, 1982)。多維多度分析(NMDS)對(duì)聚類分析進(jìn)行重復(fù)(Hunt, 2007)。聚類得到的浮游動(dòng)物群組通過相似百分度(SIMPER)來確定對(duì)群落的劃分起關(guān)鍵作用的物種。

    指示種分析(Ind Val)用于確定各站位分組之間的指示物種(Dufrêne, 1997)。

    指示種()計(jì)算公式如下:

    =A×B×100, (1)

    Bio-env分析用來解釋環(huán)境因子(表層溫度、鹽度、葉綠素、0—200 m水柱平均溫度、鹽度、葉綠素、無冰天數(shù))與浮游動(dòng)物群落劃分的關(guān)系。無冰天數(shù)是指某站位離海冰濃度低于15%的天數(shù), 海冰濃度數(shù)據(jù)由美國國家冰雪數(shù)據(jù)中心(National Snow & Ice Data Center, NSIDC)網(wǎng)站獲得。相關(guān)系數(shù)用p表示, 其數(shù)值介于0—1之間,p越接近1表示生物數(shù)據(jù)與環(huán)境因子的相關(guān)性越大,p=0則表示生物數(shù)據(jù)與環(huán)境因子無相關(guān)性。

    溫度和鹽度數(shù)據(jù)由Seabird 911型CTD現(xiàn)場(chǎng)測(cè)定。葉綠素濃度測(cè)定水樣采集自0、25、50、150和200 m, 取500 mL水樣用GF/F玻璃纖維濾膜過濾, 經(jīng)90%丙酮溶液萃取12—24 h后用Turner 熒光計(jì)測(cè)定。

    多元統(tǒng)計(jì)分析主要使用Primer軟件完成。

    2 結(jié)果

    2.1 環(huán)境因子

    調(diào)查期間, 大部分站位海冰均已消退, 少數(shù)靠近大陸邊緣的站位仍有較多浮冰。在整個(gè)采樣區(qū)域, 0—200 m水體平均溫度范圍為-1.76—0.64 °C, 東北和西南區(qū)域水溫較高, 而西北和東南區(qū)域水溫較低(圖2); 0—200 m水體平均鹽度范圍為34.11—34.46, 中部區(qū)域鹽度最高(圖2); 此外, 0—200 m水體平均葉綠素濃度范圍為0.02—1.00 mg/m3, 葉綠素濃度最大值出現(xiàn)在調(diào)查區(qū)域的西北部C4-03站位。

    圖2 0—200 m水柱平均溫度和鹽度

    2.2 浮游動(dòng)物群落

    浮游動(dòng)物總豐度在8.93—102 ind./m3之間, C6-01站位豐度最高, C3-03站位最低(圖3)。浮游動(dòng)物群落主要由橈足類組成, 其豐度占總豐度的90%以上。聚類分析結(jié)果顯示在相似度為71.3%和62.7%時(shí)浮游動(dòng)物可分為5個(gè)群落(圖4)。

    群落1主要位于采樣區(qū)域的東南部, 0—200 m水體平均溫度最低, 為-1.24 °C, 鹽度處于平均水平為34.32, 0—200 m水體平均葉綠素濃度為0.30 mg/m3(圖4c, 表1)。群落1浮游動(dòng)物平均豐度相對(duì)較低, 為(27.09低, 為平為南) ind./m3, 其指示種為、、sp.和spp. (表2)群落2主要位于采樣區(qū)域離大陸較遠(yuǎn)的北部, 0—200 m水體平均溫度相對(duì)較高, 為0.11 °C, 鹽度平均水平為34.29, 平均葉綠素濃度最高, 為0.36 mg/m3(圖4c, 表1)。群落2浮游動(dòng)物平均豐度最高, 為(80.23均豐度最高, ) ind./m3, 其指示種為、、、、、、、、、和(表2)。

    圖3 2019/2020年夏季和2006年夏季宇航員海浮游動(dòng)物的總豐度及空間分布(改自Swadling et al, 2010)

    圖4 基于浮游動(dòng)物豐度數(shù)據(jù)的聚類分析(a)、多維多度分析(b)及地圖中浮游動(dòng)物群落分組展示(c)

    注: 環(huán)流大致位置參照Orsi等(1995)及Williams等(2010)

    表1 群落1—5所處區(qū)域的環(huán)境特征及浮游動(dòng)物平均豐度(平均值及范圍)

    Tab.1 Environmental variables and average abundance of zooplankton (the mean and range) in communities 1—5

    表2 SIMPER分析顯示對(duì)各群落劃分累積貢獻(xiàn)超過90%的物種組成及豐度

    Tab.2 Average abundances (ind./m3) of zooplankton species that contributed >90% to the similarity within each group based on SIMPER analysis

    注: *表示該物種為此群落的指示種; 橈足類CⅠ—Ⅲ表示橈足類的橈足幼體階段; 磷蝦CⅠ—Ⅲ表示磷蝦原蚤狀幼體階段, FⅠ—Ⅵ表示磷蝦蚤狀幼體階段

    群落3主要位于采樣區(qū)域東北部, 0—200 m水體平均溫度和鹽度均相對(duì)較低, 分別為-0.66 °C和34.28, 平均葉綠素濃度相對(duì)略低, 為0.27 mg/m3(圖4c, 表1)。群落3浮游動(dòng)物平均豐度相對(duì)較高, 為(56.36±10.37)ind./m3, 其指示種為、、、、、、、、、sp.和spp. (表2)。群落4主要位于采樣區(qū)域的西部, 0—200 m水體平均溫度和鹽度分別為-0.30 °C和34.29, 平均葉綠素濃度相對(duì)較低, 為0.24 mg/m3(圖4c, 表1)。群落4浮游動(dòng)物平均豐度為(47.26±19.43) ind./m3, 其指示種為、、和(表2)。群落5主要位于采樣區(qū)域的中部, 0—200 m水體平均溫度和鹽度均最高, 分別為0.30 °C和34.33, 而平均葉綠素濃度最低, 為0.23 mg/m3(圖4c, 表1)。群落5浮游動(dòng)物平均豐度最低, 為(17.20±5.45) ind./m3。群落5中浮游動(dòng)物各物種豐度均相對(duì)較低, 這也是該群落無指示種的原因(表2)。群落5存在兩個(gè)亞組(圖4b, 4c), 其中5b主要包括站位C2-5、C2-9、C3-3和C4-5 (圖4c), 兩個(gè)亞組在物種組成和豐度上無明顯差異, 因此本文只對(duì)群落5作為整體進(jìn)行了物種組成和豐度展示(表2)。

    2.3 橈足類優(yōu)勢(shì)物種的種群結(jié)構(gòu)

    、、和在南大洋不同扇區(qū)廣泛分布, 也在本次所有調(diào)查站位中均被發(fā)現(xiàn)。平均豐度最高(圖5), 為(6.32±5.95)ind./m3, 從種群結(jié)構(gòu)來看, 在各群組成體階段均占有較高的比例, 幼體期CI—CIII、CⅣ—CⅤ次之(圖6);平均豐度最低, 為(1.14±1.47) ind./m3, 從種群結(jié)構(gòu)來看, 在各群組幼體期CI—CIII占有較高的比例, 成體、幼體期CⅣ—CⅤ次之;平均豐度較低, 為(1.48±1.73) ind./m3, 從種群結(jié)構(gòu)來看, 在各群組成體階段均占有較高比例, 幼體期CI—CIII、CⅣ—CⅤ次之;平均豐度較高, 為(4.90±5.15) ind./m3, 從種群結(jié)構(gòu)來看, 在各群組幼體期CI—CIII占有較高的比例, 幼體期CⅣ—CⅤ、成體次之。

    圖5 Calanoides acutus, Metridia gerlachei, Calanus propinquus和Rhincalanus gigas平均豐度(誤差棒為標(biāo)準(zhǔn)差)

    2.4 環(huán)境因子對(duì)群落結(jié)構(gòu)的影響

    生物數(shù)據(jù)與環(huán)境數(shù)據(jù)的Bio-Env分析顯示, 溫度是對(duì)群落聚類結(jié)果起關(guān)鍵作用的單因子(p= 0.330, 表3)。表層溫度、0—200 m水體平均鹽度及0—200 m水體平均葉綠素濃度三個(gè)環(huán)境因子的組合與生物數(shù)據(jù)的相關(guān)性最優(yōu)(p= 0.375, 表3), 較好解釋了群落劃分的合理性。無冰天數(shù)并沒有對(duì)浮游動(dòng)物的群落劃分起到主要解釋作用(表3)。

    3 討論

    3.1 宇航員海浮游動(dòng)物組成及群落結(jié)構(gòu)

    本研究使用的網(wǎng)具(北太平洋網(wǎng), 網(wǎng)口面積0.5 m2, 網(wǎng)目330 μm)與2006年夏季澳大利亞“BROKE-West”航次調(diào)查使用的網(wǎng)具(RMT1網(wǎng), 網(wǎng)口面積1 m2, 網(wǎng)目315 μm)具有相近的網(wǎng)目, 因此數(shù)據(jù)具有較好的可比性來研究宇航員海10余年前后夏季浮游動(dòng)物群落的組成和結(jié)構(gòu)差異。由于本次采樣位點(diǎn)未涉及冰間湖等近岸區(qū)域, 因此以往南大洋印度洋扇區(qū)經(jīng)常報(bào)道的以晶磷蝦為指示種、高豐度的近岸群落在本研究中未出現(xiàn)(Hunt, 2007; Yang, 2011)。與以往南大洋不同海域浮游動(dòng)物群落研究(Hosie, 1994; Hunt, 2007; Pakhomov, 2020)相似的是, 群落劃分很大程度上歸因于浮游動(dòng)物優(yōu)勢(shì)物種豐度尤其是橈足類豐度的差異而不是浮游動(dòng)物物種組成的變化(表2)。對(duì)本研究各群落劃分累積起到90%以上的優(yōu)勢(shì)物種/類群一共23種(表2), 均為南大洋較為常見且廣泛分布的物種(Atkinson, 2012; Pakhomov, 2020)。本研究浮游動(dòng)物群落總平均豐度與2006年夏季澳大利亞“BROKE-West”航次記錄在一個(gè)量級(jí)上(圖3), 但數(shù)值略高于其研究結(jié)果[本研究: 44.41 ind./m3對(duì)比Swadling 等(2010): 23.16 ind./m3], 這種差異可能更大程度上源于南大洋“季節(jié)性”的差異而非“年際間”的差異。本次采樣集中在海冰消退初期的夏季初期(12月), 而“BROKE-West”航次集中在海冰消退中后期的1—2月(Swadling, 2010)。由于南大洋海冰等環(huán)境因子強(qiáng)烈的季節(jié)性差異及上行效應(yīng)在食物網(wǎng)中的主導(dǎo)地位(Atkinson, 2014), 浮游動(dòng)物群落的年內(nèi)差異(季節(jié)差異)通常高于年際間的差異(Pinkerton, 2020)。在物種組成上, 本研究與“BROKE-West”航次有很大的相似性: 橈足類通常較占優(yōu)勢(shì), 其豐度對(duì)該海域浮游動(dòng)物總豐度的貢獻(xiàn)可達(dá)90%以上(表2, Swadling, 2010)。盡管我們采樣網(wǎng)具的網(wǎng)目(330 μm)對(duì)、、、等小型橈足類存在顯著的低估(Makabe, 2012), 它們?cè)趦蓚€(gè)航次均呈現(xiàn)極高的豐度(表2)。隨著近些年南大洋采樣網(wǎng)具的精細(xì)化, 小型橈足類在初級(jí)生產(chǎn)的攝食控制、顆粒有機(jī)碳向深水的傳輸?shù)壬鷳B(tài)系統(tǒng)過程中的貢獻(xiàn)越來越受到人們的關(guān)注(Calbet, 2008; Steinberg, 2017)。

    圖6 Calanoides acutus, Metridia gerlachei, Calanus propinquus和Rhincalanus gigas各發(fā)育期占比

    表3 Bio-env分析結(jié)果

    Tab.3 Results of Bio-Env analysis

    、、和四種優(yōu)勢(shì)大中型橈足類是南大洋浮游動(dòng)物生物量的主要貢獻(xiàn)類群, 在某些海域其產(chǎn)量甚至遠(yuǎn)高于南極磷蝦(Conover, 1991; Voronina, 1998; Pakhomov, 2020)。它們存在不同程度的晝夜垂直移動(dòng)和季節(jié)性垂直遷移(Atkinson, 1992)。儲(chǔ)存脂質(zhì)是高緯度地區(qū)橈足類的典型特征, 大多數(shù)物種在夏季富含浮游植物的表層水體攝食后進(jìn)入深水區(qū)進(jìn)行越冬(Atkinson, 1998), 其中冬季滯育物種的脂質(zhì)含量相對(duì)最高(Donnelly, 1994; Hagen, 1996)。這種由橈足類季節(jié)性垂直移動(dòng)主導(dǎo)的碳向深海傳輸?shù)闹|(zhì)泵過程首次在北大西洋提出(Jónasdóttir, 2015), 并在近些年的生物地球化學(xué)循環(huán)過程研究中廣泛體現(xiàn)(Boyd, 2019)。在本研究中, 0—200 m水體中四種優(yōu)勢(shì)橈足類早期幼體CI—CIII數(shù)量較多(圖6), 說明它們正處于夏季浮游植物暴發(fā)后的種群補(bǔ)充期。種群呈現(xiàn)最高的早期幼體比例(圖6), 與以往研究認(rèn)為生活史過程中具有延長的發(fā)育期的觀點(diǎn)相一致(Schnack-Schiel,1994)。

    本次研究我們?cè)谟詈絾T海主要采集到兩種磷蝦: 南極磷蝦和長臂櫻磷蝦(表2)。與“BROKE-West”航次結(jié)果相似的是, 長臂櫻磷蝦在宇航員海具有更高的豐度(Swadling, 2010)。需要指出的是由于成體磷蝦較強(qiáng)的游泳能力, 我們的調(diào)查網(wǎng)具與Swadling等(2010)的調(diào)查網(wǎng)具均不能有效地對(duì)其進(jìn)行獲取。南極磷蝦在位于大洋區(qū)的群落2、4相對(duì)較高的豐度(表2)驗(yàn)證了Atkinson等(2008)提出的南大洋印度-太平洋扇區(qū)南極磷蝦主要分布在陸架坡折附近大洋區(qū)的觀點(diǎn)。近期基于單季節(jié)環(huán)南極調(diào)查航次和KRILLBASE數(shù)據(jù)庫(1926—2016年)的研究結(jié)果顯示隨著南大洋大西洋扇區(qū)磷蝦生境的破壞, 其他扇區(qū)成為磷蝦生境的避難所(Yang, 2021)。宇航員海將來是否可以成為南極磷蝦新的生物量熱點(diǎn)區(qū)還需要今后持續(xù)的調(diào)查和采樣。

    3.2 環(huán)境因子對(duì)浮游動(dòng)物群落的影響

    水團(tuán)、鋒面、溫度、餌料供應(yīng)是南大洋浮游動(dòng)物群落及生態(tài)系統(tǒng)結(jié)構(gòu)的主要影響因子(Nicol, 2000; Chapman, 2020)。本研究中5組群落的劃分與水團(tuán)有很好的吻合(圖4)?;?006年澳大利亞“BROKE-West”航次(Meijers, 2010)及本航次的溫鹽數(shù)據(jù), 群落1和群落3靠近順時(shí)針循環(huán)的普里茲灣渦流(Prydz Bay Gyre); 群落2靠近南極繞極流并且在南極繞極流南部邊界以南, 位于更深、更溫暖的水域; 群落4主要受到順時(shí)針循環(huán)進(jìn)入宇航員海并隨南極陸坡流(Antarctic Slope Current)返回的威德爾渦流(Weddell Gyre)的影響; 群落5主要位于采樣區(qū)域的中部, 該群落0—200 m水體溫鹽均最高, 可能是受到威德爾環(huán)流和南極繞極流南部邊界的影響。Hunt等(2007)對(duì)該海域大型浮游動(dòng)物群落年際變化的研究也顯示此區(qū)域浮游動(dòng)物群落與環(huán)南極流、近岸流等水團(tuán)分布有顯著的關(guān)系。位于順時(shí)針循環(huán)普里茲灣渦流內(nèi)的群落1具有最高的南極磷蝦豐度(表2), Kawaguchi等(2010)認(rèn)為此渦流體系存在獨(dú)立的、自我維持的磷蝦種群。

    海冰是極地海域重要的環(huán)境因子之一, 南大洋海冰的變化通過餌料供應(yīng)引起的食物網(wǎng)生態(tài)學(xué)效應(yīng)受到人們的廣泛關(guān)注(Atkinson, 2004, 2019; Constable, 2014)。Swadling等(2010)發(fā)現(xiàn)海冰消退的天數(shù)很好地解釋了宇航員浮游動(dòng)物群落的劃分。但本文研究結(jié)果顯示海冰相關(guān)的參數(shù)(無冰天數(shù)等)在Bio-Env分析結(jié)果中沒有得以體現(xiàn), 我們分析浮游動(dòng)物群落結(jié)構(gòu)主要受到環(huán)流、鋒面等相對(duì)大尺度的物理海洋學(xué)過程控制。另外海冰消退通過對(duì)餌料供應(yīng)的效應(yīng)影響浮游動(dòng)物的種群補(bǔ)充需要一定的時(shí)間積累(Yang, 2011), 本航次調(diào)查期間大部分站位均處于海冰消退初期, 而BROKE-WEST航次相比此航次滯后一個(gè)月, 且絕大部分站位均處于海冰消退后期(Swadling, 2010), 海冰消退通過對(duì)餌料供應(yīng)的上行效應(yīng)已經(jīng)對(duì)浮游動(dòng)物產(chǎn)生了重要影響。這也是海冰(無冰天數(shù))在本研究和以往研究(Swadling, 2010)在對(duì)浮游動(dòng)物群落劃分解釋上存在差異的原因之一。盡管海冰對(duì)浮游動(dòng)物影響的信號(hào)在本研究中沒有被捕捉到, 餌料仍然對(duì)夏季初期宇航員海浮游動(dòng)物的種群補(bǔ)充有一定的解釋: 浮游動(dòng)物平均豐度最高的群落2區(qū)域具有最高的葉綠素濃度(表1, 2), 同時(shí)橈足類在葉綠素濃度相對(duì)較高的群落1、2、3所在區(qū)域呈現(xiàn)相對(duì)更“年輕”的種群結(jié)構(gòu)(圖6)。此外Bio-env分析結(jié)果也表明溫度、鹽度與葉綠素三者的結(jié)合更好地解釋了浮游動(dòng)物群落的劃分(表3)。這些都潛在地說明除了環(huán)流等大尺度的物理學(xué)過程, 餌料供應(yīng)也對(duì)宇航員海域浮游動(dòng)物群落結(jié)構(gòu)產(chǎn)生了重要的影響。

    4 結(jié)語

    本研究顯示宇航員海域夏季浮游動(dòng)物群落主要由橈足類、磷蝦、毛顎類、浮游被囊類等類群組成, 其中橈足類豐度對(duì)浮游動(dòng)物總豐度的貢獻(xiàn)可達(dá)90%以上。普里茲灣渦流、威德爾渦流、南極繞極流南界等物理海洋學(xué)過程是宇航員海浮游動(dòng)物群落分布模式的主要控制因素, 另外, 餌料(葉綠素)也對(duì)該海區(qū)浮游動(dòng)物的種群和群落結(jié)構(gòu)具有重要影響。

    本研究是繼10余年前2006年夏季“BROKE-West”航次以來全球唯一的一次對(duì)宇航員海較為系統(tǒng)的調(diào)查, 我們對(duì)此海域浮游動(dòng)物群落劃分及其與環(huán)境影響因子的關(guān)系進(jìn)行了探討。研究結(jié)果對(duì)數(shù)據(jù)相對(duì)匱乏的宇航員海進(jìn)行了數(shù)據(jù)補(bǔ)充, 為該海域今后進(jìn)行連續(xù)觀測(cè)的站位設(shè)計(jì)提供了基礎(chǔ)參考, 同時(shí)為南大洋變化的大背景下更好地監(jiān)測(cè)中上層浮游生態(tài)系統(tǒng)的長期變化及其對(duì)氣候環(huán)境變化的響應(yīng)奠定基礎(chǔ)。

    致謝 感謝中國第36次南極科學(xué)考察雪龍2號(hào)船全體船員對(duì)現(xiàn)場(chǎng)樣品采集工作的幫助。

    Atkinson A, 1998. Life cycle strategies of epipelagic copepods in the Southern Ocean. Journal of Marine Systems, 15(1—4): 289—311

    Atkinson A, Hill S L, Barange M, 2014. Sardine cycles, krill declines, and locust plagues: revisiting wasp-waist food webs. Trends in Ecology & Evolution, 29(6): 309—316

    Atkinson A, Hill S L, Pakhomov E A, 2019. Krill () distribution contracts southward during rapid regional warming. Nature Climate Change, 9(2): 142—147

    Atkinson A, Siegel V, Pakhomov E, 2004. Long-term decline in krill stock and increase in salps within the Southern Ocean. Nature, 432(7013): 100—103

    Atkinson A, Siegel V, Pakhomov E A, 2008. Oceanic circumpolar habitats of Antarctic krill. Marine Ecology Progress, 362(1): 1—23

    Atkinson A, Ward P, Hunt B P V, 2012. An overview of Southern Ocean zooplankton data: Abundance, biomass, feeding and functional relationships. CCAMLR Science, 19: 171—218

    Atkinson A, Ward P, Williams R, 1992. Diel vertical migration and feeding of copepods at an oceanic site near South Georgia. Marine Biology, 113(4): 583—593

    Beaugrand G, Reid P C, Iba?ez F, 2002. Reorganization of North Atlantic marine copepod biodiversity and climate. Science, 296(5573): 1692—1694

    Boyd P W, Claustre H, Levy M, 2019. Multi-faceted particle pumps drive carbon sequestration in the ocean. Nature, 568(7752): 327—335

    Calbet A, 2008. The trophic roles of microzooplankton in marine systems. ICES Journal of Marine Science, 65(3): 325—331

    Chapman C C, Lea M A, Meyer A, 2020. Defining Southern Ocean fronts and their influence on biological and physical processes in a changing climate. Nature Climate Change, 10(3): 209—219

    Conover R J, Huntley M, 1991. Copepods in ice-covered seas—Distribution, adaptations to seasonally limited food, metabolism, growth patterns and life cycle strategies in polar seas. Journal of Marine Systems, 2(1/2): 1—41

    Constable A J, Melbourne-Thomas J, Corney S P, 2014. Climate change and Southern Ocean ecosystems I: how changes in physical habitats directly affect marine biota. Global Change Biology, 20(10): 3004—3025

    Convey P, Peck L S, 2019. Antarctic environmental change and biological responses. Science Advances, 5(11): eaaz0888

    Donnelly J, Torres J J, Hopkins T L, 1994. Chemical composition of antarctic zooplankton during austral fall and winter. Polar Biology, 14(3): 171—183

    Dufrêne M, Legendre P, 1997. Species assemblages and indicator species: The need for a flexible asymmetrical approach. Ecological Monographs, 67(3): 345—366

    Estrada R, Harvey M, Gosselin M, 2012. Late-summer zooplankton community structure, abundance, and distribution in the Hudson Bay system (Canada) and their relationships with environmental conditions, 2003-2006. Progress in Oceanography, 101(1): 121—145

    Field J G, Clarke K R, Warwick R M, 1982. A Practical Strategy for Analysing Multispecies Distribution Patterns. Marine Ecology Progress Series, 8: 37—52

    Hagen W, Schnack-Schiel S B, 1996. Seasonal lipid dynamics in dominant Antarctic copepods: Energy for overwintering or reproduction?. Deep Sea Research Part I Oceanographic Research Papers, 43(2): 139—158

    Henley S F, Schofield O M, Hendry K R, 2019. Variability and change in the west Antarctic Peninsula marine system: Research priorities and opportunities. Progress in Oceanography, 173: 208—237

    Hosie G W, Cochran T G, 1994. Mesoscale distribution patterns of macrozooplankton communities in Prydz Bay, Antarctica - January to February 1991. Marine Ecology Progress Series, 106: 21—39

    Hunt B P V, Pakhomov E A, Trotsenko B G, 2007. The macrozooplankton of the Cosmonaut Sea, east Antarctica (30°E–60°E), 1987—1990. Deep Sea Research Part I: Oceanographic Research Papers, 54(7): 1042—1069

    Jónasdóttir S H, Visser A W, Richardson K, 2015. Seasonal copepod lipid pump promotes carbon sequestration in the deep North Atlantic. Proceedings of the National Academy of Sciences of the United States of America, 112(39): 12122—12126

    Kawaguchi S, Nicol S, Virtue P, 2010. Krill demography and large-scale distribution in the Western Indian Ocean sector of the Southern Ocean (CCAMLR Division 58.4.2) in Austral summer of 2006. Deep Sea Research Part II: Topical Studies in Oceanography, 57(9/10): 934—947

    Lenaerts J T M, Lhermitte S, Drews R, 2017. Meltwater produced by wind–albedo interaction stored in an East Antarctic ice shelf. Nature Climate Change, 7(1): 58—62

    Makabe R, Tanimura A, Fukuchi M, 2014. Comparison of mesh size effects on mesozooplankton collection efficiency in the Southern Ocean.Journal of Plankton Research, 34( 5) : 432—436

    Mayzaud P, Pakhomov E A, 2014. The role of zooplankton communities in carbon recycling in the ocean: the case of the Southern Ocean. Journal of Plankton Research, 36(6): 1543—1556

    Meijers A J S, Klocker A, Bindoff N L, 2010. The circulation and water masses of the Antarctic shelf and continental slope between 30 and 80°E. Deep Sea Research Part II: Topical Studies in Oceanography, 57(9/10): 723—737

    Montes-Hugo M, Doney S C, Ducklow H W, 2009. Recent changes in phytoplankton communities associated with rapid regional climate change along the western Antarctic Peninsula. Science, 323(5920): 1470—1473

    Moriarty R, O'Brien T D, 2013. Distribution of mesozooplankton biomass in the global ocean. Earth System Science Data, 5(1): 45—55

    Nicol S, Pauly T, Bindoff N L, 2000. Ocean circulation off East Antarctica affects ecosystem structure and sea-ice extent. Nature, 406(6795): 504—507

    Orsi A H, Whitworth III T, Nowlin Jr W D, 1995. On the meridional extent and fronts of the Antarctic Circumpolar Current. Deep Sea Research Part I: Oceanographic Research Papers, 42(5): 641—673

    Pakhomov E A, Pshenichnov L K, Krotet A, 2020. Zooplankton distribution and community structure in the Pacific and Atlantic sectors of the Southern Ocean during austral summer 2017—18: a pilot study conducted from Ukrainian long-liners. Journal of Marine Science and Engineering, 7(8): 488

    Pinkerton M H, Décima M, Kitchener J A, 2020. Zooplankton in the Southern Ocean from the continuous plankton recorder: Distributions and long-term change. Deep Sea Research Part I: Oceanographic Research Papers, 162: 103303

    Rintoul S R, Chown S L, Deconto R M, 2018. Choosing the future of Antarctica. Nature, 558(7709): 233—241

    Saba G K, Fraser W R, Saba V S, 2014. Winter and spring controls on the summer food web of the coastal West Antarctic Peninsula. Nature Communications, 5(1): 4318

    Schmidt K, Schlosser C, Atkinson A, 2016. Zooplankton gut passage mobilizes lithogenic iron for ocean productivity. Current Biology, 26(19): 2667—2673

    Schnack-Schiel S B, Hagen W, 1994. Life cycle strategies and seasonal variations in distribution and population structure of four dominant calanoid copepod species in the eastern Weddell Sea, Antarctica. Journal of Plankton Research, 16(11): 1543—1566

    Steinberg D K, Landry M R, 2017. Zooplankton and the ocean carbon cycle. Annual Review of Marine Science, (1): 413—444

    Swadling K M, Kawaguchi S, Hosie G W, 2010. Antarctic mesozooplankton community structure during BROKE-West (30°E–80°E), January–February 2006. Deep Sea Research Part II: Topical Studies in Oceanography, 57(9/10): 887—904

    Voronina N M, 1998. Comparative abundance and distribution of major filter-feeders in the Antarctic pelagic zone. Journal of Marine Systems, 17(1/4): 375—390

    Williams G D, Nicol S, Aoki S, 2010. Surface oceanography of BROKE-West, along the Antarctic margin of the south-west Indian Ocean (30–80°E). Deep Sea Research Part II: Topical Studies in Oceanography, 57(9/10): 738—757

    Yang G, Atkinson A, Hill L S, 2021. Changing circumpolar distributions and isoscapes of Antarctic krill: Indo-Pacific habitat refuges counter long-term degradation of the Atlantic sector. Limnology and Oceanography, 66(1): 272—287

    Yang G, Li C L, Sun S, 2011. Inter-annual variation in summer zooplankton community structure in Prydz Bay, Antarctica, from 1999 to 2006. Polar Biology, 34(6): 921—932

    THE ZOOPLANKTON COMMUNITY IN COSMONAUT SEA: COMMUNITY STRUCTURE AND ENVIRONMENTAL FACTORS

    MOU Wen-Xiu1, 2, 3, YANG Guang1, 2, 4, HAO Qiang5, XU Zhi-Qiang1, LI Chao-Lun1, 2, 3, 4

    (1. Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 2. Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China; 3. University of Chinese Academy of Sciences, Beijing 100049, China; 4. Center for Ocean Mega Science, Chinese Academy of Sciences, Qingdao 266071, China; 5.Second Institute of Oceanography, Ministry of Natural Resources, Hangzhou 310012, China)

    Based on samples collected using a Norpac net (0.5 m2, 330 μm) in the Cosmonaut Sea, Antarctic during the 36th Chinese National Antarctic Research Expedition (December 2019 to January 2020), the zooplankton community structure, abundance, population structure of dominant species and their relationship with environmental factors were studied. Our results show that total abundance of zooplankton in Cosmonaut Sea in the summer of 2019/2020 varied from 8.93 to 102.00 ind./m3. Zooplankton was mainly composed of copepods, krill, chaetognaths, and tunicates. Zooplankton abundance in each station was dominated by copepods (>90%). The large-sized copepods, such as,,and, were more widely distributed in the survey area, while small-sized ones, such as,,, andwere more abundant. Five groups of zooplankton community could be recognized based on cluster analysis and they were closely linked to local physical oceanographic processes (the Prydz Bay Gyre, Weddell Gyre, and the Southern boundary of the Antarctic Circumpolar Current). Food availability (chlorophyll) also had an important impact on the population structure and community structure of zooplankton.

    community structure; zooplankton; Cosmonaut Sea; Southern Ocean

    Q958.8

    10.11693/hyhz20210200043

    * 國家自然科學(xué)基金項(xiàng)目, 41876217號(hào); 自然資源部南大洋專項(xiàng), IRASCC2020-2022-No. 01-02-02B, 02-04-02。牟文秀, 碩士研究生, E-mail: 17863804640@163.com

    李超倫, 博士生導(dǎo)師, 教授, E-mail: lcl@qdio.ac.cn

    2021-02-08,

    2021-03-17

    猜你喜歡
    南大洋橈足類磷蝦
    磷蝦真是“蝦無敵”
    南極磷蝦粉在水產(chǎn)飼料中的應(yīng)用
    湖南飼料(2021年4期)2021-10-13 07:32:46
    基于漂流浮標(biāo)的南大洋衛(wèi)星高度計(jì)有效波高研究
    “南大洋”來了,世界地圖要改?
    白洋淀枝角類和橈足類初步調(diào)查
    世界環(huán)境(2018年6期)2019-01-12 07:45:02
    “美味”的磷蝦
    淺析南大洋海流、海冰對(duì)磷蝦分布的影響
    “美味”的磷蝦
    海洋橈足類生態(tài)學(xué)研究進(jìn)展
    河南科技(2015年21期)2015-03-26 21:21:21
    久久人妻熟女aⅴ| 日本vs欧美在线观看视频| 精品人妻熟女毛片av久久网站| 精品少妇一区二区三区视频日本电影| 一本—道久久a久久精品蜜桃钙片| 丝袜人妻中文字幕| 亚洲精品日韩在线中文字幕| 成人免费观看视频高清| 在线观看免费视频网站a站| 我要看黄色一级片免费的| 这个男人来自地球电影免费观看| 91麻豆精品激情在线观看国产 | 亚洲av片天天在线观看| 欧美黄色淫秽网站| 大香蕉久久成人网| 久久鲁丝午夜福利片| 亚洲成人免费av在线播放| 欧美精品亚洲一区二区| 啦啦啦视频在线资源免费观看| 精品久久蜜臀av无| 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻人人爽人人添夜夜欢视频| 亚洲欧美色中文字幕在线| 国产日韩欧美视频二区| 国产精品熟女久久久久浪| 亚洲欧美精品综合一区二区三区| 无遮挡黄片免费观看| 欧美精品av麻豆av| 在线观看免费日韩欧美大片| 天天操日日干夜夜撸| 久久精品人人爽人人爽视色| 捣出白浆h1v1| 欧美老熟妇乱子伦牲交| 欧美精品高潮呻吟av久久| 自线自在国产av| 欧美 日韩 精品 国产| 午夜视频精品福利| 日本欧美国产在线视频| 一区二区三区四区激情视频| 国产麻豆69| 精品熟女少妇八av免费久了| 日韩av不卡免费在线播放| 丝袜美腿诱惑在线| 久久久久国产一级毛片高清牌| 一本综合久久免费| 国产深夜福利视频在线观看| 久久人妻熟女aⅴ| 极品人妻少妇av视频| 可以免费在线观看a视频的电影网站| 首页视频小说图片口味搜索 | 国产麻豆69| 男的添女的下面高潮视频| 久久女婷五月综合色啪小说| 午夜福利免费观看在线| 99re6热这里在线精品视频| 亚洲人成网站在线观看播放| 一区二区日韩欧美中文字幕| 久久国产精品大桥未久av| 亚洲av成人精品一二三区| 高潮久久久久久久久久久不卡| 男男h啪啪无遮挡| 一级毛片我不卡| 看免费av毛片| 欧美久久黑人一区二区| 久久午夜综合久久蜜桃| 日本wwww免费看| 超色免费av| 免费看十八禁软件| 国产精品av久久久久免费| 18禁国产床啪视频网站| tube8黄色片| 女人高潮潮喷娇喘18禁视频| 成年av动漫网址| 美女视频免费永久观看网站| 手机成人av网站| 国产精品国产三级国产专区5o| 伊人久久大香线蕉亚洲五| 日韩欧美一区视频在线观看| 成人午夜精彩视频在线观看| 国产一区有黄有色的免费视频| 十八禁高潮呻吟视频| 97在线人人人人妻| 亚洲国产成人一精品久久久| 老司机靠b影院| 国产野战对白在线观看| av视频免费观看在线观看| 国产真人三级小视频在线观看| 在线观看一区二区三区激情| 在线亚洲精品国产二区图片欧美| 建设人人有责人人尽责人人享有的| 欧美精品亚洲一区二区| 午夜激情av网站| 一二三四在线观看免费中文在| 精品一区二区三区av网在线观看 | 午夜免费鲁丝| 午夜精品国产一区二区电影| 亚洲精品美女久久av网站| 99精国产麻豆久久婷婷| 国产高清不卡午夜福利| 人妻一区二区av| 丝瓜视频免费看黄片| 精品人妻1区二区| 亚洲欧美日韩另类电影网站| 亚洲专区国产一区二区| 欧美人与性动交α欧美软件| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品人妻久久久影院| 日韩伦理黄色片| 波多野结衣一区麻豆| 欧美日韩一级在线毛片| 成人黄色视频免费在线看| 亚洲美女黄色视频免费看| 汤姆久久久久久久影院中文字幕| 大香蕉久久成人网| 国产97色在线日韩免费| 啦啦啦在线观看免费高清www| 一本—道久久a久久精品蜜桃钙片| 一区二区av电影网| 日本猛色少妇xxxxx猛交久久| 婷婷丁香在线五月| 777米奇影视久久| √禁漫天堂资源中文www| 50天的宝宝边吃奶边哭怎么回事| 丝袜美足系列| 亚洲国产欧美日韩在线播放| 久久久久久人人人人人| 18禁观看日本| 中文字幕制服av| 欧美av亚洲av综合av国产av| 久久精品aⅴ一区二区三区四区| 看免费成人av毛片| 夫妻性生交免费视频一级片| 少妇裸体淫交视频免费看高清 | 亚洲国产精品999| 男女无遮挡免费网站观看| avwww免费| 久久av网站| 日本黄色日本黄色录像| 天堂中文最新版在线下载| 欧美成狂野欧美在线观看| 各种免费的搞黄视频| 两人在一起打扑克的视频| 一级毛片 在线播放| 国产人伦9x9x在线观看| 激情五月婷婷亚洲| 日本av手机在线免费观看| 中文字幕av电影在线播放| 丰满迷人的少妇在线观看| 黄色毛片三级朝国网站| 丁香六月天网| 欧美亚洲 丝袜 人妻 在线| 亚洲,一卡二卡三卡| 极品少妇高潮喷水抽搐| 又大又爽又粗| 国产免费一区二区三区四区乱码| 80岁老熟妇乱子伦牲交| 青草久久国产| 91精品国产国语对白视频| 黄网站色视频无遮挡免费观看| 巨乳人妻的诱惑在线观看| 18禁国产床啪视频网站| 亚洲精品美女久久久久99蜜臀 | 一级片'在线观看视频| 久久精品久久精品一区二区三区| 国产片特级美女逼逼视频| 女人高潮潮喷娇喘18禁视频| 男的添女的下面高潮视频| 黑人欧美特级aaaaaa片| 精品亚洲乱码少妇综合久久| 1024视频免费在线观看| 精品国产国语对白av| 男女无遮挡免费网站观看| 超色免费av| 国产成人a∨麻豆精品| videos熟女内射| 欧美激情高清一区二区三区| 日本av免费视频播放| 国产亚洲欧美精品永久| 桃花免费在线播放| xxx大片免费视频| 久久青草综合色| 亚洲欧美清纯卡通| 亚洲精品一卡2卡三卡4卡5卡 | 老司机午夜十八禁免费视频| 国产精品一区二区精品视频观看| 视频区欧美日本亚洲| 可以免费在线观看a视频的电影网站| 国产精品久久久久久精品电影小说| 新久久久久国产一级毛片| kizo精华| 国产免费现黄频在线看| 国产成人精品久久二区二区91| 成年人黄色毛片网站| 成人手机av| 又大又黄又爽视频免费| 日韩制服骚丝袜av| 成人三级做爰电影| 国产精品久久久久成人av| 波多野结衣av一区二区av| 日韩一本色道免费dvd| 久久久久久久精品精品| avwww免费| 波多野结衣一区麻豆| 精品久久久精品久久久| 男女床上黄色一级片免费看| 青草久久国产| 久久国产精品人妻蜜桃| 亚洲精品国产av成人精品| 欧美 亚洲 国产 日韩一| 丝袜在线中文字幕| 巨乳人妻的诱惑在线观看| 90打野战视频偷拍视频| 亚洲人成电影观看| 日韩中文字幕视频在线看片| 精品高清国产在线一区| 欧美日韩黄片免| 国产深夜福利视频在线观看| 亚洲国产看品久久| 亚洲国产日韩一区二区| 亚洲国产精品一区二区三区在线| 久久精品久久久久久噜噜老黄| 爱豆传媒免费全集在线观看| 黄色a级毛片大全视频| a级毛片黄视频| 精品人妻1区二区| www.av在线官网国产| 午夜福利在线免费观看网站| 婷婷色麻豆天堂久久| 国产男女内射视频| 欧美日韩视频高清一区二区三区二| 欧美性长视频在线观看| 久久av网站| 色婷婷av一区二区三区视频| 亚洲国产欧美网| 999精品在线视频| 日韩电影二区| 赤兔流量卡办理| 亚洲国产日韩一区二区| 久久精品国产亚洲av高清一级| 99久久精品国产亚洲精品| 久久国产亚洲av麻豆专区| 蜜桃国产av成人99| 国产淫语在线视频| 亚洲七黄色美女视频| 七月丁香在线播放| 一区二区三区四区激情视频| 丁香六月天网| 又大又黄又爽视频免费| 国产亚洲一区二区精品| 国产精品 国内视频| 大片免费播放器 马上看| 国产真人三级小视频在线观看| 久久久久久久久久久久大奶| 操美女的视频在线观看| 2021少妇久久久久久久久久久| 亚洲精品在线美女| 一区福利在线观看| 亚洲欧美一区二区三区黑人| 在线观看免费午夜福利视频| 亚洲国产精品一区三区| 国产精品久久久久成人av| tube8黄色片| 国产高清不卡午夜福利| 两人在一起打扑克的视频| 女性被躁到高潮视频| 久久精品熟女亚洲av麻豆精品| 亚洲一区二区三区欧美精品| 一级片'在线观看视频| 国产99久久九九免费精品| av不卡在线播放| 国产成人欧美| 老汉色av国产亚洲站长工具| 日韩av不卡免费在线播放| 国产一区二区 视频在线| 亚洲欧美日韩另类电影网站| 在线观看一区二区三区激情| 亚洲人成77777在线视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人手机| 自线自在国产av| 欧美黑人欧美精品刺激| 成人亚洲精品一区在线观看| 99久久精品国产亚洲精品| 亚洲,一卡二卡三卡| 十八禁高潮呻吟视频| 亚洲人成77777在线视频| 精品第一国产精品| 久久天堂一区二区三区四区| 国产精品欧美亚洲77777| videos熟女内射| 精品熟女少妇八av免费久了| 日本欧美国产在线视频| a级片在线免费高清观看视频| 久久精品亚洲av国产电影网| 欧美精品av麻豆av| 欧美日韩av久久| 女人高潮潮喷娇喘18禁视频| 一二三四社区在线视频社区8| 亚洲国产精品一区二区三区在线| 婷婷色综合大香蕉| 在线观看免费高清a一片| 久久精品国产a三级三级三级| 捣出白浆h1v1| 国产极品粉嫩免费观看在线| 又黄又粗又硬又大视频| 男女无遮挡免费网站观看| 日韩一本色道免费dvd| 亚洲精品国产一区二区精华液| 亚洲精品av麻豆狂野| 欧美精品亚洲一区二区| 亚洲国产精品999| 国产一区二区三区av在线| 黄色怎么调成土黄色| 国产爽快片一区二区三区| 一区二区三区精品91| 欧美日韩精品网址| 大型av网站在线播放| 久久精品国产综合久久久| 日本黄色日本黄色录像| 婷婷丁香在线五月| 黄网站色视频无遮挡免费观看| 18禁国产床啪视频网站| 如日韩欧美国产精品一区二区三区| 国产高清国产精品国产三级| 亚洲欧美日韩高清在线视频 | 久久精品国产亚洲av涩爱| 亚洲免费av在线视频| 久久久久久人人人人人| 别揉我奶头~嗯~啊~动态视频 | kizo精华| av天堂久久9| 欧美黑人精品巨大| 午夜激情久久久久久久| 欧美97在线视频| 成人国产av品久久久| 亚洲欧美成人综合另类久久久| 国产成人a∨麻豆精品| 国产成人系列免费观看| 亚洲,欧美精品.| 人人妻人人澡人人看| 看免费av毛片| 国产深夜福利视频在线观看| 久久精品成人免费网站| 欧美国产精品一级二级三级| 亚洲五月色婷婷综合| 国产欧美日韩精品亚洲av| 黄色片一级片一级黄色片| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品久久成人aⅴ小说| 后天国语完整版免费观看| 欧美乱码精品一区二区三区| 亚洲第一av免费看| 亚洲人成电影观看| 久久影院123| 亚洲专区中文字幕在线| 国产一区二区激情短视频 | 国产精品国产三级专区第一集| 日日夜夜操网爽| 青春草视频在线免费观看| 如日韩欧美国产精品一区二区三区| 免费在线观看日本一区| 十八禁高潮呻吟视频| 少妇裸体淫交视频免费看高清 | 九草在线视频观看| 黄频高清免费视频| 人人妻人人澡人人看| 日韩一区二区三区影片| 大型av网站在线播放| 99久久综合免费| 亚洲精品一二三| 国产国语露脸激情在线看| 欧美日韩精品网址| 日本黄色日本黄色录像| 最黄视频免费看| 午夜免费鲁丝| videos熟女内射| 日韩人妻精品一区2区三区| 中文字幕高清在线视频| 国产精品麻豆人妻色哟哟久久| 亚洲国产欧美一区二区综合| 人人妻人人添人人爽欧美一区卜| 亚洲精品中文字幕在线视频| 久久女婷五月综合色啪小说| 欧美中文综合在线视频| 男人舔女人的私密视频| 性色av乱码一区二区三区2| 日韩av免费高清视频| 亚洲专区国产一区二区| 欧美老熟妇乱子伦牲交| 大香蕉久久成人网| 欧美日韩视频高清一区二区三区二| 波多野结衣av一区二区av| 丝袜人妻中文字幕| 天堂8中文在线网| 国产精品香港三级国产av潘金莲 | 日日摸夜夜添夜夜爱| 欧美在线黄色| 日本一区二区免费在线视频| 老司机影院成人| 大香蕉久久网| 欧美成狂野欧美在线观看| 久久精品久久精品一区二区三区| 天天影视国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 高清欧美精品videossex| 国产成人一区二区在线| 国产精品秋霞免费鲁丝片| 老熟女久久久| xxx大片免费视频| 青春草视频在线免费观看| 男女边吃奶边做爰视频| 亚洲,欧美精品.| 色视频在线一区二区三区| 看免费成人av毛片| 欧美精品av麻豆av| 国产免费又黄又爽又色| 久久精品久久久久久噜噜老黄| 亚洲欧美精品综合一区二区三区| 五月开心婷婷网| 男女边吃奶边做爰视频| 乱人伦中国视频| 亚洲av日韩在线播放| 久久热在线av| 人妻 亚洲 视频| 美女主播在线视频| 好男人视频免费观看在线| 考比视频在线观看| 亚洲伊人久久精品综合| 亚洲欧美激情在线| 天天添夜夜摸| 国产高清国产精品国产三级| 日韩视频在线欧美| 国产成人啪精品午夜网站| 国产午夜精品一二区理论片| 久久99热这里只频精品6学生| 国产一区二区三区综合在线观看| 精品国产国语对白av| 男女之事视频高清在线观看 | 久久亚洲国产成人精品v| 色综合欧美亚洲国产小说| 亚洲情色 制服丝袜| 在线观看免费视频网站a站| 欧美大码av| 国产一区二区三区综合在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲一码二码三码区别大吗| 亚洲熟女毛片儿| 高清黄色对白视频在线免费看| 亚洲国产最新在线播放| 人人妻人人爽人人添夜夜欢视频| 亚洲欧洲日产国产| 热re99久久精品国产66热6| 大话2 男鬼变身卡| 中文字幕高清在线视频| 1024香蕉在线观看| 两人在一起打扑克的视频| 男人舔女人的私密视频| 日韩一本色道免费dvd| 久久 成人 亚洲| 91精品三级在线观看| 高清av免费在线| 亚洲国产精品国产精品| 制服诱惑二区| 亚洲成人免费电影在线观看 | 欧美日韩亚洲高清精品| 国产成人欧美在线观看 | 自拍欧美九色日韩亚洲蝌蚪91| 国产爽快片一区二区三区| 亚洲精品日韩在线中文字幕| 国产片特级美女逼逼视频| av有码第一页| 国产精品一区二区在线不卡| 午夜免费成人在线视频| 日本av免费视频播放| 亚洲欧美精品综合一区二区三区| 欧美黑人精品巨大| 首页视频小说图片口味搜索 | 大型av网站在线播放| 咕卡用的链子| 捣出白浆h1v1| 欧美大码av| 国产又爽黄色视频| 两人在一起打扑克的视频| 韩国精品一区二区三区| 免费在线观看黄色视频的| 午夜精品国产一区二区电影| 天堂中文最新版在线下载| 黄色毛片三级朝国网站| 成人免费观看视频高清| 国产精品 欧美亚洲| 少妇粗大呻吟视频| 免费女性裸体啪啪无遮挡网站| 国产1区2区3区精品| 欧美激情 高清一区二区三区| 亚洲精品乱久久久久久| 国产爽快片一区二区三区| 欧美日韩一级在线毛片| 国产亚洲一区二区精品| a级毛片在线看网站| 亚洲午夜精品一区,二区,三区| 国产有黄有色有爽视频| 久久ye,这里只有精品| 久热这里只有精品99| 国产成人欧美在线观看 | 日韩大片免费观看网站| 欧美激情 高清一区二区三区| 午夜日韩欧美国产| 人成视频在线观看免费观看| 午夜两性在线视频| 国产精品偷伦视频观看了| 亚洲精品一二三| 亚洲av国产av综合av卡| 亚洲五月色婷婷综合| 久久久久久久久久久久大奶| 黄色毛片三级朝国网站| 亚洲欧美色中文字幕在线| 少妇人妻 视频| 久久久久久久国产电影| 1024视频免费在线观看| 一级黄色大片毛片| 亚洲激情五月婷婷啪啪| 一本综合久久免费| 国产av精品麻豆| 夜夜骑夜夜射夜夜干| 少妇裸体淫交视频免费看高清 | 老熟女久久久| 久久这里只有精品19| 国产亚洲精品第一综合不卡| 两人在一起打扑克的视频| 亚洲伊人色综图| 亚洲人成电影免费在线| 欧美 日韩 精品 国产| 欧美日韩一级在线毛片| 好男人电影高清在线观看| 一区二区三区精品91| 侵犯人妻中文字幕一二三四区| 男的添女的下面高潮视频| 99国产精品一区二区三区| 成人免费观看视频高清| 国产在视频线精品| 免费观看av网站的网址| 美女国产高潮福利片在线看| 日韩人妻精品一区2区三区| 中文乱码字字幕精品一区二区三区| 精品国产乱码久久久久久男人| av电影中文网址| 久久99一区二区三区| 一级a爱视频在线免费观看| 只有这里有精品99| 少妇粗大呻吟视频| 亚洲国产中文字幕在线视频| av视频免费观看在线观看| 久久午夜综合久久蜜桃| 搡老岳熟女国产| 高清av免费在线| videos熟女内射| 成年人黄色毛片网站| 侵犯人妻中文字幕一二三四区| 交换朋友夫妻互换小说| 91九色精品人成在线观看| 国产熟女欧美一区二区| 亚洲国产精品成人久久小说| 亚洲av国产av综合av卡| 亚洲专区国产一区二区| 在线观看免费日韩欧美大片| 韩国高清视频一区二区三区| 久久精品亚洲av国产电影网| 国产视频一区二区在线看| 日本色播在线视频| 国产爽快片一区二区三区| 久9热在线精品视频| 欧美 日韩 精品 国产| 波多野结衣一区麻豆| 亚洲伊人色综图| 又大又黄又爽视频免费| 精品国产一区二区三区四区第35| 黄色毛片三级朝国网站| 久久精品成人免费网站| 伊人久久大香线蕉亚洲五| 日韩大片免费观看网站| h视频一区二区三区| 9色porny在线观看| h视频一区二区三区| 久久午夜综合久久蜜桃| 精品一品国产午夜福利视频| 亚洲精品国产色婷婷电影| 91麻豆精品激情在线观看国产 | 国产日韩欧美亚洲二区| 十八禁人妻一区二区| 亚洲专区国产一区二区| 波多野结衣一区麻豆| 黄色 视频免费看| 操美女的视频在线观看| 国产淫语在线视频| 啦啦啦在线观看免费高清www| 久久久精品免费免费高清| 久久人人爽av亚洲精品天堂| 国产精品.久久久| 又黄又粗又硬又大视频| 亚洲国产欧美一区二区综合| 热re99久久精品国产66热6| 成人黄色视频免费在线看| 自线自在国产av| 免费不卡黄色视频| 黑人巨大精品欧美一区二区蜜桃| 国产成人91sexporn| 乱人伦中国视频| 亚洲成国产人片在线观看| 午夜福利乱码中文字幕| 久久ye,这里只有精品| 亚洲国产精品999| av在线播放精品| 99国产综合亚洲精品| av不卡在线播放| 母亲3免费完整高清在线观看|