• <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.中國海洋大學(xué)海洋與大氣學(xué)院,山東 青島 266100;2.中國海洋大學(xué)物理海洋教育部重點(diǎn)實(shí)驗(yàn)室,山東 青島 266003; 3. 青島海洋科學(xué)與技術(shù)國家實(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è)資源為中國、日本、韓國等國共同開發(fā)利用。藍(lán)點(diǎn)馬鮫是我國產(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)馬鮫群體無遺傳分化的重要原因[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];邱盛堯等討論了我國漁業(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)級(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)層次的關(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ù)和美國國家海洋大氣局(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ì)我國漁業(yè)管理和藍(lán)點(diǎn)馬鮫資源的可持續(xù)利用具有一定的數(shù)據(jù)支撐和指導(dǎo)作用。

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

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

    我國近海藍(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月,但由于我國沿??缭蕉鄠€(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我國沿海藍(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án)點(diǎn)馬鮫作為我國近海最重要的漁業(yè)資源種類之一,雖然已有學(xué)者針對(duì)其生物學(xué)特性及漁業(yè)資源變化做了許多工作,但是在對(duì)環(huán)境的響應(yīng)方面缺乏深入研究。近年來越來越多的研究證明,無論是從單一的海洋生物到整個(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ù)來源

    1992—2012年黃渤海海域藍(lán)點(diǎn)馬鮫捕獲量數(shù)據(jù)及年齡組成資料引自農(nóng)業(yè)部漁業(yè)局編制的《中國漁業(yè)統(tǒng)計(jì)年鑒》和邱盛堯、孫本曉等[1,12-13,34-35]的歷史資料。1992—2012年的海表面溫度采用美國國家海洋大氣局的最優(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)馬鮫群體。我國是藍(lán)點(diǎn)馬鮫漁業(yè)資源的主要開發(fā)國,藍(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),左下角為月份。受太陽輻射影響,在北半球深水大洋中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改變的原因較多。從定義來看,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月)為我國黃渤海藍(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án)點(diǎn)馬鮫洄游模式、年齡與生長(zhǎng)以及資源評(píng)估和管理的研究,開展藍(lán)點(diǎn)馬鮫生活史參數(shù)化方案的研究,建立應(yīng)對(duì)氣候變化的早期預(yù)警機(jī)制等研究。

    国产永久视频网站| 亚洲精品国产色婷婷电影| 成人无遮挡网站| a级毛色黄片| 内地一区二区视频在线| 99re6热这里在线精品视频| 日韩欧美 国产精品| 女人精品久久久久毛片| 视频中文字幕在线观看| 久久人人爽人人片av| 国产深夜福利视频在线观看| 国产成人精品福利久久| 精品国产一区二区久久| 欧美国产精品一级二级三级 | 女人久久www免费人成看片| 少妇人妻一区二区三区视频| 最近中文字幕高清免费大全6| 国产美女午夜福利| 国产午夜精品一二区理论片| 国产永久视频网站| 久久影院123| 亚洲人与动物交配视频| av专区在线播放| 人妻夜夜爽99麻豆av| 丰满人妻一区二区三区视频av| 国产成人免费无遮挡视频| 丝袜在线中文字幕| 高清视频免费观看一区二区| 日日啪夜夜撸| 少妇裸体淫交视频免费看高清| 少妇熟女欧美另类| 尾随美女入室| 亚洲欧美精品专区久久| 国产av国产精品国产| 男人狂女人下面高潮的视频| 亚洲国产日韩一区二区| 免费播放大片免费观看视频在线观看| 在线观看免费高清a一片| 久热久热在线精品观看| 久久久久久久久久久免费av| 看非洲黑人一级黄片| 天堂8中文在线网| 秋霞在线观看毛片| 久久毛片免费看一区二区三区| 两个人的视频大全免费| 91久久精品国产一区二区三区| 婷婷色综合www| 97在线人人人人妻| 国产免费一级a男人的天堂| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产色片| 国产av一区二区精品久久| 人妻夜夜爽99麻豆av| 国产亚洲91精品色在线| 免费黄色在线免费观看| 少妇被粗大的猛进出69影院 | 国产精品一区www在线观看| 欧美精品高潮呻吟av久久| 在线精品无人区一区二区三| 少妇丰满av| 国产av国产精品国产| 丝袜在线中文字幕| 黄色日韩在线| 桃花免费在线播放| 老司机影院成人| 97在线人人人人妻| 中文天堂在线官网| 9色porny在线观看| 国产精品久久久久久精品古装| 国产精品国产av在线观看| 亚洲精品亚洲一区二区| 久久久久久久精品精品| 黑人巨大精品欧美一区二区蜜桃 | 欧美变态另类bdsm刘玥| 91aial.com中文字幕在线观看| 国产视频首页在线观看| 少妇人妻精品综合一区二区| 国产淫语在线视频| 高清黄色对白视频在线免费看 | 另类亚洲欧美激情| 黄片无遮挡物在线观看| 成人18禁高潮啪啪吃奶动态图 | av线在线观看网站| 亚洲欧美一区二区三区黑人 | 欧美+日韩+精品| 国产一区二区三区av在线| 天美传媒精品一区二区| 亚洲精品国产av蜜桃| 免费大片黄手机在线观看| av免费观看日本| 国精品久久久久久国模美| 国产精品伦人一区二区| 成人亚洲精品一区在线观看| 一本—道久久a久久精品蜜桃钙片| 最近手机中文字幕大全| a级一级毛片免费在线观看| 久久久亚洲精品成人影院| 有码 亚洲区| 国产片特级美女逼逼视频| av专区在线播放| 最新的欧美精品一区二区| 我要看日韩黄色一级片| 丝瓜视频免费看黄片| 色网站视频免费| 久久久久视频综合| 狂野欧美激情性bbbbbb| 中文资源天堂在线| 精品国产一区二区久久| 欧美精品亚洲一区二区| 日韩电影二区| 91aial.com中文字幕在线观看| 亚洲,一卡二卡三卡| 99久久人妻综合| 99精国产麻豆久久婷婷| 欧美国产精品一级二级三级 | 亚洲怡红院男人天堂| 亚洲综合精品二区| 另类亚洲欧美激情| 精品国产一区二区久久| 最近2019中文字幕mv第一页| 大香蕉久久网| 色婷婷av一区二区三区视频| 亚洲欧美一区二区三区国产| 国产精品.久久久| www.色视频.com| 天堂中文最新版在线下载| 亚洲国产欧美在线一区| a级毛色黄片| 黄色配什么色好看| 国产男人的电影天堂91| 亚洲精品中文字幕在线视频 | 久久影院123| 中文字幕亚洲精品专区| 国产精品一区二区性色av| 2018国产大陆天天弄谢| 晚上一个人看的免费电影| 夫妻性生交免费视频一级片| av一本久久久久| 婷婷色av中文字幕| 夫妻性生交免费视频一级片| 中文精品一卡2卡3卡4更新| 在线免费观看不下载黄p国产| 伦精品一区二区三区| 搡女人真爽免费视频火全软件| 亚洲欧美清纯卡通| 免费少妇av软件| 亚洲精品久久久久久婷婷小说| 色视频在线一区二区三区| 久久久久久人妻| 全区人妻精品视频| 欧美少妇被猛烈插入视频| 性色av一级| 国产精品一区二区三区四区免费观看| 一边亲一边摸免费视频| 国产91av在线免费观看| 日韩av不卡免费在线播放| 男女边吃奶边做爰视频| 久久精品国产亚洲网站| 国产黄片美女视频| av一本久久久久| 三上悠亚av全集在线观看 | 少妇熟女欧美另类| 夜夜看夜夜爽夜夜摸| 一级毛片电影观看| 国产日韩欧美亚洲二区| a级毛片在线看网站| 色视频在线一区二区三区| av黄色大香蕉| 亚洲欧美精品专区久久| 老司机影院成人| 尾随美女入室| 大码成人一级视频| 亚洲精品国产av蜜桃| 亚洲av日韩在线播放| 久久午夜综合久久蜜桃| 99热这里只有是精品50| 婷婷色麻豆天堂久久| 日日摸夜夜添夜夜爱| av在线老鸭窝| av在线观看视频网站免费| 18禁在线无遮挡免费观看视频| 男女无遮挡免费网站观看| 在线观看免费日韩欧美大片 | 丁香六月天网| 久久国产乱子免费精品| 最近手机中文字幕大全| 成年人免费黄色播放视频 | 免费黄网站久久成人精品| 在线看a的网站| 日韩强制内射视频| 国产精品蜜桃在线观看| 男女边摸边吃奶| 卡戴珊不雅视频在线播放| 夜夜爽夜夜爽视频| 人人妻人人澡人人看| 久久久久网色| 老司机影院毛片| 精品久久久久久久久av| 91成人精品电影| 久久99精品国语久久久| 亚洲国产精品成人久久小说| 一个人看视频在线观看www免费| 女人精品久久久久毛片| 在线亚洲精品国产二区图片欧美 | 久久人人爽人人爽人人片va| 免费黄网站久久成人精品| 成人毛片a级毛片在线播放| 日本爱情动作片www.在线观看| 亚洲av成人精品一区久久| a级毛片在线看网站| 亚洲精华国产精华液的使用体验| 亚洲av成人精品一二三区| 国产精品一区www在线观看| 欧美激情极品国产一区二区三区 | 亚洲综合色惰| 妹子高潮喷水视频| 日本黄大片高清| 中文精品一卡2卡3卡4更新| 午夜福利影视在线免费观看| 丝袜在线中文字幕| 秋霞伦理黄片| freevideosex欧美| 大话2 男鬼变身卡| 成人国产av品久久久| 国产精品三级大全| 简卡轻食公司| 久久久久久久久久久免费av| 日韩成人av中文字幕在线观看| 99热这里只有是精品50| 99久久精品热视频| 国产一级毛片在线| 汤姆久久久久久久影院中文字幕| 少妇裸体淫交视频免费看高清| 欧美精品亚洲一区二区| 国产欧美日韩精品一区二区| 国产精品.久久久| av在线观看视频网站免费| 欧美少妇被猛烈插入视频| 欧美精品亚洲一区二区| h日本视频在线播放| 多毛熟女@视频| 伊人久久国产一区二区| 日本与韩国留学比较| 婷婷色综合www| 国产成人freesex在线| 国产成人精品婷婷| 老女人水多毛片| 久久国产乱子免费精品| 麻豆乱淫一区二区| 伊人亚洲综合成人网| 男男h啪啪无遮挡| 亚洲精华国产精华液的使用体验| 国产精品免费大片| 乱码一卡2卡4卡精品| 深夜a级毛片| av女优亚洲男人天堂| av有码第一页| 久久99热这里只频精品6学生| h视频一区二区三区| 精品久久久久久电影网| 国产日韩欧美视频二区| 日韩一本色道免费dvd| 插阴视频在线观看视频| 免费av不卡在线播放| 18禁在线播放成人免费| 在线观看免费视频网站a站| 日韩精品有码人妻一区| 蜜臀久久99精品久久宅男| 成人影院久久| 国产日韩欧美在线精品| 日韩制服骚丝袜av| 丝袜在线中文字幕| 免费人妻精品一区二区三区视频| a 毛片基地| 人体艺术视频欧美日本| 高清不卡的av网站| 久久国产精品男人的天堂亚洲 | 成年人免费黄色播放视频 | 夜夜爽夜夜爽视频| 18禁裸乳无遮挡动漫免费视频| 女的被弄到高潮叫床怎么办| 亚洲美女黄色视频免费看| 我要看黄色一级片免费的| 国产精品一区二区三区四区免费观看| 精品少妇久久久久久888优播| 精品一品国产午夜福利视频| 久久人人爽人人爽人人片va| 免费少妇av软件| 亚洲怡红院男人天堂| 另类精品久久| 91午夜精品亚洲一区二区三区| 欧美 亚洲 国产 日韩一| 久久久国产精品麻豆| 国产日韩一区二区三区精品不卡 | 99re6热这里在线精品视频| 欧美少妇被猛烈插入视频| 国产成人精品一,二区| 日韩不卡一区二区三区视频在线| 中文字幕制服av| 汤姆久久久久久久影院中文字幕| 亚洲欧美日韩另类电影网站| 一级,二级,三级黄色视频| 一本大道久久a久久精品| 亚洲中文av在线| 在线观看美女被高潮喷水网站| 另类亚洲欧美激情| 伦理电影免费视频| 18+在线观看网站| 欧美精品国产亚洲| 国产色婷婷99| 汤姆久久久久久久影院中文字幕| 热re99久久国产66热| 涩涩av久久男人的天堂| 高清黄色对白视频在线免费看 | 日本欧美视频一区| 91久久精品电影网| 五月玫瑰六月丁香| 春色校园在线视频观看| 国产免费又黄又爽又色| freevideosex欧美| 视频中文字幕在线观看| 性色av一级| 精品酒店卫生间| 亚洲av福利一区| 免费在线观看成人毛片| 伊人久久精品亚洲午夜| 国产一区二区在线观看日韩| 一区二区三区免费毛片| 精品视频人人做人人爽| 黄片无遮挡物在线观看| 精品卡一卡二卡四卡免费| av在线老鸭窝| 国产亚洲5aaaaa淫片| 美女福利国产在线| 国产日韩欧美在线精品| 国产精品福利在线免费观看| 国产免费一区二区三区四区乱码| 久久这里有精品视频免费| 搡女人真爽免费视频火全软件| 国产成人aa在线观看| 黄色一级大片看看| 美女xxoo啪啪120秒动态图| av视频免费观看在线观看| 2018国产大陆天天弄谢| a级一级毛片免费在线观看| 国产精品免费大片| 午夜视频国产福利| 一级毛片电影观看| 91精品一卡2卡3卡4卡| 亚洲欧洲日产国产| 国产成人午夜福利电影在线观看| 亚洲欧洲精品一区二区精品久久久 | 麻豆成人午夜福利视频| 国内少妇人妻偷人精品xxx网站| 五月天丁香电影| 寂寞人妻少妇视频99o| 高清欧美精品videossex| 在线天堂最新版资源| 插阴视频在线观看视频| 高清午夜精品一区二区三区| 我要看黄色一级片免费的| 欧美日韩综合久久久久久| 精品卡一卡二卡四卡免费| 欧美最新免费一区二区三区| 欧美日韩国产mv在线观看视频| 看免费成人av毛片| 91久久精品国产一区二区成人| 在线精品无人区一区二区三| 日日爽夜夜爽网站| 夫妻性生交免费视频一级片| 亚洲国产最新在线播放| 国产精品久久久久久久电影| 亚洲四区av| 欧美三级亚洲精品| 久久久久精品性色| 91久久精品电影网| 人妻系列 视频| 日产精品乱码卡一卡2卡三| 在线免费观看不下载黄p国产| 亚洲精品视频女| 亚洲精品成人av观看孕妇| 国内少妇人妻偷人精品xxx网站| 水蜜桃什么品种好| 自拍偷自拍亚洲精品老妇| 国产淫片久久久久久久久| 久久女婷五月综合色啪小说| 中文字幕精品免费在线观看视频 | 人妻制服诱惑在线中文字幕| 久久 成人 亚洲| 国产精品.久久久| 日本爱情动作片www.在线观看| av专区在线播放| 亚洲天堂av无毛| 嫩草影院入口| 涩涩av久久男人的天堂| 一级,二级,三级黄色视频| 色视频在线一区二区三区| 观看免费一级毛片| 男女免费视频国产| 91在线精品国自产拍蜜月| 在线观看免费视频网站a站| 国产真实伦视频高清在线观看| 美女福利国产在线| 亚洲欧美成人精品一区二区| 99久久中文字幕三级久久日本| 99热全是精品| 天堂中文最新版在线下载| 卡戴珊不雅视频在线播放| 午夜福利,免费看| 午夜激情久久久久久久| 国产无遮挡羞羞视频在线观看| 午夜福利网站1000一区二区三区| 一级毛片我不卡| 少妇高潮的动态图| 成人漫画全彩无遮挡| 欧美精品高潮呻吟av久久| 日本wwww免费看| 欧美成人精品欧美一级黄| 日韩av在线免费看完整版不卡| 我要看黄色一级片免费的| 亚洲中文av在线| 欧美日韩视频高清一区二区三区二| 日韩电影二区| 99久久中文字幕三级久久日本| 狠狠精品人妻久久久久久综合| 国产精品福利在线免费观看| 曰老女人黄片| 欧美三级亚洲精品| av国产久精品久网站免费入址| 一个人看视频在线观看www免费| 91久久精品国产一区二区三区| 精品亚洲乱码少妇综合久久| 亚洲av男天堂| 22中文网久久字幕| 亚洲情色 制服丝袜| 国产视频内射| 在线看a的网站| 夜夜骑夜夜射夜夜干| av天堂中文字幕网| 香蕉精品网在线| 男人舔奶头视频| 亚洲精品自拍成人| 国产美女午夜福利| 国产永久视频网站| 亚洲国产精品成人久久小说| 99精国产麻豆久久婷婷| 成人黄色视频免费在线看| 秋霞在线观看毛片| 亚洲av日韩在线播放| 我的老师免费观看完整版| 日本欧美国产在线视频| 免费观看av网站的网址| 国产免费一区二区三区四区乱码| 久久99热6这里只有精品| 欧美日韩视频精品一区| 在线免费观看不下载黄p国产| 边亲边吃奶的免费视频| 亚洲欧美清纯卡通| 人体艺术视频欧美日本| 精品少妇内射三级| 久久鲁丝午夜福利片| 久久综合国产亚洲精品| a级片在线免费高清观看视频| 丰满乱子伦码专区| 在线观看免费视频网站a站| 日本猛色少妇xxxxx猛交久久| 国产免费又黄又爽又色| 日韩大片免费观看网站| 亚洲av在线观看美女高潮| 人人妻人人澡人人看| 18禁动态无遮挡网站| 3wmmmm亚洲av在线观看| 在线观看www视频免费| 午夜日本视频在线| 麻豆乱淫一区二区| 国产黄色免费在线视频| 少妇人妻精品综合一区二区| 国内少妇人妻偷人精品xxx网站| 国产片特级美女逼逼视频| 成年美女黄网站色视频大全免费 | 亚洲欧洲日产国产| 国产欧美日韩精品一区二区| 少妇的逼好多水| 秋霞伦理黄片| 亚洲欧洲日产国产| 99久久精品一区二区三区| 九色成人免费人妻av| 日本欧美视频一区| 亚洲av福利一区| 久久久a久久爽久久v久久| 亚洲精品乱久久久久久| 亚洲国产毛片av蜜桃av| 99热网站在线观看| 久久久久久人妻| 免费播放大片免费观看视频在线观看| 妹子高潮喷水视频| 免费高清在线观看视频在线观看| 免费黄色在线免费观看| 免费人妻精品一区二区三区视频| 国产乱来视频区| 亚洲av中文av极速乱| 国产淫片久久久久久久久| 一级毛片黄色毛片免费观看视频| 久久青草综合色| 欧美老熟妇乱子伦牲交| 国产亚洲欧美精品永久| 中文在线观看免费www的网站| 26uuu在线亚洲综合色| 一级av片app| 69精品国产乱码久久久| 国产精品久久久久久久久免| 少妇人妻一区二区三区视频| 特大巨黑吊av在线直播| 男女边摸边吃奶| 亚洲综合色惰| 亚洲不卡免费看| 3wmmmm亚洲av在线观看| 国产视频内射| 国产午夜精品久久久久久一区二区三区| 黄色视频在线播放观看不卡| 全区人妻精品视频| 91精品伊人久久大香线蕉| 人妻 亚洲 视频| 深夜a级毛片| 99热全是精品| 天天躁夜夜躁狠狠久久av| 国模一区二区三区四区视频| 久久久久久久久久成人| 日韩av免费高清视频| 黄色欧美视频在线观看| 美女主播在线视频| 国产精品蜜桃在线观看| 国产免费福利视频在线观看| 丰满少妇做爰视频| 亚洲成人av在线免费| 六月丁香七月| 日韩一区二区视频免费看| 亚洲精品色激情综合| 久久久久久久国产电影| 精品久久久久久电影网| 晚上一个人看的免费电影| 精品熟女少妇av免费看| av天堂久久9| 久久这里有精品视频免费| 18+在线观看网站| 国产精品秋霞免费鲁丝片| 日韩一区二区三区影片| 伊人亚洲综合成人网| 日韩免费高清中文字幕av| 看免费成人av毛片| 青青草视频在线视频观看| 少妇的逼水好多| 国产精品国产三级国产av玫瑰| 国产亚洲一区二区精品| 亚洲无线观看免费| 各种免费的搞黄视频| av线在线观看网站| 亚洲欧美精品专区久久| 国产片特级美女逼逼视频| 高清午夜精品一区二区三区| 亚洲成人手机| 五月开心婷婷网| 大片电影免费在线观看免费| 国产在线一区二区三区精| 久久精品国产自在天天线| 国产视频内射| 夜夜看夜夜爽夜夜摸| 久久午夜福利片| 亚洲av日韩在线播放| 六月丁香七月| 女性生殖器流出的白浆| 青青草视频在线视频观看| av视频免费观看在线观看| 少妇的逼好多水| 大又大粗又爽又黄少妇毛片口| 少妇被粗大猛烈的视频| 亚洲精品成人av观看孕妇| 国产精品一区二区性色av| 卡戴珊不雅视频在线播放| 人人妻人人爽人人添夜夜欢视频 | 男人舔奶头视频| 我要看黄色一级片免费的| av天堂中文字幕网| 啦啦啦中文免费视频观看日本| 亚洲,欧美,日韩| 青春草亚洲视频在线观看| 少妇的逼水好多| 欧美另类一区| 久久狼人影院| 亚洲激情五月婷婷啪啪| 国产黄色视频一区二区在线观看| a级毛片在线看网站| 久久精品夜色国产| 成人毛片60女人毛片免费| 夜夜骑夜夜射夜夜干| 国产精品久久久久久久电影| 午夜福利视频精品| 国产成人精品福利久久| 在线观看免费高清a一片| 伦理电影大哥的女人| 卡戴珊不雅视频在线播放| 人妻制服诱惑在线中文字幕| 人人妻人人添人人爽欧美一区卜| 高清视频免费观看一区二区| 成人18禁高潮啪啪吃奶动态图 | 日产精品乱码卡一卡2卡三| 男人爽女人下面视频在线观看| 大话2 男鬼变身卡| 热99国产精品久久久久久7| 国产有黄有色有爽视频| 亚洲国产精品成人久久小说| 免费人妻精品一区二区三区视频| 纯流量卡能插随身wifi吗| 亚州av有码| 久久亚洲国产成人精品v|