• <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)

    最近中文字幕高清免费大全6| 国产精品久久久久久精品电影小说 | 欧美xxxx黑人xx丫x性爽| av播播在线观看一区| 免费少妇av软件| 国产精品嫩草影院av在线观看| 亚洲真实伦在线观看| 成人午夜精彩视频在线观看| 麻豆精品久久久久久蜜桃| 男女国产视频网站| 亚洲国产高清在线一区二区三| av线在线观看网站| 日韩av在线大香蕉| 少妇丰满av| 国产伦一二天堂av在线观看| 青春草国产在线视频| 少妇裸体淫交视频免费看高清| 午夜福利网站1000一区二区三区| 国产成人精品久久久久久| av线在线观看网站| 91精品伊人久久大香线蕉| 男女下面进入的视频免费午夜| 久久鲁丝午夜福利片| 午夜福利成人在线免费观看| 国产精品久久久久久av不卡| 中文欧美无线码| kizo精华| 秋霞在线观看毛片| 青青草视频在线视频观看| 欧美bdsm另类| 亚洲性久久影院| 国产亚洲一区二区精品| 久久久久久国产a免费观看| 久久久久久久亚洲中文字幕| 亚洲欧美精品专区久久| 亚洲一区高清亚洲精品| 22中文网久久字幕| 大香蕉97超碰在线| 亚洲在线观看片| 成人毛片a级毛片在线播放| 亚洲欧洲日产国产| 中文字幕久久专区| 国产精品伦人一区二区| 99久久中文字幕三级久久日本| 在线a可以看的网站| 高清在线视频一区二区三区| 老司机影院成人| 亚洲最大成人av| 街头女战士在线观看网站| 亚洲激情五月婷婷啪啪| 永久网站在线| 别揉我奶头 嗯啊视频| 性色avwww在线观看| 亚洲怡红院男人天堂| 女人久久www免费人成看片| 丝袜喷水一区| 精品一区二区三区视频在线| 欧美性猛交╳xxx乱大交人| 午夜日本视频在线| 中国美白少妇内射xxxbb| 午夜爱爱视频在线播放| 激情 狠狠 欧美| 亚洲成人中文字幕在线播放| 免费av不卡在线播放| 午夜久久久久精精品| 午夜激情欧美在线| 女人十人毛片免费观看3o分钟| 国产精品国产三级国产av玫瑰| 久久精品人妻少妇| 成年版毛片免费区| 国产色婷婷99| 中文字幕制服av| 精品人妻熟女av久视频| 亚洲欧洲日产国产| 自拍偷自拍亚洲精品老妇| 国产白丝娇喘喷水9色精品| 亚洲天堂国产精品一区在线| 国产精品久久久久久精品电影小说 | 国产成人91sexporn| 成年女人看的毛片在线观看| 啦啦啦韩国在线观看视频| 久久草成人影院| 欧美成人a在线观看| 激情 狠狠 欧美| 日韩不卡一区二区三区视频在线| 成人毛片a级毛片在线播放| 免费看不卡的av| 欧美极品一区二区三区四区| 亚洲av一区综合| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品国产三级普通话版| 日日撸夜夜添| 国产精品不卡视频一区二区| 国产精品麻豆人妻色哟哟久久 | 少妇高潮的动态图| 国产午夜精品一二区理论片| 国产精品蜜桃在线观看| 国产精品一二三区在线看| 亚洲精品日本国产第一区| 日韩欧美精品免费久久| 亚洲高清免费不卡视频| 欧美激情在线99| 日韩强制内射视频| 久久久久久伊人网av| eeuss影院久久| 黄片无遮挡物在线观看| 成年女人看的毛片在线观看| 国产一区二区三区综合在线观看 | 久久人人爽人人片av| 亚洲婷婷狠狠爱综合网| 青青草视频在线视频观看| 精品久久久久久久久av| 黑人高潮一二区| 在线a可以看的网站| 国产亚洲av片在线观看秒播厂 | 韩国av在线不卡| 97热精品久久久久久| av在线播放精品| 我要看日韩黄色一级片| 亚洲av免费在线观看| 日韩不卡一区二区三区视频在线| 肉色欧美久久久久久久蜜桃 | 极品少妇高潮喷水抽搐| 欧美高清性xxxxhd video| 亚洲精品456在线播放app| 日韩一区二区三区影片| 高清欧美精品videossex| 久久精品久久精品一区二区三区| 插阴视频在线观看视频| 亚洲国产欧美人成| 深夜a级毛片| 成人毛片a级毛片在线播放| 国产女主播在线喷水免费视频网站 | 国产精品99久久久久久久久| 天堂av国产一区二区熟女人妻| 久久久久久久国产电影| 国产精品嫩草影院av在线观看| 97热精品久久久久久| 丰满少妇做爰视频| 国产成人a∨麻豆精品| 一个人观看的视频www高清免费观看| 欧美日韩视频高清一区二区三区二| 国语对白做爰xxxⅹ性视频网站| 亚洲精品成人av观看孕妇| 高清毛片免费看| 成年女人在线观看亚洲视频 | 啦啦啦韩国在线观看视频| 国产精品麻豆人妻色哟哟久久 | 日本-黄色视频高清免费观看| 久久久久国产网址| 国产成人精品福利久久| 搞女人的毛片| 特级一级黄色大片| 99久久精品一区二区三区| 国产成人精品一,二区| 亚洲丝袜综合中文字幕| 小蜜桃在线观看免费完整版高清| 久久久久久久亚洲中文字幕| 久久久久久久亚洲中文字幕| 精品一区二区三卡| 最近最新中文字幕大全电影3| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久国产电影| 国产一级毛片在线| 一级黄片播放器| 能在线免费观看的黄片| 成年人午夜在线观看视频 | 久99久视频精品免费| 亚洲国产精品成人综合色| 国精品久久久久久国模美| 在线观看美女被高潮喷水网站| 嫩草影院精品99| 久久久成人免费电影| 亚洲av中文av极速乱| 亚洲欧美成人精品一区二区| 一级毛片aaaaaa免费看小| 精品不卡国产一区二区三区| 国产成人a区在线观看| 高清在线视频一区二区三区| 久久精品久久久久久久性| 人妻少妇偷人精品九色| 老司机影院成人| av卡一久久| 少妇猛男粗大的猛烈进出视频 | 草草在线视频免费看| 久久国产乱子免费精品| av播播在线观看一区| av网站免费在线观看视频 | 22中文网久久字幕| 午夜福利网站1000一区二区三区| 日本-黄色视频高清免费观看| 99视频精品全部免费 在线| 亚洲欧美日韩东京热| 日韩在线高清观看一区二区三区| 最近最新中文字幕大全电影3| 国产 一区 欧美 日韩| 色尼玛亚洲综合影院| 亚洲不卡免费看| 韩国av在线不卡| 晚上一个人看的免费电影| 青春草国产在线视频| 又爽又黄无遮挡网站| 精品一区二区三区视频在线| 国产成人a∨麻豆精品| 亚洲av一区综合| 91精品伊人久久大香线蕉| 欧美日韩国产mv在线观看视频 | 欧美成人午夜免费资源| 人妻系列 视频| 三级国产精品片| 成人二区视频| 国产亚洲一区二区精品| 成人鲁丝片一二三区免费| 日本三级黄在线观看| 久久精品熟女亚洲av麻豆精品 | 成人毛片60女人毛片免费| 看非洲黑人一级黄片| 久久97久久精品| 亚洲精品成人av观看孕妇| 亚洲人与动物交配视频| 乱系列少妇在线播放| 日韩 亚洲 欧美在线| 建设人人有责人人尽责人人享有的 | 能在线免费看毛片的网站| 能在线免费观看的黄片| 久久精品夜色国产| 啦啦啦啦在线视频资源| 精品久久久久久成人av| 亚洲天堂国产精品一区在线| 青春草视频在线免费观看| 中国美白少妇内射xxxbb| 亚洲成色77777| 亚洲欧美一区二区三区黑人 | 一个人观看的视频www高清免费观看| 美女xxoo啪啪120秒动态图| 91狼人影院| 99视频精品全部免费 在线| 夜夜看夜夜爽夜夜摸| 久久国产乱子免费精品| 日韩一本色道免费dvd| 91aial.com中文字幕在线观看| 国内揄拍国产精品人妻在线| 日本猛色少妇xxxxx猛交久久| 国产高清国产精品国产三级 | 国产在线一区二区三区精| 亚洲图色成人| 成人一区二区视频在线观看| 亚洲精华国产精华液的使用体验| 国产黄片美女视频| 欧美最新免费一区二区三区| 欧美三级亚洲精品| 22中文网久久字幕| 欧美bdsm另类| 亚洲内射少妇av| 国精品久久久久久国模美| 国模一区二区三区四区视频| 久久午夜福利片| 啦啦啦啦在线视频资源| 国产亚洲av片在线观看秒播厂 | 国产老妇女一区| 亚洲精品自拍成人| 亚洲av男天堂| 天天躁日日操中文字幕| 插阴视频在线观看视频| 国产成人aa在线观看| 亚洲av电影不卡..在线观看| 国产久久久一区二区三区| 乱码一卡2卡4卡精品| 一级毛片aaaaaa免费看小| 久久久久久久久中文| 久久久久免费精品人妻一区二区| www.av在线官网国产| 插逼视频在线观看| 午夜福利在线观看免费完整高清在| 亚洲av日韩在线播放| 亚洲成人中文字幕在线播放| 街头女战士在线观看网站| 国产 亚洲一区二区三区 | 深夜a级毛片| 国产亚洲5aaaaa淫片| 亚洲综合色惰| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产成人aa在线观看| 一级黄片播放器| 在线观看人妻少妇| 欧美人与善性xxx| 天堂中文最新版在线下载 | 成年av动漫网址| 真实男女啪啪啪动态图| 国产探花在线观看一区二区| 国产淫语在线视频| 国产探花在线观看一区二区| 又爽又黄a免费视频| 大又大粗又爽又黄少妇毛片口| 亚洲精华国产精华液的使用体验| 国产亚洲av片在线观看秒播厂 | 18禁动态无遮挡网站| 国产一区有黄有色的免费视频 | 最近中文字幕2019免费版| freevideosex欧美| 婷婷六月久久综合丁香| 男女国产视频网站| 美女xxoo啪啪120秒动态图| 亚洲va在线va天堂va国产| 97精品久久久久久久久久精品| 又爽又黄a免费视频| 日韩欧美 国产精品| 久久久成人免费电影| 国产av不卡久久| 搡老妇女老女人老熟妇| 午夜福利视频1000在线观看| 91精品国产九色| 小蜜桃在线观看免费完整版高清| 国产精品久久久久久久电影| 亚洲欧美精品专区久久| 精品午夜福利在线看| 国产视频首页在线观看| 午夜福利在线观看吧| 51国产日韩欧美| 欧美高清性xxxxhd video| 白带黄色成豆腐渣| 久久精品人妻少妇| 色吧在线观看| 一级二级三级毛片免费看| 能在线免费看毛片的网站| 最近2019中文字幕mv第一页| 欧美高清性xxxxhd video| 国产成人一区二区在线| 午夜福利网站1000一区二区三区| 精品一区二区免费观看| 久久精品久久精品一区二区三区| 三级毛片av免费| 99热这里只有是精品50| 大陆偷拍与自拍| 亚洲精品,欧美精品| 高清av免费在线| 亚洲欧美清纯卡通| 丰满少妇做爰视频| 国产精品嫩草影院av在线观看| 99久国产av精品| 91狼人影院| 亚洲av国产av综合av卡| 精品一区二区免费观看| 国产高清有码在线观看视频| 99热这里只有是精品在线观看| 美女被艹到高潮喷水动态| 亚洲av日韩在线播放| 亚洲综合色惰| 91久久精品国产一区二区成人| 国产精品久久久久久久久免| 欧美高清成人免费视频www| 国产精品综合久久久久久久免费| 国产成人91sexporn| 美女内射精品一级片tv| 国产乱人偷精品视频| 色综合色国产| 中文欧美无线码| 国产视频首页在线观看| 国产淫片久久久久久久久| 直男gayav资源| 国产一区二区三区综合在线观看 | 十八禁网站网址无遮挡 | 国产毛片a区久久久久| 久久久亚洲精品成人影院| 能在线免费观看的黄片| 99re6热这里在线精品视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久久久久大av| 超碰av人人做人人爽久久| 欧美xxxx性猛交bbbb| 亚洲三级黄色毛片| 国产女主播在线喷水免费视频网站 | 中文精品一卡2卡3卡4更新| 欧美成人一区二区免费高清观看| 精品久久久久久久久久久久久| 波野结衣二区三区在线| 亚洲伊人久久精品综合| 嫩草影院新地址| 久久草成人影院| av在线亚洲专区| 精品人妻视频免费看| 综合色丁香网| 中文字幕亚洲精品专区| 白带黄色成豆腐渣| 婷婷六月久久综合丁香| 国产精品久久久久久av不卡| 高清毛片免费看| 肉色欧美久久久久久久蜜桃 | 免费播放大片免费观看视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 又爽又黄无遮挡网站| 精品酒店卫生间| 精品久久久久久电影网| 最近手机中文字幕大全| 美女被艹到高潮喷水动态| 日韩精品青青久久久久久| 免费大片18禁| 亚洲美女搞黄在线观看| 老司机影院毛片| 爱豆传媒免费全集在线观看| 国产亚洲精品久久久com| av专区在线播放| 91aial.com中文字幕在线观看| 亚洲aⅴ乱码一区二区在线播放| 真实男女啪啪啪动态图| 一区二区三区高清视频在线| 狂野欧美激情性xxxx在线观看| 成人二区视频| 亚洲av二区三区四区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产亚洲5aaaaa淫片| 久久99热这里只频精品6学生| 国产91av在线免费观看| 一个人看的www免费观看视频| 亚洲精品乱久久久久久| 亚洲欧美一区二区三区国产| 久久精品熟女亚洲av麻豆精品 | 我的女老师完整版在线观看| 免费观看的影片在线观看| 99久久精品一区二区三区| 深夜a级毛片| 国产淫片久久久久久久久| 搞女人的毛片| 亚洲一级一片aⅴ在线观看| 午夜精品在线福利| 国产在线男女| 少妇人妻一区二区三区视频| 性插视频无遮挡在线免费观看| 丝袜喷水一区| 22中文网久久字幕| 亚洲久久久久久中文字幕| 成年免费大片在线观看| 精品少妇黑人巨大在线播放| 只有这里有精品99| av国产久精品久网站免费入址| 国产探花在线观看一区二区| 日本午夜av视频| 婷婷色综合www| 精品久久久久久久人妻蜜臀av| av免费观看日本| 亚洲va在线va天堂va国产| 黄色欧美视频在线观看| 国产黄频视频在线观看| 18禁在线无遮挡免费观看视频| 国产日韩欧美在线精品| 午夜日本视频在线| 国产免费视频播放在线视频 | 欧美日韩精品成人综合77777| 18禁裸乳无遮挡免费网站照片| 成人综合一区亚洲| 日产精品乱码卡一卡2卡三| 好男人视频免费观看在线| a级一级毛片免费在线观看| av国产免费在线观看| 亚洲最大成人手机在线| 日日摸夜夜添夜夜添av毛片| 亚洲美女视频黄频| 亚洲欧美中文字幕日韩二区| 中国美白少妇内射xxxbb| 99热网站在线观看| 国产精品一二三区在线看| 精品酒店卫生间| 亚洲内射少妇av| 欧美成人午夜免费资源| 国产永久视频网站| 嫩草影院新地址| 精品久久久噜噜| 男人舔女人下体高潮全视频| 日韩亚洲欧美综合| 欧美日本视频| 久久这里只有精品中国| 最近手机中文字幕大全| av在线亚洲专区| 国产成年人精品一区二区| 欧美日韩亚洲高清精品| www.av在线官网国产| av一本久久久久| 高清毛片免费看| 22中文网久久字幕| 国产av不卡久久| 亚洲精品日本国产第一区| 国产淫片久久久久久久久| 建设人人有责人人尽责人人享有的 | 一个人免费在线观看电影| 成年av动漫网址| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看| 免费观看av网站的网址| 最近中文字幕高清免费大全6| 日日啪夜夜爽| 国产精品无大码| 久久久久精品性色| 色综合色国产| 舔av片在线| 十八禁国产超污无遮挡网站| 精品久久久久久久久亚洲| 欧美xxxx性猛交bbbb| 国产av码专区亚洲av| 噜噜噜噜噜久久久久久91| 十八禁国产超污无遮挡网站| 看非洲黑人一级黄片| 亚洲精品第二区| 亚洲av二区三区四区| 亚洲最大成人av| 国产精品麻豆人妻色哟哟久久 | 最近手机中文字幕大全| 高清视频免费观看一区二区 | 人妻系列 视频| 麻豆久久精品国产亚洲av| 在线观看免费高清a一片| 日韩在线高清观看一区二区三区| 亚州av有码| 国产精品久久久久久久久免| 国产乱人视频| 国语对白做爰xxxⅹ性视频网站| 亚洲在久久综合| 精品一区在线观看国产| 国产色爽女视频免费观看| 亚洲成人中文字幕在线播放| 亚洲最大成人手机在线| 亚洲色图av天堂| 青春草亚洲视频在线观看| 午夜精品一区二区三区免费看| 少妇人妻一区二区三区视频| 在线免费观看不下载黄p国产| 国产一级毛片七仙女欲春2| 少妇猛男粗大的猛烈进出视频 | 国产日韩欧美在线精品| 亚洲av成人精品一二三区| 国产美女午夜福利| 岛国毛片在线播放| 国产视频首页在线观看| 免费少妇av软件| av免费在线看不卡| 亚洲精品国产av蜜桃| 国产精品福利在线免费观看| 精品人妻熟女av久视频| 大陆偷拍与自拍| 亚洲国产精品成人综合色| 国产探花极品一区二区| 国产精品精品国产色婷婷| 国产欧美日韩精品一区二区| 亚洲精品色激情综合| 免费看日本二区| 在线观看一区二区三区| 亚洲图色成人| 蜜臀久久99精品久久宅男| 免费黄色在线免费观看| 九九在线视频观看精品| 综合色丁香网| kizo精华| 男女国产视频网站| 免费观看精品视频网站| 日韩大片免费观看网站| 亚洲精品成人av观看孕妇| 亚洲精品aⅴ在线观看| 亚洲欧美日韩卡通动漫| 中国美白少妇内射xxxbb| 国产一区二区亚洲精品在线观看| 精品人妻一区二区三区麻豆| 又大又黄又爽视频免费| 91精品国产九色| 街头女战士在线观看网站| 免费在线观看成人毛片| 中文字幕av成人在线电影| 在线 av 中文字幕| 肉色欧美久久久久久久蜜桃 | 中国国产av一级| 18+在线观看网站| 日本猛色少妇xxxxx猛交久久| 欧美97在线视频| 五月玫瑰六月丁香| 99热这里只有是精品50| 亚洲av.av天堂| 亚洲图色成人| 精品一区二区免费观看| 女的被弄到高潮叫床怎么办| 日韩欧美精品v在线| 亚洲成人一二三区av| 人人妻人人澡欧美一区二区| 国产乱来视频区| 日本与韩国留学比较| 韩国高清视频一区二区三区| 国产极品天堂在线| 日韩一区二区三区影片| 只有这里有精品99| 美女高潮的动态| 亚洲欧美精品专区久久| 欧美丝袜亚洲另类| av网站免费在线观看视频 | 最近最新中文字幕大全电影3| 在线免费观看的www视频| videos熟女内射| 国产成人免费观看mmmm| 真实男女啪啪啪动态图| 国产男人的电影天堂91| 夜夜爽夜夜爽视频| 国产av不卡久久| 国产激情偷乱视频一区二区| 久久精品综合一区二区三区| a级一级毛片免费在线观看| 狂野欧美激情性xxxx在线观看| 亚洲欧美成人精品一区二区| 最近视频中文字幕2019在线8| 国产精品人妻久久久久久| 亚洲av电影在线观看一区二区三区 | av国产久精品久网站免费入址| 成人亚洲欧美一区二区av| 久久久久九九精品影院| 99视频精品全部免费 在线| 日本一二三区视频观看| 身体一侧抽搐| 人妻少妇偷人精品九色| 亚洲精品视频女| 免费看av在线观看网站|