• <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| 欧美精品一区二区大全| 黄色配什么色好看| 波多野结衣巨乳人妻| 特大巨黑吊av在线直播| 亚洲欧美中文字幕日韩二区| 精品久久久精品久久久| 色哟哟·www| 夜夜爽夜夜爽视频| 蜜桃亚洲精品一区二区三区| 国产亚洲精品久久久com| 国产av国产精品国产| 中国国产av一级| 亚洲欧洲日产国产| 国产一区二区在线观看日韩| 男的添女的下面高潮视频| 久久久久久久久久久免费av| 亚洲,欧美,日韩| 国产又色又爽无遮挡免| 亚洲精品一区蜜桃| 又爽又黄无遮挡网站| 精品一区二区三区视频在线| 午夜日本视频在线| 亚洲人成网站在线播| 日韩在线高清观看一区二区三区| 91午夜精品亚洲一区二区三区| videossex国产| 中文精品一卡2卡3卡4更新| 国内揄拍国产精品人妻在线| 国产黄频视频在线观看| 午夜日本视频在线| 男女边摸边吃奶| 久久久久久九九精品二区国产| 国内精品美女久久久久久| 日韩 亚洲 欧美在线| 美女内射精品一级片tv| 国产高清国产精品国产三级 | 亚洲,欧美,日韩| 日韩电影二区| 久久精品国产亚洲av天美| 久久午夜福利片| 色视频在线一区二区三区| 亚洲av男天堂| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 夫妻午夜视频| 久久影院123| 国产亚洲最大av| 久久久精品欧美日韩精品| a级毛片免费高清观看在线播放| 国产成人精品婷婷| 国产成年人精品一区二区| 亚洲一级一片aⅴ在线观看| 亚洲一级一片aⅴ在线观看| 搡女人真爽免费视频火全软件| 男女下面进入的视频免费午夜| 99热国产这里只有精品6| 精品国产乱码久久久久久小说| 亚洲av成人精品一二三区| 久久精品久久久久久噜噜老黄| 97人妻精品一区二区三区麻豆| 久久国产乱子免费精品| 91aial.com中文字幕在线观看| 日产精品乱码卡一卡2卡三| 中文字幕av成人在线电影| 亚洲欧美日韩无卡精品| 91午夜精品亚洲一区二区三区| 亚洲国产高清在线一区二区三| 伊人久久精品亚洲午夜| 九色成人免费人妻av| 国产精品成人在线| 亚洲欧美日韩东京热| 老师上课跳d突然被开到最大视频| 男女边吃奶边做爰视频| 免费高清在线观看视频在线观看| 三级经典国产精品| 亚洲精品一二三| 亚洲人成网站在线播| 色婷婷久久久亚洲欧美| av一本久久久久| 91aial.com中文字幕在线观看| 久久鲁丝午夜福利片| 又爽又黄a免费视频| 亚洲欧美中文字幕日韩二区| 一级毛片久久久久久久久女| 日韩中字成人| 国内精品美女久久久久久| 中文字幕人妻熟人妻熟丝袜美| 大片免费播放器 马上看| 欧美国产精品一级二级三级 | 国产永久视频网站| 日韩人妻高清精品专区| 国产伦在线观看视频一区| 99热这里只有是精品50| 久久99热6这里只有精品| 国产精品精品国产色婷婷| 美女主播在线视频| 黄片无遮挡物在线观看| 日韩,欧美,国产一区二区三区| 99久久人妻综合| 日韩 亚洲 欧美在线| 日本爱情动作片www.在线观看| 日韩av在线免费看完整版不卡| 在线播放无遮挡| 欧美丝袜亚洲另类| 简卡轻食公司| av在线播放精品| 久久久久久国产a免费观看| 爱豆传媒免费全集在线观看| 亚洲一级一片aⅴ在线观看| 亚洲伊人久久精品综合| 国产有黄有色有爽视频| 亚洲性久久影院| 日韩电影二区| 亚洲图色成人| 人妻少妇偷人精品九色| 久久6这里有精品| 欧美日韩精品成人综合77777| 熟女电影av网| 黄色一级大片看看| 91精品一卡2卡3卡4卡| 极品少妇高潮喷水抽搐| 午夜爱爱视频在线播放| 亚洲精品影视一区二区三区av| 97在线视频观看| 日韩精品有码人妻一区| 亚洲精品日韩在线中文字幕| 国产精品.久久久| 精品久久久久久久人妻蜜臀av| 久久久a久久爽久久v久久| 又爽又黄a免费视频| 麻豆精品久久久久久蜜桃| 插阴视频在线观看视频| 91精品伊人久久大香线蕉| 精品少妇久久久久久888优播| 又粗又硬又长又爽又黄的视频| 高清日韩中文字幕在线| 久久久成人免费电影| 日产精品乱码卡一卡2卡三| 国产精品一区二区性色av| 人体艺术视频欧美日本| 亚洲无线观看免费| 午夜视频国产福利| 亚洲不卡免费看| 特大巨黑吊av在线直播| 免费黄色在线免费观看| 欧美成人a在线观看| 欧美激情在线99| 欧美少妇被猛烈插入视频| 亚洲av日韩在线播放| 国产精品国产三级专区第一集| eeuss影院久久| 黄片无遮挡物在线观看| 一区二区三区精品91| 日日啪夜夜撸| 午夜福利高清视频| 久久久久久久久久人人人人人人| 成年版毛片免费区| 99re6热这里在线精品视频| 真实男女啪啪啪动态图| 国产一级毛片在线| 欧美xxxx性猛交bbbb| 女人十人毛片免费观看3o分钟| 久久久精品欧美日韩精品| 国产精品爽爽va在线观看网站| 欧美另类一区| 亚洲国产欧美人成| 中文精品一卡2卡3卡4更新| 日本-黄色视频高清免费观看| 在线a可以看的网站| 欧美日本视频| 亚洲国产精品国产精品| 女的被弄到高潮叫床怎么办| 国产视频内射| 一级黄片播放器| 国产 精品1| 久久午夜福利片| 丝袜美腿在线中文| 国产高清国产精品国产三级 | 欧美激情久久久久久爽电影| 性色av一级| 麻豆乱淫一区二区| 免费看av在线观看网站| 18禁动态无遮挡网站| 日本三级黄在线观看| 可以在线观看毛片的网站| 亚洲国产精品999| 一个人看的www免费观看视频| 国产熟女欧美一区二区| 神马国产精品三级电影在线观看| 日本色播在线视频| 99久久中文字幕三级久久日本| 日本爱情动作片www.在线观看| 一区二区三区乱码不卡18| 99热6这里只有精品| 亚洲精品一二三| 91午夜精品亚洲一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 91精品国产九色| 亚洲最大成人中文| 一级片'在线观看视频| 日韩欧美精品免费久久| 国产一区有黄有色的免费视频| 伦理电影大哥的女人| 高清日韩中文字幕在线| 在线观看免费高清a一片| 建设人人有责人人尽责人人享有的 | 国产老妇伦熟女老妇高清| 又爽又黄a免费视频| 老司机影院毛片| 精品国产三级普通话版| 国产免费视频播放在线视频| 精品亚洲乱码少妇综合久久| 波野结衣二区三区在线| 91aial.com中文字幕在线观看| 最近手机中文字幕大全| 国产久久久一区二区三区| 看免费成人av毛片| 亚洲天堂av无毛| 99九九线精品视频在线观看视频| 免费看日本二区| 国产成人午夜福利电影在线观看| 日产精品乱码卡一卡2卡三| av一本久久久久| 国产成人aa在线观看| 成人午夜精彩视频在线观看| 男女边摸边吃奶| 亚洲精品国产成人久久av| 亚洲婷婷狠狠爱综合网| 春色校园在线视频观看| 国产久久久一区二区三区| 欧美潮喷喷水| 久久热精品热| 男女国产视频网站| 色视频www国产| 特级一级黄色大片| 哪个播放器可以免费观看大片| 久久6这里有精品| av在线app专区| 精品国产一区二区三区久久久樱花 | 狂野欧美白嫩少妇大欣赏| 可以在线观看毛片的网站| 搡女人真爽免费视频火全软件| 亚洲国产精品成人综合色| 国产在线一区二区三区精| 免费观看在线日韩| 自拍欧美九色日韩亚洲蝌蚪91 | av又黄又爽大尺度在线免费看| 国产精品熟女久久久久浪| 免费少妇av软件| 91久久精品电影网| 久久久久久久午夜电影| 少妇熟女欧美另类| 高清视频免费观看一区二区| 91aial.com中文字幕在线观看| 久久人人爽av亚洲精品天堂 | 免费黄色在线免费观看| 亚洲无线观看免费| 制服丝袜香蕉在线| 欧美97在线视频| 久久精品国产亚洲av涩爱| 亚洲欧美日韩无卡精品| 久久97久久精品| 国产欧美日韩一区二区三区在线 | 熟女电影av网| 各种免费的搞黄视频| 国内精品宾馆在线| av专区在线播放| 两个人的视频大全免费| 国产黄a三级三级三级人| 黄色视频在线播放观看不卡| 秋霞伦理黄片| 久久热精品热| 美女xxoo啪啪120秒动态图| 人妻少妇偷人精品九色| 97超碰精品成人国产| 婷婷色综合大香蕉| 超碰av人人做人人爽久久| 干丝袜人妻中文字幕| 成人免费观看视频高清| 亚洲怡红院男人天堂| 最后的刺客免费高清国语| 亚洲精品一区蜜桃| 国产永久视频网站| 亚洲,一卡二卡三卡| 麻豆精品久久久久久蜜桃| 亚洲欧美日韩东京热| 日韩伦理黄色片| 色哟哟·www| 在现免费观看毛片| 久久久久久久亚洲中文字幕| 久久精品人妻少妇| 啦啦啦在线观看免费高清www| 久久久久国产精品人妻一区二区| 久久久久久久久久成人| 99久久精品热视频| 日本av手机在线免费观看| 在线免费观看不下载黄p国产| 另类亚洲欧美激情| 国产国拍精品亚洲av在线观看| 熟妇人妻不卡中文字幕| 国产午夜福利久久久久久| 亚洲精品国产成人久久av| 欧美日韩亚洲高清精品| 亚洲欧美精品自产自拍| 亚洲欧美日韩卡通动漫| 日韩视频在线欧美| 成人黄色视频免费在线看| 美女内射精品一级片tv| 亚洲最大成人中文| 国产亚洲91精品色在线| 国产 精品1| 99精国产麻豆久久婷婷| 欧美少妇被猛烈插入视频| 日韩一区二区三区影片| 日韩国内少妇激情av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品一区蜜桃| 欧美xxⅹ黑人| 狂野欧美激情性xxxx在线观看| 在线观看一区二区三区激情| 黄片无遮挡物在线观看| 人妻系列 视频| 22中文网久久字幕| 久久久久精品性色| 1000部很黄的大片| 亚洲婷婷狠狠爱综合网| 国产黄片美女视频| 在线观看免费高清a一片| 中文资源天堂在线| 亚洲一级一片aⅴ在线观看| 久久精品国产亚洲网站| 色综合色国产| 国产一区亚洲一区在线观看| 99热6这里只有精品| 99久久九九国产精品国产免费| 亚洲自偷自拍三级| 91精品国产九色| 精品国产一区二区三区久久久樱花 | 久久精品夜色国产| 天天一区二区日本电影三级| 久久久久久伊人网av| 国产淫片久久久久久久久| 国产成人午夜福利电影在线观看| 日韩中字成人| 国产精品熟女久久久久浪| 中文资源天堂在线| 丰满人妻一区二区三区视频av| 久久精品国产a三级三级三级| av网站免费在线观看视频| 色视频在线一区二区三区| 九九爱精品视频在线观看| 边亲边吃奶的免费视频| 免费大片18禁| 国产日韩欧美在线精品| 亚洲精品色激情综合| 国产精品蜜桃在线观看| 永久免费av网站大全| 日韩强制内射视频| 日日啪夜夜爽| 男人舔奶头视频| av国产久精品久网站免费入址| 熟女av电影| 色哟哟·www| 别揉我奶头 嗯啊视频| 亚洲av国产av综合av卡| 精品一区二区三区视频在线| 亚洲色图av天堂| 91久久精品国产一区二区三区| av女优亚洲男人天堂| 麻豆成人av视频| av又黄又爽大尺度在线免费看| 成人鲁丝片一二三区免费| 亚洲最大成人中文| 国产一区亚洲一区在线观看| 美女高潮的动态| 免费观看av网站的网址| 久久99热这里只频精品6学生| 97在线人人人人妻| 国产亚洲5aaaaa淫片| 插阴视频在线观看视频| 看免费成人av毛片| 少妇的逼好多水| 亚洲在线观看片| 婷婷色综合大香蕉| 亚洲欧美中文字幕日韩二区| 人妻制服诱惑在线中文字幕| 中文精品一卡2卡3卡4更新| 女人被狂操c到高潮| 国产有黄有色有爽视频| 搡老乐熟女国产| 丝瓜视频免费看黄片| 亚洲欧美精品自产自拍| 免费观看无遮挡的男女| 日韩一区二区视频免费看| 国产成人精品婷婷| 日本三级黄在线观看| 欧美精品人与动牲交sv欧美| 欧美xxxx黑人xx丫x性爽| a级一级毛片免费在线观看| 久久精品久久久久久噜噜老黄| 一级黄片播放器| 欧美xxⅹ黑人| 有码 亚洲区| 国产精品国产av在线观看| 久久久久久久久大av| 成年免费大片在线观看| 亚洲欧美精品自产自拍| videossex国产| 黄色配什么色好看| 美女xxoo啪啪120秒动态图| 国产一区有黄有色的免费视频| 亚洲国产最新在线播放| 久热这里只有精品99| 欧美另类一区| 在线精品无人区一区二区三 | 亚洲最大成人中文| 激情五月婷婷亚洲| 亚洲av成人精品一区久久| 蜜臀久久99精品久久宅男| 高清欧美精品videossex| 午夜福利网站1000一区二区三区| 国产免费一区二区三区四区乱码| 久久精品国产亚洲av涩爱| 亚洲精品色激情综合| 久久久久久久大尺度免费视频| 亚洲av欧美aⅴ国产| 日本爱情动作片www.在线观看| 成人欧美大片| 国产一区亚洲一区在线观看| 性色av一级| 欧美精品一区二区大全| 亚洲国产精品成人综合色| 国产精品人妻久久久久久| 自拍欧美九色日韩亚洲蝌蚪91 | 97人妻精品一区二区三区麻豆| 好男人在线观看高清免费视频| 国产黄片美女视频| 日韩成人av中文字幕在线观看| 少妇的逼水好多| 国产亚洲91精品色在线| 久久精品国产a三级三级三级| 18禁裸乳无遮挡动漫免费视频 | 欧美日韩视频精品一区| av女优亚洲男人天堂| 成人亚洲精品一区在线观看 | 精品久久久久久久久亚洲| 少妇 在线观看| 亚洲欧洲日产国产| 国产成人午夜福利电影在线观看| 色吧在线观看| 久久精品夜色国产| 国产成人freesex在线| www.av在线官网国产| 九草在线视频观看| 精品熟女少妇av免费看| av免费观看日本| 免费看a级黄色片| 欧美激情在线99| 久久久a久久爽久久v久久| av又黄又爽大尺度在线免费看| 国产免费视频播放在线视频| 简卡轻食公司| 成人一区二区视频在线观看| 国产精品99久久久久久久久| 内射极品少妇av片p| 视频区图区小说| 毛片女人毛片| 国产淫语在线视频| 少妇人妻一区二区三区视频| 91久久精品国产一区二区三区| 亚洲精品第二区| 大香蕉久久网| 777米奇影视久久| 秋霞在线观看毛片| 波野结衣二区三区在线| 18禁裸乳无遮挡动漫免费视频 | 国产乱来视频区| 国产成人精品福利久久| 亚洲天堂av无毛| 97精品久久久久久久久久精品| av在线观看视频网站免费| 丰满少妇做爰视频| 成人综合一区亚洲| 国产精品久久久久久精品古装| 色综合色国产| 免费不卡的大黄色大毛片视频在线观看| 少妇人妻精品综合一区二区| 岛国毛片在线播放| 亚洲av免费在线观看| 久久久久国产精品人妻一区二区| 国产黄片美女视频| 亚洲国产欧美在线一区| 午夜精品国产一区二区电影 | 亚洲三级黄色毛片| 国内少妇人妻偷人精品xxx网站| 综合色丁香网| 老师上课跳d突然被开到最大视频| 两个人的视频大全免费| 91久久精品国产一区二区三区| 中文字幕免费在线视频6| 国产爱豆传媒在线观看| 日日摸夜夜添夜夜爱| 国产成人a区在线观看| 深夜a级毛片| 男人爽女人下面视频在线观看| 岛国毛片在线播放| 国产亚洲午夜精品一区二区久久 | 色综合色国产| 大片电影免费在线观看免费| 国产一级毛片在线| 亚洲精品aⅴ在线观看| 国产免费视频播放在线视频| 联通29元200g的流量卡| 久久99蜜桃精品久久| 性色avwww在线观看| 国产成人a区在线观看| 亚洲一区二区三区欧美精品 | 国产精品久久久久久av不卡| 成人美女网站在线观看视频| 午夜福利在线观看免费完整高清在| 欧美日韩视频高清一区二区三区二| 亚洲av在线观看美女高潮| 1000部很黄的大片| 亚洲av一区综合| 免费看日本二区| 亚洲精品久久久久久婷婷小说| 亚洲天堂国产精品一区在线| 免费不卡的大黄色大毛片视频在线观看| 观看免费一级毛片| 日本wwww免费看| 亚洲av在线观看美女高潮| 亚洲伊人久久精品综合| 又黄又爽又刺激的免费视频.| 我的女老师完整版在线观看| 成人美女网站在线观看视频| 一级片'在线观看视频| 亚洲av成人精品一区久久| 欧美bdsm另类| 狂野欧美白嫩少妇大欣赏| 成年版毛片免费区| 有码 亚洲区| 真实男女啪啪啪动态图| 国产精品成人在线| 国产片特级美女逼逼视频| 欧美人与善性xxx| 综合色av麻豆| 国产乱来视频区| 久久久成人免费电影| 久久久久精品性色| 国产精品伦人一区二区| 熟女人妻精品中文字幕| 成人美女网站在线观看视频| 成年女人看的毛片在线观看| 97在线视频观看| 韩国高清视频一区二区三区| 亚洲av二区三区四区| 亚洲自偷自拍三级| 黄色视频在线播放观看不卡| 99热这里只有是精品50| 亚洲欧美日韩东京热| 亚洲内射少妇av| 亚洲国产成人一精品久久久| 少妇高潮的动态图| 日韩av免费高清视频| 在线免费观看不下载黄p国产| 国产在线一区二区三区精| 国产精品久久久久久久电影| 在线观看国产h片| 亚洲av不卡在线观看| 成人毛片60女人毛片免费| 免费观看性生交大片5| 免费黄频网站在线观看国产| 亚洲精品亚洲一区二区| 91在线精品国自产拍蜜月| 免费看av在线观看网站| 在线a可以看的网站| 亚洲av电影在线观看一区二区三区 | 久久久a久久爽久久v久久| 亚洲精品,欧美精品| 国内精品宾馆在线| 黄色配什么色好看| 啦啦啦在线观看免费高清www| 色视频www国产| 成人毛片60女人毛片免费| 久久午夜福利片| 亚洲国产精品999| xxx大片免费视频| 啦啦啦啦在线视频资源| 天堂中文最新版在线下载 | av在线天堂中文字幕| 精品亚洲乱码少妇综合久久| 免费高清在线观看视频在线观看| 久久久久九九精品影院| 国产成人午夜福利电影在线观看| 视频区图区小说| 水蜜桃什么品种好| 女人久久www免费人成看片| 国产亚洲午夜精品一区二区久久 | 国产白丝娇喘喷水9色精品| 国产亚洲午夜精品一区二区久久 | av天堂中文字幕网| 亚洲,一卡二卡三卡| 久久久a久久爽久久v久久| 国产亚洲5aaaaa淫片| 99久国产av精品国产电影| 麻豆精品久久久久久蜜桃| 男女啪啪激烈高潮av片| videossex国产| av在线亚洲专区| 国产一区二区亚洲精品在线观看| 久久精品国产亚洲av天美| 国产免费福利视频在线观看| 国精品久久久久久国模美| 熟女av电影|