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

    膠州灣菲律賓蛤仔個(gè)體生長(zhǎng)模型的構(gòu)建與驗(yàn)證*

    2020-07-26 10:35:50董世鵬姜娓娓杜美榮高亞平王曉芹李鳳雪蔣增杰
    漁業(yè)科學(xué)進(jìn)展 2020年4期
    關(guān)鍵詞:蛤仔膠州灣軟體

    董世鵬 藺 凡 姜娓娓 杜美榮 高亞平 王曉芹 李鳳雪 侯 興 蔣增杰,3①

    膠州灣菲律賓蛤仔個(gè)體生長(zhǎng)模型的構(gòu)建與驗(yàn)證*

    董世鵬1,2藺 凡2姜娓娓2杜美榮2高亞平2王曉芹2李鳳雪1,2侯 興1,2蔣增杰2,3①

    (1. 上海海洋大學(xué) 水產(chǎn)科學(xué)國(guó)家級(jí)實(shí)驗(yàn)教學(xué)示范中心 上海 201306;2. 中國(guó)水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所 農(nóng)業(yè)農(nóng)村部海洋漁業(yè)可持續(xù)發(fā)展重點(diǎn)實(shí)驗(yàn)室 青島 266071;3. 青島海洋科學(xué)與技術(shù)試點(diǎn)國(guó)家實(shí)驗(yàn)室海洋漁業(yè)科學(xué)與食物產(chǎn)出過(guò)程功能實(shí)驗(yàn)室 青島 266071)

    為更好地掌握膠州灣菲律賓蛤仔()的動(dòng)態(tài)生長(zhǎng)情況,本研究基于動(dòng)態(tài)能量收支(Dynamic energy budget, DEB)理論,利用python軟件構(gòu)建了菲律賓蛤仔的個(gè)體生長(zhǎng)模型,以2018年4月24日~2019年1月9日觀測(cè)的膠州灣海域葉綠素和水溫為強(qiáng)制函數(shù),通過(guò)現(xiàn)場(chǎng)實(shí)驗(yàn)和已有文獻(xiàn)報(bào)道獲取模型參數(shù),模擬了菲律賓蛤仔軟體組織質(zhì)量和殼長(zhǎng)的生長(zhǎng)情況,并根據(jù)膠州灣海域菲律賓蛤仔生長(zhǎng)的實(shí)測(cè)數(shù)據(jù)對(duì)模型進(jìn)行了驗(yàn)證。結(jié)果表明,構(gòu)建的個(gè)體生長(zhǎng)模型能夠很好地模擬膠州灣海域菲律賓蛤仔軟體組織干重和殼長(zhǎng)的生長(zhǎng),軟體組織干重和殼長(zhǎng)的模擬值與實(shí)測(cè)值呈顯著線性相關(guān)關(guān)系(<0.01),2分別為0.9374和0.9168。敏感度指數(shù)最高的是阿倫紐斯溫度A和參考溫度1,如果A和1分別改變10%,菲律賓蛤仔軟體組織干重增加高達(dá)8.86%。研究結(jié)果為后續(xù)開展基于生態(tài)系統(tǒng)動(dòng)力學(xué)模型的養(yǎng)殖容量動(dòng)態(tài)評(píng)估提供了基礎(chǔ)模塊和數(shù)據(jù)支撐。

    菲律賓蛤仔;動(dòng)態(tài)能量收支;軟體組織干重;殼長(zhǎng);膠州灣

    菲律賓蛤仔()是世界性養(yǎng)殖貝類,也是我國(guó)單種產(chǎn)量最高的養(yǎng)殖貝類,年產(chǎn)量超過(guò)300萬(wàn)t,占世界總產(chǎn)量的 90%(閆喜武等, 2017)。近年來(lái),菲律賓蛤仔養(yǎng)殖產(chǎn)業(yè)發(fā)展迅速,在提高漁民收入、提升海域生態(tài)系統(tǒng)服務(wù)功能等方面發(fā)揮了重要作用(閻希柱, 2009; Chauvaud, 2003; Yang, 2013; 吳亞林等, 2018)。在規(guī)?;B(yǎng)殖情況下,菲律賓蛤仔的生長(zhǎng)情況是關(guān)系到經(jīng)濟(jì)效益高低的關(guān)鍵因素,生長(zhǎng)模型是反映貝類個(gè)體生長(zhǎng)規(guī)律性變化的有效工具。依據(jù)模型建立的方法,貝類生長(zhǎng)模型主要分為純粹生長(zhǎng)曲線模型、經(jīng)驗(yàn)統(tǒng)計(jì)模型、生理化學(xué)模型和生態(tài)生理模型四大類。純粹生長(zhǎng)曲線模型是采用數(shù)學(xué)方程或數(shù)學(xué)模型來(lái)描述貝類殼高、殼長(zhǎng)、殼寬、肉重等生長(zhǎng)相關(guān)性狀隨時(shí)間和年齡的變化規(guī)律,根據(jù)方程在坐標(biāo)圖上畫出生長(zhǎng)相關(guān)性狀隨時(shí)間變化的生長(zhǎng)曲線模型,主要包括Brody、von Bertalanffy、Gompertz、Logistic和Richards等生長(zhǎng)模型。純粹生長(zhǎng)曲線模型未將貝類生長(zhǎng)與環(huán)境條件關(guān)聯(lián),因此,不適用于自然海區(qū)養(yǎng)殖生物的研究(Urban, 2002)。經(jīng)驗(yàn)統(tǒng)計(jì)模型是通過(guò)統(tǒng)計(jì)學(xué)方法建立的貝類生長(zhǎng)速率與餌料濃度、質(zhì)量、溫度等環(huán)境條件的回歸方程,其有助于解釋影響貝類生長(zhǎng)的主要環(huán)境要素,但未考慮貝類的生態(tài)生理學(xué)和適應(yīng)力,也無(wú)法反映貝類生理活動(dòng)與環(huán)境要素之間的反饋?zhàn)饔?Goulletquer, 1988)。生理生態(tài)模型適用于許多食物質(zhì)量起決定作用的生物個(gè)體模型,比如,Powell等(2002)模擬了長(zhǎng)牡蠣()浮游幼體體內(nèi)的碳水化合物、脂類、蛋白質(zhì)和灰分的含量變化,分別模擬了幼蟲體重和殼長(zhǎng)的生長(zhǎng),考慮了遺傳差異和卵的質(zhì)量對(duì)發(fā)育的影響。生態(tài)生理學(xué)模型是對(duì)貝類生長(zhǎng)增量或能量收支方程中的各項(xiàng)生理活動(dòng)分別進(jìn)行模擬后組合成的模型,對(duì)攝食行為進(jìn)行較細(xì)致的刻化,可以動(dòng)態(tài)地反映貝類攝食與餌料供應(yīng)(豐度和質(zhì)量)的相互反饋?zhàn)饔?。目前,?yīng)用最多的生態(tài)生理學(xué)模型是基于生長(zhǎng)余力(Scope for growth)概念建立的SFG模型和基于個(gè)體動(dòng)態(tài)能量收支(Dynamic energy budget)理論的DEB模型(Bayne, 1978; Kooijman, 2000)。SFG模型能夠較全面地反映出能量分配及各因素對(duì)能量分配的影響,可以較好地預(yù)測(cè)貝類的生長(zhǎng),但未能將貝類軟體組織和貝殼有機(jī)質(zhì)之間的能量分開。DEB模型整合了物種所采取的能量分配策略,可以根據(jù)水溫和食物密度等環(huán)境變量的變化預(yù)測(cè)生長(zhǎng)情況(如個(gè)體生長(zhǎng)和繁殖潛力)(van der Meer, 2006; Filgueira, 2011; Sarà, 2013)。DEB模型的重要假設(shè)是-原則,儲(chǔ)存能量一部分()用于結(jié)構(gòu)物質(zhì)的生長(zhǎng)和維持生命活動(dòng),剩余部分(1?)用于繁殖和維持繁育活動(dòng),更適合預(yù)測(cè)雙殼貝類的生長(zhǎng)和繁殖狀態(tài)。

    膠州灣是我國(guó)北方重要的菲律賓蛤仔規(guī)?;B(yǎng)殖海灣,菲律賓蛤仔年產(chǎn)量高達(dá)32.5萬(wàn)t,占膠州灣海水養(yǎng)殖總產(chǎn)量的90%以上(張明亮, 2008; 過(guò)鋒等, 2012),菲律賓蛤仔養(yǎng)殖產(chǎn)業(yè)的健康可持續(xù)發(fā)展對(duì)于充分發(fā)揮膠州灣的生態(tài)服務(wù)功能至關(guān)重要。據(jù)了解,該灣的菲律賓蛤仔播種規(guī)格為2000~3000粒/kg,播種密度為750~1000 kg/畝,播苗時(shí)間為每年的4~5月份,經(jīng)7~8個(gè)月的養(yǎng)成后采收。構(gòu)建個(gè)體生長(zhǎng)模型,對(duì)于準(zhǔn)確掌握膠州灣菲律賓蛤仔的動(dòng)態(tài)生長(zhǎng)情況,更好地指導(dǎo)生產(chǎn)實(shí)踐具有重要的理論和現(xiàn)實(shí)意義。本研究采用python軟件構(gòu)建基于動(dòng)態(tài)能量收支理論的菲律賓蛤仔個(gè)體生長(zhǎng)模型,重點(diǎn)關(guān)注食物(以葉綠素表征)和水溫對(duì)菲律賓蛤仔能量獲取(攝食)和能量消耗(維持、生長(zhǎng)和繁殖)過(guò)程的影響,根據(jù)菲律賓蛤仔野外生長(zhǎng)實(shí)測(cè)數(shù)據(jù)與模擬數(shù)據(jù)的互較進(jìn)行部分參數(shù)的優(yōu)化調(diào)整,旨在模擬膠州灣菲律賓蛤仔個(gè)體的動(dòng)態(tài)生長(zhǎng)過(guò)程,為膠州灣菲律賓蛤仔的養(yǎng)殖容量評(píng)估和適應(yīng)性管理提供決策支持。

    1 材料與方法

    1.1 研究海域

    膠州灣(35°38¢~36°18¢N, 120°04¢~120°23¢E)位于山東半島南岸的西部,瀕臨南黃海西部,為中型半封閉的淺海灣,呈扇形,南北長(zhǎng)約為33.3 km,東西寬約為27.8 km,總面積約為362 km2,潮間帶灘涂面積約為125 km2,平均水深為7 m。

    1.2 環(huán)境參數(shù)

    于2018年4月24日~2019年1月9日,對(duì)膠州灣菲律賓的蛤仔生長(zhǎng)情況和環(huán)境參數(shù)進(jìn)行了跟蹤監(jiān)測(cè)。在蛤仔養(yǎng)殖區(qū)的水溫和葉綠素濃度連續(xù)觀測(cè),數(shù)據(jù)采用自容式濁度葉綠素儀(INFINITY-CLW, 日本)采集。自動(dòng)記錄采樣模式設(shè)定為“脈沖”模式,每隔4 h進(jìn)行1次數(shù)據(jù)監(jiān)測(cè),每次監(jiān)測(cè)3組數(shù)據(jù),每組數(shù)據(jù)間隔時(shí)間為10 min,將3組數(shù)據(jù)取平均值,記為該時(shí)間點(diǎn)的數(shù)據(jù)。王雪(2014)研究表明,膠州灣海域表底層主要環(huán)境因子差異不顯著,表底溫差小。使用3 m的尼龍繩將葉綠素儀拴在船的尾柱上,船只通過(guò)船錨固定在蛤仔養(yǎng)殖區(qū)的預(yù)定位置(36°11¢14.02N, 120°15¢2.92E),葉綠素儀的放置水層為水下1.5 m。

    1.3 菲律賓蛤仔生長(zhǎng)數(shù)據(jù)

    菲律賓蛤仔來(lái)源于膠州灣沿岸農(nóng)場(chǎng)底播養(yǎng)殖的同一批苗種。菲律賓蛤仔苗種的初始平均殼長(zhǎng)為(0.93±0.10) cm,平均軟體組織干重為0.02 g,2018年4月播種到海區(qū)底部,水深為4~7 m。實(shí)驗(yàn)期間,每 30 d左右隨機(jī)取30~50個(gè)蛤仔進(jìn)行生長(zhǎng)數(shù)據(jù)測(cè)定,生長(zhǎng)數(shù)據(jù)主要包括軟體組織干重(g)、軟體組織濕重(g)、殼長(zhǎng)(cm)、殼寬(cm)和殼髙(cm)。

    1.4 模型函數(shù)和參數(shù)

    本研究使用的動(dòng)態(tài)能量收支模型是基于動(dòng)態(tài)能量收支理論建立的,生長(zhǎng)和繁殖的動(dòng)態(tài)可以通過(guò)3個(gè)部分進(jìn)行描述:結(jié)構(gòu)體積()、儲(chǔ)備能量()、繁育能量(R)。假定同化能量首先存儲(chǔ)于體內(nèi),然后在遵循-原則的情況下進(jìn)行維持、生長(zhǎng)和繁殖的能量分配。模型的主要函數(shù)及參數(shù)的含義見表1和表2。Yokogawa (1998)研究表明,不同海域養(yǎng)殖菲律賓蛤仔的形狀差異顯著。本研究采用殼長(zhǎng)和軟體組織濕重回歸法計(jì)算形狀系數(shù)M。半飽和常數(shù)H基于模型校準(zhǔn)獲得,其他參數(shù)的取值主要來(lái)源于文獻(xiàn)(Deslous-Paoli, 1988; Laruelle, 1999; Park, 2004; 張繼紅等, 2005; Kooijman, 2006; van der Veer, 2006; Flye-Sainte-Marie, 2008; Kooijman, 2008)。

    表1 DEB模型中所用的主要關(guān)系式

    Tab.1 Main functions in the DEB model

    1.5 模型運(yùn)行與驗(yàn)證

    使用python 2.7軟件進(jìn)行模型模擬。狀態(tài)變量的初始值由第1次采樣的測(cè)量值估算,繁育儲(chǔ)能R的初始值設(shè)為0。水溫(℃)和葉綠素濃度(μg/L)是模型運(yùn)行的強(qiáng)制函數(shù)。膠州灣菲律賓蛤仔的生長(zhǎng)模擬從2018年4月24日~2019年1月9日,步長(zhǎng)為1 d,共計(jì)運(yùn)行260 d。以逐月獲得的菲律賓蛤仔生長(zhǎng)數(shù)據(jù)(殼長(zhǎng)和軟體組織干重)作為模型的驗(yàn)證。

    1.6 敏感性分析

    為了探索模型在不同環(huán)境中的適應(yīng)性,使用Majkowski(1982)的方法進(jìn)行參數(shù)敏感性分析。每次運(yùn)行模型時(shí),增加或減少與模型輸出結(jié)果相關(guān)的單個(gè)參數(shù)的固定百分比,然后比較模型輸出結(jié)果的相對(duì)變化。敏感性指數(shù)的計(jì)算公式如下:

    對(duì)膠州灣菲律賓蛤仔260 d的軟體組織干重進(jìn)行敏感性分析,增加或減少參數(shù)10%后模型輸出結(jié)果變化的平均值為敏感度指標(biāo)。比較各模型參數(shù)敏感性指標(biāo)的相對(duì)值,確定各參數(shù)的相對(duì)影響程度。

    2 結(jié)果與分析

    2.1 膠州灣表層水體水溫和葉綠素a的連續(xù)性變化情況

    2018年4月24日~2019年1月9日膠州灣菲律賓蛤仔養(yǎng)殖區(qū)水溫和葉綠素濃度的變化如圖1所示。水溫的變化范圍為–0.5℃~31.4℃,最高值出現(xiàn)在8月,最低值出現(xiàn)在12月。葉綠素濃度的變化范圍為0.53~10.89 μg/L,最高值出現(xiàn)在6月,最低值出現(xiàn)在1月。

    表2 菲律賓蛤仔DEB模型的參數(shù)取值

    Tab.2 Parameters of the DEB model of clam Ruditapes philippinarum

    2.2 菲律賓蛤仔軟體組織干重和殼長(zhǎng)的模擬與驗(yàn)證

    2018年4月24日~2019年1月9日膠州灣菲律賓蛤仔的生長(zhǎng)情況和模擬結(jié)果,生長(zhǎng)的實(shí)測(cè)值與模擬值得線性回歸結(jié)果,如圖2所示。結(jié)果顯示,軟體組織干重和殼長(zhǎng)的模擬值與實(shí)測(cè)值呈顯著線性相關(guān)(< 0.01),2分別為0.9374和0.9168。根據(jù)模擬結(jié)果可知,模擬前30 d,菲律賓蛤仔軟體組織干重和殼長(zhǎng)增長(zhǎng)迅速,模擬值略高于實(shí)測(cè)值;30~150 d,菲律賓蛤仔進(jìn)入緩慢增長(zhǎng)期,本段時(shí)間模擬值與實(shí)測(cè)值基本吻合;150~210 d為菲律賓蛤仔的高速增長(zhǎng)期,模擬增速高于實(shí)測(cè)增速;210 d后,增速急劇下降,甚至殼長(zhǎng)增長(zhǎng)出現(xiàn)停滯,軟體組織干重下降。模擬值與實(shí)測(cè)值增長(zhǎng)趨勢(shì)基本一致,表明該模型可以較好地模擬并反映菲律賓蛤仔軟體組織干重和殼長(zhǎng)隨養(yǎng)殖時(shí)間的變化情況。

    圖1 膠州灣葉綠素a濃度和水溫的監(jiān)測(cè)數(shù)據(jù)

    圖2 膠州灣菲律賓蛤仔軟體組織干重、殼長(zhǎng)實(shí)際值與模擬值結(jié)果的比較

    2.3 敏感性分析

    從敏感性分析結(jié)果可以看出,蛤仔軟體組織干重增長(zhǎng)對(duì)大多數(shù)相關(guān)模型參數(shù)的變化不敏感(圖3),表明,該模型在模擬蛤仔軟體組織干重增長(zhǎng)方面是合理的。敏感度指數(shù)最高的是A和1,對(duì)菲律賓蛤仔生長(zhǎng)的影響最大。如果A和1分別改變10%,菲律賓蛤仔軟體組織干重均改變8.86%。

    圖3 DEB模型參數(shù)的敏感性分析結(jié)果

    A:阿倫紐斯溫度;1:參考溫度;FH:半飽和常數(shù);[G]:?jiǎn)挝惑w積結(jié)構(gòu)物質(zhì)所需能量;:結(jié)構(gòu)維持以及生長(zhǎng)能量比例;[M]:最大單位體積儲(chǔ)能;[m]:?jiǎn)挝惑w積維持耗能率;E:儲(chǔ)備能量含量;R:繁殖能量分配比例

    A: Arrhenius temperature;1: Reference temperature; FH: Half saturation coefficient; [G]: Volume-specific costs for structure;: Fraction of catabolic flux to growth and maintenance; [M]: Maximum storage density; [m]: Volume-specific maintenance rate;E: Energy content of reserves;R: Fraction of reproductive reserves

    3 討論

    生長(zhǎng)模擬結(jié)果表明,構(gòu)建的個(gè)體生長(zhǎng)模型能夠較好地模擬膠州灣海區(qū)菲律賓蛤仔的生長(zhǎng)情況。DEB模型在用于牡蠣、貽貝和扇貝等多種濾食性貝類的能量學(xué)研究中,環(huán)境因子的監(jiān)測(cè)多是以月份為單位(Ren, 2001; Filgueira, 2011; 張繼紅等, 2017)。隨著對(duì)于生物生理生態(tài)過(guò)程認(rèn)識(shí)的不斷深入以及觀測(cè)技術(shù)的不斷進(jìn)步,對(duì)數(shù)據(jù)時(shí)間尺度的需求逐漸趨向于長(zhǎng)期性和連續(xù)性。本實(shí)驗(yàn)中,葉綠素和水溫在模擬時(shí)間段內(nèi)每天測(cè)量18次,更加精細(xì)化地模擬了菲律賓蛤仔生長(zhǎng)對(duì)環(huán)境的反饋。

    諸多的研究表明,合適的參數(shù)值對(duì)于DEB模型的應(yīng)用至關(guān)重要(Stéphane, 2006; Troost, 2010)。多數(shù)參數(shù)可以在實(shí)驗(yàn)室內(nèi)進(jìn)行估計(jì)或室外進(jìn)行監(jiān)測(cè),包括溫度和食物。雖然在菲律賓蛤仔進(jìn)行過(guò)許多生理實(shí)驗(yàn),但不同研究中的攝食率和代謝率有很大不同,模型的參數(shù)值也不同(van der Veer, 2006; Flye-Sainte-Marie, 2008)。這種生理速率的變化可能是由于不同種群的物種之前存在遺傳差異導(dǎo)致的(Ren, 2000)。在不改變某些或所有參數(shù)值的情況下,直接將DEB模型應(yīng)用于不同生態(tài)系統(tǒng)中的物種是不合理的。理想情況下,DEB模型的參數(shù)化應(yīng)該基于來(lái)自本地總體的數(shù)據(jù)集。因此,形狀系數(shù)選于膠州灣本地種群測(cè)定。半飽和系數(shù)決定模擬值和實(shí)測(cè)值的擬合度,其值可以通過(guò)模型模擬情況而調(diào)整給出(Stéphane, 2006)。其他模型參數(shù)主要參考van der Meer (2006)和Flye-Sainte-Marie(2008)中的參數(shù)。另外,從本研究的敏感性分析結(jié)果可以看出,阿倫紐斯溫度A、參考溫度1、半飽和系數(shù)FH和單位體積結(jié)構(gòu)物質(zhì)所需能量[G]對(duì)菲律賓蛤仔生長(zhǎng)模型的影響較大。例如,如果將[G]提高10%,模擬出的貝類生長(zhǎng)速率可提高6%。因此,這些參數(shù)應(yīng)謹(jǐn)慎賦值。

    底播養(yǎng)殖的菲律賓蛤仔生長(zhǎng)情況受食物和水溫的雙重調(diào)控。根據(jù)葉綠素、水溫與模擬結(jié)果比較,膠州灣海域模擬過(guò)程中大部分時(shí)間水溫對(duì)蛤仔生長(zhǎng)控制更明顯。4、5月,水溫開始回升,蛤仔幼苗快速生長(zhǎng)。6~8月,食物到達(dá)一年中的峰值,但夏季高溫使得菲律賓蛤仔幾乎停止生長(zhǎng)。9~11月中旬,食物與水溫都相對(duì)適宜,蛤仔生長(zhǎng)情況良好。進(jìn)入冬季后,伴隨著食物減少和水溫降低,蛤仔軟體組織出現(xiàn)停止增長(zhǎng)甚至萎縮的情況。因此,冬季低溫和食物限制可能會(huì)導(dǎo)致菲律賓蛤仔存儲(chǔ)能量耗盡,影響其生理狀態(tài)或造成結(jié)構(gòu)物質(zhì)損傷。

    本研究建立的菲律賓蛤仔動(dòng)態(tài)能量收支模型,在給定的海水溫度和食物濃度等環(huán)境條件下,較好地模擬了蛤仔在自然條件下的生長(zhǎng)狀況。但是,菲律賓蛤仔的生長(zhǎng)情況受自然海域環(huán)境的影響很大,而且這種環(huán)境不是人為可控的,所以該模型的周年穩(wěn)定性需要進(jìn)一步驗(yàn)證。菲律賓蛤仔的食物來(lái)源復(fù)雜,該模型的食物來(lái)源僅僅采用了葉綠素。據(jù)報(bào)道,對(duì)食物來(lái)源的精細(xì)化甄別可以進(jìn)一步提升模型模擬的準(zhǔn)確性,如Troost等(2010)使用葉綠素和顆粒有機(jī)物(POM)計(jì)算食物密度,更好地模擬了鳥蛤()在食物條件較差的Oosterschelde海灣的生長(zhǎng)情況。此外,在其他海洋動(dòng)物動(dòng)態(tài)模型研究過(guò)程中,雌雄模型具有一定差異,這也可能是本模型模擬值和實(shí)測(cè)值存在差異的原因(王新安等, 2014)。因此,后續(xù)研究將在現(xiàn)有的基礎(chǔ)上,對(duì)食物來(lái)源、環(huán)境影響和雌雄差異等進(jìn)行細(xì)致解析,進(jìn)一步提高模型的準(zhǔn)確性和穩(wěn)定性,以期更好地為生態(tài)系統(tǒng)動(dòng)力學(xué)模型的構(gòu)建以及養(yǎng)殖容量的評(píng)估提供基礎(chǔ)數(shù)據(jù)。

    Bayne BL, Widdows J. The physiological ecology of two populations ofL. Oecologia, 1978, 37(2): 137–162

    Chauvaud L, Thompson JK, Cloern JE,. Clams as CO2generators: Theexample in San Francisco Bay. Limnology and Oceanography, 2003, 48(6): 2086–2092

    Deslous-Paoli JM, Héral M. Biochemical composition and energy value of(Thunberg) cultured in the bay of Merennes-Oléron. Aquatic Living Resources, 1988, 1(4): 239–249

    Filgueira R, Rosland R, Grant J. A comparison of scope for growth (SFG) and dynamic energy budget (DEB) models applied to the blue mussel (). Journal of Sea Research, 2011, 66(4): 403–410

    Flye-Sainte-Marie JG. Ecophysiology of brown ring disease in the Manila clam, experimental and modelling approaches. Amsterdam: Vrije Universiteit, 2008

    Goulletquer P, Bacher C. Empirical modelling of the growth ofby means of non linear regression on factorial coordinates. Aquatic living resources, 1988, 1(3): 141–154

    Guo F, Zhao J, Chen JF,. Nitrogen and phosphorous pollution in shellfish culture areas of Jiaozhou Bay. Progress in Fishery Sciences, 2012, 33(5): 116–122 [過(guò)鋒, 趙俊, 陳聚法, 等. 膠州灣貝類養(yǎng)殖區(qū)氮、磷污染現(xiàn)狀及動(dòng)態(tài)變化. 漁業(yè)科學(xué)進(jìn)展, 2012, 33(5): 116–122]

    Kooijman SA, Sousa T, Pecquerie L,. From food-dependent statistics to metabolic parameters, a practical guide to the use of dynamic energy budget theory. Biological Reviews, 2008, 83(4): 533–552

    Kooijman SALM. Dynamic energy and mass budgets in biological systems. Cambridge: Cambridge University Press, 2000

    Kooijman SALM. Pseudo-faeces production in bivalves. Journal of Sea Research, 2006, 56(2): 103–106

    Laruelle F. Phénologie et déterminisme de la reproduction chez(L.) et(Adams et Reeve) en Bretagne. Doctoral Dissertation of Université de Bretagne Occidentale, 1999

    Majkowski JV. Usefulness and applicability of sensitivity analysis in a multispecies approach to fisheries management. Theory and management of tropical fisheries. Research Publications Repository, 1982, 149–165

    Park K, Choi K. Application of enzyme-linked immunosorbent assay for studying of reproduction in the Manila clam(Mollusca: Bivalvia): I. Quantifying eggs. Aquaculture, 2004, 241(1–4): 667–687

    Powell EN, Bochenek EA, Klinck JM,. Influence of food quality and quantity on the growth and development of, larvae: A modeling approach. Aquaculture, 2002, 210(1): 89–117

    Ren JS, Ross AH, Schiel DR. Functional descriptions of feeding and energetics of the Pacific oysterin New Zealand. Marine Ecology Progress Series, 2000, 208(3): 119–130

    Ren JS, Ross AH. A dynamic energy budget model of the Pacific oyster. Ecological Modelling, 2001, 142(1–2): 105–120

    Sarà G, Palmeri V, Montalto V,. Parameterisation of bivalve functional traits for mechanistic eco-physiological dynamic energy budget (DEB) models. Marine Ecology Progress Series, 2013, 480: 99–117

    Stéphane P, Bourles Y, Sébastien L,. Application of a dynamic energy budget model to the Pacific oyster,, reared under various environmental conditions. Journal of Sea Research, 2006, 56(2): 156–167

    Troost TA, Wijsman JWM, Saraiva S,. Modelling shellfish growth with dynamic energy budget models: An application for cockles and mussels in the Oosterschelde (southwest Netherlands). Philosophical Transactions of the Royal Society of London, 2010, 365(1557): 3567–3577

    Urban HJ. Modeling growth of different developmental stages in bivalves. Marine Ecology Progress Series, 2002, 238: 109– 114

    van der Meer J. An introduction to dynamic energy budget (DEB) models with special emphasis on parameter estimation. Journal of Sea Research, 2006, 56(2): 85–102

    van der Veer HW, Cardoso JFMF, van der Meer J. The estimation of DEB parameters for various Northeast Atlantic bivalve species. Journal of Sea Research, 2006, 56(2): 107–124

    Wang XA, Ma AJ, Huang ZH,. Developmental differences between female and male groups in turbot () breeding families. Journal of Fisheries of China, 2014, 38(4): 465–470 [王新安, 馬愛(ài)軍, 黃智慧, 等. 大菱鲆選育家系雌、雄群體的生長(zhǎng)發(fā)育差異. 水產(chǎn)學(xué)報(bào), 2014, 38(4): 465–470]

    Wang X. Numerical study to hydrodynamic conditions of the Jiaozhou Bay. Master′s Thesis of Ocean University of China, 2014 [王雪. 膠州灣三維水動(dòng)力數(shù)值模擬研究. 中國(guó)海洋大學(xué)碩士研究生學(xué)位論文, 2014]

    Wu YL, Gao YP, Lü XN,. Ecological contribution ofin the seagrass meadow of Sanggou Bay. Progress in Fishery Sciences, 2018, 39(6): 126–133 [吳亞林, 高亞平, 呂旭寧, 等. 桑溝灣楮島大葉藻床區(qū)域菲律賓蛤仔的生態(tài)貢獻(xiàn). 漁業(yè)科學(xué)進(jìn)展, 2018, 39(6): 126–133]

    Yan XW, Huo ZM, Zhang YH,.white zebra. China Fisheries, 2017(12): 72–74 [閆喜武, 霍忠明, 張躍環(huán), 等. 菲律賓蛤仔白斑馬蛤. 中國(guó)水產(chǎn), 2017(12): 72–74]

    Yan XZ. Study on the filtration rates ofin Zhaoan Gulf. Hubei Agricultural Sciences, 2009, 48(11): 2808–2810 [閻希柱. 詔安灣菲律賓蛤仔的濾水率研究. 湖北農(nóng)業(yè)科學(xué), 2009, 48(11): 2808–2810]

    Yang F, Zhao L, Yan X,. Bioaccumulation of trace elements infrom China: Public health risk assessment implications. International Journal of Environmental Research and Public Health, 2013, 10(4): 1392–1405

    Yokogawa K. Morphological variability and genetic features in Japanese common clam. Venus (Japanese Journal of Malacology), 1998, 57: 121–131

    Zhang JH, Fang JG, Sun S,. Clearance rate, ingestion rate and absorption efficiency of cultivated clamin Jiaozhou Bay. Oceanologia et Limnologia Sinica, 2005, 36(6): 548–555 [張繼紅, 方建光, 孫松, 等. 膠州灣養(yǎng)殖菲律賓蛤仔的清濾率、攝食率、吸收效率的研究. 海洋與湖沼, 2005, 36(6): 548–555]

    Zhang JH, Wu WG, Liu Y,. A dynamic energy budget (DEB) growth model for Japanese scallopcultured in China. Journal of Fishery Sciences of China, 2017, 24(3): 497–506 [張繼紅, 吳文廣, 劉毅, 等. 蝦夷扇貝動(dòng)態(tài)能量收支生長(zhǎng)模型. 中國(guó)水產(chǎn)科學(xué), 2017, 24(3): 497–506]

    Zhang ML. The research of carrying capacity of shellfish in the Jiaozhou Bay. Master′s Thesis of the First Institute of Oceanography, National Oceanic Administration, 2008 [張明亮. 膠州灣貝類養(yǎng)殖容量研究與分析. 國(guó)家海洋局第一海洋研究所碩士研究生學(xué)位論文, 2008]

    Simulation and Validation of Manila ClamGrowth with a DEB-Based Individual Growth Model in Jiaozhou Bay

    DONG Shipeng1,2, LIN Fan2, JIANG Weiwei2, DU Meirong2, GAO Yaping2, WANG Xiaoqin2, LI Fengxue1,2, HOU Xing1,2, JIANG Zengjie2,3①

    (1. National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai 201306; 2. Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Key Laboratory of Sustainable Development of Marine Fisheries, Ministry of Agriculture and Rural Affairs, Qingdao 266071; 3. Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology(Qingdao), Qingdao 266071)

    In order to better understand the dynamic growth condition of the Manila clamin Jiaozhou Bay, an individual growth model forwas developed based on the dynamic energy budget (DEB) theory. Seawater chlorophyll-concentration and temperature were considered as influencing variables and data were collected from April 24, 2018 to January 9, 2019. The DEB parameters for the clam were calculated based on field experiment measurements and previous reported data from literatures. The growth of dry tissue weight and shell length ofwere simulated in the model. Field observations (dry tissue weight and shell length) were used to validate the sensitivities of the clam-DEB model. The simulated results concur with the observed results. Regression analysis showed that there was a significant linear correlation between the simulated and observed values of dry tissue weight and shell length (2=0.9374 and 0.9168), respectively (<0.01). The model was relatively sensitive to the variations ofAand1. IfAand1increased by 10%, the dry tissue weight of the clam showed an increase of 8.86%. The results in this study established through the DEB model provide a helpful basic module and data support for the subsequent dynamic assessment of aquaculture carrying capacity based on the ecosystem model.

    ; Dynamic energy budget; Dry tissue weight; Shell length; Jiaozhou Bay

    JIANG Zengjie, E-mail: jiangzj@ysfri.ac.cn

    S967.2

    A

    2095-9869(2020)04-0094-08

    10.19663/j.issn2095-9869.20190409001

    http://www.yykxjz.cn/

    董世鵬, 藺凡, 姜娓娓, 杜美榮, 高亞平, 王曉芹, 李鳳雪, 侯興, 蔣增杰. 膠州灣菲律賓蛤仔個(gè)體生長(zhǎng)模型的構(gòu)建與驗(yàn)證. 漁業(yè)科學(xué)進(jìn)展, 2020, 41(4): 94–101

    Dong SP, Lin F, Jiang WW, Du MR, Gao YP, Wang XQ, Li FX, Hou X, Jiang ZJ. Simulation and validation of Manila clamgrowth with a DEB-based individual growth model in Jiaozhou Bay. Progress in Fishery Sciences, 2020, 41(4): 94–101

    * 科技部政府間國(guó)際科技創(chuàng)新合作重點(diǎn)專項(xiàng)(2017YFE0118300)、國(guó)家自然科學(xué)基金國(guó)際(地區(qū))合作與交流項(xiàng)目(41761134052)、青島海洋科學(xué)與技術(shù)試點(diǎn)國(guó)家實(shí)驗(yàn)室海洋漁業(yè)科學(xué)與食物產(chǎn)出過(guò)程功能實(shí)驗(yàn)室漁業(yè)科技青年人才計(jì)劃項(xiàng)目(2018-MFS-T13)和現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系(CARS-49)共同資助[This work was supported by Key Program for International Cooperation on Scientific and Technological Innovation, Ministry of Science and Technology (2017YFE0118300); International Cooperation and Exchange of the National Natural Science Foundation of China (41761134052); Youth Talent Program Supported by Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology (Qingdao)(2018-MFS-T13), and China Agriculture Research System (CARS-49)]. 董世鵬,E-mail: 742119697@qq.com

    蔣增杰,研究員,E-mail: jiangzj@ysfri.ac.cn

    2019-04-09,

    2019-05-05

    (編輯 馬璀艷)

    猜你喜歡
    蛤仔膠州灣軟體
    蚯蚓軟體機(jī)器人來(lái)到火星變身“農(nóng)夫”
    軍事文摘(2024年4期)2024-01-09 09:07:58
    菲律賓蛤仔清洗分級(jí)整機(jī)設(shè)計(jì)及參數(shù)優(yōu)化
    遼寧蛤仔產(chǎn)業(yè)發(fā)展初探
    晶格型模塊化軟體機(jī)器人自重構(gòu)序列
    會(huì)捉蒼蠅的高速軟體機(jī)器人問(wèn)世
    軍事文摘(2020年24期)2020-11-29 06:35:54
    平流霧罩,海上蓬萊膠州灣
    不同海藻餌料對(duì)菲律賓蛤仔生長(zhǎng)的影響
    軟體家具檢測(cè)國(guó)標(biāo)將出臺(tái)
    膠州灣夏季鹽度長(zhǎng)期輸運(yùn)機(jī)制分析
    兩種水溫條件下鎘對(duì)菲律賓蛤仔的急性毒性
    亚洲成人免费电影在线观看| 精品久久久久久久久久久久久| 欧美成人一区二区免费高清观看| 欧美性猛交╳xxx乱大交人| 久久久久久久久大av| 99久久精品热视频| 国产精品久久久久久精品电影| 亚洲精品成人久久久久久| 欧美+亚洲+日韩+国产| av在线蜜桃| 九色国产91popny在线| 亚洲不卡免费看| 国产精品伦人一区二区| 国产免费男女视频| 亚洲国产精品999在线| 90打野战视频偷拍视频| 97超级碰碰碰精品色视频在线观看| 美女被艹到高潮喷水动态| 动漫黄色视频在线观看| 欧美精品国产亚洲| 久久性视频一级片| 免费观看精品视频网站| 每晚都被弄得嗷嗷叫到高潮| 欧美一区二区国产精品久久精品| 一级作爱视频免费观看| 成人欧美大片| 99国产极品粉嫩在线观看| 床上黄色一级片| 精品欧美国产一区二区三| 精品人妻一区二区三区麻豆 | 午夜福利在线观看吧| 午夜福利成人在线免费观看| 亚洲最大成人中文| eeuss影院久久| 一级毛片久久久久久久久女| 特大巨黑吊av在线直播| 激情在线观看视频在线高清| 九九久久精品国产亚洲av麻豆| 成人午夜高清在线视频| www.熟女人妻精品国产| 国产精品亚洲一级av第二区| 亚洲国产欧洲综合997久久,| 可以在线观看的亚洲视频| 中文字幕久久专区| 亚洲 国产 在线| 国产一区二区在线av高清观看| 日韩亚洲欧美综合| 久久久久久久久中文| 欧美+亚洲+日韩+国产| 亚洲av日韩精品久久久久久密| 欧美色视频一区免费| 高潮久久久久久久久久久不卡| 日本黄大片高清| 久久久久久久久久黄片| 久久久色成人| 老女人水多毛片| 日本成人三级电影网站| 男女床上黄色一级片免费看| 18禁黄网站禁片午夜丰满| 丁香六月欧美| 亚洲成人精品中文字幕电影| 日本免费a在线| 国产69精品久久久久777片| 国产亚洲精品av在线| 久久久久亚洲av毛片大全| 欧美性猛交黑人性爽| 日韩 亚洲 欧美在线| 久久久久久久久中文| 婷婷精品国产亚洲av在线| 亚洲最大成人av| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品av在线| 久久人人精品亚洲av| 俄罗斯特黄特色一大片| 精品99又大又爽又粗少妇毛片 | 男人舔女人下体高潮全视频| 午夜福利视频1000在线观看| 久久久精品大字幕| 此物有八面人人有两片| 超碰av人人做人人爽久久| 免费看日本二区| 久久久久国产精品人妻aⅴ院| 中文字幕人成人乱码亚洲影| 十八禁人妻一区二区| 最近最新中文字幕大全电影3| 亚洲电影在线观看av| 波多野结衣高清作品| 亚洲国产色片| netflix在线观看网站| av天堂在线播放| or卡值多少钱| 国产精品亚洲av一区麻豆| 欧美激情国产日韩精品一区| 草草在线视频免费看| 久久天躁狠狠躁夜夜2o2o| 久久久久精品国产欧美久久久| 亚洲成人免费电影在线观看| a级毛片a级免费在线| 老司机午夜十八禁免费视频| 成人国产综合亚洲| 亚洲av.av天堂| 午夜免费激情av| 久久久精品大字幕| 亚洲专区中文字幕在线| 成人午夜高清在线视频| 国内久久婷婷六月综合欲色啪| 成人特级黄色片久久久久久久| 日本 av在线| 国产成人福利小说| 久久精品综合一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 欧美黑人巨大hd| 欧美3d第一页| 午夜视频国产福利| 午夜两性在线视频| 少妇丰满av| 国产精品日韩av在线免费观看| 91久久精品国产一区二区成人| 三级毛片av免费| 亚洲人与动物交配视频| 好男人在线观看高清免费视频| av黄色大香蕉| 国产一区二区亚洲精品在线观看| 人人妻人人澡欧美一区二区| 88av欧美| 久久性视频一级片| 亚洲久久久久久中文字幕| 国产成+人综合+亚洲专区| 国产乱人视频| 久久中文看片网| 色精品久久人妻99蜜桃| 性欧美人与动物交配| 亚洲在线观看片| 亚洲最大成人中文| ponron亚洲| 亚洲精品色激情综合| 偷拍熟女少妇极品色| 给我免费播放毛片高清在线观看| 国产单亲对白刺激| 好看av亚洲va欧美ⅴa在| 国产高清视频在线播放一区| 欧美又色又爽又黄视频| 亚洲av.av天堂| 毛片一级片免费看久久久久 | 搡老岳熟女国产| 天堂动漫精品| 一区二区三区四区激情视频 | 亚洲人成电影免费在线| 成年女人毛片免费观看观看9| 国产精品久久久久久久久免 | 最好的美女福利视频网| 久久婷婷人人爽人人干人人爱| 国产黄色小视频在线观看| 中出人妻视频一区二区| 久久精品国产亚洲av涩爱 | 美女xxoo啪啪120秒动态图 | 久久精品久久久久久噜噜老黄 | 欧美一区二区国产精品久久精品| 简卡轻食公司| 99久久99久久久精品蜜桃| 自拍偷自拍亚洲精品老妇| 一个人免费在线观看电影| 国产69精品久久久久777片| 黄色一级大片看看| 久久这里只有精品中国| 中国美女看黄片| 免费黄网站久久成人精品 | 一个人观看的视频www高清免费观看| x7x7x7水蜜桃| 国产欧美日韩精品一区二区| 日本黄大片高清| 免费高清视频大片| 美女cb高潮喷水在线观看| 九色国产91popny在线| 熟女电影av网| 热99在线观看视频| 一夜夜www| 国产日本99.免费观看| 欧美高清成人免费视频www| 国产免费一级a男人的天堂| 一二三四社区在线视频社区8| 最后的刺客免费高清国语| 亚洲av电影不卡..在线观看| 午夜日韩欧美国产| 精品人妻偷拍中文字幕| 高清日韩中文字幕在线| 18美女黄网站色大片免费观看| 午夜福利成人在线免费观看| 免费av不卡在线播放| 国产私拍福利视频在线观看| 亚洲五月天丁香| 精品久久久久久久久久久久久| 欧美不卡视频在线免费观看| 身体一侧抽搐| 黄片小视频在线播放| 制服丝袜大香蕉在线| 国产精品美女特级片免费视频播放器| 观看免费一级毛片| 国产毛片a区久久久久| 久久久久久九九精品二区国产| 亚洲av美国av| 99国产精品一区二区蜜桃av| 中出人妻视频一区二区| 国产黄片美女视频| 日韩欧美一区二区三区在线观看| 亚洲av第一区精品v没综合| 欧美黄色淫秽网站| 免费一级毛片在线播放高清视频| 午夜免费成人在线视频| 黄色视频,在线免费观看| 欧美丝袜亚洲另类 | 国产v大片淫在线免费观看| 亚洲午夜理论影院| 97碰自拍视频| 亚洲av成人av| 青草久久国产| 国内精品久久久久精免费| 熟女电影av网| 久久这里只有精品中国| 日日干狠狠操夜夜爽| 在线看三级毛片| 亚洲成av人片在线播放无| 人人妻,人人澡人人爽秒播| a级毛片免费高清观看在线播放| 国产真实乱freesex| 国产又黄又爽又无遮挡在线| 欧美日韩瑟瑟在线播放| 国产av在哪里看| 黄片小视频在线播放| 亚洲狠狠婷婷综合久久图片| 一区二区三区免费毛片| 99热这里只有是精品在线观看 | 深爱激情五月婷婷| 精品免费久久久久久久清纯| 亚洲精品久久国产高清桃花| 女人被狂操c到高潮| 国产亚洲精品久久久com| 在线免费观看不下载黄p国产 | 在线免费观看不下载黄p国产 | 亚洲精品乱码久久久v下载方式| 中文字幕久久专区| 在线十欧美十亚洲十日本专区| 精品国产亚洲在线| 大型黄色视频在线免费观看| 色综合亚洲欧美另类图片| 国产一区二区激情短视频| 毛片一级片免费看久久久久 | 高清日韩中文字幕在线| 亚洲av日韩精品久久久久久密| 国产精品久久久久久人妻精品电影| 真人一进一出gif抽搐免费| 亚洲av不卡在线观看| 小蜜桃在线观看免费完整版高清| 亚洲色图av天堂| 悠悠久久av| 精品国产三级普通话版| 听说在线观看完整版免费高清| 97热精品久久久久久| 亚洲一区高清亚洲精品| 国产在线精品亚洲第一网站| 日韩国内少妇激情av| 中文字幕av成人在线电影| 欧美午夜高清在线| 观看美女的网站| 在现免费观看毛片| 黄色丝袜av网址大全| 日韩国内少妇激情av| 国产精品综合久久久久久久免费| 久久这里只有精品中国| 欧美性感艳星| 男人和女人高潮做爰伦理| 亚洲av电影不卡..在线观看| 波多野结衣巨乳人妻| 国产一级毛片七仙女欲春2| 亚洲精品成人久久久久久| 丰满乱子伦码专区| 亚洲欧美清纯卡通| 91麻豆av在线| 午夜日韩欧美国产| 国产一区二区三区视频了| 精品国产三级普通话版| а√天堂www在线а√下载| 老司机午夜十八禁免费视频| 麻豆国产av国片精品| 日韩欧美国产在线观看| 国产免费一级a男人的天堂| 亚洲av日韩精品久久久久久密| 最新在线观看一区二区三区| 国产精品亚洲美女久久久| 美女免费视频网站| 毛片一级片免费看久久久久 | 亚洲黑人精品在线| 欧美色欧美亚洲另类二区| 亚洲av成人精品一区久久| 欧美日韩综合久久久久久 | 中文字幕av成人在线电影| 精品一区二区三区人妻视频| 一a级毛片在线观看| 中文字幕熟女人妻在线| 日韩有码中文字幕| 久久人妻av系列| 日本熟妇午夜| 日本五十路高清| 有码 亚洲区| 亚洲国产色片| 精品一区二区三区av网在线观看| av福利片在线观看| av国产免费在线观看| 90打野战视频偷拍视频| 国产人妻一区二区三区在| 成人美女网站在线观看视频| 又紧又爽又黄一区二区| 搡女人真爽免费视频火全软件 | 欧美成人免费av一区二区三区| 人人妻人人看人人澡| 精品人妻偷拍中文字幕| 亚洲精品粉嫩美女一区| 男人和女人高潮做爰伦理| 亚洲av成人精品一区久久| 欧美区成人在线视频| 久久精品夜夜夜夜夜久久蜜豆| 51午夜福利影视在线观看| 国产亚洲欧美98| 久久久久久久久中文| 又粗又爽又猛毛片免费看| 精品久久久久久久久av| 男人和女人高潮做爰伦理| 久久国产精品人妻蜜桃| 欧美日韩综合久久久久久 | 国产白丝娇喘喷水9色精品| 国产一区二区在线观看日韩| а√天堂www在线а√下载| 免费看日本二区| 天天一区二区日本电影三级| 一级a爱片免费观看的视频| 91在线精品国自产拍蜜月| 成年女人毛片免费观看观看9| 国产免费男女视频| 日韩欧美精品v在线| 欧美日韩乱码在线| 性欧美人与动物交配| 亚洲美女视频黄频| 亚洲狠狠婷婷综合久久图片| 在线看三级毛片| 在线观看舔阴道视频| 2021天堂中文幕一二区在线观| 免费看美女性在线毛片视频| 最新在线观看一区二区三区| 国产精品99久久久久久久久| 日韩欧美精品v在线| 制服丝袜大香蕉在线| 欧美中文日本在线观看视频| 老司机深夜福利视频在线观看| 97人妻精品一区二区三区麻豆| 日韩欧美免费精品| 亚洲一区高清亚洲精品| 国产伦一二天堂av在线观看| 97人妻精品一区二区三区麻豆| 女人十人毛片免费观看3o分钟| 亚洲av不卡在线观看| 欧美最新免费一区二区三区 | 我要搜黄色片| 欧美3d第一页| 黄色配什么色好看| 天堂动漫精品| 色综合欧美亚洲国产小说| 亚洲第一区二区三区不卡| 制服丝袜大香蕉在线| 精品久久久久久久久久免费视频| 亚洲精品一区av在线观看| 亚洲精品久久国产高清桃花| 赤兔流量卡办理| www.色视频.com| 三级毛片av免费| 国产免费男女视频| 一a级毛片在线观看| 国模一区二区三区四区视频| 老女人水多毛片| 老熟妇仑乱视频hdxx| 99国产综合亚洲精品| 日本一本二区三区精品| 动漫黄色视频在线观看| 男女床上黄色一级片免费看| 亚洲成av人片免费观看| 精品久久久久久久人妻蜜臀av| 精品久久久久久,| 欧美激情在线99| 精品久久久久久久久久久久久| 永久网站在线| 久久久久久大精品| 国产精品久久久久久久久免 | 日韩成人在线观看一区二区三区| 日本a在线网址| 中文字幕人妻熟人妻熟丝袜美| 精品一区二区三区视频在线观看免费| 欧美黑人欧美精品刺激| 国产真实乱freesex| 久久久久免费精品人妻一区二区| 一级a爱片免费观看的视频| 97热精品久久久久久| 脱女人内裤的视频| 午夜免费成人在线视频| 日日干狠狠操夜夜爽| av黄色大香蕉| 中文字幕av成人在线电影| 日本免费a在线| 高清在线国产一区| 两个人的视频大全免费| 午夜福利免费观看在线| 又爽又黄无遮挡网站| 国语自产精品视频在线第100页| 中文字幕熟女人妻在线| 成人美女网站在线观看视频| 日韩欧美一区二区三区在线观看| 亚洲第一欧美日韩一区二区三区| 日本五十路高清| 非洲黑人性xxxx精品又粗又长| 人妻夜夜爽99麻豆av| 久久久久国内视频| 成人高潮视频无遮挡免费网站| 亚洲18禁久久av| bbb黄色大片| 91九色精品人成在线观看| 丰满的人妻完整版| 一级作爱视频免费观看| 欧美日韩国产亚洲二区| 男人舔女人下体高潮全视频| 国产成人影院久久av| 免费一级毛片在线播放高清视频| 亚洲精品一区av在线观看| 久久午夜福利片| 看免费av毛片| 色在线成人网| 99riav亚洲国产免费| 18+在线观看网站| 一a级毛片在线观看| 麻豆国产av国片精品| 一本综合久久免费| 国产欧美日韩一区二区三| 51国产日韩欧美| 伦理电影大哥的女人| 久久精品影院6| 国产v大片淫在线免费观看| 一区二区三区激情视频| 黄色视频,在线免费观看| 久久久久精品国产欧美久久久| 国产乱人视频| 亚洲国产精品成人综合色| 天堂网av新在线| 91av网一区二区| 久久久色成人| 国产免费男女视频| 亚洲一区二区三区不卡视频| 美女免费视频网站| 午夜激情欧美在线| 国产精品久久久久久久电影| 免费观看人在逋| 日韩大尺度精品在线看网址| 3wmmmm亚洲av在线观看| 久9热在线精品视频| 网址你懂的国产日韩在线| 露出奶头的视频| 国产又黄又爽又无遮挡在线| 人人妻人人澡欧美一区二区| 欧美乱妇无乱码| 伦理电影大哥的女人| 90打野战视频偷拍视频| 十八禁国产超污无遮挡网站| 日韩成人在线观看一区二区三区| 好男人在线观看高清免费视频| 激情在线观看视频在线高清| 一个人免费在线观看的高清视频| 真实男女啪啪啪动态图| eeuss影院久久| x7x7x7水蜜桃| 琪琪午夜伦伦电影理论片6080| 精品午夜福利视频在线观看一区| 国产极品精品免费视频能看的| 亚洲欧美日韩卡通动漫| 亚洲中文日韩欧美视频| 香蕉av资源在线| 亚洲精品456在线播放app | 亚洲经典国产精华液单 | 99精品在免费线老司机午夜| 中文字幕人成人乱码亚洲影| 国产成人av教育| 变态另类成人亚洲欧美熟女| 99国产综合亚洲精品| 午夜精品一区二区三区免费看| 一级a爱片免费观看的视频| 亚洲国产色片| 在线观看66精品国产| 一级毛片久久久久久久久女| 99热精品在线国产| 婷婷六月久久综合丁香| 99久久精品国产亚洲精品| 午夜福利在线观看免费完整高清在 | 亚洲av中文字字幕乱码综合| 欧美午夜高清在线| 男人狂女人下面高潮的视频| 在线观看av片永久免费下载| 久久久久国内视频| 黄色一级大片看看| 国产精品自产拍在线观看55亚洲| av在线老鸭窝| 国产高清视频在线观看网站| 亚洲人成网站高清观看| 亚洲熟妇中文字幕五十中出| 亚洲美女搞黄在线观看 | 免费在线观看成人毛片| 日韩中文字幕欧美一区二区| 久久久久久久精品吃奶| 日日摸夜夜添夜夜添av毛片 | 露出奶头的视频| 看黄色毛片网站| 国产一区二区三区在线臀色熟女| 久久久久久久久中文| 亚洲av电影在线进入| 99久久99久久久精品蜜桃| 久久久久久大精品| 色综合婷婷激情| 最近视频中文字幕2019在线8| 小蜜桃在线观看免费完整版高清| 色综合站精品国产| 中文字幕av在线有码专区| 亚洲18禁久久av| 中文在线观看免费www的网站| 男女下面进入的视频免费午夜| 成人鲁丝片一二三区免费| 久久久久亚洲av毛片大全| 麻豆成人午夜福利视频| 成年女人看的毛片在线观看| 午夜福利在线观看吧| .国产精品久久| 午夜福利视频1000在线观看| 白带黄色成豆腐渣| 免费看美女性在线毛片视频| 国产欧美日韩精品一区二区| 亚州av有码| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 黄片小视频在线播放| 日韩欧美在线二视频| 少妇丰满av| 免费在线观看成人毛片| 日本免费一区二区三区高清不卡| 一本久久中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 亚洲七黄色美女视频| 久久久久久久久久成人| 我要看日韩黄色一级片| 国产精品野战在线观看| 国产av不卡久久| 怎么达到女性高潮| 欧美黑人欧美精品刺激| 最近最新中文字幕大全电影3| 欧美午夜高清在线| 别揉我奶头 嗯啊视频| 久久精品国产自在天天线| 久久精品国产99精品国产亚洲性色| 成人国产一区最新在线观看| 久9热在线精品视频| 国产一级毛片七仙女欲春2| 亚洲欧美清纯卡通| 99视频精品全部免费 在线| 中国美女看黄片| 免费av不卡在线播放| 老鸭窝网址在线观看| 国产精品国产高清国产av| 91在线精品国自产拍蜜月| 国产探花极品一区二区| 在线观看美女被高潮喷水网站 | 麻豆久久精品国产亚洲av| 久久九九热精品免费| 一区二区三区激情视频| a在线观看视频网站| 99视频精品全部免费 在线| 亚洲精品影视一区二区三区av| 国产精品国产高清国产av| 欧美激情在线99| 黄色一级大片看看| 69av精品久久久久久| 黄色一级大片看看| 女人十人毛片免费观看3o分钟| 欧美激情久久久久久爽电影| 69av精品久久久久久| 网址你懂的国产日韩在线| 日本 欧美在线| 久久国产精品影院| 亚洲欧美清纯卡通| 亚洲中文字幕一区二区三区有码在线看| 午夜福利视频1000在线观看| 亚洲电影在线观看av| 久久午夜亚洲精品久久| 精品一区二区三区人妻视频| 久9热在线精品视频| ponron亚洲| 国产中年淑女户外野战色| 最好的美女福利视频网| 亚洲一区高清亚洲精品| 婷婷色综合大香蕉| 内地一区二区视频在线| 久久久久久久久大av| 国产成人a区在线观看| 亚洲精品粉嫩美女一区| 久久久久久久久中文| 日韩欧美免费精品| 757午夜福利合集在线观看| 91久久精品国产一区二区成人| 国产乱人伦免费视频| www.www免费av| 小说图片视频综合网站| 亚洲国产色片| 少妇人妻一区二区三区视频| 天堂影院成人在线观看| 欧美日韩中文字幕国产精品一区二区三区|