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

    近岸典型生態(tài)系統(tǒng)大型底棲動(dòng)物的次級(jí)生產(chǎn)力
    ——以甌江口、洞頭列島為例

    2017-11-23 02:33:38賈勝華廖一波曾江寧高愛根湯雁濱
    生態(tài)學(xué)報(bào) 2017年21期
    關(guān)鍵詞:列島洞頭甌江

    賈勝華,廖一波,曾江寧,*,壽 鹿,黃 偉,高愛根,湯雁濱

    1 國(guó)家海洋局海洋生態(tài)系統(tǒng)與生物地球化學(xué)重點(diǎn)實(shí)驗(yàn)室, 杭州 310012 2 國(guó)家海洋局第二海洋研究所, 杭州 310012

    近岸典型生態(tài)系統(tǒng)大型底棲動(dòng)物的次級(jí)生產(chǎn)力
    ——以甌江口、洞頭列島為例

    賈勝華1,2,廖一波1,2,曾江寧1,2,*,壽 鹿1,2,黃 偉1,2,高愛根1,2,湯雁濱1,2

    1 國(guó)家海洋局海洋生態(tài)系統(tǒng)與生物地球化學(xué)重點(diǎn)實(shí)驗(yàn)室, 杭州 310012 2 國(guó)家海洋局第二海洋研究所, 杭州 310012

    為更好了解和估算河口及島嶼生態(tài)系統(tǒng)大型底棲動(dòng)物次級(jí)生產(chǎn)力的情況,以近岸兩個(gè)典型生態(tài)系統(tǒng)河口-甌江口、島嶼-洞頭列島為例,于2015年4月和10月在該海域進(jìn)行大型底棲動(dòng)物調(diào)查,采用3種Brey經(jīng)驗(yàn)?zāi)P蛯?duì)大型底棲動(dòng)物次級(jí)生產(chǎn)力進(jìn)行估算,并分析次級(jí)生產(chǎn)力與環(huán)境因子的關(guān)系。結(jié)果表明:(1)以去灰干重(ash-free dry weight, AFDW)計(jì),Brey經(jīng)驗(yàn)公式估算甌江口和洞頭列島海域年均次級(jí)生產(chǎn)力分別為0.48 g(AFDW) m-2a-1和0.70 g(AFDW) m-2a-1,Brey模型估算甌江口和洞頭列島海域年均次級(jí)生產(chǎn)力分別為0.51 g(AFDW) m-2a-1和1.55 g(AFDW) m-2a-1,Brey模型估算甌江口和洞頭列島海域年均次級(jí)生產(chǎn)力分別為0.25 g(AFDW) m-2a-1和0.99 g(AFDW) m-2a-1。(2)3種經(jīng)驗(yàn)?zāi)P头从炒笮偷讞珓?dòng)物年均次級(jí)生產(chǎn)力空間分布趨勢(shì)一致:甌江口海域有1高值區(qū),位于甌江口南部;洞頭列島海域有2高值區(qū),分別位于洞頭本島東北部的三盤島與花崗島之間和東黃島東南部的島嶼之間。(3)甌江口海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力主要貢獻(xiàn)種為雙鰓內(nèi)卷齒蠶Aglaophamusdibranchis、長(zhǎng)吻沙蠶Glycerachirori、異蚓蟲Heteromastusfiliforms、焦河籃蛤Potamocorbulaustulata、薄云母蛤Yoldiasimilis、圓筒原盒螺Eocylichnabraunsi,它們對(duì)次級(jí)生產(chǎn)力的貢獻(xiàn)率超54.2%;洞頭列島海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力主要貢獻(xiàn)種為雙鰓內(nèi)卷齒蠶、異蚓蟲、薄云母蛤、絨毛細(xì)足蟹Raphidopusciliatus、隆線強(qiáng)蟹Eucratecrenata、腦紐蟲Cerebratulinasp.、紅狼牙鰕虎魚Odontamblyopusrubicundus,它們對(duì)次級(jí)生產(chǎn)力的貢獻(xiàn)率超57.1%。(4)次級(jí)生產(chǎn)力與環(huán)境因子的相關(guān)性分析顯示化學(xué)需氧量、懸浮物和表層沉積物中值粒徑是影響甌江口海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力的重要環(huán)境因子,而環(huán)境因子與洞頭列島海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力相關(guān)關(guān)系不顯著。(5)甌江口海域Brey經(jīng)驗(yàn)公式與Brey模型估算結(jié)果基本一致,洞頭列島海域Brey經(jīng)驗(yàn)公式與Brey模型估算結(jié)果基本一致。

    次級(jí)生產(chǎn)力;大型底棲動(dòng)物;洞頭列島;甌江口;Brey經(jīng)驗(yàn)?zāi)P?/p>

    大型底棲動(dòng)物次級(jí)生產(chǎn)力是生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動(dòng)量化的基礎(chǔ),是表征生態(tài)系統(tǒng)功能的一個(gè)重要參數(shù)[1- 2]。估算次級(jí)生產(chǎn)力不僅可幫助人們了解大型底棲動(dòng)物的功能變化,而且對(duì)研究生態(tài)系統(tǒng)的物流、能流過程和海洋生物資源的持續(xù)利用以及評(píng)價(jià)海洋環(huán)境壓力也具有重要意義[3- 5]。河口與島嶼是近岸海域兩個(gè)典型的生態(tài)系統(tǒng),河口區(qū)位于海洋與陸地交界處,受入海河流攜帶的陸源沉積物、有機(jī)物以及無(wú)機(jī)營(yíng)養(yǎng)物質(zhì)影響,海洋環(huán)境具有高營(yíng)養(yǎng)化、水體層化、環(huán)流的特點(diǎn),歷來(lái)是研究的熱點(diǎn)區(qū)域[6];海島處于海岸帶的前沿,在發(fā)展沿海經(jīng)濟(jì)、推動(dòng)國(guó)防建設(shè)、保護(hù)海洋環(huán)境等方面具有重要意義。目前,國(guó)內(nèi)報(bào)道島嶼與河口大型底棲動(dòng)物次級(jí)生產(chǎn)力的研究較少[7],并且鮮有比較島嶼與河口兩個(gè)典型生態(tài)系統(tǒng)次級(jí)生產(chǎn)力的報(bào)道。而且國(guó)內(nèi)大型底棲動(dòng)物次級(jí)生產(chǎn)力估算基本采用Brey經(jīng)驗(yàn)公式[7- 12],估算結(jié)果比較粗略,僅有膠州灣潮間帶[13]和遼東灣西部海域[14]兩篇報(bào)道采用Brey估算模型[15],而采用Brey估算模型[1]計(jì)算大型底棲動(dòng)物次級(jí)生產(chǎn)力國(guó)內(nèi)尚無(wú)報(bào)道。

    國(guó)外學(xué)者提出許多大型底棲動(dòng)物次級(jí)生產(chǎn)力估算模型[3- 4,15- 18],并對(duì)不同模型進(jìn)行了比較[4,19],這些模型或是建模數(shù)據(jù)范圍較小,不一定適用于低緯度海域[18];或是采用生物壽命等不易確定的參數(shù),難以應(yīng)用到實(shí)踐[4,16]。目前,Brey提出的經(jīng)驗(yàn)?zāi)P驮趪?guó)內(nèi)外大型底棲動(dòng)物群落次級(jí)生產(chǎn)力估算中應(yīng)用較為廣泛[5,8,10,12- 14]。國(guó)外次級(jí)生產(chǎn)力估算模型已由Brey經(jīng)驗(yàn)公式[3]發(fā)展到Brey估算模型[15],再到最近的Brey估算模型[1]。Brey估算模型[1]考慮的參數(shù)多達(dá)20個(gè),模型更為精細(xì)、復(fù)雜,雖其計(jì)算表單設(shè)計(jì)有輸入層、隱藏層和結(jié)果層,但只需輸入數(shù)據(jù),結(jié)果則自動(dòng)顯示,因此操作非常方便。

    研究海域位于27°41′—28°13′N,120°52′—121°16′E的浙江南部,包括洞頭列島海域、甌江口海域,兼有島嶼、河口2個(gè)典型生態(tài)系統(tǒng)。目前關(guān)于該海域大型底棲動(dòng)物次級(jí)生產(chǎn)力研究尚無(wú)報(bào)道,并且沒有比較兩個(gè)區(qū)域次級(jí)生產(chǎn)力的研究,本文擬采用Brey經(jīng)驗(yàn)公式[3]、Brey估算模型[1]以及Brey估算模型[15]分別計(jì)算甌江口海域及洞頭列島海域大型底棲動(dòng)物次級(jí)生產(chǎn)力,比較兩個(gè)研究區(qū)域的次級(jí)生產(chǎn)力大小和3種估算方法產(chǎn)生的結(jié)果差距,并分析其與環(huán)境因子的關(guān)系,以期為河口、島嶼兩個(gè)典型生態(tài)系統(tǒng)大型底棲動(dòng)物功能研究、生態(tài)環(huán)境保護(hù)提供理論基礎(chǔ)。

    1 材料與方法

    1.1 站位分布與樣品采集處理

    大型底棲動(dòng)物采樣于2015年4月(春季)和10月(秋季)進(jìn)行,共布設(shè)35個(gè)調(diào)查站位(圖1)。其中,S8、S14—S17、S25—S30、S32—S34調(diào)查站位位于甌江口海域;S1—S7、S9—S13、S18—S24、S31、S35調(diào)查站位位于洞頭列島海域。采用0.1 m2的Van Veen型抓斗式底質(zhì)采樣器取樣,每站重復(fù)取樣2次,合并為一個(gè)樣品。使用孔徑為0.5 mm的網(wǎng)篩分選樣品,采用5%的福爾馬林溶液固定所獲生物樣品,帶回實(shí)驗(yàn)室進(jìn)行種類鑒定。樣品的采集、保存、計(jì)數(shù)、稱量和資料整理均按照《海洋調(diào)查規(guī)范》(GB/T 12763.6——2007)[20]進(jìn)行。

    圖1 洞頭列島及鄰近海域大型底棲動(dòng)物采樣站位Fig.1 Sampling stations of macrozoobenthos in Dongtou Islands and adjacent sea area

    1.2 數(shù)據(jù)處理

    根據(jù)經(jīng)驗(yàn)?zāi)P蛯?duì)數(shù)據(jù)的要求,需將生物量濕重?cái)?shù)據(jù)轉(zhuǎn)換為去灰干重生物量或能值。Brey經(jīng)驗(yàn)公式[3]中生物量濕重(WW)的9/50為去灰干重(AFDW)[8]。Brey估算模型[1]和Brey估算模型[15]中,生物量濕重、去灰干重、能值的轉(zhuǎn)換關(guān)系[13-15]見表1。

    1.3 次級(jí)生產(chǎn)力估算經(jīng)驗(yàn)?zāi)P?/p>

    1.3.1Brey經(jīng)驗(yàn)公式[3]

    計(jì)算所有物種的年均去灰干重生物量和年均個(gè)體去灰干重,按照公式(1)逐種計(jì)算大型底棲動(dòng)物年均次級(jí)生產(chǎn)力,將所有物種年均次級(jí)生產(chǎn)力相加得到該站位大型底棲動(dòng)物群落年均次級(jí)生產(chǎn)力。

    lgP=a+b1lgB+b2lgW

    (1)

    表1 大型底棲動(dòng)物濕重生物量、去灰干重和能值的轉(zhuǎn)換關(guān)系[13-15]

    大型底棲動(dòng)物不同類群的經(jīng)驗(yàn)公式系數(shù)a,b1和b2見表2。其中:P為大型底棲動(dòng)物年均次級(jí)生產(chǎn)力,單位為g(AFDW) m-2a-1;B為大型底棲動(dòng)物年均去灰干重,單位為g(AFDW)/m2;W為大型底棲動(dòng)物個(gè)體的平均去灰干重,單位為g(AFDW)/個(gè)。

    表2 Brey經(jīng)驗(yàn)公式中不同類群的3個(gè)系數(shù)值[3]

    a:系數(shù)a,Coefficient a;b1:系數(shù)b1,Coefficient b1;b2:系數(shù)b2,Coefficient b2

    1.3.2Brey估算模型[4- 5,15,22]

    與Brey經(jīng)驗(yàn)公式[3]相比,Brey估算模型[15]除考慮大型底棲動(dòng)物的密度、生物量等生物因素,還考慮到水深、水溫、生境等環(huán)境因素,見公式(2)。

    lgP=7.947+lgB-2.294lgM-2409.856/(T+273)+0.168/D+0.194SubT+

    0.180InEpi+0.277MoEpi+0.174Taxon1-0.188Taxon2+0.33Taxon3-0.062Habitat1+582.851lgM/(T+273)

    (2)

    其中,P為大型底棲動(dòng)物年均次級(jí)生產(chǎn)力(kJ m-2a-1);B為年平均生物量對(duì)應(yīng)的能值(kJ/m2);M為年平均個(gè)體體重能值(kJ);T為平均底層水溫(℃);D為平均水深(m);SubT:潮下帶值為1,潮間帶為0;InEpi:底內(nèi)生物值為1,底表生物為0;MoEpi:移動(dòng)生物值為1,固著生物為0;Taxon1:環(huán)節(jié)動(dòng)物或甲殼動(dòng)物值為1,其他類群為0;Taxon2,棘皮動(dòng)物值為1,其他類群為0;Taxon3:昆蟲值為1,其他類群為0;Habitat1:湖泊生境值為1,其他生境值為0。

    1.3.3Brey估算模型[1,23]

    在Brey估算模型[15]多元線性回歸估算模型基礎(chǔ)上,Brey估算模型[1]增加大型底棲動(dòng)物食性、移動(dòng)性等參數(shù)提出多參數(shù)人工神經(jīng)網(wǎng)絡(luò)(Artificial Neural Network, ANN)估算模型,可對(duì)大型底棲動(dòng)物次級(jí)生產(chǎn)力進(jìn)行更好的預(yù)測(cè)[23],見公式(3)、(4)、(5)。由于模型參數(shù)較多,各參數(shù)對(duì)應(yīng)系數(shù)值可參考2012-01版本計(jì)算表單[1]。模型計(jì)算表單具有3個(gè)層次:輸入層、隱藏層和結(jié)果層,輸入層包含平均個(gè)體體重能值、水深、水溫、不同類群、移動(dòng)性、食性、生境等20個(gè)參數(shù),結(jié)果層為P/B的對(duì)數(shù)值。只需輸入每個(gè)物種的參數(shù),結(jié)果自動(dòng)顯示,操作比較簡(jiǎn)便。計(jì)算時(shí)先把20個(gè)待輸入數(shù)據(jù)放入其他表格中,再拷貝到計(jì)算表單中。

    log(P/B)=a0+a1H1+a2H2

    (3)

    H1=TanH(b0+b1logM+b2/T+b3logD+b4Mollusca+b5Annelida+b6Crustacea+b7Echinodermata+b8Insecta+b9Infauna+b10Sessile+b11Crawler+b12FacSwim+b13Herbivor+b14Omnivor+b15Carnivor+b16Lake+b17River+b18Marine+b19Subtidal+b20Exploited)

    (4)

    H2=TanH(c0+c1logM+c2/T+c3logD+c4Mollusca+c5Annelida+c6Crustacea+c7Echinodermata+c8Insecta+c9Infauna+c10Sessile+c11Crawler+c12FacSwim+c13Herbivor+c14Omnivor+c15Carnivor+c16Lake+c17River+c18Marine+c19Subtidal+c20Exploited)

    (5)

    其中,P為大型底棲動(dòng)物年均次級(jí)生產(chǎn)力(kJ m-2a-1);B為年平均生物量對(duì)應(yīng)的能值(kJ/m2);M為年平均個(gè)體體重能值(J);T為平均底層水溫(℃);D為平均水深(m);Mollusca:軟體動(dòng)物值為1,其他為0;Annelida:環(huán)節(jié)動(dòng)物值為1,其他為0;Crustacea:甲殼動(dòng)物值為1,其他為0;Echinodermata:棘皮動(dòng)物值為1,其他為0;Insecta:昆蟲值為1,其他為0;Infauna:底內(nèi)生物值為1,其他為0;Sessile:固著生物值為1,其他為0;Crawler:爬行動(dòng)物值為1,其他為0;FacSwim:兼性游泳生物值為1,其他為0;Herbivor:草食動(dòng)物動(dòng)物值為1,其他為0;Omnivor:雜食動(dòng)物值為1,其他為0;Carnivor:肉食動(dòng)物值為1,其他為0;Lake:湖泊生境值為1,其他生境值為0;River:河流生境值為1,其他生境值為0;Marine:海洋生境值為1,其他生境值為0;Subtidal:潮下帶值為1,其他生境值為0;Exploited指魚類,通常設(shè)置為0。

    本文將春、秋兩季大型底棲動(dòng)物棲息密度和生物量取平均值,然后分別采用3種經(jīng)驗(yàn)?zāi)P椭鸱N計(jì)算大型底棲動(dòng)物年均次級(jí)生產(chǎn)力。Brey經(jīng)驗(yàn)公式[3]運(yùn)用Excel軟件進(jìn)行數(shù)據(jù)處理,Brey估算模型[1]和Brey估算模型[15]從網(wǎng)站(http://www.thomas-brey.de/science/virtualhandbook/)下載計(jì)算表單,按要求進(jìn)行數(shù)據(jù)處理。為便于經(jīng)驗(yàn)?zāi)P陀?jì)算結(jié)果間比較,將Brey估算模型[1]和Brey估算模型[15]計(jì)算結(jié)果按照表1轉(zhuǎn)換關(guān)系換算為去灰干重(g(AFDW) m-2a-1)。采用SPSS 16.0軟件進(jìn)行相關(guān)性分析。

    2 結(jié)果

    2.1 種類組成

    甌江口海域兩航次共采集到大型底棲動(dòng)物81種。其中,多毛類39種,占48.1%;甲殼動(dòng)物23種,占28.4%;軟體動(dòng)物9種,占11.1%;棘皮動(dòng)物4種,占4.9%;其他類6種(包括魚類、螠蟲動(dòng)物、紐形動(dòng)物、海葵),占7.4%。洞頭列島海域兩航次共采集到大型底棲動(dòng)物111種。其中,多毛類50種,占45.0%;甲殼動(dòng)物30種,占27.0%;軟體動(dòng)物18種,占16.2%;棘皮動(dòng)物5種,占4.5%;其他類8種(包括魚類、星蟲動(dòng)物、紐形動(dòng)物、曳鰓動(dòng)物),占7.2%。可以看出,甌江口海域與洞頭列島大型底棲動(dòng)物種類組成占比例最多的均為多毛類、甲殼動(dòng)物和軟體動(dòng)物。

    2.2 不同海域大型底棲動(dòng)物棲息密度、生物量及次級(jí)生產(chǎn)力

    甌江口海域大型底棲生物年均棲息密度為150.18 個(gè)/m2,年均生物量為2.06 g/m2。洞頭列島海域大型底棲生物年均棲息密度為214.29 個(gè)/m2,年均生物量為4.96 g/m2。以去灰干重(AFDW)計(jì),Brey經(jīng)驗(yàn)公式[3]估算甌江口和洞頭列島海域年均次級(jí)生產(chǎn)力分別為0.48 g(AFDW) m-2a-1和0.70 g(AFDW) m-2a-1,Brey模型[15]估算甌江口和洞頭列島海域年均次級(jí)生產(chǎn)力分別為0.51 g(AFDW) m-2a-1和1.55 g(AFDW) m-2a-1,Brey模型[1]估算甌江口和洞頭列島海域年均次級(jí)生產(chǎn)力分別為0.25 g(AFDW) m-2a-1和0.99 g(AFDW) m-2a-1(表3)。Brey經(jīng)驗(yàn)公式[3]與Brey估算模型[15]對(duì)甌江口海域大型底棲動(dòng)物次級(jí)生產(chǎn)力估算結(jié)果基本一致,Brey模型[1]估算結(jié)果比前兩者的平均值低49.5%;Brey經(jīng)驗(yàn)公式[3]與Brey估算模型[1]對(duì)洞頭列島海域大型底棲動(dòng)物次級(jí)生產(chǎn)力估算結(jié)果基本一致,Brey模型[15]估算結(jié)果比前兩者的平均值高45.5%。3種模型估算結(jié)果均表明,洞頭列島海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力高于甌江口海域。

    表3 甌江口海域和洞頭列島海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力

    AFDW:去灰干重,ash-free dry weight

    圖2 大型底棲動(dòng)物年均生物量(a)和Brey經(jīng)驗(yàn)公式[3](b)、Brey估算模型[15](c)、Brey估算模型[1](d)年均次級(jí)生產(chǎn)力空間分布Fig.2 Distribution of the mean biomass (a) and productivity of macrobenthos estimated by Brey′s empirical formula[3] (b), Brey′s empirical model[15] (c), and Brey′s empirical model[1] (d)

    3種經(jīng)驗(yàn)?zāi)P头从炒笮偷讞珓?dòng)物年均次級(jí)生產(chǎn)力空間分布趨勢(shì)一致,如圖2:甌江口海域有1個(gè)高值區(qū),出現(xiàn)在甌江口南部海域;洞頭列島海域有2個(gè)高值區(qū),分別出現(xiàn)在洞頭本島東北部的三盤島與花崗島之間海域和東黃島東南部的島嶼之間海域,與大型底棲動(dòng)物平均生物量分布規(guī)律相似。根據(jù)Brey經(jīng)驗(yàn)公式[3]、Brey估算模型[15]和Brey估算模型[1]計(jì)算結(jié)果,年均次級(jí)生產(chǎn)力最低值均出現(xiàn)在甌江口北部的S14站位,分別為0.01、0.04 g(AFDW) m-2a-1和0.01 g(AFDW) m-2a-1,原因是該站位大型底棲動(dòng)物物種數(shù)少,生物量極低,僅有少量多毛類。Brey經(jīng)驗(yàn)公式[3]計(jì)算結(jié)果表明,甌江口南部海域S27站位大型底棲動(dòng)物年均次級(jí)生產(chǎn)力最高,為2.62 g(AFDW) m-2a-1,原因是該站位出現(xiàn)了生物量和棲息密度均較大的焦河籃蛤Potamocorbulaustulata。該站位焦河籃蛤個(gè)體大,數(shù)量多,3種經(jīng)驗(yàn)?zāi)P陀?jì)算表明該種對(duì)該站位年均次級(jí)生產(chǎn)力貢獻(xiàn)率均達(dá)到90%以上。Brey估算模型[15]和Brey估算模型[1]計(jì)算結(jié)果均表明S11站位大型底棲動(dòng)物年均次級(jí)生產(chǎn)力最高,分別為6.25 g(AFDW) m-2a-1和4.97 g(AFDW) m-2a-1,原因是該站位大型底棲動(dòng)物物種數(shù)多,達(dá)28種,棲息密度和生物量均較大。

    2.4 大型底棲動(dòng)物年均次級(jí)生產(chǎn)力的主要貢獻(xiàn)物種

    3種經(jīng)驗(yàn)?zāi)P凸浪愦笮偷讞珓?dòng)物年均次級(jí)生產(chǎn)力排前5位的物種基本一致(表4和表5),包括多毛類、軟體動(dòng)物、甲殼動(dòng)物、魚類和紐形動(dòng)物。甌江口海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力主要貢獻(xiàn)種為雙鰓內(nèi)卷齒蠶Aglaophamusdibranchis、長(zhǎng)吻沙蠶Glycerachirori、異蚓蟲Heteromastusfiliforms、焦河籃蛤P.ustulata、薄云母蛤Yoldiasimilis、圓筒原盒螺Eocylichnabraunsi,它們對(duì)該海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力的貢獻(xiàn)率超54.2%;洞頭列島海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力主要貢獻(xiàn)種為雙鰓內(nèi)卷齒蠶、異蚓蟲、薄云母蛤、絨毛細(xì)足蟹Raphidopusciliatus、隆線強(qiáng)蟹Eucratecrenata、腦紐蟲Cerebratulinasp.、紅狼牙鰕虎魚Odontamblyopusrubicundus,它們對(duì)該海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力的貢獻(xiàn)率超57.1%。

    表43種經(jīng)驗(yàn)?zāi)P凸浪惝T江口海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力排前5位的物種

    Table4Top5speciescontributetosecondaryproductivityofmacrobenthosinOujiangRiverestuaryestimatedbythreeempiricalmodels,respectively

    經(jīng)驗(yàn)?zāi)P虴mpiricalmodels次級(jí)生產(chǎn)力排前5的物種Thetop5speciesofsecondaryproductivity合計(jì)貢獻(xiàn)率ThetotalcontributionrateBrey經(jīng)驗(yàn)公式[3]Brey'sempiricalformula[3]焦河籃蛤、薄云母蛤、長(zhǎng)吻沙蠶、圓筒原盒螺、異蚓蟲82.9%Brey估算模型[15]Brey'sempiricalmodel[15]焦河籃蛤、長(zhǎng)吻沙蠶、薄云母蛤、異蚓蟲、雙鰓內(nèi)卷齒蠶58.1%Brey估算模型[1]Brey'sempiricalmodel[1]焦河籃蛤、長(zhǎng)吻沙蠶、薄云母蛤、異蚓蟲、雙鰓內(nèi)卷齒蠶54.2%

    其中,雙鰓內(nèi)卷齒蠶、異蚓蟲、薄云母蛤、絨毛細(xì)足蟹、圓筒原盒螺、長(zhǎng)吻沙蠶屬于分布范圍較廣、棲息密度和生物量均較高的物種,而隆線強(qiáng)蟹、焦河籃蛤?qū)儆趥€(gè)別站位棲息密度和生物量極高,在其他站位棲息密度和生物量偏低的物種。紅狼牙鰕虎魚、腦紐蟲則屬于出現(xiàn)頻率和棲息密度較低,但個(gè)體大、生物量高的物種。

    表53種經(jīng)驗(yàn)?zāi)P凸浪愣搭^列島海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力排前5位的物種

    Table5Top5speciescontributetosecondaryproductivityofmacrobenthosinDongtouIslandsestimatedbythreeempiricalmodels,respectively

    經(jīng)驗(yàn)?zāi)P虴mpiricalmodels次級(jí)生產(chǎn)力排前5的物種Thetop5speciesofsecondaryproductivity合計(jì)貢獻(xiàn)率/%ThetotalcontributionrateBrey經(jīng)驗(yàn)公式[3]Brey'sempiricalformula[3]薄云母蛤、異蚓蟲、腦紐蟲、紅狼牙鰕虎魚、絨毛細(xì)足蟹57.1Brey估算模型[15]Brey'sempiricalmodel[15]異蚓蟲、腦紐蟲、雙鰓內(nèi)卷齒蠶、絨毛細(xì)足蟹、紅狼牙鰕虎魚57.7Brey估算模型[1]Brey'sempiricalmodel[1]腦紐蟲、紅狼牙鰕虎魚、絨毛細(xì)足蟹、異蚓蟲、隆線強(qiáng)蟹59.2

    2.5 大型底棲動(dòng)物年均次級(jí)生產(chǎn)力與環(huán)境因子的關(guān)系

    將3種經(jīng)驗(yàn)?zāi)P凸浪愦笮偷讞珓?dòng)物年均次級(jí)生產(chǎn)力與環(huán)境因子進(jìn)行Pearson相關(guān)性分析。由表6和表7可知,3種經(jīng)驗(yàn)?zāi)P凸浪隳昃渭?jí)生產(chǎn)力之間相關(guān)性極顯著(P<0.01)。甌江口海域Brey模型[15]估算結(jié)果與化學(xué)需氧量、懸浮物顯著負(fù)相關(guān)(P<0.05),Brey模型[1]估算結(jié)果與懸浮物顯著負(fù)相關(guān)(P<0.05)、與沉積物中值粒徑顯著正相關(guān)(P<0.05)。洞頭列島海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力與環(huán)境因子相關(guān)關(guān)系不顯著。

    表6 3種經(jīng)驗(yàn)?zāi)P凸浪惝T江口海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力與環(huán)境因子的相關(guān)性

    表7 3種經(jīng)驗(yàn)?zāi)P凸浪愣搭^列島海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力與環(huán)境因子的相關(guān)性

    3 討論

    3.1 次級(jí)生產(chǎn)力分布格局分析及與其他海域的比較

    次級(jí)生產(chǎn)力的空間分布主要受大型底棲動(dòng)物群落結(jié)構(gòu)和棲息地環(huán)境條件影響[23]。甌江口除南部以外的其他海域,年均次級(jí)生產(chǎn)力較低。河口海域由于受入海河流的影響,沉積速率較高,使底質(zhì)處于強(qiáng)烈的擾動(dòng)之中,在一定程度限制某些大型底棲動(dòng)物的生存,如多毛類和部分甲殼類,活動(dòng)能力強(qiáng),對(duì)底質(zhì)變化不敏感的物種才能生存,造成該區(qū)域大型底棲動(dòng)物種類較少,因此次級(jí)生產(chǎn)力一般較低[7,14]。如遼東灣西部海域[14]受六股河影響,大型底棲動(dòng)物的次級(jí)生產(chǎn)力有由近岸海域向遠(yuǎn)岸海域遞增的趨勢(shì);受長(zhǎng)江徑流影響,河口區(qū)海底環(huán)境極不穩(wěn)定,而距離河口較遠(yuǎn)的海域沉積環(huán)境趨于穩(wěn)定,因此長(zhǎng)江口[7]大型底棲動(dòng)物次級(jí)生產(chǎn)力具有自長(zhǎng)江入??谙驏|遞增的趨勢(shì);越靠近深圳河河口[9],大型底棲動(dòng)物的次級(jí)生產(chǎn)力越低。另外,近年來(lái)實(shí)施的溫州淺灘圍墾促淤工程在靈昆島東側(cè)到洞頭霓嶼島西側(cè)的灘涂上,建成南北兩條圍涂大堤,把靈昆島與霓嶼島連成一片,該海域大部分為新淤出底質(zhì),大型底棲動(dòng)物物種數(shù)少、生物量低。長(zhǎng)江口灘涂圍墾對(duì)大型底棲動(dòng)物的影響研究[24]表明,圍墾促淤可導(dǎo)致大型底棲動(dòng)物物種減少,群落結(jié)構(gòu)改變,生物量明顯降低。因此,圍墾促淤可能導(dǎo)致該海域次級(jí)生產(chǎn)力降低。另外,焦河籃蛤是咸淡水雙殼類動(dòng)物,為河口常見種,對(duì)環(huán)境適應(yīng)能力較強(qiáng)。甌江口南部海域年均次級(jí)生產(chǎn)力高的原因與廈門鼓浪嶼附近海域[25]較為相似,鼓浪嶼位于九龍江河口區(qū),附近海域因出現(xiàn)大量光滑河籃蛤Potamocorbulalaevis導(dǎo)致次級(jí)生產(chǎn)力很高。

    洞頭本島東北部的三盤島與花崗島之間海域和東黃島東南部的島嶼之間海域因距離河口較遠(yuǎn),沉積環(huán)境相對(duì)穩(wěn)定,適宜于大型底棲動(dòng)物生存。另外,島嶼生態(tài)系統(tǒng)有較高的生產(chǎn)力,海流流過島嶼時(shí)產(chǎn)生海底局部上升流使有機(jī)質(zhì)均勻混合,為底棲動(dòng)物提供豐富餌料,且該局部上升流或海流分支可為底棲動(dòng)物提供多樣性的棲息環(huán)境[6],島嶼之間的海域,大型底棲動(dòng)物物種數(shù)、生物量一般較高[26]。因此,該海域大型底棲動(dòng)物次級(jí)生產(chǎn)力較高。

    與其他海域比較,甌江口和洞頭列島海域大型底棲動(dòng)物次級(jí)生產(chǎn)力處于較低水平,低于象山港港口海域[27]0.99 g(AFDW) m-2a-1、膠州灣海域[8]3.41 g(AFDW) m-2a-1、長(zhǎng)江口及鄰近海域[7]3.52 g(AFDW) m-2a-1及遼東灣西部海域[14]3.82 g(AFDW) m-2a-1的調(diào)查結(jié)果。歷史調(diào)查資料[28]顯示,洞頭列島海域?qū)俚蜕锪?、低密度海?該海域是浙江大型底棲動(dòng)物年平均生物量最低的海域,棲息密度僅高于六橫島海域。河口海域由于底質(zhì)不穩(wěn)定可能次級(jí)生產(chǎn)力本身較低,如象山港口門附近站位[27]次級(jí)生產(chǎn)力僅為0.32 g(AFDW) m-2a-1。另外,甌江口和洞頭列島海域次級(jí)生產(chǎn)力低可能與近年來(lái)圍墾促淤有關(guān),如研究表明[25]九龍江口海域底質(zhì)受疏浚擾動(dòng),不利于大型底棲動(dòng)物棲息,年次級(jí)生產(chǎn)力僅為0.31 g(AFDW) m-2a-1。而且相對(duì)依據(jù)全年資料估算的次級(jí)生產(chǎn)力數(shù)值,本文依據(jù)春秋兩個(gè)季節(jié)資料估算的結(jié)果可能與實(shí)際次級(jí)生產(chǎn)力水平存在偏差。

    3.2 次級(jí)生產(chǎn)力與環(huán)境因子的關(guān)系

    環(huán)境因子如水深、水溫、溶解氧等通過影響動(dòng)物代謝速率和物種生活史,進(jìn)而影響次級(jí)生產(chǎn)力[23],分析影響次級(jí)生產(chǎn)力的環(huán)境因子對(duì)全面理解海洋生態(tài)系統(tǒng)的功能具有重要意義[18]。由大型底棲動(dòng)物次級(jí)生產(chǎn)力估算經(jīng)驗(yàn)?zāi)P蛥?shù)可知,除棲息密度和生物量外,次級(jí)生產(chǎn)力還與水深、底溫密切相關(guān)[1,17-18],二者可分別影響食物質(zhì)量和動(dòng)物代謝速率[18]。許多研究表明[3- 4,8,10-11],大型底棲動(dòng)物次級(jí)生產(chǎn)力隨水深增加而降低。此外,影響大型底棲動(dòng)物次級(jí)生產(chǎn)力的環(huán)境因子還有很多,如底質(zhì)類型、鹽度、溶解氧、沉積物中有機(jī)質(zhì)含量、葉綠素a含量、總氮、總磷[10,12,18]等。袁偉等[10]對(duì)膠州灣西部海域大型底棲動(dòng)物次級(jí)生產(chǎn)力與環(huán)境因子的相關(guān)性分析研究表明,葉綠素a是影響該海域次級(jí)生產(chǎn)力的重要因子。嚴(yán)娟等[12]采用典范對(duì)應(yīng)分析長(zhǎng)江口潮間帶大型底棲動(dòng)物群落次級(jí)生產(chǎn)力與環(huán)境變量的關(guān)系,發(fā)現(xiàn)水體總氮、總磷、鹽度、沉積物有機(jī)質(zhì)、溶解氧、葉綠素a均影響次級(jí)生產(chǎn)力,但不同斷面影響因子不盡一致。

    本研究發(fā)現(xiàn)化學(xué)需氧量、懸浮物、沉積物中值粒徑是影響甌江口海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力的重要環(huán)境因子,環(huán)境因子與洞頭列島海域大型底棲動(dòng)物年均次級(jí)生產(chǎn)力相關(guān)關(guān)系不顯著。河口生態(tài)系統(tǒng)處于海洋與陸地交接處,陸源污染物可通過河口流入海洋。甌江口海域可能由于近年來(lái)沿岸工業(yè)發(fā)展導(dǎo)致該海域有機(jī)污染物含量增加,有機(jī)污染通過微生物生化作用分解和氧化,大量消耗水中溶解氧,使水質(zhì)變黑發(fā)臭[29],進(jìn)而影響大型底棲動(dòng)物的生存。且河口海域的懸浮物顆粒量可影響海水透明度和真光層厚度,從而影響浮游生物光合作用和初級(jí)生產(chǎn)力[30],進(jìn)而影響次級(jí)生產(chǎn)力。通常認(rèn)為沉積物粒度是影響大型底棲動(dòng)物物種分布、棲息密度和生物量的重要環(huán)境因子[31],本研究發(fā)現(xiàn)甌江口海域次級(jí)生產(chǎn)力與沉積物中值粒徑顯著正相關(guān),證實(shí)了這一結(jié)論。島嶼間可能由于地理環(huán)境復(fù)雜,相關(guān)性分析結(jié)果表明環(huán)境因子與洞頭列島海域大型底棲動(dòng)物次級(jí)生產(chǎn)力相關(guān)關(guān)系不顯著。不同海域環(huán)境因素不同,且各個(gè)環(huán)境因子對(duì)底棲動(dòng)物不同類群的作用也不同,要深刻理解影響大型底棲動(dòng)物次級(jí)生產(chǎn)力的關(guān)鍵因素,需要對(duì)特定研究海域進(jìn)行長(zhǎng)期數(shù)據(jù)積累和對(duì)比分析,該領(lǐng)域有待進(jìn)一步深入研究。

    3.3 次級(jí)生產(chǎn)力估算經(jīng)驗(yàn)?zāi)P偷谋容^

    次級(jí)生產(chǎn)力估算經(jīng)驗(yàn)?zāi)P褪菍?duì)單個(gè)種群次級(jí)生產(chǎn)力研究的總結(jié),極大地促進(jìn)該領(lǐng)域研究的進(jìn)步,使大型底棲動(dòng)物群落次級(jí)生產(chǎn)力估算成為可能。由Brey經(jīng)驗(yàn)公式[3]到Brey估算模型[15]再到Brey估算模型[1],模型越來(lái)越復(fù)雜、精細(xì),表現(xiàn)在兩方面:一是模型考慮的參數(shù)由2個(gè)增加到20個(gè),并考慮環(huán)境因子;二是濕重生物量、去灰干重、能值轉(zhuǎn)換由所有底棲動(dòng)物相同轉(zhuǎn)化比例精細(xì)到不同類群不同轉(zhuǎn)換比例。另外,從建模方法看,Brey經(jīng)驗(yàn)公式[3]和Brey估算模型[15]建模采用多元線性回歸法,而Brey估算模型[1]采用人工神經(jīng)網(wǎng)絡(luò)法(Artificial Neural Network, ANN)。多元線性回歸法對(duì)獨(dú)立變量間的交互性非常敏感,而人工神經(jīng)網(wǎng)絡(luò)法可以模擬自變量和因變量間復(fù)雜的、非線性的、間斷性的關(guān)系[1]。在預(yù)測(cè)大型底棲動(dòng)物P/B值上,人工神經(jīng)網(wǎng)絡(luò)法優(yōu)于多元線性回歸法[32]。

    經(jīng)驗(yàn)?zāi)P椭皇菍?duì)群落次級(jí)生產(chǎn)力的估算,并非真實(shí)值。3種經(jīng)驗(yàn)?zāi)P退么渭?jí)生產(chǎn)力估算值存在一定差異,但在反映空間分布趨勢(shì)和計(jì)算大型底棲動(dòng)物次級(jí)生產(chǎn)力排前5的物種上大致相同。由于群落次級(jí)生產(chǎn)力的真實(shí)值較難測(cè)定,可將不同估算模型結(jié)果相吻合的數(shù)據(jù)認(rèn)為是次級(jí)生產(chǎn)力的可靠估算值[13],本研究中Brey經(jīng)驗(yàn)公式[3]與Brey估算模型[15]對(duì)甌江口海域大型底棲動(dòng)物次級(jí)生產(chǎn)力估算結(jié)果基本一致,Brey經(jīng)驗(yàn)公式[3]與Brey估算模型[1]對(duì)洞頭列島海域大型底棲動(dòng)物次級(jí)生產(chǎn)力估算結(jié)果基本一致,可認(rèn)為Brey經(jīng)驗(yàn)公式[3]和Brey估算模型[15]對(duì)甌江口海域的估算結(jié)果比較可靠,Brey經(jīng)驗(yàn)公式[3]和Brey估算模型[1]對(duì)洞頭列島海域的估算結(jié)果比較可靠。

    致謝:感謝劉小涯、于培松、張海峰、蔡小霞、趙晨輝、陳文深等協(xié)助采樣并提供相關(guān)理化參數(shù)。

    [1] Brey T. A multi-parameter artificial neural network model to estimate macrobenthic invertebrate productivity and production. Limnology and Oceanography, 2012, 10(8): 581- 589.

    [2] Dolbeth M, Cusson M, Sousa R, Pardal M A. Secondary production as a tool for better understanding of aquatic ecosystems. Canadian Journal of Fisheries and Aquatic Sciences, 2012, 69(7): 1230- 1253.

    [3] Brey T. Estimating productivity of macrobenthic invertebrates from biomass and mean individual weight. Meeresforsch, 1990, 32(4): 329- 343.

    [4] Cusson M, Bourget E. Global patterns of macroinvertebrate production in marine benthic habitats. Marine Ecology Progress Series, 2005, 297: 1- 14.

    [5] Petracco M, Cardoso R S, Corbisier T N, Turra A. Secondary production of sandy beach macrofauna: an evaluation of predictive models. Estuarine, Coastal and Shelf Science, 2012, 115: 359- 365.

    [6] 壽鹿, 曾江寧, 廖一波, 趙永強(qiáng), 江志兵, 陳全震, 高愛根, 楊家新. 甌江口海域大型底棲動(dòng)物分布及其與環(huán)境的關(guān)系. 應(yīng)用生態(tài)學(xué)報(bào), 2009, 20(8): 1958- 1964.

    [7] 劉勇, 線薇薇, 孫世春, 吳耀泉. 長(zhǎng)江口及其鄰近海域大型底棲動(dòng)物生物量、豐度和次級(jí)生產(chǎn)力的初步研究. 中國(guó)海洋大學(xué)學(xué)報(bào), 2008, 38(5): 749- 756.

    [8] 李新正, 王洪法, 張寶琳. 膠州灣大型底棲動(dòng)物次級(jí)生產(chǎn)力初探. 海洋與湖沼, 2005, 36(6): 527- 533.

    [9] 金亮, 蔡立哲, 周細(xì)平, 徐華林. 深圳灣北岸泥灘大型底棲動(dòng)物次級(jí)生產(chǎn)力研究. 臺(tái)灣海峽, 2007, 26(3): 415- 421.

    [10] 袁偉, 張志南, 于子山. 膠州灣西部海域大型底棲動(dòng)物次級(jí)生產(chǎn)力初步研究. 應(yīng)用生態(tài)學(xué)報(bào), 2007, 18(1): 145- 150.

    [11] 王宗興, 孫丕喜, 劉彩霞, 張新軍, 代欣欣, 周穎敏. 桑溝灣大型底棲動(dòng)物的次級(jí)生產(chǎn)力. 應(yīng)用與環(huán)境生物學(xué)報(bào), 2011, 17(4): 495- 498.

    [12] 嚴(yán)娟, 莊平, 侯俊利, 張濤, 劉鑒毅, 馮廣朋, 楊剛. 長(zhǎng)江口潮間帶大型底棲動(dòng)物次級(jí)生產(chǎn)力及其影響因子. 應(yīng)用與環(huán)境生物學(xué)報(bào), 2012, 18(6): 935- 942.

    [13] 張崇良, 徐賓鐸, 任一平, 薛瑩, 紀(jì)毓鵬. 膠州灣潮間帶大型底棲動(dòng)物次級(jí)生產(chǎn)力的時(shí)空變化. 生態(tài)學(xué)報(bào), 2011, 31(17): 5071- 5080.

    [14] 紀(jì)瑩璐, 趙寧, 楊傳平, 季相星, 王振鐘, 隋吉星, 曲方圓, 于子山. 遼東灣西部海域大型底棲動(dòng)物群落次級(jí)生產(chǎn)力初探. 中國(guó)海洋大學(xué)學(xué)報(bào), 2015, 45(4): 53- 58.

    [15] Brey T. Population dynamics in benthic invertebrates. A virtual handbook. 2001. http://www.thomas-brey.de/science/virtualhandbook/navlog/index.html.

    [16] Robertson A I. The relationship between annual production: biomass ratios and lifespans for marine macrobenthos. Oecologia, 1979, 38(2): 193- 202.

    [17] Schwinghamer P, Hargrave B, Peer D, Hawkins C M. Partitioning of production and respiration among size groups of organisms in an intertidal benthic community. Marine Ecology Progress Series, 1986, 31(2): 131- 142.

    [18] Tumbiolo M L, Downing J A. An empirical model for the prediction of secondary production in marine benthic invertebrate populations. Marine Ecology Progress Series, 1994, 114: 165- 174.

    [19] Dolbeth M, Lilleb? A I, Cardoso P G, Ferreira S M, Pardal M A. Annual production of estuarine fauna in different environmental conditions: an evaluation of the estimation methods. Journal of Experimental Marine Biology and Ecology, 2005, 326(2): 115- 127.

    [20] 國(guó)家海洋局. GB/T 12763.6—2007 海洋調(diào)查規(guī)范 第6部分: 海洋生物調(diào)查. 北京: 中國(guó)標(biāo)準(zhǔn)出版社, 2008.

    [21] 國(guó)家海洋局. GB 17378.4—2007 海洋監(jiān)測(cè)規(guī)范 第4部分: 海水分析. 北京: 中國(guó)標(biāo)準(zhǔn)出版社, 2008.

    [22] Soliman Y S, Rowe G T. Secondary production ofAmpeliscamississippianaSoliman and Wicksten 2007 (Amphipoda, Crustacea) in the head of the Mississippi Canyon, northern Gulf of Mexico. Deep Sea Research Part II: Topical Studies in Oceanography, 2008, 55(24/26): 2692- 2698.

    [23] Fuhrmann M M, Pedersen T, Ramasco V, Nilssen E M. Macrobenthic biomass and production in a heterogenic subarctic fjord after invasion by the red king crab. Journal of Sea Research, 2015, 106: 1- 13.

    [24] 馬長(zhǎng)安, 徐霖林, 田偉, 呂巍巍, 趙云龍. 圍墾對(duì)南匯東灘濕地大型底棲動(dòng)物的影響. 生態(tài)學(xué)報(bào), 2012, 32(4): 1007- 1015.

    [25] 周細(xì)平, 蔡立哲, 梁俊彥, 傅素晶, 林和山, 黃昆. 廈門海域大型底棲動(dòng)物次級(jí)生產(chǎn)力的初步研究. 廈門大學(xué)學(xué)報(bào): 自然科學(xué)版, 2008, 47(6): 902- 906.

    [26] 高愛根, 楊俊毅, 陳全震, 胡錫鋼, 董永庭, 胡月妹. 達(dá)山島、平島、車牛山島鄰近海域大型底棲生物分布特征. 海洋學(xué)報(bào), 2003, 25(6): 135- 141.

    [27] 王金輝, 楊春旺, 孫亞偉, 楊元利. 象山港大型底棲動(dòng)物的生物多樣性和次級(jí)生產(chǎn)力研究. 天津農(nóng)學(xué)院學(xué)報(bào), 2006, 13(2): 24- 28.

    [28] 國(guó)家海洋局第二海洋研究所. 浙江省海島海域生物資源綜合調(diào)查報(bào)告. 杭州: 國(guó)家海洋局第二海洋研究所, 1993.

    [29] 鄭青, 韓海波, 周保學(xué), 李金花, 白晶, 蔡偉民. 化學(xué)需氧量(COD)快速測(cè)定新方法研究進(jìn)展. 科學(xué)通報(bào), 2009, 54(21): 3241- 3250.

    [30] 趙冬梅, 高愛國(guó), 李超, 郭占榮, 張延頗. 閩江下游及河口海域水體中懸浮物的分布特征. 海洋地質(zhì)前沿, 2012, 28(6): 35- 39.

    [31] Nanami A, Saito H, Akita T, Motomatsu K I, Kuwahara H. Spatial distribution and assemblage structure of macrobenthic invertebrates in a brackish lake in relation to environmental variables. Estuarine, Coastal and Shelf Science, 2005, 63(1/2): 167- 176.

    [32] Brey T, Jarre-Teichmann A, Borlich O. Artificial neural network versus multiple linear regression: predictingP/Bratios from empirical data. Marine Ecology Progress Series, 1996, 140: 251- 256.

    Secondaryproductivityofmacrobenthosintypicalcoastalecosystems:OujiangRiverestuaryandDongtouIslandsasexamples

    JIA Shenghua1,2, LIAO Yibo1,2, ZENG Jiangning1,2,*, SHOU Lu1,2, HUANG Wei1,2, GAO Aigen1,2, TANG Yanbin1,2

    1KeylaboratoryofMarineEcosystemandBiogeochemistry,StateOceanicAdministration,Hangzhou310012,China2SecondInstituteofOceanography,StateOceanicAdministration,Hangzhou310012,China

    To better understand and estimate the secondary productivity of macrobenthos in estuary ecosystems and island ecosystems, we conducted investigations on the macrobenthos in two typical coastal ecosystems, Oujiang River estuary and Dongtou Islands, in the Eastern China sea, in April and October 2015. In Oujiang River estuary, 81 species were found, including 39 polychaete species, 23 crustacean species, 9 mollusk species, 4 echinodermata species, and 6 other species. In Dongtou Islands, 111 species were found, including 50 polychaete species, 30 crustacean species, 18 mollusk species, 5 echinodermata species, and 8 other species. Three different empirical formulas were applied to calculate the secondary productivity of these macrobenthos. The relationships between the secondary productivity and environmental factors were also analyzed. The results showed that: (1) The average annual secondary productivity of macrobenthos, measured as ash-free dry weight and calculated using three empirical formulas, i.e., Brey′s empirical formula (1990) and Brey′s empirical models (2001 and 2012), were 0.48, 0.51, and 0.25 g m-2a-1, respectively, in the Oujiang River estuary and 0.70, 1.55, and 0.99 g m-2a-1, respectively, in the Dongtou Islands. (2) The spatial distribution of the average annual secondary productivity estimated by the three empirical formulas was consistent. A zone with high average annual secondary productivity of macrobenthos was found in the southern waters of the Oujiang River estuary. Two such zones were found in the Dongtou Islands, one located in the sea area between Sanpan Island and Huagang Island, northeast of the main island of Dongtou, and the other in the sea area between the islands southeast of the Donghuang Island. (3) Six species,Aglaophamusdibranchis,Glycerachirori,Heteromastusfiliformis,Potamocorbulaustulata,Yoldiasimilis, andEocylichnabraunsi, mainly contributed—more than 54.2% of the total—to the average annual secondary productivity of macrobenthos, in the Oujiang River estuary. In the Dongtou Islands,A.dibranchis,H.filiformis,Y.similis,Raphidopusciliatus,Eucratecrenata,Cerebratulinasp., andOdontamblyopusrubicunduscontributed the most—more than 57.1% of the total—to the average annual secondary productivity of macrobenthos. (4) Correlation analysis indicated that chemical oxygen demand, suspended solids, and median size of surface sediments were the important environmental factors affecting the average annual secondary productivity of macrobenthos in the Oujiang River estuary. In contrast, environmental factors had no measurable impact on the average annual secondary productivity of macrobenthos in the Dongtou Islands. (5) The average annual secondary productivity of macrobenthos estimated by Brey′s empirical formula (1990) was consistent with that estimated by Brey′s empirical model (2001), for the Oujiang River estuary, and with that estimated by Brey′s empirical model (2012), for the Dongtou Islands.

    secondary productivity; macrobenthos; Dongtou Islands; Oujiang River estuary; Brey′s empirical model

    海洋公益性行業(yè)科研專項(xiàng)(201405007,201505004- 3,201505027- 4);國(guó)家海洋局第二海洋研究所基本科研業(yè)務(wù)費(fèi)專項(xiàng) (JG1311,JG1312,JG1412,JG1616);浙江省自然科學(xué)基金(LY14D060007);國(guó)家科技支撐計(jì)劃(2015BAD08B01)

    2016- 09- 18; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 07- 11

    *通訊作者Corresponding author.E-mail: jiangningz@126.com

    10.5846/stxb201609181872

    賈勝華,廖一波,曾江寧,壽鹿,黃偉,高愛根,湯雁濱.近岸典型生態(tài)系統(tǒng)大型底棲動(dòng)物的次級(jí)生產(chǎn)力——以甌江口、洞頭列島為例.生態(tài)學(xué)報(bào),2017,37(21):7140- 7151.

    Jia S H, Liao Y B, Zeng J N, Shou L, Huang W, Gao A G, Tang Y B.Secondary productivity of macrobenthos in typical coastal ecosystems: Oujiang River estuary and Dongtou Islands as examples.Acta Ecologica Sinica,2017,37(21):7140- 7151.

    猜你喜歡
    列島洞頭甌江
    A new species of the gudgeon genus Microphysogobio Mori, 1934 (Cypriniformes: Cyprinidae) from Zhejiang Province, China
    洞頭:三無(wú)村變身“桃花源”
    溫州甌江北口大橋綠色公路創(chuàng)建實(shí)踐
    上海公路(2019年3期)2019-11-25 07:39:26
    《甌江叢韻》在浙首發(fā)
    “意外”的插曲撞擊出數(shù)學(xué)思維“火花”
    東方教育(2017年9期)2017-07-19 14:12:36
    日本列島封土石室墓的展開與連云港石室墓
    東南文化(2016年4期)2016-08-22 09:35:42
    洞頭近海養(yǎng)殖水域沉積物碳含量空間分布格局
    舟山風(fēng)物志
    排水固結(jié)法在甌江口新區(qū)起步區(qū)道路軟基處理中的應(yīng)用
    長(zhǎng)山列島國(guó)家地質(zhì)公園
    国产精品综合久久久久久久免费| 亚洲人成网站在线播放欧美日韩| 97碰自拍视频| 免费在线观看完整版高清| 香蕉丝袜av| 国产熟女午夜一区二区三区| 高清在线国产一区| 在线十欧美十亚洲十日本专区| 午夜福利在线在线| 变态另类丝袜制服| 91国产中文字幕| 午夜福利成人在线免费观看| 99国产精品99久久久久| 国产免费男女视频| 日本一本二区三区精品| 波多野结衣高清无吗| 久久久久久免费高清国产稀缺| 国产午夜精品论理片| 欧美一级a爱片免费观看看 | 巨乳人妻的诱惑在线观看| 人人妻,人人澡人人爽秒播| 久99久视频精品免费| 青草久久国产| 亚洲欧美一区二区三区黑人| 久久中文字幕人妻熟女| 99久久国产精品久久久| 精品国内亚洲2022精品成人| 色哟哟哟哟哟哟| 亚洲,欧美精品.| 亚洲精品久久国产高清桃花| www.熟女人妻精品国产| 久久久久国产精品人妻aⅴ院| 日韩欧美国产一区二区入口| 亚洲成人精品中文字幕电影| 亚洲人成网站在线播放欧美日韩| 男男h啪啪无遮挡| 啦啦啦韩国在线观看视频| 精品久久久久久久人妻蜜臀av| 99久久久亚洲精品蜜臀av| 亚洲男人天堂网一区| 亚洲五月天丁香| 99久久综合精品五月天人人| 国产熟女xx| 男人舔奶头视频| 哪里可以看免费的av片| 最近在线观看免费完整版| 最新美女视频免费是黄的| 国产精品1区2区在线观看.| 啦啦啦免费观看视频1| 香蕉丝袜av| 亚洲成a人片在线一区二区| 真人做人爱边吃奶动态| 我的老师免费观看完整版| 宅男免费午夜| 嫁个100分男人电影在线观看| 欧洲精品卡2卡3卡4卡5卡区| 成人国产综合亚洲| 99re在线观看精品视频| 亚洲欧美一区二区三区黑人| 精品午夜福利视频在线观看一区| 亚洲专区字幕在线| 午夜影院日韩av| 欧美日韩中文字幕国产精品一区二区三区| 久99久视频精品免费| 啦啦啦韩国在线观看视频| 精品国产亚洲在线| 免费在线观看完整版高清| 久久天躁狠狠躁夜夜2o2o| 欧美日韩亚洲国产一区二区在线观看| 又黄又爽又免费观看的视频| 一区福利在线观看| 国产成人精品无人区| 久久草成人影院| 久久精品亚洲精品国产色婷小说| 亚洲成人精品中文字幕电影| 亚洲七黄色美女视频| 国产v大片淫在线免费观看| 成人永久免费在线观看视频| av国产免费在线观看| 一夜夜www| 欧美中文日本在线观看视频| 欧美日韩瑟瑟在线播放| 亚洲成人中文字幕在线播放| 色播亚洲综合网| 老司机福利观看| 日韩国内少妇激情av| 久久久国产欧美日韩av| cao死你这个sao货| 日韩欧美国产在线观看| 亚洲中文字幕一区二区三区有码在线看 | 精品久久久久久久人妻蜜臀av| 国内少妇人妻偷人精品xxx网站 | 一进一出抽搐动态| 黄色毛片三级朝国网站| 久久精品人妻少妇| 亚洲18禁久久av| 国产成人精品无人区| 午夜免费观看网址| 免费av毛片视频| 亚洲,欧美精品.| 少妇裸体淫交视频免费看高清 | 国产精品一区二区免费欧美| 日韩欧美国产在线观看| 欧美激情久久久久久爽电影| 亚洲成人中文字幕在线播放| 亚洲国产欧美一区二区综合| 亚洲第一电影网av| 久久精品91蜜桃| 久久精品影院6| 在线观看日韩欧美| 精品高清国产在线一区| 国产三级在线视频| 制服丝袜大香蕉在线| 少妇的丰满在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩高清综合在线| 91老司机精品| ponron亚洲| a级毛片a级免费在线| 国产区一区二久久| 亚洲色图av天堂| 精品国产美女av久久久久小说| 男插女下体视频免费在线播放| 夜夜躁狠狠躁天天躁| 国产三级中文精品| 一个人免费在线观看电影 | 白带黄色成豆腐渣| 美女黄网站色视频| 亚洲色图 男人天堂 中文字幕| 人妻久久中文字幕网| 亚洲人成伊人成综合网2020| 一级片免费观看大全| 日韩高清综合在线| 欧美人与性动交α欧美精品济南到| 成年人黄色毛片网站| 国产区一区二久久| 免费在线观看亚洲国产| 韩国av一区二区三区四区| 日韩 欧美 亚洲 中文字幕| 18美女黄网站色大片免费观看| 91老司机精品| 午夜影院日韩av| 亚洲国产欧美一区二区综合| 一级毛片女人18水好多| 亚洲精品在线观看二区| 大型av网站在线播放| 级片在线观看| 国产爱豆传媒在线观看 | 国产精品98久久久久久宅男小说| www.999成人在线观看| 99久久国产精品久久久| 国内精品一区二区在线观看| 在线免费观看的www视频| 99精品久久久久人妻精品| 亚洲午夜理论影院| avwww免费| 日本一二三区视频观看| 人人妻人人看人人澡| 国产一区在线观看成人免费| 国产亚洲av嫩草精品影院| 黄色毛片三级朝国网站| 五月伊人婷婷丁香| 久久午夜亚洲精品久久| 一本大道久久a久久精品| 一进一出好大好爽视频| 黄色片一级片一级黄色片| 亚洲自偷自拍图片 自拍| 国产免费av片在线观看野外av| 麻豆成人午夜福利视频| 不卡av一区二区三区| 精品乱码久久久久久99久播| e午夜精品久久久久久久| 亚洲专区国产一区二区| 国产一区二区激情短视频| 欧美丝袜亚洲另类 | 亚洲美女视频黄频| 国产熟女午夜一区二区三区| 又黄又爽又免费观看的视频| 99在线人妻在线中文字幕| www国产在线视频色| 全区人妻精品视频| 国产精品自产拍在线观看55亚洲| 亚洲男人的天堂狠狠| 日本黄大片高清| 丝袜美腿诱惑在线| 观看免费一级毛片| 亚洲乱码一区二区免费版| 18美女黄网站色大片免费观看| 国产1区2区3区精品| 中文字幕人成人乱码亚洲影| 久9热在线精品视频| 国产伦在线观看视频一区| 日韩 欧美 亚洲 中文字幕| АⅤ资源中文在线天堂| 黄片小视频在线播放| 亚洲激情在线av| 国产一区二区三区在线臀色熟女| 免费搜索国产男女视频| 一级作爱视频免费观看| 亚洲欧美精品综合久久99| 又黄又粗又硬又大视频| 欧美黑人欧美精品刺激| 欧美中文日本在线观看视频| 999精品在线视频| 日本免费a在线| 91大片在线观看| 国产精品影院久久| 日韩高清综合在线| 香蕉久久夜色| 五月玫瑰六月丁香| 91麻豆av在线| 12—13女人毛片做爰片一| 亚洲一码二码三码区别大吗| 欧美av亚洲av综合av国产av| 啪啪无遮挡十八禁网站| 亚洲中文日韩欧美视频| 国产精品一及| 99久久精品国产亚洲精品| 国产视频一区二区在线看| 极品教师在线免费播放| 午夜精品一区二区三区免费看| 欧美av亚洲av综合av国产av| av福利片在线| 国产精品1区2区在线观看.| 91成年电影在线观看| 亚洲国产精品久久男人天堂| av天堂在线播放| 亚洲五月婷婷丁香| 亚洲va日本ⅴa欧美va伊人久久| 天堂动漫精品| www.精华液| 日韩免费av在线播放| av超薄肉色丝袜交足视频| 99国产精品99久久久久| 在线免费观看的www视频| 免费看日本二区| 亚洲熟女毛片儿| 免费看美女性在线毛片视频| 夜夜躁狠狠躁天天躁| 亚洲国产看品久久| 久久久精品国产亚洲av高清涩受| 亚洲 欧美 日韩 在线 免费| 欧美日韩乱码在线| 精品国内亚洲2022精品成人| 色综合亚洲欧美另类图片| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线美女| 国产精品 国内视频| 99久久99久久久精品蜜桃| 久久精品国产清高在天天线| 国产单亲对白刺激| 久久婷婷成人综合色麻豆| 亚洲av熟女| 黄色毛片三级朝国网站| 精品乱码久久久久久99久播| 蜜桃久久精品国产亚洲av| 婷婷精品国产亚洲av在线| 嫁个100分男人电影在线观看| 日韩国内少妇激情av| 国产一区二区激情短视频| 久久久久精品国产欧美久久久| 午夜免费成人在线视频| 香蕉国产在线看| 亚洲av五月六月丁香网| 制服丝袜大香蕉在线| 丁香六月欧美| 美女黄网站色视频| √禁漫天堂资源中文www| av中文乱码字幕在线| 国产亚洲av嫩草精品影院| 成人国产综合亚洲| 叶爱在线成人免费视频播放| 国产1区2区3区精品| 在线免费观看的www视频| 一卡2卡三卡四卡精品乱码亚洲| 此物有八面人人有两片| 最近在线观看免费完整版| 在线观看66精品国产| 久久中文字幕人妻熟女| 中文字幕高清在线视频| 国产av不卡久久| 成人18禁在线播放| 美女大奶头视频| 日韩有码中文字幕| 国产一级毛片七仙女欲春2| 色播亚洲综合网| 一本综合久久免费| 久久精品综合一区二区三区| 精品久久久久久久久久免费视频| 久久国产精品影院| 亚洲欧美一区二区三区黑人| 国产熟女xx| 久久久水蜜桃国产精品网| 熟女少妇亚洲综合色aaa.| 2021天堂中文幕一二区在线观| 1024手机看黄色片| 亚洲人成网站在线播放欧美日韩| av超薄肉色丝袜交足视频| 熟妇人妻久久中文字幕3abv| 国产精品爽爽va在线观看网站| 正在播放国产对白刺激| 在线观看一区二区三区| 老司机午夜福利在线观看视频| 国产精品久久久久久亚洲av鲁大| 99精品久久久久人妻精品| 日韩中文字幕欧美一区二区| av超薄肉色丝袜交足视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av电影在线进入| 欧美一级a爱片免费观看看 | 亚洲午夜理论影院| 九九热线精品视视频播放| 免费观看精品视频网站| 美女高潮喷水抽搐中文字幕| 日韩欧美国产在线观看| av欧美777| 首页视频小说图片口味搜索| 日本五十路高清| 欧美人与性动交α欧美精品济南到| 亚洲成人免费电影在线观看| 高清毛片免费观看视频网站| 精品少妇一区二区三区视频日本电影| 久久亚洲真实| 欧美精品亚洲一区二区| 夜夜爽天天搞| 白带黄色成豆腐渣| 国产精品久久久久久人妻精品电影| 女人爽到高潮嗷嗷叫在线视频| videosex国产| 午夜日韩欧美国产| 日韩三级视频一区二区三区| 免费观看精品视频网站| 哪里可以看免费的av片| 啦啦啦韩国在线观看视频| 男女床上黄色一级片免费看| 99热这里只有是精品50| 亚洲精品色激情综合| 高清在线国产一区| 国产欧美日韩一区二区三| 成人亚洲精品av一区二区| 国产三级在线视频| 国产精品亚洲av一区麻豆| 此物有八面人人有两片| 免费在线观看亚洲国产| 女人高潮潮喷娇喘18禁视频| 成年免费大片在线观看| 日韩欧美精品v在线| 老司机午夜十八禁免费视频| 亚洲精品一区av在线观看| 国产成人aa在线观看| 我的老师免费观看完整版| 成年人黄色毛片网站| 亚洲国产高清在线一区二区三| 18禁裸乳无遮挡免费网站照片| 久久国产精品影院| 黄片大片在线免费观看| 国产一区二区三区视频了| 给我免费播放毛片高清在线观看| 国产99白浆流出| 青草久久国产| 亚洲成人久久性| 欧美成狂野欧美在线观看| 日韩中文字幕欧美一区二区| 久久九九热精品免费| 亚洲 欧美 日韩 在线 免费| 亚洲欧美日韩高清专用| 悠悠久久av| 午夜精品久久久久久毛片777| 色哟哟哟哟哟哟| 精品欧美国产一区二区三| 妹子高潮喷水视频| 精品欧美国产一区二区三| 久久婷婷成人综合色麻豆| 亚洲欧美日韩高清专用| 亚洲专区中文字幕在线| 久99久视频精品免费| 亚洲专区中文字幕在线| 在线看三级毛片| 最好的美女福利视频网| 久久草成人影院| 五月伊人婷婷丁香| 色尼玛亚洲综合影院| 久久这里只有精品中国| 丰满人妻熟妇乱又伦精品不卡| 国产免费av片在线观看野外av| 黄色丝袜av网址大全| 1024香蕉在线观看| 怎么达到女性高潮| 九九热线精品视视频播放| 亚洲国产欧美网| 中文字幕熟女人妻在线| 成人永久免费在线观看视频| 欧美zozozo另类| 国产精品国产高清国产av| 成人午夜高清在线视频| 欧美精品亚洲一区二区| 国产亚洲av嫩草精品影院| 身体一侧抽搐| 一进一出抽搐动态| 午夜福利18| 99久久综合精品五月天人人| 在线永久观看黄色视频| 99热这里只有是精品50| 琪琪午夜伦伦电影理论片6080| 午夜日韩欧美国产| а√天堂www在线а√下载| 亚洲国产精品合色在线| 日韩大码丰满熟妇| 每晚都被弄得嗷嗷叫到高潮| 99久久国产精品久久久| 美女大奶头视频| 亚洲男人的天堂狠狠| 精品国产美女av久久久久小说| xxx96com| videosex国产| 99国产精品一区二区三区| 免费人成视频x8x8入口观看| 一二三四在线观看免费中文在| 欧美成人一区二区免费高清观看 | 中文字幕人妻丝袜一区二区| 一级毛片高清免费大全| 在线观看美女被高潮喷水网站 | 最新美女视频免费是黄的| 不卡一级毛片| 一个人免费在线观看的高清视频| 校园春色视频在线观看| 三级国产精品欧美在线观看 | 国产精品美女特级片免费视频播放器 | 久久久久亚洲av毛片大全| 亚洲免费av在线视频| 99re在线观看精品视频| 国产又黄又爽又无遮挡在线| 草草在线视频免费看| 久久草成人影院| 两个人的视频大全免费| 琪琪午夜伦伦电影理论片6080| 给我免费播放毛片高清在线观看| 国产熟女xx| 久久99热这里只有精品18| 国产亚洲欧美98| 人人妻人人澡欧美一区二区| 成人三级做爰电影| 国产欧美日韩一区二区三| 夜夜躁狠狠躁天天躁| 麻豆国产97在线/欧美 | 国产精品日韩av在线免费观看| 成人精品一区二区免费| 51午夜福利影视在线观看| 五月伊人婷婷丁香| 久久久久九九精品影院| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品美女久久av网站| 伦理电影免费视频| 99re在线观看精品视频| 成年免费大片在线观看| 成人av一区二区三区在线看| 毛片女人毛片| 国产精品免费视频内射| 亚洲一码二码三码区别大吗| 国产成+人综合+亚洲专区| 久久精品综合一区二区三区| 免费观看精品视频网站| 岛国视频午夜一区免费看| 国产片内射在线| xxxwww97欧美| 伊人久久大香线蕉亚洲五| 毛片女人毛片| 中文亚洲av片在线观看爽| 岛国在线免费视频观看| 一本大道久久a久久精品| 久久人妻福利社区极品人妻图片| 韩国av一区二区三区四区| 最近视频中文字幕2019在线8| 国产99久久九九免费精品| 日韩欧美精品v在线| 非洲黑人性xxxx精品又粗又长| 在线观看66精品国产| 18禁美女被吸乳视频| 啦啦啦观看免费观看视频高清| 五月玫瑰六月丁香| 国产aⅴ精品一区二区三区波| 欧美黑人精品巨大| 久久精品国产99精品国产亚洲性色| 男女午夜视频在线观看| x7x7x7水蜜桃| 男人舔奶头视频| 黄色片一级片一级黄色片| 国内少妇人妻偷人精品xxx网站 | 别揉我奶头~嗯~啊~动态视频| 人妻夜夜爽99麻豆av| 一本一本综合久久| 色哟哟哟哟哟哟| 一本一本综合久久| 亚洲成人久久性| 成人av在线播放网站| ponron亚洲| 一进一出抽搐动态| 97碰自拍视频| 国产成人aa在线观看| 神马国产精品三级电影在线观看 | 欧美+亚洲+日韩+国产| 日本一二三区视频观看| 999久久久精品免费观看国产| 国产真实乱freesex| 好男人电影高清在线观看| 午夜老司机福利片| 此物有八面人人有两片| 国产视频内射| 村上凉子中文字幕在线| 老司机福利观看| 久久精品aⅴ一区二区三区四区| 黄色丝袜av网址大全| 国产欧美日韩一区二区三| 黑人欧美特级aaaaaa片| 亚洲,欧美精品.| 成人一区二区视频在线观看| 国产午夜精品久久久久久| 亚洲精品av麻豆狂野| 日韩三级视频一区二区三区| 可以在线观看毛片的网站| 亚洲av成人av| 日日爽夜夜爽网站| 哪里可以看免费的av片| 可以在线观看毛片的网站| 亚洲人成77777在线视频| 欧美日本视频| 国产精品自产拍在线观看55亚洲| 香蕉久久夜色| 久久精品国产亚洲av高清一级| 欧美不卡视频在线免费观看 | 欧美性长视频在线观看| 国内毛片毛片毛片毛片毛片| 国产亚洲精品一区二区www| 视频区欧美日本亚洲| 一级毛片女人18水好多| 亚洲精品色激情综合| 国产1区2区3区精品| 国产精品av久久久久免费| 成人国产综合亚洲| 欧美日韩乱码在线| 欧美日韩瑟瑟在线播放| 亚洲aⅴ乱码一区二区在线播放 | 欧美在线一区亚洲| 欧美精品亚洲一区二区| 亚洲av第一区精品v没综合| 国产亚洲精品一区二区www| 舔av片在线| 窝窝影院91人妻| or卡值多少钱| 国产伦在线观看视频一区| 亚洲aⅴ乱码一区二区在线播放 | 亚洲av熟女| 熟女电影av网| 男插女下体视频免费在线播放| 国产一区二区三区视频了| 两性夫妻黄色片| 欧美成人午夜精品| 国产成人一区二区三区免费视频网站| 精品人妻1区二区| 最新在线观看一区二区三区| 欧美av亚洲av综合av国产av| 欧美丝袜亚洲另类 | 伦理电影免费视频| 午夜免费观看网址| 亚洲在线自拍视频| 变态另类丝袜制服| 午夜日韩欧美国产| 日本免费a在线| 日韩精品中文字幕看吧| 精品久久久久久久毛片微露脸| 看片在线看免费视频| 亚洲熟妇熟女久久| 日韩大尺度精品在线看网址| 一个人免费在线观看的高清视频| 国产精品乱码一区二三区的特点| 亚洲精华国产精华精| 久久国产精品影院| 久久九九热精品免费| 又粗又爽又猛毛片免费看| 欧美丝袜亚洲另类 | 亚洲欧美日韩高清在线视频| 黄色视频不卡| 亚洲中文日韩欧美视频| 国产黄色小视频在线观看| 欧美日韩瑟瑟在线播放| 欧美zozozo另类| 欧美精品啪啪一区二区三区| 亚洲一码二码三码区别大吗| 欧美乱码精品一区二区三区| 国产精品自产拍在线观看55亚洲| 欧美不卡视频在线免费观看 | 国语自产精品视频在线第100页| a级毛片在线看网站| 一二三四在线观看免费中文在| 国产真实乱freesex| 五月伊人婷婷丁香| 最新美女视频免费是黄的| 久久久水蜜桃国产精品网| 成人永久免费在线观看视频| 免费人成视频x8x8入口观看| 亚洲成人久久性| 久久精品国产99精品国产亚洲性色| 午夜福利欧美成人| 久久午夜综合久久蜜桃| 久久精品综合一区二区三区| 国产伦一二天堂av在线观看| 大型av网站在线播放| 久久香蕉精品热| 成人手机av| 熟女电影av网| 午夜免费激情av| 欧美激情久久久久久爽电影| 国产精品久久视频播放| 精品少妇一区二区三区视频日本电影| 高清毛片免费观看视频网站|