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

    海洋熱含量對太陽總輻射11年周期變化的響應(yīng)

    2014-03-02 05:25:06王剛顏雙喜林敏
    關(guān)鍵詞:太陽活動太陽輻射高值

    王剛顏雙喜林敏

    (1 國家海洋局數(shù)據(jù)分析與應(yīng)用重點實驗室,國家海洋局第一海洋研究所, 青島 266061;2 國家海洋局海洋環(huán)境科學(xué)與數(shù)值模擬重點實驗室,國家海洋局第一海洋研究所,青島 266061;3 中國海洋大學(xué)數(shù)學(xué)科學(xué)學(xué)院數(shù)學(xué)系,青島 266100)

    海洋熱含量對太陽總輻射11年周期變化的響應(yīng)

    王剛1,2顏雙喜3林敏3

    (1 國家海洋局數(shù)據(jù)分析與應(yīng)用重點實驗室,國家海洋局第一海洋研究所, 青島 266061;2 國家海洋局海洋環(huán)境科學(xué)與數(shù)值模擬重點實驗室,國家海洋局第一海洋研究所,青島 266061;3 中國海洋大學(xué)數(shù)學(xué)科學(xué)學(xué)院數(shù)學(xué)系,青島 266100)

    地球接收到的太陽總輻射(TSI)存在一個約11年的微小周期變化,該變化可能對地球氣候系統(tǒng)造成影響。眾多研究表明,局部海域海表面溫度、海洋熱含量的平均時間序列與TSI的11年周期變化有顯著的相關(guān)性。但海洋對太陽輻射11年周期的響應(yīng)并非空間均勻的,響應(yīng)機制也不確定。這里采用合成平均差方法,根據(jù)太陽活動的高值年和低值年,對美國國家海洋和大氣管理局(NOAA)和日本氣象廳(JMA)的全球上700m熱含量進行合成分析。兩套數(shù)據(jù)產(chǎn)品得到的海洋熱含量對太陽輻射的響應(yīng)空間模態(tài)基本一致。在太平洋分別選取其中一個正響應(yīng)區(qū)域和一個負響應(yīng)區(qū)域,發(fā)現(xiàn)這兩個區(qū)域的上700m熱含量確實存在較明顯的11年周期變化。

    太陽總輻射,海表面溫度,海洋熱含量

    1 引言

    太陽輻射是地球氣候系統(tǒng)的主要驅(qū)動力。太陽總輻射(TSI)指位于地球大氣層之外一個標準地—日距離處的單位圓盤所接收到的太陽輻射通量。衛(wèi)星觀測時代開始之前,TSI的測量在地面進行,科學(xué)家普遍認為它是一個常量(約1367W·m-2),因而也稱其為“太陽常數(shù)”[1-4]。只有少數(shù)科學(xué)家根據(jù)一些間接資料推測,TSI可能隨太陽黑子等太陽活動而發(fā)生變化[5]。有了衛(wèi)星觀測以后,TSI的觀測精確度大大提升,TSI與太陽活動的密切關(guān)系也被逐漸證實。例如,TSI隨太陽黑子數(shù)有一個11年左右的周期變化[6-7],太陽活動極大年(峰年)與太陽活動極小年(谷年)之間的TSI相差約1W·m-2。通常所說的太陽活動周(或太陽黑子周期)就是指太陽活動(太陽輻射和高能粒子等)和表象(黑子數(shù)、耀斑等)的平均11年左右(通常在9~13.6年之間)的周期變化。此外,在年代際尺度上,TSI還存在22年的Hale周期和80~90年的Gleissberg周期[8]。根據(jù)TSI與太陽活動的關(guān)系,天文學(xué)家可以重構(gòu)上百年的TSI時間序列[9-10]。

    Kopp等[11]根據(jù)衛(wèi)星觀測數(shù)據(jù)將目前TSI的平均值修正到約1361W·m-2。

    早在200多年前,Herschel[12]就推測,太陽活動可能對地球氣候系統(tǒng)產(chǎn)生影響。作為太陽活動的一個主要指標,TSI在年代至千年時間尺度上的變化都可能會影響到氣候系統(tǒng)[13-14]。在有衛(wèi)星觀測之前,Eddy[5]就提出,TSI隨黑子數(shù)變化的準10年周期可能對地球氣候系統(tǒng)產(chǎn)生影響。Beer等[15]采用一個非線性回歸模型估計,1850—1990年的140年中,全球變暖約40%源于太陽輻射。Scafetta等[16]以太陽—氣候耦合模型估算,太陽對1900—2000年全球變暖的貢獻率達到45%~50%,對1980—2000年全球變暖貢獻率約為25%~35%。Gleiser等[17]在低緯度和中緯度對流層的速度、熱力變量中都發(fā)現(xiàn)了顯著的太陽11年周期。Friis-Christensen等[18]發(fā)現(xiàn)太陽黑子周期長度的變化與北半球氣候距平時間序列有很高相關(guān)性,因而可以作為太陽活動與氣候系統(tǒng)的一個相關(guān)指數(shù)。

    海洋是地球氣候系統(tǒng)的重要組成部分。本文主要探討海洋熱含量對TSI的11年周期的響應(yīng)。第二部分,將回顧以往研究中給出的TSI的11年周期變化影響上層海洋溫度的證據(jù)及可能的物理機制;第三部分,給出兩個海洋熱含量具有明顯11年周期的海區(qū),這兩個海區(qū)可能是對TSI的11年周期具有顯著響應(yīng)的區(qū)域;第四部分是總結(jié)和展望。

    2 海洋溫度對太陽11年周期的響應(yīng)

    Reid[19]發(fā)現(xiàn),1985年之前130年的全球平均海表面溫度(SST)異常與太陽黑子數(shù)的變化,在11年滑動平均之后具有非常一致的變化趨勢。根據(jù)1980—1985年間TSI的變率以及130年的全球平均SST數(shù)據(jù)推測,SST可能與TSI的11年周期變化的包絡(luò)(即80~90年的Gleissberg周期)具有相同的相位[20]。Reid建立的海洋熱結(jié)構(gòu)模式計算表明,驅(qū)動海洋變化需要TSI約0.1%~1%的變化量。從而,TSI在一個太陽黑子周期內(nèi)的變化幅度(約0.1%)可能會引起氣候系統(tǒng)的變化[20]。不過,也有研究認為其還不足以對氣候造成顯著影響[15]。White等[21]計算了1955—1991年共38年的全球平均SST異常的功率譜,發(fā)現(xiàn)各海盆和全球平均SST時間序列在9~13年以及18~25年的時間尺度上對TSI的變化有顯著響應(yīng)。1955—1996年間,全球海洋上層熱含量也有8~15年及15~30年周期的顯著變化,并與TSI的變化同步[22]。曲維政等[23]對1955—1999年間北太平洋中部400m深海洋溫度的異常進行了功率譜分析,發(fā)現(xiàn)北太平洋中部和南太平洋中部海水溫度具有明顯11年周期。

    上層海洋對太陽11年周期活動的響應(yīng)在空間上也是分布不均的。通常采用合成平均差(Composite Mean Difference)方法,即根據(jù)上層海洋溫度在太陽活動峰年(或高值年)和谷年(或低值年)的平均差,來確定海洋對太陽響應(yīng)的空間結(jié)構(gòu)。眾多研究表明,太平洋SST中發(fā)現(xiàn)的太陽11年周期信號類似于ENSO信號。van Loon等[24-25]給出太平洋SST對TSI的11年周期峰年的響應(yīng),發(fā)現(xiàn)熱帶北太平洋區(qū)域?qū)SI的11年周期的響應(yīng)為La Ni?a型信號,即在太陽活動峰年,熱帶西、中太平洋SST異常偏暖,而熱帶東太平洋異常偏冷。Meehl等[26-27]在熱帶太平洋的SST時間序列中發(fā)現(xiàn)了類似的強La Ni?a型冷信號。然而White等[21-22]、Tung等[28]以及Roy等[29]則在同一區(qū)域發(fā)現(xiàn)了弱El Ni?o型的暖信號。不同學(xué)者得到的海洋對TSI的11年周期響應(yīng)的空間分布在相位上的差異,可能是他們采用不同的數(shù)據(jù)濾波方法造成的[25,30];也有可能太陽對海洋的影響并不顯著,因而被ENSO信號所掩蓋[31]。

    在地球表面監(jiān)測TSI的11年周期變化對氣候系統(tǒng)的影響并不容易,這一方面是由于TSI本身的變化非常小,另一方面是由于近些年來溫室氣體、火山、氣溶膠等的變化對氣候系統(tǒng)的影響更加顯著,很難從中分離出太陽的影響。太陽活動影響地球氣候系統(tǒng)的觀測證據(jù),似乎都只能從時間序列的統(tǒng)計特征來證明它們的相關(guān)性,因而受到一些學(xué)者的質(zhì)疑[32]。此外,從能量角度來說,TSI的11年周期變化只能使全球平均SST產(chǎn)生約0.1K的變化,這與觀測(至少是局部觀測)所顯示的太陽11年尺度的氣候變化的量值相差較大。一種可能的解釋是,太陽與地球氣候系統(tǒng)存在一個復(fù)雜的相互作用機制,使得地球表面只有局部地區(qū)對TSI的11年周期變化有顯著響應(yīng)[33]。

    學(xué)者們所提出的太陽對氣候系統(tǒng)的影響機制大致可以分為兩類:一是太陽通過高能粒子的調(diào)制作用對氣候系統(tǒng)產(chǎn)生影響[34];二是TSI的變化直接作用于氣候系統(tǒng)[35]。無論是高能粒子還是TSI本身,它們的變化對氣候系統(tǒng)造成的影響都非常小,但是其長期累積效應(yīng)或者非線性效應(yīng)導(dǎo)致的放大反饋過程可能會對氣候系統(tǒng)起到明顯作用。

    太陽通過高能粒子的調(diào)制作用對氣候系統(tǒng)的影響主要是在高緯大氣中引發(fā)化學(xué)反應(yīng),進而作用到氣候系統(tǒng)。目前尚未有高能粒子直接影響海洋溫度的研究和報道。TSI直接作用機制包括由底至上機制(海洋吸收太陽輻射,然后傳遞到大氣中)和自頂向下機制(大氣對太陽輻射的響應(yīng)驅(qū)動了海洋的響應(yīng))。由底

    至上機制主要強調(diào)海洋的熱吸收能力。海洋熱容量巨大,因而其直接吸收TSI造成的長期累積效應(yīng)可能會比較顯著[36]。Meehl等[26]提出由底至上機制解釋了熱帶太平洋SST對TSI的11年周期響應(yīng)的La Ni?a型空間結(jié)構(gòu)形成過程(圖1):太陽活動峰年時,熱帶太平洋東部無云區(qū)海表面接收到的太陽輻射增加(圖1a),從而增加了海洋的感熱通量;蒸發(fā)導(dǎo)致大氣中的水汽增多,更多的水汽隨著信風(fēng)被輸送到中西太平洋多云區(qū),加強了降水量和向?qū)μ栞椛涞姆瓷洌▓D1b);東太平洋無云區(qū)大氣沉降的增強,使得云量進一步減少,然后導(dǎo)致更強的信風(fēng)和更強的赤道輻合帶(圖1c)。海氣相互作用過程使得太陽輻射在太平洋東部增加而西部較少。自頂向下機制則主要來解釋太陽輻射的紫外部分如何通過同溫層自上而下影響氣候系統(tǒng)。White[35]對海洋SST和熱含量等的熱收支的診斷認為,SST中的11年周期信號不能僅僅以由底至上機制來解釋,TSI中的紫外輻射部分通過大氣同溫層自頂向下影響氣候系統(tǒng)的作用也可能存在。Misios等[37]的數(shù)值模擬實驗也說明,大氣對TSI的11年周期的響應(yīng)驅(qū)動了熱帶太平洋中該尺度的信號,并通過?!獨夥答佭^程放大。熱帶和副熱帶海區(qū)可能同時存在由底至上和自頂向下兩種響應(yīng)機制[30]。

    圖1 熱帶太平洋SST對太陽輻射11年周期響應(yīng)的機制(a)全球平均太陽輻射增加0.2W·m-2即可使東太平洋無云區(qū)增加約2W·m-2的熱通量;(b)無云區(qū)海面溫度的增加導(dǎo)致海洋潛熱通量增加,大氣中的水汽增多。信風(fēng)將更多的水汽帶到西太平洋的赤道輻合帶,從而輻合帶云量增加,太陽輻射反射增強且降水增多,海表面溫度反而降低;(c)隨著東部大氣沉降的增強,無云區(qū)云量進一步減少,海面接收到的太陽輻射進一步增強,導(dǎo)致更強的信風(fēng)和赤道輻合帶(據(jù)Meehl等[26]改編)

    3 海洋熱含量對TSI的11年周期的響應(yīng)

    目前關(guān)于海洋對太陽11年周期的響應(yīng)研究,大致有如下3種方法:一是根據(jù)海洋觀測數(shù)據(jù)計算局部海區(qū)(海盆或更大尺度)平均溫度(SST或海洋熱含量等)的時間序列和太陽活動(TSI或黑子數(shù))時間序列的相關(guān)性;二是通過對海洋觀測數(shù)據(jù)的合成平均或合成平均差來確定海洋對太陽11年周期響應(yīng)的空間分布;三是采用?!獨怦詈夏J剑谀J街懈淖兣c太陽活動有關(guān)的參數(shù)來做敏感性試驗,從而推測海洋對太陽11年周期的響應(yīng)過程和機制。

    本研究探討了海洋熱含量(如果沒有特別說明,海洋熱含量就是指海洋熱含量異常)與TSI的11年周期的關(guān)系。TSI數(shù)據(jù)基于Wang等[38]的通量傳輸模型模擬得到(http://lasp.colorado.edu/home/sorce/data/tsi-data/),然后進行平移,以便和Kopp等[11]校正后的TSI觀測值銜接。海洋熱含量采用源自美國國家海洋和大氣管理局(NOAA)和日本氣象廳(JMA)的共兩套數(shù)據(jù)。NOAA數(shù)據(jù)由Levitus等[39]基于WOD09構(gòu)建,該數(shù)據(jù)集提供全球1°×1°格點海洋熱含量數(shù)據(jù),時間從1955—2011年,包括季節(jié)平均、年平均和5年平均,垂向包括0~700m和0~2000m兩個水深。本文采用的是0~700m層的5年平均數(shù)據(jù)。JMA數(shù)據(jù)由Ishii等[40]基于WOD05構(gòu)建,空間分辨率也是1°×1°,時間從1950—2011年,該數(shù)據(jù)為年平均數(shù)據(jù)。為與NOAA數(shù)據(jù)保持一致,我們對JMA數(shù)據(jù)進行了5年滑動平均。

    首先,根據(jù)TSI的量值(圖2)將1955—2012年這段時間分為太陽活動高值年和太陽活動低值年兩類。前者為包括太陽活動峰年在內(nèi)的TSI的11年周期中輻射值較高的年份(大約5年),后者為包括太陽活動谷年在內(nèi)的TSI的11年周期中輻射值較低的年份。我們注意到,如果以一個常數(shù)(如1955—2012年的TSI平均值,圖2藍線)來劃分高值年和低值年,可能會造成一個太陽黑子周期內(nèi)有多個高值年時段或多個低值年時段的情況(如1965—1976年間,1971年為低值

    年,而1970年和1972年為高值年)。為此,本文采用經(jīng)驗正交分解(EMD)方法,提取TSI中的11年周期的本征模態(tài)函數(shù),以此模態(tài)代替TSI來確定高值年和低值年:高于平均TSI的EMD分量所對應(yīng)的年份為太陽活動高值年,其余年份為太陽活動低值年。在做EMD分解時,為了減小EMD分解中邊界效應(yīng)所引起的誤差,選取的時間段從1947年開始。在做合成平均差時,實際時段根據(jù)海洋熱含量的時段截取。

    太陽活動高值年和低值年確定之后,我們采用合成平均差(將海洋中每個點的0~700m熱含量按照太陽活動高值年和低值年分為兩組,然后計算兩組的平均值之差)來確定海洋熱含量對太陽活動的響應(yīng)。圖3是兩套熱含量數(shù)據(jù)得到的熱含量對太陽11年周期響應(yīng)的空間模態(tài)。

    圖2 年平均TSI(黑線)以及它的11年周期的EMD分量(紅線)。圖中高于1955—2012年的TSI平均值(藍線)的EMD分量所對應(yīng)的年份為太陽活動高值年,其余為太陽活動低值年

    圖3 合成平均差得到的0~700m全球海洋熱含量對太陽11年周期的響應(yīng)模態(tài)。(a)和(b)分別是NOAA和JMA數(shù)據(jù)得到的結(jié)果。如果某個點的海洋熱含量具有與TSI同相位的11年周期振蕩,則該點表現(xiàn)為暖色調(diào);反之,如果某點的海洋熱含量具有同TSI反相位的11年周期振蕩,則該點在表現(xiàn)為冷色調(diào)。矩形區(qū)域A和B對應(yīng)著太平洋海洋熱含量對TSI正響應(yīng)顯著和負響應(yīng)顯著的兩個區(qū)域

    很明顯,TSI的11年周期只是圖3中空間各點熱含量對太陽輻射響應(yīng)為正或負的充分條件,而不是必要條件。換言之,圖3中為正響應(yīng)的點,其熱含量未必具有與TSI同相位的11年周期振蕩;同樣,負響應(yīng)的點也未必具有與TSI反相位的11年周期振蕩。為此,我們在熱帶太平洋選擇兩個海洋熱含量對TSI響應(yīng)比較顯著的海區(qū),來進一步觀察其是否具有11年周期振蕩。圖3中A區(qū)域和B區(qū)域分別對應(yīng)太平洋海洋熱含量對TSI正響應(yīng)或負響應(yīng)比較顯著的兩個區(qū)域。圖4和圖5分別給出這兩個區(qū)域經(jīng)向(子午向)平均海洋熱含量隨時間的變化。在NOAA數(shù)據(jù)(圖4a)和JMA數(shù)據(jù)(圖4b)中,A區(qū)域的海洋熱含量都表現(xiàn)出明顯的準10年振蕩。B區(qū)域的海洋熱含量在1980年代之后表現(xiàn)出明顯的準10年周期振蕩。如果認為1980年左右存在一個冷相位(圖5a中比較明顯,圖5b中則不明顯),則從1950年代至2011年會出現(xiàn)一個連續(xù)的10年周期振蕩。這個冷相位的弱化或缺失,可能與1975年左右熱帶太平洋的穩(wěn)態(tài)轉(zhuǎn)換(regime shift)[41]有關(guān)。以Levitus數(shù)據(jù)為例,我們對這兩個區(qū)域進行了小波分析。功率譜(圖6)顯示,兩個區(qū)域的平均熱含量都在11年左右周期上有明顯譜峰,并通過了紅噪聲95%的置信度檢驗。對于A區(qū)域(圖6a),11年周期為最顯著周期;對于B區(qū)域(圖6b),25年左右周期最顯著,11年周期次之。

    圖4 赤道中太平洋(圖3中A區(qū)域)海洋熱含量隨時間的變化。海洋熱含量具有明顯的11年周期振蕩,且與太陽11年周期信號同相位。(a)和(b)分別是NOAA和JMA數(shù)據(jù)得到的結(jié)果

    圖5 西太平洋暖池(圖3中B區(qū)域)海洋熱含量隨時間的變化。海洋熱含量具有明顯的11年周期振蕩,且與太陽11年周期信號反相位。(a)和(b)分別是NOAA和JMA數(shù)據(jù)得到的結(jié)果

    然后,采用一個簡單的辦法對圖3中給出的海洋熱含量對TSI的響應(yīng)模態(tài)做顯著性檢驗。1955—2012年間,TSI大致經(jīng)歷了5個黑子周期,這里將其分為9段(每段長度為一個太陽黑子周期,相鄰兩段有半個周期的重疊),在每段時間采用合成平均差做海洋熱含量的響應(yīng)投影。若某點海洋熱含量對TSI的9段響應(yīng)中同為正或同為負的出現(xiàn)7次以上,則認為該點的海洋熱含量對TSI的11年周期響應(yīng)是顯著的。結(jié)果顯示,圖3中A區(qū)域的大部分格點上對TSI的11年周期的響應(yīng)是顯著的,而B區(qū)域的顯著響應(yīng)區(qū)域集中在該區(qū)域顏色較深的位置(圖略)。這也解釋了為何圖5中B區(qū)域的11年周期不如圖4中A區(qū)域的11年周期明顯。

    此外,在黑潮、灣流以及阿加勒斯海流這幾個西邊界流的延伸區(qū),也分別看到了熱含量對太陽11年周期的顯著響應(yīng)。其中,黑潮延伸區(qū)為顯著的正響應(yīng),后兩者為負響應(yīng)。這些區(qū)域的海洋熱含量變化對局地氣候系統(tǒng)都有重要影響。

    4 總結(jié)與展望

    太陽是地球氣候系統(tǒng)的主要能量來源。太陽活動必然會影響到地球氣候系統(tǒng)。然而,太陽總輻射在一個黑子周期(11年)內(nèi)的微小變化是否能夠顯著影響到地球氣候系統(tǒng),還存在廣泛爭議。目前已有的觀測只能證明,在大氣和海洋的某些區(qū)域,存在與TSI的11年周期相關(guān)的周期信號。時間序列的相關(guān)性無法解釋TSI的微小能量變化如何引起地球氣候系統(tǒng)(局部)地區(qū)顯著的11年振蕩的問題。TSI信號通過非線性作用在氣候系統(tǒng)中的放大過程,還需要復(fù)雜的太陽—氣候系統(tǒng)的相互作用機制來解釋。目前數(shù)值模式在機制分析上起重要作用,理論證明還需要進一步發(fā)展。

    我們近期的工作從海洋熱含量對TSI響應(yīng)的空間模態(tài)中挑選了兩個響應(yīng)顯著的區(qū)域。在這兩個區(qū)域中,我們看到海洋熱含量具有明顯的準10年周期振蕩。Hasegawa等[42]曾在熱帶太平洋的海洋熱含量中發(fā)現(xiàn)了13年的振蕩信號,并認為它和ENSO事件有關(guān)。與他們發(fā)現(xiàn)的信號不同的是,我們給出的準10年信號周期更接近太陽11年周期,并且沒有發(fā)現(xiàn)明顯的東西向傳播現(xiàn)象,因而可能是不同的信號。考慮到得到該信號的過程,我們認為,它是海洋熱含量對TSI的11年周期的響應(yīng)。此外,黑潮、灣流、阿加勒斯海流這幾個對區(qū)域氣候影響顯著的西邊界流,它們延伸區(qū)的海洋熱含量對TSI的11年周期的響應(yīng)也很顯著,只是相位上有所不同。不同海區(qū)對TSI響應(yīng)的強度或相位上的差異,可能是海—氣非線性相互作用的結(jié)果,也可能是海洋本身動力過程的結(jié)果,響應(yīng)過程及作用機理有待深入探討。

    [1]Johnson F S. The solar constant. J Meteor, 1954, 11: 431-439.

    [2]Drummond A J, Hickey J R, Scholes W J, et al. New value for the solar constant of radiation. Nature, 1968, 218: 259-261.

    [3]Thekaekara M P, Drummond A J. Standard values for the solar constant and its spectral components. Nature, 1971, 229: 6-9.

    [4]Laue E G, Drummond A J. Solar constant: First direct measurements. Science, 1968, 161(3844): 888-891.

    [5]Eddy J A. The Maunder minumum. Science, 1976, 192(4245): 1189-1202.

    [6]Kuhn J R, Libbrecht K G, Dicke R H. The surface temperature of the Sun and changes in the solar constant. Science, 1988, 242(4880): 908-911.

    [7] Willson R C, Hudson H S. The Sun’s luminosity over a complete solar cycle. Nature, 1991, 351(42): 42-44.

    [8]Tsiropoula G. Signatures of solar activity variability in meteorological parameters. Journal of Atmospheric and Solar-Terrestrial Physics, 2003, 65: 469-482.

    [9] Lean J. Evolution of the Sun’s spectral irradiance since the Maunde r Minimum. Geophysical Research Letters, 2000, 27(16): 2425-2428.

    [10]Hoyt D V, Schatten K H. A discussion of plausible solar irradiance variations, 1700-1992. Journal of Geophysical Research, 1993, 98(A11): 18895-18906.

    [11]Kopp G, Lean J L. A new, lower value of total solar irradiance: Evidence and climate significance. Geophys Res Lett, 2011, 38: L01706.

    [12]Herschel W. Observations tending to investigate the nature of the Sun, in order to find the causes or symptoms of its variable emission of light and heat: With remarks on the use that may possibly be drawn from solar observations. Philos Trans R Soc London, 1801, 91: 265-318.

    [13]Njau E C. Sun-weather/climate relationships: A review—Part Ⅱ. PINSA, 2000, 66,A(5): 451-466.

    [14]Njau E C. Sun-weather/climate relationships: A review—Part I. PINSA, 2000, 66,A(3&4): 415-441.

    [15]Beer J, Mende W, Stellmacher R. The role of the sun in climate forcing. Quaternary Science Reviews, 2000, 19: 403-415.

    [16]Scafetta N, West B J. Phenomenological solar contribution to the 1900—2000 global surface warming. Geophysical Research Letters, 2006, 33: L05708.

    [17]Gleisner H, Thejll P. Patterns of tropospheric response to solar variability. Geophysical Research Letters, 2003, 30(13): 1711.

    [18]Friis-Christensen E, Lassen K. Length of the solar cycle: An indicator of solar activity closely associated with climate. Science, 1991, 254: 698-700.

    [19]Reid G C. Influence of solar variability on global sea surface temperatures. Nature, 1987, 329: 142-143.

    [20]Reid G C. Solar total irradiance variations and the global sea surface temperature record. Journal of Geophysical Research, 1991, 96(D2): 2835-2844.

    [21]White W B, Lean J, Cayan D R, et al. Response of global upper ocean temperature to changing solar irradiance. Journal of Geophysical Research, 1997, 102(C2): 3255-3266.

    [22]White W B, Cayan D R, Lean J. Global upper ocean heat storage response to radiative forcing from changing solar irradiance and increasing greenhouse gas/aerosol concentrations. Journal of Geophysical Research, 1998, 103(C10): 21355-21366. 67: 2649-2657.

    [23]曲維政, 鄧聲貴, 黃菲, 等. 深海溫度變化對太陽活動的響應(yīng). 第四紀研究, 2004, (3): 285-292.

    [24]van Loon H, Meehl G A, Shea D J. Coupled air sea response to solar forcing in the Pacific region during northern winter. J Geophys Res, 2007, 112(D2): D02108.

    [25]van Loon H, Meehl G A. The response in the Pacific to the Sun’s decadal peaks and contrasts to cold events in the Southern Oscillation. J Atmos Sol-Terr Phy, 2008, 70: 1046-1055.

    [26]Meehl G A, Arblaster J M, Branstator G, et al. A coupled air–sea response mechanism to solar forcing in the Pacific region. J Climate, 2008, 21: 2883-2897.

    [27]Meehl G A, Arblaster J M. A lagged warm event–like response to peaks in solar forcing in the Paci fi c region. J Climate, 2009, 22: 3647-3660.

    [28] Tung K K, Zhou J S. The Paci fi c’s response to surface heating in 130 yr of SST: La Ni?a-like or El Ni?o-like? J Atmos Sci, 2010,

    [29]Roy I, Haigh J D. Solar cycle signals in sea level pressure and sea surface temperature. Atmos Chem Phys, 2010, 10: 3147-3153.

    [30]Meehl G A, Arblaster J M, Matthes K, et al. Amplifying the Paci fi c climate system response to a small 11-year solar cycle forcing. Science, 2009, 325: 1114-1118.

    [31]Roy I, Haigh J D. Solar cycle signals in the Paci fi c and the issue of timings. Journal of the Atmospheric Sciences, 2012, 69: 1446-1451.

    [32]Legras B, Mestre O, Bard E, et al. A critical look at solar-climate relationships from long temperature series. Clim Past, 2010, 6: 745-758.

    [33]Christoforou P, Hameed S. Solar cycle and the Paci fi c“centers of action”. Geophys Res Lett, 1997, 24(3): 293-296.

    [34]Dickinson R E. Solar variability and the lower atmosphere. Bull Am Meteorol Soc, 1975, 56: 1240-1248.

    [35]White W B. Response of tropical global ocean temperature to the Sun’s quasi-decadal UV radiative forcing of the stratosphere. J Geophys Res, 2006, 111: C09020.

    [36]Gray L J, Beer J, Geller M, et al. Solar influences on climate. Reviews of Geophysics, 2010, 48: RG4001.

    [37]Misios S, Schmidt H. Mechanisms involved in the ampli fi cation of the 11-yr solar cycle signal in the tropical Paci fi c Ocean. Journal of Climate, 2012, 25: 5102-5118.

    [38]Wang Y M, Lean J L, Sheeley J N R. Modeling the Sun’s Magnetic Field and Irradiance since 1713. The Astrophysical Journal, 2005, 625: 522-538.

    [39]Levitus S, Antonov J I, Boyer T P, et al. World ocean heat content and thermosteric sea level change (0–2000m), 1955–2010 Geophysical Research Leters, 2012, 39: L10603.

    [40]Ishii M, Kimoto M. Reevaluation of historical ocean heat content variations with timevarying XBT and MBT depth bias corrections. J Oceanogr, 2009, 65: 287- 299.

    [41]Stephens C, Levitus S, Antonov J, et al. On the Pacific Ocean regime shift. Geophysical Research Leters, 2001, 28(19): 3721-3724.

    [42]Hasegawa T, Hanawa K. Decadal-scale variability of upper ocean heat content in the tropical Paci fi c. Geophysical Research Leters, 2003, 30(6): 1272.

    Response of Ocean Heat Content to the 11-year Variation of Total Solar Irradiance

    Wang Gang1,2, Yan Shuangxi3, Lin Min3
    (1 Key Laboratory of Data Analysis and Applications (LDAA), the First Institute of Oceanography (FIO), State Oceanic Administration (SOA), Qingdao 266061 2 Key Laboratory of Marine Science and Numerical Modeling (MASNUM), FIO, SOA, Qingdao 266061 3 Department of Mathematics, School of Mathematical Sciences, Ocean University of China, Qingdao 266100)

    The total solar irradiance (TSI) impinging on the Earth has a slight variation of approximate 11-year period, which is expected to affect the variability of the Earth climate system. Researches have found correlations between sea surface temperature or ocean heat content of local ocean basins and the 11-year TSI time series. However, the response of the upper ocean to the 11-year solar cycle is not spatially uniform, and the response mechanism is still open to debate. No region has even been identif i ed to have a clear 11-year variation in its sea surface temperature or ocean heat content yet. We analysed the ocean heat content (OHC) in upper-700-m of the global ocean from NOAA and JMA datasets, respectively. Composite mean difference method is used to obtain the response pattern of OHC to the 11-year solar activity. The two datasets give rather similar results. We select a positive response area and a negative one in the Pacif i c and fi nd a signif i cant 11-year variation in the upper 700-m OHC in these two regions.

    total solar irradiance (TSI), sea surface temperature (SST), ocean heat content (OHC)

    10.3969/j.issn.2095-1973.2014.04.003

    2013年10月16日;

    2014年3月10日

    王剛(1977— ),Email: wangg@fio.org.cn

    資助信息:國家重大科學(xué)研究計劃項目(2012CB957803);國家重點基礎(chǔ)研究發(fā)展計劃項目(2011CB403502);國家自然科學(xué)基金資助項目(41006018,41476024);山東省優(yōu)秀中青年科學(xué)家科研獎勵基金項目(BS2011HZ019)

    猜你喜歡
    太陽活動太陽輻射高值
    養(yǎng)殖廢棄物快速發(fā)酵及高值轉(zhuǎn)化土壤修復(fù)生物肥料關(guān)鍵技術(shù)
    麻文化發(fā)展與高值利用前景展望
    第24太陽活動周中國地區(qū)電離層閃爍統(tǒng)計特性研究
    第23和24太陽活動周高緯地磁感應(yīng)電流分布特性
    邯鄲太陽輻射時空分布特征
    基于PCA 的太陽輻射觀測算法研究
    太陽輻射作用下鋼筒倉結(jié)構(gòu)溫度場分析研究
    高值無害化利用 廢白土大有可為
    洛陽地區(qū)太陽輻射變化特征及影響因子分析
    河南科技(2015年11期)2015-03-11 16:25:00
    榜單
    天堂网av新在线| 国产美女午夜福利| 午夜激情欧美在线| 成人精品一区二区免费| 1000部很黄的大片| 国产精品,欧美在线| 九九热线精品视视频播放| 午夜福利视频1000在线观看| 观看免费一级毛片| 日本黄色视频三级网站网址| 久久精品国产自在天天线| 欧美一区二区国产精品久久精品| 国产成人aa在线观看| av在线老鸭窝| 国语自产精品视频在线第100页| 免费看a级黄色片| 国产精品三级大全| 美女cb高潮喷水在线观看| 国产成人啪精品午夜网站| 搡老熟女国产l中国老女人| 亚洲专区中文字幕在线| 国产欧美日韩精品亚洲av| 18美女黄网站色大片免费观看| 女人被狂操c到高潮| 中文字幕高清在线视频| .国产精品久久| 九九在线视频观看精品| a级毛片a级免费在线| 国产精品久久久久久精品电影| 免费电影在线观看免费观看| 欧美成狂野欧美在线观看| 日韩欧美国产在线观看| 日本黄大片高清| 亚洲激情在线av| 天天躁日日操中文字幕| 一本精品99久久精品77| 麻豆久久精品国产亚洲av| 欧美黄色片欧美黄色片| 日韩 亚洲 欧美在线| 又黄又爽又刺激的免费视频.| 亚洲精品一卡2卡三卡4卡5卡| 最新中文字幕久久久久| 狠狠狠狠99中文字幕| 国产精品亚洲av一区麻豆| 18+在线观看网站| 亚洲精品亚洲一区二区| 18禁在线播放成人免费| 特大巨黑吊av在线直播| 欧美激情久久久久久爽电影| 国产精品一及| 精品一区二区三区视频在线| 天堂av国产一区二区熟女人妻| 两个人的视频大全免费| 哪里可以看免费的av片| 国产麻豆成人av免费视频| 亚洲av中文字字幕乱码综合| av在线观看视频网站免费| 中文字幕熟女人妻在线| 女人被狂操c到高潮| 中文资源天堂在线| 免费看a级黄色片| 成人特级av手机在线观看| 国产真实乱freesex| 免费在线观看日本一区| 亚洲精品亚洲一区二区| 亚洲av熟女| 中文字幕免费在线视频6| 在线十欧美十亚洲十日本专区| 精品人妻视频免费看| 观看美女的网站| 国产单亲对白刺激| 亚洲欧美日韩东京热| 精品久久久久久久久亚洲 | 美女xxoo啪啪120秒动态图 | 久久久久久久久大av| 欧美成人一区二区免费高清观看| 黄色丝袜av网址大全| 看十八女毛片水多多多| 久久国产乱子伦精品免费另类| 亚洲成人精品中文字幕电影| 我要搜黄色片| 亚洲av免费在线观看| 偷拍熟女少妇极品色| 亚洲av一区综合| 亚洲最大成人手机在线| 精品久久久久久久久久久久久| 免费一级毛片在线播放高清视频| 午夜两性在线视频| 草草在线视频免费看| 亚洲内射少妇av| 久久久久久久精品吃奶| 国内精品久久久久久久电影| 久久久久精品国产欧美久久久| 网址你懂的国产日韩在线| 亚洲国产精品sss在线观看| avwww免费| 少妇的逼好多水| 久久午夜福利片| 1024手机看黄色片| 日韩 亚洲 欧美在线| 国产成人欧美在线观看| 国产高清三级在线| 亚洲va日本ⅴa欧美va伊人久久| 网址你懂的国产日韩在线| 老司机深夜福利视频在线观看| 国产亚洲精品av在线| 一区二区三区激情视频| 欧美bdsm另类| 亚洲最大成人av| 男女视频在线观看网站免费| 午夜激情欧美在线| 国产精品一区二区免费欧美| 一本综合久久免费| 欧美精品啪啪一区二区三区| 国产高清视频在线播放一区| 国产精品久久久久久亚洲av鲁大| 久久精品久久久久久噜噜老黄 | 亚洲一区二区三区不卡视频| 小说图片视频综合网站| 久久欧美精品欧美久久欧美| 国产毛片a区久久久久| 超碰av人人做人人爽久久| 亚洲自偷自拍三级| 午夜两性在线视频| 亚洲精品成人久久久久久| 亚洲av.av天堂| 亚洲av成人av| 欧美日韩乱码在线| 在现免费观看毛片| 亚洲av电影不卡..在线观看| 国产精品电影一区二区三区| 午夜福利在线观看免费完整高清在 | 久久午夜亚洲精品久久| 欧美激情国产日韩精品一区| 久久久久久久亚洲中文字幕 | 亚洲熟妇熟女久久| 一个人免费在线观看电影| 亚洲精品影视一区二区三区av| 欧美黄色淫秽网站| 成人精品一区二区免费| 亚洲va日本ⅴa欧美va伊人久久| 久久午夜福利片| 国内精品久久久久久久电影| 中文字幕久久专区| 欧美日韩瑟瑟在线播放| 亚洲,欧美,日韩| 三级毛片av免费| 精品午夜福利视频在线观看一区| 亚洲天堂国产精品一区在线| 麻豆一二三区av精品| 少妇高潮的动态图| 99国产极品粉嫩在线观看| 自拍偷自拍亚洲精品老妇| 搡女人真爽免费视频火全软件 | 少妇人妻精品综合一区二区 | 一本精品99久久精品77| 国内精品一区二区在线观看| 我的老师免费观看完整版| 91av网一区二区| 国产精品三级大全| av女优亚洲男人天堂| 欧洲精品卡2卡3卡4卡5卡区| 国产成+人综合+亚洲专区| 99久久久亚洲精品蜜臀av| 高清日韩中文字幕在线| 一个人观看的视频www高清免费观看| 两人在一起打扑克的视频| 国产探花在线观看一区二区| 久久久久国产精品人妻aⅴ院| 国产麻豆成人av免费视频| 亚洲av成人精品一区久久| 国产视频一区二区在线看| 午夜a级毛片| 我要搜黄色片| 免费在线观看亚洲国产| 日韩欧美在线乱码| 久久国产精品影院| 无遮挡黄片免费观看| 蜜桃亚洲精品一区二区三区| 又爽又黄a免费视频| 天堂网av新在线| 特级一级黄色大片| 久久精品久久久久久噜噜老黄 | 国产精品1区2区在线观看.| 人人妻,人人澡人人爽秒播| 天堂动漫精品| 老鸭窝网址在线观看| 蜜桃亚洲精品一区二区三区| 国产精品不卡视频一区二区 | 一区二区三区免费毛片| 亚洲av电影不卡..在线观看| 欧美不卡视频在线免费观看| av在线观看视频网站免费| 国内毛片毛片毛片毛片毛片| 亚洲五月婷婷丁香| 综合色av麻豆| 国产三级在线视频| 亚洲av成人av| 桃色一区二区三区在线观看| 成人国产综合亚洲| 男女下面进入的视频免费午夜| 色综合站精品国产| 午夜福利在线观看吧| 精品久久久久久成人av| www.www免费av| 99国产精品一区二区蜜桃av| 18禁黄网站禁片午夜丰满| 男人狂女人下面高潮的视频| 亚洲一区二区三区不卡视频| 亚洲成人中文字幕在线播放| 国产精品久久久久久人妻精品电影| 三级毛片av免费| 午夜日韩欧美国产| 精品人妻1区二区| 一个人免费在线观看的高清视频| 国产精品久久久久久久久免 | 嫩草影视91久久| 色综合亚洲欧美另类图片| 国产黄色小视频在线观看| 久久久国产成人免费| 午夜久久久久精精品| 亚洲精品一卡2卡三卡4卡5卡| 婷婷色综合大香蕉| 日本a在线网址| 亚洲人成电影免费在线| 搡老熟女国产l中国老女人| 自拍偷自拍亚洲精品老妇| 久久久久久国产a免费观看| 亚洲欧美激情综合另类| 成人美女网站在线观看视频| 嫩草影视91久久| АⅤ资源中文在线天堂| 一本一本综合久久| 99热这里只有精品一区| 亚洲 欧美 日韩 在线 免费| 国产av不卡久久| 男人和女人高潮做爰伦理| av黄色大香蕉| 又黄又爽又刺激的免费视频.| 日本黄色视频三级网站网址| 丰满人妻一区二区三区视频av| 99久久无色码亚洲精品果冻| 免费黄网站久久成人精品 | 久久精品国产自在天天线| 制服丝袜大香蕉在线| 国产又黄又爽又无遮挡在线| h日本视频在线播放| 他把我摸到了高潮在线观看| 好男人在线观看高清免费视频| netflix在线观看网站| 精华霜和精华液先用哪个| 桃红色精品国产亚洲av| 在线国产一区二区在线| 亚洲国产精品sss在线观看| 一个人观看的视频www高清免费观看| 国产精品永久免费网站| 欧美区成人在线视频| 亚洲,欧美精品.| 99在线视频只有这里精品首页| 成人一区二区视频在线观看| 人人妻,人人澡人人爽秒播| 亚洲欧美日韩卡通动漫| 国产成人啪精品午夜网站| 身体一侧抽搐| 露出奶头的视频| 亚洲成a人片在线一区二区| 亚洲无线观看免费| avwww免费| 一本一本综合久久| 亚洲午夜理论影院| 男女做爰动态图高潮gif福利片| 俺也久久电影网| 中文字幕av在线有码专区| av天堂中文字幕网| 男女做爰动态图高潮gif福利片| 亚洲最大成人手机在线| 脱女人内裤的视频| 国产亚洲精品久久久久久毛片| 热99在线观看视频| 18美女黄网站色大片免费观看| 99热这里只有是精品50| 神马国产精品三级电影在线观看| 国产黄a三级三级三级人| 少妇人妻精品综合一区二区 | 亚洲av电影不卡..在线观看| 激情在线观看视频在线高清| 九九久久精品国产亚洲av麻豆| 青草久久国产| 国产色婷婷99| 日本免费一区二区三区高清不卡| 成年女人永久免费观看视频| 成年人黄色毛片网站| 久久人妻av系列| 三级国产精品欧美在线观看| 99热只有精品国产| 男人的好看免费观看在线视频| 亚洲在线观看片| 给我免费播放毛片高清在线观看| 国产中年淑女户外野战色| 国产精品一区二区免费欧美| 国产三级在线视频| 天天躁日日操中文字幕| 国产综合懂色| 2021天堂中文幕一二区在线观| 啪啪无遮挡十八禁网站| 淫妇啪啪啪对白视频| 男人舔奶头视频| 成人无遮挡网站| 成人精品一区二区免费| 亚洲美女黄片视频| 性色av乱码一区二区三区2| 久久国产精品人妻蜜桃| 一区二区三区高清视频在线| 九九热线精品视视频播放| 成人亚洲精品av一区二区| 亚洲第一区二区三区不卡| 欧美黄色片欧美黄色片| 狂野欧美白嫩少妇大欣赏| av福利片在线观看| 免费高清视频大片| 久久久久九九精品影院| 真人做人爱边吃奶动态| 无遮挡黄片免费观看| 18禁裸乳无遮挡免费网站照片| a级一级毛片免费在线观看| 97碰自拍视频| 欧美成人性av电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲成人久久性| 天美传媒精品一区二区| 亚洲第一电影网av| 嫩草影视91久久| 狠狠狠狠99中文字幕| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 午夜福利在线观看免费完整高清在 | av在线观看视频网站免费| 久久久久久久精品吃奶| 级片在线观看| 国产免费一级a男人的天堂| av在线天堂中文字幕| 99久国产av精品| 精品一区二区三区人妻视频| 久久6这里有精品| 欧美一级a爱片免费观看看| 午夜免费成人在线视频| 亚洲精华国产精华精| 成人欧美大片| 婷婷色综合大香蕉| 在线观看午夜福利视频| 成人美女网站在线观看视频| 一a级毛片在线观看| 国产黄色小视频在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久久久久久久大av| 性欧美人与动物交配| 免费av不卡在线播放| 最近视频中文字幕2019在线8| 久久伊人香网站| 中文亚洲av片在线观看爽| 90打野战视频偷拍视频| 国产成人啪精品午夜网站| 午夜福利高清视频| 69av精品久久久久久| 久久国产精品人妻蜜桃| 亚洲国产精品999在线| 99riav亚洲国产免费| 国产亚洲精品综合一区在线观看| 日韩中字成人| 精品人妻熟女av久视频| 少妇熟女aⅴ在线视频| 高清日韩中文字幕在线| 国产麻豆成人av免费视频| 婷婷亚洲欧美| 亚洲熟妇中文字幕五十中出| 亚洲av电影不卡..在线观看| 国产免费一级a男人的天堂| 日本免费a在线| 国产v大片淫在线免费观看| 九九在线视频观看精品| 国产精品日韩av在线免费观看| 亚洲人成网站在线播| 成人毛片a级毛片在线播放| 99久久无色码亚洲精品果冻| a级一级毛片免费在线观看| 精品国产三级普通话版| 一级作爱视频免费观看| 国产乱人伦免费视频| 欧美一区二区国产精品久久精品| 中文字幕高清在线视频| 欧美一区二区国产精品久久精品| 国产午夜精品论理片| 极品教师在线免费播放| 日韩国内少妇激情av| 亚洲av第一区精品v没综合| 熟妇人妻久久中文字幕3abv| а√天堂www在线а√下载| 亚洲最大成人中文| 天天躁日日操中文字幕| 亚洲第一区二区三区不卡| 1024手机看黄色片| 大型黄色视频在线免费观看| 国产亚洲欧美98| 亚洲自拍偷在线| 精品午夜福利视频在线观看一区| 国产主播在线观看一区二区| 真实男女啪啪啪动态图| 免费高清视频大片| 久久久精品大字幕| 色噜噜av男人的天堂激情| 国产精品久久电影中文字幕| 婷婷六月久久综合丁香| 淫妇啪啪啪对白视频| 国产一区二区在线av高清观看| 免费一级毛片在线播放高清视频| 亚洲av电影在线进入| 国产亚洲精品av在线| 天美传媒精品一区二区| 变态另类成人亚洲欧美熟女| 青草久久国产| 97碰自拍视频| 亚洲精品456在线播放app | 搡老妇女老女人老熟妇| 人妻制服诱惑在线中文字幕| 一个人观看的视频www高清免费观看| 欧美日韩综合久久久久久 | 热99re8久久精品国产| 日韩免费av在线播放| 日韩 亚洲 欧美在线| 人人妻人人澡欧美一区二区| 麻豆av噜噜一区二区三区| 丁香欧美五月| 久久久久精品国产欧美久久久| 男女床上黄色一级片免费看| 国产精品野战在线观看| 男女之事视频高清在线观看| 免费av不卡在线播放| 欧美成狂野欧美在线观看| 成人毛片a级毛片在线播放| 日本撒尿小便嘘嘘汇集6| 亚洲精品456在线播放app | 国产精品久久久久久人妻精品电影| 亚洲专区国产一区二区| 禁无遮挡网站| 亚洲av成人不卡在线观看播放网| 老司机午夜福利在线观看视频| 欧美bdsm另类| 国产色爽女视频免费观看| 国产一区二区亚洲精品在线观看| 中文资源天堂在线| 国产精品久久久久久人妻精品电影| 一级黄片播放器| 亚洲成人久久性| 97超级碰碰碰精品色视频在线观看| 久久国产乱子伦精品免费另类| 亚洲美女搞黄在线观看 | 国产午夜精品论理片| 免费av毛片视频| 中文字幕熟女人妻在线| 亚洲天堂国产精品一区在线| 国产不卡一卡二| 国产亚洲av嫩草精品影院| 在线十欧美十亚洲十日本专区| 亚洲av美国av| 亚洲av不卡在线观看| 日韩欧美国产一区二区入口| 可以在线观看的亚洲视频| 在线观看66精品国产| 国产精品自产拍在线观看55亚洲| 国产精品一区二区免费欧美| 亚洲激情在线av| 亚洲成a人片在线一区二区| 不卡一级毛片| 久久午夜福利片| 搞女人的毛片| АⅤ资源中文在线天堂| 99热这里只有是精品50| 免费av毛片视频| 日本在线视频免费播放| 美女免费视频网站| 夜夜夜夜夜久久久久| 亚洲自拍偷在线| xxxwww97欧美| 神马国产精品三级电影在线观看| 99精品久久久久人妻精品| 日本黄大片高清| 成人无遮挡网站| 免费一级毛片在线播放高清视频| 亚洲18禁久久av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲av.av天堂| 乱码一卡2卡4卡精品| 国产高清激情床上av| 精品午夜福利视频在线观看一区| 一级a爱片免费观看的视频| 欧美丝袜亚洲另类 | 午夜a级毛片| www.熟女人妻精品国产| 男女做爰动态图高潮gif福利片| 内地一区二区视频在线| 亚洲av免费在线观看| 午夜免费男女啪啪视频观看 | 舔av片在线| 成年版毛片免费区| 毛片女人毛片| 国产三级在线视频| 亚洲天堂国产精品一区在线| 午夜福利免费观看在线| 人妻制服诱惑在线中文字幕| 午夜精品在线福利| 亚洲国产精品999在线| 一本精品99久久精品77| 欧美性猛交╳xxx乱大交人| 麻豆成人av在线观看| 亚洲欧美日韩高清专用| 亚洲精品成人久久久久久| 国产又黄又爽又无遮挡在线| 国产精品爽爽va在线观看网站| 一个人看的www免费观看视频| 国产精华一区二区三区| 色在线成人网| 日本三级黄在线观看| 成年女人永久免费观看视频| 国产乱人伦免费视频| 国产淫片久久久久久久久 | 亚洲成a人片在线一区二区| 一a级毛片在线观看| 日日夜夜操网爽| 日本a在线网址| 三级国产精品欧美在线观看| 精品日产1卡2卡| 亚洲国产欧洲综合997久久,| 在线十欧美十亚洲十日本专区| 国产精华一区二区三区| 亚洲欧美日韩无卡精品| 亚洲天堂国产精品一区在线| 深爱激情五月婷婷| 女同久久另类99精品国产91| 免费在线观看成人毛片| av欧美777| 国产探花在线观看一区二区| 丝袜美腿在线中文| 欧美一区二区精品小视频在线| 99热这里只有精品一区| 国产精品三级大全| 亚洲精品一卡2卡三卡4卡5卡| 在线免费观看不下载黄p国产 | 97热精品久久久久久| 国产三级中文精品| 亚洲男人的天堂狠狠| 国产精品三级大全| .国产精品久久| 51国产日韩欧美| 国产主播在线观看一区二区| 麻豆一二三区av精品| 午夜福利在线观看吧| 日韩欧美国产一区二区入口| 乱码一卡2卡4卡精品| 极品教师在线免费播放| 国产极品精品免费视频能看的| 国产精品久久久久久久久免 | 级片在线观看| av在线天堂中文字幕| 亚洲精品在线美女| 老司机午夜福利在线观看视频| 亚洲熟妇熟女久久| 激情在线观看视频在线高清| 非洲黑人性xxxx精品又粗又长| 久久精品国产亚洲av香蕉五月| 99久久99久久久精品蜜桃| 亚洲精品久久国产高清桃花| 亚洲自拍偷在线| 啦啦啦观看免费观看视频高清| 一级毛片久久久久久久久女| 人妻夜夜爽99麻豆av| 成熟少妇高潮喷水视频| 精品久久国产蜜桃| 亚洲经典国产精华液单 | 他把我摸到了高潮在线观看| 韩国av一区二区三区四区| 亚洲五月婷婷丁香| 中文字幕人成人乱码亚洲影| 国产av麻豆久久久久久久| 亚洲欧美日韩卡通动漫| 国产三级中文精品| 午夜久久久久精精品| 一个人免费在线观看电影| 一进一出抽搐gif免费好疼| 蜜桃久久精品国产亚洲av| 舔av片在线| 精品国产亚洲在线| 中文字幕免费在线视频6| 久久久久久久午夜电影| 三级男女做爰猛烈吃奶摸视频| 日本黄色片子视频| 欧美精品国产亚洲| .国产精品久久| 夜夜夜夜夜久久久久| 国产av不卡久久| av天堂中文字幕网| 最新在线观看一区二区三区| 国产高清视频在线观看网站| 成人国产综合亚洲| 欧美色欧美亚洲另类二区| 欧美精品国产亚洲| 伦理电影大哥的女人| 一级作爱视频免费观看| 99久久精品热视频| 久久精品国产亚洲av涩爱 | 搡老岳熟女国产| 丁香六月欧美| 亚洲精品久久国产高清桃花| 搡老岳熟女国产| 精品久久久久久久人妻蜜臀av| 日本一二三区视频观看| 久久精品夜夜夜夜夜久久蜜豆|