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

    西太平洋暖池海域上層海洋熱鹽含量初步研究*

    2016-05-12 03:03:37楊小欣吳曉芬劉增宏

    楊小欣, 吳曉芬, 劉增宏

    (國家海洋局第二海洋研究所衛(wèi)星海洋環(huán)境動力學國家重點實驗室,浙江 杭州 310012)

    ?

    西太平洋暖池海域上層海洋熱鹽含量初步研究*

    楊小欣, 吳曉芬, 劉增宏

    (國家海洋局第二海洋研究所衛(wèi)星海洋環(huán)境動力學國家重點實驗室,浙江 杭州 310012)

    摘要:基于2001年1月—2014年7月期間的Argo溫鹽剖面資料,利用循環(huán)平穩(wěn)經(jīng)驗正交函數(shù)(CSEOF)分解、最大熵譜分析和相關分析等方法,研究了西太平洋暖池海域上層海洋熱鹽含量的空間分布、季節(jié)和年際變化特征,并探討了其可能的變化原因。結果表明,暖池海域次表層熱含量的變化要遠遠大于表層,鹽含量以表層變化為主,次表層變化不大,且熱含量受ENSO事件影響較大,而鹽含量受ENSO事件的影響并不大,以氣候態(tài)變化為主;CSEOF分析表明,暖池海域熱含量第一模態(tài)空間場具有顯著的東-西反相位年際振蕩,鹽含量第一模態(tài)則呈正-負-正的三極子變化模態(tài),時間序列顯示,熱含量在2007年以后伴隨ENSO事件經(jīng)過了3次位相調整,鹽含量在2007年以后只經(jīng)過一次位相調整。進一步分析表明,熱含量變化主要受到ENSO、局地風和緯向流的影響,而引起鹽含量變化的原因更為復雜,其中緯向流和淡水通量對研究海域鹽含量的影響均為正反饋效應。

    關鍵詞:熱含量;鹽含量;CSEOF分析;Argo資料;西太平洋暖池

    引用格式:楊小欣, 吳曉芬, 劉增宏. 西太平洋暖池海域上層海洋熱鹽含量初步研究[J]. 中國海洋大學學報(自然科學版), 2016, 46(4): 1-12.

    YANG Xiao-Xin, WU Xiao-Fen, LIU Zeng-Hong. Preliminary research on upper ocean heat and salt content of the West Pacific Warm Pool[J]. Periodical of Ocean University of China, 2016, 46(4): 1-12.

    西太平洋暖池海域是全球大洋表層溫度最高的海域,也是全球大氣運動的主要熱源之一,其變化對全球大氣和海洋起著十分重要的作用,歷來是大氣和海洋科學研究的熱點。1980年代以前,海洋和大氣環(huán)境要素的現(xiàn)場觀測資料相當有限,人們無法可靠地對發(fā)生在海洋上的眾多物理海洋現(xiàn)象進行監(jiān)測并及時提供預警預報。自1985年開始,國際上組織實施了為期10年的熱帶海洋-全球大氣計劃(TOGA)及“熱帶海洋全球大氣耦合響應試驗(TOGA-COARE),以及1990年開始的世界大洋環(huán)流實驗計劃(WOCE),同樣進行了為期10年的全球準同步觀測。期間,一些沿海國家也相繼開展了許多區(qū)域性的專題調查,取得了一批寶貴的海上第一手資料,從而使人們對發(fā)生在全球海洋上的熱點問題,特別是西太平洋暖池及其上層海洋熱鹽含量問題(諸如暖池和上層海洋熱含量的分布、變化及其相互之間的聯(lián)系,上層海洋熱含量與黑潮、ENSO等海氣相互作用事件間的關系)的分析、研究工作開始活躍起來。

    同時,大洋鹽含量的年際變化和ENSO事件也有著很強的相關性[12],ENSO導致的降雨異常會引起鹽含量變化,而大的鹽含量變化會引起熱鹽環(huán)流的顯著變化,熱鹽環(huán)流的變化則勢必會造成顯著的熱通量和氣候變化。有研究指出,在一定的溫度范圍內增加1g/L鹽含量對密度的影響相當于溫度降低3~4℃造成的影響[13]。對于鹽含量的變化機制,研究表明,在開闊的大洋內,亞熱帶海域較大的蒸發(fā)量以及兩極海域較大的降水量會導致大洋鹽含量的變化,同時在風、溫度等因素影響下的環(huán)流運動造成鹽含量的水平運輸,這種水平運輸在上層幾百米是十分迅速且顯著的,因此蒸發(fā)、降水以及洋流運輸都是造成大洋鹽含量分布的主要因素[14]。

    盡管目前國內外對海洋熱含量已經(jīng)有了豐富的研究成果,但對熱含量的估算仍然面臨著很大的挑戰(zhàn),如有關觀測空白區(qū)提出的不同mapping方案會帶來熱含量估算的不確定性[15-16],觀測儀器偏差(尤其是1966—2000年XBT資料出現(xiàn)偏差)帶來熱含量估算的不確定性,以及氣候態(tài)的選取(過去幾十年基于船舶觀測數(shù)據(jù)的時空分辨率不足以構造一個好的氣候態(tài))所帶來的熱含量分析誤差[17]等問題。同樣,由于歷史上大面積鹽度剖面資料的缺乏,以往對鹽含量研究的數(shù)據(jù)大都采用如SODA等海洋同化數(shù)據(jù)或者采用時空分布不均勻的走航觀測擬合所得[18],且目前對鹽含量的研究中都將鹽度對體積的積分值近似為鹽含量,而未考慮密度的影響等,故鹽含量的估算精度、對大范圍海洋鹽含量分布及其變化特征的研究還十分有限。

    2000年開始的全球Argo實時海洋觀測網(wǎng)建設給準確估算海盆尺度的海洋熱、鹽含量帶來了難得的機遇。Argo剖面浮標可以提供全球大洋0~2000dbar水深內的高時空分辨率的溫、鹽度資料,并于2007年11月初就已獲取了100萬條溫鹽度剖面。目前,正以每年12萬條剖面的速度在不斷遞增。Argo資料有效地補充了地球觀測衛(wèi)星所無法獲取的海洋次表層信息。Lyman等研究指出1955—2002年期間,由于觀測手段的提高,對熱含量的不確定值在Argo計劃實施之后達到了最低值[20]。在研究全球氣候變化,特別是深入探討西太平洋海域上層海洋熱鹽含量的分布、變化及其與暖池、ENSO等海洋現(xiàn)象和海氣相互作用事件之間的聯(lián)系,Argo觀測網(wǎng)有著比以往任何一個觀測系統(tǒng)更好的時空范圍的一致分布[21],為更加準確估算海洋上層熱鹽含量,探討熱鹽含量的季節(jié)和年際變化特征,深入研究年際尺度上ENSO事件與熱含量異常的關系及其海氣相互作用過程等問題提供了廣闊的前景。因此,本文擬利用已經(jīng)積累了14a的新穎Argo剖面資料,探討西太平洋暖池海域熱含量,特別是鹽含量的空間分布及季節(jié)和年際變化特征。

    1資料來源與分析方法

    1.1 資料來源

    (1)溫、鹽度剖面數(shù)據(jù)來自于日本海洋科學與技術廳(Japan Marine-Earth Science and Technology Center, JMESTC)提供的2001年1月—2014年7月間利用Argo浮標資料制作的逐月客觀分析網(wǎng)格化數(shù)據(jù)集(Grid Point Value of the Monthly Objective Analysis using Argo float data,MOAA GPV)。除了Argo浮標數(shù)據(jù)外,還采用了TRITON數(shù)據(jù)以及CTD觀測獲得的實測數(shù)據(jù),并使用最優(yōu)插值法得到水平分辨率為1(°)×1(°)經(jīng)、緯度的網(wǎng)格化數(shù)據(jù),垂向為10~2000dbar共25層,本文選取了西太平洋暖池海域(30°N~20°S,120°E~140°W),120dbar深度以上7層(10,20,30,50,75,100,120dbar)的溫鹽度數(shù)據(jù),并計算了該海域的熱鹽含量[22]。

    (3)海流資料為HYCOM模式海流數(shù)據(jù),本文選取的研究范圍為西太平洋暖池海域(30°N~20°S,120°E~140°W),120dbar深度以上21層,水平分辨率為1(°)×1(°)[23]。

    (4)降水量數(shù)據(jù)CMAP(CPC Merged Analysis of Precipitation)由NOAA地球系統(tǒng)研究實驗室(Earth System Research Laboratory)提供,主要基于5顆衛(wèi)星(即GPI,OPI,SSM/I scattering, SSM/I emission and MSU)觀測提供的實測數(shù)據(jù)反演所獲得的月平均降水量數(shù)據(jù),水平分辨率為2.5(°)×2.5(°)[24]。

    (5)海表面蒸發(fā)量數(shù)據(jù)來自于伍茲霍爾海洋研的oaflux數(shù)據(jù)集,本文選取的研究范圍為西太平洋暖池海域(30°N~20°S,120°E~140°W),水平分辨率為1(°)×1(°)[25]

    (6)海表面凈熱通量數(shù)據(jù)由美國國家環(huán)境預報中心(NCEP)和美國國家大氣研究中心(NCAR)提供的NECP/NCAR再分析數(shù)據(jù),包括海表面感熱、潛熱、長波輻射、短波輻射等參數(shù),空間分辨率為2(°)×2(°)[26]。

    1.2 計算方法

    熱鹽含量水平積分暖池各標準層熱含量(OHC)和鹽含量(OSC)的面積積分為:

    OHCh=∫ρCpθ′dA;OSCh=0.001∫ρS′dA。

    (1)

    其中:ρ為海水密度;Cp為定壓比熱;S′、θ′分別代表溫度異常值,研究范圍為西太平洋暖池海域(30°N~20°S,120°E~140°W,0~120dbar)。本文將研究范圍內各個層次上每個經(jīng)、緯度網(wǎng)格點上的熱鹽含量值求和并除以相應層次上經(jīng)、緯度網(wǎng)格點數(shù)來近似計算該層次上的熱鹽含量積分均值,獲得西太平洋暖池海域各層次的熱鹽含量水平積分值,并制作各個層次上熱鹽含量水平積分均值的逐月時間序列。

    熱鹽含量的垂直積分西太平洋暖池海域上層海洋熱含量和鹽含量的垂向積分為:

    (2)

    本文將研究范圍內每個經(jīng)、緯度網(wǎng)格點上的數(shù)據(jù)沿著深度進行積分,表層熱鹽含量用10dbar層數(shù)據(jù)近似,積分下界取暖池的最大深度120dbar。

    CSEOF分解法CSEOF分解法為:

    (3)

    式中:PC為時間系數(shù);LV為空間模態(tài);i為空間模特的數(shù)目。與傳統(tǒng)EOF分解有所不同的是LVi不獨立于時間變量,而是限制在一個嵌套周期(nested period)內,即:

    LVi(t)=LVi(t+d)

    (4)

    式中d為嵌套周期(本文的d取為12個月)。該方法的特點是可以體現(xiàn)長周期變化背景下的短周期變化[27-29]。本文利用該方法探討西太平洋暖池海域熱鹽含量場的時空變化特征。

    2暖池海域熱鹽含量的基本特征

    2.1 ENSO期間表層溫、鹽度分布

    圖1(b)為ENSO期間暖池海域10dbar層上34.8等鹽線分布(其中線條顏色同圖1(a))。由圖可見,鹽度隨ENSO事件的變化沒有明顯的規(guī)律。在日界線以東海域,34.8等鹽度線基本呈緯向分布,且無論是在暖事件或冷事件中,其走向基本不變。鹽度變化較大的區(qū)域位于日界線以西,且南部變化比北部區(qū)域更為復雜,可能與該區(qū)域位于赤道輻合帶附近,有著較強的降水量以及較大的鹽含量緯向流運輸有關[30-31]。

    圖1暖池海域10 dbar層28.5℃等溫線(a)和

    34.8等鹽線(b)分布

    Fig.1Positions of the 28.5℃ isotherms (a) and

    of the 34.8 isohalines (b) at 10 dbar

    2.2 各標準層熱鹽含量變化

    (經(jīng)過12個月滑動平均;圖中黑色實線表示氣候態(tài)均值。Analyzed by moving average in 12 months; Black solid line is climatic average.)

    圖2熱含量(無量綱化)在3個標準層(10、50和120dbar)上

    的區(qū)域平均逐年變化

    Fig.2Yearly variations of area-averaged OHC

    at 10、50 and 120dbar

    各標準層鹽含量的逐年變化相對比較單一,且以短期氣候態(tài)變化為主要特征(見圖3)。與熱含量變化明顯不同的是,暖池海域表層(10dbar)鹽含量的變化幅度要遠大于次表層(50及120dbar),即由表層向下,鹽含量變化的振幅逐漸縮小。由圖3還可以看出,表層鹽含量2008年相對于多年平均值有一個大的增長,而且該增長自2009—2013年底一直維持在一個較高的水平,50dbar層鹽含量的增長要滯后于表層,且幅度要小得多,而120dbar鹽含量在氣候態(tài)均值左右波動,變化不大。

    (經(jīng)過12個月滑動平均;圖中黑色實線表示氣候態(tài)均值。Analyzed by moving average in 12 months; Black solid line is climatic average.)

    圖3鹽含量(無量綱化)在3個標準層(10、50和120dbar)上

    的區(qū)域平均逐年變化

    Fig.3Yearly variations of area-averaged OSC at

    10、50 and 120dbar

    3CSEOF分析

    3.1 熱含量

    采用(2)式計算獲得2001年1月—2014年7月間西太平洋暖池海域熱含量場,并對其進行CSEOF分解,得到前三個模態(tài)的方差貢獻率分別為46.57%、10.91%和8.26%,前兩個模態(tài)的特征向量已經(jīng)占據(jù)西太平洋暖池海域熱含量距平場的57.49%,體現(xiàn)了該海域的大部分信息,因篇幅所限,本文主要討論第一模態(tài)(見圖4)。

    圖4 熱含量CSEOF分析第一模態(tài)空間場(1—12月)

    西太平洋暖池海域熱含量變化的第一模態(tài)距平場存在著明顯的東-西反相位特征,即在0~12°N,130°E~155°E附近海域存在一舌狀變異高值區(qū),其主軸約在8°N左右,而在日界線以東海域存在另一舌狀變異高值區(qū),其主軸位置約在赤道附近,且二者的極值絕對值大小相當。這種緯向的反相位變化特點表明,當西太平洋暖池海域西部熱含量增多時,東部的熱含量則減少,反之亦然。此外,圖4給出的西太平洋暖池海域0~120dbar層海洋熱含量1~12個月變化表明,這種東-西反相位變化還具有明顯的季節(jié)變化特征。冬季(11月—次年2月)強度最大,夏季(5—8月)最小,冬季兩變異區(qū)極值絕對值增加(強度增大),變異區(qū)面積增大,且西變異區(qū)向東擴展,東部變異區(qū)向西擴展,夏季兩變異區(qū)極值絕對值減小,變異區(qū)面積也縮小,且西變異區(qū)向西收縮,東部變異區(qū)向東收縮。

    CSEOF第一模態(tài)時間序列如圖5所示,結合空間距平場分布(見圖4)可以看出,第一模態(tài)時間序列達到峰值時空間距平場表現(xiàn)為東高西低的分布特點,而達到谷值時則表現(xiàn)為西高東低的特點,說明西太平洋暖池海域0~120dbar層熱含量在2007年以后經(jīng)過了3次調整。此外,利用最大熵譜分析方法對第一模態(tài)時間序列進行分析,可知該時間序列具有3年的振蕩周期。

    圖5 熱含量CSEOF分析第一模態(tài)時間序列及其最大熵譜分析

    3.2鹽含量

    鹽含量經(jīng)CSEOF分析得到的前三個模態(tài)方差貢獻分別為33.14%,21.96%和10.43%,其中前兩個模態(tài)的特征向量之和為55.1%。圖6給出了是西太平洋暖池海域0~120dbar層鹽含量CSEOF分析第一模態(tài)分布。由圖可知,西太平洋暖池海域鹽含量由南向北基本呈“正-負-正”的三極子變化模態(tài)。其中,研究區(qū)域西北海域及南部海域為正值變異區(qū),而在18°N~8°S為負值變異區(qū),且該負值變異區(qū)范圍自西向東逐漸擴大,并向東北和東南方向延伸,幾乎跨越整個熱帶太平洋海域。此外,負值變異區(qū)出現(xiàn)兩個極值中心,一個位于赤道偏南海域,另一個位于12°N附近。由CSEOF分解得到的西太平洋暖池海域上層海洋鹽含量1~12個月的變化可見,這種三極子變化模態(tài)還具有明顯的季節(jié)性變異,其中,無論是該海域鹽含量的增加和減少,還是負值變異區(qū)2個極值中心的變化,在冬春季(1—4月)均達到最大強度,夏秋季(9—12月)強度最小。冬春季,負值變異區(qū)兩個極值絕對值達到最大,南部極值中心逐漸向東南方向延伸直到144°W,因而其極值絕對值減小,面積也縮?。欢辈康臉O值中心則向東南方向延伸,且極值中心也會逐漸向東移動,只是跨度不大,在日界線以西來回振蕩。北部正值變異區(qū)和南部正值變異區(qū)同樣在冬春季強度最大,夏秋季強度最小,北部正值變異區(qū)面積隨季節(jié)變化并不明顯,但極值中心位置在冬春季向南移動,而在夏秋季向北移動;南部正值變異區(qū)面積在冬春季向西北方向擴展,在夏秋季向東南方向收縮。

    圖6 鹽含量CSEOF分析第一模態(tài)空間場(1—12月)

    圖7 鹽含量CSEOF分析第一模態(tài)時間序列及其最大熵譜分析

    CSEOF得到的鹽含量場第一模態(tài)時間系數(shù)如圖7所示,結合空間距平場分布(見圖6)可以看出,鹽含量在第一模態(tài)時間序列的峰值時空間距平場表現(xiàn)為南北部高、中部低的分布特點,而在谷值則表現(xiàn)為南北部低、中部高的分布特點,而時間序列表明2007年以后鹽含量只經(jīng)過一次位相調整。利用最大熵譜分析方法對第一模態(tài)時間序列進行分析,可知該時間序列同樣具有3年的震蕩周期。

    4討論

    從以上分析可以看到,西太平洋暖池海域熱鹽含量都有著明顯的年際變化,其主要周期均為3年,且存在著明顯的季節(jié)變化,只不過空間場變化模態(tài)有所不同,熱含量呈東-西反相位變化,而鹽含量呈三極子變化模態(tài)。前人的研究表明,熱含量的主要變化原因包括海氣熱通量、水平平流、垂直對流等因素[35],但在年際變化上主要受到水平輻合輻散的影響[36];而引起鹽含量變化的主要因素有區(qū)域性降水、水平平流、垂直對流引起的混合等因素[37];同時,熱鹽含量都會受到ENSO事件的影響[32,12]。接下來將主要探討引起西太平洋暖池海域熱鹽含量變化的可能原因。

    4.1 影響熱含量變化的可能原因

    圖8 熱含量CSEOF1時間序列(藍色)和指數(shù)(紅色)

    圖9 緯向風異常與熱含量CSEOF第一模態(tài)時間序列相關場

    (上圖為緯向流輸送量CSEOF第一模態(tài)1—12月均值;下圖藍色實線為其時間序列,紅色實線為Nino3.4指數(shù)。Monthly average spaual distributional patter of CSEOF mode 1 is in the upper picture; In the bottom picture, the blue solid line is for time series of it and the red line is for Nino3.4 Index..)

    圖10緯向流輸送量CSEOF第一模態(tài)1-12月平均態(tài)及其時間序列

    Fig.10Monthly average spatial distributional pattern of CSEOF mode 1 for OHC and Time series of CSEOF1 for OHC (blue)

    4.2 影響鹽含量變化的可能原因

    目前國內外對大尺度海盆鹽含量變化的理論機制研究還十分有限,但已有的針對中小尺度海盆鹽含量的研究發(fā)現(xiàn),海流和海表面淡水通量是影響鹽含量變化的主要因素[37]。本文給出的鹽含量場是通過深度積分計算獲得的,因此這里不討論垂向流的影響。前人研究[42]指出,南太平洋輻合帶海表面鹽度變化最大(本文分析同樣顯示該區(qū)域是鹽含量的主要變異區(qū)),且在這片海域內海表面淡水通量主要依賴于降水量的影響,故作者將主要從緯向流和降水兩個方面來討論引起鹽含量變化的可能機制,同時探討鹽含量與ENSO循環(huán)之間可能的聯(lián)系。

    進一步分析西太平洋暖池海域淡水通量(E-P)的年際變異情況(見圖11)。淡水通量EOF第一模態(tài)時間系數(shù)顯示,淡水通量超前于鹽含量6個月二者的時間系數(shù)相關最大,相關系數(shù)可達-0.71。由圖11可知,淡水通量最大變異區(qū)位于西太平洋暖池海域中部,并呈現(xiàn)由北向南的負正負三級模式,其空間分布特征與鹽含量相似,其中淡水通量呈現(xiàn)正異常時代表海表面失去淡水,負異常則相反。結合淡水通量第一模態(tài)空間及時間序列可知,當研究海域中部失去淡水,西北部和南部獲得淡水時,而鹽含量滯后于淡水通量6個月在研究海域中部出現(xiàn)升高,在西北部和南部出現(xiàn)降低。因此可以推測,海表面淡水通量可能對鹽含量有著正反饋作用。

    (紅線為E-P第一模態(tài)時間函數(shù);藍線為超前鹽含量第一模態(tài)時間系數(shù)6個月。The red line is time series of EOF1 for E-P; The blue line is six months ahead of time series of CSEOF1 for OCS.)

    圖11海表面淡水通量(E-P)EOF第一模態(tài)空間分布及時間函數(shù)

    Fig.11EOF1 pattern and timeseries of the Fresh water flux

    5結論

    本文基于Argo資料,利用CSEOF分析、譜分析等方法研究了西太平洋暖池海域上層海洋熱鹽含量的時空變化特征,并探討了熱鹽含量的可能變化機制。得到結果如下:

    (1)暖池海域次表層熱含量的變化要遠遠大于表層,而鹽含量恰好相反,以表層變化為主,次表層變化不大,且熱含量受ENSO事件影響較大,而鹽含量受ENSO時間的影響并不大,以氣候態(tài)變化為主。

    (2)西太平洋暖池海域上層海洋熱鹽含量空間場有著明顯的季節(jié)和年際變化,其年際變化主要周期都為3年;熱含量空間場呈現(xiàn)東-西反相位振蕩特征,鹽含量空間場呈現(xiàn)正-負-正三極子變化模式,且時間序列顯示,熱含量在2007年以后經(jīng)過3次位相調整,且這種年際變化都與ENSO事件有關,而鹽含量2007年以后只經(jīng)過一次位相調整。

    (4)西太平洋暖池海域鹽含量變化受到緯向流和淡水通量異常正反饋式的影響。

    致謝:感謝許建平老師對本文的指導和幫助,感謝兩位審稿人提出的修改意見,極大地幫助了文章的完善,感謝中國海洋大學和國家海洋局第二海洋研究所對筆者的培養(yǎng)。本文所使用的Argo資料來自于日本海洋科學與技術廳,是本文研究內容不可或缺的數(shù)據(jù)基礎,在此感謝。

    參考文獻:

    [1]林傳蘭. 1964-1982年熱帶西北太平洋海洋上層熱含量的變化特征[J]. 熱帶海洋, 1990, 9(2): 78-85.

    Lin Chuanlan. Some features of heat content changes of the oceanic upper layer in north west pacific during 1964-1982[J]. Tropic Oceanology, 1990, 9(2): 78-85.

    [2]于淑秋, 林學椿. 北太平洋海溫的氣候躍度及其對中國汛期降水的影響[J]. 熱帶氣象學報, 1997, 13(3): 265-275.

    Yu Shuqiu, Lin Xuechun. Climatic jump of north pacific SST and its effect on precipitation of floods season in China[J]. Journal of tropical Meteorology, 1997, 13(3): 265-275.

    [3]Suarez M J, Schopf P S. A delayed action oscillator for ENSO[J]. Journal of Atmosphere Science, 1988, 45(21): 3283-3287.

    [5]Picaut J, Masia F, Penhoat Y. An advective reflective conceptual model for the oscillatory nature of the ENSO[J]. Science, 1997, 277(1): 663-666.

    [6]Jin F F. An equatorial ocean recharge paradigm for ENSO. Part I: Conceptualmodel[J]. Journal of Atmosphere Science, 1997, 54(7): 811-829.

    [7]Jin F F. An equatorial ocean recharge paradigm for ENSO. Part II: A stripped-down coupled model[J]. Journal of Atmosphere. Science, 1997, 54(7): 830-847.

    [9]巢紀平. 對“厄爾尼諾”、“拉尼娜”發(fā)展的新認識[J]. 中國科學院院刊, 2001, 6: 412-417.

    [10]張人禾, 巢紀平. 對ENSO循環(huán)機理的一些新認識[J]. 氣候與環(huán)境研究, 2002, 7(2): 175-183.

    Zhang Renhe, Chao Jiping. Some new aspects in understanding of ENSO cycle[J]. Climatic and Environmental Research, 2002, 7(2): 175-183.

    [11]齊慶華, 侯一筠, 張啟龍. 赤道太平洋緯向風和流異常與西太平洋暖池緯向運移[J]. 海洋與湖沼, 2010, 41(3): 469-476.

    Qi Qinghua, Hou Yijun, Zhang Qilong. Zonal wind stress and current anomalies in equatorial pacific ocean and the zonal displacement of western pacific warm pool[J]. Oceanologia et Limnologia Sinica, 2010, 41(3): 469-476.

    [12]謝強, 李海洋, 王東曉. 熱帶太平洋鹽含量的年際變化[J]. 海洋科學進展, 2009, 27(2): 155-165.

    Xie Qiang, Li Haiyang, Wang Dongxiao. Interannual variations of salt content in the tropical pacific ocean[J]. Advances in Marine Science, 2009, 27(2): 155-165.

    [13]Broecker W S. The great ocean conveyor [J]. Oceanography, 1991, 4(2): 79-89.

    [14]Harvey H W. Recent Advances in the Chemistry and Biology of Sea Water [M]. Britain: Cambridge University Press, 1945.

    [15]Palmer M D, Hanines K. Estimating oceanic heat content change using isotherms[J]. Journal of Climate, 2009, 22(19): 4953-4969.

    [16]Lyman J M, Johnson G C. Estimating annual global upper-ocean heat content anomalies despite irregular in situ ocean sampling[J]. Climate, 2008, 21(21): 5629-5641.

    [17]Abraham J, Baringer M, Bindoff N et al. A review of global ocean temperature observations: Implications for ocean heat content estimates and climate change[J]. Reviews of Geophysics, 2013, 51(3): 450-483.

    [18]Levitus S. Annual cycle of temperature and heat storage in the World Ocean[J]. Journal of Physical Oceanography, 1986, 16(2): 322-343.

    [19]吳曉芬, 許建平, 張啟龍. 基于Argo資料的熱帶西太平洋上層熱含量初步研究 [J]. 海洋預報, 2001, 28(4): 76-86.

    Wu Xiaofen, Xu Jianping, Zhang Qilong et al. A preliminary study on upper ocean heat content of tropical Western Pacific[J]. Marine Forecasts, 2001, 28(4): 76-86.

    [20]Lyman J M, Willis J K, Johnson G C. Recent cooling of the upper ocean[J]. Geophysical Research Letters, 2006, 33: L18604.

    [21]許建平, 劉增宏, 孫朝輝. 全球Argo實時海洋觀測網(wǎng)全面建成[J]. 海洋技術, 2008, 27(1): 68-70.

    Xu Jianping, Liu Zenghong Sun Chaohui. Entire achievement of global argo real-time ocean observing network[J]. Ocean Technology, 2008, 27(1): 68-70.

    [22]Shigeki Hosoda, Tsuyoshi Ohira, Tomoaki Nakamura. A monthly mean dataset of global oceanic temperature and salinity derived from Argo float observations[J]. JAMSTEC Rep Res Dev, 2008, 8: 47-59.

    [23]Wallcraft A. Hybrid Coordinate Ocean Model, User's Guide[M]. Miami: University of Miami Press, 2003: 76.

    [24]Xie P P, Arkin P A. Global precipitation: A 17-year monthly analysis based on gauge observations, satellite estimates, and numerical model outputs[J]. Bulletin of the American Meteorological Society, 1997, 78(11): 2539 - 2558.

    [25]Yu L, Jin X, Weller R A. Multidecade Global Flux Datasets from the Objectively Analyzed Air-sea Fluxes (OAFlux) Project: Latent and sensible heat fluxes, ocean evaporation, and related surface meteorological variables[M]. Massachusetts: Woods Hole Oceanographic Institution, 2008: 64.

    [26]Kalnay E, Kanamitsu M, Kistler R et al. The NCEP/NCAR 40-year reanalysis project[J]. Bulletin of the American Meteorological Society, 1996, 77: 437-470.

    [27]Kwang Y K, Gerald R N. EOFs of Harmonizable cyc-lostationary process[J]. Journal of the Atmospheric Sciences, 1997, 54: 2416-2427.

    [28]Kwang Y K, Wu Q G. A comparision study of EOF techniques: Analysis of nonstationary data with periodic sta-tistics[J]. Journal of the Atmospheric Sciences, 1999, 12: 185-199.

    [29]Kwang Y K, Cheng C L. Notes and correspondence on the evolution of the annual cycle in the tropical Pacific[J]. Journal of the Atmospheric Sciences, 2001, 14: 991-994.

    [30]Delcroix T, Henin C, Porte V et al. Precipitation and sea-surface salinity in the tropical Pacific[J]. Deep Sea Research, 1996, 43(7): 1123-1141.

    [31]Reverdin G, Frankignoul C, Kestenare E, et al. Seasonal variability in the surface currents of the equatorial Pacific[J]. Journal of Geophysical Research, 1994, 99(C10): 20323- 20344.

    [32]張啟龍, 翁學傳. 熱帶西太平洋暖池表層熱含量分析[J]. 高原氣象, 1999, 18(4): 584-589.

    Zhang Qilong, Weng Xuefan. Analysis of heat content of the tropical western pacific warm pool[J]. Plateau Meteorology, 1999, 18(4): 584-589.

    [33]巢紀平, 袁紹宇, 巢清塵,等. 熱帶西太平洋暖池次表層暖水的起源——對1997/1998年ENSO事件的分析[J]. 大氣科學, 2003, 27(2): 145-151.

    [34]Weihong Qian, Haoran Hu, Yafen Zhu. Thermocline oscillation and warming event in the tropical Indian Ocean[J]. Atmosphere-Ocean, 2003, 41(3): 241-258.

    [35]王凡, 胡敦欣, 穆穆, 等. 熱帶太平洋海洋環(huán)流與暖池的結構特征、變異機理和氣候效應 [J]. 地球科學進展, 2012, 27(6): 585-602.

    Wang Fan, Hu Dunxin, Mu Mu et al. Structure, variations and climatic impacts of ocean circulation and the warm pool in the tropical pacific ocean[J]. Advances in Earth Science, 2012, 27(6): 585-602.

    [36]Wu Xiaofen, Liu Zenghong, Liao Guanghong. Variation of Indo-Pacific upper ocean heat content during 2001—2012 revealed by Argo[J]. Acta Oceanologica Sinica, 2015, 34(5): 29-38.

    [38]穆明權, 李崇銀. 西太平洋暖池次表層海溫異常與ENSO循環(huán)的相互作用 [J]. 大氣科學, 2000, 24(4): 447-460.

    Mu Mingquan, Li Chongyin. Interactions between subsurface ocean temperature anomalies in the western pacific warm pool and ENSO cycle[J]. Chinese Journal of Atmospheric Sciences, 2000, 24(4): 447-460.

    [39]Cayan D R. Latent and sensible heat flux anomalies over the northern oceans: Driving the sea surface temperature[J]. Journal of Physical Oceanography, 1992, 22(8): 859-881.

    [40]Kelly K A. The relationship between oceanic heat transport and surface fluxes in the western North Pacific: 1970-2000[J]. Journal of Climate, 2004, 17(3): 573-588.

    [41]Moon J H, Song Y T. Sea level and heat content changes in the western north Pacific[J]. Journal of Geophysical Research: Oceans, 2013, 118(4): 2014-2022.

    [42]Gouriou Y, Delcroix T. Seasonal and ENSO variations of sea surface salinity and temperature in the South Pacific Convergence Zone during 1976-2000[J]. Journal of Geophysical Research: Oceans, 2002, 107(C12): SRF-12.

    責任編輯龐旻

    Preliminary Research on Upper Ocean Heat & Salt Content of the West Pacific Warm Pool

    YANG Xiao-Xin, WU Xiao-Fen, LIU Zeng-Hong

    (State Key Lab of Satellite Ocean Environment Dynamics, Second Institute of Oceanography, SOA, Hangzhou 310012, China)

    Abstract:The knowledge of heat and salt budget is crucial to our understanding of global climate change and the coupling between atmosphere and ocean. But, due to the highly disagreement of temporary and spatial distributions of history ship-based observational data in subsurface ocean, biases of different oceanographic instruments, the choice of climatology on heat and salt content anomaly calculation and so on, estimation of the upper ocean heat and salt content had much uncertainty, particularly, the seasonal and inter-annual variations of heat and salt content and their related physical mechanisms have not been researched deeply. Fortunately, after the implementation of global Argo observation program which began in the year 2000, more than 1 million profiles of temperature and salinity with high measurement accuracy and high resolution have been accumulated. Such excellence makes the precisely calculation of upper ocean heat and salt content and the analysis of their variations to be possible. So, based on gridded Argo profile data from January 2001 to July 2014, together with Cyclostationary Empirical Orthogonal Function (CSEOF), Maximum Entropy Method (MEM) and correlation analysis, we studied the spatial distribution and temporal variations of upper ocean heat and salt content (OHC & OSC) in the Western Pacific Warm Pool (WPWP), and the probable reasons of the variations were also discussed. The main results show that variations of the OHC in subsurface region are stronger than that at the surface, which are significantly influenced by ENSO events, whereas the great change of OSC occurs at the surface and displays a quasi-decadal variation. Moreover, CSEOF analysis shows that the main mode of the OHC in the WPWP region has an east-west anti-phase oscillation on inter-annual time scale. while the OSC anomaly has a "positive-negative-positive" triple mode from north to south. According to time series analysis, we found that, after 2007, the OHC had three times phase adjustments which was highly related to ENSO (El o-Southern Oscillation) events, but only once for OSC. Finally, the variations of the OHC are closely related to ENSO events, locally zonal wind and zonal currents, while the causes of OSC variation are more complicated. Among the factors, zonal currents and fresh water flux has a positive feedback on OSC change.

    Key words:heat content; salt content; CSEOF analysis; Argo data; the Western Pacific Warm Pool

    DOI:10.16441/j.cnki.hdxb.20150320

    中圖法分類號:X87

    文獻標志碼:A

    文章編號:1672-5174(2016)04-001-12

    作者簡介:楊小欣(1990-),男,碩士生,主要從事物理海洋學調查研究。E-mail:yangxiaoxin1@126.com

    收稿日期:2015-09-19;

    修訂日期:2015-10-26

    *基金項目:衛(wèi)星海洋環(huán)境動力學國家重點實驗室自主項目(SOEDZZ1522);國家自然科學基金委員會青年科學基金項目(41406022);科技基礎性工作專項項目(2012FY112300);國家自然科學基金青年科學基金項目(41206022)資助

    Supported by the Project of State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography(No.SOEDZZ1522); National Natural Science Foundation of China (41406022); Special Program for National Basic Research(2012FY112300); National Natural Science Foundation of China(41206022)

    中文在线观看免费www的网站| 一卡2卡三卡四卡精品乱码亚洲| 身体一侧抽搐| 午夜爱爱视频在线播放| 亚洲美女黄片视频| 国产精品不卡视频一区二区| 免费看av在线观看网站| 亚洲精品久久国产高清桃花| АⅤ资源中文在线天堂| 欧美人与善性xxx| 国产麻豆成人av免费视频| 国产精品人妻久久久影院| 波多野结衣高清无吗| a级一级毛片免费在线观看| 美女高潮的动态| 久久久久国内视频| 深爱激情五月婷婷| 亚洲av中文av极速乱| 国产午夜精品论理片| 国产色婷婷99| 日本黄色片子视频| 成人无遮挡网站| 亚洲无线在线观看| 久久精品久久久久久噜噜老黄 | 永久网站在线| 真实男女啪啪啪动态图| 51国产日韩欧美| 成人漫画全彩无遮挡| 在线免费观看的www视频| 1000部很黄的大片| av天堂中文字幕网| 男女那种视频在线观看| 亚洲成人久久爱视频| а√天堂www在线а√下载| 精品久久久久久久久av| 精品一区二区免费观看| 欧美日韩精品成人综合77777| 日韩人妻高清精品专区| 日本成人三级电影网站| 91在线观看av| 真人做人爱边吃奶动态| 啦啦啦观看免费观看视频高清| 国产人妻一区二区三区在| 69av精品久久久久久| 国产高清视频在线观看网站| 国产成人精品久久久久久| 精品福利观看| 国语自产精品视频在线第100页| 国产高清三级在线| 高清毛片免费看| 国产精品久久久久久久久免| 可以在线观看的亚洲视频| 久久久国产成人免费| 国产激情偷乱视频一区二区| 99久国产av精品国产电影| 女生性感内裤真人,穿戴方法视频| 久久这里只有精品中国| 日韩人妻高清精品专区| 变态另类成人亚洲欧美熟女| 午夜日韩欧美国产| 久久精品国产鲁丝片午夜精品| 午夜激情欧美在线| 久久午夜福利片| 噜噜噜噜噜久久久久久91| 久久婷婷人人爽人人干人人爱| 一区二区三区高清视频在线| 伊人久久精品亚洲午夜| 特级一级黄色大片| 中国美白少妇内射xxxbb| 亚洲欧美日韩高清在线视频| 97超级碰碰碰精品色视频在线观看| 可以在线观看的亚洲视频| 男人和女人高潮做爰伦理| 男女视频在线观看网站免费| 亚洲无线在线观看| 国产成年人精品一区二区| 一进一出好大好爽视频| 久久亚洲精品不卡| 我的女老师完整版在线观看| 国产精品一区二区三区四区久久| 噜噜噜噜噜久久久久久91| 亚洲国产欧美人成| 国产黄a三级三级三级人| 黑人高潮一二区| 一a级毛片在线观看| 深夜a级毛片| 看非洲黑人一级黄片| 99久久久亚洲精品蜜臀av| 日韩欧美 国产精品| 国产精品久久久久久久电影| 成人永久免费在线观看视频| 国产伦在线观看视频一区| 久久鲁丝午夜福利片| 欧美日韩国产亚洲二区| 国产午夜精品论理片| 国产精品99久久久久久久久| 中国美白少妇内射xxxbb| 亚洲经典国产精华液单| 久久久欧美国产精品| 啦啦啦韩国在线观看视频| 久久久久久久久久久丰满| 精品99又大又爽又粗少妇毛片| 蜜桃久久精品国产亚洲av| 欧美性感艳星| 亚洲一级一片aⅴ在线观看| 婷婷色综合大香蕉| 亚洲成人久久性| av.在线天堂| 精品熟女少妇av免费看| 永久网站在线| 91久久精品电影网| 久久久久久九九精品二区国产| 免费观看在线日韩| 精品熟女少妇av免费看| 一区二区三区高清视频在线| 日日啪夜夜撸| av国产免费在线观看| 自拍偷自拍亚洲精品老妇| 亚洲欧美日韩无卡精品| 国内少妇人妻偷人精品xxx网站| 永久网站在线| 国产乱人视频| 国产三级在线视频| 国产高潮美女av| 亚洲真实伦在线观看| 亚洲乱码一区二区免费版| 97人妻精品一区二区三区麻豆| 日韩精品有码人妻一区| 2021天堂中文幕一二区在线观| 亚洲精品亚洲一区二区| 日韩人妻高清精品专区| 天天躁夜夜躁狠狠久久av| 国产爱豆传媒在线观看| 国产午夜福利久久久久久| 在线天堂最新版资源| av.在线天堂| 在线天堂最新版资源| 国产亚洲精品久久久久久毛片| 在线观看66精品国产| 国产精品亚洲一级av第二区| av.在线天堂| av在线蜜桃| 亚洲第一区二区三区不卡| 亚洲电影在线观看av| 老熟妇仑乱视频hdxx| 大型黄色视频在线免费观看| 日本-黄色视频高清免费观看| 欧美日韩一区二区视频在线观看视频在线 | 久久国内精品自在自线图片| 免费搜索国产男女视频| 性插视频无遮挡在线免费观看| 一个人观看的视频www高清免费观看| 日本爱情动作片www.在线观看 | 又爽又黄a免费视频| 国产精品伦人一区二区| 内地一区二区视频在线| 成熟少妇高潮喷水视频| 精品久久久久久久久久免费视频| 欧美潮喷喷水| 两性午夜刺激爽爽歪歪视频在线观看| 国产探花极品一区二区| 一个人看的www免费观看视频| 日本黄色视频三级网站网址| 又粗又爽又猛毛片免费看| 成人三级黄色视频| 久久人人爽人人爽人人片va| 国产一区亚洲一区在线观看| 久久婷婷人人爽人人干人人爱| 欧美bdsm另类| 九九久久精品国产亚洲av麻豆| 国产精品一及| 国产欧美日韩精品亚洲av| 国产高清有码在线观看视频| 成人性生交大片免费视频hd| 国产v大片淫在线免费观看| 国产av麻豆久久久久久久| 人妻夜夜爽99麻豆av| 成人国产麻豆网| 亚洲乱码一区二区免费版| 少妇被粗大猛烈的视频| 日韩三级伦理在线观看| 综合色丁香网| 免费搜索国产男女视频| 日韩 亚洲 欧美在线| 亚洲av电影不卡..在线观看| 一本精品99久久精品77| 亚洲人成网站在线观看播放| 国产精品乱码一区二三区的特点| 精品福利观看| 不卡一级毛片| 精品久久久久久久久亚洲| 偷拍熟女少妇极品色| 哪里可以看免费的av片| 舔av片在线| 国产精品久久视频播放| 亚洲成人久久爱视频| 高清毛片免费看| 精品人妻一区二区三区麻豆 | 中文资源天堂在线| 日本成人三级电影网站| 国产精品久久久久久亚洲av鲁大| 亚洲av免费在线观看| 日韩欧美三级三区| 露出奶头的视频| 国产成人影院久久av| 在线看三级毛片| 国产在视频线在精品| 老女人水多毛片| 国产av麻豆久久久久久久| 午夜日韩欧美国产| 国产精品一区二区三区四区免费观看 | 日本黄大片高清| av在线观看视频网站免费| 国产精品,欧美在线| 国产一区二区激情短视频| 国产成人a∨麻豆精品| 国产真实乱freesex| 99久久成人亚洲精品观看| 观看免费一级毛片| 97超碰精品成人国产| 免费无遮挡裸体视频| 免费av不卡在线播放| 两个人的视频大全免费| 一a级毛片在线观看| 国产成人福利小说| 丝袜美腿在线中文| 亚洲av免费高清在线观看| 91精品国产九色| 国产精品野战在线观看| 欧美区成人在线视频| 国产午夜精品久久久久久一区二区三区 | 看片在线看免费视频| 亚洲,欧美,日韩| 亚洲欧美精品自产自拍| www.色视频.com| 噜噜噜噜噜久久久久久91| 国产精华一区二区三区| 特大巨黑吊av在线直播| 联通29元200g的流量卡| 日韩欧美在线乱码| 熟女电影av网| 精品国产三级普通话版| a级毛片a级免费在线| 国模一区二区三区四区视频| 国产中年淑女户外野战色| 最近最新中文字幕大全电影3| 亚洲最大成人中文| 国产乱人视频| av免费在线看不卡| 黄色视频,在线免费观看| 日本三级黄在线观看| 国产精品美女特级片免费视频播放器| 国产大屁股一区二区在线视频| 亚洲国产日韩欧美精品在线观看| 波野结衣二区三区在线| av天堂在线播放| 淫秽高清视频在线观看| 男女之事视频高清在线观看| 国产片特级美女逼逼视频| av黄色大香蕉| 大香蕉久久网| 久久久午夜欧美精品| 亚洲熟妇熟女久久| 午夜福利在线观看吧| 国产单亲对白刺激| 欧美潮喷喷水| 成人av一区二区三区在线看| 草草在线视频免费看| 中国国产av一级| 老司机影院成人| 精品久久久久久久久久久久久| 人人妻人人看人人澡| 国产精品久久久久久久电影| 国产一区二区在线观看日韩| 欧美一区二区亚洲| 一进一出抽搐动态| 美女被艹到高潮喷水动态| 亚洲av五月六月丁香网| 国产三级在线视频| 国产一区二区在线av高清观看| 午夜福利在线观看吧| 精品午夜福利视频在线观看一区| 国产伦一二天堂av在线观看| 中文字幕av成人在线电影| 国国产精品蜜臀av免费| 日本欧美国产在线视频| 久久久久久久亚洲中文字幕| 国产亚洲精品av在线| 一进一出抽搐gif免费好疼| 久久亚洲精品不卡| 一级毛片久久久久久久久女| 午夜激情福利司机影院| 国产精品野战在线观看| 我的老师免费观看完整版| 欧美精品国产亚洲| 久久午夜福利片| 亚洲不卡免费看| 看免费成人av毛片| 日本色播在线视频| 亚洲最大成人av| av卡一久久| 国产淫片久久久久久久久| 国产91av在线免费观看| 久久精品国产亚洲网站| 国内精品久久久久精免费| 国产av麻豆久久久久久久| 国产精品亚洲一级av第二区| 香蕉av资源在线| 亚洲aⅴ乱码一区二区在线播放| 欧美3d第一页| 亚洲无线观看免费| 国产精品福利在线免费观看| 91久久精品国产一区二区三区| 最近的中文字幕免费完整| 国产在视频线在精品| 伊人久久精品亚洲午夜| 午夜福利18| 一级毛片久久久久久久久女| 国产日本99.免费观看| 亚洲欧美成人精品一区二区| 日日啪夜夜撸| 欧美性猛交黑人性爽| 国内精品宾馆在线| 黄色欧美视频在线观看| 成人亚洲欧美一区二区av| 亚洲真实伦在线观看| 五月玫瑰六月丁香| 午夜福利在线观看吧| 亚洲美女搞黄在线观看 | 99热只有精品国产| 毛片女人毛片| 亚洲在线自拍视频| 成年女人看的毛片在线观看| 日韩三级伦理在线观看| 在现免费观看毛片| 久久久久久久亚洲中文字幕| 亚洲电影在线观看av| 18禁黄网站禁片免费观看直播| 日本爱情动作片www.在线观看 | 人人妻人人澡人人爽人人夜夜 | 婷婷精品国产亚洲av| 最好的美女福利视频网| 偷拍熟女少妇极品色| 人妻久久中文字幕网| 久久久久国产网址| 在线观看一区二区三区| 在线观看午夜福利视频| 久久精品国产鲁丝片午夜精品| 久久久久国产网址| 亚洲精品日韩av片在线观看| 免费大片18禁| 欧美3d第一页| 九九在线视频观看精品| 久久精品国产亚洲av香蕉五月| 少妇猛男粗大的猛烈进出视频 | 看免费成人av毛片| 人人妻人人澡欧美一区二区| 搡老妇女老女人老熟妇| 蜜桃亚洲精品一区二区三区| 日本熟妇午夜| 日本色播在线视频| 超碰av人人做人人爽久久| 人妻久久中文字幕网| 免费av不卡在线播放| 高清日韩中文字幕在线| 午夜精品在线福利| 免费一级毛片在线播放高清视频| 蜜臀久久99精品久久宅男| 欧美zozozo另类| 人妻制服诱惑在线中文字幕| 2021天堂中文幕一二区在线观| 俄罗斯特黄特色一大片| 亚洲一区高清亚洲精品| 久久久久国产精品人妻aⅴ院| 国产不卡一卡二| 男人的好看免费观看在线视频| 级片在线观看| 不卡一级毛片| 成人特级av手机在线观看| 99热这里只有是精品在线观看| 内射极品少妇av片p| avwww免费| 毛片一级片免费看久久久久| 久久久久久久亚洲中文字幕| 免费看a级黄色片| videossex国产| av黄色大香蕉| 日本黄大片高清| 我的女老师完整版在线观看| 亚洲aⅴ乱码一区二区在线播放| 成年av动漫网址| 九九在线视频观看精品| 有码 亚洲区| 午夜激情福利司机影院| 国产伦在线观看视频一区| 国产一区二区三区av在线 | 五月玫瑰六月丁香| 国产美女午夜福利| 91狼人影院| 深夜a级毛片| 99九九线精品视频在线观看视频| 亚洲av不卡在线观看| 亚洲性久久影院| 99国产极品粉嫩在线观看| 夜夜爽天天搞| 在线播放国产精品三级| 国产女主播在线喷水免费视频网站 | 亚洲性久久影院| www日本黄色视频网| 精品免费久久久久久久清纯| 国产大屁股一区二区在线视频| 免费在线观看影片大全网站| 俄罗斯特黄特色一大片| 日本三级黄在线观看| 国产乱人视频| 成人无遮挡网站| 国产探花在线观看一区二区| 久久久久久九九精品二区国产| 精品久久久久久久久av| 99热这里只有是精品在线观看| 伦理电影大哥的女人| 黑人高潮一二区| 在线看三级毛片| 国产白丝娇喘喷水9色精品| 深夜精品福利| 精品无人区乱码1区二区| 国产精品一区二区三区四区免费观看 | 久久久久久久久大av| 99久国产av精品国产电影| 国产伦精品一区二区三区四那| avwww免费| 久久午夜福利片| 午夜视频国产福利| 日本 av在线| 乱系列少妇在线播放| 国产精品不卡视频一区二区| 成人国产麻豆网| 国产欧美日韩一区二区精品| 亚洲中文日韩欧美视频| 久久欧美精品欧美久久欧美| 偷拍熟女少妇极品色| 99久久成人亚洲精品观看| 悠悠久久av| 淫秽高清视频在线观看| 99热这里只有精品一区| 熟女电影av网| 久久精品综合一区二区三区| 日本在线视频免费播放| 国产精品1区2区在线观看.| 亚洲成人精品中文字幕电影| 美女xxoo啪啪120秒动态图| 亚洲精品456在线播放app| 日本黄大片高清| 欧美一区二区亚洲| 内地一区二区视频在线| 国产女主播在线喷水免费视频网站 | 国产一级毛片七仙女欲春2| 激情 狠狠 欧美| 国产aⅴ精品一区二区三区波| av天堂在线播放| 一本久久中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 天堂影院成人在线观看| 内地一区二区视频在线| av.在线天堂| 嫩草影院新地址| 成人欧美大片| 亚洲成人中文字幕在线播放| 精品一区二区三区视频在线| .国产精品久久| 亚洲性夜色夜夜综合| 成熟少妇高潮喷水视频| 日韩精品青青久久久久久| or卡值多少钱| 99久久成人亚洲精品观看| 人人妻人人澡欧美一区二区| 毛片女人毛片| 人人妻人人澡人人爽人人夜夜 | 日韩欧美 国产精品| 两个人视频免费观看高清| 联通29元200g的流量卡| 国产 一区精品| 久久久国产成人免费| 国语自产精品视频在线第100页| 久久精品国产亚洲av香蕉五月| 最近的中文字幕免费完整| 国产熟女欧美一区二区| 国产黄色小视频在线观看| 一级毛片aaaaaa免费看小| 亚洲欧美成人综合另类久久久 | 日本-黄色视频高清免费观看| 久久精品综合一区二区三区| 国产91av在线免费观看| 天天躁夜夜躁狠狠久久av| 一本久久中文字幕| 日韩,欧美,国产一区二区三区 | 欧美中文日本在线观看视频| 在线a可以看的网站| 在线免费观看的www视频| 99热这里只有是精品在线观看| 免费在线观看影片大全网站| 久久韩国三级中文字幕| 欧美+日韩+精品| 国产精品一区二区三区四区久久| 精品99又大又爽又粗少妇毛片| 日日啪夜夜撸| 毛片一级片免费看久久久久| 赤兔流量卡办理| 99热网站在线观看| 女的被弄到高潮叫床怎么办| 日韩人妻高清精品专区| 午夜福利视频1000在线观看| 看免费成人av毛片| 日本精品一区二区三区蜜桃| 别揉我奶头 嗯啊视频| 搡老岳熟女国产| 最近的中文字幕免费完整| 夜夜看夜夜爽夜夜摸| 嫩草影院入口| 国产精品乱码一区二三区的特点| 九色成人免费人妻av| 亚洲天堂国产精品一区在线| 国产精品无大码| 国产中年淑女户外野战色| 亚洲一区二区三区色噜噜| 99久国产av精品国产电影| av天堂中文字幕网| 一边摸一边抽搐一进一小说| 91狼人影院| 九九热线精品视视频播放| 啦啦啦观看免费观看视频高清| 国产蜜桃级精品一区二区三区| 夜夜夜夜夜久久久久| 国内少妇人妻偷人精品xxx网站| 国产精华一区二区三区| 老熟妇乱子伦视频在线观看| 国产精品女同一区二区软件| 青春草视频在线免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲在线自拍视频| 一进一出好大好爽视频| 精品99又大又爽又粗少妇毛片| 成人三级黄色视频| 麻豆精品久久久久久蜜桃| 久久久午夜欧美精品| 人妻少妇偷人精品九色| 99久久久亚洲精品蜜臀av| 18禁在线无遮挡免费观看视频 | 九九热线精品视视频播放| 成人综合一区亚洲| 色吧在线观看| 国产av一区在线观看免费| 男女下面进入的视频免费午夜| 亚洲一区高清亚洲精品| 国产三级中文精品| 亚洲欧美中文字幕日韩二区| 亚洲中文字幕日韩| 欧美成人一区二区免费高清观看| 日韩精品中文字幕看吧| 一个人免费在线观看电影| 国产在线男女| 国产亚洲91精品色在线| 天堂网av新在线| 久久久久久久久久黄片| 天堂影院成人在线观看| 日韩欧美一区二区三区在线观看| 免费黄网站久久成人精品| 精品国产三级普通话版| 亚洲人成网站在线播| 精华霜和精华液先用哪个| 国产v大片淫在线免费观看| 亚洲美女黄片视频| 国产69精品久久久久777片| 搡老岳熟女国产| 中出人妻视频一区二区| 亚洲精品一区av在线观看| 中文字幕精品亚洲无线码一区| 亚洲激情五月婷婷啪啪| 麻豆国产av国片精品| 99久久久亚洲精品蜜臀av| 欧美xxxx性猛交bbbb| 神马国产精品三级电影在线观看| 国产麻豆成人av免费视频| 97碰自拍视频| 日本撒尿小便嘘嘘汇集6| 午夜福利18| 国产精品久久久久久久久免| 欧美日韩精品成人综合77777| 午夜精品国产一区二区电影 | 99riav亚洲国产免费| 你懂的网址亚洲精品在线观看 | 欧美日韩综合久久久久久| 一区二区三区免费毛片| 欧美区成人在线视频| 国模一区二区三区四区视频| 亚洲成人av在线免费| 久久天躁狠狠躁夜夜2o2o| 欧美日韩综合久久久久久| 高清毛片免费观看视频网站| 国产极品精品免费视频能看的| 亚洲欧美日韩高清专用| 国产爱豆传媒在线观看| 欧美xxxx性猛交bbbb| av中文乱码字幕在线| 成年女人毛片免费观看观看9| 日本色播在线视频| 成人二区视频| 国产三级在线视频| 如何舔出高潮| 国产伦精品一区二区三区视频9| 精品久久久久久久久久久久久| 免费电影在线观看免费观看| 天堂网av新在线| 一级毛片电影观看 | 国产精品福利在线免费观看|