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

    黃渤海藍(lán)點(diǎn)馬鮫種群變動(dòng)與海表面溫度的關(guān)系?

    2019-11-07 01:46:18王凱迪于華明于海慶伊藤進(jìn)一

    王凱迪, 于華明??, 于海慶, 伊藤進(jìn)一

    (1.中國(guó)海洋大學(xué)海洋與大氣學(xué)院,山東 青島 266100;2.中國(guó)海洋大學(xué)物理海洋教育部重點(diǎn)實(shí)驗(yàn)室,山東 青島 266003; 3. 青島海洋科學(xué)與技術(shù)國(guó)家實(shí)驗(yàn)室,山東 青島 266237;4.東京大學(xué)大氣海洋研究所,日本 千葉縣 277-8564)

    藍(lán)點(diǎn)馬鮫(Scomberomorusniphonius, Japanese Spanish mackerel)隸屬于鱸形目(Perciformes)鯖科(Scombridae)馬鮫屬(Scomberomorus),俗稱藍(lán)點(diǎn)鮫、鲅魚、條燕、板鲅、竹鮫、尖頭馬加、馬鮫、青箭,屬海洋中上層暖水肉食性魚類,具有攝食量大、游速快、游程遠(yuǎn)等特征,有洄游習(xí)性,主要分布于渤海、黃海、東海、日本海以及日本太平洋沿岸海域,其漁業(yè)資源為中國(guó)、日本、韓國(guó)等國(guó)共同開發(fā)利用。藍(lán)點(diǎn)馬鮫是我國(guó)產(chǎn)量較大的重要經(jīng)濟(jì)魚種之一,具有較高的經(jīng)濟(jì)價(jià)值,是目前黃渤海唯一年漁獲量超過100 000 t的大型經(jīng)濟(jì)魚類[1]。

    前人針對(duì)藍(lán)點(diǎn)馬鮫的生物學(xué)與漁業(yè)資源變動(dòng)特征進(jìn)行了一定研究,如Shoji J等發(fā)現(xiàn)藍(lán)點(diǎn)馬鮫以其他魚類仔魚為食[2-3],且攝食有明顯的晝夜變化[4]等;韋晟等對(duì)黃渤海藍(lán)點(diǎn)馬鮫種群鑒別進(jìn)行了研究, 認(rèn)為黃海南部、中部、北部及渤海諸海域間藍(lán)點(diǎn)馬鮫均屬同一種群,而黃海南部藍(lán)點(diǎn)馬鮫可視作同一種群的不同群體[5];水柏年等提出個(gè)體較強(qiáng)的擴(kuò)散能力和長(zhǎng)距離的洄游可能是藍(lán)點(diǎn)馬鮫群體無(wú)遺傳分化的重要原因[6];袁楊洋等分析了黃海南部水域藍(lán)點(diǎn)馬鮫洄游行為及與漁場(chǎng)形成、漁期更迭關(guān)系,指出水溫是藍(lán)點(diǎn)馬鮫漁場(chǎng)形成、漁期更迭的決定因素[7];宋超、王宇壇、林楠、楊林林等對(duì)象山港藍(lán)點(diǎn)馬鮫魚卵、仔稚魚的分布、生長(zhǎng)、攝食、耳石形態(tài)等進(jìn)行了詳細(xì)研究[8-11];邱盛堯等討論了我國(guó)漁業(yè)管理對(duì)藍(lán)點(diǎn)馬鮫漁業(yè)資源的貢獻(xiàn)[12];孫本曉對(duì)2006—2008年黃渤海藍(lán)點(diǎn)馬鮫資源特性與變動(dòng)進(jìn)行了研究,并提出種質(zhì)資源保護(hù)的措施和建議[13]。

    20世紀(jì)末,隨著全球變暖,海洋溫度整體上升,強(qiáng)ENSO事件頻繁發(fā)生,物理環(huán)境的變化逐漸成為漁業(yè)種類組成和種群結(jié)構(gòu)改變的重要原因之一。在此背景下,作為大型、高營(yíng)養(yǎng)級(jí)的重要經(jīng)濟(jì)魚類,藍(lán)點(diǎn)馬鮫漁獲量卻仍在波動(dòng)中上升,呈現(xiàn)出與其他主要魚種截然不同的變動(dòng)特征[12-15],這說明藍(lán)點(diǎn)馬鮫的生活史或許具有應(yīng)對(duì)物理環(huán)境變化多重壓力的強(qiáng)韌性和柔軟性[16]。藍(lán)點(diǎn)馬鮫作為高營(yíng)養(yǎng)層次的關(guān)鍵種,位于海洋食物鏈的頂端,其資源變動(dòng)不僅受到小型中上層魚類等餌料生物的影響,也會(huì)影響到餌料魚種和海洋生態(tài)系統(tǒng)變動(dòng),在海洋生態(tài)系統(tǒng)中發(fā)揮著不可替代的關(guān)鍵作用。因此,研究藍(lán)點(diǎn)馬鮫對(duì)海洋環(huán)境變化的響應(yīng)有助于我們進(jìn)一步揭示藍(lán)點(diǎn)馬鮫種群變動(dòng)機(jī)理和生態(tài)系統(tǒng)的調(diào)控機(jī)制,具有一定的科學(xué)意義。然而,由于藍(lán)點(diǎn)馬鮫是高度洄游性的大型魚類,生活史過程與種群結(jié)構(gòu)相對(duì)復(fù)雜,數(shù)據(jù)在時(shí)間、空間尺度較為離散,連續(xù)多年、大尺度研究比較困難;藍(lán)點(diǎn)馬鮫種群數(shù)量受人類活動(dòng)影響較大,難以建立與海洋水文環(huán)境之間的對(duì)應(yīng)關(guān)系;且海洋漁業(yè)、物理海洋學(xué)科均具有較強(qiáng)專業(yè)性,學(xué)科之間的界限很難打破。因此,針對(duì)藍(lán)點(diǎn)馬鮫種群動(dòng)態(tài)對(duì)環(huán)境變化響應(yīng)仍缺乏深入的研究。

    在傳統(tǒng)漁業(yè)資源評(píng)估中,單位捕獲努力量漁獲量(CPUE,Catch Per Unit Effort)常用作為資源豐度的標(biāo)志量,但研究表明,受社會(huì)需求、捕撈成本、捕撈技術(shù)、捕撈區(qū)域變化等因素影響,CPUE與資源量的正比關(guān)系并不完全成立,且CPUE標(biāo)準(zhǔn)化可能隨魚類種群、時(shí)空尺度、模型方法選擇的不同存在偏差[17-18]。而實(shí)際種群分析方法(VPA, Virtual Population Analysis)將魚類資源總量視作漁獲量(Catches)、自然死亡量(Natural losses)與幸存量(Survivors)三者的總和[19],把自然與社會(huì)對(duì)魚類資源變動(dòng)的影響分開,將人類捕撈活動(dòng)等社會(huì)影響納入漁獲量中,消除社會(huì)因素對(duì)其他變量的干擾。該方法如今已被廣泛應(yīng)用于漁業(yè)資源評(píng)估與管理工作中[20-22]。再生產(chǎn)成功率(Recruitment Per Spawner, RPS)是子魚與親魚幸存量的比值[23-24],已剔除人為影響,只與自然環(huán)境有關(guān),相比CPUE,RPS能夠更加客觀、準(zhǔn)確地描繪魚類種群變動(dòng)對(duì)自然環(huán)境的響應(yīng)。

    本文根據(jù)1992—2012年黃渤海海域藍(lán)點(diǎn)馬鮫捕獲量數(shù)據(jù)和美國(guó)國(guó)家海洋大氣局(National Oceanic and Atmospheric Administration,NOAA)的最優(yōu)插值海表面溫度(Sea Surface Temperature ,SST)數(shù)據(jù),利用實(shí)際種群分析方法推算了黃渤海各年份的RPS,分析了藍(lán)點(diǎn)馬鮫RPS與SST的關(guān)系??蔀楹Q鬂O業(yè)與海洋水動(dòng)力研究的學(xué)科交叉與學(xué)科融合提供新的思路,為進(jìn)一步深入研究藍(lán)點(diǎn)馬鮫對(duì)物理環(huán)境的響應(yīng)提供必要參考,對(duì)我國(guó)漁業(yè)管理和藍(lán)點(diǎn)馬鮫資源的可持續(xù)利用具有一定的數(shù)據(jù)支撐和指導(dǎo)作用。

    1 我國(guó)近海藍(lán)點(diǎn)馬鮫洄游及分布規(guī)律

    藍(lán)點(diǎn)馬鮫從出生到死亡,產(chǎn)卵、索餌、越冬等生活史特性具有明顯的時(shí)空尺度和地域特性。

    我國(guó)近海藍(lán)點(diǎn)馬鮫產(chǎn)卵場(chǎng)主要有渤海的萊州灣、渤海灣、遼東灣,黃海的煙威近海、海洋島附近水域、青島近海、海州灣和江蘇呂泗沿海,東海的長(zhǎng)江口海域、舟山近海、溫州近海、官井洋附近水域等[5-6,25-26]。越冬場(chǎng)主要分布在黃海東南部濟(jì)州島西南部海域、東海中南部舟山與溫州外海海域及東海南部至南海北部海域[5-6,13]。

    藍(lán)點(diǎn)馬鮫的越冬期為1—2 月。每年初春(3月中旬),海域水溫升高,促進(jìn)了藍(lán)點(diǎn)馬鮫性腺的加速發(fā)育,促使其游離越冬場(chǎng)進(jìn)行生殖洄游[27]。藍(lán)點(diǎn)馬鮫的產(chǎn)卵期大致為4—6月,但由于我國(guó)沿??缭蕉鄠€(gè)緯度,南北海區(qū)水溫梯度大,故不同海區(qū)產(chǎn)卵期略有差別,大致由南至北相應(yīng)推遲[25](見表1)。藍(lán)點(diǎn)馬鮫產(chǎn)卵之后,親魚和幼魚均在產(chǎn)卵場(chǎng)附近海域分散索餌。至8月下旬,隨著近岸水溫下降,魚群陸續(xù)向較深水域進(jìn)行適溫洄游,繼續(xù)攝食,生長(zhǎng)育肥,并于12月中下旬陸續(xù)回到各個(gè)越冬場(chǎng)。藍(lán)點(diǎn)馬鮫的產(chǎn)卵場(chǎng)、越冬場(chǎng)及洄游路線見圖1,粉色為主要產(chǎn)卵場(chǎng),藍(lán)色為主要越冬場(chǎng),實(shí)心箭頭為產(chǎn)卵洄游路線,空心箭頭為越冬洄游路線,箭頭上的月份為部分洄游發(fā)生時(shí)間。

    圖1 藍(lán)點(diǎn)馬鮫主要產(chǎn)卵場(chǎng)、越冬場(chǎng)分布及洄游路線Fig.1 Main spawning grounds, wintering grounds and migration trajectory for Scomberomorus niphonius

    表1我國(guó)沿海藍(lán)點(diǎn)馬鮫主要產(chǎn)卵期
    Table 1 Main spawning season in coastal areas of China

    海域Sea area產(chǎn)卵期Spawning period產(chǎn)卵盛期Peak of spawning period最適溫度Optimum temperature/℃渤海The Bohai Sea4月下旬至6月下旬5月下旬至6月中旬14~17黃海北部North of the Yellow Sea4月下旬至6月上旬5月中旬至5月下旬10~15黃海中部Middle of the Yellow Sea4月中下旬至5月下旬5月上旬至5月下旬11~14黃海南部South of the Yellow Sea4月中旬至5月中旬4月下旬至5月上旬9~13福建沿海Coastal area of Fujian Province3月至6月4月上旬至5月下旬

    2 數(shù)據(jù)來(lái)源與處理

    藍(lán)點(diǎn)馬鮫作為我國(guó)近海最重要的漁業(yè)資源種類之一,雖然已有學(xué)者針對(duì)其生物學(xué)特性及漁業(yè)資源變化做了許多工作,但是在對(duì)環(huán)境的響應(yīng)方面缺乏深入研究。近年來(lái)越來(lái)越多的研究證明,無(wú)論是從單一的海洋生物到整個(gè)海洋生態(tài)系統(tǒng),都或多或少受環(huán)境變化的脅迫[31-32]。而藍(lán)點(diǎn)馬鮫從出生到死亡,受海洋環(huán)境變化影響很大[7]。海水溫度與鹽度是描述海洋環(huán)境狀況中最重要的因子,常作為研究水文性質(zhì),描述水體運(yùn)動(dòng)的基本指標(biāo),由于藍(lán)點(diǎn)馬鮫對(duì)鹽度適應(yīng)性較強(qiáng)[2,28-29],對(duì)海表面溫度要求嚴(yán)格[6,26-27,29-30],且早期生活史過程主要發(fā)生在海水表層,故本文討論藍(lán)點(diǎn)馬鮫種群變動(dòng)與SST的關(guān)系。藍(lán)點(diǎn)馬鮫種群變動(dòng)的標(biāo)志量選用再生產(chǎn)成功率(RPS),由于在計(jì)算過程中剔除了人類捕撈活動(dòng)的影響,RPS只與自然環(huán)境有關(guān)[23-24]。

    2.1 數(shù)據(jù)來(lái)源

    1992—2012年黃渤海海域藍(lán)點(diǎn)馬鮫捕獲量數(shù)據(jù)及年齡組成資料引自農(nóng)業(yè)部漁業(yè)局編制的《中國(guó)漁業(yè)統(tǒng)計(jì)年鑒》和邱盛堯、孫本曉等[1,12-13,34-35]的歷史資料。1992—2012年的海表面溫度采用美國(guó)國(guó)家海洋大氣局的最優(yōu)插值海表面溫度(Optimum Interpolation Sea Surface Temperature,OISST)的日平均資料,空間分辨率為0.25°×0.25°,具有較高的準(zhǔn)確性[33]。

    2.2 處理方法

    本章采用實(shí)際種群分析方法(VPA),基于年齡結(jié)構(gòu)數(shù)據(jù)的漁獲量(Catch)數(shù)據(jù),概算各年齡組的幸存量(Survivors) 由于歷史資料中僅有年總捕獲量的數(shù)據(jù),沒有各年齡層捕獲量,因此需要對(duì)年齡組成做出假設(shè),將藍(lán)點(diǎn)馬鮫群體按年齡劃分為5組,分別為1齡、2齡、3齡、4齡和4+齡組,其中4+齡組為5齡及以上的藍(lán)點(diǎn)馬鮫群體。我國(guó)是藍(lán)點(diǎn)馬鮫漁業(yè)資源的主要開發(fā)國(guó),藍(lán)點(diǎn)馬鮫資源的開發(fā)利用可分為 5 個(gè)階段: 1962 年以前的原始階段,1962—1976 年的發(fā)展階段,1977—1989 年的充分利用階段,1990—1997 年的衰退階段,1997 年以后的管理階段[12-13](見表2)。本文主要討論1992—2012年的藍(lán)點(diǎn)馬鮫資源數(shù)量變動(dòng),因此將該時(shí)間段分為衰退時(shí)期(1992—1996年)和管理階段(1997年之后)。衰退階段的藍(lán)點(diǎn)馬鮫直接引用1992—1996年黃渤海藍(lán)點(diǎn)馬鮫資源數(shù)量變動(dòng)監(jiān)測(cè)數(shù)據(jù)的年齡組成比例[34-35];管理階段由于缺少逐年年齡組成資料,因此取2006—2008年黃渤海藍(lán)點(diǎn)馬鮫資源數(shù)量變動(dòng)監(jiān)測(cè)數(shù)據(jù)的平均年齡組成,1—4+齡的年齡組成比例分別為32.0%、20.6%、12.9%、10.5%和24.0%[13]。

    根據(jù)以上假設(shè)可得到具有年齡結(jié)構(gòu)的藍(lán)點(diǎn)馬鮫漁獲量數(shù)據(jù),然后利用R軟件的TropFishR(Tropical Fisheries Analysis)工具包的VPA函數(shù),對(duì)年齡結(jié)構(gòu)化的漁獲量數(shù)據(jù)進(jìn)行處理,求得每年各年齡層的幸存量。其中,根據(jù)Ven-Bertalanffy[36]生長(zhǎng)方程對(duì)藍(lán)點(diǎn)馬鮫的年齡與叉長(zhǎng)、體重的關(guān)系描述:

    Lt=746(1-e-0.460 6(t+0.733)),

    Wt=3.157(1-e-0.460 6(t+0.733))3。

    計(jì)算可得1、2、3、4齡魚的體重分別0.52、1.16、1.76、2.20 kg,4+齡魚取5~7齡魚平均體重2.72 kg,并建立叉長(zhǎng)-體重關(guān)系式W=c×Ld,c、d為常量,根據(jù)生長(zhǎng)方程分別取0.000 023和2.815 9。

    再生產(chǎn)成功率(RPS)為次年1齡魚與當(dāng)年性成熟成魚幸存量的比值,根據(jù)歷史資料,1990年代以后1齡魚性成熟比例已達(dá)90%以上[34-35],設(shè)1齡魚的性成熟比例為0.9,2、3、4、4+齡魚默認(rèn)已全部性成熟。由于藍(lán)點(diǎn)馬鮫的繁殖力隨年齡增長(zhǎng)而增加,不同年齡層藍(lán)點(diǎn)馬鮫對(duì)補(bǔ)充群體的繁殖貢獻(xiàn)率不同,因此引入繁殖力修正系數(shù)對(duì)RPS結(jié)果進(jìn)行修正。根據(jù)邱盛堯等[37]的研究,個(gè)體絕對(duì)生殖力(R)與年齡(t)之間的關(guān)系為:

    Rt=89.01t1.817 3。

    表2 藍(lán)點(diǎn)馬鮫資源開發(fā)的不同階段Table 2 Development stages of Scomberomorus niphonius resource

    各個(gè)年齡層的繁殖力修正系數(shù)分別為nt為:

    nt=Rt/Ravg。

    其中,Ravg為該年平均年齡個(gè)體絕對(duì)生殖力,衰退階段平均年齡為2.65齡,管理階段的藍(lán)點(diǎn)馬鮫平均年齡為2.78齡[13]。以1992年為例,RPS計(jì)算公式如下,其中RPS(1992)為1992年的RPS值,survivors(1993,1)為1993年1齡魚的幸存量,以此類推:

    RPS(1992)=[survivors(1993,1)]/[ 0.9×n1×survivors(1992,1)+n2×survivors(1992,2)+n3×survivors(1992,3)+n4×survivors(1992,4)+n4+×survivors(1992,4+)]。

    由于本文更加關(guān)心年、月時(shí)間尺度的SST變化,因此對(duì)日平均SST資料做年、月平均處理。在比較不同海域的海表面溫度變化趨勢(shì)時(shí),采用極差標(biāo)準(zhǔn)化方法處理:

    Xt′=(Xt-Xmin)/(Xmax-Xmin)。

    3 結(jié)果

    根據(jù)歷史資料,1992—2012年黃渤海海域藍(lán)點(diǎn)馬鮫捕獲量見圖2,最低值出現(xiàn)在1992年,最高值出現(xiàn)在2005年,整體呈波動(dòng)中上升趨勢(shì)。

    圖2 1992—2012年黃渤海藍(lán)點(diǎn)馬鮫漁獲量Fig.2 Catch of Scomberomorus niphonius in the Bohai and Yellow Sea from 1992 to 2012

    基于漁獲量數(shù)據(jù),我們采用實(shí)際種群分析方法(VPA)概算各年齡組的幸存量(Survivors)(見圖3)。由于人類捕撈活動(dòng)與魚類的自然死亡,同年的幸存量隨年齡逐漸遞減,但相同年齡層的幸存量均大致隨年份增加。

    根據(jù)圖3所示各年齡層幸存量,可以求得每年的RPS值。圖4為黃渤海海域藍(lán)點(diǎn)馬鮫RPS隨年份的變化,實(shí)線為RPS值,虛線為線性擬合后的變化趨勢(shì)。由圖3可知,1992—1994年1齡魚幸存量迅速增長(zhǎng),而其他年齡層增幅較小, RPS于1992—1993年有小幅上升;1994—1997年1齡魚幸存量增長(zhǎng)速度放緩,而其他年齡層幸存量增長(zhǎng)速度提高,導(dǎo)致RPS于1993年后開始下降;1997年后各年齡層幸存量增速較小,均呈平穩(wěn)上升狀態(tài),RPS在0.45~0.40之間波動(dòng)下降??傮w而言,RPS在1992—2012年間呈下降趨勢(shì)。

    圖3 1992~2012年各年齡層藍(lán)點(diǎn)馬鮫幸存量

    圖4 黃渤海再生產(chǎn)成功率隨年份的變化Fig.4 Variation of average RPS from 1992 to 2012

    圖5為研究區(qū)域內(nèi)年平均SST隨年份的變化,實(shí)線為SST,虛線為線性擬合后的變化趨勢(shì)。由圖5可知,黃渤海平均溫度總體變化不大,變化范圍為15~16 ℃,最低溫度出現(xiàn)在1993年,最高溫度出現(xiàn)在距離現(xiàn)在最近的2012年及強(qiáng)ENSO事件發(fā)生的1998年,且總體在波動(dòng)中有上升趨勢(shì)。這可能與“溫室效應(yīng)”所導(dǎo)致的全球溫度上升有關(guān)。

    圖6展示了1992—2012年間黃渤海氣候態(tài)月平均海表面溫度(SST),左下角為月份。受太陽(yáng)輻射影響,在北半球深水大洋中SST基本呈“南高北低”狀態(tài);在近岸海域受地形及人類活動(dòng)影響,SST總體上略低于遠(yuǎn)海。在夏季(6、7、8月)溫度最高,最高溫度可達(dá)29 ℃,出現(xiàn)在8月的黃海呂泗沿海及東南部海域,在冬季(12、1、2月)溫度最低,最低溫出現(xiàn)在1—2月渤海海域,約為5 ℃左右。

    圖5 黃渤海平均海表面溫度隨年份的變化Fig.5 Variation of average SST from 1992 to 2012

    圖6 黃渤海1—12月海表面溫度Fig.6 SST from January to December in Bohai and Yellow Sea

    4 討論與分析

    藍(lán)點(diǎn)馬鮫為大范圍洄游性魚種,據(jù)韋晟等[5]研究,黃渤海海域藍(lán)點(diǎn)馬鮫存在兩個(gè)群體,即黃渤海群體與黃海南部群體,不同群體因其生活史、洄游路徑等不同,可能會(huì)對(duì)海洋環(huán)境變動(dòng)產(chǎn)生不同的響應(yīng)。本文以36°N為界劃分研究區(qū)域,區(qū)域a為渤海、黃海北部、黃海中部海域,存在藍(lán)點(diǎn)馬鮫黃渤海群體,區(qū)域b為黃海南部海域,存在藍(lán)點(diǎn)馬鮫黃海南部群體。圖7為1992—2012年平均海表面溫度,黑色虛線為36°N所在位置,全場(chǎng)溫度變化范圍為10~22 ℃,區(qū)域a溫度略低于區(qū)域b。

    圖7 1992—2012年平均海表面溫度Fig.7 Average SST of from 1992 to 2012

    為比較不同區(qū)域溫度變化趨勢(shì),本文對(duì)兩區(qū)域海表面溫度做標(biāo)準(zhǔn)化處理。圖8為標(biāo)準(zhǔn)化后的海表面溫度隨年份的變化,實(shí)線為區(qū)域a,虛線為區(qū)域b。由圖8可知兩區(qū)域最低溫度均出現(xiàn)在1993年,區(qū)域a最高溫度出現(xiàn)在2012年,區(qū)域b最高溫度出現(xiàn)在1998年,兩區(qū)域的海表面溫度變化趨勢(shì)基本相同,相關(guān)系數(shù)達(dá)0.91(p<0.001)。因此我們認(rèn)為海表面溫度變化對(duì)藍(lán)點(diǎn)馬鮫2個(gè)群體的影響是一致的,以下討論將不再劃分區(qū)域。

    圖8 標(biāo)準(zhǔn)化后的平均海表面溫度隨年份的變化Fig.8 Variation of standardized average SST of from 1992 to 2012

    將RPS隨年份變化的時(shí)間序列與計(jì)算海域中各點(diǎn)SST隨年份變化的時(shí)間序列作相關(guān)性分析,可得圖9的結(jié)果,左下角為月份,色標(biāo)中紅色系代表正相關(guān)較強(qiáng),藍(lán)色系代表負(fù)相關(guān)較強(qiáng),黑線包圍的區(qū)域?yàn)橥ㄟ^顯著性檢驗(yàn)的區(qū)域,顯著性水平p<0.05。由圖9可以看出,渤海渤海灣、萊州灣于6月出現(xiàn)明顯負(fù)相關(guān)性區(qū)域,黃海長(zhǎng)江口附近海域4—5月負(fù)相關(guān)性較強(qiáng),海州灣5—6月負(fù)相關(guān)性較強(qiáng),黃海西南海域8—11月負(fù)相關(guān)性較強(qiáng),以上相關(guān)系數(shù)均在0.6~0.8之間。因缺少漁業(yè)數(shù)據(jù),在此對(duì)黃海東部海域不做深入討論。

    圖9 再生產(chǎn)成功率與海表面溫度相關(guān)性分布圖(p<0.05)Fig.9 Correlation of SST and RPS (p<0.05)

    由于藍(lán)點(diǎn)馬鮫具有復(fù)雜的生活史過程與種群結(jié)構(gòu),導(dǎo)致RPS改變的原因較多。從定義來(lái)看,RPS值的大小由補(bǔ)充量與親魚資源量的比值決定,而決定補(bǔ)充量的兩個(gè)過程,即成魚產(chǎn)卵與幼魚生長(zhǎng),均與SST密切相關(guān)。將出現(xiàn)明顯負(fù)相關(guān)性的月份劃為三個(gè)時(shí)期,第一時(shí)期為4—6月,第二時(shí)期為8—11月,分別對(duì)應(yīng)黃渤海藍(lán)點(diǎn)馬鮫的產(chǎn)卵時(shí)期與攝食洄游時(shí)期。由于第一時(shí)期負(fù)相關(guān)系數(shù)較大區(qū)域處存在藍(lán)點(diǎn)馬鮫產(chǎn)卵場(chǎng),第二時(shí)期系數(shù)較大區(qū)域有藍(lán)點(diǎn)馬鮫在此攝食溯游,因此推測(cè)SST可通過影響性成熟成魚產(chǎn)卵與當(dāng)年幼魚的早期攝食過程,導(dǎo)致次年補(bǔ)充量的改變,從而引發(fā)RPS的變化。以下將分時(shí)期進(jìn)行詳細(xì)討論。

    第一時(shí)期(4—6月)為我國(guó)黃渤海藍(lán)點(diǎn)馬鮫的主要產(chǎn)卵期,長(zhǎng)江口附近海域、海州灣、渤海灣與萊州灣均存在藍(lán)點(diǎn)馬鮫產(chǎn)卵場(chǎng),產(chǎn)卵最適溫度為9~17 ℃(見表1)。由于藍(lán)點(diǎn)馬鮫產(chǎn)卵過程對(duì)SST較敏感,過高的溫度會(huì)限制藍(lán)點(diǎn)馬鮫的產(chǎn)卵量與產(chǎn)卵范圍,而這四處海域的實(shí)際溫度均已超過本地最適產(chǎn)卵溫度(見表3),因此隨著溫度的升高,藍(lán)點(diǎn)馬鮫的產(chǎn)卵過程受到抑制,從而影響了第二年藍(lán)點(diǎn)馬鮫群體的補(bǔ)充,使得SST與RPS呈現(xiàn)明顯的負(fù)相關(guān)。

    表3 負(fù)相關(guān)性較強(qiáng)海域的藍(lán)點(diǎn)馬鮫最適溫度與實(shí)際溫度Table 3 Optimum and actual temperature in the area with strong negative correlation

    Note:①Period;②Month with strong negative correlation;③Area with strong negative correlation;④Optimum temperature;⑤Actual temperature;⑥The first period(spawning);⑦The second period(feeding)

    第二時(shí)期(8—11月)是藍(lán)點(diǎn)馬鮫攝食階段,SST除了能直接影響藍(lán)點(diǎn)馬鮫成魚產(chǎn)卵過程外,還可通過影響食物分布間接作用于藍(lán)點(diǎn)馬鮫。8、9月已進(jìn)入藍(lán)點(diǎn)馬鮫當(dāng)年幼魚的主要生長(zhǎng)期[25],一年中藍(lán)點(diǎn)馬鮫在該時(shí)的攝食強(qiáng)度達(dá)到最高[13]。在此階段,食物是決定仔稚魚生長(zhǎng)、存活以及資源量補(bǔ)充的決定因素[38]。藍(lán)點(diǎn)馬鮫在呂泗外海產(chǎn)卵后,在附近及水深較深處索餌,以覓食豐盛的成鳀和幼鳀[27],而鳀魚分布與水溫關(guān)系密切,8月中旬至10月上旬為鳀魚的產(chǎn)卵末期,黃海西南部最適水溫為12~19℃[39],但此時(shí)負(fù)相關(guān)性較強(qiáng)的海域海溫超過20℃,溫度過高會(huì)使得鳀魚幼魚減少,致使藍(lán)點(diǎn)馬鮫幼魚攝食受抑制,影響藍(lán)點(diǎn)馬鮫幼魚成活率,間接致使RPS下降。隨著溫度的進(jìn)一步降低, 11月鳀魚均開始越冬洄游,此時(shí)鳀魚會(huì)向著11~13℃聚集[39],而黃海西南部的溫度仍為17℃以上,此處的藍(lán)點(diǎn)馬鮫因食物缺乏而影響成活率,導(dǎo)致資源補(bǔ)充量的減少與RPS下降。

    5 結(jié)論與研究展望

    本文基于1992—2012年黃渤海海域藍(lán)點(diǎn)馬鮫漁獲量數(shù)據(jù),利用實(shí)際種群分析方法推算了藍(lán)點(diǎn)馬鮫再生產(chǎn)成功率(RPS),并進(jìn)一步探究了其與同時(shí)期海表面溫度(SST)的關(guān)系。結(jié)果表明,藍(lán)點(diǎn)馬鮫的RPS與SST于4月在長(zhǎng)江口附近海域負(fù)相關(guān)性較強(qiáng),于5月在海州灣與長(zhǎng)江口附近負(fù)相關(guān)性較強(qiáng),于6月在渤海灣、萊州灣、海州灣負(fù)相關(guān)性較強(qiáng),于8—11月在黃海西南部海域負(fù)相關(guān)性較強(qiáng),以上相關(guān)性均達(dá)0.6以上(P<0.05)。根據(jù) SST與RPS呈現(xiàn)明顯的負(fù)相關(guān)性的海域位置與時(shí)間,推測(cè)SST對(duì)RPS的影響路徑有兩條:(1)產(chǎn)卵時(shí)期,過高的SST通過限制藍(lán)點(diǎn)馬鮫產(chǎn)卵范圍與產(chǎn)卵量,直接降低次年補(bǔ)充量,從而導(dǎo)致RPS的下降。(2)攝食時(shí)期,過高的SST通過影響藍(lán)點(diǎn)馬鮫的主要食物(鳀魚)的分布,抑制藍(lán)點(diǎn)馬鮫幼魚的攝食過程,降低幼魚成活率,間接降低次年補(bǔ)充量,從而導(dǎo)致RPS的下降。但由于藍(lán)點(diǎn)馬鮫生活史過程與種群結(jié)構(gòu)較為復(fù)雜,魚類資源變動(dòng)是多尺度氣候變化影響的結(jié)果,藍(lán)點(diǎn)馬鮫對(duì)海洋環(huán)境響應(yīng)的更多機(jī)制還需開展進(jìn)一步有針對(duì)性的研究。

    針對(duì)目前藍(lán)點(diǎn)馬鮫研究面臨的問題,未來(lái)可開展藍(lán)點(diǎn)馬鮫洄游模式、年齡與生長(zhǎng)以及資源評(píng)估和管理的研究,開展藍(lán)點(diǎn)馬鮫生活史參數(shù)化方案的研究,建立應(yīng)對(duì)氣候變化的早期預(yù)警機(jī)制等研究。

    国产不卡一卡二| 成年女人永久免费观看视频| 永久网站在线| 国内精品一区二区在线观看| 亚洲精华国产精华液的使用体验| 乱码一卡2卡4卡精品| 国产视频内射| 亚洲天堂国产精品一区在线| 久久人妻av系列| 欧美成人精品欧美一级黄| 最近2019中文字幕mv第一页| 亚洲婷婷狠狠爱综合网| 欧美一区二区亚洲| 亚洲国产精品专区欧美| 国产不卡一卡二| 日韩欧美精品v在线| 日韩欧美 国产精品| .国产精品久久| 特级一级黄色大片| 亚洲国产日韩欧美精品在线观看| 又黄又爽又刺激的免费视频.| 女的被弄到高潮叫床怎么办| 国产精品蜜桃在线观看| 国产免费一级a男人的天堂| 国产精品国产高清国产av| 午夜免费男女啪啪视频观看| 美女被艹到高潮喷水动态| 国产在视频线在精品| 国产精品蜜桃在线观看| 国产淫片久久久久久久久| 国产成人a∨麻豆精品| 国产高清视频在线观看网站| 国产午夜精品久久久久久一区二区三区| 国产成人午夜福利电影在线观看| 老女人水多毛片| 极品教师在线视频| 精品无人区乱码1区二区| 毛片女人毛片| 校园人妻丝袜中文字幕| 精品久久久久久久久av| 我的女老师完整版在线观看| 免费电影在线观看免费观看| 我的女老师完整版在线观看| 最近最新中文字幕免费大全7| 久久久成人免费电影| 99久国产av精品| 日日干狠狠操夜夜爽| 国产黄色小视频在线观看| 毛片一级片免费看久久久久| 日本欧美国产在线视频| 精品人妻熟女av久视频| 亚洲最大成人手机在线| 最后的刺客免费高清国语| 亚洲va在线va天堂va国产| 美女cb高潮喷水在线观看| 久久国内精品自在自线图片| 中文乱码字字幕精品一区二区三区 | 91在线精品国自产拍蜜月| av卡一久久| 最近手机中文字幕大全| 亚洲av免费高清在线观看| 日韩亚洲欧美综合| 午夜激情福利司机影院| 麻豆乱淫一区二区| 欧美极品一区二区三区四区| 18禁动态无遮挡网站| 亚洲av电影不卡..在线观看| 蜜桃久久精品国产亚洲av| av线在线观看网站| 黄色欧美视频在线观看| 看片在线看免费视频| 亚洲无线观看免费| 国产午夜精品论理片| 天天一区二区日本电影三级| 国产精品久久久久久久电影| 天天躁夜夜躁狠狠久久av| 免费播放大片免费观看视频在线观看 | 美女大奶头视频| 久久精品久久久久久噜噜老黄 | 亚洲av成人精品一二三区| 成人欧美大片| 午夜福利在线观看免费完整高清在| 国产单亲对白刺激| 少妇猛男粗大的猛烈进出视频 | 久久精品国产鲁丝片午夜精品| 欧美一级a爱片免费观看看| 色综合色国产| 神马国产精品三级电影在线观看| 国产精品一及| 国产高清有码在线观看视频| 午夜激情福利司机影院| 欧美另类亚洲清纯唯美| 99热精品在线国产| 色综合色国产| 国产极品天堂在线| 亚洲欧洲日产国产| 亚洲国产最新在线播放| 热99re8久久精品国产| 十八禁国产超污无遮挡网站| 看免费成人av毛片| 国产成人精品久久久久久| eeuss影院久久| 国产一区二区在线av高清观看| 六月丁香七月| 人妻少妇偷人精品九色| 日本三级黄在线观看| 免费av观看视频| 天堂中文最新版在线下载 | 女人十人毛片免费观看3o分钟| 最近的中文字幕免费完整| 亚洲av成人精品一区久久| 亚洲国产精品合色在线| 只有这里有精品99| 精品久久久久久久久av| 亚洲欧美清纯卡通| 成人亚洲精品av一区二区| 伦精品一区二区三区| 亚洲激情五月婷婷啪啪| 色综合亚洲欧美另类图片| 九色成人免费人妻av| 男女视频在线观看网站免费| 国产乱来视频区| 日本三级黄在线观看| videos熟女内射| 少妇熟女aⅴ在线视频| 日韩av在线免费看完整版不卡| 菩萨蛮人人尽说江南好唐韦庄 | 韩国av在线不卡| 成人综合一区亚洲| 欧美日韩综合久久久久久| 最新中文字幕久久久久| 青青草视频在线视频观看| 99热6这里只有精品| 亚洲三级黄色毛片| 白带黄色成豆腐渣| 美女大奶头视频| 亚洲av电影在线观看一区二区三区 | 国产精品一区二区三区四区免费观看| 91午夜精品亚洲一区二区三区| 天堂√8在线中文| 男插女下体视频免费在线播放| 极品教师在线视频| 午夜免费男女啪啪视频观看| 成年女人看的毛片在线观看| 国产精品av视频在线免费观看| 听说在线观看完整版免费高清| 国产成人aa在线观看| 级片在线观看| 欧美最新免费一区二区三区| 精品久久久久久久久亚洲| 综合色av麻豆| 亚洲成人久久爱视频| 老女人水多毛片| 国产黄片美女视频| 日本wwww免费看| av在线播放精品| 国产精品久久久久久久久免| 国产伦一二天堂av在线观看| 欧美日韩在线观看h| 一级二级三级毛片免费看| 国产人妻一区二区三区在| 日韩视频在线欧美| 我的老师免费观看完整版| 又粗又硬又长又爽又黄的视频| 男人的好看免费观看在线视频| 亚洲图色成人| 亚洲经典国产精华液单| 少妇人妻一区二区三区视频| 成人漫画全彩无遮挡| 午夜激情欧美在线| 啦啦啦韩国在线观看视频| 身体一侧抽搐| 99热网站在线观看| 少妇猛男粗大的猛烈进出视频 | 男女边吃奶边做爰视频| 一本久久精品| 国产又黄又爽又无遮挡在线| 国产白丝娇喘喷水9色精品| 精品国产三级普通话版| 久久6这里有精品| 欧美丝袜亚洲另类| 欧美97在线视频| 夜夜看夜夜爽夜夜摸| 亚洲av男天堂| 午夜免费激情av| 在线播放国产精品三级| 人人妻人人澡人人爽人人夜夜 | 国产高潮美女av| 精品久久久久久久久av| 亚洲国产色片| 啦啦啦韩国在线观看视频| 成人毛片60女人毛片免费| 国产精品国产三级国产av玫瑰| 插阴视频在线观看视频| 嘟嘟电影网在线观看| 国产午夜精品论理片| 精华霜和精华液先用哪个| 欧美一区二区国产精品久久精品| 国产在视频线精品| 国产色爽女视频免费观看| 精品国产露脸久久av麻豆 | 在线观看av片永久免费下载| 变态另类丝袜制服| 国产精品美女特级片免费视频播放器| 精品久久久久久久人妻蜜臀av| 日韩一本色道免费dvd| 色综合站精品国产| 99热全是精品| 欧美激情久久久久久爽电影| 五月伊人婷婷丁香| 搡老妇女老女人老熟妇| 日本黄色片子视频| 亚洲欧美成人精品一区二区| 麻豆久久精品国产亚洲av| 亚洲精品乱久久久久久| 成人无遮挡网站| 亚洲国产精品久久男人天堂| 黄色配什么色好看| 黄片无遮挡物在线观看| 日韩强制内射视频| 麻豆乱淫一区二区| 亚洲,欧美,日韩| 2021天堂中文幕一二区在线观| 91久久精品国产一区二区三区| 精品人妻熟女av久视频| 久久久精品大字幕| 你懂的网址亚洲精品在线观看 | 日本黄色视频三级网站网址| 亚洲美女视频黄频| 91久久精品国产一区二区成人| 久热久热在线精品观看| 亚洲在久久综合| 国产一区二区在线观看日韩| 日本爱情动作片www.在线观看| 91精品国产九色| 少妇熟女aⅴ在线视频| 亚洲精品成人久久久久久| 免费观看人在逋| 乱系列少妇在线播放| 欧美高清成人免费视频www| 国语自产精品视频在线第100页| 免费观看a级毛片全部| 久久99热这里只频精品6学生 | 久久久久国产网址| 国产不卡一卡二| 一级毛片我不卡| 成人一区二区视频在线观看| 午夜福利高清视频| 亚洲欧美日韩卡通动漫| 精品人妻一区二区三区麻豆| 精品一区二区三区人妻视频| 亚洲欧洲日产国产| 啦啦啦啦在线视频资源| 岛国毛片在线播放| 久久久成人免费电影| 精品欧美国产一区二区三| 亚洲人成网站在线观看播放| 亚洲精品,欧美精品| 狠狠狠狠99中文字幕| 99久久成人亚洲精品观看| 亚洲av日韩在线播放| 欧美性猛交黑人性爽| 国产69精品久久久久777片| 伦理电影大哥的女人| 亚洲精品乱久久久久久| 成人av在线播放网站| 中国美白少妇内射xxxbb| 91av网一区二区| 能在线免费看毛片的网站| 青春草视频在线免费观看| 久久久亚洲精品成人影院| 99国产精品一区二区蜜桃av| videos熟女内射| 91午夜精品亚洲一区二区三区| 校园人妻丝袜中文字幕| 高清在线视频一区二区三区 | 国产成人福利小说| 久久精品国产99精品国产亚洲性色| 色综合色国产| 免费不卡的大黄色大毛片视频在线观看 | 国产视频首页在线观看| 免费看av在线观看网站| 久久精品国产亚洲av涩爱| 26uuu在线亚洲综合色| 一级毛片我不卡| 国产淫语在线视频| 日韩精品青青久久久久久| 日韩成人av中文字幕在线观看| 亚洲美女搞黄在线观看| 欧美三级亚洲精品| 久久人人爽人人爽人人片va| 少妇人妻精品综合一区二区| 久久久久免费精品人妻一区二区| 五月伊人婷婷丁香| 国语对白做爰xxxⅹ性视频网站| 久久精品国产亚洲av天美| av天堂中文字幕网| 99在线人妻在线中文字幕| 1024手机看黄色片| 高清视频免费观看一区二区 | 国产欧美另类精品又又久久亚洲欧美| 1024手机看黄色片| 美女被艹到高潮喷水动态| 国产精品福利在线免费观看| 色播亚洲综合网| 久久久久久久亚洲中文字幕| 国产爱豆传媒在线观看| 偷拍熟女少妇极品色| 级片在线观看| 成人毛片a级毛片在线播放| 久久久久久久午夜电影| 蜜桃亚洲精品一区二区三区| 成人亚洲欧美一区二区av| 亚洲18禁久久av| 亚洲人与动物交配视频| 高清视频免费观看一区二区 | 国产精品av视频在线免费观看| 中文乱码字字幕精品一区二区三区 | 人妻少妇偷人精品九色| 精品少妇黑人巨大在线播放 | 久久久久久久国产电影| 亚洲中文字幕日韩| 国产国拍精品亚洲av在线观看| 99热6这里只有精品| 人妻制服诱惑在线中文字幕| 爱豆传媒免费全集在线观看| 99热这里只有是精品在线观看| 白带黄色成豆腐渣| 一边亲一边摸免费视频| 色视频www国产| 国产黄片美女视频| 亚洲内射少妇av| 国产v大片淫在线免费观看| 国产精品久久久久久久电影| 亚洲精品国产成人久久av| 欧美极品一区二区三区四区| 日本-黄色视频高清免费观看| 国产亚洲午夜精品一区二区久久 | 美女黄网站色视频| 亚洲欧美中文字幕日韩二区| 欧美人与善性xxx| 99久久精品热视频| 1024手机看黄色片| 日韩成人伦理影院| 精品一区二区免费观看| 久久久色成人| 看十八女毛片水多多多| 亚洲av成人av| 亚洲欧美精品综合久久99| 午夜福利在线观看免费完整高清在| 国产高清国产精品国产三级 | 精品少妇黑人巨大在线播放 | 日韩视频在线欧美| 精华霜和精华液先用哪个| 亚洲国产最新在线播放| 三级国产精品片| 日韩国内少妇激情av| 亚洲av福利一区| 国产91av在线免费观看| 亚洲av日韩在线播放| 99国产精品一区二区蜜桃av| 午夜精品一区二区三区免费看| 国产一区二区三区av在线| 国产亚洲av嫩草精品影院| 日日摸夜夜添夜夜添av毛片| 亚洲综合色惰| 亚洲人与动物交配视频| 欧美最新免费一区二区三区| 亚洲精品国产成人久久av| 亚洲美女搞黄在线观看| 国内少妇人妻偷人精品xxx网站| 综合色丁香网| 亚洲在线自拍视频| 欧美日本亚洲视频在线播放| 午夜免费激情av| 天堂网av新在线| 成年女人看的毛片在线观看| 深爱激情五月婷婷| 老师上课跳d突然被开到最大视频| 欧美又色又爽又黄视频| 最后的刺客免费高清国语| 中国美白少妇内射xxxbb| 国产成人免费观看mmmm| 成年免费大片在线观看| 两个人的视频大全免费| 天天躁日日操中文字幕| 91精品一卡2卡3卡4卡| 国产av不卡久久| 国产精品永久免费网站| 国产精品福利在线免费观看| 男插女下体视频免费在线播放| 熟妇人妻久久中文字幕3abv| 美女黄网站色视频| 国产亚洲一区二区精品| 久久精品91蜜桃| 一二三四中文在线观看免费高清| 成人综合一区亚洲| 国产一区二区在线观看日韩| 熟女人妻精品中文字幕| 网址你懂的国产日韩在线| 亚洲精品国产成人久久av| 久久久久久久亚洲中文字幕| 久久久久精品久久久久真实原创| 最近2019中文字幕mv第一页| 亚洲av二区三区四区| 最近视频中文字幕2019在线8| 亚洲人成网站高清观看| 日本午夜av视频| 中文字幕精品亚洲无线码一区| 国产免费一级a男人的天堂| av卡一久久| 人人妻人人看人人澡| АⅤ资源中文在线天堂| 日本三级黄在线观看| 亚洲欧美日韩无卡精品| 毛片女人毛片| 国产老妇女一区| 校园人妻丝袜中文字幕| 国产激情偷乱视频一区二区| 高清日韩中文字幕在线| 成年免费大片在线观看| 1000部很黄的大片| 日本一二三区视频观看| 白带黄色成豆腐渣| 日韩成人av中文字幕在线观看| 日韩欧美在线乱码| 国产黄a三级三级三级人| 国产综合懂色| 久久亚洲精品不卡| 亚洲一区高清亚洲精品| 在线观看av片永久免费下载| 国产精华一区二区三区| 成年版毛片免费区| 国产精品女同一区二区软件| 成人一区二区视频在线观看| 七月丁香在线播放| 欧美成人免费av一区二区三区| 亚洲欧美精品专区久久| 免费黄色在线免费观看| 精品国产一区二区三区久久久樱花 | 人妻夜夜爽99麻豆av| 欧美成人一区二区免费高清观看| 欧美激情国产日韩精品一区| 男女边吃奶边做爰视频| 又黄又爽又刺激的免费视频.| 美女国产视频在线观看| 国产精品伦人一区二区| 日本五十路高清| 国产中年淑女户外野战色| 大话2 男鬼变身卡| 久久久久九九精品影院| 三级经典国产精品| 日本免费a在线| 国产午夜精品一二区理论片| 国产av不卡久久| 国产色爽女视频免费观看| 草草在线视频免费看| 久热久热在线精品观看| 午夜福利在线在线| 1024手机看黄色片| av.在线天堂| 性色avwww在线观看| 亚洲色图av天堂| 岛国在线免费视频观看| 国产激情偷乱视频一区二区| 国产男人的电影天堂91| 久久久a久久爽久久v久久| 亚洲怡红院男人天堂| 亚洲av电影不卡..在线观看| 少妇人妻一区二区三区视频| 建设人人有责人人尽责人人享有的 | 亚洲性久久影院| 狂野欧美激情性xxxx在线观看| 美女被艹到高潮喷水动态| 欧美精品国产亚洲| 国产成人精品一,二区| 亚洲av.av天堂| 级片在线观看| 草草在线视频免费看| 免费黄网站久久成人精品| 在线观看一区二区三区| 欧美成人免费av一区二区三区| 日本午夜av视频| 中文字幕亚洲精品专区| 免费观看a级毛片全部| 国产精品伦人一区二区| 自拍偷自拍亚洲精品老妇| 国产麻豆成人av免费视频| 男的添女的下面高潮视频| 熟女人妻精品中文字幕| 国产精品一区二区性色av| 国产老妇伦熟女老妇高清| 精品熟女少妇av免费看| 一本久久精品| 人体艺术视频欧美日本| 丰满少妇做爰视频| 老司机影院毛片| 嫩草影院入口| 一级黄色大片毛片| 一区二区三区免费毛片| 老司机影院毛片| 久久午夜福利片| 国产熟女欧美一区二区| 好男人视频免费观看在线| 3wmmmm亚洲av在线观看| 中文字幕精品亚洲无线码一区| 亚洲精品自拍成人| 国产黄色视频一区二区在线观看 | 国产片特级美女逼逼视频| 身体一侧抽搐| 精品一区二区三区视频在线| 国产极品精品免费视频能看的| 性色avwww在线观看| 人体艺术视频欧美日本| 日韩欧美精品免费久久| 你懂的网址亚洲精品在线观看 | 日本免费在线观看一区| 国产精华一区二区三区| 国产一级毛片七仙女欲春2| 好男人在线观看高清免费视频| 精品人妻偷拍中文字幕| 成人性生交大片免费视频hd| 日韩一本色道免费dvd| 国产黄色视频一区二区在线观看 | 桃色一区二区三区在线观看| 人人妻人人澡欧美一区二区| 久久午夜福利片| 国产白丝娇喘喷水9色精品| 男女视频在线观看网站免费| 亚洲成色77777| 久久久久九九精品影院| 精品一区二区三区人妻视频| av女优亚洲男人天堂| 在线观看av片永久免费下载| 免费黄色在线免费观看| 免费看a级黄色片| 一级av片app| 在线观看一区二区三区| 内地一区二区视频在线| 久久久久久久久大av| 日本一本二区三区精品| 高清av免费在线| 亚洲欧美成人精品一区二区| 69人妻影院| 老司机影院毛片| 2021少妇久久久久久久久久久| 国产乱人视频| 国产精品国产三级国产专区5o | 亚洲精品456在线播放app| 久久精品国产自在天天线| 亚洲不卡免费看| 亚州av有码| 欧美变态另类bdsm刘玥| 少妇裸体淫交视频免费看高清| 少妇熟女aⅴ在线视频| 亚洲高清免费不卡视频| 在线天堂最新版资源| 99在线视频只有这里精品首页| 男女啪啪激烈高潮av片| 国产乱来视频区| 亚洲国产精品sss在线观看| 亚洲美女搞黄在线观看| 色吧在线观看| 麻豆一二三区av精品| 国产精品精品国产色婷婷| 亚洲,欧美,日韩| 亚洲国产精品成人综合色| 一级黄色大片毛片| 国产片特级美女逼逼视频| 日产精品乱码卡一卡2卡三| 欧美xxxx性猛交bbbb| 亚洲成人精品中文字幕电影| 老司机福利观看| 亚洲精品乱久久久久久| 美女被艹到高潮喷水动态| 亚洲av电影在线观看一区二区三区 | 国产精品久久久久久精品电影| 久久人妻av系列| 色综合站精品国产| 国产精品乱码一区二三区的特点| 亚洲五月天丁香| 午夜福利在线观看免费完整高清在| 国产精品美女特级片免费视频播放器| av播播在线观看一区| 蜜桃久久精品国产亚洲av| 亚洲av免费在线观看| 日韩一区二区视频免费看| 超碰av人人做人人爽久久| 亚洲经典国产精华液单| 久久国内精品自在自线图片| 91午夜精品亚洲一区二区三区| 嘟嘟电影网在线观看| 免费av毛片视频| 天堂av国产一区二区熟女人妻| 日本免费a在线| 欧美日韩综合久久久久久| 精品一区二区三区视频在线| 国产午夜精品久久久久久一区二区三区| 国产91av在线免费观看| 国产欧美另类精品又又久久亚洲欧美| 波多野结衣巨乳人妻| 成人午夜高清在线视频| 欧美日韩一区二区视频在线观看视频在线 | 国产又黄又爽又无遮挡在线| 亚洲欧洲日产国产| 国产成人freesex在线| 国语对白做爰xxxⅹ性视频网站| 午夜福利高清视频| 婷婷色av中文字幕| 色综合亚洲欧美另类图片| 欧美潮喷喷水| 免费观看人在逋| 色网站视频免费| 最近的中文字幕免费完整| 男人舔女人下体高潮全视频|