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

    捕撈壓力和氣候變化對黃渤海小黃魚漁獲量的影響*

    2017-06-19 19:25:27劉笑笑徐賓鐸任一平
    關(guān)鍵詞:黃渤海漁獲量小黃魚

    劉笑笑, 王 晶, 徐賓鐸, 薛 瑩, 任一平

    (中國海洋大學(xué)水產(chǎn)學(xué)院,山東 青島 266003)

    捕撈壓力和氣候變化對黃渤海小黃魚漁獲量的影響*

    劉笑笑, 王 晶, 徐賓鐸, 薛 瑩**, 任一平

    (中國海洋大學(xué)水產(chǎn)學(xué)院,山東 青島 266003)

    捕撈壓力和氣候變化是影響漁業(yè)資源數(shù)量變動的主要因素,會對漁獲量產(chǎn)生較大影響。本文通過分析長時(shí)間序列的漁業(yè)統(tǒng)計(jì)資料和氣候變化數(shù)據(jù),研究了1962—2012年間捕撈壓力和氣候變化對黃渤海小黃魚(Pseudosciaenapolyactis)漁獲量的影響,并應(yīng)用Fox模型擬合捕撈壓力對黃渤海小黃魚漁獲量的影響,同時(shí)根據(jù)多個(gè)氣候因子及氣候指數(shù)的年間變動數(shù)據(jù),分析每個(gè)顯著相關(guān)的氣候變量對黃渤海小黃魚漁獲量的影響。研究表明,黃渤海小黃魚漁獲量在1962—1971年呈波動下降趨勢,在1972—1990年保持平穩(wěn)狀態(tài),在1991—2012年急劇增加。黃渤海小黃魚漁獲量與渤海冬季季風(fēng)、北太平洋指數(shù)(North pacific index,NPI)呈顯著負(fù)相關(guān)(P<0.05),與黃海夏季季風(fēng)、黃海海表溫度呈顯著正相關(guān)(P<0.05)。加入氣候變量進(jìn)行優(yōu)化的Fox模型考慮了捕撈壓力和氣候變化的影響,該模型擬合的漁獲量與統(tǒng)計(jì)的漁獲量呈極顯著相關(guān)(P<0.01)。研究表明:捕撈壓力的增大是導(dǎo)致黃渤海小黃魚漁獲量變動的主要原因,氣候變化也會對小黃魚漁獲量產(chǎn)生顯著影響。

    黃渤海;小黃魚漁獲量;Fox模型;捕撈壓力;氣候變化

    小黃魚(Pseudosciaenapolyactis)隸屬于鱸形目(Perciformes)石首魚科(Sciaenidae)黃魚屬(Pseudosciaena)。分布范圍較廣,主要集中在渤海、黃海、東海和朝鮮半島西岸海區(qū),水深小于100 m的海域,特別是受到長江徑流影響的黃海南部和東海北部水深為40~80 m的海域[1]。自1850年代以來,小黃魚在我國底拖網(wǎng)漁業(yè)中始終占有重要地位。然而,隨著漁業(yè)裝備和捕撈技術(shù)的不斷發(fā)展,小黃魚資源受到過度捕撈的嚴(yán)重威脅,加之環(huán)境惡化等因素的影響,致使小黃魚漁獲量在1952—1962年期間急劇下降[2-4],漁業(yè)種群也呈現(xiàn)諸如低齡化、個(gè)體生長加快、性成熟提前等生物學(xué)特征的變化[5]。在1962—2012年間,黃渤海小黃魚漁獲量經(jīng)歷了急劇衰退,進(jìn)而平穩(wěn)發(fā)展,隨后有所恢復(fù),并快速上升的過程,小黃魚漁獲量的巨大波動已經(jīng)引起了越來越多科學(xué)家的重視。

    聯(lián)合國糧食及農(nóng)業(yè)組織(Food and Agriculture Organization of the United,F(xiàn)AO)有關(guān)專家認(rèn)為,氣候變化不僅影響海洋物種的分布格局,甚至改變了海洋生物的生理過程,這些變化終將對漁業(yè)資源產(chǎn)生無法預(yù)估的影響[6]。Chang等[7]研究表明,海表水溫、降雨以及氣壓等氣候環(huán)境的變化對韓國近岸海域的漁業(yè)資源量以及海洋生態(tài)系統(tǒng)都產(chǎn)生了相當(dāng)明顯的影響。王躍中等[8]在研究中發(fā)現(xiàn),我國東海帶魚漁獲量變動不僅與捕撈壓力的變化有關(guān),還與海表水溫、熱帶氣旋、冬夏季季風(fēng)、海表徑流以及陸地降雨等氣候因素的變化相關(guān)。Tian等[9]研究發(fā)現(xiàn),太平洋秋刀魚的資源波動受ENSO的影響并不明顯,卻與北太平洋指數(shù)NPI(North pacific index)顯著相關(guān)。上述研究均表明,分析漁獲量的變動規(guī)律不僅要考慮捕撈壓力的影響,還應(yīng)結(jié)合相應(yīng)時(shí)間段內(nèi)的環(huán)境和氣候因子來解釋[10]。

    本研究基于1962—2012年《中國漁業(yè)統(tǒng)計(jì)年鑒》數(shù)據(jù),并結(jié)合長時(shí)間序列的氣候變化數(shù)據(jù),分析捕撈壓力和氣候變化對黃渤海小黃魚漁獲量的影響,旨在探討黃渤海小黃魚漁獲量變動的原因和機(jī)制,為我國小黃魚漁業(yè)資源的管理和可持續(xù)利用提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 數(shù)據(jù)來源

    1.1.1 漁業(yè)統(tǒng)計(jì)資料 漁業(yè)統(tǒng)計(jì)資料主要來源于中國農(nóng)業(yè)部漁業(yè)局編制的《中國漁業(yè)統(tǒng)計(jì)年鑒》(1949—2012年)[11-12]。分析的數(shù)據(jù)主要包括1962—2012年間我國黃渤海沿岸三省一市(遼寧省、河北省、山東省和天津市)的小黃魚漁獲量和捕撈努力量數(shù)據(jù)。本研究采用的捕撈努力量為漁船功率,由機(jī)動漁船和非機(jī)動漁船相加而得,單位統(tǒng)一用千瓦(kW)來表示,其中非機(jī)動漁船捕撈努力量是按其單位捕撈努力量漁獲量(CPUE)的年間變動趨勢與該海區(qū)相同功率的機(jī)動漁船比例標(biāo)準(zhǔn)化而來[13]。

    1.1.2 氣候數(shù)據(jù)資料 根據(jù)國內(nèi)外有關(guān)氣候因素對漁獲量影響的研究,選擇黃渤海冬季季風(fēng)、黃渤海夏季季風(fēng)、黃渤海海表溫度、黃渤海沿岸降雨量、太平洋年代際濤動(Pacific Decadal Oscillation,PDO)、北太平洋指數(shù)等氣候數(shù)據(jù)擬合小黃魚漁獲量。其中,降雨量數(shù)據(jù)來源于Climate Explorer(http://climexp.knmi.nl)網(wǎng)站,基于Global Precipitation Climatology Center(GPCC)V4,1°×1°的網(wǎng)格化數(shù)據(jù)。季風(fēng)數(shù)據(jù)來自于Climate Explorer(http://climexp.knmi.nl)網(wǎng)站,獲取每月2°×2°COADS(Comprehensive Ocean-Atmospheric Data Set)海面風(fēng)速,夏季季風(fēng)數(shù)據(jù)來自于每年6—8月海表平均風(fēng)速,冬季季風(fēng)數(shù)據(jù)來自于每年10月—次年3月海表平均風(fēng)速[14]。海表溫度(SST)同樣來源于Climate Explorer(http://climexp.knmi.nl)網(wǎng)站,基于每月的2°×2°COADS(Comprehensive Ocean-Atmospheric Data Set)海表水溫序列數(shù)據(jù)。NPI和PDO數(shù)據(jù)均來源于日本氣象廳發(fā)布的氣象數(shù)據(jù)。

    為了提高氣候數(shù)據(jù)的精確度,將黃海(32°N~38°N,120°E~126°E)和渤海(36°N~40°N,115°E~122°E)的氣候數(shù)據(jù)分別計(jì)算,并把氣候數(shù)據(jù)全部經(jīng)過Z-score標(biāo)準(zhǔn)化進(jìn)行無量綱化處理,之后再對氣候數(shù)據(jù)進(jìn)行K-S檢驗(yàn)(Kolmogorov-Smirnov Test)顯著性水平均大于0.05,即所得數(shù)據(jù)均符合正態(tài)分布,可以直接應(yīng)用于小黃魚漁獲量的回歸性分析及其他相關(guān)計(jì)算。

    1.2 分析方法

    1.2.1 分析捕撈壓力對小黃魚漁獲量的影響 利用Fox模型分析捕撈壓力對黃渤海小黃魚漁獲量的影響[15]。

    Fox模型為:

    Ye=Cfee-dfe。

    (1)

    式中:Ye為黃渤海小黃魚漁獲量;fe為黃渤海漁船的捕撈努力量;C、d為參數(shù)。

    Fox模型中計(jì)算CPUE公式如下:

    (2)

    公式(2)兩邊取自然對數(shù):

    (3)

    參數(shù)C、d的計(jì)算和置信區(qū)間采用paired data bootstrap方法估計(jì)。具體步驟:(1)對1962—2012年小黃魚漁獲量和捕撈努力量41對數(shù)據(jù)進(jìn)行有放回式隨機(jī)取樣,每次抽取1 000組數(shù)據(jù),得到新的1組1 000對數(shù)據(jù);(2)根據(jù)公式(3)估計(jì)模型參數(shù)C和d;3)重復(fù)上述1、2步驟10 000次,得到參數(shù)C、d的分布情況[16]。得到參數(shù)分布后,代入公式(1),計(jì)算捕撈漁獲量Fox模型估計(jì)值Ye。

    1.2.2 分析氣候變化對小黃魚漁獲量的影響 應(yīng)用優(yōu)化的Fox模型分析捕撈壓力對黃渤海小黃魚漁獲量的影響[16]:

    Yo=Cfee-dfe+a1x1+a2x2……+anxn+b。

    (4)

    式中:Yo為黃渤海小黃魚擬合漁獲量;x1,x2,……,xn為氣候因子;a1,a2,……,an為氣候因子的偏相關(guān)系數(shù),b為常數(shù)。

    首先移除捕撈效應(yīng)引起的漁獲量變動ΔY=Y-Ye,ΔY式中為氣候變化引起的小黃魚漁獲量變動,Y為實(shí)際捕撈量,Ye為Fox模型計(jì)算所得小黃魚漁獲量。

    ΔY的年際變化與降雨、季風(fēng)、海表溫度、氣候指數(shù)等因子相關(guān)。將氣候數(shù)據(jù)代入公式ΔY=a1x1+a2x2+……+anxn+b,利用SPSS軟件對漁獲量變動對各氣候數(shù)據(jù)進(jìn)行線性回歸分析。在回歸分析中,各氣候因子之前的偏相關(guān)系數(shù)用來說明其對小黃魚漁獲量變動的影響程度。

    考慮到氣候因子對漁獲量的影響具有時(shí)滯性,在線性回歸分析中,通過對小黃魚漁獲量前0~5年的氣候因子逐年分析,找到每個(gè)因子對應(yīng)的最適時(shí)滯年限,從而建立優(yōu)化的Fox模型,擬合小黃魚漁獲量的變動。

    2 結(jié)果

    2.1 漁獲量變動與捕撈努力量的關(guān)系

    1962—2012年黃渤海小黃魚漁獲量和漁船功率變化如圖1所示。從黃渤海小黃魚歷年漁獲量變化可以看出,其漁獲量存在明顯的年際變化。1962—1971年小黃魚漁獲量呈波動下降趨勢,1972年漁獲量僅有0.41×104t;1972—1990年小黃魚漁獲量保持平穩(wěn)狀態(tài);在1991—2012年間隨著漁船功率的持續(xù)增長,小黃魚漁獲量也開始上升,2010年小黃魚漁獲量最高達(dá)26.7×104t。從黃渤海小黃魚CPUE(t/kW)年際變化可以看出(見圖2),雖然捕撈壓力增加導(dǎo)致黃渤海小黃魚漁獲量增加,但黃渤海小黃魚CPUE在1962—1972年間急劇下降,在1973—1989年間保持平穩(wěn),隨后有所回升。

    圖1 1962—2012年黃渤海小黃魚漁獲量和捕撈努力量的年間變化

    圖2 1962—2012年黃渤海小黃魚CPUE的年際變化

    采用Bootstrap方法對41組漁獲量和捕撈努力量數(shù)據(jù)在R軟件(R-3.3.1)中進(jìn)行分析得到C=0.035 16,d=-0.000 121,所以小黃魚的Fox模型為:

    Ye=0.035 16fee0.000 121fe。

    在SPSS軟件中,用Fox模型根據(jù)黃渤海捕撈努力量擬合小黃魚漁獲量估計(jì)值,擬合結(jié)果顯示,黃渤海小黃魚漁獲量與捕撈努力量的Fox模型回歸系數(shù)達(dá)0.87,且模型具有顯著性(P<0.01),說明黃渤海小黃魚漁獲量與捕撈努力量之間具有顯著的相關(guān)性。

    2.2 漁獲量變動與氣候因子的關(guān)系

    應(yīng)用Fox模型移除捕撈壓力引起的變動,由相關(guān)氣候因素引起的黃渤海小黃魚漁獲量年間變動見圖3。將黃渤海小黃魚漁獲量年間變動與相關(guān)氣候變量進(jìn)行多元線形回歸分析,相關(guān)氣候變量包括黃渤海沿岸降雨量、黃海海表溫度、渤海海表溫度、黃海冬季季風(fēng)、黃海夏季季風(fēng)、渤海冬季季風(fēng)、渤海夏季季風(fēng)、NPI指數(shù)和PDO指數(shù)(見圖4)。

    圖3 移除捕撈壓力影響的黃渤海小黃魚漁獲量的年間變動

    考慮到氣候影響的時(shí)滯性,對移除捕撈壓力影響后漁獲量的前0~5年的氣候因子變動趨勢逐年進(jìn)行線性回歸分析,選出回歸模型中各因子具有顯著相關(guān)性,并且優(yōu)化后的模型與漁獲量統(tǒng)計(jì)值具最大方差值的氣候因子。結(jié)果表明,黃渤海小黃魚漁獲量與渤海冬季季風(fēng)、北太平洋指數(shù)NPI呈顯著負(fù)相關(guān)(P<0.05),與黃海夏季季風(fēng)、黃海海表溫度呈顯著正相關(guān)(P<0.05),與黃渤海沿岸降雨、黃海冬季季風(fēng)、渤海夏季季風(fēng)、渤海表面溫度、太平洋年代際濤動PDO沒有顯著相關(guān)性(P>0.05)。

    比較發(fā)現(xiàn),時(shí)滯為1年時(shí)具有最大方差值(R2=0.537,P<0.01)(見表1),因此選擇時(shí)滯1年時(shí)的氣候因子進(jìn)行線性回歸分析,擬合出氣候因子對黃渤海小黃魚漁獲量的影響。通過計(jì)算各氣候因子的偏相關(guān)系數(shù)和顯著性水平(見表2),得到優(yōu)化后的黃渤海小黃魚漁獲量Fox模型為:

    Yo=0.035 16fee0.000 121fe+1.453xai-1+1.403xbi-1-

    0.126xci-1-1.695xdi-1+1.099。

    式中:Yo為黃渤海小黃魚擬合漁獲量;fe為黃渤海海洋捕撈漁船捕撈努力量;xa為黃海夏季季風(fēng);xb為黃海海表溫度;xc為NPI;xd為渤海冬季季風(fēng)。

    2.3 漁獲量擬合結(jié)果

    黃渤海小黃魚漁獲量的年間變動,一方面由于捕撈壓力引起,另一方面由于氣候變化引起。根據(jù)優(yōu)化后的Fox模型代入黃渤海漁船捕撈努力量和氣候變量數(shù)據(jù),計(jì)算出擬合的黃渤海區(qū)小黃魚漁獲量(見圖5),擬合結(jié)果表示,黃渤海區(qū)小黃魚漁獲量擬合值與統(tǒng)計(jì)值之間的回歸系數(shù)為0.90,其置信水平達(dá)99%以上。

    圖4 1957—2012年黃渤海主要?dú)夂蛞蜃拥淖兓厔?/p>

    時(shí)滯/aTimelagsig.R200.0010.49010.0000.53720.0010.48230.0020.44040.0020.43950.0150.371

    表2 移除捕撈壓力影響的小黃魚漁獲量變動與氣候因子多元線性回歸方程的偏相關(guān)性分析(時(shí)滯為1年)

    Note:①Climatic variables;②Partial correlation coefficient;③North Pacific Index NPI;④Winter wend speed in the Bohai Sea;⑤Summer wend speed in the Yellow Sea;⑥Sea surface temperature in the Yellow Sea.

    圖5 1962—2012年黃渤海小黃魚實(shí)際漁獲量與擬合漁獲量

    3 討論

    3.1 捕撈壓力對黃渤海小黃魚漁獲量的影響

    本研究應(yīng)用Fox模型擬合捕撈壓力對黃渤海小黃魚漁獲量的影響,發(fā)現(xiàn)小黃魚漁獲量與捕撈努力量的Fox模型回歸系數(shù)達(dá)0.87,且具有顯著相關(guān)性,說明捕撈壓力的增大是導(dǎo)致黃渤海小黃魚漁獲量變動的主要原因。黃渤海是我國開發(fā)較早且?guī)砭薮蠛Q蠼?jīng)濟(jì)效益的海域,1950年黃渤海的海洋捕撈產(chǎn)量已經(jīng)占據(jù)中國大陸捕撈總產(chǎn)量的1/2以上[17]。隨著我國海洋漁業(yè)的迅猛發(fā)展,黃渤海捕撈漁船的捕撈努力量不斷增加,漁船功率從1962年的14.1×104kW增加到2003年的382.9×104kW[12]。黃渤海小黃魚在1960—1980年代由于過度捕撈導(dǎo)致漁獲量緩慢下降,1972年降至最低,僅為0.41×104t[12]。自1980年代后期以來,漁船捕撈壓力急劇增加,漁獲量迅速上升,造成嚴(yán)重的捕撈過度和資源衰退,導(dǎo)致小黃魚個(gè)體生長速度升高[18],出現(xiàn)了生命周期縮短、性成熟提前、個(gè)體繁殖力上升等生物學(xué)特性的演變[19],以適應(yīng)高強(qiáng)度的捕撈壓力,維持種群的延續(xù)。2010年黃渤海小黃魚的漁獲量達(dá)到最高值,為26.7×104t[12]。本研究表明,如果僅考慮捕撈努力量的影響,根據(jù)Fox模型擬合出的黃渤海小黃魚漁獲量與其實(shí)際漁獲量的變動趨勢存在一定差異(見圖1),這說明小黃魚漁獲量的變動不僅要考慮捕撈壓力的影響,還應(yīng)結(jié)合相應(yīng)時(shí)間段內(nèi)的環(huán)境氣候因子的變化來解釋[8,10,13]。

    3.2 氣候變化對黃渤海小黃魚漁獲量的影響

    本研究發(fā)現(xiàn),加入氣候因子的Fox模型與小黃魚實(shí)際漁獲量擬合度更優(yōu),回歸系數(shù)達(dá)0.90,且置信水平達(dá)99%以上,說明黃渤海小黃魚漁獲量的變動與捕撈壓力和氣候變化均密切相關(guān)。影響黃渤海小黃魚漁獲量的氣候因子主要有:黃海海表溫度、渤海冬季季風(fēng)、黃海夏季季風(fēng)和NPI指數(shù)。

    黃海海表面溫度與小黃魚漁獲量變動呈顯著正相關(guān)。研究表明,小黃魚屬于暖溫性魚類,主要分布在暖水域或其邊緣地帶[1]。在適溫范圍內(nèi),隨著溫度適當(dāng)升高,小黃魚達(dá)到性成熟所需的有效積溫時(shí)間提前,使得繁殖盛期提前[20],有利于促進(jìn)小黃魚的性腺發(fā)育,并增加小黃魚產(chǎn)卵數(shù)量[21]。另外,溫度適當(dāng)升高有利于浮游植物的光合作用,提高海域的初級生產(chǎn)力[22],進(jìn)而為小黃魚提供充足的食物供給,對增加補(bǔ)充群體數(shù)量的有著重要影響。

    渤海冬季季風(fēng)、黃海夏季季風(fēng)與小黃魚漁獲量變動分別呈顯著負(fù)相關(guān)和顯著正相關(guān)。渤海冬季季風(fēng)產(chǎn)生負(fù)面影響可能與渤海季風(fēng)造成海域中營養(yǎng)鹽流失有關(guān),渤海的季風(fēng)會使渤海海域產(chǎn)生風(fēng)生環(huán)流,在該環(huán)流驅(qū)使下,渤海高營養(yǎng)鹽的海水會從南部峽口處流失,而缺乏營養(yǎng)鹽的外海水又從峽口處流入,從而造成渤海初級生產(chǎn)力下降,進(jìn)而對漁獲量產(chǎn)生負(fù)面影響[23]。黃海夏季季風(fēng)產(chǎn)生正面影響可能由于黃海海域夏季盛行西南風(fēng),阻礙了南下的黃海沿岸流[24],減弱了黃海中部冷水團(tuán)與呂泗至長江口水域海水混合導(dǎo)致的溫度變化,為小黃魚提供穩(wěn)定的水溫環(huán)境[25]。王躍中等[13]研究發(fā)現(xiàn),與小黃魚同屬于暖溫種的東海馬面鲀(Navodonspp.),其漁獲量的年間變動也與海表溫度及黃海夏季季風(fēng)呈顯著正相關(guān),與本文的研究結(jié)論相似。

    北太平洋指數(shù)NPI與黃渤海小黃魚漁獲量變動呈顯著負(fù)相關(guān)。NPI是反映阿留申低壓變化和異常的環(huán)流指數(shù)之一。阿留申低壓是北半球主要的半永久大氣活動中心之一,其強(qiáng)度和位置異常不僅導(dǎo)致“下游”美國西海岸和北美地區(qū)的氣候異常,也會影響“上游”西太平洋和東亞地區(qū)的天氣和氣候。當(dāng)NPI降低時(shí),說明黃渤海氣溫偏高,導(dǎo)致黃渤海區(qū)域初級生產(chǎn)力升高,且有利于小黃魚性腺發(fā)育和產(chǎn)卵。這一結(jié)論,也與上文中黃渤海水溫升高使?jié)O獲量增加的結(jié)論相一致。根據(jù)相關(guān)報(bào)道,在1962—2012年間共出現(xiàn)過6次強(qiáng)厄爾尼諾事件[26],在其發(fā)生的年份黃渤海小黃魚的漁獲量均呈現(xiàn)出不同程度的上升趨勢。例如,在1986—1988年發(fā)生強(qiáng)厄爾尼諾事件期間[26],黃渤海小黃魚的平均漁獲量為1.70×104t,均高于之前一年及之后一年的漁獲量(分別為1.26×104和1.29×104t)。進(jìn)一步證明海表面溫度升高有助于小黃魚漁獲量增加。

    本研究根據(jù)黃渤海小黃魚漁獲量變動趨勢選擇Fox模型進(jìn)行擬合[27],經(jīng)典的剩余產(chǎn)量模型假設(shè)漁業(yè)資源種群是處于平衡狀態(tài)的,可是在現(xiàn)實(shí)漁業(yè)中,漁業(yè)資源很少有處于平衡狀態(tài)的情形[28]。近年來,在持續(xù)高強(qiáng)度捕撈狀態(tài)下[29],小黃魚漁業(yè)面臨資源崩潰的風(fēng)險(xiǎn),小黃魚生物學(xué)表型特征的進(jìn)化演變幾乎已達(dá)極限的危險(xiǎn)預(yù)警[19],因此當(dāng)前的剩余產(chǎn)量模型可能無法完全準(zhǔn)確反映近年來小黃魚漁獲量的動態(tài)變化。在今后的研究中需要考慮根據(jù)小黃魚不同漁業(yè)狀態(tài)選擇多種模型進(jìn)行擬合。

    致謝:感謝中國海洋大學(xué)的田永軍教授對本文提出了寶貴的修改意見,特此致謝!

    [1] 李建生, 林龍山, 程家驊. 東海北部秋季小黃魚分布特征及其與底層溫度和鹽度的關(guān)系[J]. 中國水產(chǎn)科學(xué), 2009, 16(3): 348-356. Li Jian-Sheng, Lin Long-Shan, Cheng Jia-Hua. Distribution characteristic of small yellow croaker (LarimichthyspolyactisBleeker) and its relationship with bottom water temperature and salinity in the northern East China Sea in autumn[J]. Journal of Fishery Sciences of China, 2009, 16(3): 348-356.

    [2] Chikuni S. The fish resources of the northwest Pacific: FAO Fish[J]. Tech Pap, 1985, 266: 190.

    [3] Wang S D, Zhan B. Marine fishery resource management in P. R. China[J]. Marine Policy, 1992, 16: 197-209.

    [4] Yu H. Marine fishery management in P.R. China[J]. Marine Policy, 1991, 15: 23-32.

    [5] 林龍山, 程家驊, 任一平, 等. 東海區(qū)小黃魚種群生物學(xué)特性的分析[J]. 中國水產(chǎn)科學(xué), 2004, 11(4): 333-338. Lin Lin-Shan, Cheng Jia-Hua, Ren Yi-Ping, et al. Analysis of population biology of small yellow croakerPseudosciaenapolyactisin the East China Sea region[J]. Journal of Fishery Sciences of China, 2004, 11(4): 333-338.

    [6] FAO. The State of World Fisheries and Aquaculture 2010[R]. Rome: FAO, 2010: 115-117.

    [7] Chang I Z, Jae B L, Suam K, et al. Climatic regime shifts and their impacts on marine ecosystem and fisheries resources in Korean waters[J]. Prog Oceanogr, 2000, 47: 171-190.

    [8] 王躍中, 孫典榮, 林昭進(jìn), 等. 捕撈壓力和氣候因素對黃渤海帶魚漁獲量變化的影響[J]. 中國水產(chǎn)科學(xué), 2012, 19(6): 1043-1050. Wang Yue-Zhong, Sun Dian-Rong, Lin Zhao-Jin, et al. Analysis on responses of hairtail catches to fishing and climate factors in the Yellow Sea and Bohai Sea, China[J]. Journal of Fishery Sciences of China, 2012, 19(6): 1043-1050.

    [9] Tian Y, Uchikawa K, Ueda Y, et al. Comparison of fluctuations in fish communities and trophic structures of ecosystems from three currents around Japan: Synchronies and differences[J]. ICES Journal of Marine Science, 2014, 71(1): 19-34.

    [10] 劉允芬. 氣候變化對我國沿海漁業(yè)生產(chǎn)影響的評價(jià)[J]. 中國農(nóng)業(yè)氣象, 2000, 21(4): 1-5. Liu Yun-Fen. Study of the impact of climatic change on fishery production at coastal areas in China[J]. Chinese Journal of Agrometeorology, 2000, 21(4): 1-5.

    [11] 農(nóng)業(yè)部水產(chǎn)司. 中國漁業(yè)統(tǒng)計(jì)四十年[M]. 北京: 海洋出版社, 1991: 7-185. Department of Fisheries in the Ministry of Agriculture of China. China Fishery Statistical Yearbook in Forty Years[M]. Beijing: China Ocean Press, 1991: 7-185.

    [12] 農(nóng)業(yè)部漁業(yè)局. 中國漁業(yè)統(tǒng)計(jì)年鑒[M]. 北京: 中國農(nóng)業(yè)出版社, 1962-2012: 10-80. Ministry of agriculture fisheries bureau. China Fishery Statistical Yearbook[M]. Beijing: China Agriculture Press, 1962—2012: 10-80.

    [13] 王躍中, 孫典榮, 賈曉平, 等. 捕撈壓力和氣候變化對東海馬面鲀漁獲量的影響[J]. 南方水產(chǎn)科學(xué), 2013, 9(1): 8-15. Wang Yue-Zhong, Sun Dian-Rong, Jia Xiao-Ping, et al. Influence of fishing pressure and climate change on file fishcatches in East China Sea[J]. South China Fisheries Science, 2013, 9(1): 8-15.

    [14] 管秉賢. 中國沿岸的表面海流與風(fēng)的關(guān)系的初步研究[J]. 海洋與湖沼, 1957, 1(1): 95-116. Guan Bing-Xian. A preliminary study on the relations between the surface currents and the winds along China coast[J]. Oceanologia et Limnologia Sinica, 1957, 1(1): 95-116.

    [15] 詹秉義. 漁業(yè)資源評估[M]. 北京: 中國農(nóng)業(yè)出版社, 1995: 176-178. Chan Bing-Yi. Fish Stock Assessment[M]. Beijing: China Agriculture Press, 1995: 176-178.

    [16] 王繼隆, 李繼龍, 楊文波, 等. 利用氣候因子對Fox 模型計(jì)算東??偨?jīng)濟(jì)魚類CPUE的優(yōu)化[J]. 中國水產(chǎn)科學(xué), 2011, 18(1): 136-144. Wang Ji-Long, Li Ji-Long, Yang Wen-Bo, et al. Using climatic factors to optimize Fox surplus yield model for estimating CPUE of total commercial fishes of East China Sea[J]. Journal of Fishery Sciences of China, 2011, 18(1): 136-144.

    [17] 王建文, 張洪濤, 張合成, 等. 我國專屬經(jīng)濟(jì)區(qū)和大陸架勘測專項(xiàng)綜合報(bào)告[R]. 北京: 海洋出版社, 2002: 393-406. Wang Jian-Wen, Zhang Hong-Tao, Zhang He-Cheng, et al. Report of Special Investigations in the Exclusive Economic Zone and Continental Shelf of China[R]. Beijing: China Ocean Press, 2002: 393-406.

    [18] 李忠爐, 金顯仕, 單秀娟, 等. 小黃魚體長-體質(zhì)量關(guān)系和肥滿度的年際變化[J]. 中國水產(chǎn)科學(xué), 2011, 18(3): 602-610. Li Zhong-Lu, Jin Xian-Shi, Shan Xiu-Juan, et al. Inter-annual changes on body weight-length relationship and relative fatness of small yellow croaker (Larimichthyspolyactis) [J]. Journal of Fishery Sciences of China, 2011, 18(3): 602-610.

    [19] 嚴(yán)利平, 劉尊雷, 張輝, 等. 小黃魚生物學(xué)特征與資源數(shù)量的演變[J]. 海洋漁業(yè), 2014, 36(6): 481-488. Yan Li-Ping, Liu Zun-Lei, Zhang Hui, et al. On the evolution of biological characteristics and resources of small yellow croaker[J]. Marine Fisheries, 2014, 36(6): 481-488.

    [20] 單秀娟, 李忠爐, 戴芳群, 等. 黃海中南部小黃魚種群生物學(xué)特征的季節(jié)變化和年際變化[J]. 漁業(yè)科學(xué)進(jìn)展, 2011, 32(6): 7-16. Shan Xiu-Juan, Li Zhong-Lu, Dai Fang-qun, et al. Seasonal and annual variations in biological characteristics of small yellow croakerLarimichthyspolyactisin the central and southern Yellow Sea[J]. Progress in Fishery Sciences, 2011, 32(6): 7-16.

    [21] 金顯仕. 黃海小黃魚生態(tài)和種群動態(tài)的研究[J]. 中國水產(chǎn)科學(xué), 1996, 3(1): 32-46. Jin Xian-Shi. Ecology and population dynamics of small yellow croaker in the yellow sea[J]. Journal of Fishery Sciences of China, 1996, 3(1): 32-46.

    [22] 陳新軍. 漁業(yè)資源與漁場學(xué)[M]. 北京: 海洋出版社, 2004: 116-134. Chen Xin-Jun. Fishery Resources and Fishing Grounds[M]. Beijing: China Ocean Press, 2004: 116-134.

    [23] Lin C L, Su J L, Xu B R, et al. Long-term variations of temperature and salinity of Bohai Sea and their influence on its ecosystem[J]. Prog Oceanogr, 2001, 49: 7-19.

    [24] 陳冠賢, 顧新根, 郜洪仁, 等. 中國海洋漁業(yè)環(huán)境[M]. 杭州: 浙江科學(xué)出版社, 1991: 3-112. Chen Guan-Xian, Gu Xin-Gen, Gao Hong-Ren, et al. Marine Fishery Environment in China[M]. Hangzhou: Zhejiang Science Press, 1991: 3-112.

    [25] 劉尊雷, 袁興偉, 楊林林, 等. 氣候變化對東海北部外海越冬場漁業(yè)群落格局的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2015, 26(3): 901-911. Liu Zun-Lei, Yuan Xing-Wei, Yang Lin-Lin, et al. Effect of climate change on the fisheries community pattern in the overwintering ground of open waters of northern East China Sea[J]. Chinese Journal of Applied Ecology, 2015, 26(3): 901-911.

    [26] 唐曉春, 袁中友. 近60年來厄爾尼諾事件對廣東省旱災(zāi)的影響[J]. 地理研究, 2010, 29(11): 1932-1939. Tang Xiao-Chun, Yuan Zhong-You. The influence of the drought disaster caused by El Ninoevents in Guangdong Province in recent 60 years[J]. Geographical Research, 2010, 29(11): 1932-1939.

    [27] 崔鶴. 連續(xù)形 Fox 產(chǎn)量模型在模擬和實(shí)際漁業(yè)評估中的應(yīng)用[D]. 青島: 中國海洋大學(xué), 2008. Cui He. Investigation of a Continuous Fox-Form Production Model for Fishery Stock Assessment[D]. Qin Dao: Ocean University of China, 2008.

    [28] 王玉. 剩余產(chǎn)量模型的研究[D]. 青島: 中國海洋大學(xué), 2013. Wang Yu. Study of the Surplus Production Models[D]. Qing dao: Ocean University of China, 2013.

    [29] 林龍山, 程家驊. 東海區(qū)小黃魚漁業(yè)生物學(xué)現(xiàn)狀[J]. 中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2004, 34(4): 565-570. Lin Long-Shan, Cheng Jia-Hua. An analysis of the current situation of fishery biology of smallyellow croaker in the east China sea[J]. Periodical of Ocean University of China, 2004, 34(4): 565-570.

    責(zé)任編輯 朱寶象

    Impacts of Fishing Pressure and Climate Change on Catches of Small Yellow Croaker in the Yellow Sea and the Bohai Sea

    LIU Xiao-Xiao, WANG Jing, XU Bin-Duo, XUE Ying, REN Yi-Ping

    (College of Fisheries, Ocean University of China, Qingdao 266003, China)

    The influence of fishing pressure and climate change on the catch of small yellow croaker from Yellow Sea and Bohai Sea was studied according to the data collected from 1962 to 2012. We used the Fox surplus production model to fit the influence. The partial correlation between the variation of small yellow croaker catch and the climatic variables wasdetermined using a multivariate linear regression equation according to the inter-annual variation of multiple climatic factors and climatic indexes. We found that small yellow croaker catch from Yellow Sea and Bohai Sea was relatively stable from 1962 to 1989, but soared from 1990 to 2010. After removing variation trend of catch caused by increasing fishing pressure, the variation of catch was significantly partial correlated with sea surface temperature and summer monsoon in Yellow Sea(P<0.05), while negatively partial correlated with North Pacific Index and winter monsoon in Bohai Sea(P<0.05). The Fox model was optimized by considering the influence of climatic factors. The catch of small yellow croaker fitted by fishing effort and climatic variables was significantly correlated to the actual small yellow croaker catch from Yellow Sea and Bohai Sea(R=0.90,P<0.01), which indicated that variation of small yellow croaker catch was affected by fishing pressure and climate change. According to the future scenario of climate, small yellow croaker catch will experience greater fluctuation in future.

    Yellow Sea and Bohai Sea; small yellow croaker catch; Fox model; fishing pressure; climate change

    青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室鰲山科技創(chuàng)新計(jì)劃項(xiàng)目(2015ASKJ02);高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金項(xiàng)目(20120132130001);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(201262004,201562030)資助 Supported by The Scientific and Technological Innovation Project Financially Supported by Qingdao National Laboratory for Marine Science and Technology(2015ASKJ02);the Specialized Research Fund for the Doctoral Program of Higher Education(20120132130001);the Fundamental Research Funds for the Central Universities(201262004,201562030)

    2016-07-21;

    2016-12-10

    劉笑笑(1992-),女,碩士生。E-mail:476166776@qq.com

    ** 通訊作者:E-mail:xueying@ouc.edu.cn

    S932.4

    A

    1672-5174(2017)08-058-07

    10.16441/j.cnki.hdxb.20160263

    劉笑笑, 王晶, 徐賓鐸, 等. 捕撈壓力和氣候變化對黃渤海小黃魚漁獲量的影響[J]. 中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2017, 47(8): 58-64.

    LIU Xiao-Xiao, WANG Jing, XU Bin-Duo, et al.Impacts of fishing pressure and climate change on catches of small yellow croaker in the Yellow Sea and Bohai Sea[J]. Periodical of Ocean University of China, 2017, 47(8): 58-64.

    猜你喜歡
    黃渤海漁獲量小黃魚
    黃渤海不同水層中浮游植物對灰霾添加的響應(yīng)
    小黃魚
    東南太平洋公海魷魚釣漁業(yè)狀況分析
    世界竹筴魚資源開發(fā)現(xiàn)狀及其建議
    炸小黃魚
    小讀者(2020年4期)2020-06-16 03:34:08
    出海
    休漁補(bǔ)網(wǎng)忙
    鄱陽湖湖口刀鱭幼魚出現(xiàn)的時(shí)間及漁獲量
    黃渤海區(qū)域衛(wèi)星云圖反演海面溫度設(shè)計(jì)初探夏光濱
    小心染色小黃魚
    女人爽到高潮嗷嗷叫在线视频| 99热只有精品国产| 久久久国产精品麻豆| 亚洲黑人精品在线| 一进一出抽搐动态| 国产激情久久老熟女| 成人永久免费在线观看视频| 日本五十路高清| 岛国毛片在线播放| 在线观看免费午夜福利视频| xxx96com| 黄网站色视频无遮挡免费观看| 亚洲精品久久成人aⅴ小说| 黄频高清免费视频| 人人妻人人澡人人看| 黄色怎么调成土黄色| 亚洲欧美精品综合一区二区三区| 日韩欧美一区视频在线观看| 色精品久久人妻99蜜桃| 人妻久久中文字幕网| 欧美乱色亚洲激情| 涩涩av久久男人的天堂| 黑丝袜美女国产一区| 99精品久久久久人妻精品| 国产精品欧美亚洲77777| 免费在线观看视频国产中文字幕亚洲| 黑人巨大精品欧美一区二区蜜桃| 国产精品一区二区在线观看99| 在线观看日韩欧美| 午夜日韩欧美国产| 91字幕亚洲| 老汉色av国产亚洲站长工具| 黄色片一级片一级黄色片| 99精品在免费线老司机午夜| 老熟女久久久| 91av网站免费观看| 国产男女内射视频| 国产亚洲精品久久久久久毛片 | 两性夫妻黄色片| www.熟女人妻精品国产| 日韩中文字幕欧美一区二区| 亚洲精品美女久久av网站| 久久性视频一级片| 中文字幕av电影在线播放| 飞空精品影院首页| 国产一区二区三区在线臀色熟女 | 黄色毛片三级朝国网站| 亚洲精品在线观看二区| 亚洲av电影在线进入| 久久久精品国产亚洲av高清涩受| 日韩免费高清中文字幕av| 精品人妻1区二区| 午夜福利在线免费观看网站| 中出人妻视频一区二区| 亚洲欧美精品综合一区二区三区| 亚洲成av片中文字幕在线观看| 青草久久国产| 丁香六月欧美| 免费黄频网站在线观看国产| 久久人妻熟女aⅴ| 国产精品亚洲av一区麻豆| 欧美日韩国产mv在线观看视频| 日韩制服丝袜自拍偷拍| 国产三级黄色录像| 精品久久蜜臀av无| 曰老女人黄片| 黄色女人牲交| 一区在线观看完整版| 人人妻人人添人人爽欧美一区卜| 欧美 亚洲 国产 日韩一| 国产精品99久久99久久久不卡| 婷婷精品国产亚洲av在线 | 下体分泌物呈黄色| 精品乱码久久久久久99久播| 一区在线观看完整版| 日日摸夜夜添夜夜添小说| 中文字幕制服av| 国产精品.久久久| 久久这里只有精品19| 黄色怎么调成土黄色| 12—13女人毛片做爰片一| 啦啦啦 在线观看视频| 18禁黄网站禁片午夜丰满| 精品高清国产在线一区| 国产精品久久久av美女十八| a在线观看视频网站| xxx96com| 成人av一区二区三区在线看| 国产三级黄色录像| 色综合欧美亚洲国产小说| 女性被躁到高潮视频| 国产精品香港三级国产av潘金莲| 国产精品久久久久久人妻精品电影| 欧美另类亚洲清纯唯美| 少妇 在线观看| 热99久久久久精品小说推荐| 操出白浆在线播放| 99国产精品一区二区三区| 日韩 欧美 亚洲 中文字幕| 69av精品久久久久久| 日韩欧美在线二视频 | 欧美日韩一级在线毛片| 王馨瑶露胸无遮挡在线观看| 国产又色又爽无遮挡免费看| 如日韩欧美国产精品一区二区三区| 黑人欧美特级aaaaaa片| 黄色 视频免费看| 女人久久www免费人成看片| 久99久视频精品免费| 99国产精品99久久久久| 久久狼人影院| 中亚洲国语对白在线视频| 最新美女视频免费是黄的| 妹子高潮喷水视频| 亚洲中文av在线| 国产欧美日韩一区二区精品| 免费一级毛片在线播放高清视频 | 精品国内亚洲2022精品成人 | 国产精品一区二区在线观看99| 久9热在线精品视频| 50天的宝宝边吃奶边哭怎么回事| 老汉色∧v一级毛片| 精品人妻熟女毛片av久久网站| 大香蕉久久成人网| 一区二区三区国产精品乱码| 国产在线观看jvid| 亚洲人成电影观看| 久久人人爽av亚洲精品天堂| 三上悠亚av全集在线观看| 亚洲视频免费观看视频| 欧美色视频一区免费| 精品少妇久久久久久888优播| 不卡一级毛片| 色94色欧美一区二区| 日本wwww免费看| 热re99久久国产66热| 999久久久精品免费观看国产| 免费在线观看完整版高清| 如日韩欧美国产精品一区二区三区| 一级毛片女人18水好多| 婷婷精品国产亚洲av在线 | 欧美激情久久久久久爽电影 | 亚洲午夜精品一区,二区,三区| 国产一区在线观看成人免费| 国产黄色免费在线视频| 亚洲精品久久午夜乱码| 亚洲av成人不卡在线观看播放网| 免费观看人在逋| 欧美成人免费av一区二区三区 | av有码第一页| 国产97色在线日韩免费| 97人妻天天添夜夜摸| 国产免费av片在线观看野外av| 高清毛片免费观看视频网站 | 久久精品国产99精品国产亚洲性色 | 亚洲全国av大片| 91成年电影在线观看| e午夜精品久久久久久久| netflix在线观看网站| 国产男靠女视频免费网站| 午夜亚洲福利在线播放| 欧美久久黑人一区二区| 乱人伦中国视频| 国产精品一区二区在线观看99| 99热国产这里只有精品6| 精品国产国语对白av| 丁香欧美五月| 欧美日韩成人在线一区二区| 色播在线永久视频| 亚洲精品久久成人aⅴ小说| 18在线观看网站| 9色porny在线观看| 精品国产亚洲在线| 欧美激情极品国产一区二区三区| 国产成人免费观看mmmm| 在线观看66精品国产| 黄色视频,在线免费观看| 免费在线观看亚洲国产| 午夜激情av网站| 国产精品二区激情视频| 中出人妻视频一区二区| 黄色 视频免费看| 国产伦人伦偷精品视频| 高潮久久久久久久久久久不卡| 乱人伦中国视频| 最近最新免费中文字幕在线| 美国免费a级毛片| 国产av一区二区精品久久| 别揉我奶头~嗯~啊~动态视频| 好看av亚洲va欧美ⅴa在| 午夜福利乱码中文字幕| 免费在线观看亚洲国产| 久久久久久免费高清国产稀缺| 黄色视频,在线免费观看| 伦理电影免费视频| 日本欧美视频一区| 欧美日韩一级在线毛片| 老司机在亚洲福利影院| 国产男女超爽视频在线观看| 夫妻午夜视频| 精品亚洲成a人片在线观看| 亚洲精品自拍成人| 9191精品国产免费久久| 丰满的人妻完整版| 日韩一卡2卡3卡4卡2021年| 水蜜桃什么品种好| 欧美日韩亚洲综合一区二区三区_| svipshipincom国产片| 法律面前人人平等表现在哪些方面| 国产精品久久久久久精品古装| 日韩欧美免费精品| 久久青草综合色| 国产在线观看jvid| 中文字幕人妻熟女乱码| 黄色成人免费大全| 日本a在线网址| 淫妇啪啪啪对白视频| 亚洲美女黄片视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品国产高清国产av | 午夜福利影视在线免费观看| 亚洲熟妇熟女久久| 成人三级做爰电影| 嫩草影视91久久| 51午夜福利影视在线观看| 久久久久久久久免费视频了| 亚洲第一av免费看| 日本一区二区免费在线视频| 精品少妇一区二区三区视频日本电影| 国内毛片毛片毛片毛片毛片| 99精品在免费线老司机午夜| 亚洲,欧美精品.| 国产97色在线日韩免费| 悠悠久久av| 欧美在线一区亚洲| 韩国av一区二区三区四区| 亚洲第一青青草原| 777米奇影视久久| 国产高清视频在线播放一区| 国产精品 国内视频| 国产精品电影一区二区三区 | 99精品在免费线老司机午夜| 黄片小视频在线播放| 欧美丝袜亚洲另类 | 曰老女人黄片| 黄频高清免费视频| 韩国av一区二区三区四区| 深夜精品福利| av一本久久久久| 极品少妇高潮喷水抽搐| 99热国产这里只有精品6| 一级黄色大片毛片| 亚洲人成77777在线视频| xxx96com| 亚洲国产看品久久| 国产区一区二久久| 美女视频免费永久观看网站| 这个男人来自地球电影免费观看| 欧美人与性动交α欧美软件| 亚洲av成人av| av天堂久久9| 欧美av亚洲av综合av国产av| 在线观看免费午夜福利视频| 婷婷成人精品国产| 高清毛片免费观看视频网站 | 一边摸一边做爽爽视频免费| 高清黄色对白视频在线免费看| 热re99久久国产66热| 日韩欧美三级三区| 亚洲一区高清亚洲精品| 久久人人97超碰香蕉20202| 免费看十八禁软件| 欧美 亚洲 国产 日韩一| 日韩熟女老妇一区二区性免费视频| 大片电影免费在线观看免费| 国产精品影院久久| 欧美精品一区二区免费开放| 一边摸一边做爽爽视频免费| 日本a在线网址| 色播在线永久视频| 视频区欧美日本亚洲| 国内久久婷婷六月综合欲色啪| 国产在线一区二区三区精| 亚洲色图 男人天堂 中文字幕| 国产精品国产高清国产av | 黑人巨大精品欧美一区二区mp4| 成在线人永久免费视频| 真人做人爱边吃奶动态| 国产精品九九99| 黑人欧美特级aaaaaa片| 久久中文看片网| 免费久久久久久久精品成人欧美视频| 欧美精品高潮呻吟av久久| 欧美日韩一级在线毛片| 女人被躁到高潮嗷嗷叫费观| 日韩欧美免费精品| 法律面前人人平等表现在哪些方面| 亚洲欧美激情在线| 搡老岳熟女国产| 他把我摸到了高潮在线观看| 午夜福利在线观看吧| √禁漫天堂资源中文www| 欧美精品啪啪一区二区三区| 午夜福利欧美成人| 99riav亚洲国产免费| 91精品国产国语对白视频| 飞空精品影院首页| 欧美日韩中文字幕国产精品一区二区三区 | 精品国产乱子伦一区二区三区| 搡老乐熟女国产| 激情视频va一区二区三区| 国产精品九九99| 久久久精品国产亚洲av高清涩受| 我的亚洲天堂| 久久婷婷成人综合色麻豆| 国产精品自产拍在线观看55亚洲 | 精品久久久久久,| 欧美精品人与动牲交sv欧美| 国产成人av激情在线播放| 亚洲伊人色综图| 中亚洲国语对白在线视频| 国产1区2区3区精品| 亚洲片人在线观看| 久久久国产成人免费| 国产免费男女视频| 宅男免费午夜| aaaaa片日本免费| 下体分泌物呈黄色| 黄色视频不卡| 欧美另类亚洲清纯唯美| 91字幕亚洲| 又大又爽又粗| 午夜91福利影院| 人人妻人人爽人人添夜夜欢视频| 久久热在线av| 热99久久久久精品小说推荐| 制服人妻中文乱码| 亚洲国产精品sss在线观看 | 国产成人精品无人区| 热re99久久国产66热| 亚洲男人天堂网一区| 国产高清videossex| 国产精品香港三级国产av潘金莲| 午夜久久久在线观看| 亚洲国产看品久久| 欧美一级毛片孕妇| 国产日韩欧美亚洲二区| 国产亚洲欧美98| 黑人巨大精品欧美一区二区蜜桃| 日本撒尿小便嘘嘘汇集6| 伊人久久大香线蕉亚洲五| 日韩有码中文字幕| 国产精品影院久久| 十分钟在线观看高清视频www| 女人高潮潮喷娇喘18禁视频| 国产单亲对白刺激| 一a级毛片在线观看| 国产亚洲欧美精品永久| 久久国产亚洲av麻豆专区| 亚洲,欧美精品.| 中文字幕最新亚洲高清| 99热只有精品国产| 欧美久久黑人一区二区| 久久精品国产清高在天天线| 91老司机精品| 亚洲成av片中文字幕在线观看| 国产精品一区二区免费欧美| 国产精品二区激情视频| 怎么达到女性高潮| 岛国毛片在线播放| 亚洲午夜理论影院| 精品免费久久久久久久清纯 | 久久性视频一级片| 女人被躁到高潮嗷嗷叫费观| 操出白浆在线播放| 极品少妇高潮喷水抽搐| 国产97色在线日韩免费| 精品国产一区二区三区四区第35| 成人18禁在线播放| 深夜精品福利| xxxhd国产人妻xxx| 露出奶头的视频| 午夜91福利影院| 国内毛片毛片毛片毛片毛片| 色老头精品视频在线观看| 啪啪无遮挡十八禁网站| 久久国产精品影院| tocl精华| 三上悠亚av全集在线观看| 黄色片一级片一级黄色片| 99久久人妻综合| 精品无人区乱码1区二区| 乱人伦中国视频| 久久久久久亚洲精品国产蜜桃av| 亚洲国产精品sss在线观看 | 黄色毛片三级朝国网站| 极品教师在线免费播放| 色婷婷久久久亚洲欧美| 国产精品 国内视频| 亚洲欧美日韩另类电影网站| 男人的好看免费观看在线视频 | 天堂√8在线中文| 欧美最黄视频在线播放免费 | 亚洲av欧美aⅴ国产| 黄色成人免费大全| 欧美成狂野欧美在线观看| 国产精品久久久久久精品古装| 久9热在线精品视频| 亚洲专区中文字幕在线| 成人国语在线视频| av在线播放免费不卡| 成熟少妇高潮喷水视频| 他把我摸到了高潮在线观看| 日韩欧美在线二视频 | 18禁裸乳无遮挡动漫免费视频| 在线观看66精品国产| 亚洲国产精品合色在线| 久久精品国产清高在天天线| 精品欧美一区二区三区在线| 在线观看www视频免费| 国产精品av久久久久免费| 久久久久久亚洲精品国产蜜桃av| 国产乱人伦免费视频| 亚洲av美国av| 首页视频小说图片口味搜索| ponron亚洲| 日韩大码丰满熟妇| 老熟妇乱子伦视频在线观看| 少妇裸体淫交视频免费看高清 | 亚洲中文av在线| 精品第一国产精品| 99re在线观看精品视频| 99热网站在线观看| 亚洲中文字幕日韩| av网站免费在线观看视频| 精品久久久精品久久久| 国产三级黄色录像| 亚洲在线自拍视频| 两个人免费观看高清视频| 久久久国产一区二区| 亚洲国产精品sss在线观看 | 女警被强在线播放| 丝袜在线中文字幕| 免费人成视频x8x8入口观看| av不卡在线播放| 一级,二级,三级黄色视频| 精品国产超薄肉色丝袜足j| 女警被强在线播放| 欧美在线黄色| 国产在线精品亚洲第一网站| 五月开心婷婷网| 老司机深夜福利视频在线观看| 日韩人妻精品一区2区三区| 中文亚洲av片在线观看爽 | 国产一区在线观看成人免费| 91大片在线观看| 大香蕉久久网| 又黄又粗又硬又大视频| 午夜福利在线观看吧| 黄色 视频免费看| av片东京热男人的天堂| 国产av一区二区精品久久| 精品少妇久久久久久888优播| 亚洲人成77777在线视频| 久久青草综合色| 国产亚洲欧美在线一区二区| 中文字幕精品免费在线观看视频| 亚洲精华国产精华精| 亚洲精品中文字幕在线视频| 十分钟在线观看高清视频www| 一区二区日韩欧美中文字幕| 99精品在免费线老司机午夜| 天堂中文最新版在线下载| 精品福利观看| av福利片在线| 亚洲国产精品sss在线观看 | 精品久久久精品久久久| 久久精品91无色码中文字幕| 999久久久国产精品视频| 十分钟在线观看高清视频www| 久久久精品免费免费高清| 免费观看精品视频网站| 男女之事视频高清在线观看| 一边摸一边做爽爽视频免费| 精品一区二区三区四区五区乱码| 日韩欧美一区视频在线观看| 国产精品电影一区二区三区 | 丰满饥渴人妻一区二区三| 国产欧美日韩综合在线一区二区| 日本欧美视频一区| 50天的宝宝边吃奶边哭怎么回事| 精品国产乱码久久久久久男人| 午夜激情av网站| 国产熟女午夜一区二区三区| 丁香欧美五月| 50天的宝宝边吃奶边哭怎么回事| 少妇裸体淫交视频免费看高清 | 婷婷精品国产亚洲av在线 | av有码第一页| 欧美不卡视频在线免费观看 | 成人精品一区二区免费| 男女下面插进去视频免费观看| 丰满饥渴人妻一区二区三| aaaaa片日本免费| 黄色女人牲交| 在线国产一区二区在线| 一级片免费观看大全| 婷婷成人精品国产| 久久久久国产一级毛片高清牌| a级片在线免费高清观看视频| 美女福利国产在线| 一级片免费观看大全| 三上悠亚av全集在线观看| 大片电影免费在线观看免费| 在线天堂中文资源库| 色婷婷久久久亚洲欧美| 69精品国产乱码久久久| 在线观看一区二区三区激情| 精品国产乱子伦一区二区三区| 欧美日韩乱码在线| 国产无遮挡羞羞视频在线观看| 国产精品亚洲av一区麻豆| 黄片小视频在线播放| 欧美乱码精品一区二区三区| 黄片大片在线免费观看| 极品少妇高潮喷水抽搐| 成人18禁高潮啪啪吃奶动态图| 久久久国产精品麻豆| 亚洲精品一卡2卡三卡4卡5卡| 一进一出抽搐动态| 午夜福利一区二区在线看| 一级片'在线观看视频| 精品人妻熟女毛片av久久网站| 不卡av一区二区三区| 男女下面插进去视频免费观看| 一级a爱视频在线免费观看| a在线观看视频网站| 女同久久另类99精品国产91| 桃红色精品国产亚洲av| 一级a爱视频在线免费观看| 男女下面插进去视频免费观看| 一本综合久久免费| 天天操日日干夜夜撸| 97人妻天天添夜夜摸| 久久天躁狠狠躁夜夜2o2o| 亚洲欧美激情在线| 男人操女人黄网站| 亚洲专区中文字幕在线| 色婷婷av一区二区三区视频| 日本欧美视频一区| 法律面前人人平等表现在哪些方面| 亚洲精品久久成人aⅴ小说| 丰满人妻熟妇乱又伦精品不卡| 午夜日韩欧美国产| 黄色女人牲交| 成人18禁在线播放| 欧美激情 高清一区二区三区| 精品一区二区三区av网在线观看| 日韩三级视频一区二区三区| 99久久99久久久精品蜜桃| 国产欧美日韩精品亚洲av| 久久国产亚洲av麻豆专区| 叶爱在线成人免费视频播放| 欧美不卡视频在线免费观看 | 午夜影院日韩av| 成熟少妇高潮喷水视频| 国产主播在线观看一区二区| 中文字幕最新亚洲高清| 91精品三级在线观看| 美女福利国产在线| 亚洲国产精品一区二区三区在线| 一边摸一边做爽爽视频免费| 国产极品粉嫩免费观看在线| 精品久久久久久,| 国产精品 欧美亚洲| 99精品欧美一区二区三区四区| 啦啦啦在线免费观看视频4| 久久国产精品大桥未久av| 无限看片的www在线观看| 久久精品国产亚洲av高清一级| 一本大道久久a久久精品| 亚洲av日韩在线播放| 黄色成人免费大全| 久久精品亚洲av国产电影网| 国产亚洲精品久久久久久毛片 | 久久婷婷成人综合色麻豆| 国产成人一区二区三区免费视频网站| 国产野战对白在线观看| 可以免费在线观看a视频的电影网站| 在线观看免费视频网站a站| 久久ye,这里只有精品| 十八禁网站免费在线| 中文字幕高清在线视频| 90打野战视频偷拍视频| 免费在线观看视频国产中文字幕亚洲| а√天堂www在线а√下载 | 免费人成视频x8x8入口观看| 又紧又爽又黄一区二区| 成人黄色视频免费在线看| 一进一出抽搐gif免费好疼 | 亚洲成国产人片在线观看| 日韩 欧美 亚洲 中文字幕| 国内久久婷婷六月综合欲色啪| 一本一本久久a久久精品综合妖精| 日日夜夜操网爽| 香蕉丝袜av| 亚洲午夜精品一区,二区,三区| 久久久久国内视频| 999精品在线视频| av片东京热男人的天堂| 日韩 欧美 亚洲 中文字幕| 国产精品一区二区在线不卡| 久久久久精品人妻al黑| 99久久99久久久精品蜜桃| 亚洲七黄色美女视频|