• <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ì)菲律賓蛤仔的急性毒性
    观看免费一级毛片| 国产黄频视频在线观看| 亚洲精品乱久久久久久| 丝袜美腿在线中文| 亚洲精品国产av蜜桃| 精品一区二区三卡| 精品国产一区二区三区久久久樱花 | 亚洲丝袜综合中文字幕| 网址你懂的国产日韩在线| 2021天堂中文幕一二区在线观| 在线观看三级黄色| 日韩亚洲欧美综合| 中文字幕久久专区| 国内精品美女久久久久久| 少妇人妻一区二区三区视频| 成人漫画全彩无遮挡| 国产精品久久久久久精品电影小说 | 国产69精品久久久久777片| 男人添女人高潮全过程视频| 麻豆久久精品国产亚洲av| www.色视频.com| 嫩草影院新地址| 久久99热6这里只有精品| 国产在线一区二区三区精| 免费看不卡的av| 身体一侧抽搐| 午夜激情久久久久久久| 热99国产精品久久久久久7| 亚洲精品乱码久久久v下载方式| 欧美人与善性xxx| 亚洲精品自拍成人| 国产亚洲最大av| 大香蕉97超碰在线| videos熟女内射| 女人久久www免费人成看片| 亚洲精品,欧美精品| 久久精品综合一区二区三区| 黄色视频在线播放观看不卡| 国产在线男女| 国产极品天堂在线| 精品久久久久久久久av| 免费黄色在线免费观看| 亚洲人成网站在线观看播放| 汤姆久久久久久久影院中文字幕| 舔av片在线| 亚洲欧美一区二区三区国产| 亚洲国产精品999| 亚洲精品乱码久久久久久按摩| 最近最新中文字幕免费大全7| 黄色日韩在线| 精品人妻视频免费看| 亚洲最大成人av| 久久久久精品性色| 美女被艹到高潮喷水动态| 国产高清不卡午夜福利| 国产精品秋霞免费鲁丝片| 超碰av人人做人人爽久久| 99热全是精品| 久久久久久久久久成人| 欧美日韩精品成人综合77777| 午夜亚洲福利在线播放| 亚洲精品国产av蜜桃| 国产毛片a区久久久久| 亚洲自拍偷在线| 欧美日韩视频高清一区二区三区二| 国产精品国产三级国产专区5o| 国产伦精品一区二区三区四那| 久久久久久九九精品二区国产| 日韩av免费高清视频| 久久久色成人| 一级爰片在线观看| 久久久久久九九精品二区国产| 有码 亚洲区| 男插女下体视频免费在线播放| 99热这里只有是精品50| 亚洲精品亚洲一区二区| 欧美日韩视频高清一区二区三区二| 秋霞在线观看毛片| 免费大片黄手机在线观看| 精华霜和精华液先用哪个| 午夜免费鲁丝| av线在线观看网站| 在线a可以看的网站| 精品久久久久久电影网| 亚洲欧美中文字幕日韩二区| 国产精品99久久久久久久久| 日本wwww免费看| 性插视频无遮挡在线免费观看| 亚洲国产成人一精品久久久| 在线天堂最新版资源| 国产精品一及| 丰满少妇做爰视频| av线在线观看网站| av在线蜜桃| 一边亲一边摸免费视频| 麻豆乱淫一区二区| 人妻制服诱惑在线中文字幕| 亚洲一级一片aⅴ在线观看| 亚洲精品国产av成人精品| 老女人水多毛片| 日本免费在线观看一区| 国产亚洲最大av| 美女高潮的动态| 成人高潮视频无遮挡免费网站| 日韩中字成人| 久久韩国三级中文字幕| 精品亚洲乱码少妇综合久久| 日韩av不卡免费在线播放| 激情 狠狠 欧美| 久久99热6这里只有精品| 日韩成人伦理影院| 成人国产麻豆网| 一区二区三区免费毛片| 可以在线观看毛片的网站| 男女啪啪激烈高潮av片| 成人特级av手机在线观看| 久久精品国产亚洲网站| 亚洲无线观看免费| 日本色播在线视频| 大香蕉97超碰在线| 欧美变态另类bdsm刘玥| 亚洲欧洲国产日韩| 好男人在线观看高清免费视频| 秋霞伦理黄片| 只有这里有精品99| 久久韩国三级中文字幕| 在线观看人妻少妇| 寂寞人妻少妇视频99o| 久久人人爽人人爽人人片va| 超碰97精品在线观看| 日韩,欧美,国产一区二区三区| 精品国产三级普通话版| 在线 av 中文字幕| 国产午夜精品一二区理论片| 亚洲精品日韩在线中文字幕| 亚洲国产精品成人综合色| 国产在线一区二区三区精| 欧美xxⅹ黑人| 午夜亚洲福利在线播放| 99热这里只有是精品50| 亚洲欧美精品专区久久| 国模一区二区三区四区视频| 成人美女网站在线观看视频| 免费高清在线观看视频在线观看| 国产综合精华液| 欧美 日韩 精品 国产| 女人久久www免费人成看片| 日韩大片免费观看网站| 日韩一区二区三区影片| 看十八女毛片水多多多| 有码 亚洲区| 成人鲁丝片一二三区免费| 国产毛片a区久久久久| 我的女老师完整版在线观看| 免费黄色在线免费观看| 九九在线视频观看精品| 国产视频首页在线观看| 成人亚洲精品av一区二区| 日本免费在线观看一区| 狂野欧美激情性xxxx在线观看| 美女视频免费永久观看网站| 永久免费av网站大全| 国产精品女同一区二区软件| 极品少妇高潮喷水抽搐| 亚洲欧美精品专区久久| 丰满少妇做爰视频| 欧美+日韩+精品| 99久久精品热视频| 午夜爱爱视频在线播放| 大香蕉久久网| 亚洲精品色激情综合| 麻豆乱淫一区二区| 成年女人看的毛片在线观看| 69人妻影院| av在线app专区| 男的添女的下面高潮视频| 日本av手机在线免费观看| 直男gayav资源| av在线蜜桃| 超碰av人人做人人爽久久| 精品人妻一区二区三区麻豆| 国内少妇人妻偷人精品xxx网站| 在线观看美女被高潮喷水网站| 亚洲美女搞黄在线观看| 91午夜精品亚洲一区二区三区| 免费大片18禁| 成人综合一区亚洲| 日韩伦理黄色片| 毛片女人毛片| 亚洲最大成人中文| 麻豆国产97在线/欧美| 久久精品久久久久久久性| 成人亚洲精品一区在线观看 | 国产 一区精品| av国产精品久久久久影院| tube8黄色片| 亚洲精品视频女| 免费不卡的大黄色大毛片视频在线观看| 亚洲欧美清纯卡通| 亚洲最大成人手机在线| 一区二区三区四区激情视频| 成人鲁丝片一二三区免费| 岛国毛片在线播放| 国产精品.久久久| 国内精品美女久久久久久| 国产精品人妻久久久影院| 大香蕉97超碰在线| 亚洲三级黄色毛片| 肉色欧美久久久久久久蜜桃 | 国产精品嫩草影院av在线观看| 又爽又黄a免费视频| 亚洲第一区二区三区不卡| 高清日韩中文字幕在线| 欧美激情在线99| 亚洲精品中文字幕在线视频 | 在线观看国产h片| 国产成人aa在线观看| 日韩电影二区| 中文天堂在线官网| 国产av不卡久久| 春色校园在线视频观看| 精品人妻一区二区三区麻豆| 夫妻性生交免费视频一级片| 亚洲精品成人久久久久久| 91久久精品国产一区二区三区| 国内精品宾馆在线| 国产综合懂色| 成人亚洲精品一区在线观看 | 欧美三级亚洲精品| 国内揄拍国产精品人妻在线| 亚洲aⅴ乱码一区二区在线播放| 色视频www国产| 18+在线观看网站| 午夜福利高清视频| 一边亲一边摸免费视频| 色婷婷久久久亚洲欧美| 最近最新中文字幕大全电影3| 三级经典国产精品| 在线观看一区二区三区| 搡老乐熟女国产| 午夜激情福利司机影院| 人妻 亚洲 视频| 国产永久视频网站| 一边亲一边摸免费视频| 亚洲欧美中文字幕日韩二区| 亚洲在线观看片| 黄色配什么色好看| 婷婷色综合大香蕉| 18禁在线播放成人免费| 久久久久国产精品人妻一区二区| 男女那种视频在线观看| 国产精品一区二区在线观看99| 26uuu在线亚洲综合色| 午夜亚洲福利在线播放| 男人狂女人下面高潮的视频| 久久精品国产亚洲av涩爱| 亚洲精品aⅴ在线观看| 国产国拍精品亚洲av在线观看| 最近手机中文字幕大全| 亚洲天堂国产精品一区在线| 欧美三级亚洲精品| 国产成年人精品一区二区| 日本欧美国产在线视频| 国产精品99久久久久久久久| 成人亚洲欧美一区二区av| 亚洲熟女精品中文字幕| 国产精品女同一区二区软件| 精品人妻熟女av久视频| 尾随美女入室| 日本与韩国留学比较| 日本猛色少妇xxxxx猛交久久| 毛片女人毛片| 香蕉精品网在线| 成人亚洲精品一区在线观看 | 午夜免费男女啪啪视频观看| 国产熟女欧美一区二区| 日韩不卡一区二区三区视频在线| 色播亚洲综合网| 欧美极品一区二区三区四区| 免费大片黄手机在线观看| 亚洲欧洲国产日韩| 午夜福利视频精品| 亚洲久久久久久中文字幕| 黄色配什么色好看| 精品国产三级普通话版| 天天躁日日操中文字幕| 一级av片app| 可以在线观看毛片的网站| 国内少妇人妻偷人精品xxx网站| 精品午夜福利在线看| 麻豆成人午夜福利视频| 中文字幕久久专区| 日本黄大片高清| 一级毛片aaaaaa免费看小| 亚洲av二区三区四区| 国产精品国产三级国产av玫瑰| 精品久久久精品久久久| 直男gayav资源| 国产精品久久久久久久电影| 国产成人精品福利久久| 天天一区二区日本电影三级| 日本午夜av视频| 国产高清有码在线观看视频| 欧美成人精品欧美一级黄| 亚洲最大成人中文| 精品久久久久久久久亚洲| 美女内射精品一级片tv| 观看免费一级毛片| 在线看a的网站| 亚洲精品乱久久久久久| 男人和女人高潮做爰伦理| 国产综合懂色| av在线蜜桃| 欧美国产精品一级二级三级 | 1000部很黄的大片| 在线观看国产h片| 婷婷色av中文字幕| 伦理电影大哥的女人| 国产精品av视频在线免费观看| videossex国产| 成人一区二区视频在线观看| 久久久午夜欧美精品| 2021天堂中文幕一二区在线观| 国产精品成人在线| 久久人人爽av亚洲精品天堂 | 毛片一级片免费看久久久久| 欧美xxxx黑人xx丫x性爽| 久久精品夜色国产| av专区在线播放| 男的添女的下面高潮视频| 2018国产大陆天天弄谢| 一本一本综合久久| 18禁动态无遮挡网站| 成人高潮视频无遮挡免费网站| 哪个播放器可以免费观看大片| av在线天堂中文字幕| 国产亚洲一区二区精品| 少妇高潮的动态图| 亚洲怡红院男人天堂| 下体分泌物呈黄色| tube8黄色片| 亚洲精品视频女| 亚洲aⅴ乱码一区二区在线播放| 久久久久国产精品人妻一区二区| 黄色一级大片看看| 丰满人妻一区二区三区视频av| 亚洲精品日本国产第一区| 一级爰片在线观看| av在线蜜桃| 亚洲色图综合在线观看| 狠狠精品人妻久久久久久综合| 国产有黄有色有爽视频| 国国产精品蜜臀av免费| 成年免费大片在线观看| 成人免费观看视频高清| 免费不卡的大黄色大毛片视频在线观看| 午夜激情福利司机影院| 国产黄片美女视频| 高清av免费在线| 男人狂女人下面高潮的视频| 插逼视频在线观看| 国产高清国产精品国产三级 | 女人被狂操c到高潮| .国产精品久久| a级毛片免费高清观看在线播放| 又爽又黄无遮挡网站| av天堂中文字幕网| 国产一区二区亚洲精品在线观看| 美女高潮的动态| 99视频精品全部免费 在线| 国产一区二区亚洲精品在线观看| 蜜桃亚洲精品一区二区三区| 高清午夜精品一区二区三区| 亚洲av电影在线观看一区二区三区 | 久久久久久久久久久丰满| 免费观看在线日韩| 91久久精品国产一区二区三区| 别揉我奶头 嗯啊视频| 成人毛片60女人毛片免费| 久久影院123| 又黄又爽又刺激的免费视频.| 亚洲天堂av无毛| 国产黄色视频一区二区在线观看| 国产精品.久久久| 成人特级av手机在线观看| 日日啪夜夜爽| 国产精品国产三级专区第一集| 一级毛片我不卡| 特级一级黄色大片| 免费看不卡的av| 少妇人妻 视频| 国产精品国产三级国产专区5o| 夜夜看夜夜爽夜夜摸| 69av精品久久久久久| 男女那种视频在线观看| 日韩中字成人| 国产精品一区www在线观看| 国产高清不卡午夜福利| 免费大片18禁| 日本一本二区三区精品| 久久精品久久久久久噜噜老黄| 亚洲人与动物交配视频| 精品国产露脸久久av麻豆| 少妇裸体淫交视频免费看高清| 真实男女啪啪啪动态图| 午夜福利网站1000一区二区三区| 欧美激情在线99| 中国美白少妇内射xxxbb| 久久人人爽人人爽人人片va| 又爽又黄a免费视频| 久久99精品国语久久久| 国内精品宾馆在线| 国产成人精品福利久久| 日韩强制内射视频| 老女人水多毛片| 美女主播在线视频| 久久久久国产精品人妻一区二区| 你懂的网址亚洲精品在线观看| 天美传媒精品一区二区| 免费人成在线观看视频色| 一区二区av电影网| 色视频www国产| 国产免费福利视频在线观看| 最近中文字幕高清免费大全6| 精品视频人人做人人爽| 精品国产三级普通话版| 美女高潮的动态| 午夜福利高清视频| 亚洲精品国产av成人精品| 免费看光身美女| 激情五月婷婷亚洲| 亚洲精品中文字幕在线视频 | 日本av手机在线免费观看| 亚洲精品456在线播放app| 在线a可以看的网站| 国产亚洲午夜精品一区二区久久 | 国产精品国产av在线观看| 国产精品偷伦视频观看了| 亚洲国产精品成人综合色| 免费黄频网站在线观看国产| 亚洲欧美精品自产自拍| 亚洲精品国产成人久久av| 女的被弄到高潮叫床怎么办| 高清视频免费观看一区二区| 亚洲欧洲日产国产| 国产亚洲av片在线观看秒播厂| 婷婷色麻豆天堂久久| 亚洲四区av| 不卡视频在线观看欧美| 午夜激情久久久久久久| 午夜激情福利司机影院| 看十八女毛片水多多多| 精品久久国产蜜桃| 欧美极品一区二区三区四区| 国产精品福利在线免费观看| 国语对白做爰xxxⅹ性视频网站| 搞女人的毛片| 黑人高潮一二区| 熟妇人妻不卡中文字幕| 亚洲精品自拍成人| 亚洲经典国产精华液单| 日日撸夜夜添| 麻豆精品久久久久久蜜桃| 亚洲欧美成人精品一区二区| 国产精品人妻久久久影院| 婷婷色av中文字幕| 男人舔奶头视频| 最后的刺客免费高清国语| 色综合色国产| 国产一区亚洲一区在线观看| 亚洲国产精品国产精品| 精品人妻偷拍中文字幕| 国产成年人精品一区二区| 最近的中文字幕免费完整| 欧美最新免费一区二区三区| 久久精品综合一区二区三区| 亚洲人成网站在线播| 亚洲欧美精品专区久久| 在线 av 中文字幕| 午夜福利在线观看免费完整高清在| 一级毛片久久久久久久久女| 国产精品国产av在线观看| 日韩精品有码人妻一区| 在线 av 中文字幕| 搞女人的毛片| 免费av毛片视频| 人妻系列 视频| 久久久久久久久久久免费av| 91在线精品国自产拍蜜月| 熟女av电影| 日韩三级伦理在线观看| 欧美成人午夜免费资源| 亚洲一区二区三区欧美精品 | 天堂网av新在线| 99热网站在线观看| 久久99热6这里只有精品| 美女脱内裤让男人舔精品视频| 欧美xxxx黑人xx丫x性爽| 日日摸夜夜添夜夜添av毛片| 五月伊人婷婷丁香| 日韩 亚洲 欧美在线| 白带黄色成豆腐渣| 青春草亚洲视频在线观看| 欧美变态另类bdsm刘玥| 亚洲av.av天堂| 高清在线视频一区二区三区| 91aial.com中文字幕在线观看| 国产视频内射| 日韩在线高清观看一区二区三区| 国内揄拍国产精品人妻在线| 在线 av 中文字幕| 成人免费观看视频高清| 欧美日韩视频高清一区二区三区二| 国产精品国产av在线观看| 国产精品人妻久久久影院| 真实男女啪啪啪动态图| 欧美日韩视频精品一区| 国产精品久久久久久精品古装| 人体艺术视频欧美日本| 毛片一级片免费看久久久久| 国产伦精品一区二区三区四那| 国产国拍精品亚洲av在线观看| 直男gayav资源| 九九久久精品国产亚洲av麻豆| 免费黄色在线免费观看| 国产永久视频网站| 插阴视频在线观看视频| 欧美区成人在线视频| 午夜亚洲福利在线播放| 午夜精品国产一区二区电影 | 成人亚洲精品一区在线观看 | 色视频www国产| 国产极品天堂在线| 99热这里只有是精品在线观看| 日本黄大片高清| 久久国内精品自在自线图片| 人人妻人人澡人人爽人人夜夜| 99九九线精品视频在线观看视频| 91狼人影院| 搡老乐熟女国产| 亚洲人成网站在线观看播放| 久久精品国产自在天天线| a级毛色黄片| 国产 一区精品| 久久久久精品性色| 亚洲第一区二区三区不卡| 深爱激情五月婷婷| 黄色欧美视频在线观看| 久久精品久久久久久噜噜老黄| 2021少妇久久久久久久久久久| 尤物成人国产欧美一区二区三区| 街头女战士在线观看网站| 婷婷色综合www| 黄色日韩在线| 亚洲国产最新在线播放| 伦理电影大哥的女人| 色播亚洲综合网| 免费av观看视频| 欧美一区二区亚洲| 亚洲在线观看片| 亚洲av福利一区| 又爽又黄a免费视频| 免费电影在线观看免费观看| 熟妇人妻不卡中文字幕| 亚洲怡红院男人天堂| 嘟嘟电影网在线观看| 欧美日韩视频精品一区| 免费在线观看成人毛片| 国产精品偷伦视频观看了| 精品人妻熟女av久视频| freevideosex欧美| 熟女人妻精品中文字幕| 亚洲精品日韩在线中文字幕| 2021少妇久久久久久久久久久| 十八禁网站网址无遮挡 | 国产精品一及| 国产高清三级在线| 亚洲欧美成人综合另类久久久| 国产在线一区二区三区精| 97人妻精品一区二区三区麻豆| 美女视频免费永久观看网站| 2022亚洲国产成人精品| 天堂俺去俺来也www色官网| 欧美精品人与动牲交sv欧美| 51国产日韩欧美| 午夜福利在线观看免费完整高清在| 国内精品宾馆在线| 国产有黄有色有爽视频| 高清欧美精品videossex| 亚洲av二区三区四区| 午夜福利视频1000在线观看| 中文资源天堂在线| av在线天堂中文字幕| 午夜免费男女啪啪视频观看| av在线天堂中文字幕| 国产色爽女视频免费观看| 精品人妻视频免费看| 韩国高清视频一区二区三区| 亚洲欧美精品专区久久| 久久人人爽人人片av| 免费观看a级毛片全部| 日韩一区二区三区影片| 国产一区二区三区av在线| 亚洲精品成人av观看孕妇| 日本与韩国留学比较| 亚洲欧美日韩另类电影网站 | 精品99又大又爽又粗少妇毛片| 看非洲黑人一级黄片| 国产乱来视频区| 黄片wwwwww| 亚洲精品中文字幕在线视频 | 国产视频首页在线观看| 久久精品国产亚洲av涩爱| 美女内射精品一级片tv| 1000部很黄的大片|