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

    中國遠(yuǎn)洋作業(yè)漁場(chǎng)海表溫度異常年際變動(dòng)分析

    2016-07-06 09:42:20紀(jì)世建周為峰靳少非崔雪森農(nóng)業(yè)部東海與遠(yuǎn)洋漁業(yè)資源開發(fā)利用重點(diǎn)實(shí)驗(yàn)室上海200090上海海洋大學(xué)海洋科學(xué)學(xué)院上海2006中國科學(xué)院大氣物理研究所東亞區(qū)域氣候環(huán)境重點(diǎn)實(shí)驗(yàn)室北京00029
    海洋科學(xué) 2016年1期
    關(guān)鍵詞:海表海溫漁場(chǎng)

    紀(jì)世建, 周為峰, 樊 偉, 靳少非, 崔雪森(. 農(nóng)業(yè)部東海與遠(yuǎn)洋漁業(yè)資源開發(fā)利用重點(diǎn)實(shí)驗(yàn)室, 上海 200090; 2. 上海海洋大學(xué) 海洋科學(xué)學(xué)院, 上海2006; . 中國科學(xué)院 大氣物理研究所, 東亞區(qū)域氣候環(huán)境重點(diǎn)實(shí)驗(yàn)室, 北京 00029)

    ?

    中國遠(yuǎn)洋作業(yè)漁場(chǎng)海表溫度異常年際變動(dòng)分析

    紀(jì)世建1, 2, 周為峰1, 樊 偉1, 靳少非3, 崔雪森1
    (1. 農(nóng)業(yè)部東海與遠(yuǎn)洋漁業(yè)資源開發(fā)利用重點(diǎn)實(shí)驗(yàn)室, 上海 200090; 2. 上海海洋大學(xué) 海洋科學(xué)學(xué)院, 上海201306; 3. 中國科學(xué)院 大氣物理研究所, 東亞區(qū)域氣候環(huán)境重點(diǎn)實(shí)驗(yàn)室, 北京 100029)

    基于美國國家海洋大氣局(NOAA)氣候預(yù)測(cè)中心月平均SST資料, 利用時(shí)間序列的統(tǒng)計(jì)學(xué)特征分析了中國7個(gè)主要遠(yuǎn)洋作業(yè)漁場(chǎng)1982~2011年海表溫度異常(SSTA)年際變動(dòng), 用功率譜方法計(jì)算時(shí)間序列的顯著變動(dòng)周期, 并用相關(guān)分析探討了去除趨勢(shì)項(xiàng)后的SSTA與南方濤動(dòng)指數(shù)(Southern Oscillation Index, SOI)的相關(guān)性。結(jié)果表明, 1982~2011年中國主要遠(yuǎn)洋作業(yè)漁場(chǎng)SSTA為–0.3~0.3℃, 波動(dòng)周期約為3~4 a, 平均SSTA總體呈現(xiàn)上升趨勢(shì), 與SOI存在顯著的相關(guān)性(r=–0.509), 說明漁場(chǎng)的SSTA 與ENSO事件有著密切聯(lián)系; 從漁場(chǎng)SSTA升降趨勢(shì)來看, 除了東太平洋和東南太平洋SSTA出現(xiàn)下降趨勢(shì), 其余漁場(chǎng)的SSTA均有一定程度的上升, 其中西北太平洋SSTA上升最為顯著; 從漁場(chǎng)SSTA的變化周期來看, 東太平洋和西南大西洋的SSTA變化周期為3~4 a, 東南太平洋為4 a, 西太平洋為5 a,其余的短期周期性較不明顯, 約為10 a; 與SOI時(shí)間序列進(jìn)行相關(guān)分析得到, 東太平洋、中大西洋以及西南大西洋均與SOI存在顯著的相關(guān)性, 相關(guān)系數(shù)分別為–0.895、0.471和–0.598, 其余漁場(chǎng)與SOI無顯著相關(guān)。通過各漁場(chǎng)間的對(duì)比得到以下特征: 赤道附近海域東太平洋SSTA變化往往與西太平洋和中大西洋反相, 而與印度洋同相; 中緯度海域的3個(gè)漁場(chǎng)中, 南半球中緯度漁場(chǎng)溫度變化要比北半球中緯度漁場(chǎng)??; 東、西印度洋SSTA存在顯著相關(guān)性, 印度洋內(nèi)部SSTA正負(fù)變化情況基本一致。

    遠(yuǎn)洋漁場(chǎng); 海表溫度; 海表溫度異常

    [Foundation: National Key Technology Support Program (2013BAD13B06);Scientific Research Program Supported by Shanghai Science and Technology Commission (15DZ1202201)]

    近半個(gè)世紀(jì)以來, 全球氣溫有明顯上升。作為世界最大的熱能儲(chǔ)存庫, 海洋在全球氣候變化過程中占據(jù)著核心地位, 全球溫室氣體釋放的能量中約93%被海洋所吸收[1]。因此, 全球氣溫上升勢(shì)必會(huì)造成海溫變動(dòng)。聯(lián)合國政府間氣候變化專門委員會(huì)(IPCC)第5次評(píng)估報(bào)告指出, 1970年以來海洋在持續(xù)升溫, 75 m以上的上層海水溫度每10年升溫幅度超過0.11℃, 1950 ~2009年間印度洋、大西洋和太平洋的平均海表溫度(SST)分別上升了0.65℃、0.41℃和0.31℃[2]。

    在海洋漁業(yè)領(lǐng)域, 海表溫度是研究海洋漁場(chǎng)環(huán)境的重要因子, 海表溫度異常(sea surface temperature anomaly, SSTA), 又稱海表溫度距平, 常用于捕捉海表溫度在年際間的異常變化特征。海表溫度的異常變動(dòng)如ENSO事件會(huì)影響海中魚類的繁殖、生長、發(fā)育及洄游等行為, 進(jìn)而引起漁場(chǎng)分布以及資源量豐度等發(fā)生一系列變化[3-4]。曹曉怡等[5]在對(duì)印度洋大眼金槍魚和黃鰭金槍魚延繩釣漁場(chǎng)的重心變化分析時(shí), 發(fā)現(xiàn)漁場(chǎng)重心的年間變化受到了ENSO事件的影響, 在發(fā)生厄爾尼諾的當(dāng)年或次年,年漁獲量分布重心向西、南偏移; 在發(fā)生拉尼娜的年份漁獲量分布重心會(huì)向東、北偏移。汪金濤和陳新軍[6]通過研究1990~2010年間ENSO事件對(duì)中西太平洋鰹魚漁場(chǎng)重心變化的影響, 發(fā)現(xiàn)厄爾尼諾現(xiàn)象發(fā)生時(shí), 鰹魚漁獲量重心明顯東移; 在拉尼娜現(xiàn)象發(fā)生時(shí), 鰹魚漁獲量重心有整體西移趨勢(shì)。拉尼娜事件會(huì)導(dǎo)致柔魚補(bǔ)充量下降, 而厄爾尼諾事件會(huì)導(dǎo)致補(bǔ)充量增加, 并且ENSO事件還會(huì)影響柔魚索餌場(chǎng)分布。

    中國遠(yuǎn)洋作業(yè)漁場(chǎng)按海域位置可以分為西北太平洋漁場(chǎng)、東太平洋漁場(chǎng)、東南太平洋漁場(chǎng)、西太平洋漁場(chǎng)、印度洋漁場(chǎng)、中大西洋漁場(chǎng)和西南大西洋漁場(chǎng)[7]。研究中國遠(yuǎn)洋作業(yè)漁場(chǎng)的海溫異常變動(dòng),有助于把握遠(yuǎn)洋漁場(chǎng)時(shí)空分布以及漁業(yè)資源量的變化, 為漁業(yè)資源管理和開發(fā)提供環(huán)境參考, 實(shí)現(xiàn)海洋漁業(yè)資源的可持續(xù)利用。

    1 數(shù)據(jù)與方法

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

    SST數(shù)據(jù)來自美國國家海洋大氣局(NOAA)氣候預(yù)報(bào)中心提供的海溫最優(yōu)插值再分析數(shù)據(jù)(optimum interpolation SST, OISST), 該數(shù)據(jù)經(jīng)過最優(yōu)插值處理并以月為單位, 其空間分辨率為1°×1°。文章截選了1982年1月~2011年12月共30 a的SST數(shù)據(jù)集。為了探究海表溫度異常與ENSO事件的相關(guān)性, 文章還采用了美國氣候中心網(wǎng)站(http: //www.bom.gov. au/climate/current/soi2.shtml)提供的南方濤動(dòng)指數(shù)(Southern Oscillation Index, SOI)1982 ~2011年的時(shí)間序列。SOI是南太平洋塔希提島(Tahiti)與達(dá)爾文港(Darwin)兩地的標(biāo)準(zhǔn)化海平面氣壓差, 它與ENSO關(guān)系密切, 能夠反映ENSO事件的發(fā)展?fàn)顩r及劇烈程度。圖1為中國遠(yuǎn)洋作業(yè)漁場(chǎng)的分布。

    圖1 中國遠(yuǎn)洋作業(yè)漁場(chǎng)分布圖Fig.1 Distribution map of Chinese pelagic operation fishing ground

    1.2 研究方法

    首先, 根據(jù)漁場(chǎng)的區(qū)域范圍對(duì)SST數(shù)據(jù)集進(jìn)行截取, 計(jì)算出所有漁場(chǎng)區(qū)域每年的平均SST, 然后,利用公式(1)計(jì)算海表溫度異常, 式中SSTAi為第i年的海表溫度異常值, SSTi為第i年的海表溫度值,為30 a的平均海表溫度值。

    然后, 計(jì)算所有漁場(chǎng)的SSTA時(shí)間序列的統(tǒng)計(jì)學(xué)特征, 包括最大值、最小值和標(biāo)準(zhǔn)差, 可以得到各漁場(chǎng)SSTA在30 a間的變化范圍和分析各漁場(chǎng)SSTA的長期的升降幅度。最后, 將所有漁場(chǎng)的SSTA時(shí)間序列在EXCEL中用帶直線的散點(diǎn)圖顯示, 可以得到中國遠(yuǎn)洋作業(yè)漁場(chǎng)年平均海表溫度異常時(shí)間序列。用同樣的方法可以得到SOI的時(shí)間序列以及每個(gè)漁場(chǎng)各自年平均海表溫度異常時(shí)間序列,通過這些時(shí)間序列可以直觀的看出其年際間變化趨勢(shì)。

    作者采用matlab中的detrend函數(shù)將時(shí)間序列中的線性趨勢(shì)剔除, 得到平穩(wěn)時(shí)間序列。其原理就是,將待去除趨勢(shì)項(xiàng)的時(shí)間序列用最小二乘法擬合出一個(gè)線性回歸方程, 然后用原來的時(shí)間序列值減去這條回歸線上對(duì)應(yīng)的值, 得到去除趨勢(shì)項(xiàng)后的時(shí)間序列。時(shí)間序列趨勢(shì)項(xiàng)的存在會(huì)導(dǎo)致相關(guān)性分析和功率譜密度分析的結(jié)果不可靠, 去趨勢(shì)項(xiàng)分析可以檢測(cè)非平穩(wěn)時(shí)間序列中長周期波動(dòng)所造成的偽相關(guān)現(xiàn)象[8]。進(jìn)一步利用經(jīng)典功率譜分析方法計(jì)算各漁場(chǎng)SSTA時(shí)間序列的顯著周期, 并用相關(guān)分析方法探討去除趨勢(shì)項(xiàng)后的SSTA和SOI之間的關(guān)系, 分析各漁區(qū)對(duì)ENSO事件的響應(yīng)。

    2 結(jié)果與分析

    2.1 遠(yuǎn)洋作業(yè)漁場(chǎng)SSTA的總體變化分析

    全球氣候變暖引起海洋溫度發(fā)生相應(yīng)變化, 為了呈現(xiàn)出中國遠(yuǎn)洋作業(yè)漁場(chǎng)海表溫度異常在1982~2011年間的總體變化情況, 本文將所有漁場(chǎng)作為一個(gè)整體, 計(jì)算其年平均海表溫度異常, 得到了中國遠(yuǎn)洋作業(yè)漁場(chǎng)的年平均SSTA時(shí)間序列(圖2)。

    圖2 中國遠(yuǎn)洋作業(yè)漁場(chǎng)年平均海表溫度異常SSTA及SOIFig. 2 Annual mean SSTA and SOI of Chinese pelagic fishing grounds

    從總體變化來看, 漁場(chǎng)平均海表溫度異常呈現(xiàn)上升趨勢(shì)。30 a來SSTA在–0.3~0.3℃范圍內(nèi)上下波動(dòng), 波動(dòng)周期約為3~4 a, 最大的負(fù)異常出現(xiàn)在1985年, 最大正異常出現(xiàn)在1998年。1998年以前主要表現(xiàn)為負(fù)異常, 1998年以后主要表現(xiàn)為正異常, 并且在2001~2006年間有一個(gè)持續(xù)穩(wěn)定的正異常特征。與同期SOI時(shí)間序列(圖2)對(duì)比, 可以發(fā)現(xiàn), 在年際尺度上兩者波動(dòng)周期一致, 長期的升降趨勢(shì)基本相反, 對(duì)二者去除趨勢(shì)項(xiàng)后進(jìn)行相關(guān)分析, 發(fā)現(xiàn)二者在0.01的顯著性水平上呈負(fù)相關(guān), 相關(guān)系數(shù)為–0.509, 說明漁場(chǎng)的海表溫度異常與ENSO事件有著密切的聯(lián)系。

    2.2 不同海域漁場(chǎng)SSTA的變化分析

    為了便于分析, 將印度洋漁場(chǎng)分為東印度洋和西印度洋, 各海域漁場(chǎng)的SSTA變化范圍、標(biāo)準(zhǔn)差(變化幅度)、波動(dòng)周期及趨勢(shì)斜率見表1。

    表1 各海域漁場(chǎng)的SSTA變化范圍、標(biāo)準(zhǔn)差、波動(dòng)周期、趨勢(shì)及斜率Tab.1 Variation range, standard deviation, variation period, trend, and slope of pelagic fishing grounds’ SSTA

    位于赤道附近海域的漁場(chǎng)包括東太平洋、西太平洋以及中大西洋3個(gè)漁場(chǎng), 它們的SSTA的時(shí)序變化見圖3。東太平洋SSTA變化非常顯著, 標(biāo)準(zhǔn)差為0.4, 其值在–0.58~0.8℃上下波動(dòng), 1997年出現(xiàn)最大值, 1999年出現(xiàn)最小值, 30 a的總體趨勢(shì)表現(xiàn)為下降,功率譜分析得到的周期約為3~4 a, 與SOI的變化周期一致, 其去除趨勢(shì)項(xiàng)后與SOI時(shí)間序列在0.01的顯著性水平上呈負(fù)相關(guān), 相關(guān)系數(shù)達(dá)–0.895, 具有極高的相關(guān)性, 可見ENSO是東太平洋SSTA年際變化的最強(qiáng)信號(hào); 西太平洋SSTA變化稍顯平緩, 標(biāo)準(zhǔn)差為0.2, 變化范圍在–0.33~0.34℃, 最小值出現(xiàn)在1993年, 最大值出現(xiàn)在2003年, 并且在2001~2007年間保持了較長的正異常狀態(tài), 30 a總體處于上升趨勢(shì), 波動(dòng)周期為5 a, 與SOI時(shí)間序列無顯著的相關(guān)性; 中大西洋年際變化振幅較大, 標(biāo)準(zhǔn)差為0.22, 變化范圍在–0.6~0.52℃, 分別出現(xiàn)在1992年和2010年, 其周期性變化較不明顯, 總體有上升趨勢(shì), 主要的波動(dòng)周期為10 a, 與SOI在0.01的顯著性水平上呈正相關(guān), 相關(guān)系數(shù)為0.471, 說明中大西洋與ENSO事件也存在一定的聯(lián)系。

    圖3 赤道海域漁場(chǎng)的SSTA時(shí)間序列Fig. 3 Time serials of SSTA in pelagic fishing grounds near the equator

    對(duì)比發(fā)現(xiàn), 東太平洋與西太平洋及中大西洋變化基本反相, 尤其在1998年以前。也就是說, 在東太平洋處于極大值的年份中, 西太平洋和中大西洋往往表現(xiàn)為極小值, 反之亦然。這種差異可能與太平洋所發(fā)生的厄爾尼諾和拉尼娜2種現(xiàn)象交替出現(xiàn)有關(guān)。正常情況下, 東太平洋的大量暖水會(huì)被信風(fēng)吹到西太平洋地區(qū), 只能靠海面以下的冷水進(jìn)行補(bǔ)充,東太平洋SST一般要比西太平洋偏低。厄爾尼諾現(xiàn)象發(fā)生時(shí), 東太平洋冷水上翻現(xiàn)象消失, 表層暖水向東回流, 導(dǎo)致太平洋東部SST異常增溫。而在西太平洋出現(xiàn)異常下沉氣流, 海表溫度降低。這種大范圍的對(duì)流異常會(huì)引起熱帶海洋的大氣環(huán)流的變化[9-10]。與此同時(shí)在中大西洋東風(fēng)也會(huì)得到加強(qiáng), 中大西洋溫躍層變淺, 使得中大西洋海表溫度變冷。拉尼娜現(xiàn)象則使東太平洋深層海水上翻現(xiàn)象比正常情況更加劇烈,東部氣流下沉而西部氣流加劇上升, 加劇太平洋東部冷水發(fā)展, SST、大氣環(huán)流以及短波輻射的變化與厄爾尼諾相反。一般拉尼娜現(xiàn)象會(huì)隨著厄爾尼諾現(xiàn)象而來, 出現(xiàn)厄爾尼諾現(xiàn)象的第二年, 往往都會(huì)出現(xiàn)拉尼娜現(xiàn)象, 而且有時(shí)拉尼娜現(xiàn)象會(huì)持續(xù)2~3年。1984年底~1985年、1988~1989年、1995 ~1996年、1998年底~2000年初、2007年底~2008年、2010年底~2011年都發(fā)生了拉尼娜現(xiàn)象, 令太平洋東部、東南部的海水溫度比正常低了1~2℃。東太平洋SSTA 與SOI具有顯著的相關(guān)性, 這正是由于東太平洋是ENSO事件的產(chǎn)生區(qū)域, 每當(dāng)ENSO暖相位(厄爾尼諾現(xiàn)象)和ENSO冷相位(拉尼娜現(xiàn)象)發(fā)生時(shí), 都會(huì)導(dǎo)致太平洋東部SST發(fā)生顯著變動(dòng), 這也解釋了太平洋東部SSTA變化幅度大這一特征。

    位于中緯度海域的漁場(chǎng)包括西北太平洋、東南太平洋和西南大西洋3個(gè)漁場(chǎng), 它們的時(shí)間序列變化見圖4。西北太平洋SSTA變化較為劇烈, 標(biāo)準(zhǔn)差為0.39, 變化范圍在–0.65~0.7℃, 1986年出現(xiàn)最小值,1999年出現(xiàn)最大值, 并且在1998年后基本表現(xiàn)為正異常, 30年的總體變化具有明顯的上升趨勢(shì), 其主要的變化周期為10年, 與SOI相關(guān)性不顯著; 東南太平洋SSTA相較于其他漁場(chǎng)來說, 其變化最為平緩,標(biāo)準(zhǔn)差為0.12, 變化范圍在–0.17~0.22℃, 最小值出現(xiàn)在1986年, 最大值出現(xiàn)在1997年, 其SSTA變化趨勢(shì)線與零值線幾乎重合, 但總體上仍存在輕微的下降趨勢(shì), 變化周期為4 a, 與SOI相關(guān)性不顯著;西南大西洋變化在–0.43~0.25℃, 標(biāo)準(zhǔn)差為0.15, 在1988年有一個(gè)明顯的谷值, 在2005年出現(xiàn)最大值,總體上也有緩慢的上升趨勢(shì), 變化周期3~4 a, 與SOI 在0.01的顯著性水平上呈負(fù)相關(guān), 相關(guān)系數(shù)為–0.598,說明西南大西洋與ENSO事件存在一定的聯(lián)系。

    圖4 中緯度海域漁場(chǎng)的SSTA時(shí)間序列Fig. 4 Time serials of SSTA in mid-attitude pelagic fishing grounds

    中緯度海域的3個(gè)漁場(chǎng)中, 西北太平洋漁場(chǎng)SSTA變化最為顯著, 而東南太平洋和西南大西洋這2個(gè)漁場(chǎng)SSTA變化較為平緩, 這可能是因?yàn)槲鞅碧窖鬂O場(chǎng)位于北半球的中緯度海域, 而后2個(gè)漁場(chǎng)位于南半球中緯度海域, 南半球的中緯度海域面積較大, 海水的熱容量也相對(duì)較大, 其吸收相同的熱量所引起的海表溫度變化就必然比北半球要小。

    印度洋漁場(chǎng)在赤道和中緯度均有分布, 將印度洋漁場(chǎng)按經(jīng)度劃分為東、西印度洋, 分析發(fā)現(xiàn)東印度洋和西印度洋變化趨勢(shì)相似, 在30年間都有一定程度的升溫。其中東印度洋SSTA的變化范圍在–0.31~0.48℃, 標(biāo)準(zhǔn)差為0.18, 最小值出現(xiàn)在1986年,最大值出現(xiàn)在1998年, 波動(dòng)周期為10年, 與SOI時(shí)間序列無顯著相關(guān); 西印度洋SSTA的變化范圍在–0.39~0.26℃, 標(biāo)準(zhǔn)差為0.16, 最小值出現(xiàn)在1984年,最大值出現(xiàn)在2010年(圖5), 波動(dòng)周期為10年, 與東印度洋同期一致, 與SOI時(shí)間序列亦無顯著相關(guān)。對(duì)東印度洋和西印度洋二者進(jìn)行相關(guān)分析, 得到二者在0.01的顯著性水平上呈正相關(guān), 相關(guān)系數(shù)為0.483, 說明東西印度洋的SSTA變化在時(shí)間尺度上具有相關(guān)性。

    對(duì)比發(fā)現(xiàn), 印度洋的SSTA變化與東太平洋變化同相。ENSO通過大氣對(duì)印度洋海溫的影響主要有熱力學(xué)和動(dòng)力學(xué)兩方面因素。(1)KLEIN等[11]認(rèn)為熱帶太平洋和印度洋之間會(huì)產(chǎn)生異常的Walker環(huán)流, 其會(huì)在印度尼西亞和西太平洋地區(qū)形成顯著的海平面氣壓場(chǎng)異常, 進(jìn)而通過海氣相互影響印度洋(尤其是熱帶印度洋)海表云量和風(fēng)場(chǎng)等發(fā)生異常變化, 導(dǎo)致表面凈熱通量變化從而引起海溫的異常。(2)異常風(fēng)場(chǎng)會(huì)使東赤道印度洋沿岸上升流發(fā)生調(diào)制, 從而引起海面溫度異常和溫躍層發(fā)生起伏現(xiàn)象[12]。劉巖松等[13]構(gòu)建了一個(gè)同時(shí)考慮動(dòng)力和熱力的印-太海氣耦合模式, 其實(shí)驗(yàn)結(jié)果表明, 熱帶太平洋暖異常中心更加接近觀測(cè)值, 熱帶印度洋出現(xiàn)了海盆尺度的海溫正異常, 通過模擬的方式說明了ENSO事件會(huì)通過海氣界面的動(dòng)力作用和熱力交換導(dǎo)致印度洋海盆尺度增暖。

    3 討論

    統(tǒng)計(jì)分析結(jié)果表明, 在1982~2011年的SSTA變化中, 中國遠(yuǎn)洋作業(yè)漁場(chǎng)總體SSTA在–0.3~0.3℃范圍內(nèi)上下波動(dòng), 波動(dòng)周期約為3~4 a。與同期SOI時(shí)間序列對(duì)比, 可以發(fā)現(xiàn), 在年際尺度上兩者波動(dòng)周期一致, 長期的升降趨勢(shì)基本相反, 對(duì)二者去除趨勢(shì)項(xiàng)后進(jìn)行相關(guān)分析, 發(fā)現(xiàn)二者在0.01的顯著性水平上呈負(fù)相關(guān), 相關(guān)系數(shù)為–0.509。從各個(gè)漁場(chǎng)SSTA升降趨勢(shì)來看, 除了東太平洋和東南太平洋出現(xiàn)下降趨勢(shì), 其余漁場(chǎng)的SSTA都有一定程度的上升, 其中西北太平洋SSTA上升最為顯著。BELKIN[14]對(duì)1957~2006年全球63個(gè)大海洋生態(tài)系的海表溫度及其異常進(jìn)行長時(shí)序的統(tǒng)計(jì)和分析, 發(fā)現(xiàn)全球的SSTA具有較強(qiáng)的局部變化特征, 所有大海洋生態(tài)系的SSTA的年際變化表現(xiàn)為除了位于美洲西海岸的2個(gè)大海洋生態(tài)系具有海溫下降的趨勢(shì)外, 其余均有一定程度的上升, 這與本文的結(jié)果一致。東太平洋和西南大西洋的SSTA變化周期為3~4 a, 東南太平洋為4 a,西太平洋為5 a, 其余的短期周期性較不明顯, 約為10 a。與SOI時(shí)間序列進(jìn)行相關(guān)分析得到, 東太平洋、中大西洋以及西南大西洋均與SOI存在顯著的相關(guān)性, 相關(guān)系數(shù)分別為–0.895、0.471和–0.598。東太平洋處于正異常的年份中, 西太平洋往往表現(xiàn)為負(fù)異常, 印度洋與東太平洋同相, 中大西洋與西太平洋同相, 表明在赤道附近海域存在著一定的聯(lián)系。事實(shí)上, LAU和NATH[15]提出了“熱帶大氣橋”理論來解釋赤道附近海域之間的聯(lián)系, 他們認(rèn)為與ENSO有關(guān)的大氣環(huán)流異常會(huì)導(dǎo)致蒸發(fā)和云量的異常, 進(jìn)而影響其他海盆的凈熱通量變化并形成局部的海表溫度異常。在分析中緯度漁場(chǎng)時(shí)發(fā)現(xiàn), 南半球中緯度漁場(chǎng)溫度變化要比北半球中緯度漁場(chǎng)更加穩(wěn)定。這一結(jié)果與前人的研究結(jié)果相符, 如江志紅等[16]在分析南北半球的月平均溫度序列時(shí)發(fā)現(xiàn)在全球變暖的背景下, 南半球溫度變化要比北半球更穩(wěn)定。東、西印度洋具有相同的周期, 并且二者存在顯者的相關(guān)性,相關(guān)系數(shù)為0.483, 說明熱帶印度洋內(nèi)部SSTA變化一致。這與前的研究結(jié)果也基本相同, 譚言科[17]等發(fā)現(xiàn)印度洋海溫的變化主要存在全區(qū)符號(hào)一致的特點(diǎn),也有學(xué)者把熱帶印度洋呈現(xiàn)的“海盆一致”稱為印度洋海盆模態(tài)(Indian Ocean Basin Mode, IOBM), 它是印度洋對(duì)太平洋的ENSO響應(yīng)模態(tài)[18-19]。

    總的來說, 在全球變暖背景的驅(qū)動(dòng)下, 中國遠(yuǎn)洋作業(yè)漁場(chǎng)SST也逐漸變暖。而引起SST變動(dòng)的因素較為復(fù)雜, 不同區(qū)域會(huì)受到該區(qū)域的海流、盛行風(fēng)以及大氣環(huán)流等因素的相互作用, 從而使得不同區(qū)域海表溫度變動(dòng)的程度有所不同, 厄爾尼諾和拉尼娜現(xiàn)象就是一個(gè)典型的例子。印度洋的SST變動(dòng)對(duì)ENSO事件響應(yīng)及關(guān)系目前還存在一些爭論, 有一些學(xué)者認(rèn)為全球海氣是一個(gè)整體, ENSO與印度洋海溫異常密切相關(guān)[20-21]。與此相反, 也有學(xué)者認(rèn)為兩者相互獨(dú)立, 是各自的動(dòng)力作用和熱力交換的產(chǎn)物[22]。事實(shí)上印度洋發(fā)生較強(qiáng)的SST異常確實(shí)可以發(fā)生在非厄爾尼諾的年份, 其異常變動(dòng)并不完全依賴于太平洋的ENSO事件, 意味著印度洋和太平洋之間的海溫異常變動(dòng)既存在一定的獨(dú)立性, 又存在一定的相關(guān)性。此外, 通過對(duì)比大陸東西海岸的漁場(chǎng)SSTA變化, 可以發(fā)現(xiàn)大陸東海岸漁場(chǎng)(西北太平洋漁場(chǎng)、西印度洋漁場(chǎng)以及西南大西洋漁場(chǎng))往往要比大陸西海岸漁場(chǎng)(東太平洋漁場(chǎng)和東南太平洋漁場(chǎng))SSTA上升更為顯著。這可能與洋流的分布有關(guān), 大陸的東海岸多以暖流為主, 而西海岸則多以寒流為主, 如北美洲西海岸的加利福尼亞寒流與南美洲西海岸的秘魯寒流使得東太平洋和東南太平洋漁場(chǎng)SSTA下降,而位于亞洲東海岸的日本暖流和北太平洋暖流、南美東海岸的巴西暖流以及非洲東海岸的厄加勒斯暖流則使得西北太平洋, 西印度洋和西南大西洋的SSTA上升。

    海表溫度是漁場(chǎng)環(huán)境的重要指標(biāo)之一, 海表溫度異常用來表征海表溫度的變化狀態(tài)和趨勢(shì)。海表溫度及其變化對(duì)魚類資源的分布、洄游和集群等行為都有著直接或間接的影響, 漁場(chǎng)SST變暖的趨勢(shì)一般會(huì)使暖水性中上層魚類漁獲量增多, 但也存在特殊的情況, 如厄爾尼諾現(xiàn)象使SST溫度異常升高,也會(huì)導(dǎo)致海水中浮游生物明顯減少, 從而導(dǎo)致以浮游生物為餌料的魚類大量死亡。因此, 海溫的變動(dòng)對(duì)漁場(chǎng)的形成關(guān)鍵而復(fù)雜, 需要進(jìn)一步考慮魚類的生活習(xí)性和生理特征。本文從漁場(chǎng)的角度來分析各海域的SSTA年際變動(dòng)及其對(duì)ENSO的響應(yīng), 為漁業(yè)資源管理和開發(fā)提供環(huán)境參考。BELKIN[14]文中的大海洋生態(tài)系統(tǒng)雖然也是從漁場(chǎng)經(jīng)營的角度劃分的, 但它大多靠近海岸洋流和上涌流或者是被陸地包圍的海域, 因此這種劃分方式忽略了一些離陸地較遠(yuǎn)的大洋漁場(chǎng)的海溫變化情況, 而文章按海域位置將中國遠(yuǎn)洋作業(yè)漁場(chǎng)分為西北太平洋漁場(chǎng)、東太平洋漁場(chǎng)、東南太平洋漁場(chǎng)、西太平洋漁場(chǎng)、印度洋漁場(chǎng)、中大西洋漁場(chǎng)和西南大西洋漁場(chǎng), 則充分考慮了非陸緣漁場(chǎng)的海溫變化情況。文章中海表溫度只是一個(gè)平均化的結(jié)果, 每個(gè)漁場(chǎng)的SST只用該空間范圍的算術(shù)平均值來表示, 因此只能說明整體情況而不能揭示出局部或空間尺度上的變化特征。此外, 本文只探討了各漁場(chǎng)SSTA的年際變化與ENSO的聯(lián)系,但并未深入研究各漁場(chǎng)SSTA的年際變異機(jī)理, 關(guān)于大氣(風(fēng)、太陽輻射、降水等)和海洋要素(水平平流、垂向流、溫躍層等)對(duì)各漁場(chǎng)SSTA年際變化的影響等方面的研究有待今后擬另文分別敘述, 為不同的漁場(chǎng)中長期預(yù)報(bào)提供參考。

    [1] 秦大河. 氣候變化科學(xué)與人類可持續(xù)發(fā)展[J]. 地理科學(xué)進(jìn)展, 2014, 33(7): 874-883. Qin Dahe. Climate change science and sustainable development[J]. Progress in Geography, 2014, 33(7):874-883.

    [2] Intergovernmental Panel on Climate Change. Working Group II. Climate Change 2014: Impacts, Adaptation,and Vulnerability[M]. New York: Cambridge University Press, 2014: 411-484.

    [3] 周甦芳. 厄爾尼諾-南方濤動(dòng)現(xiàn)象對(duì)中西太平洋鰹魚圍網(wǎng)漁場(chǎng)的影響[J]. 中國水產(chǎn)科學(xué), 2006, 12(6):739-744. Zhou Sufang. Impacts of the El Ni(n)o Southern Oscillation on skipjack tuna purse-seine fishing grounds in the Western and Central Pacific Ocean[J]. Journal of Fishery Sciences of China, 2006, 12(6): 739-744.

    [4] 崔科, 陳新軍. 東黃海鮐魚資源豐度與表溫關(guān)系[J].南方水產(chǎn)科學(xué), 2007, 3(4): 20-25. Cui Ke, Chen Xinjun. Study of the relationships between SST and mackerel abundances in the Yellow and East China Seas[J]. South China Fisheries Science,2007, 3(4): 20-25.]

    [5] 曹曉怡, 周為峰, 樊偉, 等. 印度洋大眼金槍魚, 黃鰭金槍魚延繩釣漁場(chǎng)重心變化分析[J]. 上海海洋大學(xué)學(xué)報(bào), 2009, 18(4): 466-471. Cao Xiaoyi, Zhou Weifeng, Fan Wei, et al. Analysis on barycenter of fishing ground for bigeye tuna and ellowfin tuna by longlining in the Indian Ocean[J]. Journal of Shanghai University, 2009, 18(4): 466-471.

    [6] 汪金濤, 陳新軍. 中西太平洋鰹魚漁場(chǎng)的重心變化及其預(yù)測(cè)模型建立[J]. 中國海洋大學(xué)學(xué)報(bào): 自然科學(xué)版, 2013, 43(8): 44-48. Wang Jintao, Chen Xinjun. Changes and Prediction of the Fishing Ground Gravity of Skipjack (Katsuwonus pelamis) in Western-Central Pacific[J]. Periodical of Ocean University of China, 2013, 43(8): 44-48.

    [7] FAO 1990-2015. CWP Handbook of Fishery Statistical Standards. Section H: FISHING AREAS FOR STATISTICAL PURPOSES. CWP Data Collection. In:FAO Fisheries and Aquaculture Department [online]. Rome. Updated 10 January 2002. [Cited 4 February 2015]. http: //www.fao.org/fishery/cwp/handbook/H/en.

    [8] 王國杰, 姜彤, 陳桂亞. 長江干流徑流的時(shí)序結(jié)構(gòu)與長期記憶[J]. 地理學(xué)報(bào), 2006, 61(1): 47-56. Wang Guojie, Jiang Tong, Chen Guiya. Structure and Long-term Memory of Discharge Series in Yangtze River[J]. Acta Geographica Sinica, 2006, 61(1): 47-56.]

    [9] Alexander M A, Blade I, Newman M, et al. The atmospheric bridge: The influence of ENSO teleconnections on air-sea interaction over the global oceans[J]. J Climate, 2002, 15(16): 2205-2231.

    [10] Liu Z Y, Alexander M. Atmospheric Bridge, oceanic tunnel, and global climatic teleconnections[J]. Reviews of Geophysics, 2007, 45(2): 1-34.

    [11] Klein S A, Soden b J, Lau N C. Remote sea surface temperature variations during ENSO: Evidence for a tropical atmospheric bridge[J]. J Climate, 1999, 12(4):917-932.

    [12] Susanto R D, Gordon A L, Zhang Q A. Upwelling along the coasts of Java and Sumatra and its relation to ENSO[J]. Geophysical Research Letters, 2001, 28(8):1599-1602.

    [13] 劉巖松, 王法明. 一個(gè)簡單的印-太海氣耦合模式[J].海洋與湖沼, 2013, 44(6): 1462-1468. Liu Yansong, Wang Faming. An Intermedia Air-Sea Coupled Model for the Indo-Pacific Ocean[J]. Oceanologia et Limnologia Sinica, 2013, 44(6): 1462-1468.

    [14] Belkin I M. Rapid warming of large marine ecosystems[J]. Progress in Oceanography, 2009, 81(1): 207-213.

    [15] Lau N C, Nath m J. The role of the “atmospheric bridge” in linking tropical Pacific ENSO events to extra tropical SST anomalies[J]. J Climate, 1996, 9(9):2036-2057.

    [16] 江志紅, 張強(qiáng), 屠其璞, 等. 南北半球平均溫度年代際振蕩的非均衡性及其對(duì)變暖顯著性的影響[J]. 熱帶氣象學(xué)報(bào), 2003, 19(2): 122-136. Jiang Zhihong, Zhang Qiang, Tu Qipu, et al. Differences of Interdacadal Oscillation between Mean Temperature of Northern and Southern Hemisphere and Their Influences on Warming Signal[J]. Journal of Tropical Meteorology, 2003, 19(2): 122-136.

    [17] 譚言科, 張人禾, 何金海, 等. 熱帶印度洋海溫的年際變化與ENSO[J]. 氣象學(xué)報(bào), 2004, 62(6): 831-840. Tan Yan-ke, Zhang Ren-he, He Jin-hai, et al. Relationship of the Interannual Variations of Sea Surface Temperature in Tropical Indian Ocean to Enso[J]. Acta Meteorologica Sinica, 2004, 62(6): 831-840.

    [18] Vinayachandran P N, Goswami B N, Saji N H, et al. A dipole mode in the tropical Indian Ocean[J]. Nature International Weekly Journal of Science, 1999, 401(6751):360-363.

    [19] Yang J L, Liu Q Y, Xie S P, et al. Impact of the Indian Ocean SST basin mode on the Asian summer monsoon[J]. Geophysical Research Letters, 2007, 34(2):155-164.

    [20] Du Y, Xie S P, Huang G, et al. Role of Air-Sea Interaction in the Long Persistence of El Ni?o-Induced North Indian Ocean Warming[J]. J Climate, 2009,22(8): 2023-2038.

    [21] 范伶俐, 郭品文. 熱帶太平洋-印度洋相互關(guān)系的年代際改變[J]. 熱帶氣象學(xué)報(bào), 2005, 21(2): 172-180. FAN Ling-li, Guo Pin-wen. The Interdecadal Variation of the Interrelation between the Tropical Pacific Ocean and Indian Ocean [J]. Journal of Tropical Meteorology,2005, 21(2): 172-180.

    [22] Valsala V K, Ikeda M. Pathways and effects of the Indonesian Through flow water in the Indian Ocean using particle trajectory and tracers in an OGCM[J]. J Climate, 2007, 20(13): 2994-3017.

    (本文編輯: 梁德海)

    Sea surface temperature anomaly’s interannual variability in pelagic fishing grounds of China

    JI Shi-jian1, 2, ZHOU Wei-feng1, FAN Wei1, JIN Shao-fei3, CUI Xue-sen1
    (1. Key Laboratory of East China Sea & Oceanic Fishery Resources Exploitation and Utilization, Ministry of Agriculture, Shanghai 200090, China; 2. College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China; 3. Key Laboratory of Regional Climate-Environment for Temperate East Asia, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China)

    May, 25, 2015

    pelagic fishing ground; sea surface temperature; sea surface temperature anomaly

    Based on monthly averaged sea surface temperature data from the U.S. National Oceanic and Atmospheric Administration’s Climate Prediction Center, we conducted a statistical analysis of sea surface temperature anomalies (SSTA) of China’s seven main pelagic fishing grounds, the East Pacific, West Pacific, Middle Atlantic,North Pacific, Southeast Pacific, Southwest Atlantic, and Indian Ocean fishing grounds over a 30-year time period (1982–2011) using the power spectrum method to calculate remarkable variation periods in the time series. We also performed a correlation analysis to determine the correlation between trend-removed SSTA and the southern oscillation index (SOI). The results show that the SSTA of China’s main pelagic fishing grounds fluctuate between –0.3℃and 0.3℃ with a 3–4-year period and demonstrate a general rising trend over the 30-year period. A significant correlation (r = ?0.509) exists between the trend-removed SSTA and SOI, indicating that the SSTA of the fishing grounds are closely related to ENSO events. In terms of the SSTA fluctuation trend of the fishing grounds, we found that except for the downward SSTA trends in the East Pacific and Southeast Pacific, the other fishing grounds’SSTA increases to a certain degree, especially in the North Pacific. In terms of the fishing ground SSTA fluctuation period, the SSTA variation period is 34 years in the East Pacific as well as in the Southwest Atlantic, 4 years in the Southeast Pacific, 5 years in the West Atlantic, and 10 years in the other areas. The correlation analysis between the trend-removed SSTA of each fishing ground and the SOI indicates that the SOI is correlated to the SSTA of the East Pacific, middle Atlantic, and Southwest Atlantic, with correlation coefficients of –0.895, 0.471, and –0.598, respectively, and no significant correlation was observed for the others. A comparison of the fishing grounds reveals the following features. At the equator, the SSTA change of the East Pacific fishing ground is always opposite to that of the West Pacific and Middle Atlantic fishing grounds but is in-phase with that of the Indian Ocean fishing ground. Among the three fishing grounds in the mid-latitudes, SSTA changes in the Southern Hemisphere latitudes are more stable than those in the northern latitudes. A significant correlation exists between the East Indian and West Indian Oceans, and positive and negative changes in the SSTAs of the entire Indian Ocean are basically the same.

    S931.3

    A

    1000-3096(2016)01-0085-09

    10.11759/hykx20150313003

    2015-05-25;

    2015-10-11

    國家科技支撐計(jì)劃項(xiàng)目(2013BAD13B06); 上海市科學(xué)技術(shù)委員會(huì)科研計(jì)劃項(xiàng)目(15DZ1202201)

    紀(jì)世建(1991-), 男, 碩士研究生, 主要從事南海漁情預(yù)報(bào)研究, E-mail: jackjsj@foxmail.com;周為峰(1978-), 通信作者, 女, 副研究員, 主要從事漁業(yè)遙感與漁場(chǎng)海洋學(xué)研究, E-mail: zhwfzhwf@163.com

    猜你喜歡
    海表海溫漁場(chǎng)
    《關(guān)于加強(qiáng)遠(yuǎn)洋魷釣漁船作業(yè)管理的通知》公布
    基于深度學(xué)習(xí)的海溫觀測(cè)數(shù)據(jù)質(zhì)量控制應(yīng)用研究
    基于Argo、XBT數(shù)據(jù)的蘇拉威西海溫鹽特征分析
    埃及超大規(guī)模養(yǎng)魚場(chǎng)
    基于無人機(jī)的海表環(huán)境智能監(jiān)測(cè)系統(tǒng)設(shè)計(jì)與應(yīng)用
    南印度洋偶極型海溫與中國西南地區(qū)初秋降水的關(guān)系
    2016與1998年春季北大西洋海表溫度異常的差異及成因
    融合海表溫度產(chǎn)品在渤黃東海的對(duì)比分析及初步驗(yàn)證
    太陽總輻照度對(duì)熱帶中太平洋海表溫度年代際變化的可能影響
    連續(xù)三年保持50%增長,入駐通威、新希望養(yǎng)殖基地,這家漁機(jī)企業(yè)將打造未來無人漁場(chǎng)
    18禁黄网站禁片免费观看直播| 免费看光身美女| 国产亚洲精品综合一区在线观看| 精品国产亚洲在线| av中文乱码字幕在线| 亚洲国产看品久久| 日韩成人在线观看一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 桃红色精品国产亚洲av| 首页视频小说图片口味搜索| 午夜精品一区二区三区免费看| 天天一区二区日本电影三级| 国产三级中文精品| 99在线人妻在线中文字幕| 99国产综合亚洲精品| 午夜福利视频1000在线观看| 99久久综合精品五月天人人| 99国产精品99久久久久| 亚洲成人免费电影在线观看| 日韩成人在线观看一区二区三区| 国产成年人精品一区二区| 国产精品av久久久久免费| 亚洲五月婷婷丁香| 亚洲第一电影网av| 亚洲av五月六月丁香网| 精品人妻1区二区| 99热只有精品国产| 18美女黄网站色大片免费观看| xxxwww97欧美| 精品国内亚洲2022精品成人| 黄色日韩在线| 男女下面进入的视频免费午夜| 90打野战视频偷拍视频| 成人av在线播放网站| 国产久久久一区二区三区| 午夜日韩欧美国产| 国产乱人视频| 十八禁网站免费在线| 日本a在线网址| 欧美日韩一级在线毛片| 最近最新中文字幕大全电影3| 亚洲专区字幕在线| 2021天堂中文幕一二区在线观| 99精品欧美一区二区三区四区| 天天添夜夜摸| 可以在线观看的亚洲视频| 国产精品 欧美亚洲| 99国产精品99久久久久| 一区二区三区激情视频| 99re在线观看精品视频| 黄色视频,在线免费观看| 毛片女人毛片| 亚洲精品乱码久久久v下载方式 | 欧美乱码精品一区二区三区| 亚洲成人久久性| 中文字幕人妻丝袜一区二区| 日韩有码中文字幕| 国产1区2区3区精品| 麻豆久久精品国产亚洲av| 深夜精品福利| 欧美日韩一级在线毛片| 国产熟女xx| 精品人妻1区二区| 国产精品自产拍在线观看55亚洲| 性欧美人与动物交配| 亚洲无线观看免费| 99re在线观看精品视频| 久久精品91无色码中文字幕| 欧美精品啪啪一区二区三区| 青草久久国产| 法律面前人人平等表现在哪些方面| 日韩欧美在线乱码| 一二三四在线观看免费中文在| 最好的美女福利视频网| 淫秽高清视频在线观看| 久久久色成人| 国产在线精品亚洲第一网站| 日韩国内少妇激情av| 国产精品亚洲美女久久久| 在线观看免费视频日本深夜| 国产极品精品免费视频能看的| 草草在线视频免费看| 精品久久久久久久毛片微露脸| 精品欧美国产一区二区三| 成年女人毛片免费观看观看9| 天天一区二区日本电影三级| 亚洲熟女毛片儿| 男人舔女人下体高潮全视频| 午夜福利在线观看吧| 成人亚洲精品av一区二区| 丁香六月欧美| 日韩免费av在线播放| 久久久久久久久久黄片| 中亚洲国语对白在线视频| 婷婷丁香在线五月| 欧美不卡视频在线免费观看| 国产伦精品一区二区三区视频9 | 免费一级毛片在线播放高清视频| 亚洲人成网站在线播放欧美日韩| 香蕉av资源在线| 99国产精品一区二区蜜桃av| 中文字幕精品亚洲无线码一区| 99热这里只有是精品50| 久久香蕉精品热| 成人三级黄色视频| 18禁美女被吸乳视频| 欧美3d第一页| 亚洲专区字幕在线| 一a级毛片在线观看| 免费观看人在逋| aaaaa片日本免费| 国产探花在线观看一区二区| 伊人久久大香线蕉亚洲五| 一二三四社区在线视频社区8| 久久天躁狠狠躁夜夜2o2o| 婷婷丁香在线五月| 国产成人精品久久二区二区91| 久久草成人影院| 韩国av一区二区三区四区| 午夜福利高清视频| 美女高潮喷水抽搐中文字幕| 亚洲五月婷婷丁香| 美女被艹到高潮喷水动态| 我的老师免费观看完整版| 日本免费一区二区三区高清不卡| 日本 av在线| 日日干狠狠操夜夜爽| 欧美中文日本在线观看视频| 久久久久国产精品人妻aⅴ院| 久久国产乱子伦精品免费另类| or卡值多少钱| 国产69精品久久久久777片 | 最新中文字幕久久久久 | 精品午夜福利视频在线观看一区| 久久久久国产精品人妻aⅴ院| 一级a爱片免费观看的视频| 亚洲性夜色夜夜综合| 午夜激情欧美在线| 亚洲国产日韩欧美精品在线观看 | 伊人久久大香线蕉亚洲五| a级毛片a级免费在线| 欧美成狂野欧美在线观看| 特大巨黑吊av在线直播| 99国产精品一区二区三区| 亚洲精品在线观看二区| 天天添夜夜摸| 亚洲男人的天堂狠狠| 岛国在线观看网站| 免费av毛片视频| 麻豆av在线久日| 欧美极品一区二区三区四区| 欧美激情久久久久久爽电影| 久久精品国产清高在天天线| 嫩草影院入口| 淫妇啪啪啪对白视频| 每晚都被弄得嗷嗷叫到高潮| 成年女人永久免费观看视频| 成人鲁丝片一二三区免费| 国产精品99久久99久久久不卡| 在线免费观看的www视频| 在线观看免费午夜福利视频| 免费无遮挡裸体视频| 精华霜和精华液先用哪个| 日日干狠狠操夜夜爽| 亚洲男人的天堂狠狠| av福利片在线观看| 一区福利在线观看| 欧美激情在线99| 校园春色视频在线观看| 一个人看的www免费观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲五月天丁香| 国产欧美日韩精品亚洲av| 色老头精品视频在线观看| 亚洲中文字幕日韩| 变态另类丝袜制服| 性欧美人与动物交配| 欧美精品啪啪一区二区三区| 亚洲av美国av| 一级作爱视频免费观看| 99国产精品99久久久久| 精品久久久久久,| 国产精品爽爽va在线观看网站| 亚洲国产色片| 亚洲中文日韩欧美视频| 精品国产三级普通话版| 亚洲va日本ⅴa欧美va伊人久久| 国产极品精品免费视频能看的| 视频区欧美日本亚洲| 欧美一级毛片孕妇| 身体一侧抽搐| 99久久无色码亚洲精品果冻| 国产1区2区3区精品| av国产免费在线观看| 男人的好看免费观看在线视频| 看黄色毛片网站| 午夜福利欧美成人| 一进一出抽搐动态| 九九在线视频观看精品| 久久久国产精品麻豆| 美女扒开内裤让男人捅视频| 国产精品影院久久| 啪啪无遮挡十八禁网站| 亚洲成av人片在线播放无| 国产精品av视频在线免费观看| 国产欧美日韩一区二区精品| 亚洲专区国产一区二区| 久久中文字幕人妻熟女| 国产激情欧美一区二区| 热99在线观看视频| 在线观看免费午夜福利视频| 欧美日韩综合久久久久久 | 亚洲色图 男人天堂 中文字幕| 欧美乱色亚洲激情| 国产高清视频在线观看网站| 日本三级黄在线观看| 999久久久国产精品视频| 亚洲精品中文字幕一二三四区| 国产成人av激情在线播放| av中文乱码字幕在线| 日韩高清综合在线| 91在线观看av| 欧美黄色片欧美黄色片| 亚洲av电影不卡..在线观看| 人妻夜夜爽99麻豆av| 91在线精品国自产拍蜜月 | 999精品在线视频| 制服丝袜大香蕉在线| 国产 一区 欧美 日韩| 99国产极品粉嫩在线观看| 欧美xxxx黑人xx丫x性爽| 欧美高清成人免费视频www| 男女做爰动态图高潮gif福利片| 欧美日韩福利视频一区二区| 国产精品1区2区在线观看.| 欧美一区二区国产精品久久精品| 精品一区二区三区av网在线观看| 亚洲精品一区av在线观看| 亚洲五月婷婷丁香| 动漫黄色视频在线观看| 中文字幕熟女人妻在线| 中文字幕久久专区| 美女黄网站色视频| 狂野欧美激情性xxxx| 两个人的视频大全免费| 亚洲精品久久国产高清桃花| 757午夜福利合集在线观看| 亚洲精品粉嫩美女一区| 久久热在线av| 久久婷婷人人爽人人干人人爱| 无人区码免费观看不卡| 色哟哟哟哟哟哟| 亚洲欧美日韩卡通动漫| 国产亚洲精品av在线| 免费电影在线观看免费观看| 日本三级黄在线观看| av在线蜜桃| 99久久无色码亚洲精品果冻| 黑人操中国人逼视频| 国产精品免费一区二区三区在线| 久久草成人影院| 18禁观看日本| 国产精品久久视频播放| 成人国产综合亚洲| 香蕉丝袜av| 一级毛片精品| 色吧在线观看| 久久久久国产一级毛片高清牌| 精品电影一区二区在线| 欧美成狂野欧美在线观看| 不卡av一区二区三区| 国产一区二区在线av高清观看| 国产成人啪精品午夜网站| 日韩欧美国产一区二区入口| 一级作爱视频免费观看| 999久久久精品免费观看国产| 亚洲欧美日韩卡通动漫| 香蕉久久夜色| 亚洲精品色激情综合| 亚洲第一欧美日韩一区二区三区| 嫩草影视91久久| 毛片女人毛片| 一区二区三区激情视频| 啪啪无遮挡十八禁网站| 欧美激情久久久久久爽电影| 日本a在线网址| 看黄色毛片网站| 欧美黑人欧美精品刺激| 天堂影院成人在线观看| 久久99热这里只有精品18| 久久久国产成人精品二区| 国产高清视频在线播放一区| 亚洲成av人片免费观看| 熟女电影av网| 久久天躁狠狠躁夜夜2o2o| 老司机深夜福利视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日日夜夜操网爽| 国产精品,欧美在线| 成人av一区二区三区在线看| 久久久国产精品麻豆| xxx96com| 久久久久久久久久黄片| 国产 一区 欧美 日韩| 热99在线观看视频| www.精华液| 久久久国产欧美日韩av| 欧美黄色淫秽网站| 2021天堂中文幕一二区在线观| 久久久国产成人免费| 亚洲中文字幕日韩| 99riav亚洲国产免费| 成年版毛片免费区| 99国产精品一区二区蜜桃av| 久久久久国产精品人妻aⅴ院| av在线天堂中文字幕| 天堂av国产一区二区熟女人妻| 国产精品免费一区二区三区在线| 一本精品99久久精品77| 欧美午夜高清在线| 国产精品永久免费网站| 久久婷婷人人爽人人干人人爱| 亚洲第一欧美日韩一区二区三区| 搡老岳熟女国产| 桃色一区二区三区在线观看| 国产野战对白在线观看| 亚洲成a人片在线一区二区| 91字幕亚洲| 午夜精品久久久久久毛片777| 午夜两性在线视频| 在线观看午夜福利视频| 久久国产精品影院| 18美女黄网站色大片免费观看| 亚洲午夜精品一区,二区,三区| 91av网站免费观看| 国产伦一二天堂av在线观看| 91av网站免费观看| 99国产精品一区二区三区| 极品教师在线免费播放| 成人精品一区二区免费| 午夜福利免费观看在线| 亚洲av片天天在线观看| 综合色av麻豆| 欧美高清成人免费视频www| 成年女人永久免费观看视频| 一区二区三区国产精品乱码| 麻豆国产97在线/欧美| 99久久成人亚洲精品观看| 日韩免费av在线播放| 亚洲精华国产精华精| 欧美绝顶高潮抽搐喷水| 国产高清视频在线观看网站| 人妻久久中文字幕网| 少妇人妻一区二区三区视频| 久久久国产成人精品二区| 亚洲精品在线观看二区| 精品无人区乱码1区二区| 精品电影一区二区在线| 久99久视频精品免费| 九九在线视频观看精品| 欧美乱妇无乱码| 全区人妻精品视频| 给我免费播放毛片高清在线观看| 亚洲精品在线美女| 亚洲18禁久久av| 色av中文字幕| 男女做爰动态图高潮gif福利片| 91麻豆av在线| 男女下面进入的视频免费午夜| 校园春色视频在线观看| av欧美777| 国产精品久久电影中文字幕| 他把我摸到了高潮在线观看| 国产 一区 欧美 日韩| 黄片小视频在线播放| 亚洲天堂国产精品一区在线| 男女床上黄色一级片免费看| 窝窝影院91人妻| 国产午夜精品久久久久久| 人妻丰满熟妇av一区二区三区| 男女床上黄色一级片免费看| 成人欧美大片| 久久久久久大精品| 好看av亚洲va欧美ⅴa在| 亚洲九九香蕉| 免费看光身美女| 岛国在线免费视频观看| 国语自产精品视频在线第100页| 国产精品久久久久久人妻精品电影| 国产精品乱码一区二三区的特点| 免费搜索国产男女视频| 在线a可以看的网站| 亚洲精华国产精华精| 在线a可以看的网站| www.www免费av| 九九在线视频观看精品| 99在线视频只有这里精品首页| 成人av一区二区三区在线看| 很黄的视频免费| 人人妻人人澡欧美一区二区| 色综合欧美亚洲国产小说| 一级a爱片免费观看的视频| 精品日产1卡2卡| 五月伊人婷婷丁香| 欧美国产日韩亚洲一区| 在线免费观看的www视频| 国产成人av激情在线播放| 老汉色av国产亚洲站长工具| 在线国产一区二区在线| 精品免费久久久久久久清纯| 亚洲色图av天堂| 免费看十八禁软件| 色哟哟哟哟哟哟| 操出白浆在线播放| 欧美日本亚洲视频在线播放| 一卡2卡三卡四卡精品乱码亚洲| 久久精品aⅴ一区二区三区四区| 精品久久久久久久久久久久久| 97超视频在线观看视频| 综合色av麻豆| 搡老熟女国产l中国老女人| 男女之事视频高清在线观看| 亚洲精品一卡2卡三卡4卡5卡| 18禁裸乳无遮挡免费网站照片| 精品国产超薄肉色丝袜足j| 99在线视频只有这里精品首页| 大型黄色视频在线免费观看| 国产单亲对白刺激| 久久精品国产亚洲av香蕉五月| 丰满人妻熟妇乱又伦精品不卡| 在线免费观看的www视频| 中文字幕久久专区| 天天添夜夜摸| 曰老女人黄片| 给我免费播放毛片高清在线观看| 欧美+亚洲+日韩+国产| 最新在线观看一区二区三区| 精品国产三级普通话版| 观看免费一级毛片| 19禁男女啪啪无遮挡网站| 亚洲av成人一区二区三| 一个人免费在线观看的高清视频| 日韩国内少妇激情av| 老鸭窝网址在线观看| 全区人妻精品视频| 99久久国产精品久久久| 亚洲成人中文字幕在线播放| 国产精品永久免费网站| 哪里可以看免费的av片| netflix在线观看网站| 精品免费久久久久久久清纯| 黑人巨大精品欧美一区二区mp4| 精品久久久久久久人妻蜜臀av| 国产av不卡久久| 日韩欧美在线乱码| 久久久久久久精品吃奶| 国产精品久久久久久久电影 | 亚洲av成人精品一区久久| 亚洲av成人av| 色综合婷婷激情| 日本免费a在线| 老熟妇乱子伦视频在线观看| 男人的好看免费观看在线视频| 一区二区三区高清视频在线| 日本黄色片子视频| 婷婷六月久久综合丁香| 99精品在免费线老司机午夜| 91av网站免费观看| 黑人巨大精品欧美一区二区mp4| 精品一区二区三区视频在线 | 噜噜噜噜噜久久久久久91| 久久精品亚洲精品国产色婷小说| 中国美女看黄片| 一区二区三区高清视频在线| 欧美丝袜亚洲另类 | 色尼玛亚洲综合影院| 岛国视频午夜一区免费看| 天堂av国产一区二区熟女人妻| 国产黄色小视频在线观看| 黄色成人免费大全| 别揉我奶头~嗯~啊~动态视频| 久久久久亚洲av毛片大全| 首页视频小说图片口味搜索| 国产欧美日韩一区二区三| 精品国内亚洲2022精品成人| 男女之事视频高清在线观看| 国产精品久久久久久久电影 | 亚洲av成人精品一区久久| 一个人免费在线观看的高清视频| 亚洲色图av天堂| 欧美色视频一区免费| 村上凉子中文字幕在线| 丰满人妻熟妇乱又伦精品不卡| 变态另类丝袜制服| 亚洲最大成人中文| 一本一本综合久久| 午夜免费成人在线视频| 国产精品 欧美亚洲| 久久久水蜜桃国产精品网| 国产熟女xx| 淫秽高清视频在线观看| 国产精品99久久99久久久不卡| 欧美色视频一区免费| 亚洲成人精品中文字幕电影| 最近最新中文字幕大全电影3| 90打野战视频偷拍视频| 国产精品一区二区精品视频观看| 亚洲国产色片| 国产一区二区在线观看日韩 | 欧美日本视频| 最近最新免费中文字幕在线| 国产乱人视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品九九99| 老司机深夜福利视频在线观看| 高清毛片免费观看视频网站| 99国产极品粉嫩在线观看| 国产精品久久电影中文字幕| 久久香蕉国产精品| 亚洲国产精品合色在线| 真人一进一出gif抽搐免费| 午夜影院日韩av| 级片在线观看| 国产精品精品国产色婷婷| 精品不卡国产一区二区三区| 很黄的视频免费| 国产精品久久久久久精品电影| 成人欧美大片| 色播亚洲综合网| 免费在线观看影片大全网站| 精品午夜福利视频在线观看一区| 久久精品影院6| 亚洲欧美日韩东京热| 美女cb高潮喷水在线观看 | 级片在线观看| 亚洲性夜色夜夜综合| 在线观看免费视频日本深夜| 国产日本99.免费观看| 欧美午夜高清在线| 91九色精品人成在线观看| a在线观看视频网站| 精品人妻1区二区| 中文字幕最新亚洲高清| 最近视频中文字幕2019在线8| 小说图片视频综合网站| 桃色一区二区三区在线观看| 亚洲av免费在线观看| 欧美绝顶高潮抽搐喷水| 麻豆国产av国片精品| 国产一区在线观看成人免费| 美女被艹到高潮喷水动态| 亚洲黑人精品在线| 欧美一区二区精品小视频在线| 亚洲avbb在线观看| 视频区欧美日本亚洲| 在线观看66精品国产| 亚洲人成伊人成综合网2020| 久久中文看片网| 欧美乱色亚洲激情| 亚洲色图 男人天堂 中文字幕| 啦啦啦免费观看视频1| 成人欧美大片| 亚洲五月天丁香| 午夜免费成人在线视频| 最新美女视频免费是黄的| 免费大片18禁| 国产极品精品免费视频能看的| 国产v大片淫在线免费观看| 九色成人免费人妻av| 曰老女人黄片| 国产成人一区二区三区免费视频网站| 久久99热这里只有精品18| 韩国av一区二区三区四区| 久久性视频一级片| 美女大奶头视频| 日韩大尺度精品在线看网址| av视频在线观看入口| 大型黄色视频在线免费观看| 级片在线观看| 国产伦人伦偷精品视频| 免费在线观看影片大全网站| 欧洲精品卡2卡3卡4卡5卡区| 国产精品自产拍在线观看55亚洲| 精品一区二区三区视频在线观看免费| 成年女人毛片免费观看观看9| 久久久水蜜桃国产精品网| 国产精品一区二区精品视频观看| 久久久久性生活片| 亚洲,欧美精品.| 少妇熟女aⅴ在线视频| 别揉我奶头~嗯~啊~动态视频| 久久这里只有精品中国| 在线国产一区二区在线| 午夜福利视频1000在线观看| www日本黄色视频网| 真人做人爱边吃奶动态| 俺也久久电影网| 欧美在线黄色| 在线国产一区二区在线| 国产亚洲精品久久久com| 少妇丰满av| 在线国产一区二区在线| 久久九九热精品免费| 五月伊人婷婷丁香| 成人精品一区二区免费| 免费观看的影片在线观看| 啦啦啦免费观看视频1| 中文字幕av在线有码专区| 青草久久国产| 午夜a级毛片| 国产熟女xx| 午夜精品久久久久久毛片777| 久久久成人免费电影| 欧美日韩国产亚洲二区| 成人性生交大片免费视频hd|