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

    魚類粒徑譜研究進展

    2017-11-07 02:34:58郭建忠陳作志徐姍楠
    海洋漁業(yè) 2017年5期
    關(guān)鍵詞:正態(tài)斜率魚類

    郭建忠,陳作志,徐姍楠

    魚類粒徑譜研究進展

    郭建忠1,2,陳作志1,徐姍楠1

    (1.中國水產(chǎn)科學(xué)研究院南海水產(chǎn)研究所,農(nóng)業(yè)部南海漁業(yè)資源開發(fā)利用重點實驗室,廣東省漁業(yè)生態(tài)環(huán)境重點實驗室,廣東廣州 510300;2.上海海洋大學(xué)水產(chǎn)與生命學(xué)院,上海 201306)

    粒徑譜表示的是生物量或者生物數(shù)量與粒徑大小關(guān)系的曲線,能夠反映生態(tài)系統(tǒng)結(jié)構(gòu)和功能以及系統(tǒng)內(nèi)部的動態(tài)聯(lián)系,是當(dāng)前水生生態(tài)學(xué)研究的熱點。綜述了魚類粒徑譜的理論、研究方法及影響因子等,重點介紹了魚類粒徑譜在魚類生產(chǎn)量評估,開展生態(tài)系統(tǒng)監(jiān)測,以及預(yù)測各種人類活動對生態(tài)系統(tǒng)的影響等方面的研究進展。針對目前研究中所存在的不足與問題,提出了新觀測和采樣技術(shù)的應(yīng)用以及基于粒徑譜評估魚類生物量的通用理論模式研發(fā)等是今后魚類粒徑譜研究的趨勢。

    粒徑譜;魚類粒徑譜;生態(tài)系統(tǒng)

    粒徑譜表示的是生物量或者生物數(shù)量與粒徑大小關(guān)系的曲線,是研究生態(tài)系統(tǒng)特征的一個重要途徑。它不僅能夠從宏觀上研究海洋生態(tài)系統(tǒng)的結(jié)構(gòu)、動態(tài)和影響機制,而且能夠比較不同生態(tài)系統(tǒng)的特征、評估生物群落的生產(chǎn)力與生物量,對生態(tài)系統(tǒng)的研究起著重要的指示作用[1]。在水生生態(tài)系統(tǒng)中[2],魚類粒徑譜被廣泛應(yīng)用于評估不同營養(yǎng)水平的生產(chǎn)狀況,尤其是魚類的生產(chǎn)、豐度和物種的豐度[3-9],以及預(yù)測各種人類活動對生態(tài)系統(tǒng)的影響等[10-11]。因此,研究魚類粒徑譜對于漁業(yè)的管理、開發(fā)和生態(tài)系統(tǒng)的可持續(xù)健康發(fā)展有著重要的意義。

    1 粒徑譜

    1.1 粒徑譜概念及發(fā)展

    1967年 SHELDON等[12-13]首次提出粒徑譜的概念,即粒徑譜是在相等的對數(shù)粒徑間隔上測量生物體總的豐度或生物量;將某一特定的生態(tài)系統(tǒng)中的生物按照統(tǒng)一顆粒粒徑大小(equivalent spherical diameter,ESD)劃分為不同的粒徑范圍,計算各粒徑范圍上的生物量,并將其對數(shù)轉(zhuǎn)化,那么這種生物量在對數(shù)粒徑范圍上的分布稱為生物量粒徑譜。之后又提出了粒徑譜理論,即相等對數(shù)化粒徑間隔的顆粒物質(zhì)的生物量大致相同。之后的40年間,該理論得到了迅速的發(fā)展[14]。20世紀(jì) 70年代,SHELDON等[8]以捕食者-被捕食者相互作用關(guān)系建立模型來研究海洋生態(tài)系統(tǒng);KERR[15]在 SHELDON模型基礎(chǔ)上以能量流動與營養(yǎng)級關(guān)系來構(gòu)建模型;THIEBAUX等[16-17]以雙對數(shù)坐標(biāo)繪出非正態(tài)化生物量譜,并提出正態(tài)化生物量的概念[13];SILVERT等[18]以時間為基礎(chǔ)建立模型。80年代,BORGMANN[3]提出粒徑轉(zhuǎn)換效率“ε”概念,并用來評估生物體粒徑組的生產(chǎn)力[10]。90年代,BOUDREAU等[19]用生態(tài)動力學(xué)理論來分析生物量譜穹頂與營養(yǎng)級、粒徑的關(guān)系;THIEBAUX等[16]在P-D模型基礎(chǔ)上建立完全異速生長模型;ZHOU等[20]建立正態(tài)化生物量譜來平衡生物體出生率、生長率、死亡率與捕食的關(guān)系。21世紀(jì)初期,ROSEN[21]建立粒徑模型來預(yù)測任意粒級范圍內(nèi)不同粒徑的物種比例;LAW等[22]以特定捕食者粒徑為基礎(chǔ)建立參數(shù)化模型;ROSSBERG[23]通過攝食、新陳代謝和生長的能量消耗等關(guān)系將不同規(guī)格的性成熟魚類物種連接起來建立模型。

    1.2 粒徑譜研究方法

    對于水生生態(tài)系統(tǒng),粒徑譜有兩種主要的研究方法[24]。第一種是假定生物體由小型生長到大型的生物量是連續(xù)的,集中分析生物量粒徑譜的斜率,采用動態(tài)方程式描述水生生物群落特征[25]。第二種是在以粒徑為基礎(chǔ)的生態(tài)系統(tǒng)結(jié)構(gòu)一致的前提下分析生物量粒徑譜的不規(guī)則性和分布結(jié)構(gòu)[26]。兩種方法形成對照,第一種強調(diào)生長,第二種強調(diào)捕食行為,且第二種已經(jīng)廣泛用于分析粒徑譜形成的“峰型”曲線上穹頂形成的原因和研究海洋底棲無脊椎動物、魚類粒徑譜的斜率[27]。

    2 魚類粒徑譜

    2.1 魚類粒徑譜基本原理

    1991年,SPRULES等[9]對美國密歇根湖海洋生物進行研究,首次完成了生物量粒徑譜從細(xì)菌到大型魚類的應(yīng)用。自此,國內(nèi)外對于魚類粒徑譜的研究逐步開展起來。魚類粒徑譜曲線所形成的“峰型”形狀反映魚類群落的結(jié)構(gòu)特征,如圖1安大略湖魚類粒徑譜呈現(xiàn)“單峰”型。正態(tài)化魚類粒徑譜表面上反映魚類群落生物量的轉(zhuǎn)化分布,而本質(zhì)上代表魚類群落的粒級豐度;方法上強調(diào)異速生長比例,其理論斜率為-1,且斜率與海域的生產(chǎn)狀況、魚類粒徑大小、捕撈強度、棲息地環(huán)境等因素有關(guān),表示生物量與粒徑之間的關(guān)系,如圖1安大略湖魚類正態(tài)化粒徑譜斜率為-1.249 8,小于-1,說明生物量隨粒徑的增大而減小。

    2.2 魚類粒徑譜類型

    2.2.1 SHELDON型生物量粒徑譜

    SHELDON型生物量粒徑譜是用以2為底對數(shù)轉(zhuǎn)換的個體干重生物量劃分的粒徑間隔為橫坐標(biāo),用以2為底對數(shù)轉(zhuǎn)換的單位面積上對應(yīng)粒徑間隔的干重生物量為縱坐標(biāo)[26]。SHELDON型生物量粒徑譜曲線上穹頂與波谷構(gòu)成的“峰型”能夠反映不同粒徑范圍內(nèi)魚類群落的結(jié)構(gòu)特征。

    2.2.2 正態(tài)化生物量粒徑譜

    正態(tài)化生物量粒徑譜(normalized biomass size spectrum,NBSS)是用以2為底對數(shù)轉(zhuǎn)換的個體干重生物量劃分的粒級為橫坐標(biāo),用以2為底對數(shù)轉(zhuǎn)換后的單位面積上對應(yīng)粒級總干重生物量除以該粒級上個體干重的變化幅度為縱坐標(biāo)[28]。正態(tài)化生物量粒徑譜包含幾個重要參數(shù)信息:斜率(a)、截距(b)和回歸系數(shù)(R2)。

    正態(tài)化生物量粒徑譜參數(shù)能夠描述魚類群落的結(jié)構(gòu)特征、比較不同魚類群落間的差異和預(yù)測魚類群落未來的發(fā)展?fàn)顩r[13]。a表示魚類群落生物量與粒徑之間的關(guān)系:a=-1,魚類群落生物量隨著粒徑的增大而保持為一個不變的常數(shù);a<-1,魚類群落生物量隨著粒徑的增大而減??;a>-1,魚類群落生物量隨著粒徑的增大而增大。b表示魚類群落的豐度和生產(chǎn)力水平,b值不同則代表不同魚類群落生產(chǎn)力水平之間的差異[26];一般而言,生產(chǎn)力水平的大小與b值的大小呈正相關(guān)[29]。R2表示魚類群落偏離穩(wěn)定狀態(tài)的程度。

    圖1 安大略湖魚類生物量譜,虛線為正態(tài)化生物量譜線(線性項)[13]Fig.1 Fish biomass size spectrum for Lake Ontario with fitted parabolas and the linear“backbone”[13]

    2.3 魚類粒徑譜主要研究方法

    2.3.1 庫爾特計數(shù)法

    1967年庫爾特發(fā)明計數(shù)器以后,庫爾特計數(shù)法開始應(yīng)用于海洋生物粒徑譜的研究[14],根據(jù)小孔電阻原理來測量顆粒大小。通過測量每個顆粒大小,統(tǒng)計出粒徑分布,是最早研究海洋生態(tài)系統(tǒng)粒徑的一種方法。

    2.3.2 拖網(wǎng)捕撈法

    對于魚類群落的研究一般是海上現(xiàn)場調(diào)查,根據(jù)所需魚種類型選擇網(wǎng)型和網(wǎng)目直徑,通過拖網(wǎng)捕撈獲取魚類數(shù)據(jù)開展研究。該方法較為直觀,但也具有一定的局限性。如網(wǎng)目過大會導(dǎo)致小型魚類遺漏,網(wǎng)目過小會導(dǎo)致過度捕撈,都不利于魚類群落的研究和魚類資源合理開發(fā)。MACPHERSON等[30]用拖網(wǎng)調(diào)查對地中海西北部魚類群落研究表明,受保護海域,魚類群落粒徑結(jié)構(gòu)較平穩(wěn),生物粒徑譜斜率較平緩;未受保護海域,魚類群落粒徑結(jié)構(gòu)急劇變化,生物粒徑譜斜率較陡峭。

    2.3.3 粒徑分選法

    1983年,SHIRAYAMA[31]最早將此方法用于海洋底棲生物粒徑結(jié)構(gòu)的研究,之后逐步應(yīng)用于魚類的研究。該方法主要是根據(jù)魚類大小來選取網(wǎng)目直徑,獲得不同粒級范圍內(nèi)魚類物種的生物量與豐度狀況,是研究魚類粒徑譜最常用的方法。

    2.3.4 異速生長 P∶B比率法

    異速生長P∶B比率法即魚類生產(chǎn)量與生物量之間比率(production:biomass,P∶B),最早由BANSE等[32]提出被廣泛應(yīng)用于評估魚類群落的生產(chǎn)力。JUNG等[33]于1995~2000年對切薩皮克海灣魚類群落研究,用P∶B比率來評估魚類群落的生產(chǎn)力狀況,并提出公式:P(w)=β(w)×100.44-0.26lg(w)。其中,P(w)代表粒級 w的年際生產(chǎn)量;β(w)代表粒級w的年均現(xiàn)存生物量;w單位是kcal,1 g濕重=1 269 kcal;年際總魚類生產(chǎn)量 P=ΣP(w)。

    2.3.5 PPMR和 TE

    平均捕食者規(guī)格與被捕食者規(guī)格比值(the ratio of the mean predator size to mean prey size,PPMR)和營養(yǎng)轉(zhuǎn)化效率(trophic transfer efficiency,TE)對于魚類群落豐度-質(zhì)量關(guān)系粒徑譜的斜率有顯著的約束性,且二者在不同的生態(tài)系統(tǒng)中具有顯著的一致性[34]。CYR[35]和BROWN等[36]指出,在一個完整的粒徑譜中,PPMR和TE可以預(yù)測不同粒徑魚類物種的能量變化和斜率變化。JENNINGS等[37]為評估捕撈對北海魚類群落的影響,計算PPMR和TE來預(yù)測未開發(fā)海域魚類生物粒徑譜的斜率,與現(xiàn)有數(shù)據(jù)估算的豐度-質(zhì)量關(guān)系進行比較,得出未開發(fā)海域到開發(fā)海域的魚類生物量隨生物體質(zhì)量的縮放范圍為 M-0.1~M-1.0(M:body mass)。

    2.3.6 生化分析法

    用生物、化學(xué)原理和方法來研究魚類粒徑譜。主要包括化合物分析、元素分析(常量元素、微量元素[38]、同位素)、C/N比值[39](主要是脂類與蛋白質(zhì)的組成成分)和能量分析等方法。CHANDLER等[40]于1998年利用穩(wěn)定同位素方法(δ15N)對美國安大略湖東北部湖泊魚類群落研究表明:魚類群落營養(yǎng)結(jié)構(gòu)呈現(xiàn)“金字塔”模式[41]。

    2.4 影響魚類粒徑譜的主要因素

    2.4.1 生產(chǎn)力

    魚類正態(tài)化生物粒徑譜的斜率與海域生產(chǎn)力相關(guān)[42]。生產(chǎn)力高的海域(上升流中心),最小粒級的豐度非常高,魚類群落生物量隨著粒徑的增大而明顯降低。生產(chǎn)力低的海域,魚類群落生物量隨著粒徑的增加而增加或保持為一個不變的常數(shù),且魚類粒徑譜的斜率往往較低[42]。EMMRICH等[43]于2011年對德國北部湖泊魚類群落研究表明,生產(chǎn)力低的湖泊,大型食肉魚類的豐度較高,魚類生物粒徑譜的斜率在水體深處比較平緩。

    2.4.2 棲息地環(huán)境

    棲息地環(huán)境影響魚類的生長,導(dǎo)致魚類生物粒徑譜上各粒級范圍內(nèi)魚類物種結(jié)構(gòu)的不同。MACPHERSON等[30]于1997~1998年對地中海西北部沿海區(qū)域魚類群落研究表明,巖石基底海域,魚類群落生物量隨著粒徑的增大而保持為一個不變的常數(shù)。海草床基底海域,保護區(qū)魚類群落生物量隨著粒徑的增大而增大;未保護區(qū)魚類群落生物量隨著粒徑的增大而保持為一個不變的常數(shù)。

    2.4.3 捕撈強度

    捕撈強度影響魚類的粒徑結(jié)構(gòu)與豐度[44],導(dǎo)致魚類群落平均粒徑下降和大型魚類比例減少,且魚類生物粒徑譜的斜率隨著捕撈強度的增大而更加陡峭。MACPHERSON等[30]對地中海西北部魚類群落研究表明,未受捕撈海域,魚類群落粒徑結(jié)構(gòu)較平穩(wěn)且粒徑譜斜率較平緩;受捕撈海域,魚類群落粒徑結(jié)構(gòu)急劇變化且粒徑譜斜率較陡峭。ZWANENBURG[45]也表明,與捕撈強度低的魚類群落相比,捕撈強度高的魚類群落粒徑譜的斜率往往比較陡峭。

    2.4.4 競爭

    在生態(tài)系統(tǒng)食物網(wǎng)中,生物之間存在捕食、競爭、捕食和競爭三種關(guān)系。競爭影響魚類群落的粒徑結(jié)構(gòu)。ROSEN[21]研究表明,正是由于生態(tài)系統(tǒng)中存在競爭關(guān)系,群落中大型物種生物量減少促使小型物種生物量的增殖,從而導(dǎo)致魚類生物粒徑譜中任意粒級范圍內(nèi)不同粒徑物種比例的差異。

    2.4.5 營養(yǎng)級水平

    強大的粒徑結(jié)構(gòu)在海洋生態(tài)系統(tǒng)中是普遍存在的[7],且營養(yǎng)水平與生物體粒徑大小有關(guān)[46]。在大多數(shù)水生生態(tài)群落中,營養(yǎng)級水平隨著生物體粒徑的增大而不斷地上升[47]。PAHL[48]于1993年對生態(tài)系統(tǒng)食物網(wǎng)時空格局研究表明,食物網(wǎng)中營養(yǎng)級水平不僅可以比較不同生態(tài)系統(tǒng)的結(jié)構(gòu)特征,而且可以作為預(yù)測食物網(wǎng)內(nèi)部動態(tài)結(jié)構(gòu)的指標(biāo)。JENNINGS等[47]于2002年對英國北海魚類群落研究表明,在不考慮物種差異前提下,魚類生物體粒徑大小與其對應(yīng)的營養(yǎng)級水平呈顯著性相關(guān);在相似的魚類群落中,營養(yǎng)級水平與生物體質(zhì)量呈線性關(guān)系,且這種線性關(guān)系可用來預(yù)測魚類群落的營養(yǎng)級結(jié)構(gòu)特征。

    2.4.6 漁具選擇

    漁業(yè)具有選擇性,往往體型大的、營養(yǎng)價值高的魚類物種是消費者所熱衷的,因此需要根據(jù)所需魚類物種來選擇漁具,常見的主要漁具有底拖網(wǎng)、流刺網(wǎng)和張網(wǎng)[49]等,同時漁具的物理影響會導(dǎo)致魚類選擇性死亡。JENNINGS等[47]于2002年對英國北海魚類群落研究表明,當(dāng)拖網(wǎng)拖過海床時,脆弱的大型魚類物種可能被捕獲,而小型魚類物種可能被漁具前面的壓力波推到一邊,從而導(dǎo)致拖網(wǎng)海域底層的群落結(jié)構(gòu)由小型魚類物種組成,且粒徑譜的斜率在拖網(wǎng)嚴(yán)重的海域較大。

    2.4.7 季節(jié)變化

    季節(jié)變化是影響魚類群落粒徑結(jié)構(gòu)的主要因子,不同的季節(jié),魚類群落的粒徑結(jié)構(gòu)組成不同。當(dāng)魚類物種由小型向大型轉(zhuǎn)變時、魚類群落動態(tài)變化跨越大的湖泊和海洋時,可以通過季節(jié)變化來研究給定魚類群落粒徑結(jié)構(gòu)的時空差異。SPRULES等[9]于1991年對美國密歇根湖魚類群落研究表明,魚類群落粒徑結(jié)構(gòu)在季節(jié)上存在顯著的差異性。

    2.4.8 補償比率

    種群間魚類物種的補償比率影響魚類群落的粒徑結(jié)構(gòu),且年際補償中大的波動影響魚類的密度,導(dǎo)致魚類群落生物量與粒徑結(jié)構(gòu)的變化[50]。補償比率大,魚類群落的粒徑結(jié)構(gòu)就相對完整;補償比率小,魚類群落的粒徑結(jié)構(gòu)就相對混亂。JENNINGS等[50]于2007年對海洋生物粒徑開發(fā)和保護的研究中指出,魚類群落密度依賴于調(diào)控,產(chǎn)卵-補充關(guān)系中現(xiàn)存的補償比率對魚類生物粒徑譜的斜率起著關(guān)鍵的機制作用。

    2.5 魚類粒徑譜的應(yīng)用

    2.5.1 環(huán)境擾動監(jiān)測

    通過魚類正態(tài)化粒徑譜斜率可以捕獲到人為和自然干擾破壞生態(tài)的過程。DUPLISEA等[51]用魚類正態(tài)化粒徑譜參數(shù)來說明自然擾動對加拿大司考田沙洲深層魚類的影響。BLANCHARD等[52]以北海4個魚類種群為例建立粒徑譜模型來評估種群和群落對捕撈的響應(yīng),進行環(huán)境監(jiān)測。GRAHAM等[53]用魚類粒徑譜來研究捕撈對珊瑚礁魚類的影響,得出結(jié)論,魚類粒徑譜的斜率隨著捕撈強度的增加而逐步減小,且在高捕撈強度海域斜率變得顯著增加;魚類粒徑譜的中點高度(粒徑間隔中點值所對應(yīng)的y軸生物量值)隨著捕撈強度的增加而顯著地下降。在粒級范圍內(nèi),魚類群落的豐度與捕撈強度呈負(fù)相關(guān)性。在捕撈海域,大型魚類豐度的下降速度比小型魚類豐度的下降速度更快(圖2、圖3)。同時魚類粒徑譜可以用來監(jiān)測水體的營養(yǎng)狀態(tài)。如:SPRULES等[28]研究表明,在富營養(yǎng)化生態(tài)系統(tǒng)中,魚類生物粒徑譜的斜率較小、截距較高或曲率較大、峰值較高。EMMRICH等[43]對德國北部湖泊魚類群落研究表明:在生產(chǎn)力較低的湖泊,大型食肉魚類豐度較高,魚類生物粒徑譜的斜率在水體深處較小。

    2.5.2 反映魚類群落特征

    魚類粒徑譜曲線上穹頂與波谷所形成的“峰型”能夠反映魚類群落的結(jié)構(gòu)特征。JUNG等[33]對切薩皮克海灣魚類群落研究得出,魚類生物粒徑譜呈現(xiàn)“雙峰”型,且第一峰穹頂對應(yīng)粒徑物種結(jié)構(gòu)由小型食浮游動物的魚類組成,主要是鳀魚(Engraulis japonicus);第二峰穹頂對應(yīng)粒徑物種結(jié)構(gòu)由大型雜食性的魚類組成。季節(jié)上,第一峰穹頂對應(yīng)的魚類生物量水平和生物個體平均粒徑從4月到10月呈逐步增加趨勢;第二峰穹頂對應(yīng)的魚類生物量水平相對穩(wěn)定。地域上,上、下海灣的第二峰穹頂?shù)纳锪克礁哂诘谝环羼讽?,但?dāng)季節(jié)性和間接性缺氧同時發(fā)生時,中海灣區(qū)的第一、二峰穹頂?shù)纳锪克阶钚。▓D3)。同時,在水生生態(tài)系統(tǒng)中,魚類粒徑譜能夠反映魚類群落的生態(tài)過程、生物量分布和能量流動[54]。重要的生態(tài)作用和生理過程都依賴于生物體粒徑的大小,如:死亡率、轉(zhuǎn)換效率、新陳代謝和 生 長 等[24]。PETERSON等[55]在 SILVERT等[18]基礎(chǔ)上以粒徑為基礎(chǔ)建立粒徑模型來估算魚類群落的死亡率。

    圖2 捕撈強度與粒徑譜斜率關(guān)系[53]Fig.2 Relationship between slope and fishing intensity[53]

    圖3 捕撈強度與中點高度關(guān)系[53]Fig.3 Relationship between mid-point height and fishing intensity[53]

    2.5.3 預(yù)測魚類生物量和生產(chǎn)量

    BORGMANN[3]和 THIEBAUX等[16-17]用捕食者-被捕食者異速生長關(guān)系建立模型來評估不同營養(yǎng)組的生物量和生產(chǎn)量。BORGMANN認(rèn)為,如果一個粒徑組的生產(chǎn)量或生物量、粒徑轉(zhuǎn)換效率ε、以粒徑為依賴的指數(shù)流動率n是已知的,那么任意粒徑組生物體的粒徑范圍、生物量和生產(chǎn)量都可以被預(yù)測。SPRULES等[9]于1991年用已公布的ε和n值來預(yù)測美國密歇根湖浮游植物、浮游動物和魚類的生物量,且與實際觀測密切符合。THYGESEN等[56]用水生生物群落豐度粒徑譜模型來計算給定粒徑的魚類的生物量大小。

    圖4 1995~2000年切薩皮克海灣中層拖網(wǎng)收集的魚類平均正態(tài)化生物量粒徑譜,由第一穹頂?shù)膬?yōu)勢種和第二穹頂物種雜食性魚類組成(虛線為線性項)[33]Fig.4 Mean normalized biomass size spectrum for fish species,the dominant species in the first dome and the piscivorous fishes component of the second dome collected by mid-water trawl in Chesapeake Bay from 1995 to 2000(dotted line is linear term)[33]

    2.5.4 反映生態(tài)系統(tǒng)的結(jié)構(gòu)和功能

    魚類正態(tài)化粒徑譜參數(shù)能夠反映生態(tài)系統(tǒng)的結(jié)構(gòu)和功能[18]。MACPHERSON等[30]對地中海西北部魚類群落研究表明,受保護海域,斜率均大于-1,則魚類生物量隨著粒徑的增大而增大;大的破碎巖石基底海域的截距大于海藻床基底海域的截距,則其魚類群落生產(chǎn)力水平高于海藻床基底海域的生產(chǎn)力水平。未保護海域,輕微傾斜巖石基底海域的斜率小于-1,則魚類生物量隨著粒徑的增大而減??;海藻床基底海域的斜率大于-1,則魚類生物量隨著粒徑的增大而增大;輕微傾斜巖石基底海域的截距大于海藻床基底海域的截距,則其魚類群落的生產(chǎn)力水平高于海藻床基底海域的生產(chǎn)力水平。R2值不同,則偏離群落穩(wěn)定狀態(tài)的程度有所差異,如表1。

    表1 地中海西北部沿海保護海域與未保護海域魚類群落NBSS參數(shù)值[30]Tab.1 Statistics for normalized biomass size spectra by area(protected and unprotected areas)of littoral fishes in the north western Mediterranean[30]

    3 研究展望

    自從粒徑譜理論被引入水生生態(tài)系統(tǒng)研究以來,我國對于浮游生物、海洋底棲生物粒徑譜的研究較為廣泛,而對于魚類粒徑譜的研究甚少,幾乎為空白。粒徑譜理論引入為研究魚類乃至海洋生態(tài)學(xué)提供了新的方法。與經(jīng)典的依靠分類階元進行魚類生物群落特征的描述相比,粒徑譜方法可以不依賴于分類學(xué)知識,避免因分類鑒定不準(zhǔn)確而引起的計算誤差,可以更好地增強對海洋生態(tài)系統(tǒng)群落特征以及魚類物種生態(tài)和生理過程的認(rèn)識,同時預(yù)測、評估魚類群落的生產(chǎn)量、生產(chǎn)力、分布以及對環(huán)境進行監(jiān)測。

    當(dāng)前大多數(shù)魚類粒徑譜理論的研究均是基于海上調(diào)查的生態(tài)過程,且許多理論預(yù)測的線性正態(tài)化生物量粒徑譜的斜率為-1,同時也產(chǎn)生了一些新的問題需要我們今后去解決[57]:第一,與淡水生物量粒徑譜相比,為什么海洋生物量粒徑譜更多呈現(xiàn)線性關(guān)系;第二,普遍觀察到的正態(tài)化生物量粒徑譜的穹頂是否是穩(wěn)定的;第三,在受到外界干擾恢復(fù)之后,魚類群落能否有一個平衡狀態(tài)的粒徑結(jié)構(gòu)[58]。

    魚類粒徑譜的研究極大依賴于魚類群落生物個體的粒徑大小和所抽取的樣本大小。研究樣品的獲得以及測量是目前學(xué)術(shù)界對魚類粒徑譜的研究存在差異性的主要原因。首先,不同的研究者所采取的研究方法、條件有所不同,必然造成結(jié)果出現(xiàn)差異;其次,多數(shù)研究者均是人工采樣及測量數(shù)據(jù),避免不了出現(xiàn)誤差,不能真實反映魚類群落特征。比如:測量與換算誤差、劃分粒徑依據(jù)、樣本大小、儀器設(shè)備新舊、研究區(qū)域差異、魚類生活習(xí)性(產(chǎn)卵期、捕食期、競爭)以及季節(jié)差異等因素,都會影響研究結(jié)果;最后,因為研究對象和地域的時空差異,先前多數(shù)學(xué)者的研究成果不能很好地進行比較研究,給后人研究帶來一定的困難。因此,在抽樣調(diào)查時應(yīng)抽取盡可能大且粒徑范圍廣泛的生物體樣本,確保所有粒徑間隔上的生物體豐度都能被嚴(yán)格的估算到,并加大對新采樣技術(shù)的研發(fā)與應(yīng)用,如水聲學(xué)和光學(xué)等,以取代費時費力的拖網(wǎng)采樣,以便更好地反映整個魚類群落的粒徑結(jié)構(gòu)和生態(tài)特征。同時,研發(fā)通用的基于粒徑譜的魚類生物量評估理論模型,并規(guī)范衡量標(biāo)準(zhǔn),應(yīng)該是未來魚類粒徑譜研究的趨勢之一。

    [1] 沈國英.海洋生態(tài)學(xué)(第3版)[M].北京:科學(xué)出版社,2010.

    SHEN G Y.Marine ecology(The third edition)[M].Beijing:Science press,2010.

    [2] BENOIT E,ROCHET M J.A continuous model of biomass size spectra governed by predation and the effects of fishing on them[J].Journal of Theoretical Biology,2004,226(1):9-21.

    [3] BORGMANN U.Particle-size-conversion efficiency and total animal production in pelagic ecosystems[J].Canadian Journal of Fisheries and Aquatic Sciences,2011,39(5):668-674.

    [4] BOUDREAU PR,DICKIE L M.Biomass spectra of aquatic ecosystems in relation to fisheries yield[J].Canadian Journal of Fisheries and Aquatic Sciences,2011,49(8):1528-1538.

    [5] CYR H,PETERS R H.Biomass-size spectra and the prediction of fish biomass in lakes[J].Canadian Journal of Fisheries and Aquatic Sciences,2011,53(5):994-1006.

    [6] LEACH J H,DICKIE L M,SHUTER B J,et al.A review of methods for prediction of potential fish production with application to the Great Lakes and Lake Winnipeg[J].Canadian Journal of Fisheries and Aquatic Sciences,2011,44(S2):471-485.

    [7] SHELDON R W,PRAKASH A,SUTCLIFFE W H.The size distribution of particles in the ocean[J].Limnology and Oceanography,1972,17(3):327-340.

    [8] SHELDON R W,SUTCLIFFE W H,PARANJAPE M A.Structure of pelagic food chain and relationship between plankton and fish production[J].Journal of the Fisheries Research Board of Canada,1977,34(12):2344-2353.

    [9] SPRULESW G,BRANDT SB,STEWART D J,et al.Biomass size spectrum of the Lake Michigan Pelagic food web[J].Canadian Journal of Fisheries and Aquatic Sciences,1991,48(1):105-115.

    [10] BORGMANN U,WHITTLE D M.Particle-sizeconversion efficiency and contaminant concentrations in Lake Ontario Biota[J].Canadian Journal of Fisheries and Aquatic Sciences,2011,40(3):328-336.

    [11] COTTINGHAM K L.Nutrients and zooplankton as multiple stressors of phytoplankton communities:Evidence from size structure[J].Limnology and Oceanography,1999,44(3):810-827.

    [12] SHELDON R W,PARSONST R.A continuous size spectrum for particulate matter in the sea[J].Journal of the Fisheries Research Board of Canada,2011,24(5):909-915.

    [13] 周林濱,譚燁輝,黃良民,等.水生生物粒徑譜/生物量譜研究進展[J].生態(tài)學(xué)報,2010,30(12):3319-3333.

    ZHOU L B,TAN Y H,HUANG L M,et al.The advances in the aquatic particle/biomass size spectra study[J].Acta Ecologica Sinica,2010,30(12):3319-3333.

    [14] 王睿照,張志南.海洋底棲生物粒徑譜的研究[J].海洋湖沼通報,2003(4):61-68.WANG R Z,ZHANG Z N.A study on the size spectra of benthos[J].Transactions of Oceanology and Limnology,2003(4):61-68.

    [15] KERR SR.Theory of size distribution in ecological communities[J].Journal of the Fisheries Research Board of Canada,2011,31(12):1859-1862.

    [16] THIEBAUX M L,DICKIE L M.Models of aquatic biomass size spectra and the common structure of their solutions[J].Journal of Theoretical Biology,1992,159(2):147-161.

    [17] THIEBAUX M L,DICKIE L M.Structure of the body-size spectrum of the biomass in aquatic ecosystems:A consequence of allometry in predatorprey interactions[J].Canadian Journal of Fisheries and Aquatic Sciences,1993,50(6):1308-1317.

    [18] SILVERT W,PLATT T.Energy flux in the pelagic ecosystem: A time-dependent equation[J].Limnology and Oceanography,1978,23(4):813-816.

    [19] BOUDREAU PR,DICKIE L M,KERR SR.Bodysize spectra of production and biomass as systemlevel indicators of ecological dynamics[J].Journal of Theoretical Biology,1991,152(3):329-339.

    [20] ZHOU M,HUNTLEY M E.Population dynamics theory of plankton based on biomass spectra[J].Marine Ecology Progress Series,1997,159(159):61-73.

    [21] ROSEN C J.Modelling an exploited marine fish community with 15 parameters-results from a simple size-based model[J]. Ices Journal of Marine Science,2006(63):1029-1044.

    [22] LAWR,PLANK M J,JAMESA,et al.Size-spectra dynamics from stochastic predation and growth of individuals[J].Ecology,2009,90(3):802-811.

    [23] ROSSBERG A G.6-a complete analytic theory for structure and dynamics of populations and communities spanning wide ranges in body size[J].Advances in Ecological Research,2012(46):427-521.

    [24] KERR SR,DICKIE L M.The biomass spectrum:A predator-prey theory of aqautic production[M].Japan:Columbia University Press,2001.

    [25] PLATT T,DENMAN K.The structure of pelagic marine ecosystems[J].Rapports Et Proces Verbaux Des Reunions,2010(173):60-65.

    [26] 鄧 可,張志南,黃 勇,等.南黃海典型站位底棲動物粒徑譜及其應(yīng)用[J].中國海洋大學(xué)學(xué)報(自然科學(xué)版),2005,35(6):1005-1010.DENG K, ZHANG Z N, HUANG Y, et al.

    Application of benthic biomass size spectra of macroand microfuana at typical stations in the Southern Yellow Sea, China[J]. Periodical of Ocean University of China(Natural Science Edition),2005,35(6):1005-1010.

    [27] WITEK Z,KRAJEWSKASOLTYSA.Some examples of the epipelagic plankton size structure in high latitude oceans[J].Journal of Plankton Research,1989,11(6):1143-1155.

    [28] SPRULESW,MUNAWAR M.Plankton size spectra in relation to ecosystem productivity, size, and perturbation[J].Canadian Journal of Fisheries and Aquatic Sciences,1986,43(9):1789-1794.

    [29] PLATT T,DENMAN K.Organisation in the pelagic ecosystem [J]. Helgol?nder Wissenschaftliche Meeresunter suchungen,1977,30(1-4):575-581.

    [30] MACPHERSON E,GORDOA A,GARCIA-RUBIES A.Biomass size spectra in littoral fishes in protected and unprotected areas in the NWMediterranean[J].Estuarine Coastal and Shelf Science,2002,55(5):777-788.

    [31] SHIRAYAMA Y.Ecology of deep-sea meiobenthos in the western Pacific[J]. Journal of the Oceanographical Society of Japan,1989,45(45):83-93.

    [32] BANSE K,MOSHER S.Adult body mass and annual production/biomass relationships of field populations[J].Ecological Monographs,1980,50(3):355-379.

    [33] JUNG S,HOUDE E D.Fish biomass size spectra in Chesapeake Bay[J].Estuaries and Coasts,2005,28(2):226-240.

    [34] JENNINGS S.Size-based analysis of aquatic food webs[J].Oxford University Press,2005:86-97.

    [35] CYR H.Individual energy use and the allometry of population density:Scaling in biology[M].New York:Oxford University Press,2000.

    [36] BROWN JH,GILLOOLY JF.Ecological food webs:High-quality data facilitate theoretical unification[J]. Proceedings of the National Academy of Sciences,2003,100(4):1467-1468.

    [37] JENNINGS S,BLANCHARD J L.Fish abundance with no fishing:predictions based on macroecological theory[J].Journal of Animal Ecology,2004,73(4):632-642.

    [38] 劉 云,王際英,李寶山,等.海水魚類微量元素需求研究進展[J].海洋漁業(yè),2015,37(4):378-385.

    LIU Y,WANG J Y,LI B S,et al.Research progress on requirements of marine fishes for trace elements[J].Marine Fisheries,2015,37(4):378-385.

    [39] 趙書燕,林黑著,黃 忠,等.不同蛋白質(zhì)水平下添加小肽對石斑魚生長、消化酶、血清生化和抗氧化能力的影響[J].南方水產(chǎn)科學(xué),2016,12(3):15-23.

    ZHAO SY,LIN H Z,HUANG Z,et al.Effect of small peptide supplementation at different protein levels on growth performance,digestive enzymes activities,serum biochemical indices and antioxidant abilities of grouper(Epinephelus akaara)[J].South China Fisheries Science,2016,12(3):15-23.

    [40] CHANDLER M,PETERS R.Mapping trophic continua of benthic foodwebs: Body size-δ15N relationships[J].Marine Ecology Progress Series,1998(174):301-306.

    [41] 紀(jì)煒煒,李圣法,陳雪忠.魚類營養(yǎng)級在海洋生態(tài)系統(tǒng)研究中的應(yīng)用[J].中國水產(chǎn)科學(xué),2010,17(4):878-887.

    JI W W,LI S F,CHEN X Z.Application of fish trophic level in marine ecosystem[J].Journal of Fishery Sciences of China,2010,17(4):878-887.

    [42] MACPHERSON E,GORDOA A.Biomass spectra in benthic fish assemblages in the Benguela System[J].Marine Ecology Progress Series,1996,138(1-3):27-32.

    [43] EMMRICH M,BRUCET S,RITTERBUSCH D,et al.Size spectra of lake fish assemblages:Responses along gradients of general environmental factors and intensity of lake-use[J].Freshwater Biology,2011,56(11):2316-2333.

    [44] 王躍中,孫典榮,賈曉平,等.捕撈壓力和氣候變化對東海馬面鲀漁獲量的影響[J].南方水產(chǎn)科學(xué),2013,9(1):8-15.

    WANG Y Z,SUN D R,JIA X P,et al.Influence of fishing pressure and climate change on filefish catches in East China Sea[J].South China Fisheries Science,2013,9(1):8-15.

    [45] ZWANENBURG K C T.The effects of fishing on demersal fish communities of the Scotian Shelf[J].Ices Journal of Marine Science,2000,57(3):503-509.

    [46] FRY B,QUINONESR.Biomass spectra and stable isotope indicators of trophic level in zooplankton of the northwest Atlantic[J].Marine Ecology Progress Series,1994,112(1-2):201-204.

    [47] JENNINGS S,PINNEGAR J K,POLUNIN N V C,et al.Linking size-based and trophic analyses of benthic community structure[J].Marine Ecology Progress Series,2002,226(1):77-85.

    [48] PAHL WC.Food webs and ecological networks across temporal and spatial scales[J].Oikos,1993,67(3):415-432.

    [49] 晏 磊,譚永光,楊炳忠,等.基于張網(wǎng)漁業(yè)休漁前后的黃茅海河口漁業(yè)資源群落比較[J].南方水產(chǎn)科學(xué),2016,12(6):1-8.

    YAN L,TAN Y G,YANG B Z,et al.Comparison on resources community of stow-net fishery before and after fishing off season in Huangmaohai Estuary[J].South China Fisheries Science,2016,12(6):1-8.

    [50] JENNINGS S,REYNOLDS J D.Body size,exploitation and conservation of marine organisms[M]//Body size:The Structure and Function of Aquatic Ecosystems.England:Cambridge University Press,2007.

    [51] DUPLISEA D E,KERR S R.Application of a biomass size spectrum model to demersal fish data from the Scotian Shelf[J].Journal of Theoretical Biology,1995,177(3):263-269.

    [52] BLANCHARD J L,ANDERSEN K H,SCOTT F,et al.Evaluating targets and trade-offs among fisheries and conservation objectives using a multispecies size spectrum model[J].Journal of Applied Ecology,2014,51(3):612-622.

    [53] GRAHAM N,DULVY N K,JENNINGS S,et al.Size-spectra as indicators of the effects of fishing on coral reef fish assemblages[J].Coral Reefs,2005,24(1):118-124.

    [54] TREBILCO R,BAUM J K,SALOMON A K,et al.Ecosystem ecology:Size-based constraints on the pyramids of life[J]. Trends in Ecology and Evolution,2013,28(7):423-431.

    [55] PETERSON I,WROBLEWSKI J S.Mortality rate of fishes in the pelagic ecosystem[J]. Canadian Journal of Fisheries and Aquatic Sciences,2011,41(7):1117-1120.

    [56] THYGESEN U H,F(xiàn)ARNSWORTH K D,ANDERSEN K H,et al.How optimal life history changes with the community size-spectrum[J].Proceedings Biological Sciences,2005,272(1570):1323-1331.

    [57] SPRULES W G,BARTH L E.Surfing the biomass size spectrum:Some remarks on history,theory,and application[J].Canadian Journal of Fisheries and Aquatic Sciences,2015,12(9):3459-3466.

    [58] ANDERSEN K H,JACOBSEN NS,F(xiàn)ARNSWORTH K D.The theoretical foundations for size spectrum models of fish communities[J].Canadian Journal of Fisheries and Aquatic Sciences,2016,73(4):575-588.

    Advances in fish particle size spectra study

    GUO Jian-zhong1,2,CHEN Zuo-zhi1,XU Shan-nan1
    (1.Key Laboratory of South China Sea Fishery Resources Exploitation&Utilization of Ministry of Agriculture;Key Laboratory of Fishery Ecology and Environment of Guangdong Province,South China Sea Fisheries Research Institute,Chinese Academy of Fishery Sciences,Guangzhou 510300,China;2.College of Fisheries and Life Science,Shanghai Ocean University,Shanghai 201306,China)

    Particle size spectrum is a kind of curve graph that represents the relationship between biomass and particle size.It can reflect the structure and function of the ecosystem and its dynamic connection inside the system,which is currently becoming a hotspot issue of the aquatic ecology.At present,the research on fish particle size spectrum in our country is quite little,but there are more researches in this field and more advanced theory and methods in foreign country,which can provide worthy reference for domestic scholars.In this paper,a large number of fish particle size spectrum research literatures and the development of recent years were summarized.The paper introduced theoretical principles,research methods,influence factors and application of particle size spectrum,which referred to fish production and biomass assessment,ecological monitoring,reflection on the characteristics of fish community and ecosystem structure and function,and the prediction of human activities and natural disturbance impact on marine ecosystem.Fishing intensity,marine productivity level,the level of nutrition,fish community habitat environment,fish community species competition and compensation efficiency,seasonal changes and the types of fishing gears and meshes are main factors affecting the fish particle size spectrum.Most of the current fish particle size spectrum theory studies were based on the investigation of marine ecological processes,and many theoretical predictions of the linear slope of normalized biomass particle size spectra as-1.The slope,intercept and regression coefficient respectively affected the characteristics of fish community.The slope represented the relationship between the biomass and the size.The intercept represented the magnitude of productivity level.The regression coefficient indicated the degree of community deviation from the steady state.It also brought some new problems needed to solve in the future:firstly,compared with fresh biomass particle size spectra,why did the marine biomass particle size spectra show more linear relationships?Secondly,was the generally observed normalized biomass particle size spectrum of dome is stable?Thirdly,after the resumption from the interference by environment,did the fish community have a balanced state of the particle size structure?Based on the deficiencies and problems in current researches,this paper proposes the application of new observation and sampling technology,and the generally theoretical model research based on fish particle size spectrum.

    particle size spectrum;fish particle size spectrum;ecosystem

    S 931

    A

    1004-2490(2017)05-0582-10

    2016-10-08

    國家重點基礎(chǔ)研究發(fā)展計劃(2015CB452904);廣東省自然科學(xué)基金項目(2015A030313785);國家自然科學(xué)基金項目(31100362);廣州市科技計劃項目珠江科技新星專項(2014J2200020)

    郭建忠(1990-),男,河南省南召縣人,碩士研究生,主要從事漁業(yè)生態(tài)研究。

    E-mail:1115982121@qq.com

    徐姍楠,副研究員。Tel:020-89108303,E-mail:xushannan@scsfri.ac.cn

    猜你喜歡
    正態(tài)斜率魚類
    基于MFCC和ResNet的魚類行為識別
    物理圖像斜率的變化探討
    物理之友(2020年12期)2020-07-16 05:39:16
    奇妙的古代動物 泥盆紀(jì)的魚類
    雙冪變換下正態(tài)線性回歸模型參數(shù)的假設(shè)檢驗
    求斜率型分式的取值范圍
    基于子孔徑斜率離散采樣的波前重構(gòu)
    魚類是怎樣保護自己的
    基于泛正態(tài)阻抗云的諧波發(fā)射水平估計
    MMC-MTDC輸電系統(tǒng)新型直流電壓斜率控制策略
    電測與儀表(2016年6期)2016-04-11 12:05:54
    半?yún)?shù)EV模型二階段估計的漸近正態(tài)性
    久久影院123| 六月丁香七月| 日韩强制内射视频| 亚洲av免费在线观看| 在线a可以看的网站| 亚洲欧美日韩东京热| 老师上课跳d突然被开到最大视频| 国语对白做爰xxxⅹ性视频网站| 在线观看免费高清a一片| 少妇裸体淫交视频免费看高清| 国产亚洲午夜精品一区二区久久 | 热99国产精品久久久久久7| 卡戴珊不雅视频在线播放| 国产成人免费无遮挡视频| 亚洲成人中文字幕在线播放| 久久久久久久久久久免费av| 最近中文字幕2019免费版| 1000部很黄的大片| 日韩人妻高清精品专区| 少妇被粗大猛烈的视频| 国产真实伦视频高清在线观看| 99热国产这里只有精品6| 一级毛片 在线播放| 亚洲内射少妇av| 色视频在线一区二区三区| 一区二区三区乱码不卡18| 色婷婷久久久亚洲欧美| 久久99热这里只频精品6学生| 国产亚洲午夜精品一区二区久久 | 免费看不卡的av| 久久久久久久久久成人| 麻豆成人av视频| 国产精品.久久久| 国产精品.久久久| 亚洲精品一二三| 国产午夜精品一二区理论片| 高清视频免费观看一区二区| 亚洲成色77777| 中文字幕av成人在线电影| 免费观看a级毛片全部| 我的女老师完整版在线观看| 国产黄a三级三级三级人| 观看美女的网站| 亚洲国产最新在线播放| 久久久久精品性色| 超碰97精品在线观看| 男女那种视频在线观看| 亚洲无线观看免费| 精品国产露脸久久av麻豆| 天天躁夜夜躁狠狠久久av| 国内揄拍国产精品人妻在线| 免费人成在线观看视频色| 青春草亚洲视频在线观看| 日本一二三区视频观看| 国产黄色视频一区二区在线观看| 国产老妇伦熟女老妇高清| 如何舔出高潮| 国产综合懂色| 尾随美女入室| 你懂的网址亚洲精品在线观看| 五月伊人婷婷丁香| 久久久成人免费电影| 亚洲av电影在线观看一区二区三区 | 99久久中文字幕三级久久日本| 激情 狠狠 欧美| 一本一本综合久久| 久久久久久久国产电影| 精品人妻一区二区三区麻豆| 亚洲丝袜综合中文字幕| 国产成人精品福利久久| 九九在线视频观看精品| 欧美精品一区二区大全| 欧美日韩在线观看h| 日韩大片免费观看网站| 各种免费的搞黄视频| 久久综合国产亚洲精品| 亚洲四区av| 国产精品久久久久久精品古装| 国产大屁股一区二区在线视频| 亚洲经典国产精华液单| 亚洲自拍偷在线| xxx大片免费视频| 18禁在线无遮挡免费观看视频| 国产精品一区www在线观看| 中文字幕免费在线视频6| 国产伦精品一区二区三区四那| 精品酒店卫生间| 亚洲精品日韩在线中文字幕| 精华霜和精华液先用哪个| 久久精品久久久久久噜噜老黄| 一级毛片电影观看| 亚洲,一卡二卡三卡| 性色av一级| 日韩视频在线欧美| 三级男女做爰猛烈吃奶摸视频| 久久99精品国语久久久| 老师上课跳d突然被开到最大视频| 一级黄片播放器| 99热全是精品| 国产综合懂色| 伦理电影大哥的女人| 汤姆久久久久久久影院中文字幕| 免费av毛片视频| 各种免费的搞黄视频| 18禁在线播放成人免费| 在现免费观看毛片| 超碰97精品在线观看| 男女无遮挡免费网站观看| 国产精品.久久久| 亚洲av中文av极速乱| 亚洲av免费在线观看| 女人十人毛片免费观看3o分钟| 欧美xxxx黑人xx丫x性爽| 少妇人妻精品综合一区二区| 中文字幕人妻熟人妻熟丝袜美| 久久鲁丝午夜福利片| 久久热精品热| 国产亚洲91精品色在线| 国产成人免费无遮挡视频| 欧美xxxx性猛交bbbb| 日本一本二区三区精品| 黄色一级大片看看| 国产精品国产三级专区第一集| 高清日韩中文字幕在线| 少妇猛男粗大的猛烈进出视频 | 日本免费在线观看一区| 少妇人妻一区二区三区视频| 亚洲国产欧美人成| 成年人午夜在线观看视频| 高清在线视频一区二区三区| 麻豆成人av视频| 午夜福利高清视频| 国产成人免费观看mmmm| 偷拍熟女少妇极品色| 国产美女午夜福利| 国产成人a区在线观看| av天堂中文字幕网| 国产 一区精品| 国产男女内射视频| 插阴视频在线观看视频| 高清毛片免费看| 国产毛片a区久久久久| av播播在线观看一区| 久久ye,这里只有精品| 最近2019中文字幕mv第一页| av专区在线播放| 99精国产麻豆久久婷婷| 精品久久久久久久久亚洲| 欧美少妇被猛烈插入视频| 麻豆成人午夜福利视频| 国产欧美日韩精品一区二区| 日本-黄色视频高清免费观看| 亚洲一区二区三区欧美精品 | 高清在线视频一区二区三区| 久久精品国产亚洲av涩爱| 亚洲精品,欧美精品| 舔av片在线| 欧美3d第一页| 成人午夜精彩视频在线观看| av.在线天堂| 最近最新中文字幕免费大全7| 国产成人a∨麻豆精品| 人妻系列 视频| 小蜜桃在线观看免费完整版高清| 亚洲国产高清在线一区二区三| 菩萨蛮人人尽说江南好唐韦庄| 最近中文字幕高清免费大全6| 国产精品一二三区在线看| 久久久欧美国产精品| 少妇高潮的动态图| 国产免费一级a男人的天堂| 国产av不卡久久| 久久精品久久精品一区二区三区| 免费看av在线观看网站| 亚洲性久久影院| 国产 一区精品| 丝袜美腿在线中文| 成人二区视频| 尾随美女入室| 日产精品乱码卡一卡2卡三| 欧美日韩视频高清一区二区三区二| 美女视频免费永久观看网站| 人妻制服诱惑在线中文字幕| 免费观看a级毛片全部| 久久久久久伊人网av| 国产伦理片在线播放av一区| 午夜免费观看性视频| 人体艺术视频欧美日本| 在线观看一区二区三区| 好男人视频免费观看在线| 精品亚洲乱码少妇综合久久| 国产成人精品福利久久| eeuss影院久久| 一区二区三区四区激情视频| 久久国内精品自在自线图片| 国产成人精品一,二区| 中文天堂在线官网| 欧美性猛交╳xxx乱大交人| 香蕉精品网在线| 又爽又黄a免费视频| 精品熟女少妇av免费看| 又粗又硬又长又爽又黄的视频| 一级黄片播放器| 日韩国内少妇激情av| 九九久久精品国产亚洲av麻豆| 一级二级三级毛片免费看| 久热久热在线精品观看| 精品国产一区二区三区久久久樱花 | xxx大片免费视频| 一级毛片电影观看| 欧美xxxx黑人xx丫x性爽| 婷婷色av中文字幕| 日韩成人伦理影院| av天堂中文字幕网| 我的老师免费观看完整版| 欧美日韩精品成人综合77777| 亚洲美女视频黄频| 国产免费视频播放在线视频| 国产中年淑女户外野战色| 高清视频免费观看一区二区| 亚洲伊人久久精品综合| 毛片女人毛片| 精品一区二区三卡| 一区二区三区精品91| 色播亚洲综合网| 中文在线观看免费www的网站| 纵有疾风起免费观看全集完整版| av一本久久久久| 国产69精品久久久久777片| a级一级毛片免费在线观看| 国产黄频视频在线观看| 在线观看免费高清a一片| 亚洲av在线观看美女高潮| 国产成人精品久久久久久| 夫妻午夜视频| 麻豆精品久久久久久蜜桃| tube8黄色片| 亚洲成人av在线免费| 免费av毛片视频| 深爱激情五月婷婷| 六月丁香七月| 亚洲精品色激情综合| 狂野欧美激情性xxxx在线观看| 国产亚洲精品久久久com| 伊人久久国产一区二区| 成人亚洲精品一区在线观看 | 男女无遮挡免费网站观看| 三级男女做爰猛烈吃奶摸视频| 国产精品成人在线| 在线观看免费高清a一片| 精品久久久噜噜| 只有这里有精品99| 亚洲国产精品成人综合色| 联通29元200g的流量卡| 91aial.com中文字幕在线观看| 毛片一级片免费看久久久久| 日本一本二区三区精品| 成人毛片60女人毛片免费| 成人国产av品久久久| 人妻一区二区av| 国产精品国产三级专区第一集| 搡老乐熟女国产| 99久国产av精品国产电影| 丝袜喷水一区| 波野结衣二区三区在线| 免费av毛片视频| 国产黄色视频一区二区在线观看| 亚洲精品一二三| 又大又黄又爽视频免费| 亚洲四区av| 亚洲人成网站在线观看播放| 亚洲成人av在线免费| 在线免费观看不下载黄p国产| 18禁动态无遮挡网站| h日本视频在线播放| 亚洲欧洲日产国产| 狂野欧美白嫩少妇大欣赏| 国产免费又黄又爽又色| 男女边摸边吃奶| 国产av国产精品国产| 夜夜爽夜夜爽视频| 男女啪啪激烈高潮av片| 亚洲国产欧美人成| 美女高潮的动态| .国产精品久久| 久久精品综合一区二区三区| 成人二区视频| 三级经典国产精品| 午夜激情久久久久久久| 狂野欧美白嫩少妇大欣赏| 91狼人影院| 国产精品无大码| 国产免费又黄又爽又色| videos熟女内射| 熟女电影av网| av又黄又爽大尺度在线免费看| 99热6这里只有精品| 色网站视频免费| 丰满乱子伦码专区| 成人二区视频| 欧美日韩一区二区视频在线观看视频在线 | 在线看a的网站| 久久99热这里只频精品6学生| videos熟女内射| 亚洲一区二区三区欧美精品 | 91狼人影院| 日韩精品有码人妻一区| 简卡轻食公司| 国产成人精品福利久久| av免费在线看不卡| 色婷婷久久久亚洲欧美| 18+在线观看网站| eeuss影院久久| 久久久久国产精品人妻一区二区| 国产高清不卡午夜福利| 麻豆国产97在线/欧美| 亚洲人成网站在线播| 亚洲av欧美aⅴ国产| 丝袜美腿在线中文| 亚洲人成网站在线观看播放| 免费av毛片视频| 三级国产精品片| 波野结衣二区三区在线| 国产亚洲91精品色在线| 啦啦啦啦在线视频资源| 日本熟妇午夜| 婷婷色综合大香蕉| 性色avwww在线观看| 另类亚洲欧美激情| 亚洲第一区二区三区不卡| 国产男女内射视频| 中文字幕制服av| 国产午夜福利久久久久久| 亚洲精品国产色婷婷电影| 亚洲精品国产av蜜桃| 精品人妻偷拍中文字幕| 舔av片在线| 国产爽快片一区二区三区| 国产免费福利视频在线观看| 日日撸夜夜添| 久久久久久久久久久免费av| 日产精品乱码卡一卡2卡三| 最近中文字幕高清免费大全6| 久久久久久久亚洲中文字幕| 精品人妻偷拍中文字幕| 国产精品女同一区二区软件| 超碰97精品在线观看| 成人亚洲欧美一区二区av| 国产一区亚洲一区在线观看| 国产精品一区www在线观看| 国产日韩欧美亚洲二区| 日本免费在线观看一区| 国产成人精品婷婷| 国产人妻一区二区三区在| 国产伦理片在线播放av一区| 亚洲精品成人久久久久久| 亚洲国产精品成人综合色| 国产一区有黄有色的免费视频| 国产黄色免费在线视频| 亚洲,一卡二卡三卡| 国产亚洲av片在线观看秒播厂| 天堂中文最新版在线下载 | 国产亚洲91精品色在线| 最近2019中文字幕mv第一页| 如何舔出高潮| 三级男女做爰猛烈吃奶摸视频| 久久精品国产自在天天线| 一本一本综合久久| a级一级毛片免费在线观看| 国产乱人视频| 免费观看无遮挡的男女| 午夜精品国产一区二区电影 | 大又大粗又爽又黄少妇毛片口| 午夜福利在线观看免费完整高清在| 精品国产乱码久久久久久小说| 天天躁日日操中文字幕| 欧美97在线视频| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 久久精品国产亚洲网站| 欧美97在线视频| 国产成人精品久久久久久| 99九九线精品视频在线观看视频| 久久久久性生活片| 人妻制服诱惑在线中文字幕| 日韩三级伦理在线观看| 三级国产精品片| 欧美日本视频| 亚洲美女视频黄频| 亚洲精品影视一区二区三区av| av线在线观看网站| 在线看a的网站| 亚洲精品国产成人久久av| 美女高潮的动态| 你懂的网址亚洲精品在线观看| 久久久精品94久久精品| 亚洲不卡免费看| 只有这里有精品99| 午夜福利视频精品| 久久久久久久亚洲中文字幕| 婷婷色综合大香蕉| 日本欧美国产在线视频| 搞女人的毛片| 久久久久久久午夜电影| 成年av动漫网址| freevideosex欧美| 少妇猛男粗大的猛烈进出视频 | 国产成人免费观看mmmm| 秋霞在线观看毛片| 亚洲经典国产精华液单| 日本猛色少妇xxxxx猛交久久| 最近的中文字幕免费完整| 欧美激情国产日韩精品一区| av免费在线看不卡| 国产久久久一区二区三区| 国产一区二区三区av在线| 亚洲精品成人av观看孕妇| 成年女人看的毛片在线观看| 久久精品久久精品一区二区三区| 久久精品熟女亚洲av麻豆精品| 欧美区成人在线视频| 成人特级av手机在线观看| 成人亚洲精品av一区二区| 国产精品一区二区性色av| 99九九线精品视频在线观看视频| 中文在线观看免费www的网站| 免费看a级黄色片| 精品久久久精品久久久| 看黄色毛片网站| 精品久久久久久久久亚洲| 男插女下体视频免费在线播放| 夜夜爽夜夜爽视频| 日本一本二区三区精品| 国产精品秋霞免费鲁丝片| 日韩一区二区三区影片| 三级国产精品欧美在线观看| 国产免费视频播放在线视频| 午夜免费男女啪啪视频观看| 久久久成人免费电影| 国产精品久久久久久精品古装| 久久久亚洲精品成人影院| 亚洲精品国产av成人精品| 成人免费观看视频高清| 欧美性猛交╳xxx乱大交人| 久久精品熟女亚洲av麻豆精品| 成人国产av品久久久| 日韩制服骚丝袜av| 一区二区三区乱码不卡18| 午夜视频国产福利| 国产一区二区在线观看日韩| 成年人午夜在线观看视频| 一个人看视频在线观看www免费| 色吧在线观看| 亚洲人与动物交配视频| 一级毛片黄色毛片免费观看视频| 91久久精品国产一区二区三区| 夜夜爽夜夜爽视频| 亚洲国产av新网站| 蜜桃亚洲精品一区二区三区| 能在线免费看毛片的网站| 国产淫语在线视频| 我的女老师完整版在线观看| 午夜精品国产一区二区电影 | av一本久久久久| 精品一区二区三卡| 亚州av有码| 欧美成人午夜免费资源| 亚洲成人精品中文字幕电影| 少妇人妻精品综合一区二区| av在线观看视频网站免费| 中文字幕制服av| 97超视频在线观看视频| 18+在线观看网站| 超碰av人人做人人爽久久| 午夜免费男女啪啪视频观看| 国产亚洲午夜精品一区二区久久 | 五月玫瑰六月丁香| 亚洲精品一区蜜桃| 一级毛片aaaaaa免费看小| 欧美最新免费一区二区三区| 久久精品久久精品一区二区三区| 亚洲国产日韩一区二区| 久久久精品欧美日韩精品| 亚洲av中文字字幕乱码综合| 黑人高潮一二区| 成年av动漫网址| 最新中文字幕久久久久| 亚洲天堂国产精品一区在线| 亚洲av免费在线观看| 亚洲丝袜综合中文字幕| 激情 狠狠 欧美| 人人妻人人爽人人添夜夜欢视频 | 亚洲av成人精品一区久久| 成人国产av品久久久| 一边亲一边摸免费视频| 国产成人精品一,二区| 七月丁香在线播放| 老师上课跳d突然被开到最大视频| 色吧在线观看| 在线观看一区二区三区激情| 日韩av免费高清视频| 国产淫语在线视频| 有码 亚洲区| av在线app专区| 99热网站在线观看| 国产爽快片一区二区三区| 哪个播放器可以免费观看大片| 国产一级毛片在线| 爱豆传媒免费全集在线观看| 国产伦理片在线播放av一区| 能在线免费看毛片的网站| 男人舔奶头视频| 亚洲精品乱码久久久v下载方式| 国产男女内射视频| 大香蕉久久网| 午夜爱爱视频在线播放| 日本午夜av视频| 国产色爽女视频免费观看| 嫩草影院精品99| 成人亚洲欧美一区二区av| 久久精品夜色国产| 九九久久精品国产亚洲av麻豆| av专区在线播放| 欧美zozozo另类| 我的老师免费观看完整版| 嘟嘟电影网在线观看| 亚洲精品色激情综合| 亚洲,欧美,日韩| 婷婷色综合www| 国产精品福利在线免费观看| 国产高清三级在线| 亚洲精品国产色婷婷电影| 久久女婷五月综合色啪小说 | 久久久久久久久久久免费av| 国产黄色视频一区二区在线观看| 精品国产三级普通话版| 80岁老熟妇乱子伦牲交| 久久午夜福利片| av天堂中文字幕网| 能在线免费看毛片的网站| 国产有黄有色有爽视频| 亚洲综合精品二区| 你懂的网址亚洲精品在线观看| 一区二区三区四区激情视频| 久久韩国三级中文字幕| 少妇的逼水好多| 狂野欧美激情性xxxx在线观看| 97超碰精品成人国产| 一级毛片电影观看| 亚洲精华国产精华液的使用体验| 国产成人福利小说| 亚洲精华国产精华液的使用体验| 精品少妇久久久久久888优播| 国产精品一区二区在线观看99| 欧美成人一区二区免费高清观看| 看十八女毛片水多多多| 国产成人a∨麻豆精品| av福利片在线观看| 一区二区三区免费毛片| 国产熟女欧美一区二区| 99久国产av精品国产电影| av网站免费在线观看视频| 人妻少妇偷人精品九色| 亚洲精品亚洲一区二区| 99热这里只有是精品在线观看| 黄片无遮挡物在线观看| 在线亚洲精品国产二区图片欧美 | 久久久国产一区二区| 成人午夜精彩视频在线观看| 听说在线观看完整版免费高清| 免费看不卡的av| 亚洲av一区综合| 99久久九九国产精品国产免费| 欧美成人一区二区免费高清观看| 亚洲电影在线观看av| 18禁裸乳无遮挡免费网站照片| 日本黄大片高清| 久久人人爽人人片av| 欧美日韩精品成人综合77777| 亚洲精品一区蜜桃| 爱豆传媒免费全集在线观看| 99久久精品一区二区三区| 卡戴珊不雅视频在线播放| 神马国产精品三级电影在线观看| 国产日韩欧美亚洲二区| 日本av手机在线免费观看| 久久久久久九九精品二区国产| 亚洲国产精品国产精品| 欧美zozozo另类| 午夜福利高清视频| 日日啪夜夜爽| 欧美日韩视频高清一区二区三区二| 国产成人精品久久久久久| 熟妇人妻不卡中文字幕| 欧美性猛交╳xxx乱大交人| 亚洲av日韩在线播放| 亚洲精品乱码久久久v下载方式| 国产精品精品国产色婷婷| 在线看a的网站| 直男gayav资源| 综合色丁香网| 精品人妻偷拍中文字幕| 国产大屁股一区二区在线视频| 91在线精品国自产拍蜜月| 亚洲精品视频女| 99热这里只有是精品50| 久久韩国三级中文字幕| 2021天堂中文幕一二区在线观| 国产探花在线观看一区二区| 六月丁香七月| 2021天堂中文幕一二区在线观| 丝袜喷水一区| 一二三四中文在线观看免费高清| 2022亚洲国产成人精品| 男人添女人高潮全过程视频|