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

    基于MaxEnt模型預(yù)測(cè)極北海帶在我國(guó)黃渤海的適生情況*

    2019-09-27 01:42:14劉福利梁洲瑞汪文俊孫修濤
    漁業(yè)科學(xué)進(jìn)展 2019年5期
    關(guān)鍵詞:黃渤海適生區(qū)海帶

    孫 昕 劉福利 梁洲瑞 汪文俊 孫修濤

    基于MaxEnt模型預(yù)測(cè)極北海帶在我國(guó)黃渤海的適生情況*

    孫 昕 劉福利①梁洲瑞 汪文俊 孫修濤

    (中國(guó)水產(chǎn)科學(xué)研究院黃海水產(chǎn)研究所 農(nóng)業(yè)農(nóng)村部海洋漁業(yè)可持續(xù)發(fā)展重點(diǎn)實(shí)驗(yàn)室 青島 266071)

    我國(guó)黃渤海區(qū)域的自然藻場(chǎng)存在不同程度的退化,建立人工藻場(chǎng)是制止海底荒漠化、恢復(fù)海底植被和建設(shè)海洋牧場(chǎng)的重要手段。本研究模擬極北海帶()在我國(guó)黃渤海區(qū)域的適生情況,使用MaxEnt生態(tài)位預(yù)測(cè)方法,針對(duì)影響極北海帶分部的21項(xiàng)關(guān)鍵水質(zhì)及水文條件以及在原產(chǎn)地的分布情況,建立我國(guó)黃渤海海域生態(tài)位分布模型,劃分相應(yīng)的適生區(qū)域,并通過(guò)刀切法篩選對(duì)模型貢獻(xiàn)度較大的環(huán)境因子,逐一分析其對(duì)極北海帶適生性的影響。結(jié)果顯示,在篩選出的9項(xiàng)對(duì)模型貢獻(xiàn)度較大的環(huán)境因子中,溫度與光照強(qiáng)度對(duì)極北海帶分布的影響最大。極北海帶在溫度7℃~18℃范圍、光強(qiáng)不超過(guò)52 E/(m2?d)的區(qū)間有較大的適生概率;氮、磷濃度對(duì)預(yù)測(cè)結(jié)果也有一定影響,但天然海水中的氮、磷濃度不足以成為極北海帶分布的限制因子??傮w來(lái)看,黃渤海區(qū)域有1.67%海域?yàn)橹羞m生區(qū),主要集中在遼東半島東岸長(zhǎng)??h附近海域;山東半島東部、北部,渤海西北部有一定的低適生范圍,占總海域面積的5.12%;非適生區(qū)和邊緣適生區(qū)占海域總面積的88.51%和4.70%。本研究表明,遼東半島東岸長(zhǎng)海縣附近海域是建立極北海帶藻場(chǎng)的理想?yún)^(qū)域。

    極北海帶;MaxEnt;人工藻場(chǎng);適生范圍;生態(tài)位

    海洋大型藻類是海洋中重要的初級(jí)生產(chǎn)者,同時(shí)也是多種魚(yú)類、貝類、底棲生物的重要棲息地,為其提供必不可少的遮蔽環(huán)境、食物來(lái)源、棲息場(chǎng)所、繁殖區(qū)域等環(huán)境條件。近年來(lái),圍繞山東半島及遼東半島藻場(chǎng)修復(fù)的研究有了很大進(jìn)展,在受損海區(qū)生態(tài)環(huán)境功能修復(fù)上取得一定成效(涂忠等, 2007、2009),并發(fā)揮出相應(yīng)生態(tài)功能(Largo, 1994; Seaman, 2007; Jiao, 2014)。然而,海藻場(chǎng)的建設(shè)更多看重其生態(tài)功能,以吸引魚(yú)蝦生物及浮游生物作為主要目的(張磊等, 2012; 袁小楠等, 2017),若能將富有工業(yè)價(jià)值的海洋大型藻類,如海帶屬(),作為主導(dǎo)藻種構(gòu)建海藻床,在起到生態(tài)恢復(fù)功能的同時(shí),也可以綜合提高藻場(chǎng)的經(jīng)濟(jì)產(chǎn)出率。極北海帶()系多年生大型海藻,源于大西洋東北海域,是當(dāng)?shù)亍昂T迳帧钡膬?yōu)勢(shì)藻種(Steneck, 2002)。極北海帶具有較高經(jīng)濟(jì)價(jià)值,其柄和葉片的褐藻膠含量分別為25%~30%和17%~33%(干重比),屬于褐藻膠含量較高的褐藻之一,且其褐藻膠凝膠強(qiáng)度大,因古羅糖醛酸(Guluronic acid, G)含量在其柄和葉片的褐藻膠中分別占比高達(dá)68%~70%和44%~55%,而在其他商用褐藻膠中的含量一般為10%~45% (Peter, 2015),工業(yè)和商業(yè)價(jià)值顯著。

    在人工藻場(chǎng)的構(gòu)建中,由于需要考慮生態(tài)恢復(fù)效果,提高底棲生物豐富度,海藻附著基通常設(shè)置在水深為15~25 m的硬質(zhì)海底(許強(qiáng)等, 2011)。然而,我國(guó)常見(jiàn)的海帶()通常的垂直適生深度小于15 m (曾呈奎等, 1981),在人工魚(yú)礁深度的海水層中,光照強(qiáng)度往往無(wú)法達(dá)到海帶正常生長(zhǎng)所需的光強(qiáng),使海帶生長(zhǎng)較慢,甚至死亡,導(dǎo)致生態(tài)功能恢復(fù)效果不理想。海帶的適生深度主要受光照強(qiáng)度的影響,自然生長(zhǎng)的極北海帶對(duì)光照的需求較低,因此,在更深的海域也有分布。Joanna(1997)研究發(fā)現(xiàn),在25 m以下深度的海底存在自然生長(zhǎng)的極北海帶。Luning(1979)對(duì)極北海帶進(jìn)行光照實(shí)驗(yàn)發(fā)現(xiàn),極北海帶生長(zhǎng)所需光照顯著低于其他海帶屬種類(,)。因此,針對(duì)藻礁建設(shè)區(qū)域深度較深、光照不足的問(wèn)題,極北海帶可能是一個(gè)理想的適應(yīng)品種。

    Phillips等(2008)基于生態(tài)位理論,用最大熵原理作為統(tǒng)計(jì)推斷工具,以物種在已知分布區(qū)的信息以及目標(biāo)區(qū)的環(huán)境變量為基礎(chǔ),通過(guò)比較該物種在已知分布區(qū)的生態(tài)環(huán)境變量來(lái)預(yù)測(cè)其占有的生態(tài)位,并編寫程序MaxEnt(Maxent entropy)。MaxEnt僅考慮目標(biāo)物種在所在區(qū)域的“出現(xiàn)”情況(Presence-only),綜合考慮各種給定環(huán)境因子下物種出現(xiàn)的最大熵,從而對(duì)物種在預(yù)測(cè)區(qū)域的適生情況進(jìn)行分析。該模型在預(yù)測(cè)物種適生性上有廣泛應(yīng)用,許多學(xué)者使用該模型對(duì)珍稀物種、入侵物種以及一些植物病蟲(chóng)害的潛在分布區(qū)域進(jìn)行預(yù)測(cè)評(píng)價(jià)(曾輝等, 2008; 張飛等, 2017; 王運(yùn)生等, 2007)。此外,該模型也被證實(shí)適合預(yù)測(cè)海洋藻類的適生范圍(Martín- García, 2017)。

    盡管是構(gòu)建海藻場(chǎng)的優(yōu)良藻種,極具開(kāi)發(fā)潛力,但作為外來(lái)物種,極北海帶在我國(guó)黃渤海海區(qū)的適應(yīng)性有待研究,是否會(huì)成為入侵物種的風(fēng)險(xiǎn)也尚不明確。因此,本研究采用MaxEnt軟件對(duì)極北海帶在我國(guó)黃渤海海區(qū)進(jìn)行適生范圍預(yù)測(cè),評(píng)價(jià)其生態(tài)入侵風(fēng)險(xiǎn),為構(gòu)建極北海帶藻場(chǎng)的可行性提供參考依據(jù)。

    1 模型與方法

    1.1 生態(tài)位預(yù)測(cè)模型

    本研究選用MaxEnt v3.3.3對(duì)我國(guó)近海海域模擬極北海帶的適生性進(jìn)行評(píng)估,選取logistic形式表達(dá),其他算法選擇模型自身默認(rèn)值,并將所選用的環(huán)境因子做刀切法(JackKnife)檢驗(yàn),判斷每個(gè)環(huán)境因子對(duì)預(yù)測(cè)結(jié)果的重要性。

    1.2 極北海帶空間位點(diǎn)

    本研究使用的極北海帶空間數(shù)據(jù)從全球生物多樣性信息網(wǎng)絡(luò)GBIF (http://www.gbif.org/)中獲取。收集1900~2018年極北海帶出現(xiàn)位點(diǎn)共11894處。除去實(shí)驗(yàn)、標(biāo)本存儲(chǔ)及其他人為培養(yǎng)位點(diǎn)822處,得到極北海帶自然擴(kuò)散位點(diǎn)11072處的經(jīng)緯度數(shù)據(jù)??紤]到模型運(yùn)算量以及許多點(diǎn)位較為集中,使用Excel在11072個(gè)點(diǎn)位中隨機(jī)選擇300個(gè)位點(diǎn)作為極北海帶的空間位點(diǎn)數(shù)據(jù),將300個(gè)空間位點(diǎn)數(shù)據(jù)轉(zhuǎn)換成*.csv格式導(dǎo)入MaxEnt模型。

    1.3 環(huán)境因子數(shù)據(jù)

    于全球海洋生物擴(kuò)散模型環(huán)境數(shù)據(jù)庫(kù)Bio- ORICLE(http://bio-oracle.org/)中篩選出影響極北海帶生長(zhǎng)的21項(xiàng)關(guān)鍵環(huán)境因子,用于適生模型構(gòu)建。其中包含12項(xiàng)水質(zhì)參數(shù),即:最高月均溫度、溶解氧、氮濃度、磷濃度;最低月均溫度、溶解氧、氮濃度、磷濃度,年平均溫度、溶解氧、氮濃度、磷濃度,以及9項(xiàng)水文參數(shù),即:最大海流流速、海區(qū)深度、光照強(qiáng)度;最小海流流速、海區(qū)深度、有效太陽(yáng)輻射,年平均海流流速、海區(qū)深度、光照強(qiáng)度。

    1.4 環(huán)境因子選擇

    本研究采用MaxEnt模型中自帶的刀切法,將上述21項(xiàng)因子進(jìn)行篩選。將全球生物多樣性信息網(wǎng)絡(luò)中所收集到的極北海帶空間位點(diǎn)經(jīng)緯度數(shù)據(jù)轉(zhuǎn)化成數(shù)據(jù)文本格式(.csv),并于Bio-ORICLE中所收集的環(huán)境數(shù)據(jù)文件(.asc)同時(shí)導(dǎo)入MaxEnt模型,以logistic格式輸出最終結(jié)果。模型共計(jì)算2次,第1次計(jì)算用于篩選對(duì)適生模型貢獻(xiàn)度為0的環(huán)境因子,舍去這些因子,將剩余環(huán)境因子再次導(dǎo)入MaxEnt,進(jìn)行第2次預(yù)測(cè)模型構(gòu)建。

    1.5 適生區(qū)評(píng)估

    模型預(yù)測(cè)結(jié)果以*.asc格式輸出。將模型輸出結(jié)果導(dǎo)入ArcGIS 10.3,利用ArcGIS工具箱系統(tǒng)中的格式轉(zhuǎn)換工具( conversion tools),將預(yù)測(cè)結(jié)果數(shù)據(jù)轉(zhuǎn)化為柵格數(shù)據(jù)(*.flt)格式。選定合適的閾值,將預(yù)測(cè)結(jié)果按適生概率等級(jí)劃分為高適生區(qū)(>0.5)、中適生區(qū)(0.3~0.5)、低適生區(qū)(0.1~0.3)、邊緣適生區(qū)(0.05~0.1)、非適生區(qū)(<0.05)(Martín-García, 2014; 齊國(guó)軍等, 2012; 柯明劍等, 2017)。

    1.6 MaxEnt模型評(píng)估

    本研究通過(guò)計(jì)算指受試者工作特征曲線(ROC)的曲線下面積(AUC)來(lái)評(píng)估模型準(zhǔn)確性(Canran, 2005)。在進(jìn)行模擬時(shí),使用模型自身功能將所選取的300個(gè)位點(diǎn)分為2部分,隨機(jī)選取75%的物種分布點(diǎn)用于建模,剩余25%用于模型的檢驗(yàn)和校正(Phillips, 2006、2008)。在“presence-only”模型中,ROC曲線表示正確預(yù)測(cè)的出現(xiàn)位點(diǎn)概率與假信號(hào)出現(xiàn)概率的比值,ROC曲線下的面積(AUC)即表示預(yù)測(cè)結(jié)果的準(zhǔn)確性。AUC數(shù)值范圍為0~1,AUC數(shù)值越大,表示預(yù)測(cè)結(jié)果越準(zhǔn)確,AUC值大于0.9時(shí),表示模型預(yù)測(cè)結(jié)果有很高的可信度。

    2 結(jié)果

    黃渤海海區(qū)統(tǒng)計(jì)共5213個(gè)柵格數(shù)據(jù),最大值為0.472,最小值為0.000(精確到小數(shù)點(diǎn)后3位)。將適生概率按適生區(qū)域閾值劃分,頻數(shù)與頻率見(jiàn)表1。

    極北海帶在中國(guó)黃渤海海域的適生性預(yù)測(cè)結(jié)果如圖1所示,適生概率從藍(lán)色到紅色依次上升。極北海帶在山東半島沿岸、渤海灣近岸海域以及遼東半島沿岸均有一定的適生概率,低適生區(qū)多分布在渤海西北部及遼東半島沿岸。山東半島南岸有少量邊緣適生區(qū),再往南均無(wú)適生區(qū)。中、高適生區(qū)域主要分布在朝鮮及日本沿海地區(qū),我國(guó)海域無(wú)分布。

    評(píng)估預(yù)測(cè)模型準(zhǔn)確性的ROC曲線如圖2所示,其中,紅色為用于模型檢驗(yàn)校正,藍(lán)色為模型建立。其AUC值分別為0.966與0.965。

    表1 極北海帶在我國(guó)黃渤海海域的適生區(qū)分布情況

    Tab.1 The occurrence zones of L. hyperborean in Bohai and Huanghai Sea

    注:適生頻數(shù)及概率范圍通過(guò)ArcGIS中的工具箱系統(tǒng)中的格式轉(zhuǎn)換工具重新劃分

    Note: Occurrence probability zones (OPZs) and frequency are reclassified by the use of conversion tool from ArcGIS

    圖1 MaxEnt模型預(yù)測(cè)極北海帶在中國(guó)黃渤海海域的適生性概率范圍

    采用MaxEnt模型中自帶的刀切法功能進(jìn)行環(huán)境變量貢獻(xiàn)度分析。第1次分析后,從所選21項(xiàng)環(huán)境變量中,將貢獻(xiàn)度為0的12項(xiàng)變量刪除,保留剩余9項(xiàng)環(huán)境變量。SPSS檢測(cè)9項(xiàng)變量?jī)蓛芍g無(wú)顯著相關(guān)性。圖3為9項(xiàng)有貢獻(xiàn)度的環(huán)境變量刀切法分析結(jié)果。納入考慮的9個(gè)環(huán)境變量分別為最大月光照強(qiáng)度Light(max)、年均光照強(qiáng)度Light(mean)、平均海流流速Current(mean)、年均溶解氧濃度DO(mean)、最低月均氮濃度N(min)、最低月均磷濃度P(min)、最大月均溫度(max)、年均溫度(mean)以及最小月均溫度(min)。光照強(qiáng)度和溫度對(duì)極北海帶適生概率值的貢獻(xiàn)度最為顯著。

    選取MaxEnt模型篩選出的對(duì)極北海帶適生性影響較大的環(huán)境因子逐個(gè)分析(圖4),僅考慮單一環(huán)境變量時(shí)的極北海帶適生概率,用來(lái)反映模型在預(yù)測(cè)適生性時(shí),對(duì)所選變量的依賴性以及所選變量和其他變量之間的相關(guān)性引起的依賴性。由圖4可見(jiàn),極北海帶適生概率隨光照強(qiáng)度的增加而增加,在最大光照強(qiáng)度52 E/(m2?d) (Photosynt. Avail. Radiation)、平均強(qiáng)度48 E/(m2?d)時(shí)達(dá)到頂點(diǎn)。極北海帶適生溫度范圍為7℃~18℃、在年均溫度為13℃時(shí)達(dá)到最大適生概率。同時(shí),極北海帶的適生溫度范圍最小月均溫度不低于–2℃,最大月均溫度不超過(guò)25℃。氮、磷濃度的單因子適生概率分布圖顯示,極北海帶在0~10 mg/L氮濃度范圍、0~0.22 mg/L磷濃度范圍內(nèi)存在適生性。

    圖2 MaxEnt預(yù)測(cè)模型的ROC曲線

    圖3 環(huán)境因子對(duì)模型的貢獻(xiàn)度

    圖4 單因子環(huán)境變量分析

    3 討論

    本研究通過(guò)MaxEnt建立的極北海帶適生性預(yù)測(cè)模型,訓(xùn)練集和測(cè)試集的AUC值均非常接近于1,說(shuō)明結(jié)果預(yù)測(cè)準(zhǔn)確(Phillips, 2006)。本研究通過(guò)刀切法篩選出對(duì)極北海帶適生性貢獻(xiàn)最大的9項(xiàng)環(huán)境因子也均符合極北海帶的生長(zhǎng)條件(Susan, 2011)。MaxEnt模型預(yù)測(cè)結(jié)果顯示,除朝鮮西岸與日本主島東部沿岸以外,在我國(guó)范圍內(nèi)僅遼東半島東部大連長(zhǎng)??h附近海域存在中適生區(qū),其余地區(qū)均為低適生區(qū)、邊緣適生區(qū)和非適生區(qū)。

    溫度是影響極北海帶適生性的重要因子。本研究結(jié)果顯示,極北海帶適生溫度年均值在6~17℃范圍內(nèi),月均溫度不超過(guò)25℃,最低不低于–2℃。大連長(zhǎng)海縣附近海域年均水溫通常為–2~24℃,平均溫度和最低溫度在極北海帶適生溫度范圍之內(nèi),但夏季月均最高溫度接近于極北海帶的最高耐受溫度范圍。因此,在夏季,溫度有可能是在該海區(qū)內(nèi)影響極北海帶最重要的環(huán)境因子,在構(gòu)建海帶場(chǎng)時(shí)應(yīng)重點(diǎn)關(guān)注。

    光照強(qiáng)度是另一項(xiàng)影響極北海帶分部主要的環(huán)境因子之一。模型自身的刀切法檢測(cè)將最小光照強(qiáng)度對(duì)極北海帶適生性貢獻(xiàn)排除。該結(jié)果與預(yù)測(cè)結(jié)果相符,說(shuō)明極北海帶對(duì)低光照環(huán)境的適應(yīng)性較好,相比與其他環(huán)境因子,自然的低光照環(huán)境對(duì)極北海帶可能不構(gòu)成影響。

    氮、磷濃度對(duì)極北海帶的適生范圍也有一定的影響,單因子適生性分析結(jié)果顯示,極北海帶在氮、磷濃度極低的水平下均能達(dá)到較高的適生性,表明極北海帶適生范圍受氮、磷濃度影響相對(duì)較小。此外,最大及平均氮、磷濃度對(duì)該預(yù)測(cè)模型沒(méi)有貢獻(xiàn)值,也說(shuō)明自然水體能達(dá)到的最高氮、磷濃度不足以成為極北海帶生長(zhǎng)的限制因子。

    海洋底質(zhì)對(duì)極北海帶分布的影響沒(méi)有納入考慮范圍。盡管海洋底質(zhì)對(duì)大型海藻的自然分布影響很大,海帶生長(zhǎng)通常需要礁巖作為附著基。但藻場(chǎng)的建設(shè)通常包含人工投放礁石和海藻附著基,人為為海帶生長(zhǎng)提供適宜的海洋底質(zhì)環(huán)境,或直接采用苗繩方法進(jìn)行培養(yǎng)。因此,將底質(zhì)類型排除在模型所考慮的環(huán)境因子之外。

    引入外來(lái)物種需要嚴(yán)格把控其適生范圍,防止其成為入侵物種。遼東半島東部海域底質(zhì)以沙質(zhì)為主,不利于極北海帶的自然生長(zhǎng)。因此,在沙質(zhì)底質(zhì)建立引入的極北海帶品種還可降低其成為入侵物種的風(fēng)險(xiǎn),通過(guò)控制海洋牧場(chǎng)附著基質(zhì)的投放區(qū)域,起到控制極北海帶生長(zhǎng)范圍的作用。另一方面,適生概率超過(guò)0.5的高適生區(qū)有較高的入侵風(fēng)險(xiǎn)(Phillips, 2006),而在我國(guó)近海范圍內(nèi)均無(wú)高適生區(qū)。因此,極北海帶在我國(guó)海域內(nèi)成為入侵物種的風(fēng)險(xiǎn)較低。需要注意的是,遼東灣西部和山東半島北部海域底質(zhì)以巖石為主,極北海帶在此有一定的適生概率。然而,MaxEnt預(yù)測(cè)此海域?qū)儆诘瓦m生區(qū),成為入侵物種的概率較低。其余低適生區(qū)例如渤海灣西岸,其底質(zhì)以淤泥為主,不適合極北海帶生長(zhǎng)。

    綜上所述,引進(jìn)極北海帶作為海洋牧場(chǎng)的主導(dǎo)藻種在我國(guó)遼東半島東岸的長(zhǎng)??h附近海域有較高的可行性。同時(shí),受到環(huán)境因子的限制,極北海帶成為入侵物種的風(fēng)險(xiǎn)較低,但在遼東灣西部和山東半島北部海域仍需要嚴(yán)格審查。

    Canran L, Pam MB, Terence PD,. Selecting thresholds of occurrence in the prediction of species distributions. Ecography, 2005, 28(3): 385–393

    Holdt SL, Kraan S. Bioactive compounds in seaweed: Functional food applications and legislation. Journal of Applied Phycology, 2011, 23(3): 543–597

    Jiao L, Chang G, Yong C,. Progress on research and construction of marine ranching along the coast of Shandong Province of China. Journal of Shipping and Ocean Engineering, 2014, 4: 338–347

    Joanna MK. The biology ofX. the effect of depth on some populations. Journal of the Marine Biological Association of the United Kingdom, 1997, 57(3): 587–607

    Ke MJ, Wang HB, Yang Z,. Prediction of the potential distribution areas of brown-banded cockroach in China based on MaxEnt model. Port Health Control, 2017, 22(2): 59–63 [柯明劍, 汪海波, 楊澤, 等. 用最大熵生態(tài)位模型預(yù)測(cè)入侵媒介生物褐帶皮蠊在中國(guó)的潛在分布區(qū). 口岸衛(wèi)生控制, 2017, 22(2): 59–63]

    Largo DB, Ohno M, Critchley AT. Seasonal changes in the growth and reproduction ofC. Ag. andJ. Ag. (Sargassaceae, Fucales) from Liloan, Cebu, in Central Philippines. Sorui, 1994, 42(1): 53–61

    Luning K. Growth strategies of threespecies (Phaeophyceae) inhabiting different depth zones in the sublittoral region of Helgoland (North Sea). Marine Ecology, 1979(1): 195–207

    Martín-García, L, Herrera R, Moro-Abad L,. Predicting the potential habitat of the harmful cyanobacteria: Ecological niche model of invasive seaweeds. Harmful Algae, 2014, 34: 76–86

    Peter S, Kenneth DB, Michele SS,. The seasonal variation in the chemical composition of the kelp species,,and. Journal of Applied Phycology, 2015, 27(1): 363–373

    Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecological Modelling, 2006, 190(3–4): 231–259

    Phillips SJ, Dudík M. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography, 2008, 31(2): 161–175

    Qi GJ, Gao Y, Huang DC,. Historical invasion, expansion process and the potential geographic distributions for the rice water weevil,in China based on MaxEnt. Acta Phytophylacica Sinica, 2012, 39(2): 129–136 [齊國(guó)君, 高燕, 黃德超, 等. 基于MaxEnt的稻水象甲在中國(guó)的入侵?jǐn)U散動(dòng)態(tài)及適生性分析. 植物保護(hù)學(xué)報(bào), 2012, 39(2): 129–136]

    Seaman W. Artificial habitats and the restoration of degraded marine ecosystems and fisheries. Hydrobiologia, 2007, 580(1): 143–155

    Steneck RS, Graham BJ, Bourque D,. Kelp forest ecosystems: Biodiversity, stability, resilience and future. Environment Conservation, 2002, 29: 436–459

    Tu Z, Zhang XM, Zhang PD,. Artificial fish reef near Lidao Island in Rongcheng—Early assessment on resource and environment in sea area for artificial fish reef. Modern Fisheries Information, 2007, 22(12): 3–8 [涂忠, 張秀梅, 張沛東, 等. 榮成俚島人工魚(yú)礁建設(shè)項(xiàng)目——人工魚(yú)礁投放海域資源環(huán)境狀況前期評(píng)價(jià). 現(xiàn)代漁業(yè)信息, 2007, 22(12): 3–8]

    Tu Z, Zhang XM, Zhang PD, Huang GQ,. Artificial reef project in Lidao, Rongcheng Ⅱ. Assessment on the effect of biological resource enhancement by the artificial reef. Modern Fisheries Information, 2009, 24(11): 16–20 [涂忠, 張秀梅, 張沛東, 等. 榮成俚島人工魚(yú)礁建設(shè)項(xiàng)目Ⅱ. 人工魚(yú)礁養(yǎng)護(hù)生物資源效果的評(píng)價(jià). 現(xiàn)代漁業(yè)信息, 2009, 24(11): 16–20]

    Wang YS, Xie BY, Wan FH,. Application of ROC curve analysis in evaluating the performance of alien species’ potential distribution models. Biodiversity Science, 2007, 15(4): 365–372 [王運(yùn)生, 謝丙炎, 萬(wàn)方浩, 等. ROC曲線分析在評(píng)價(jià)入侵物種分布模型中的應(yīng)用. 生物多樣性雜志, 2007, 15(4): 365–372]

    Xu Q, Liu SB, Xu M,. Preliminary study on sites selection of marine ranching—taking Zhoushan as an example. Fishery Modernization, 2011, 38(2): 27–31 [許強(qiáng), 劉舜斌, 許敏, 等. 海洋牧場(chǎng)建設(shè)選址的初步研究: 以舟山為例. 漁業(yè)現(xiàn)代化, 2011, 38(2): 27–31]

    Yuan XN, Liang ZL, Lü ZB,. The effects of ecological restoration of the artificial reef in Weihai offshore. Acta Oceanologica Sinica, 2017, 39(10): 54–64 [袁小楠, 梁振林, 呂振波, 等. 威海近岸人工魚(yú)礁布設(shè)對(duì)生物資源恢復(fù)效果. 海洋學(xué)報(bào), 2017, 39(10): 54–64]

    Zeng CK, Pan ZZ, Zhou BC. Comparative photosynthetic studies on benthic seaweeds Ⅱ. The effect of light intensity on photosynthesis of intertidal brown algae. Chinese Journal of Oceanology and Limnology, 1981, 12(3): 254–257 [曾呈奎, 潘忠正, 周百成. 底棲海藻比較光合作用研究Ⅱ. 潮間帶褐藻的光合作用與光強(qiáng)的關(guān)系. 海洋與湖沼, 1981, 12(3): 254–257]

    Zeng H, Huang GS, Lin W,. Prediction of potential geographic distribution ofin the world using MaxEnt. Plant Protection, 2008, 34(3): 88–92 [曾輝, 黃冠勝, 林偉, 等. 利用MaxEnt預(yù)測(cè)橡膠南美葉疫病菌在全球的潛在地理分布. 植物保護(hù), 2008, 34(3): 88–92]

    Zhang F, Chen SQ, Wang LL,. Study on ecological suitability regionalization ofbased on MaxEnt and ArcGIS model. China Journal of Chinese Materia Medica, 2017, 42(16): 3078–3083 [張飛, 陳隨清, 王利麗, 等. 基于MaxEnt和ArcGIS的山茱萸生態(tài)適宜性區(qū)劃研究. 中國(guó)中藥雜志, 2017, 42(16): 3078–3083]

    Zhang L, Zhang XM, Wu ZX,Effect of environment on benthic macro-algal communities of artificial reefs in Lidao, Rongcheng. Journal of Fishery Sciences of China, 2012, 19(1): 116–125 [張磊, 張秀梅, 吳忠鑫, 等. 榮成俚島人工魚(yú)礁區(qū)大型底棲藻類群落及其與環(huán)境因子的關(guān)系. 中國(guó)水產(chǎn)科學(xué), 2012, 19(1): 116–125]

    Potential Geographic Distribution ofin the Bohaiand Huanghai Seas of China Based on MaxEnt Model

    SUN Xin, LIU Fuli①, LIANG Zhourui, WANG Wenjun, SUN Xiutao

    (Key Laboratory of Sustainable Development of Marine Fisheries, Ministry of Agriculture and Rural Affairs, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071)

    Ecological niche model MaxEnt was applied to estimate the suitability probability ofin the Bohai and Huanghai Seas of China. In this study, the native occurrence data ofand 21 key environmental variables were selected to build the ecological niche model and then redefined into occurrence zones by suitability probabilities. The contribution of each environmental variable to the model was estimated using JackKnifing techniques. The results showed that the contributions from temperature and light intensity to the spreading ofdominated all considered environmental parameters. In the forecasted areas, suitable probabilities were significant when the temperatures range from 7~18℃ and the light intensities were below 52 E/(m2?d). Although the concentrations of N and P contributed to the prediction model, they would barely reach the limitation value forin natural sea water. Additionally, our results showed that 1.67% of the entire forecasted area had a medium distribution probability, which was mainly limited within Changhai Town of Dalian City. Of the total area, 5.12% was low occurrence probability zones, mainly spreading along the east to the Liaodong Peninsula, north to the Shandong Peninsula, and northwest of the Bohai Sea. Furthermore, unbefitting and marginal suitable distribution areas represent 88.51% and 4.70% of Chinese coast, respectively. The results indicate that Changhai Town could be an ideal location for the construction ofbeds.

    ; MaxEnt; Artificial algae field; Suitable scope; Ecological Niche

    LIU Fuli, Email: liufl@ysfri.ac.cn

    *中國(guó)水產(chǎn)科學(xué)研究院基本科研業(yè)務(wù)費(fèi)(2016HY-JC0301)、青島市民生科技計(jì)劃項(xiàng)目(17-3-3-65-nsh)和現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系專項(xiàng)(CARS-50)共同資助 [This work was supported by Special Scientific Research Funds for Central Non-Profit Institutes, Chinese Academy of Fishery Sciences (2016HY-JC0301), Qingdao Science and Technology Project of People′s Livelihood (17-3-3-65-nsh), and China Agriculture Research System (CARS-50)]. 孫 昕, E-mail: 307368070@qq.com

    劉福利,副研究員,E-mail: liufl@ysfri.ac.cn

    2018-05-24,

    2018-08-07

    Q89

    A

    2095-9869(2019)05-0071-07

    10.19663/j.issn2095-9869.20180524001

    http://www.yykxjz.cn/

    孫昕, 劉福利, 梁洲瑞, 汪文俊, 孫修濤. 基于MaxEnt模型預(yù)測(cè)極北海帶在我國(guó)黃渤海的適生情況. 漁業(yè)科學(xué)進(jìn)展, 2019, 40(5): 71–77

    Sun X, Liu FL, Liang ZR, Wang WJ, Sun XT. Potential geographic distribution ofin the Bohai and Huanghai Seas of China based on MaxEnt model. Progress in Fishery Sciences, 2019, 40(5): 71–77

    (編輯 馮小花)

    猜你喜歡
    黃渤海適生區(qū)海帶
    即興
    飛碟探索(2023年3期)2023-11-23 16:08:36
    氣候變化下中國(guó)蒟蒻薯科箭根薯的地理分布格局預(yù)測(cè)
    黃渤海不同水層中浮游植物對(duì)灰霾添加的響應(yīng)
    未來(lái)氣候條件下當(dāng)歸適生區(qū)預(yù)測(cè)及時(shí)空變化分析
    如果沒(méi)有你
    飛碟探索(2023年1期)2023-05-30 22:59:34
    夢(mèng)的解析
    氣候變化下瀕危植物半日花在中國(guó)的潛在分布
    巴拉圭瓜多竹適生區(qū)分布研究
    出海
    波多野结衣av一区二区av| 亚洲无线在线观看| 国产成年人精品一区二区| 成人三级黄色视频| 久久精品国产清高在天天线| 亚洲av美国av| 色播亚洲综合网| 大香蕉久久成人网| 成熟少妇高潮喷水视频| 在线国产一区二区在线| 国产欧美日韩一区二区三区在线| 一二三四社区在线视频社区8| 欧美日韩亚洲综合一区二区三区_| 老鸭窝网址在线观看| 国产麻豆69| 久久影院123| 亚洲自偷自拍图片 自拍| 成人特级黄色片久久久久久久| 欧美av亚洲av综合av国产av| 淫秽高清视频在线观看| 免费在线观看黄色视频的| 国产精品永久免费网站| 亚洲国产精品合色在线| 久久国产精品男人的天堂亚洲| 亚洲精品美女久久久久99蜜臀| 黄色视频不卡| 69精品国产乱码久久久| 成人18禁在线播放| 在线观看舔阴道视频| 老熟妇乱子伦视频在线观看| 国产精品98久久久久久宅男小说| 波多野结衣巨乳人妻| 亚洲国产中文字幕在线视频| 日韩有码中文字幕| 桃红色精品国产亚洲av| 成年人黄色毛片网站| 视频在线观看一区二区三区| 视频在线观看一区二区三区| 欧美成狂野欧美在线观看| 后天国语完整版免费观看| 一级作爱视频免费观看| 精品久久蜜臀av无| 国产成人av激情在线播放| 亚洲色图av天堂| 亚洲色图av天堂| 999久久久精品免费观看国产| 午夜福利18| 日本精品一区二区三区蜜桃| 可以在线观看的亚洲视频| 欧美最黄视频在线播放免费| 成人国产一区最新在线观看| 嫩草影院精品99| 很黄的视频免费| 一区二区三区国产精品乱码| 色在线成人网| cao死你这个sao货| 欧美日韩亚洲国产一区二区在线观看| 黄色成人免费大全| 99精品欧美一区二区三区四区| 禁无遮挡网站| 国产精品野战在线观看| 亚洲精品一区av在线观看| 亚洲性夜色夜夜综合| 国产97色在线日韩免费| 美女国产高潮福利片在线看| 国产精品久久久久久精品电影 | 别揉我奶头~嗯~啊~动态视频| 我的亚洲天堂| 18禁观看日本| 久久久久久国产a免费观看| 操出白浆在线播放| 精品久久久精品久久久| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放 | 999久久久精品免费观看国产| 老熟妇仑乱视频hdxx| 自线自在国产av| 亚洲国产看品久久| 亚洲精品美女久久久久99蜜臀| 狠狠狠狠99中文字幕| 欧美人与性动交α欧美精品济南到| 最新在线观看一区二区三区| 国内精品久久久久精免费| 欧美国产精品va在线观看不卡| 999精品在线视频| 搡老岳熟女国产| av有码第一页| a在线观看视频网站| 可以免费在线观看a视频的电影网站| 人妻丰满熟妇av一区二区三区| 欧美精品啪啪一区二区三区| 侵犯人妻中文字幕一二三四区| 制服人妻中文乱码| 18禁国产床啪视频网站| 大型av网站在线播放| 最新在线观看一区二区三区| 精品卡一卡二卡四卡免费| 一区二区三区激情视频| 精品久久久久久,| 欧美日本中文国产一区发布| 两性夫妻黄色片| 国产成人啪精品午夜网站| 又黄又爽又免费观看的视频| 最好的美女福利视频网| 亚洲中文av在线| 黄片大片在线免费观看| 久久天躁狠狠躁夜夜2o2o| 欧美日韩乱码在线| 满18在线观看网站| 国产亚洲精品久久久久5区| 后天国语完整版免费观看| 一级毛片女人18水好多| 国产99白浆流出| 欧美黄色片欧美黄色片| 国产成人影院久久av| av天堂久久9| 黄色成人免费大全| 久久亚洲精品不卡| 黄色视频,在线免费观看| 可以在线观看的亚洲视频| 日日夜夜操网爽| 亚洲专区国产一区二区| 久久狼人影院| 久久久久久久久久久久大奶| 国产成人欧美在线观看| 午夜福利影视在线免费观看| 亚洲成人精品中文字幕电影| 国产av精品麻豆| 国产欧美日韩一区二区精品| 熟女少妇亚洲综合色aaa.| 久久久久久久久免费视频了| 天天一区二区日本电影三级 | 长腿黑丝高跟| 日韩高清综合在线| 精品欧美一区二区三区在线| 热re99久久国产66热| 亚洲精品中文字幕在线视频| 欧美午夜高清在线| 香蕉丝袜av| 大香蕉久久成人网| 熟女少妇亚洲综合色aaa.| 91九色精品人成在线观看| 亚洲精品在线观看二区| or卡值多少钱| 欧美日韩一级在线毛片| 精品国产乱子伦一区二区三区| 黄网站色视频无遮挡免费观看| 一a级毛片在线观看| 亚洲黑人精品在线| www国产在线视频色| 精品卡一卡二卡四卡免费| 少妇 在线观看| 九色亚洲精品在线播放| 18禁观看日本| 国产又爽黄色视频| 国产欧美日韩精品亚洲av| 亚洲中文av在线| 两个人视频免费观看高清| 国产主播在线观看一区二区| 18禁观看日本| 久热这里只有精品99| 国产精品久久久久久精品电影 | 亚洲国产毛片av蜜桃av| 好男人电影高清在线观看| 一进一出好大好爽视频| 午夜精品在线福利| 色老头精品视频在线观看| 18禁国产床啪视频网站| 狂野欧美激情性xxxx| 少妇 在线观看| 国产一区二区在线av高清观看| 97超级碰碰碰精品色视频在线观看| 美女国产高潮福利片在线看| 久久久水蜜桃国产精品网| 国产成人欧美在线观看| 12—13女人毛片做爰片一| 亚洲av日韩精品久久久久久密| 一级片免费观看大全| 十分钟在线观看高清视频www| 法律面前人人平等表现在哪些方面| 久久国产亚洲av麻豆专区| 亚洲欧美日韩无卡精品| tocl精华| 成人国产综合亚洲| 成人免费观看视频高清| 精品欧美一区二区三区在线| 亚洲专区中文字幕在线| 99riav亚洲国产免费| 免费一级毛片在线播放高清视频 | svipshipincom国产片| 亚洲一区二区三区色噜噜| 91av网站免费观看| 伦理电影免费视频| 久久国产精品影院| 精品久久久久久,| 美女高潮喷水抽搐中文字幕| 日本a在线网址| 久久午夜综合久久蜜桃| 波多野结衣高清无吗| 国产精品 国内视频| 国产成人免费无遮挡视频| 国产精品香港三级国产av潘金莲| 日本a在线网址| 久久香蕉激情| 日本三级黄在线观看| 欧美乱码精品一区二区三区| 婷婷六月久久综合丁香| 欧美日韩一级在线毛片| 日韩精品中文字幕看吧| 免费高清视频大片| 多毛熟女@视频| 欧美+亚洲+日韩+国产| av网站免费在线观看视频| 丁香六月欧美| 欧美老熟妇乱子伦牲交| 色av中文字幕| 涩涩av久久男人的天堂| 成年人黄色毛片网站| 亚洲自偷自拍图片 自拍| 亚洲成av片中文字幕在线观看| 国产精品亚洲av一区麻豆| 国产精品亚洲一级av第二区| 国产精品av久久久久免费| 在线播放国产精品三级| 国产xxxxx性猛交| 午夜精品久久久久久毛片777| 国产免费男女视频| av有码第一页| 黄色a级毛片大全视频| 亚洲成人精品中文字幕电影| tocl精华| 激情在线观看视频在线高清| 久久久久久久午夜电影| 国产日韩一区二区三区精品不卡| 色在线成人网| 久久久精品国产亚洲av高清涩受| 亚洲人成网站在线播放欧美日韩| 两个人视频免费观看高清| 丁香六月欧美| 制服诱惑二区| 香蕉久久夜色| 90打野战视频偷拍视频| 久久久久久久久免费视频了| 两个人免费观看高清视频| 黄色片一级片一级黄色片| 亚洲欧美日韩另类电影网站| 亚洲国产精品成人综合色| 一级毛片高清免费大全| 宅男免费午夜| 日本在线视频免费播放| 97人妻精品一区二区三区麻豆 | 一进一出抽搐动态| 亚洲午夜精品一区,二区,三区| 国产一区二区在线av高清观看| 欧美国产日韩亚洲一区| 成人三级黄色视频| 老司机午夜十八禁免费视频| 国产高清视频在线播放一区| 香蕉国产在线看| 欧美不卡视频在线免费观看 | 91精品三级在线观看| 亚洲男人的天堂狠狠| 成人国产一区最新在线观看| 久久香蕉激情| 激情视频va一区二区三区| 狠狠狠狠99中文字幕| 免费一级毛片在线播放高清视频 | 国产1区2区3区精品| 久热这里只有精品99| 欧美丝袜亚洲另类 | 丁香六月欧美| 99精品欧美一区二区三区四区| 亚洲一区高清亚洲精品| 黄频高清免费视频| 91老司机精品| 1024香蕉在线观看| 国内久久婷婷六月综合欲色啪| 国产男靠女视频免费网站| 成人手机av| 日韩国内少妇激情av| 成人免费观看视频高清| 久久青草综合色| av在线播放免费不卡| 国产精品秋霞免费鲁丝片| 亚洲精品一卡2卡三卡4卡5卡| 久久中文字幕人妻熟女| 精品卡一卡二卡四卡免费| 男女床上黄色一级片免费看| 91大片在线观看| 国产精品久久久av美女十八| 热99re8久久精品国产| 亚洲电影在线观看av| 国产成人精品无人区| 亚洲午夜理论影院| 一个人免费在线观看的高清视频| 99国产精品99久久久久| 久久香蕉国产精品| 可以在线观看毛片的网站| 国产精品,欧美在线| 熟女少妇亚洲综合色aaa.| 午夜福利18| 无人区码免费观看不卡| 午夜精品国产一区二区电影| 波多野结衣一区麻豆| 国产视频一区二区在线看| 亚洲国产精品成人综合色| 国产精品亚洲美女久久久| 国产一区二区激情短视频| 天堂影院成人在线观看| 他把我摸到了高潮在线观看| 国产精品永久免费网站| 黄色女人牲交| 亚洲欧美精品综合久久99| 乱人伦中国视频| 国产av又大| 亚洲专区国产一区二区| 1024视频免费在线观看| 亚洲欧美精品综合一区二区三区| 一级片免费观看大全| 欧美国产精品va在线观看不卡| 大陆偷拍与自拍| 黄色a级毛片大全视频| 国产一卡二卡三卡精品| 日韩一卡2卡3卡4卡2021年| 在线观看www视频免费| 午夜亚洲福利在线播放| 给我免费播放毛片高清在线观看| 一区二区三区精品91| 啦啦啦韩国在线观看视频| 国产又色又爽无遮挡免费看| 别揉我奶头~嗯~啊~动态视频| 欧美日韩精品网址| 啦啦啦 在线观看视频| 99精品久久久久人妻精品| 91在线观看av| 动漫黄色视频在线观看| 一个人观看的视频www高清免费观看 | 欧美+亚洲+日韩+国产| 久久精品91无色码中文字幕| 脱女人内裤的视频| 美女午夜性视频免费| 欧美成狂野欧美在线观看| 亚洲专区字幕在线| 亚洲精品美女久久久久99蜜臀| 香蕉国产在线看| 热re99久久国产66热| 久久精品国产亚洲av香蕉五月| 18禁黄网站禁片午夜丰满| 成年版毛片免费区| 亚洲中文av在线| 日本在线视频免费播放| 久久久久亚洲av毛片大全| 最近最新免费中文字幕在线| 看黄色毛片网站| 两个人免费观看高清视频| 久久这里只有精品19| 麻豆成人av在线观看| 国产精品日韩av在线免费观看 | 男人舔女人下体高潮全视频| 午夜老司机福利片| 亚洲三区欧美一区| 99国产精品一区二区蜜桃av| 一级毛片高清免费大全| 免费观看精品视频网站| 久久久久亚洲av毛片大全| 亚洲人成电影观看| 欧美乱码精品一区二区三区| 99热只有精品国产| 精品一区二区三区四区五区乱码| 亚洲五月色婷婷综合| а√天堂www在线а√下载| 真人做人爱边吃奶动态| 久久狼人影院| 精品国产一区二区久久| 激情在线观看视频在线高清| 91成年电影在线观看| 制服人妻中文乱码| 搡老岳熟女国产| 无人区码免费观看不卡| 国产精品一区二区精品视频观看| 91精品三级在线观看| 色精品久久人妻99蜜桃| 成人精品一区二区免费| 婷婷六月久久综合丁香| 免费在线观看亚洲国产| 国产午夜福利久久久久久| 满18在线观看网站| 亚洲一区二区三区不卡视频| 90打野战视频偷拍视频| 久久国产亚洲av麻豆专区| 亚洲第一青青草原| 国产精品亚洲av一区麻豆| 午夜久久久在线观看| 男人舔女人的私密视频| 久久久久久久久久久久大奶| 一区二区三区国产精品乱码| 他把我摸到了高潮在线观看| 久久久久久亚洲精品国产蜜桃av| 国产一区二区三区视频了| av视频在线观看入口| 久久 成人 亚洲| 欧美激情 高清一区二区三区| 可以免费在线观看a视频的电影网站| 国产精品综合久久久久久久免费 | 亚洲国产中文字幕在线视频| 91字幕亚洲| 99国产精品一区二区蜜桃av| 一级毛片高清免费大全| 国产乱人伦免费视频| 91字幕亚洲| 18禁裸乳无遮挡免费网站照片 | 给我免费播放毛片高清在线观看| 777久久人妻少妇嫩草av网站| 麻豆av在线久日| 侵犯人妻中文字幕一二三四区| 一区二区三区精品91| 老汉色∧v一级毛片| 日韩有码中文字幕| 亚洲少妇的诱惑av| 色av中文字幕| 精品乱码久久久久久99久播| 91精品国产国语对白视频| 黑丝袜美女国产一区| 超碰成人久久| 丝袜人妻中文字幕| 午夜精品国产一区二区电影| 国产精品亚洲av一区麻豆| 久久狼人影院| 国产在线精品亚洲第一网站| 黄色视频不卡| 国产成+人综合+亚洲专区| 最近最新中文字幕大全免费视频| 麻豆成人av在线观看| 久久午夜亚洲精品久久| 韩国精品一区二区三区| 成熟少妇高潮喷水视频| 18禁国产床啪视频网站| 91精品国产国语对白视频| 欧美日韩亚洲综合一区二区三区_| 免费在线观看视频国产中文字幕亚洲| 老鸭窝网址在线观看| 亚洲 欧美 日韩 在线 免费| 日韩有码中文字幕| 国产成人系列免费观看| 丰满人妻熟妇乱又伦精品不卡| 精品乱码久久久久久99久播| 脱女人内裤的视频| 欧美黑人欧美精品刺激| 精品卡一卡二卡四卡免费| 亚洲精品中文字幕一二三四区| 在线观看66精品国产| 亚洲色图 男人天堂 中文字幕| 亚洲 国产 在线| 在线av久久热| 两个人免费观看高清视频| 我的亚洲天堂| 丝袜美腿诱惑在线| 一级黄色大片毛片| 国产熟女午夜一区二区三区| 村上凉子中文字幕在线| 级片在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 97人妻天天添夜夜摸| 18禁裸乳无遮挡免费网站照片 | 精品国产一区二区三区四区第35| 精品第一国产精品| 国产一区二区三区在线臀色熟女| 久久久久国产精品人妻aⅴ院| 久久精品成人免费网站| 成人国产一区最新在线观看| 校园春色视频在线观看| 亚洲熟妇中文字幕五十中出| 一级毛片高清免费大全| 国产成人精品无人区| 久久精品91蜜桃| 午夜老司机福利片| 亚洲精品一区av在线观看| 精品国产一区二区久久| 91成人精品电影| 国产精品久久电影中文字幕| 免费av毛片视频| 多毛熟女@视频| 国产免费av片在线观看野外av| 免费观看精品视频网站| 亚洲av第一区精品v没综合| 中文亚洲av片在线观看爽| 国产极品粉嫩免费观看在线| 国产一区在线观看成人免费| 亚洲在线自拍视频| 日韩欧美免费精品| 一区二区三区激情视频| 国产97色在线日韩免费| 宅男免费午夜| 91九色精品人成在线观看| 亚洲五月天丁香| 欧美一区二区精品小视频在线| 中文字幕色久视频| 亚洲,欧美精品.| 亚洲午夜理论影院| 久久精品成人免费网站| 日韩精品青青久久久久久| 首页视频小说图片口味搜索| 日韩成人在线观看一区二区三区| 嫁个100分男人电影在线观看| 在线av久久热| 亚洲国产精品sss在线观看| 非洲黑人性xxxx精品又粗又长| 天天躁夜夜躁狠狠躁躁| 国产av又大| 久久久久久人人人人人| 免费在线观看黄色视频的| 神马国产精品三级电影在线观看 | 日本黄色视频三级网站网址| 欧美在线一区亚洲| 麻豆成人av在线观看| 久久九九热精品免费| 九色亚洲精品在线播放| 69av精品久久久久久| 搡老熟女国产l中国老女人| 国产精品日韩av在线免费观看 | 国产一区二区三区在线臀色熟女| 国产亚洲精品综合一区在线观看 | 女警被强在线播放| 午夜精品国产一区二区电影| 老司机午夜十八禁免费视频| 欧美日本亚洲视频在线播放| 成熟少妇高潮喷水视频| 日韩成人在线观看一区二区三区| 国产亚洲欧美精品永久| 国产主播在线观看一区二区| 国产高清videossex| 亚洲男人天堂网一区| 看片在线看免费视频| 欧美黑人欧美精品刺激| 国产精品野战在线观看| 嫩草影院精品99| 制服人妻中文乱码| 99热只有精品国产| 精品国产乱子伦一区二区三区| 老汉色av国产亚洲站长工具| 丁香六月欧美| 精品久久久久久久毛片微露脸| 啦啦啦观看免费观看视频高清 | 啦啦啦韩国在线观看视频| 啦啦啦 在线观看视频| 丝袜美足系列| 久久人妻福利社区极品人妻图片| 咕卡用的链子| 亚洲七黄色美女视频| 国产精品免费视频内射| 给我免费播放毛片高清在线观看| 每晚都被弄得嗷嗷叫到高潮| 欧美成狂野欧美在线观看| 丝袜在线中文字幕| 中亚洲国语对白在线视频| 一级a爱视频在线免费观看| 亚洲精品久久国产高清桃花| 午夜两性在线视频| 香蕉国产在线看| 好男人电影高清在线观看| 两人在一起打扑克的视频| 天天躁狠狠躁夜夜躁狠狠躁| 精品午夜福利视频在线观看一区| 男人舔女人下体高潮全视频| 午夜久久久久精精品| 变态另类丝袜制服| 欧美在线一区亚洲| 在线观看免费视频日本深夜| 男人操女人黄网站| www.自偷自拍.com| 狂野欧美激情性xxxx| 日韩欧美一区视频在线观看| 国产精品一区二区免费欧美| 国产一区二区在线av高清观看| 精品国产国语对白av| 欧美在线一区亚洲| 一区二区三区激情视频| 欧美日韩一级在线毛片| 亚洲伊人色综图| 两人在一起打扑克的视频| 亚洲精品久久成人aⅴ小说| 制服丝袜大香蕉在线| 国产人伦9x9x在线观看| 99热只有精品国产| 国产99久久九九免费精品| 国产精品爽爽va在线观看网站 | 老司机午夜十八禁免费视频| 免费少妇av软件| 成人18禁在线播放| 十八禁网站免费在线| 成人欧美大片| 欧美激情高清一区二区三区| 一二三四在线观看免费中文在| 午夜精品在线福利| 亚洲伊人色综图| 亚洲精品国产一区二区精华液| 在线免费观看的www视频| 99国产精品99久久久久| 最近最新中文字幕大全免费视频| 亚洲欧美激情在线| 精品熟女少妇八av免费久了| 亚洲国产精品sss在线观看| 国产高清videossex| 久久久久久国产a免费观看| 超碰成人久久| 一区二区三区精品91| 黑人巨大精品欧美一区二区蜜桃| 人妻久久中文字幕网| 人人妻人人爽人人添夜夜欢视频| 亚洲av成人一区二区三| 国产精品久久久久久人妻精品电影| 亚洲成人免费电影在线观看| 免费在线观看亚洲国产| 一二三四社区在线视频社区8| 一级a爱视频在线免费观看|