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

    北半球大氣環(huán)流能量循環(huán)的氣候特征

    2015-02-25 12:52:08李新新管兆勇李明剛
    大氣科學(xué)學(xué)報(bào) 2015年4期
    關(guān)鍵詞:動(dòng)能

    李新新,管兆勇,李明剛

    (1.氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室(南京信息工程大學(xué)),江蘇 南京 210044;2.鹽城市氣象局,江蘇 鹽城 224005)

    李新新,管兆勇,李明剛.2015.北半球大氣環(huán)流能量循環(huán)的氣候特征[J].大氣科學(xué)學(xué)報(bào),38(4):531-539.doi:10.13878/j.cnki.dqkxxb.20130306002.

    Li Xin-xin,Guan Zhao-yong,Li Ming-gang.2015.Climatological features of atmospheric energy cycle in Northern Hemisphere[J].Trans Atmos Sci,38(4):531-539.(in Chinese).

    北半球大氣環(huán)流能量循環(huán)的氣候特征

    李新新1,2,管兆勇1,李明剛1

    (1.氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室(南京信息工程大學(xué)),江蘇 南京 210044;2.鹽城市氣象局,江蘇 鹽城 224005)

    摘要:利用1958—2011年NCEP/NCAR逐日再分析資料,根據(jù)Lorenz能量循環(huán)理論框架,分析了北半球大氣能量循環(huán)的年變化特征,在此基礎(chǔ)上給出了更具普適性的多年平均的大氣能量循環(huán)框圖。結(jié)果表明:北半球大氣能量循環(huán)的年變化特征十分明顯。大氣能量及能量轉(zhuǎn)換率均表現(xiàn)為冬季高、夏季低、春秋季過(guò)渡的演變特征;緯向平均有效位能、緯向平均動(dòng)能和渦動(dòng)動(dòng)能中有少許能量在冬季時(shí)由南半球向北半球進(jìn)行越赤道輸送,夏季時(shí)則由北半球向南半球輸送,而渦動(dòng)有效位能的輸送方向則與此相反;緯向平均有效位能的制造在秋季最大,渦動(dòng)有效位能的制造在夏季最大;動(dòng)能的耗散冬季最強(qiáng),夏季最弱。就年平均而言,相較于能量轉(zhuǎn)化過(guò)程,能量越赤道交換過(guò)程非常微弱。在經(jīng)向上,緯向平均有效位能主要分布于高緯地區(qū),緯向平均動(dòng)能主要分布于中低緯地區(qū),而渦動(dòng)能量主要貯存在中緯和高緯地區(qū);此外,能量轉(zhuǎn)化過(guò)程一般在中緯度地區(qū)較活躍。

    關(guān)鍵詞:能量循環(huán);有效位能;動(dòng)能;年變化

    0引言

    通常,人們側(cè)重于研究半球或全球范圍內(nèi)大氣運(yùn)動(dòng)的平均特征,即側(cè)重于研究環(huán)流的平衡狀態(tài),且用大氣的位勢(shì)高度、溫度以及風(fēng)速等基本要素場(chǎng)來(lái)表述(施能和朱乾根,1996;楊哲等,2012)。而大氣運(yùn)動(dòng)需要外界不斷地供給能量,因此,從能量學(xué)的觀點(diǎn)來(lái)研究維持大氣環(huán)流平衡狀態(tài)的過(guò)程—大氣能量循環(huán)是非常重要的,它能正確地反映氣候系統(tǒng)中多圈層相互作用和轉(zhuǎn)化的真實(shí)物理機(jī)制(Lorenz,1967),是診斷大氣環(huán)流的一個(gè)非常有用的工具。

    大氣通過(guò)吸收太陽(yáng)輻射使得全位能增加,進(jìn)而轉(zhuǎn)化為維持風(fēng)場(chǎng)的動(dòng)能。而實(shí)際觀測(cè)結(jié)果表明,大氣的動(dòng)能還不到全位能的千分之一(Lorenz,1967),這說(shuō)明大氣中的全位能只有很小一部分能被釋放而轉(zhuǎn)化為動(dòng)能。由此Lorenz(1955)提出了有效位能的概念,將其定義為實(shí)際大氣的全位能與溫度場(chǎng)經(jīng)絕熱調(diào)整后具有水平正壓、穩(wěn)定層結(jié)的最小全位能之差,并被廣泛應(yīng)用于實(shí)際(魏民和仇永炎,1997;羅連升和楊修群,2003)。在有效位能概念的基礎(chǔ)上,Lorenz進(jìn)一步將有效位能和動(dòng)能分別分解為緯向平均和渦旋部分,從而發(fā)展了著名的四種能量模態(tài)收支框架,它包含了有效位能的產(chǎn)生、與動(dòng)能間的轉(zhuǎn)化以及動(dòng)能的耗散過(guò)程。在該框架中,有效位能的產(chǎn)生主要是由于不同的非絕熱加熱(輻射、感熱和潛熱)造成的,平均能量與擾動(dòng)能量之間存在相互轉(zhuǎn)化,而擾動(dòng)能量之間的轉(zhuǎn)化才是真正完成有效位能和動(dòng)能之間轉(zhuǎn)化的關(guān)鍵所在(高麗和李建平,2007),對(duì)大氣環(huán)流的維持和演變起著重要作用(吳洪寶,1982;李俊,1990;袁建強(qiáng),1991)。

    Lorenz能量框圖是描述大氣環(huán)流能量學(xué)的最佳框架。然而由于早期觀測(cè)技術(shù)和觀測(cè)手段的限制,氣象資料呈現(xiàn)非均勻性分布,尤其是在赤道及南半球地區(qū),資料嚴(yán)重缺乏,使得大氣能量循環(huán)的研究受到諸多限制。Krueger et al.(1965)、Wiin-Nielsen(1967)、Oort and Peixoto(1974)以及Peixoto and Oort(1974)曾根據(jù)能量循環(huán)框圖對(duì)方程組進(jìn)行計(jì)算,試圖分析北半球大氣的能量循環(huán)過(guò)程,但由于所用的資料均局限在20°N以北,因此計(jì)算結(jié)果并不可靠。唯一的一次估算全球大氣的能量循環(huán)(Oort,1983)所采用的是10 a(1963—1973年)探空數(shù)據(jù)集,但由于南半球地區(qū)測(cè)站稀疏和不良觀測(cè)條件等的約束,資料分布不均勻,且部分測(cè)站的數(shù)據(jù)差異較大,使得研究結(jié)果的準(zhǔn)確性較低。近年來(lái)隨著資料同化技術(shù)和衛(wèi)星等觀測(cè)手段的發(fā)展,氣象資料的精度得到了很大程度的提高,同時(shí)資料也具備了均勻的水平分布條件。各種再分析資料的出現(xiàn)為氣象研究工作提供了極大的方便。Li et al.(2006,2007)、Marques et al.(2009)以及Kim and Kim(2012)利用不同的再分析資料(NCEP-R2、ERA-40和JRA-25)診斷分析了全球大氣的能量循環(huán)過(guò)程,并對(duì)幾種再分析資料對(duì)大氣運(yùn)動(dòng)的描述能力進(jìn)行了對(duì)比。

    此外,許多學(xué)者還擴(kuò)展了原始的能量循環(huán)框架研究。Saltzman(1957,1970)在波數(shù)域上對(duì)大氣能量進(jìn)行分解,得到的波數(shù)域方程更能體現(xiàn)出分波能量的貢獻(xiàn)以及波緯相互作用。Wiin-Nielsen and Brown(1962)、Smagorinsky(1963)以及Wiin-Nielsen and Chen(1993)在垂直方向上對(duì)有效位能和動(dòng)能進(jìn)一步分解,從而得到能量的正壓和斜壓分量。除了對(duì)大氣能量循環(huán)參數(shù)在空間上的分解外,Sheng and Derome(1991)在頻率域上將平均運(yùn)動(dòng)分解成3個(gè)頻帶,即季節(jié)平均、低頻和高頻渦動(dòng),得到的頻率域方程可用來(lái)分析各個(gè)季節(jié)的能量循環(huán)狀況。這些工作分別從不同角度對(duì)大氣能量進(jìn)行了劃分和深入研究。

    然而,關(guān)于能量循環(huán)年變化特征的研究,目前還沒(méi)有詳細(xì)的結(jié)果。在再分析資料出現(xiàn)之前,Oort and Peixoto(1974)以及Peixoto and Oort(1974)曾利用測(cè)站資料對(duì)能量循環(huán)的年變化特征進(jìn)行過(guò)分析,但由于資料的時(shí)間長(zhǎng)度較短(1958—1963年),且垂直速度無(wú)法直接獲得,只能對(duì)循環(huán)框架中部分參數(shù)的年變化特征給予定性說(shuō)明,因而不具備普遍性。20世紀(jì)80年代之后,Ulbrich and Speth(1991)、張韜等(2006)分別利用ECMWF和NCEP/NCAR再分析資料對(duì)全球大氣1月和7月的能量循環(huán)過(guò)程進(jìn)行了分析,但并沒(méi)有涉及能量循環(huán)年變化特征的研究。Hu et al.(2004)對(duì)比了20世紀(jì)80年代前后北半球大氣在各個(gè)季節(jié)的能量循環(huán),但對(duì)其演變過(guò)程并未作過(guò)多說(shuō)明。

    基于此,本文將根據(jù)Peixoto and Oort(1974)定義的能量循環(huán)公式,從傳統(tǒng)的能量循環(huán)框架出發(fā),利用1958—2011年NCEP/NCAR再分析資料對(duì)北半球大氣的能量循環(huán)及其年變化特征進(jìn)行定性和定量描述,從能量學(xué)的觀點(diǎn)揭示大氣環(huán)流的運(yùn)動(dòng)規(guī)律。

    1資料及方法

    采用美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心與美國(guó)國(guó)家大氣研究中心(NCEP/NCAR)提供的1958—2011年緯向風(fēng)u、經(jīng)向風(fēng)v、p坐標(biāo)垂直速度ω、溫度T和位勢(shì)高度Z的逐日再分析資料(Kalnay et al.,1996),水平分辨率均為2.5°×2.5°經(jīng)緯度網(wǎng)格,垂直方向上1 000~10 hPa共17個(gè)等壓面層,其中垂直速度1 000~100 hPa共12層。

    根據(jù)Lorenz能量循環(huán)框圖(圖1),在混合時(shí)空域上將有效位能P和動(dòng)能K分別分解為緯向平均部分(PM和KM)和渦動(dòng)部分(PE和KE),而渦動(dòng)又包括了定常渦動(dòng)(時(shí)間平均,本文取為月平均,PSE和KSE)和瞬變渦動(dòng)(與時(shí)間平均的離差,PTE和KTE)兩部分。文中對(duì)不同能量的貯存和轉(zhuǎn)換項(xiàng)的計(jì)算是根據(jù)Peixoto and Oort(1974)的能量循環(huán)公式來(lái)實(shí)現(xiàn)的,具體的能量收支方程如公式(1),方程中各項(xiàng)具體的表達(dá)式可參閱文獻(xiàn)Oort and Peixoto(1974)。

    (1)

    圖1 1958—2011年北半球多年平均大氣的能量循環(huán)框圖(方框內(nèi)為貯蓄項(xiàng),單位:105 J·m-2;箭矢旁為轉(zhuǎn)換項(xiàng),單位:W·m-2;箭矢表示正值的轉(zhuǎn)換方向)Fig.1 Diagram of the Northern Hemispheric energy cycle for monthly mean circulation averaged over the period 1958—2011(Units of reservoirs(in boxes) are in 105 J·m-2,and conversions(near arrows) are in W·m-2.The arrows indicate the orientations of the conversion from one reservoir to another)

    本文在估算不同貯存(PM和KM)的制造項(xiàng)(G(PM)和G(PE))和耗散項(xiàng)(D(KM)和D(KE))的半球平均時(shí),是根據(jù)能量收支方程并忽略不同貯存項(xiàng)的時(shí)間導(dǎo)數(shù)(如?PM/?t,數(shù)值非常小)而作為剩余項(xiàng)來(lái)得到的。Ulbrich and Speth(1991)利用ECMWF再分析資料所得的結(jié)果表明,上述對(duì)不同貯存項(xiàng)的時(shí)間導(dǎo)數(shù)的忽略是完全合理的。

    這里對(duì)文中逐月估算的月平均大氣環(huán)流的能量循環(huán)的計(jì)算過(guò)程進(jìn)行具體說(shuō)明。在任意一年的某一月份,變量(即u、v、ω、T、Z)的時(shí)間平均是指由逐日資料計(jì)算得到的月平均,時(shí)間距平則是相應(yīng)的月距平,在計(jì)算涉及瞬變渦動(dòng)過(guò)程的物理量時(shí)(如〈T′2〉和〈v′T′〉),是對(duì)單變量距平值的平方或某兩個(gè)變量距平值的乘積作月平均(即為求變量方差或協(xié)方差的過(guò)程);然后對(duì)公式中(Oort and Peixoto,1974)的各項(xiàng)在緯向上進(jìn)行分解,最后經(jīng)過(guò)簡(jiǎn)單的數(shù)學(xué)運(yùn)算得到能量循環(huán)參數(shù)在該月的月平均值。對(duì)上述結(jié)果進(jìn)行多年平均處理即可得到參數(shù)在該月的氣候平均值,再對(duì)各月的氣候平均值進(jìn)行年平均處理即得到如圖1所示的能量循環(huán)框架。對(duì)于逐年估算的年平均大氣的能量循環(huán),只需計(jì)算時(shí)將公式中變量的時(shí)間平均取為年平均,時(shí)間距平即為相應(yīng)的年距平,其他步驟不變,從而得到參數(shù)在任意一年的年平均值。

    不同于全球大氣的能量循環(huán)研究,在半球大氣中,有必要考慮赤道邊界的能量流動(dòng)過(guò)程。Peixoto and Oort(1974)、Oort and Peixoto(1974)指出,對(duì)于大氣能量的年循環(huán)而言,在忽略水汽輸送的情況下,北半球可以看作是一個(gè)孤立的閉合系統(tǒng),即赤道邊界的能量流動(dòng)過(guò)程可以忽略;但對(duì)于季節(jié)以及逐月的能量循環(huán)過(guò)程,尤其是在冬夏月份,北半球是一個(gè)開(kāi)放的系統(tǒng),其邊界上緯向平均有效位能的交換相較于整個(gè)半球的能量收支非常重要?;谶@一思想,本文對(duì)有效位能和動(dòng)能在赤道邊界上的流通率進(jìn)行了估算,以得到更為準(zhǔn)確的北半球大氣能量循環(huán)框圖。

    天氣尺度的渦動(dòng)(空間尺度達(dá)1 000 km,生命史達(dá)幾天,如氣旋、反氣旋)在大氣能量循環(huán)過(guò)程中扮演著重要角色,而那些空間和時(shí)間尺度小于NCEP/NCAR逐日再分析資料分辨率(分別為250 km和1 d)的渦動(dòng),相比之下,則是次網(wǎng)格運(yùn)動(dòng)(Lorenz,1955),暫不作考慮。對(duì)大部分天氣尺度渦動(dòng)而言,其生命史都小于1 mon,因此本文所估算的月平均大氣的能量循環(huán)中,其瞬變渦動(dòng)包括了大部分的渦動(dòng)能量貢獻(xiàn),這也是目前研究大氣能量循環(huán)最常選取的時(shí)間尺度。

    2大氣能量循環(huán)的年變化特征

    2.1 大氣能量的年變化

    對(duì)北半球大氣逐月估算的能量循環(huán)過(guò)程進(jìn)行多年氣候平均,得到的各項(xiàng)能量參數(shù)如表1所示。從表中可以看出,大氣中各能量的貯存均是在冬季最高,夏季最低,這與Peixoto and Oort(1974)的結(jié)果一致。其中,除了緯向平均有效位能(PM)外,其他三種形式的能量(PE、KM、KE)具有大致相當(dāng)?shù)牧恐?且遠(yuǎn)小于PM的貯存。有效位能與動(dòng)能各自渦動(dòng)量與緯向平均量的比值(PE/PM和KE/KM)則是在夏季高,冬季低(且夏季的比值約為冬季情況下的2倍),說(shuō)明緯向平均量的年變化幅度比渦動(dòng)量的年變化幅度相對(duì)要大。年平均而言(表中Ave1所示),PE的貯存不到PM的1/5,而KE的貯存則比KM大,其比值為1.26??偟膩?lái)說(shuō),大氣中有效位能的貯存明顯大于動(dòng)能的貯存,K/P的比值也是在夏季高于冬季,平均為0.31。

    此外,表1最后一行(Ave2)還給出了逐年估算的各項(xiàng)能量,與上述逐月估算的年平均結(jié)果相比,大氣的緯向平均能量有所減少,而渦動(dòng)能量有所增加,從而導(dǎo)致渦動(dòng)量與緯向平均量的比值顯著增大。從理論上來(lái)講,年平均環(huán)流的能量循環(huán)中瞬變渦動(dòng)項(xiàng)包含了時(shí)間尺度在1 d~1 a的所有渦動(dòng)的能量貢獻(xiàn),與月平均的計(jì)算結(jié)果相比,其中增加了時(shí)間尺度大于30 d的低頻渦動(dòng)的貢獻(xiàn),因此,渦動(dòng)能量會(huì)有所增加;而依據(jù)公式(1),在某一時(shí)段內(nèi)大氣的總能量是基本不變的,渦動(dòng)量的增加必然會(huì)伴隨著緯向平均量的減少,從而說(shuō)明了以上觀測(cè)結(jié)果。

    為了更直觀地反映各能量隨季節(jié)的演變特征,將表1的結(jié)果繪制成圖。圖2給出了54 a平均的緯向平均、渦動(dòng)有效位能和動(dòng)能以及總能量的逐月變化曲線,可以看出,PM(圖2a)、PE、KM和KE(圖2b)具有明顯的年變化特征,且變化幾乎同位相,均是冬季高,夏季低,春季和秋季處于過(guò)渡階段,全年呈一次很明顯的振動(dòng),從高能量到低能量再回到高能量。其中,PM具有最大的年變化值,1月和7月的差值高達(dá)40×105J·m-2以上,而PE、KM和KE三者的振幅相當(dāng)。另外,緯向平均能量的貯存最小出現(xiàn)在夏季的7月,而渦動(dòng)能量的貯存最小則出現(xiàn)在8月。總體來(lái)講,北半球冬季時(shí)能量最高,夏季時(shí)能量最低(圖2b)。

    表11958—2011年北半球大氣各能量的逐月平均值和年平均值

    Table 1The monthly and yearly mean energy reservoirs in the Northern Hemisphere for the period of 1958—2011

    月份PM/(105J·m-2)PE/(105J·m-2)KM/(105J·m-2)KE/(105J·m-2)P/(105J·m-2)K/(105J·m-2)PE/PMKE/KMK/P1月55.209.0210.219.3064.2219.500.160.910.302月55.188.3610.278.8063.5419.060.150.860.303月51.106.958.797.8058.0516.590.140.890.294月41.735.536.336.5947.2712.920.131.040.275月29.754.913.975.7234.669.680.171.440.286月18.894.632.755.1423.527.890.251.870.347月13.194.112.664.4817.307.140.311.680.418月15.123.762.914.3318.887.240.251.490.389月22.654.123.345.0226.778.360.181.500.3110月33.975.204.496.1739.1710.660.151.370.2711月45.306.836.637.5052.1314.130.151.130.2712月52.098.368.828.7860.4517.600.160.990.29Ave136.185.985.936.6342.1612.570.181.260.31Ave231.8612.854.478.0544.7112.520.401.800.28

    注:Ave1表示逐月估算(即公式中氣象要素的時(shí)間平均取為月平均)的能量各項(xiàng)進(jìn)行多年平均后的年平均結(jié)果;Ave2表示逐年估算(要素的時(shí)間平均取為年平均)的能量各項(xiàng)進(jìn)行多年平均后的結(jié)果.

    圖2 1958—2011年北半球大氣各能量逐月平均值的變化(單位:105 J·m-2)Fig.2 The monthly mean energy reservoirs in the Northern Hemisphere for the period of 1958—2011(units:105 J·m-2)

    2.2 能量轉(zhuǎn)化過(guò)程的年變化

    能量轉(zhuǎn)化過(guò)程也具有顯著的年變化特征(圖3)。與能量貯存類似,能量轉(zhuǎn)化率一般也是冬季高、夏季低,春秋季居中。其中,緯向平均有效位能向渦動(dòng)有效位能的轉(zhuǎn)化C(PM,PE)具有最大的年變化值,1月和7月的差值達(dá)2 W·m-2以上;渦動(dòng)有效位能向渦動(dòng)動(dòng)能的轉(zhuǎn)化C(PE,KE)其次(圖3a)。除了緯向平均有效位能向緯向平均動(dòng)能的轉(zhuǎn)化C(PM,KM)在初夏時(shí)出現(xiàn)負(fù)值外(圖3b),其他轉(zhuǎn)化項(xiàng)全年都為正值。這些說(shuō)明了這樣一個(gè)事實(shí),即能量轉(zhuǎn)換過(guò)程一般是從緯向平均有效位能開(kāi)始,通過(guò)渦動(dòng),有效位能轉(zhuǎn)化為渦動(dòng)動(dòng)能,最后轉(zhuǎn)化為緯向平均動(dòng)能,而緯向平均有效位能與緯向平均動(dòng)能之間的轉(zhuǎn)化較少,且轉(zhuǎn)化方向隨季節(jié)而變,由此可見(jiàn)渦動(dòng)能量之間的轉(zhuǎn)化才是真正完成有效位能向動(dòng)能轉(zhuǎn)化的關(guān)鍵所在。

    圖3 1958—2011年北半球大氣各能量轉(zhuǎn)換項(xiàng)逐月平均值的變化(單位:W·m-2)Fig.3 The monthly mean energy conversion rates in the Northern Hemisphere for the period of 1958—2011(units:W·m-2)

    2.3 能量越赤道交換過(guò)程的年變化

    能量越赤道交換過(guò)程亦具有顯著的年變化(圖4)。緯向平均有效位能的越赤道交換B(PM)(圖4a)在半球大氣的能量循環(huán)中起著非常重要的作用,其振幅與C(PM,PE)相當(dāng),而其他三種形式能量的越赤道交換過(guò)程(圖4b)非常微弱,比B(PM)低一個(gè)量級(jí),說(shuō)明了能量的越赤道交換以PM為主。此外,單從各自的逐月變化曲線來(lái)看,PM、KM和KE越赤道交換過(guò)程的年變化具有大致相同的位相,能量一般都是在冬季時(shí)由南半球向北半球輸送(圖中表現(xiàn)為正值),夏季時(shí)能量則由北半球向南半球輸送(表現(xiàn)為負(fù)值);而PE越赤道交換過(guò)程的年變化情況與此相反,即冬季時(shí)能量由北半球向南半球輸送,夏季時(shí)能量由南半球向北半球輸送。

    2.4 能量制造和耗散過(guò)程的年變化

    作為能量平衡方程剩余項(xiàng)計(jì)算得到的有效位能制造率和動(dòng)能耗散率的季節(jié)變化如圖5所示,可以看出,緯向平均有效位能的制造G(PM)在秋季時(shí)達(dá)到最大(圖5a),這可能是由于在該季節(jié)降水與溫度之間的緯度相關(guān)性較高的緣故;渦動(dòng)有效位能的制造G(PE)則是在夏季最大(圖5a),說(shuō)明夏季的潛熱釋放過(guò)程最為強(qiáng)烈;動(dòng)能耗散過(guò)程D(KM)和D(KE)具有相同的年變化特征(圖5b),均是冬季高,夏季低,這可以作為大氣環(huán)流強(qiáng)弱程度的體現(xiàn),冬季時(shí)環(huán)流較強(qiáng),動(dòng)能耗散率大,夏季時(shí)環(huán)流較弱,動(dòng)能耗散率也小。

    圖5 1958—2011年北半球大氣各能量制造和耗散項(xiàng)逐月平均值的變化(單位:W·m-2)Fig.5 The monthly mean rates of generation and dissipation terms of each energy in the Northern Hemisphere for the period of 1958—2011(units:W·m-2)

    3北半球氣候平均的能量循環(huán)

    將上述逐月估算的能量循環(huán)過(guò)程再進(jìn)行年平均,得到如圖1所示的經(jīng)典能量循環(huán)框架??梢钥闯?盡管逐月的PM在南北半球的交換顯著,但年平均而言,B(PM)的值與其他三種形式的能量一致,與能量間的轉(zhuǎn)化率相比都非常小。其中,KM和KE均是由北半球向南半球輸送,交換率分別為0.02和0.01 W·m-2;PM、PE的交換率均為0.03 W·m-2,PE是由南半球向北半球輸送,但PM卻是由北半球向南半球輸送。

    圖1表明,首先,由于太陽(yáng)輻射的緯度差異(低緯輻射加熱、高緯輻射冷卻)和暖區(qū)的潛熱釋放作用使得PM和PE的制造率分別為1.77和0.29 W·m-2;其次,PM通過(guò)斜壓擾動(dòng)(對(duì)感熱的輸送)過(guò)程向PE轉(zhuǎn)化,轉(zhuǎn)化率為1.53 W·m-2;而PE再以1.85 W·m-2的功率向KE進(jìn)行斜壓轉(zhuǎn)化(通過(guò)中緯度擾動(dòng)所形成的暖空氣上升和冷空氣下沉運(yùn)動(dòng));其中大部分KE由于摩擦而耗散,耗散率為1.41 W·m-2,剩余的KE則通過(guò)正壓過(guò)程向KM轉(zhuǎn)化(中緯度Rossby波結(jié)構(gòu)對(duì)西風(fēng)角動(dòng)量的輸送),轉(zhuǎn)化率為0.43 W·m-2;PM向KM的直接轉(zhuǎn)化率較小,為0.21 W·m-2;最后直接和間接轉(zhuǎn)化而來(lái)的KM以0.62 W·m-2功率被耗散掉,從而完成大氣的整個(gè)能量循環(huán)。在該過(guò)程中,PM、PE、KM、KE的貯存分別為36.18×105J·m-2、5.98×105J·m-2、5.93×105J·m-2和6.63×105J·m-2。

    與Peixoto and Oort(1974)以及Oort and Peixoto(1974)的研究結(jié)果相比,大氣能量中除了緯向平均動(dòng)能的貯存偏大外,其他三種形式能量的貯存均偏小;轉(zhuǎn)換過(guò)程中C(PM,KM)的轉(zhuǎn)化方向發(fā)生了變化,即能量由緯向平均動(dòng)能向緯向平均有效位能轉(zhuǎn)化,而其他轉(zhuǎn)化項(xiàng)的方向未發(fā)生變化。

    4大氣能量循環(huán)的緯向平均分布

    4.1 大氣能量的緯向平均分布

    為了更清晰地了解能量循環(huán)過(guò)程的季節(jié)變化,以下從空間上對(duì)能量各項(xiàng)進(jìn)行分析。圖6給出了冬季(12、1、2月的平均)、夏季(6、7、8月的平均)和年平均情況下北半球大氣各能量的緯向平均分布。從圖6a可以看出,緯向平有效位能PM在35°N附近具有最小值,這是由于該緯度上的入射太陽(yáng)輻射和地球向外長(zhǎng)波輻射趨于平衡;從數(shù)學(xué)上來(lái)講,就是該緯帶上空氣的溫度與北半球的平均溫度近似相等,從而使得PM的貯存值為零。PM的最大值出現(xiàn)在高緯地區(qū),說(shuō)明該區(qū)域的溫度擾動(dòng)(與北半球平均的偏差)最為強(qiáng)烈。冬季和夏季均是如此,但冬季PM的貯存更多,年平均的情況介于兩者之間。

    圖6 北半球大氣各能量PM(a)、PE(b)、KM(c)、KE(d)的緯向平均分布(實(shí)線表示年平均情形,短虛線表示冬季情形,長(zhǎng)虛線表示夏季情形;單位:105 J·m-2)Fig.6 The meridional contributions for different energies in the Northern Hemisphere(Solid lines indicate the annual mean,short dashed lines indicate the winter mean,and long dashed lines indicate the summer mean;units:105 J·m-2)  a.zonal mean available potential energy PM;b.eddy available potential energy PE;c.zonal mean kinetic energy KM;d.eddy kinetic energy KE

    渦動(dòng)有效位能PE的貯存(圖6b)在冬夏季節(jié)同樣具有相似的分布型。在數(shù)值上,年平均的結(jié)果介于兩者之間??梢钥闯?PE在赤道地區(qū)數(shù)值最小,從10°N開(kāi)始首先隨著緯度的增加而單調(diào)遞增,到了65°N附近達(dá)到最大值,之后又隨緯度的增加而逐漸減小。

    與PM的分布相反,緯向平均動(dòng)能KM的貯存(圖6c)則是冬季時(shí)在30°N附近達(dá)到最大值,對(duì)應(yīng)于副熱帶急流的位置,而在低緯和高緯地區(qū),其貯存非常小。到了夏季,由于急流位置的變化,KM的分布也發(fā)生了相應(yīng)的變化,極大值的位置由30°N北跳到了45°N附近,且數(shù)值明顯比冬季的小。年平均情況下KM的貯存以及極大值分布的位置位于冬夏季節(jié)之間。

    與PE類似,渦動(dòng)動(dòng)能KE的分布(圖6d)在冬夏季節(jié)較為一致,且均是先隨緯度而有所增加,在中緯度地區(qū)達(dá)到最大,到高緯度之后又有所減小。年平均情況處于兩者之間。

    4.2 能量轉(zhuǎn)化過(guò)程的緯向平均分布

    冬季、夏季和年平均情況下能量轉(zhuǎn)化過(guò)程的緯向平均分布可見(jiàn)圖7。由圖7a可看出,冬季時(shí)緯向平均有效位能向渦動(dòng)有效位能的轉(zhuǎn)換C(PM,PE)最活躍的區(qū)域出現(xiàn)在中緯度地區(qū),在該區(qū)域轉(zhuǎn)化率達(dá)到最大。在低緯和高緯地區(qū),這一轉(zhuǎn)化過(guò)程非常微弱。夏季的情形與冬季相似,能量轉(zhuǎn)化過(guò)程主要也發(fā)生在中緯度地區(qū),但其強(qiáng)度明顯減弱,年平均的情形介于兩者之間。

    圖7 北半球大氣各能量轉(zhuǎn)化項(xiàng)C(PM,PE)(a)、C(PE,KE)(b)、C(KE,KM)(c)、C(PM,KM)(d)的緯向平均分布(實(shí)線表示年平均情形,短虛線表示冬季情形,長(zhǎng)虛線表示夏季情形,點(diǎn)線表示0值;單位:W·m-2)Fig.7 The meridional contributions for energy conversions in the Northern Hemisphere(Solid lines indicate the annual mean,short dashed lines indicate the winter mean,long dashed lines indicate the summer mean,and dotted lines mean the value of zero;units:W·m-2)  a.C(PM,PE);b.C(PE,KE);c.C(KE,KM);d.C(PM,KM)

    圖7b顯示,渦動(dòng)有效位能向渦動(dòng)動(dòng)能的轉(zhuǎn)換C(PE,KE)在整個(gè)北半球都比較活躍,尤以中緯度為最,冬夏季節(jié)均是如此,但冬季時(shí)這一過(guò)程的轉(zhuǎn)化率明顯高于夏季,年平均的結(jié)果居中。從數(shù)值上來(lái)看,在整個(gè)北半球,該轉(zhuǎn)換過(guò)程的轉(zhuǎn)化率始終為正值,可見(jiàn)不管是冬季還是夏季,能量總是由渦動(dòng)有效位能向渦動(dòng)動(dòng)能轉(zhuǎn)化的,這一過(guò)程是維持大氣環(huán)流運(yùn)動(dòng)的關(guān)鍵。

    冬季渦動(dòng)動(dòng)能向緯向平均動(dòng)能的轉(zhuǎn)換C(KE,KM)(圖7c)在低緯和中緯各有一大值區(qū),以30°N為界,在其南側(cè),大值區(qū)為正值,說(shuō)明能量是由渦動(dòng)動(dòng)能向緯向平均動(dòng)能轉(zhuǎn)化的;而在其北側(cè)的大值區(qū)為負(fù)值,表示能量的轉(zhuǎn)化方向發(fā)生了反轉(zhuǎn)。夏季的情況與冬季類似,但正負(fù)極值區(qū)臨界緯度的位置發(fā)生了變化,北跳到了45°N附近,對(duì)應(yīng)于冬夏季節(jié)急流位置的變化。在強(qiáng)度上,夏季的能量轉(zhuǎn)換過(guò)程小于冬季。年平均情況下,中緯地區(qū)的負(fù)值區(qū)消失,表明能量轉(zhuǎn)化方向隨緯度不發(fā)生變化。

    緯向平均有效位能向緯向平均動(dòng)能的轉(zhuǎn)換C(PM,KM)(圖7d)反映的是三圈環(huán)流的作用??梢钥闯?不管是冬季還是夏季,在30°N以南,這一過(guò)程的轉(zhuǎn)化率為正值,低緯環(huán)流的作用使得能量由緯向平均有效位能向緯向平均動(dòng)能轉(zhuǎn)化;在30~65°N區(qū)域,其轉(zhuǎn)化率為負(fù)值,中緯環(huán)流的作用使得能量由緯向平均動(dòng)能向緯向平均有效位能轉(zhuǎn)化;而在65°N以北,高緯環(huán)流的作用與低緯環(huán)流一致,轉(zhuǎn)換率為正值,但其數(shù)值卻明顯小于低緯和中緯地區(qū)的轉(zhuǎn)化率,說(shuō)明高緯環(huán)流的作用非常弱,因此這一轉(zhuǎn)化過(guò)程的方向主要取決于低緯環(huán)流和中緯環(huán)流的凈作用。

    5結(jié)論

    1)大氣能量循環(huán)具有明顯的年變化特征。大氣能量及其轉(zhuǎn)換過(guò)程的年變化幾乎同位相,均是冬季高、夏季低,全年呈一次很明顯的波動(dòng)。緯向平均有效位能的振幅明顯大于其他三種形式能量的振幅。能量轉(zhuǎn)化過(guò)程中C(PM,PE)具有最大的年變化值。能量的越赤道交換過(guò)程主要表現(xiàn)為緯向平均有效位能在南北半球間的交換,其振幅與C(PM,PE)相當(dāng)。PM、KM和KE越赤道交換過(guò)程的年變化大致相同,冬季時(shí)能量由南半球向北半球輸送,夏季則反之;而PE越赤道交換過(guò)程的年變化特征與此相反。緯向平均有效位能的制造率在秋季時(shí)最大,這可能是由于在該季節(jié)降水與溫度之間的緯度相關(guān)性較高;渦動(dòng)有效位能的制造率在夏季最大,說(shuō)明夏季的潛熱釋放過(guò)程最為強(qiáng)烈;緯向平均動(dòng)能和渦動(dòng)動(dòng)能的耗散過(guò)程年變化特征相同,也是冬季高,夏季低,春秋季節(jié)處于過(guò)渡階段,這與大氣環(huán)流的強(qiáng)弱程度有關(guān)。

    2)月平均大氣能量循環(huán)的年平均結(jié)果表明,與能量間的轉(zhuǎn)化率相比,能量越赤道交換過(guò)程非常微弱。能量轉(zhuǎn)換過(guò)程一般是從緯向平均有效位能開(kāi)始,通過(guò)渦動(dòng),有效位能轉(zhuǎn)換成渦動(dòng)動(dòng)能,最后轉(zhuǎn)換為緯向平均動(dòng)能;50多年的平均結(jié)果說(shuō)明,Hadley環(huán)流的作用比Ferrel環(huán)流的作用更大,使得能量由緯向平均有效位能向緯向平均動(dòng)能轉(zhuǎn)化。

    3)年平均情況下,PM主要分布于高緯地區(qū),在35°N附近具有最小值,接近于零;PE在中高緯地區(qū)具有較大貯存,在赤道地區(qū)貯存最小;KM隨緯度的分布型與PM正好相反,即在35°N附近具有最大值,而在赤道和高緯地區(qū)具有最小值;KE隨緯度的分布與PE類似,在中緯和高緯地區(qū)的貯存較大。轉(zhuǎn)化項(xiàng)C(PM,PE)在中緯度地區(qū)最為活躍,轉(zhuǎn)化率達(dá)最大,在低緯地區(qū)這一轉(zhuǎn)化過(guò)程非常微弱;C(PE,KE)在整個(gè)北半球都很活躍,只有在極地附近減弱明顯;C(KE,KM)主要發(fā)生在中低緯地區(qū),高緯地區(qū)幾乎無(wú)轉(zhuǎn)化;C(PM,KM)反映的是三圈環(huán)流的貢獻(xiàn),主要取決于低緯直接環(huán)流和中緯間接環(huán)流的相對(duì)強(qiáng)度,至于高緯的直接環(huán)流非常微弱,貢獻(xiàn)很小。冬季和夏季的情況與此類似,主要在數(shù)值上有所差異。

    由于資料質(zhì)量的提高和資料時(shí)間長(zhǎng)度的增長(zhǎng),本文所估算的能量循環(huán)的年變化結(jié)果與前人的研究存在一定的差異。比如,Oort和Peixoto(1974)的結(jié)果中C(KE,KM)并沒(méi)有明顯的年變化;C(PM,KM)只有在秋季和初冬月份(10—12月)轉(zhuǎn)化率為正值,而在其他季節(jié)能量均由KM向轉(zhuǎn)化PM。另外,盡管D(KM)的年變化特征與Oort and Peixoto(1974)一致,但其在數(shù)值上差別較大,本文D(KM)始終為正值,而Oort and Peixoto(1974)中KM的耗散率在春季和夏季月份(4—9月)則為負(fù)值。由此可見(jiàn),隨著資料的不斷積累,未來(lái)對(duì)于大氣能量循環(huán)的再估算仍具有重要意義。

    要說(shuō)明的是,由能量循環(huán)年變化的研究結(jié)果表明,大氣環(huán)流的能量循環(huán)過(guò)程是非常復(fù)雜的。在半球際的能量交換過(guò)程中,是何種原因?qū)е铝薆(PM)較大,目前除了平均經(jīng)圈環(huán)流(主要是Hadley環(huán)流)的作用外,其他尚不清楚;另外,能量循環(huán)的年際、年代際變化如何,以及與全球增暖的關(guān)系如何,仍需作進(jìn)一步研究。

    參考文獻(xiàn)(References):

    高麗,李建平.2007.大氣能量有效性的研究進(jìn)展[J].地球科學(xué)進(jìn)展,22(5):486-494.Gao Li,Li Jianping.2007.Progress in the study of atmospheric energy efficiency[J].Advances in Earth Science,22(5):468-494.(in Chinese).

    Hu Q,Tawaye Y,Feng S.2004.Variations of the northern hemisphere atmospheric energetic:1948—2000[J].J Climate,17(10):1975-1986.

    Kalnay E,Kanamitsu M,Kistler R,et al.1996.The NCEP/NCAR 40-year reanalysis project[J].Bull Amer Meteor Soc,77(3):437-471.

    Kim Y H,Kim M K.2012.Examination of the global Lorenz energy cycle using MERRA and NCEP-reanalysis 2[J].Clim Dyn,40:1358.doi:10.1007/s00382-021-1358-4.

    Krueger A F,Winston J S,Haines D A.1965.Computation of atmospheric energy and its transformation for the Northern Hemisphere for a recent five-year period[J].Mon Wea Rev,93(4):227-238.

    李俊.1990.北半球冬季強(qiáng)厄爾尼諾期大尺度波動(dòng)能量收支的異常特征[J].南京氣象學(xué)院學(xué)報(bào),13(3):307-315.Li Jun.1990.Abnormal features of the large-scal wave energy budget during the Northern winter strong El Nino event[J].J Nanjing Inst Meteor,13(3):307-315.(in Chinese).

    Li L,Ingersoll A P,Jiang X,et al.2006.Variations in the mechanical energy cycle of the atmosphere[J].Geophys Res Lett,34.doi:10.1029/2007GL029985.

    Li L,Ingersoll A P,Jiang X,et al.2007.Lorenz energy cycle of the global atmosphere based on reanalysis datasets[J].Geophys Res Lett,34,L16813.doi:10.1029/2007GL029985.

    Lorenz E N.1955.Available potential energy and the maintenance of the general circulation[J].Tellus,7(2):157-167.

    Lorenz E N.1967.The nature and theory of the general circulation of the atmosphere[M].Geneva:World Meteorological Organization Publication.

    羅連升,楊修群.2003.從有效位能變化來(lái)分析El Nino的年代際變化[J].氣象科學(xué),23(1):1-11.Luo Liansheng,Yang Xiuqun.2003.Interdecadal change of El Nino as seen from variations of available potential energy[J].J Meteor Sci,23(1):1-11.(in Chinese).

    Marques C A F,Richa A,Corte-Real J,et al.2009.Global atmospheric energetics from NCEP-Reanalysis 2 and ECMWF-ERA40 Reanalysis[J].Int J Climatol,29(2):159-174.

    Oort A H.1983.Global atmospheric circulation statistics,1958—1973[M].Washington:Government Printing Office:180-226.

    Oort A H,Peixoto J P.1974.The annual cycle of the energetic of the atmosphere on a planetary scale[J].J Geophys Res,79(18):2705-2719.

    Peixoto J P,Oort A H.1974.The annual distribution of atmospheric energy on a planetary scale[J].J Geophys Res,79(15):2149-2159.

    Saltzman B.1957.Equations governing the energetic of the larger scales of atmospheric turbulence in the domain of wave number[J].J Meteor,14(6):513-523.

    Saltzman B.1970.Larger-scale atmospheric energetics in the wave-number domain[J].Rev Geophys Space Phys,8:289-302.

    Sheng J,Derome J.1991.An observational study of energy transfer between the seasonal mean flow and eddies[J].Tellus A,43(2):128-144.

    施能,朱乾根.1996.北半球大氣環(huán)流特征量的長(zhǎng)期趨勢(shì)及年代際變化[J].南京氣象學(xué)院學(xué)報(bào),19(3):283-289.Shi Neng,Zhu Qiangeng.1996.Secular trends and interdecadal variation of characteristic quantities of atmospheric circulations for Northern winter[J].J Nanjing Inst Meteor,19(3):283-289.(in Chinese).

    Smagorinsky J.1963.General circulation experiments with the primitive equationsⅠ:The basic experiment[J].Mon Wea Rev,91(3):99-164.

    Ulbrich U,Speth P.1991.The global energy cycle of stationary and transient waves:Result from ECWMF analyses[J].Meteor Atmos Phys,45:125-138.

    魏民,仇永炎.1997.全球不同緯度帶平均有效位能的季節(jié)急變[J].應(yīng)用氣象學(xué)報(bào),8(2):137-146.Wei Min,Qiu Yongyan.1997.The seasonal rapid changes of global potential energy at different latitudes[J].Quart J Appl Meteor,8(2):137-146.(in Chinese).

    Wiin-Nielsen A.1967.On the annual variation and spectral distribution of atmospheric energy[J].Tellus,19(4):540-559.

    Wiin-Nielsen A,Brown J A.1962.On diagnostic computations of atmospheric heat sources and sinks and the generation of available potential energy[R]//Proceedings of the international symposium on numerical weather prediction 1960.Tokyo:Meteorological Society of Japan:593-613.

    Wiin-Nielsen A,Chen T C.1993.Fundamentals of atmospheric energetics[M].New York:Oxford University Press:5-101.

    吳洪寶.1982.北半球500毫巴候平均環(huán)流動(dòng)能變化的動(dòng)力統(tǒng)計(jì)學(xué)分析[J].南京氣象學(xué)院學(xué)報(bào),(2):178-188.Wu Hongbao.1982.A dynamic-statistical analysis of the kinetic energy variations in the mean pentad circulation at 500 mb in the Northern Hemisphere[J].J Nanjing Inst Meteor,(2):178-188.(in Chinese).

    楊哲,管兆勇,蔡佳熙.2012.近30 a來(lái)北半球?qū)α鞒檀髿庠戮h(huán)流的渦度減弱現(xiàn)象[J].大氣科學(xué)學(xué)報(bào),35(6):702-711.Yang Zhe,Guan Zhaoyong,Cai Jiaxi.2012.Weakening trends of tropospheric enstrophy of monthly mean circulation in the Northern Hemisphere in the past three decades[J].Trans Atmos Sci,35(6):702-711.(in Chinese).

    袁建強(qiáng).1991.瞬變渦動(dòng)強(qiáng)迫對(duì)大氣環(huán)流的作用[J].南京氣象學(xué)院學(xué)報(bào),14(1):128-132.Yuan Jianqiang.1991.The role of the transient eddy forcing in the formation of the atmospheric circulation[J].J Nanjing Inst Meteor,14(1):128-132.(in Chinese).

    張韜,吳國(guó)雄,郭裕福.2006.GOALS模式中大氣能量循環(huán)的診斷分析與不同版本計(jì)算結(jié)果的比較研究[J].大氣科學(xué),30(1):38-45.Zhang Tao,Wu Guoxiong,Guo Yufu.2005.The comparison between two versions of the GOALS model on the atmospheric energy cycle diagnosis[J].J Atmos Sci,30(1):38-45.(in Chinese).

    (責(zé)任編輯:張福穎)

    Climatological features of atmospheric energy cycle in

    Northern Hemisphere

    LI Xin-xin1,2,GUAN Zhao-yong1,LI Ming-gang1

    (1.Key Laboratory of Meteorological Disaster(NUIST),Ministry of Education,Nanjing 210044,China;

    2.Yancheng Meteorological Bureau,Yancheng 224005,China)

    Abstract:According to the framework of energy cycle proposed by Lorenz in 1967,the atmospheric energy cycle is estimated on a monthly basis in the Northern Hemisphere using the NCEP/NCAR reanalysis daily data for the period of 1958—2011.The annual variations of energy reservoirs,conversions,exchanges across the equator,generations and dissipations have been investigated in the present paper.The results show that the energy cycle of the Northern Hemispheric atmosphere varies on annual time-scale evidently.Reservoirs and conversion rates change in the same phase,with values larger in winter,smaller in summer.In winter,the zonal mean energies and eddy kinetic energy are transported northward from the Southern Hemisphere into the Northern Hemisphere.However,in boreal summer,they are transported southward from the Northern Hemisphere into the Southern Hemisphere,whereas the exchange of eddy available potential energy at the equatorial boundary is in contrast.A maximum in generation of the zonal mean available potential energy takes place in the fall,and the generation of eddy available potential energy has a maximum in the summer.The dissipations of both zonal mean and eddy kinetic energy are stronger in winter but weaker in summer.In the annual mean case,the cross-equatorial fluxes are smaller than the energy conversion processes.In the meridional direction,the zonal mean available potential energy mainly appears in high-latitudes,the zonal mean kinetic energy mainly in mid-and low-latitudes,and the eddy energies mainly in mid-and high-latitudes.Further more,the conversion processes of atmospheric energies are more active in the middle latitude regions.

    Key words:energy cycle;available potential energy;kinetic energy;annual variation

    doi:10.13878/j.cnki.dqkxxb.20130306002

    中圖分類號(hào):

    文章編號(hào):1674-7097(2015)04-0531-09P434.5

    文獻(xiàn)標(biāo)志碼:A

    通信作者:管兆勇,博士,教授,博士生導(dǎo)師,研究方向?yàn)闅夂騽?dòng)力學(xué),guanzy@nuist.edu.cn.

    基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)目(41105056;41175062);江蘇省研究生科研創(chuàng)新計(jì)劃項(xiàng)目(CXZZ12-0485);江蘇省青藍(lán)工程創(chuàng)新團(tuán)隊(duì)及高校優(yōu)勢(shì)學(xué)科工程(PAPD)

    收稿日期:2013-03-06;改回日期:2013-04-25

    猜你喜歡
    動(dòng)能
    “動(dòng)能和勢(shì)能”知識(shí)鞏固
    新動(dòng)能,源自創(chuàng)新力
    金橋(2022年4期)2022-05-05 06:09:50
    斜面之上探動(dòng)能
    多種競(jìng)賽激發(fā)創(chuàng)新動(dòng)能
    新動(dòng)能,新聊城
    走向世界(2018年11期)2018-12-26 01:12:32
    高新項(xiàng)目,“領(lǐng)跑”動(dòng)能轉(zhuǎn)換
    走向世界(2018年11期)2018-12-26 01:12:32
    創(chuàng)新,動(dòng)能轉(zhuǎn)換的“加速度”
    金橋(2018年4期)2018-09-26 02:24:46
    為構(gòu)建開(kāi)放創(chuàng)新新高地增添動(dòng)能
    “金企對(duì)接”轉(zhuǎn)換舊動(dòng)能
    商周刊(2018年15期)2018-07-27 01:41:16
    金融開(kāi)放與動(dòng)能轉(zhuǎn)換
    商周刊(2018年15期)2018-07-27 01:41:16
    一级片'在线观看视频| 色网站视频免费| 亚洲av成人不卡在线观看播放网 | 自线自在国产av| a级毛片黄视频| 久久 成人 亚洲| 又粗又硬又长又爽又黄的视频| 哪个播放器可以免费观看大片| 亚洲精品aⅴ在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲第一av免费看| 国产爽快片一区二区三区| 热re99久久精品国产66热6| 女人被躁到高潮嗷嗷叫费观| 午夜免费观看性视频| 在线天堂中文资源库| 制服丝袜香蕉在线| 街头女战士在线观看网站| 99久久精品国产亚洲精品| 国产精品三级大全| 国产在视频线精品| av电影中文网址| 久久99热这里只频精品6学生| 日日摸夜夜添夜夜爱| 青春草亚洲视频在线观看| 久久国产亚洲av麻豆专区| 老司机影院成人| 久久久久久免费高清国产稀缺| 人人妻,人人澡人人爽秒播 | 美女主播在线视频| 欧美日韩亚洲国产一区二区在线观看 | 国产成人精品无人区| 国产精品 国内视频| 免费高清在线观看视频在线观看| 亚洲av国产av综合av卡| 少妇的丰满在线观看| 一级毛片电影观看| 亚洲精品日本国产第一区| 国产精品久久久久久精品古装| tube8黄色片| 丝瓜视频免费看黄片| 国产精品久久久久久人妻精品电影 | 久久精品aⅴ一区二区三区四区| 国产亚洲精品第一综合不卡| av电影中文网址| 中文字幕高清在线视频| 亚洲av电影在线观看一区二区三区| 男女高潮啪啪啪动态图| 国产男人的电影天堂91| 日本色播在线视频| 99香蕉大伊视频| 啦啦啦在线观看免费高清www| 国产成人av激情在线播放| 天天躁日日躁夜夜躁夜夜| 狂野欧美激情性bbbbbb| 国产乱来视频区| 97在线人人人人妻| 国产在线免费精品| 日韩av在线免费看完整版不卡| 亚洲欧美成人综合另类久久久| 欧美中文综合在线视频| 99久国产av精品国产电影| 免费看av在线观看网站| 五月天丁香电影| 巨乳人妻的诱惑在线观看| 免费在线观看完整版高清| 黄频高清免费视频| bbb黄色大片| 亚洲 欧美一区二区三区| 亚洲精品日本国产第一区| 熟妇人妻不卡中文字幕| 国产不卡av网站在线观看| 90打野战视频偷拍视频| 午夜av观看不卡| 在线观看三级黄色| 欧美乱码精品一区二区三区| 在线亚洲精品国产二区图片欧美| kizo精华| 视频在线观看一区二区三区| 高清av免费在线| 亚洲欧美一区二区三区久久| 亚洲精品国产区一区二| 免费少妇av软件| 九草在线视频观看| 欧美精品一区二区免费开放| 国产精品成人在线| 老司机在亚洲福利影院| 午夜福利乱码中文字幕| 久久精品国产a三级三级三级| 国产无遮挡羞羞视频在线观看| 亚洲,一卡二卡三卡| 国产成人av激情在线播放| 高清av免费在线| 成人三级做爰电影| 多毛熟女@视频| 毛片一级片免费看久久久久| 亚洲欧美中文字幕日韩二区| 国产精品久久久久久久久免| 精品人妻一区二区三区麻豆| 午夜福利乱码中文字幕| 欧美日韩成人在线一区二区| 卡戴珊不雅视频在线播放| 日韩制服丝袜自拍偷拍| av.在线天堂| 十分钟在线观看高清视频www| 无限看片的www在线观看| 波多野结衣av一区二区av| 亚洲欧美一区二区三区黑人| 色婷婷av一区二区三区视频| 国产 精品1| 久久国产精品男人的天堂亚洲| 夫妻性生交免费视频一级片| 少妇人妻久久综合中文| 久久97久久精品| 国产精品嫩草影院av在线观看| 亚洲精品在线美女| 欧美黄色片欧美黄色片| 汤姆久久久久久久影院中文字幕| 多毛熟女@视频| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产欧美日韩在线播放| 亚洲男人天堂网一区| 高清黄色对白视频在线免费看| 亚洲伊人久久精品综合| 丁香六月欧美| 国产精品蜜桃在线观看| 十八禁网站网址无遮挡| 久久性视频一级片| 两性夫妻黄色片| 午夜久久久在线观看| 美国免费a级毛片| 一本大道久久a久久精品| av卡一久久| 超碰97精品在线观看| 久久精品aⅴ一区二区三区四区| 在线看a的网站| 日本wwww免费看| 两个人看的免费小视频| av卡一久久| 国产精品蜜桃在线观看| 精品酒店卫生间| 久久精品亚洲av国产电影网| 国产亚洲最大av| 国产在视频线精品| 一级爰片在线观看| 免费高清在线观看日韩| 成人国产麻豆网| 中文字幕色久视频| svipshipincom国产片| 国产成人精品无人区| 99久久人妻综合| 一区二区三区激情视频| 久久久久久免费高清国产稀缺| 国产亚洲一区二区精品| 一区二区三区激情视频| 侵犯人妻中文字幕一二三四区| 一级黄片播放器| 狂野欧美激情性bbbbbb| 蜜桃国产av成人99| 人人妻,人人澡人人爽秒播 | 亚洲国产看品久久| 深夜精品福利| 一区在线观看完整版| 久久影院123| 人人妻人人添人人爽欧美一区卜| 亚洲欧洲日产国产| 久久综合国产亚洲精品| 国产不卡av网站在线观看| 国产一级毛片在线| 少妇人妻精品综合一区二区| 麻豆乱淫一区二区| 国产 精品1| 国精品久久久久久国模美| 精品久久久精品久久久| 无限看片的www在线观看| 亚洲精品成人av观看孕妇| 亚洲熟女毛片儿| 伊人久久国产一区二区| 极品少妇高潮喷水抽搐| h视频一区二区三区| 日本黄色日本黄色录像| 青草久久国产| av天堂久久9| 亚洲欧美精品综合一区二区三区| 日本91视频免费播放| www.av在线官网国产| 国产 一区精品| 国产高清国产精品国产三级| 国产成人欧美| 精品少妇黑人巨大在线播放| 波多野结衣一区麻豆| 国产精品一区二区在线不卡| 一边摸一边做爽爽视频免费| 99re6热这里在线精品视频| 久久99精品国语久久久| 久久人人爽av亚洲精品天堂| 麻豆乱淫一区二区| 老汉色av国产亚洲站长工具| av免费观看日本| 我的亚洲天堂| 天天躁夜夜躁狠狠躁躁| 90打野战视频偷拍视频| 精品国产一区二区久久| 亚洲精品国产av成人精品| av线在线观看网站| 亚洲第一区二区三区不卡| 国产成人免费观看mmmm| 国产在视频线精品| 午夜福利一区二区在线看| 久久99一区二区三区| 777米奇影视久久| 高清欧美精品videossex| 在线观看人妻少妇| 久久精品国产亚洲av高清一级| 免费久久久久久久精品成人欧美视频| 欧美日韩av久久| 国产精品一区二区精品视频观看| 久久鲁丝午夜福利片| 婷婷色综合大香蕉| 免费久久久久久久精品成人欧美视频| 久久综合国产亚洲精品| 精品久久久精品久久久| 亚洲国产精品成人久久小说| a 毛片基地| 国产精品久久久久久久久免| 亚洲伊人色综图| 日日啪夜夜爽| 久久国产精品男人的天堂亚洲| 亚洲国产精品999| 色网站视频免费| 91国产中文字幕| 日韩制服骚丝袜av| 欧美国产精品一级二级三级| 少妇人妻久久综合中文| 精品一区在线观看国产| 久久人人爽av亚洲精品天堂| 久久青草综合色| 丝袜美足系列| www日本在线高清视频| 国产精品久久久久久精品古装| 免费少妇av软件| 亚洲在久久综合| 制服丝袜香蕉在线| 成人手机av| 欧美日韩成人在线一区二区| 久久这里只有精品19| 午夜久久久在线观看| 亚洲成人手机| 18禁观看日本| 免费观看人在逋| 两个人看的免费小视频| 80岁老熟妇乱子伦牲交| 久久国产亚洲av麻豆专区| 老熟女久久久| 涩涩av久久男人的天堂| 欧美av亚洲av综合av国产av | 国产一卡二卡三卡精品 | 狂野欧美激情性xxxx| 男女高潮啪啪啪动态图| 王馨瑶露胸无遮挡在线观看| 天堂中文最新版在线下载| 90打野战视频偷拍视频| 午夜福利网站1000一区二区三区| 又大又黄又爽视频免费| 欧美日韩视频高清一区二区三区二| 曰老女人黄片| 国产亚洲一区二区精品| 美女中出高潮动态图| 中文字幕另类日韩欧美亚洲嫩草| 99精品久久久久人妻精品| 最新在线观看一区二区三区 | 人妻 亚洲 视频| 一区二区日韩欧美中文字幕| 18禁国产床啪视频网站| 国产免费又黄又爽又色| 1024香蕉在线观看| 国产一区二区激情短视频 | 中文精品一卡2卡3卡4更新| 欧美中文综合在线视频| 国产99久久九九免费精品| 色播在线永久视频| 在线 av 中文字幕| 中文字幕亚洲精品专区| 天天躁狠狠躁夜夜躁狠狠躁| 欧美黑人精品巨大| 美女高潮到喷水免费观看| 国产精品久久久久久人妻精品电影 | 成人三级做爰电影| 两个人看的免费小视频| 久久青草综合色| a级毛片黄视频| 晚上一个人看的免费电影| 免费黄网站久久成人精品| 男女无遮挡免费网站观看| 90打野战视频偷拍视频| 日韩一区二区视频免费看| 国产亚洲av片在线观看秒播厂| 国产精品国产av在线观看| 久久鲁丝午夜福利片| 久久久久精品久久久久真实原创| 日韩精品有码人妻一区| 久久这里只有精品19| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av电影在线观看一区二区三区| 毛片一级片免费看久久久久| 下体分泌物呈黄色| 乱人伦中国视频| 亚洲国产看品久久| 国产精品二区激情视频| 久久99一区二区三区| 日韩人妻精品一区2区三区| 国产成人91sexporn| 天美传媒精品一区二区| 午夜激情av网站| 亚洲少妇的诱惑av| 午夜福利,免费看| 国产欧美亚洲国产| 久久99一区二区三区| 一区二区av电影网| 伦理电影免费视频| 一本大道久久a久久精品| 在线观看www视频免费| 国产成人系列免费观看| 五月开心婷婷网| 一级片免费观看大全| 色综合欧美亚洲国产小说| 亚洲欧美一区二区三区久久| 国产精品一区二区精品视频观看| 国产亚洲av片在线观看秒播厂| 欧美日本中文国产一区发布| 国产黄色免费在线视频| 99久久精品国产亚洲精品| 久久久久久久久免费视频了| 下体分泌物呈黄色| 一个人免费看片子| 欧美日本中文国产一区发布| 在线观看免费日韩欧美大片| 七月丁香在线播放| 日韩大片免费观看网站| 亚洲情色 制服丝袜| 老司机在亚洲福利影院| 国产一区有黄有色的免费视频| 亚洲在久久综合| 欧美精品人与动牲交sv欧美| 亚洲伊人久久精品综合| 搡老乐熟女国产| 在线 av 中文字幕| 精品亚洲乱码少妇综合久久| 熟女av电影| 夜夜骑夜夜射夜夜干| 国产片特级美女逼逼视频| 新久久久久国产一级毛片| 国产av一区二区精品久久| 成人亚洲欧美一区二区av| 久久99精品国语久久久| 中文字幕色久视频| 91老司机精品| 日韩一区二区三区影片| 青春草视频在线免费观看| 欧美日韩亚洲综合一区二区三区_| 国产高清不卡午夜福利| 国产免费现黄频在线看| 亚洲av国产av综合av卡| 婷婷色综合大香蕉| 久久久久视频综合| 一二三四中文在线观看免费高清| 色婷婷久久久亚洲欧美| 美女视频免费永久观看网站| 国产精品久久久人人做人人爽| 国产精品国产三级专区第一集| 亚洲精品国产一区二区精华液| 少妇人妻精品综合一区二区| 亚洲男人天堂网一区| 久久99精品国语久久久| 国产精品欧美亚洲77777| 99久国产av精品国产电影| 午夜福利乱码中文字幕| 欧美亚洲 丝袜 人妻 在线| 校园人妻丝袜中文字幕| 中文字幕亚洲精品专区| 国产有黄有色有爽视频| 国产精品人妻久久久影院| 美国免费a级毛片| 一个人免费看片子| 亚洲精品久久久久久婷婷小说| 亚洲国产精品一区三区| 国产爽快片一区二区三区| 国产成人精品无人区| 男女高潮啪啪啪动态图| 亚洲精品中文字幕在线视频| 久久精品人人爽人人爽视色| 国产一区二区在线观看av| 国语对白做爰xxxⅹ性视频网站| 午夜影院在线不卡| 日韩av在线免费看完整版不卡| 久久久久人妻精品一区果冻| 在线免费观看不下载黄p国产| 国产毛片在线视频| 少妇人妻 视频| 2021少妇久久久久久久久久久| 欧美变态另类bdsm刘玥| 免费人妻精品一区二区三区视频| 人人澡人人妻人| 男人舔女人的私密视频| 精品免费久久久久久久清纯 | 80岁老熟妇乱子伦牲交| 婷婷色综合大香蕉| 国产精品久久久久久精品电影小说| 国产午夜精品一二区理论片| 天天躁夜夜躁狠狠久久av| 久久久久久久久免费视频了| 中文字幕亚洲精品专区| 精品国产露脸久久av麻豆| 国产亚洲午夜精品一区二区久久| 国产精品久久久久久精品古装| 久久av网站| e午夜精品久久久久久久| 国产精品熟女久久久久浪| 久久久精品免费免费高清| 汤姆久久久久久久影院中文字幕| 一区二区三区四区激情视频| 国产亚洲精品第一综合不卡| 亚洲av欧美aⅴ国产| 久久热在线av| 欧美日韩亚洲高清精品| 精品国产一区二区三区久久久樱花| 成人三级做爰电影| 另类精品久久| 国产有黄有色有爽视频| 男女国产视频网站| 精品亚洲乱码少妇综合久久| 亚洲av国产av综合av卡| 亚洲 欧美一区二区三区| 国产精品麻豆人妻色哟哟久久| 国产极品天堂在线| 亚洲国产欧美网| 波多野结衣av一区二区av| 十八禁人妻一区二区| 黄色怎么调成土黄色| 巨乳人妻的诱惑在线观看| 精品少妇黑人巨大在线播放| 女人精品久久久久毛片| 在线观看www视频免费| 国产伦人伦偷精品视频| 丝袜人妻中文字幕| 亚洲免费av在线视频| 久久久久网色| 女人久久www免费人成看片| 大片电影免费在线观看免费| 天堂俺去俺来也www色官网| 波野结衣二区三区在线| 亚洲欧洲精品一区二区精品久久久 | 国产 精品1| 19禁男女啪啪无遮挡网站| 一级片'在线观看视频| 9热在线视频观看99| 国产高清不卡午夜福利| 国产精品久久久久久久久免| 欧美日韩亚洲综合一区二区三区_| 国产精品 欧美亚洲| 18禁观看日本| 欧美日韩综合久久久久久| 久久综合国产亚洲精品| 亚洲第一av免费看| 亚洲av国产av综合av卡| 欧美乱码精品一区二区三区| 亚洲人成网站在线观看播放| av片东京热男人的天堂| 日韩中文字幕欧美一区二区 | 老汉色av国产亚洲站长工具| 久久久国产一区二区| 中文字幕制服av| 交换朋友夫妻互换小说| 狂野欧美激情性xxxx| 999精品在线视频| 一级,二级,三级黄色视频| 久久 成人 亚洲| 亚洲国产精品国产精品| 亚洲国产精品一区二区三区在线| 一区二区av电影网| 久久热在线av| 国产日韩欧美在线精品| 亚洲美女视频黄频| 丝袜脚勾引网站| 五月开心婷婷网| 青春草亚洲视频在线观看| 国产精品国产三级国产专区5o| 免费观看av网站的网址| 精品卡一卡二卡四卡免费| 欧美激情 高清一区二区三区| 香蕉国产在线看| 国产黄色视频一区二区在线观看| 久久久国产精品麻豆| 亚洲欧洲日产国产| 国产日韩欧美在线精品| 少妇人妻 视频| 蜜桃国产av成人99| 免费高清在线观看视频在线观看| 在线 av 中文字幕| 国产av精品麻豆| 精品卡一卡二卡四卡免费| 精品少妇黑人巨大在线播放| 成人三级做爰电影| 中文字幕高清在线视频| 亚洲中文av在线| 我的亚洲天堂| 国产精品久久久人人做人人爽| av网站在线播放免费| 国产精品麻豆人妻色哟哟久久| 日本vs欧美在线观看视频| 成年av动漫网址| 在线观看免费日韩欧美大片| 久久久久国产精品人妻一区二区| 纵有疾风起免费观看全集完整版| 人妻人人澡人人爽人人| 免费高清在线观看日韩| 久久午夜综合久久蜜桃| 涩涩av久久男人的天堂| 日韩精品有码人妻一区| 晚上一个人看的免费电影| 免费人妻精品一区二区三区视频| 国产精品99久久99久久久不卡 | 精品少妇内射三级| 亚洲美女搞黄在线观看| 在线观看免费日韩欧美大片| 黄色视频在线播放观看不卡| 男女之事视频高清在线观看 | 少妇的丰满在线观看| 新久久久久国产一级毛片| 欧美日韩福利视频一区二区| 国产日韩欧美视频二区| 国产精品 国内视频| 老司机影院成人| 欧美日韩一区二区视频在线观看视频在线| 国产高清国产精品国产三级| 久久亚洲国产成人精品v| 制服人妻中文乱码| 成年av动漫网址| 久久免费观看电影| 一二三四在线观看免费中文在| 你懂的网址亚洲精品在线观看| 欧美乱码精品一区二区三区| 国产免费现黄频在线看| 欧美日韩国产mv在线观看视频| 日本黄色日本黄色录像| 久久ye,这里只有精品| 日韩精品免费视频一区二区三区| 精品久久蜜臀av无| 新久久久久国产一级毛片| 中国三级夫妇交换| 在线观看人妻少妇| 日本vs欧美在线观看视频| 中文字幕亚洲精品专区| 久久精品久久久久久噜噜老黄| 日韩制服丝袜自拍偷拍| 男人操女人黄网站| 欧美 日韩 精品 国产| 哪个播放器可以免费观看大片| 亚洲欧美激情在线| 麻豆av在线久日| 在线免费观看不下载黄p国产| 国产一区有黄有色的免费视频| 亚洲欧洲精品一区二区精品久久久 | 日日啪夜夜爽| 成人国产麻豆网| 99久久人妻综合| av片东京热男人的天堂| 韩国高清视频一区二区三区| 免费黄频网站在线观看国产| 两个人免费观看高清视频| 成人亚洲精品一区在线观看| 操美女的视频在线观看| 国产精品久久久久久久久免| 岛国毛片在线播放| 中文字幕色久视频| a级毛片黄视频| 日韩熟女老妇一区二区性免费视频| 一级片'在线观看视频| av线在线观看网站| 哪个播放器可以免费观看大片| 水蜜桃什么品种好| 51午夜福利影视在线观看| 观看美女的网站| 免费看不卡的av| 久久久国产一区二区| 亚洲第一av免费看| 最近最新中文字幕免费大全7| 精品人妻一区二区三区麻豆| 日本wwww免费看| 久久精品国产a三级三级三级| 18在线观看网站| 色婷婷久久久亚洲欧美| 久久精品人人爽人人爽视色| 久久人人爽av亚洲精品天堂| 成人午夜精彩视频在线观看| 国产男女内射视频| 国产乱人偷精品视频| 久久久久精品久久久久真实原创| 韩国高清视频一区二区三区| 日韩av在线免费看完整版不卡| 中文字幕精品免费在线观看视频| av有码第一页| 精品国产国语对白av| 久久精品人人爽人人爽视色| 久久99热这里只频精品6学生| 尾随美女入室| 中国国产av一级| 日韩欧美一区视频在线观看| 成年女人毛片免费观看观看9 | 国产老妇伦熟女老妇高清| 亚洲欧美一区二区三区黑人| 久久久久精品性色| 美女脱内裤让男人舔精品视频| 美女福利国产在线| 岛国毛片在线播放|