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

    青藏高原對(duì)流層頂高度與臭氧總量及上升運(yùn)動(dòng)的耦合關(guān)系

    2012-09-11 02:31:40楊雙艷周順武張人禾吳萍李慧馬振鋒
    大氣科學(xué)學(xué)報(bào) 2012年4期
    關(guān)鍵詞:季節(jié)變化青藏高原

    楊雙艷,周順武,張人禾,吳萍,李慧,馬振鋒

    (1.南京信息工程大學(xué)氣象災(zāi)害省部共建教育部重點(diǎn)實(shí)驗(yàn)室,江蘇南京210044; 2.中國(guó)氣象科學(xué)研究院災(zāi)害天氣國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京100081;3.四川省氣候中心,四川成都610071)

    青藏高原對(duì)流層頂高度與臭氧總量及上升運(yùn)動(dòng)的耦合關(guān)系

    楊雙艷1,2,周順武1,張人禾2,吳萍1,李慧1,馬振鋒3

    (1.南京信息工程大學(xué)氣象災(zāi)害省部共建教育部重點(diǎn)實(shí)驗(yàn)室,江蘇南京210044; 2.中國(guó)氣象科學(xué)研究院災(zāi)害天氣國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京100081;3.四川省氣候中心,四川成都610071)

    根據(jù)1979—2008年青藏高原地區(qū)14個(gè)探空站對(duì)流層頂氣壓資料以及同期各標(biāo)準(zhǔn)等壓面上的溫度資料,分析了不同季節(jié)高原上空兩類對(duì)流層頂高度與高空各層溫度之間的關(guān)系;在此基礎(chǔ)上,結(jié)合同期的NCEP/NCAR月平均再分析資料以及NASA提供的TOMS/SBUV月平均臭氧總量資料,分別討論了高原上升運(yùn)動(dòng)以及高原臭氧總量與對(duì)流層頂高度的耦合關(guān)系。結(jié)果表明:高原第一(二)對(duì)流層頂高度全年處在300~200 hPa(100 hPa附近)高度,在季節(jié)變化、年際變化以及長(zhǎng)期變化趨勢(shì)上,兩類對(duì)流層頂高度與各自對(duì)應(yīng)高度層上的溫度存在著密切的反相變化關(guān)系,當(dāng)對(duì)流層頂高度偏高(低)時(shí),相應(yīng)高度上的溫度偏低(高)。上升運(yùn)動(dòng)有助于兩類對(duì)流層頂高度的抬升,尤其是當(dāng)高空200(100)hPa附近有上升運(yùn)動(dòng)時(shí),有利于第一(二)對(duì)流層頂高度抬升。各季節(jié)高原臭氧總量與第二對(duì)流層頂高度均呈顯著的負(fù)相關(guān)關(guān)系,當(dāng)臭氧含量減少(增加)時(shí),該對(duì)流層頂高度將偏高(偏低),近年來(lái)伴隨著高原臭氧總量的減少,高原第二對(duì)流層頂高度出現(xiàn)了明顯的抬升。關(guān)鍵詞:青藏高原;對(duì)流層頂高度;季節(jié)變化;上升運(yùn)動(dòng);臭氧總量

    Coupling relationship between tropopause height and total ozone as well as ascending motion over the Tibetan Plateau

    Abstract:Using tropopause pressure data and upper temperature data of 14 sounding stations over the Tibetan Plateau from 1979 to 2008,the relationships between two types of tropopause heights and upper air temperature in different seasons are analyzed.Based on these,the coupling connection between ascending motion(total ozone)and tropopause height is discussed by monthly NCEP/NCAR reanalysis data from NOAA and monthly TOMS/SBUV total ozone data from NASA.The results indicate that the first(second) tropopause is at 300—200 hPa(near 100 hPa)over the year.The seasonal,interannual variations and long-term trend of the two types of tropopause heights display closely opposite relation to those of temperature at their corresponding altitudes,respectively,namely the higher(lower)tropopause height,the lower (higher)temperature at the corresponding altitude.The ascending motion is conducive to the rising of thetropopause height,especially the ascending motion at about 200 hPa(100 hPa)is helpful to the rising of the first(second)tropopause height.There is an obvious negative correlation between the total ozone and the second tropopause height at each season.The second tropopause rises(reduces)with the losing(increasing)of total ozone over the Tibetan Plateau.In recent years,the second tropopause height clearly lifts with the reducing of the total ozone.

    Key words:Tibetan Plateau;tropopause height;seasonal variation;ascending motion;total ozone

    0 引言

    隨著對(duì)氣候變化研究的深入,近年來(lái)人們已將氣候變化的關(guān)注點(diǎn)延伸至整個(gè)自由大氣,其中關(guān)于對(duì)流層頂高度變化已有不少研究成果,我國(guó)氣象工作者基于探空資料對(duì)我國(guó)不同區(qū)域的對(duì)流層頂高度演變特征進(jìn)行了分析(吳香玲,1995;張廣興等,2005;陳芳等,2007;李輯等,2009;楊雙艷等,2010)。周順武等(2010)根據(jù)青藏高原14個(gè)探空站近30 a (1979—2008年)的對(duì)流層頂觀測(cè)資料,分析了青藏高原(簡(jiǎn)稱高原)上空兩類對(duì)流層頂高度的季節(jié)變化特征,發(fā)現(xiàn)高原上空全年均可觀測(cè)到第二對(duì)流層頂,其中在暖季(6—10月)第二對(duì)流層頂占絕對(duì)主導(dǎo)地位;兩類對(duì)流層頂高度在季節(jié)變化上存在著明顯的差異,第一(二)對(duì)流層頂高度的年變化曲線呈雙(單)峰型;近30 a來(lái)高原第一(二)對(duì)流層頂高度在不同季節(jié)存在著不同程度的降低(上升)。

    影響對(duì)流層頂變化的因子主要有太陽(yáng)輻射(地理緯度)、大氣環(huán)流、下墊面性質(zhì)與大地形等(Reid and Gage,1981;Maxobep,1987;鄒進(jìn)上等,1989; Shimizu and Tsuda,2000)。朱保林(2005)通過(guò)對(duì)高原對(duì)流層頂高度異常年份進(jìn)行個(gè)例分析后發(fā)現(xiàn),制約高原對(duì)流層頂高度(氣壓)變動(dòng)的因素主要有對(duì)流層和平流層中溫度異常、整層大氣加熱異常等。李鵬(2007)研究認(rèn)為,高原臭氧低值與對(duì)流層頂存在著密切的關(guān)系。

    本文基于高原地區(qū)14個(gè)探空站近30 a的對(duì)流層頂觀測(cè)資料以及同期各標(biāo)準(zhǔn)等壓面的月平均溫度資料,在分析高原高空各層溫度與兩類對(duì)流層頂高度關(guān)系的基礎(chǔ)上,討論了高原地區(qū)對(duì)流層頂高度與上升運(yùn)動(dòng)以及臭氧總量的耦合關(guān)系。

    1 資料

    本文選取由中國(guó)氣象局信息中心提供的高原地區(qū)(包括西藏、青海以及周邊高海拔山區(qū))14個(gè)探空站近30 a(1979—2008年)觀測(cè)的逐日對(duì)流層頂氣壓資料,通過(guò)壓高公式換算出高原月平均對(duì)流層頂高度(周順武等,2010)。同時(shí)還根據(jù)該中心提供的上述14個(gè)站各標(biāo)準(zhǔn)等壓面的月平均溫度資料(由于10 hPa的數(shù)據(jù)缺測(cè)嚴(yán)重,同時(shí)考慮到高原地區(qū)的海拔高度,文中只分析500~20 hPa之間11個(gè)標(biāo)準(zhǔn)等壓面溫度),將各站各層溫度作算術(shù)平均得到反映高原地區(qū)各層逐月的溫度序列(周順武和張人禾,2009)。文中用TH1(TH2)表示第一(第二)對(duì)流層頂高度。

    文中使用的資料還包括:1)NCEP/NCAR提供的同期(1968—1996年)月平均再分析資料(Kalany et al.,1996);2)NASA通過(guò)TOMS和SBUV相結(jié)合得到的最近30 a(1979—2008年)全球(10°經(jīng)度× 5°緯度)臭氧總量月平均資料(來(lái)自網(wǎng)站http:// code916.gsfc.nasa.gov/Data_services/merged/index.html)(Camp et al.,2003;Steinbrecht et al.,2003),該資料的優(yōu)點(diǎn)是時(shí)間連續(xù)性好,近年來(lái)TOMS/ SBUV臭氧總量資料已被廣泛用于相關(guān)的臭氧變化研究中(Hu and Tung,2002;Randel and Wu,2002; Zhou and Zhang,2005;Ziemke and Chandra,2005;韋惠紅和鄭有飛,2006;Zhang and Zhou,2009)。

    2 兩類對(duì)流層頂高度與溫度變化的關(guān)系

    2.1 兩類對(duì)流層頂高度與高空溫度在季節(jié)變化上的聯(lián)系

    圖1給出了高原上空各月多年平均溫度隨高度的分布??梢?jiàn),各月TH2(粗實(shí)線)在100 hPa高度附近波動(dòng),在此高度以下(上)溫度隨高度升高而下降(上升),高原上空最低溫度出現(xiàn)在第二對(duì)流層頂附近,其中6—9月還存在一個(gè)小于-70℃的低溫閉合中心,與TH2在春夏(秋冬)季偏高(低)的年變化特征(周順武等,2010)相反,150~50 hPa之間的溫度表現(xiàn)出夏低冬高的變化特征。而TH1(粗虛線)全年處在300~200 hPa高度之間,對(duì)應(yīng)的溫度年變化曲線與TH1年變化曲線(周順武等,2010)基本相反,表現(xiàn)出夏季偏高、春秋季偏低的變化特征。由此可見(jiàn),兩類對(duì)流層頂高度均與其所在高度層上的溫度存在著相反的季節(jié)變化;對(duì)應(yīng)高度層上的溫度偏高時(shí),對(duì)流層頂高度偏低,反之亦然。

    圖1 高原上空各月多年平均溫度隨高度的分布(圖中粗虛線和粗實(shí)線分別表示第一和第二對(duì)流層頂高度)Fig.1The variation of annual average monthly temperature with height over the Tibetan Plateau(Heavy dashed and solid lines indicate the first and second tropopause height,respectively)

    2.2 兩類對(duì)流層頂高度與高空溫度在年際和長(zhǎng)期變化趨勢(shì)上的聯(lián)系

    由于高原上空兩類對(duì)流層頂高度存在明顯的季節(jié)差異,所以分別計(jì)算了近30 a TH1和TH2與高空各層溫度之間的逐月相關(guān)關(guān)系(圖2),以便更清楚地理解它們隨高度的變化關(guān)系。由圖2a可見(jiàn),除了在TH1出現(xiàn)頻率最低的夏季外,其余季節(jié)TH1與300~150 hPa各層的溫度存在明顯的負(fù)相關(guān)關(guān)系(通過(guò)0.05信度的顯著性檢驗(yàn)),特別是與TH1平均高度——300~200 hPa各層溫度之間的負(fù)相關(guān)關(guān)系顯著(通過(guò)0.01信度的顯著性檢驗(yàn));而與100 hPa以上和400 hPa以下各層溫度則普遍為正相關(guān),特別是與春、秋季100~50 hPa高度的溫度的正相關(guān)關(guān)系通過(guò)了0.05信度的顯著性檢驗(yàn)。這表明TH1與該對(duì)流層頂高度附近的溫度存在反相變化關(guān)系,即相應(yīng)高度上的溫度偏低(高)則TH1偏高(低),同時(shí)與平流層下部溫度存在一定的正相關(guān)關(guān)系,其中在春、秋季兩者的正相關(guān)關(guān)系顯著。

    由圖2b可見(jiàn),各月相關(guān)系數(shù)為0的高度位于200~150 hPa附近,而在其之下(上)則存在較明顯的正(負(fù))相關(guān)關(guān)系。4—11月,TH2與對(duì)流層中上層(150 hPa以下)各層溫度之間存在明顯的正相關(guān)關(guān)系;而TH2與平流層下部各月溫度之間均為顯著的負(fù)相關(guān)關(guān)系,其中與100~50 hPa高度上溫度的相關(guān)系數(shù)普遍通過(guò)了0.01信度的顯著性檢驗(yàn)。由此可見(jiàn),TH2與相應(yīng)高度上的溫度之間存在反相變化關(guān)系;該對(duì)流層頂高度附近的溫度偏低(高)時(shí),TH2則偏高(低)。

    圖2 高原上空各層溫度與第一對(duì)流層頂高度(a)和第二對(duì)流層頂高度(b)的逐月相關(guān)系數(shù)分布(淺色和深色陰影區(qū)分別表示通過(guò)0.05和0.01信度的顯著性檢驗(yàn))Fig.2Correlation coefficients between(a)the first and(b)second tropopause heights and monthly temperature over the Tibetan Plateau(Light and heavy shaded regions with correlation exceeding 0.05 and 0.01 significance levels are shown,respectively)

    由此可見(jiàn),TH1與其平均高度所在的200~300 hPa的溫度存在密切聯(lián)系,而TH2與100~50 hPa溫度的負(fù)相關(guān)非常顯著,這與朱保林(2005)對(duì)高原對(duì)流層頂高度典型異常年份分析后得到的結(jié)論是一致的。

    為進(jìn)一步分析兩類對(duì)流層頂高度與相應(yīng)高度的溫度的逐年變化關(guān)系,圖3分別給出了近30a春季TH1與同期200 hPa溫度以及夏季TH2與同期50 hPa溫度的逐年變化曲線??梢?jiàn),對(duì)流層頂高度變化(實(shí)線,左縱坐標(biāo))與對(duì)應(yīng)高度層的溫度變化(虛線,表示溫度的右縱坐標(biāo)被顛倒)幾乎相反,二者存在顯著的負(fù)相關(guān)關(guān)系,其中春季TH1與200 hPa溫度的相關(guān)系數(shù)高達(dá)-0.913,夏季TH2與50 hPa溫度的相關(guān)系數(shù)也達(dá)到-0.866。顯著的負(fù)相關(guān)關(guān)系表明,對(duì)應(yīng)高度層的溫度升高(降低),則對(duì)流層頂高度降低(升高)。此外,由圖3中的線性變化趨勢(shì)線也可以看出,近30 a來(lái)春季TH1表現(xiàn)出明顯的下降趨勢(shì)(-110 m/(10 a)),而200 hPa溫度存在上升趨勢(shì)(0.388℃/(10 a));夏季TH2存在上升趨勢(shì)(59 m/(10 a)),而50 hPa溫度存在明顯的下降趨勢(shì)(-0.860℃/(10 a))。

    綜上所述,在季節(jié)、年際和長(zhǎng)期變化趨勢(shì)上,兩類對(duì)流層頂高度與各自對(duì)應(yīng)高度的溫度存在密切的聯(lián)系,當(dāng)對(duì)流層頂高度偏高時(shí),相應(yīng)高度上的溫度偏低,反之亦然。

    3 高原上空臭氧總量變化與對(duì)流層頂高度變化的聯(lián)系

    大氣臭氧主要產(chǎn)生于平流層,作為一種重要的溫室氣體,其含量變化將改變大氣垂直輻射平衡,從而導(dǎo)致氣候變化。自20世紀(jì)80年代發(fā)現(xiàn)南極臭氧洞(Farman et al.,1985)以來(lái),大氣臭氧變化與對(duì)流層頂高度變化之間的關(guān)系一直受到關(guān)注(王衛(wèi)國(guó)和秦芳,1994;Bojkov et al.,1995;Schubert and Munteanu,1998;李國(guó)輝等,2003;楊建和呂達(dá)仁,2003;Santer et al.,2003a,2003b)。Steinbrecht et al.(1998)分析了德國(guó)南部探空站對(duì)流層頂高度和臭氧總量的關(guān)系,發(fā)現(xiàn)當(dāng)臭氧總量偏低(高)時(shí),該站對(duì)流層頂高度偏高(低)。Chakrabarty et al.(2000)分析了印度兩個(gè)探空站的對(duì)流層頂高度和溫度資料后發(fā)現(xiàn),近幾十年來(lái)該地區(qū)熱帶對(duì)流層頂高度(溫度)的上升(下降)趨勢(shì)與平流層臭氧減少有關(guān)。春、夏季高原存在異常的大氣臭氧低值,且臭氧虧損最大值的高度位于對(duì)流層頂附近(鄒捍等,1998)。因此,臭氧含量的變化也是影響高原對(duì)流層頂高度變化的因子之一(李鵬,2007)。

    上述分析表明,TH2的變化與對(duì)流層和平流層溫度存在密切聯(lián)系,因此對(duì)流層頂?shù)慕Y(jié)構(gòu)與平流層下部的臭氧分布存在著密切聯(lián)系。Zhou and Zhang (2005)和Zhang and Zhou(2009)在對(duì)近24 a(1979—2002年)高原上空溫度變化特征進(jìn)行分析時(shí)發(fā)現(xiàn),在年際和年代際尺度上,高原上空各層溫度和臭氧總量變化之間均存在密切聯(lián)系,平流層低層和對(duì)流層上層(對(duì)流層中下層)的溫度與臭氧總量之間具有明顯(一定)的正(負(fù))相關(guān)關(guān)系。伴隨臭氧總量的減少,高原平流層低層(對(duì)流層中上層)的溫度出現(xiàn)了下降(上升)。這種由于高原臭氧總量改變導(dǎo)致的對(duì)流層中上層增溫和平流層下層降溫,勢(shì)必影響到高原對(duì)流層頂高度的變化。

    圖3 近30 a春季第一對(duì)流層頂高度與同期200 hPa溫度(a)以及夏季第二對(duì)流層頂高度與同期50 hPa溫度(b)的年際變化(左縱坐標(biāo)表示對(duì)流層頂高度,單位:km;右縱坐標(biāo)表示溫度,為了便于比較,表示溫度的左坐標(biāo)被顛倒,單位:℃)Fig.3Interannual variations of(a)the first tropopause height and 200 hPa temperature in spring,and(b)the second tropopause height and 50 hPa temperature in summer in recent 30 years(the left coordinate indicates the tropopause height,units:km; the right coordinate by reverse indicates temperature,units:℃)

    那么TH2的變化與臭氧總量變化之間又存在著怎樣的聯(lián)系呢?圖4給出了高原臭氧總量與TH2的逐月相關(guān)系數(shù)??梢?jiàn),高原臭氧總量與TH2之間普遍存在相反的變化關(guān)系,且大多數(shù)月份兩者之間的負(fù)相關(guān)系數(shù)通過(guò)了0.05信度的顯著性檢驗(yàn),特別是在春夏季,TH2與臭氧總量的相關(guān)關(guān)系更密切(通過(guò)了0.01信度的顯著性檢驗(yàn)),這表明當(dāng)高原臭氧總量減少(增加)時(shí),TH2將偏高(低)。

    圖4 高原臭氧總量與高原第二對(duì)流層頂高度逐月相關(guān)系數(shù)(短、長(zhǎng)虛線分別表示通過(guò)0.01、0.05信度的顯著性檢驗(yàn))Fig.4Monthly correlation coefficient between total ozone and the second tropopause height over the Tibetan Plateau(Short and long dashed lines indicate correlation coefficients exceeding 0.01 and 0.05 significance levels,respectively)

    為進(jìn)一步比較TH2與臭氧總量的逐年變化情況,圖5給出了近30 a夏季TH2與同期和前期(春季)高原臭氧總量的演變曲線。由圖5a可見(jiàn),近30 a來(lái)高原夏季臭氧總量大約減少6.5 DU(趨勢(shì)為-2.17 DU/(10 a)),而夏季TH2大約抬升180 m (為了便于比較,表示對(duì)流層頂高度的左坐標(biāo)被顛倒;趨勢(shì)為59 m/(10 a));TH2的變化與臭氧總量變化存在明顯的負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)達(dá)-0.678 (通過(guò)0.001信度的顯著性檢驗(yàn))。不僅如此,高原春季的臭氧總量變化(趨勢(shì)為-5.72 DU/(10 a))與夏季TH2變化之間關(guān)系密切(圖5b),相關(guān)系數(shù)達(dá)-0.477(通過(guò)0.01信度的顯著性檢驗(yàn))。由此可見(jiàn),高原臭氧總量虧損導(dǎo)致高原平流層下層明顯降溫,以及包括高原臭氧總量減少等多因素導(dǎo)致的高原對(duì)流層增溫,將改變高原上空溫度場(chǎng)的垂直結(jié)構(gòu)(Zhang and Zhou,2009),從而引起高原對(duì)流層頂高度的抬升。

    4 高原上空抬升作用與對(duì)流層頂高度的聯(lián)系

    大氣環(huán)流的變化勢(shì)必影響對(duì)流層頂高度的變化。由于青藏高原的熱力作用(Li et al.,2001)和動(dòng)力作用(黃榮輝,1985;Yanai et al.,1992;Ye and Wu,1998)對(duì)中國(guó)、東亞大氣環(huán)流和氣候變化都有極其重要的影響,因此高原本身的熱力和動(dòng)力作用造成的抬升作用對(duì)其對(duì)流層頂高度變化也必定存在著影響。那么高原上空的抬升運(yùn)動(dòng)作用對(duì)高原對(duì)流層頂高度存在怎樣的影響呢?

    圖5 夏季第二對(duì)流層頂(實(shí)線,左坐標(biāo)被顛倒,單位:km)與夏季(a)和春季(b)高原臭氧總量(虛線,右坐標(biāo),單位:DU)的年際變化(直線為線性趨勢(shì)線)Fig.5Interannual variaitions of TH2(solid curves;the left coordinate by reverse;units:km)and total ozone(dashed curves;right coordinate;units:DU)over the Tibetan Plateau in(a)summer and(b)spring(The straight line indicates the linear trend)

    鑒于再分析資料的垂直速度在高原地區(qū)存在較大的偏差(周任君和陳月娟,2005),同時(shí)考慮到高原上空對(duì)流層上部和平流層下部的水汽含量低,大氣運(yùn)動(dòng)可視為干絕熱運(yùn)動(dòng),位溫恒定不變。當(dāng)空氣上升時(shí),將下層較低的位溫帶到上層,使上層的位溫下降,局地的等位溫面抬升。因此,參考Zou and Gao(1997)的方法,采用位溫來(lái)表示高空垂直運(yùn)動(dòng),當(dāng)某層位溫明顯偏低時(shí),可視為有明顯的上升運(yùn)動(dòng)。

    以下根據(jù)兩類對(duì)流層頂高度在不同月份的異常偏高年和偏低年(表1),分別對(duì)高原上空各層位溫距平進(jìn)行合成分析。在5個(gè)TH1異常偏高年(圖6a),各月位溫在150 hPa以上(下)大致為正(負(fù))距平,尤其是冬春季300~200 hPa位溫表現(xiàn)為明顯的負(fù)距平,而在春秋季100~50 hPa位溫則表現(xiàn)出明顯的正距平;而在5個(gè)TH1異常偏低年(圖6b),位溫距平合成分布基本上與偏高年的位溫距平合成分布相反,即大致以150 hPa為界,低(高)層位溫為正(負(fù))距平。由TH1偏高年減去偏低年的位溫合成差值分布(圖6c)可見(jiàn),150 hPa以下(上)各月位溫為負(fù)(正)距平差值,除在TH1出現(xiàn)頻率最低的夏季外,其余季節(jié)300~150 hPa位溫負(fù)距平差值非常顯著,尤其在冬春季200 hPa附近,位溫表現(xiàn)為顯著的負(fù)距平差值(超過(guò)了0.01信度的顯著性檢驗(yàn)),100 hPa以上的各層位溫則為正距平差值,其中在春秋季70 hPa附近的正距平差值通過(guò)了0.05信度的顯著性檢驗(yàn)。這表明,在TH1偏高年與偏低年,高原高低層位溫距平分布相反,其中第一對(duì)流層頂所在高度的位溫差異最明顯,TH1偏高時(shí),該層位溫偏小,反之亦然。

    圖7為各月TH2偏高年和偏低年高原上空各層位溫距平合成分布,以及偏高年減去偏低年的各層位溫合成差值分布。在TH2偏高年(圖7a),各月(除9月表現(xiàn)不明顯外)大致在150 hPa以下(上)為正(負(fù))位溫距平,尤其是春夏季100~30 hPa正位溫距平明顯。在TH2偏低年(圖7b),各月位溫距平分布與偏高年基本相反,150 hPa以下低層為負(fù)位溫距平,而高層為正位溫距平,尤其是3—10月,100~30 hPa出現(xiàn)了明顯的正位溫距平,而春秋季200 hPa附近為明顯的位溫負(fù)距平。由圖7c可見(jiàn),150 hPa以上(下)為負(fù)(正)位溫距平差值,其中在平流層下層(100~30 hPa),各月負(fù)位溫距平差值顯著,而在對(duì)流層中上層(200 hPa以下),3—6月及秋季正位溫距平差值較明顯。由此可知,TH2偏高年與偏低年,高原上空平流層下層與對(duì)流層中上層的位溫距平分布存在反相變化特征,特別100~30 hPa位溫距平差值非常顯著。

    上述分析表明,兩類對(duì)流層頂高度異常年高原上空各層位溫變化存在明顯差異,特別是各自對(duì)流層頂所在高度附近的位溫存在顯著的變化特征(盡管在不同季節(jié)其密切程度存在差異),即相應(yīng)高度抬升明顯時(shí)位溫偏低,其對(duì)應(yīng)的對(duì)流層頂高度則偏高。為了進(jìn)一步比較春夏季兩類對(duì)流層頂高度與相應(yīng)高度位溫的關(guān)系,圖8給出了近30 a春季TH1與同期200 hPa位溫以及夏季TH2與同期100 hPa位溫的年際變化曲線。由圖8a可見(jiàn),TH1(實(shí)線)與位溫(虛線)的年際變化幾乎相反(表示位溫的右縱坐標(biāo)被顛倒),其線性相關(guān)系數(shù)達(dá)-0.835,表明兩者關(guān)系非常密切,與TH1的下降趨勢(shì)(趨勢(shì)為-110 m/(10 a))相反,200 hPa位溫呈現(xiàn)為上升趨勢(shì)(0.079 K/(10 a))。由圖8b可見(jiàn),TH2(實(shí)線)與位溫(虛線)存在密切的反相變化關(guān)系,其線性相關(guān)系數(shù)為-0.619(通過(guò)0.01信度的顯著性檢驗(yàn)),近30 a TH2上升近180 m(趨勢(shì)為59 m/(10 a)),而100 hPa位溫每10 a減少1.67 K。由此可知,兩類對(duì)流層頂與其所在高度附近的位溫表現(xiàn)出非常顯著的負(fù)相關(guān)關(guān)系,這意味著各自對(duì)流層頂高度上的位溫下降(即有上升運(yùn)動(dòng))時(shí),對(duì)流層頂高度抬高。

    表1 各月兩類對(duì)流層頂高度5個(gè)異常偏高年和5個(gè)異常偏低年的情況Table 1Five abnormal high years and five abnormal low years for the first and second tropopause heights in each month

    圖6 各月第一對(duì)流層頂高度偏高年(a)和偏低年(b)及偏高年減去偏低年(c)的高原上空各層的位溫距平合成分布(等值線為位溫距平,單位:K;圖中淺色和深色陰影區(qū)分別表示通過(guò)0.05和0.01信度的顯著性檢驗(yàn))Fig.6Composite monthly potential temperature anomalies in(a)the abnormal high years and(b)abnormal low years of the first tropopause height over the Tibetan Plateau,and(c)their differences(high years minus low years) (Isoline is potential temperature anomaly,units:K;Light and heavy shaded regions with difference exceeding 0.05 and 0.01 significance levels are shown,respectively)

    由此可以初步推斷,當(dāng)高原上空100 hPa位溫下降(即有上升運(yùn)動(dòng))時(shí),強(qiáng)烈的上升運(yùn)動(dòng)有利于對(duì)流層低濃度的臭氧向平流層輸送,從而導(dǎo)致臭氧總量減少(周秀驥等,1995;周秀驥和李維亮,2004;卞建春等,1997),而臭氧總量減少,又將改變高空溫度場(chǎng)的垂直結(jié)構(gòu),從而引起對(duì)流層頂高度的抬升。

    5 結(jié)論與討論

    在分析青藏高原兩類對(duì)流層頂高度季節(jié)變化特征的基礎(chǔ)上,本文分析了相應(yīng)高度上溫度的季節(jié)變化、年際變化和長(zhǎng)期趨勢(shì)變化與對(duì)流層頂高度變化的關(guān)系,并初步分析了高原上空抬升作用以及高原臭氧總量變化與高原對(duì)流層頂高度變化之間的耦合關(guān)系,得到以下結(jié)論:

    1)高原第一對(duì)流層頂高度全年處在300~200 hPa,第二對(duì)流層頂高度各月基本在100 hPa高度上下波動(dòng)。在季節(jié)變化、年際變化以及長(zhǎng)期變化趨勢(shì)上,兩類對(duì)流層頂高度與各自對(duì)應(yīng)高度上的溫度具有密切的反相變化關(guān)系;當(dāng)對(duì)流層頂高度偏高(低)時(shí),相應(yīng)高度上的溫度則偏低(高)。近30 a來(lái),高原上空平流層下層(對(duì)流層中層)存在明顯(一定)的降(增)溫,導(dǎo)致第一對(duì)流層頂高度降低,而第二對(duì)流層頂普遍升高。

    2)除在第一對(duì)流層頂出現(xiàn)頻率較低的夏季外,各季節(jié)該對(duì)流層頂高度偏高(低)與同期高原上空平流層下層增(降)溫和對(duì)流層中層降(增)溫有關(guān);各季節(jié)第二對(duì)流層頂高度偏高與平流層下層降溫以及對(duì)流層中層增溫存在密切關(guān)系,該對(duì)流層頂高度偏高年,對(duì)流層中層偏暖,整層氣柱伸長(zhǎng),特別是平流層下層偏冷,引起平流層氣柱收縮,都將使得該對(duì)

    圖7 各月第二對(duì)流層頂高度偏高年(a)和偏低年(b)及偏高年減去偏低年(c)的高原上空各層的位溫距平合成分布(等值線為位溫距平,單位:K;圖中淺色和深色陰影區(qū)分別表示通過(guò)0.05和0.01信度的顯著性檢驗(yàn))Fig.7Composite monthly potential temperature anomalies in(a)the abnormal high years and(b)abnormal low years of the second tropopause height over the Tibetan Plateau,and(c)their differences(high years minus low years) (Isoline is potential temperature anomaly,units:K;Light and heavy shaded regions with difference exceeding 0.05 and 0.01 significance levels are shown,respectively)

    圖8 近30 a春季第一對(duì)流層頂高度與同期200 hPa位溫(a)以及夏季第二對(duì)流層頂高度與同期100 hPa位溫(b)的年際變化(左縱坐標(biāo)表示對(duì)流層頂高度,單位:km;右縱坐標(biāo)表示位溫,為了便于比較,表示溫度的左坐標(biāo)被顛倒,單位: K)Fig.8Interannual variations of(a)the first tropopause height and 200 hPa potential temperature in spring,and(b)the second tropopause height and 100 hPa potential temperature in summer in recent 30 years(the left coordinate indicates the tropopause height,units:km;the right coordinate by reverse indicates potential temperature,units:K)

    流層頂上移,反之亦然。

    3)高原第二對(duì)流層頂高度與高原臭氧總量之間存在明顯的反相變化關(guān)系。當(dāng)臭氧總量減少(增加)時(shí),第二對(duì)流層頂高度升高(降低)。臭氧總量減少引起平流層下層溫度降低以及對(duì)流層中層溫度增高,改變了高原上空溫度場(chǎng)的垂直結(jié)構(gòu),從而引起第二對(duì)流層頂高度的抬升。

    4)對(duì)流層頂和平流層中下層的位溫變化可以反映該層的垂直運(yùn)動(dòng)。上升運(yùn)動(dòng)有助于高原上空兩類對(duì)流層頂高度抬升,尤其當(dāng)高空200 hPa附近有上升運(yùn)動(dòng)時(shí),有利于第一對(duì)流層頂高度的抬升;而當(dāng)100 hPa附近有上升運(yùn)動(dòng)時(shí),則有利于第二對(duì)流層頂高度的升高。

    文中討論了高原地區(qū)對(duì)流層頂與高空溫度、上升運(yùn)動(dòng)以及臭氧總量的關(guān)系,但是對(duì)流層頂與這些因子之間的聯(lián)系極有可能不是孤立的。臭氧含量可能是通過(guò)影響對(duì)流層與平流層的溫度結(jié)構(gòu),進(jìn)而對(duì)對(duì)流層頂產(chǎn)生影響,而上升運(yùn)動(dòng)可能是通過(guò)影響臭氧和溫度的分布來(lái)影響對(duì)流層頂高度的。當(dāng)然,對(duì)流層頂高度的變化也會(huì)影響到高空溫度和臭氧含量的變化(李國(guó)輝等,2003)。

    許多研究表明,影響對(duì)流層頂高度的因子眾多,如槽脊、低渦(鄒進(jìn)上等,1989),太陽(yáng)輻射(Reid and Gage,1981;鄒進(jìn)上等,1989;Shimizu and Tsuda,2000),云量(鄒進(jìn)上和江靜,1990;王旻燕和呂達(dá)仁,2007),急流、鋒(Maxobep,1987;鄒進(jìn)上和江靜,1990)以及地表溫度和海洋熱含量等(Sausen and Santer,2003)。本文僅分析了高原抬升作用以及臭氧總量變化與高原對(duì)流層頂高度的耦合關(guān)系,其分析本身還顯得比較簡(jiǎn)單,其中的具體物理機(jī)理尚待進(jìn)一步深入研究。

    致謝:中國(guó)國(guó)家氣象信息中心氣象資料室(Climatic DataCenter,NationalMeteorologicalInformation Center,CMA)提供了對(duì)流層頂探空資料和高空溫度資料,美國(guó)NASA提供臭氧總量觀測(cè)資料以及美國(guó)NOAA網(wǎng)站提供NCEP/NCAR再分析資料,在此一并表示感謝。

    卞建春,李維亮,周秀驥.1997.青藏高原及其鄰近地區(qū)流場(chǎng)結(jié)構(gòu)季節(jié)性變化的特征分析[M]//周秀驥.中國(guó)地區(qū)大氣臭氧變化及其對(duì)氣候環(huán)境的影響II.北京:氣象出版社:257-273.

    陳芳,馬英芳,朱西德.2007.青海省對(duì)流層頂若干統(tǒng)計(jì)特征[J].氣象科技,35(1):57-60.

    李國(guó)輝,呂達(dá)仁,Tie X X.2003.對(duì)流層頂變化對(duì)上對(duì)流層/下平流層臭氧分布的影響[J].空間科學(xué)學(xué)報(bào),23(4):269-277.

    李輯,蔡福,明惠青,等.2009.遼寧地區(qū)第一對(duì)流層頂高度變化特征分析[J].氣象與環(huán)境學(xué)報(bào),25(2):9-15.

    李鵬.2007.青藏高原大氣臭氧與對(duì)流層頂?shù)年P(guān)系[D].北京:中國(guó)科學(xué)院研究生院.

    王旻燕,呂達(dá)仁.2007.東亞強(qiáng)對(duì)流云的季節(jié)變化及其與對(duì)流層頂關(guān)系初探[J].大氣科學(xué),31(5):937-949.

    王衛(wèi)國(guó),秦芳.1994.臭氧層季節(jié)變化與對(duì)流層頂?shù)年P(guān)系[J].云南大學(xué)學(xué)報(bào):自然科學(xué)版,16(增刊1):29-34.

    韋惠紅,鄭有飛.2006.我國(guó)臭氧總量時(shí)空分布特征[J].南京氣象學(xué)院學(xué)報(bào),29(3):390-395.

    吳香玲.1995.北京地區(qū)極地對(duì)流層頂與地面要素之間的關(guān)系[J].氣象,21(11):42-45.

    楊健,呂達(dá)仁.2003.平流層—對(duì)流層交換研究進(jìn)展[J].地球科學(xué)進(jìn)展,18(3):380-385.

    假設(shè)Bob想給他人轉(zhuǎn)賬,Bob就得向全網(wǎng)廣播他要轉(zhuǎn)賬的消息,并需全網(wǎng)達(dá)成共識(shí),才能認(rèn)為他的消息是合法的,且每個(gè)節(jié)點(diǎn)都會(huì)保存他轉(zhuǎn)賬信息。全網(wǎng)沒(méi)有中心服務(wù)器,沒(méi)有人能擁有管理的權(quán)力,只要規(guī)則定好了,就必須照著規(guī)則做,沒(méi)有人可以改變,這其實(shí)就是區(qū)塊鏈去中心化的魅力所在。

    楊雙艷,周順武,陳鶴.2010.甘肅省對(duì)流層頂高度的季節(jié)變化特征分析[J].氣象,36(4):57-62.

    張廣興,李娟,崔彩霞,等.2005.新疆1960—1999年第一對(duì)流層頂高度變化及其突變分析[J].氣候變化研究進(jìn)展,1(3):106-110.

    周任君,陳月娟.2005.青藏高原和伊朗高原上空臭氧變化特征及其與南亞高壓的關(guān)系[J].中國(guó)科學(xué)技術(shù)大學(xué)學(xué)報(bào),35(6): 899-908.

    周順武,張人禾.2009.青藏高原地區(qū)上空NCEP/NCAR再分析溫度資料和位勢(shì)高度資料與觀察資料的比較分析[J].氣候與境研究,14(2):284-292.

    周順武,楊雙艷,張人禾,等.2010.青藏高原兩類對(duì)流層頂高度的季節(jié)變化特征分析[J].大氣科學(xué)學(xué)報(bào),33(3):307-314.

    周秀驥,李維亮.2004.青藏高原地區(qū)大氣臭氧變化的研究[J].氣象學(xué)報(bào),62(5):513-527.

    周秀驥,羅超,李維亮,等.1995.中國(guó)地區(qū)臭氧總量變化與青藏高原低值中心[J].科學(xué)通報(bào),40(15):1396-1398.

    朱保林.2005.北半球夏季對(duì)流層頂變動(dòng)的現(xiàn)象及機(jī)制研究[D].南京:南京信息工程大學(xué).

    鄒捍,郜永祺,周立波.1998.大尺度山地上空的臭氧低值及地面加熱[J].氣候與環(huán)境研究,3(3):209-217.

    鄒進(jìn)上,江靜.1990.高空氣候?qū)W[M].北京:氣象出版社.

    鄒進(jìn)上,張降秋,王炳忠.1989.我國(guó)對(duì)流層頂?shù)臅r(shí)空分布特征及其影響因子[J].氣象科學(xué),9(4):417-426.

    Bojkov R D,Bishop L,F(xiàn)ioletov V E.1995.Total ozone trends from quality-controlled ground-based data(1964—1994)[J].J Geophys Res,100:25867-25876.

    Camp C D,Roulston M S,Yung Y L.2003.Temporal and spatial patterns of the interannual variability of total ozone in the tropics[J].J Geophys Res,108(D20),4643.doi:10.1029/2001JD001504.

    Chakrabarty D K,Shah N C,Pandya K V.2000.Long-term trend of tropopause over New Delhi and Thiruvananthapuram[J].Geophys Res Lett.,27(15):2181-2184.

    Farman J C,Gardiner B G,Shanklin J D.1985.Large losses of total ozone in Antarctica reveal seasonal ClOx/NOx interaction[J].Nature,315:207-210.

    Hu Y Y,Tung K K.2002.Interannual and decadal variations of planetary wave activity,stratospheric cooling,and Northern Hemisphere annular mode[J].J Climate,15(13):1659-1673.

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

    Li Dongliang,Ji Guoliang,Lü Lanzhi.2001.Impact of Tibetan Plateau surface heating field intensity on Northern Hemispherical generalcirculations and weather and the climate of China[J].Science in China:Series D,44(3):390-399.

    Maxobep 3 M.1987.對(duì)流層頂氣候?qū)W[M].張貴銀,廖壽發(fā),譯.北京:氣象出版社.

    Randel W J,Wu F.2002.Changes in column ozone correlated with the stratospheric EP flux[J].J Meteor Soc Japan,80(4b):849-862.

    Reid G C,Gage K S.1981.On the annual variation in height of the tropical tropopause[J].J Atmos Sci,38(9):1928-1938.

    Santer B D,Wehner M F,Wigley T M L,et al.2003a.Contributions of anthropogenic and natural forcing to recent tropopause height changes[J].Science,301(7):479-483.

    Santer B D,Sausen R,Wigley T M L,et al.2003b.Behavior of tropopause height and atmospheric temperature in models,reanalyses,and observations:Decadal changes[J].J Geophys Res,108(D1),4002.doi:10.1029/2002JD002258.

    Sausen R,Santer B D.2003.Use of changes in tropopause height to detect human influences on climate[J].Meteorol Z,12(33): 131-136.

    Schubert S D,Munteanu M J.1998.An analysis of tropopause pressure and total ozone correlations[J].Mon Wea Rev,116(3):569-582.Shimizu A,Tsuda T.2000.Variations in tropical tropopause observed with radiosondes in Indonesia[J].Geophys Res Lett,27(16): 2541-2544.

    Steinbrehct W,Cluade H,Kohler U,et al.1998.Correlations between tropopause height and total ozone:Implications for long-term changes[J].J Gepohys Res,103(D15):19183-19192.

    Steinbrecht W,Haβler B,Claude H,et al.2003.Global distribution of total ozone and lower stratospheric temperature variations[J].Atmos Chem Phys,3(4):3411-3449.

    Yanai M,Li C F,Song Z.1992.Seasonal heating of the Tibetan plateau and its effects on the evolution of the Asian summer monsoon[J].J Meteor Soc Japan,70(1):319-351.

    Ye D Z,Wu G X.1998.The role of the heat source of the Tibetan Plateau in the general circulation[J].Meteor Atmos Phys,67(1): 181-198.

    Zhang Renhe,Zhou Shunwu.2009.Air Temperature changes over the Tibetan Plateau and other regions in the same latitudes and the role of ozone depletion[J].Acta Meteor Sinica,23(3):290-299.

    Zhou Shunwu,Zhang Renhe.2005.Decadal variations of temperature and geopotential height over the Tibetan Plateau and their relations with Tibet ozone depletion[J].Geophys Res Lett,32(18),18705.doi:10.1029/2005GL023496.

    Ziemke J R,Chandra S.2005.A 25-year data record of atmospheric ozone in the Pacific from Total Ozone Mapping Spectrometer (TOMS)cloud slicing:Implications for ozone trends in the stratosphere and troposphere[J].J Geophys Res,110(D15),105.doi: 10.1029/2004JD005687.

    Zou Han,Gao Yongqi.1997.Vertical ozone profile over Tibet using Sage I and II data[J].Adv Atmos Sci,14(4):505-512.

    (責(zé)任編輯:倪東鴻)

    P403

    A

    1674-7097(2012)04-0438-10

    楊雙艷,周順武,張人禾,等.2012.青藏高原對(duì)流層頂高度與臭氧總量及上升運(yùn)動(dòng)的耦合關(guān)系[J].大氣科學(xué)學(xué)報(bào),35(4):438-447.

    Yang Shuang-yan,Zhou Shun-wu,Zhang Ren-he,et al.2012.Coupling relationship between tropopause height and total ozone as well as ascending motion over the Tibetan Plateau[J].Trans Atmos Sci,35(4):438-447.(in Chinese)

    2011-08-06;改回日期:2012-05-22

    公益性行業(yè)(氣象)科研專項(xiàng)(GYHY(QX)200906014);國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃項(xiàng)目(2010CB428505);中國(guó)氣象局成都高原氣象研究所開(kāi)放實(shí)驗(yàn)室基金項(xiàng)目(LPM2011015)

    楊雙艷(1983—),女,湖北隨州人,博士生,研究方向?yàn)閰^(qū)域氣候變化,152850201@163.com;周順武(通信作者),男,博士,教授,研究方向?yàn)閰^(qū)域氣候變化,zhou@nuist.edu.cn.

    猜你喜歡
    季節(jié)變化青藏高原
    破解青藏高原氣候變化的密碼
    軍事文摘(2022年20期)2023-01-10 07:18:22
    青藏高原上的“含羞花”
    給青藏高原的班公湖量體溫
    軍事文摘(2021年20期)2021-11-10 01:58:40
    為了讓青藏高原的天更藍(lán)、水更綠、草原更美
    長(zhǎng)江師范學(xué)院校園鳥(niǎo)類的季節(jié)性變動(dòng)規(guī)律初探
    不同土地利用方式土壤酶活性的季節(jié)變化
    綠色科技(2017年10期)2017-07-05 13:06:54
    金華市小型封閉水體氮污染季節(jié)變化規(guī)律及其原因分析
    綠色科技(2016年24期)2017-03-30 17:24:11
    基于GIS的東川地區(qū)植被指數(shù)季節(jié)變化研究
    商(2016年22期)2016-07-08 14:05:14
    季節(jié)變化對(duì)蛋雞養(yǎng)殖的危害及其避抗措施
    不同類型人工濕地凈化水質(zhì)季節(jié)變化分析
    欧美潮喷喷水| 久久人人精品亚洲av| 一卡2卡三卡四卡精品乱码亚洲| 日本黄色片子视频| 国产成人影院久久av| 五月玫瑰六月丁香| 国产精品久久久久久精品电影| 自拍偷自拍亚洲精品老妇| 日本黄大片高清| 亚洲欧美日韩东京热| av天堂中文字幕网| 又黄又爽又刺激的免费视频.| 九九爱精品视频在线观看| 亚洲av免费高清在线观看| 毛片女人毛片| 成人特级av手机在线观看| 国产高清有码在线观看视频| 国产成人91sexporn| 久久久成人免费电影| 内射极品少妇av片p| 伦理电影大哥的女人| 熟妇人妻久久中文字幕3abv| 看免费成人av毛片| 女的被弄到高潮叫床怎么办| 免费观看的影片在线观看| 亚洲最大成人手机在线| 国产黄色小视频在线观看| 亚洲高清免费不卡视频| 搡老妇女老女人老熟妇| 亚洲av一区综合| av在线天堂中文字幕| 国产真实伦视频高清在线观看| 久久99热这里只有精品18| 国产高清激情床上av| 精华霜和精华液先用哪个| 国产毛片a区久久久久| 国产成人精品久久久久久| 欧美日韩一区二区视频在线观看视频在线 | 啦啦啦韩国在线观看视频| 国产av一区在线观看免费| 精品久久久久久久久久免费视频| 国产精品久久久久久久久免| 男女那种视频在线观看| 高清在线视频一区二区三区 | а√天堂www在线а√下载| 亚洲美女视频黄频| 中文字幕av成人在线电影| 久久热精品热| 一卡2卡三卡四卡精品乱码亚洲| 国产单亲对白刺激| 久久6这里有精品| 能在线免费看毛片的网站| 亚洲四区av| 51国产日韩欧美| 男插女下体视频免费在线播放| av女优亚洲男人天堂| 国产亚洲精品久久久com| 狂野欧美白嫩少妇大欣赏| 久久综合国产亚洲精品| 一个人观看的视频www高清免费观看| 亚洲精品色激情综合| 中文字幕熟女人妻在线| av在线播放精品| 午夜福利在线观看吧| 久久精品夜色国产| 好男人视频免费观看在线| 在线观看午夜福利视频| 欧美zozozo另类| 亚洲欧美日韩高清专用| 超碰av人人做人人爽久久| 欧美性感艳星| 观看美女的网站| 国产一区亚洲一区在线观看| 美女xxoo啪啪120秒动态图| 美女内射精品一级片tv| 日韩欧美精品免费久久| 高清毛片免费观看视频网站| 简卡轻食公司| 亚洲精品456在线播放app| 国产精华一区二区三区| 麻豆成人午夜福利视频| 国产精品一二三区在线看| 亚洲美女搞黄在线观看| 精品久久久久久成人av| 中文字幕人妻熟人妻熟丝袜美| 午夜激情欧美在线| 99热6这里只有精品| 日本撒尿小便嘘嘘汇集6| 亚洲在线自拍视频| 亚洲精华国产精华液的使用体验 | 国产不卡一卡二| 色哟哟·www| 色哟哟·www| 噜噜噜噜噜久久久久久91| 观看美女的网站| 好男人视频免费观看在线| а√天堂www在线а√下载| 一本一本综合久久| 中文在线观看免费www的网站| 不卡视频在线观看欧美| 直男gayav资源| 欧美一级a爱片免费观看看| 国内精品一区二区在线观看| 亚洲人成网站在线观看播放| 国语自产精品视频在线第100页| 亚洲欧美日韩东京热| 在线观看av片永久免费下载| 亚洲精品成人久久久久久| 麻豆国产97在线/欧美| 成人午夜高清在线视频| 欧美色视频一区免费| 中文字幕免费在线视频6| 99久久精品一区二区三区| 亚洲中文字幕日韩| 国产亚洲5aaaaa淫片| 久久久久久伊人网av| 国产精品精品国产色婷婷| 亚洲va在线va天堂va国产| 久久久国产成人精品二区| 精品久久久噜噜| 黄色配什么色好看| 日韩人妻高清精品专区| www.色视频.com| 免费观看人在逋| 麻豆成人午夜福利视频| 老司机福利观看| 日韩一区二区视频免费看| 国产熟女欧美一区二区| 国产精品福利在线免费观看| 床上黄色一级片| 好男人视频免费观看在线| 99热这里只有是精品在线观看| 一本精品99久久精品77| a级一级毛片免费在线观看| 插逼视频在线观看| 男女下面进入的视频免费午夜| 免费观看人在逋| 非洲黑人性xxxx精品又粗又长| 精品久久久噜噜| 成人亚洲欧美一区二区av| 精品人妻视频免费看| 亚洲无线在线观看| 国产成人影院久久av| 国产精品久久电影中文字幕| 久久精品国产亚洲av涩爱 | 国产淫片久久久久久久久| 一本精品99久久精品77| 久久久久久久久久久免费av| 亚洲欧美日韩无卡精品| 狂野欧美激情性xxxx在线观看| 久久午夜亚洲精品久久| 国产乱人偷精品视频| 久久久久国产网址| 久久99蜜桃精品久久| 国产av不卡久久| 亚洲欧美清纯卡通| 青春草视频在线免费观看| 中国美白少妇内射xxxbb| 又粗又爽又猛毛片免费看| 国产精品人妻久久久影院| 男女视频在线观看网站免费| 精品日产1卡2卡| 亚洲在久久综合| 国产三级中文精品| 中文精品一卡2卡3卡4更新| 一区二区三区免费毛片| 亚洲图色成人| 岛国毛片在线播放| 日本色播在线视频| 噜噜噜噜噜久久久久久91| 最近的中文字幕免费完整| 久久久精品大字幕| 人体艺术视频欧美日本| 久久精品久久久久久噜噜老黄 | 一边摸一边抽搐一进一小说| 如何舔出高潮| 日日摸夜夜添夜夜添av毛片| 欧美精品一区二区大全| 99在线人妻在线中文字幕| 亚洲精品国产av成人精品| 99久久人妻综合| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美成人综合另类久久久 | 国产人妻一区二区三区在| 亚洲人成网站高清观看| 国产亚洲精品久久久久久毛片| 美女国产视频在线观看| 99久久无色码亚洲精品果冻| 亚洲精品456在线播放app| 内射极品少妇av片p| 国产精品无大码| 美女xxoo啪啪120秒动态图| 99热这里只有是精品在线观看| 久久亚洲国产成人精品v| 亚洲欧美精品专区久久| 青青草视频在线视频观看| 久久99热6这里只有精品| 中国美白少妇内射xxxbb| 亚洲乱码一区二区免费版| 最近视频中文字幕2019在线8| 免费在线观看成人毛片| 人人妻人人看人人澡| 蜜臀久久99精品久久宅男| 亚洲最大成人av| 国产中年淑女户外野战色| 免费一级毛片在线播放高清视频| 1000部很黄的大片| 黑人高潮一二区| 亚洲四区av| 国产欧美日韩精品一区二区| 99久久久亚洲精品蜜臀av| АⅤ资源中文在线天堂| 人妻少妇偷人精品九色| 亚洲不卡免费看| 99久久九九国产精品国产免费| 精品国产三级普通话版| av卡一久久| 亚洲美女搞黄在线观看| 国产黄片美女视频| 亚洲美女视频黄频| 国产精品久久视频播放| 日本在线视频免费播放| 国产精品1区2区在线观看.| 色噜噜av男人的天堂激情| 两个人视频免费观看高清| 如何舔出高潮| 不卡视频在线观看欧美| 热99在线观看视频| 天天躁日日操中文字幕| 别揉我奶头 嗯啊视频| 变态另类成人亚洲欧美熟女| 亚洲综合色惰| av卡一久久| 精品人妻一区二区三区麻豆| 欧美区成人在线视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人无遮挡网站| 亚洲国产精品国产精品| 夫妻性生交免费视频一级片| 精品一区二区免费观看| 插阴视频在线观看视频| 国产高清三级在线| 12—13女人毛片做爰片一| 国产一区二区在线av高清观看| 非洲黑人性xxxx精品又粗又长| 久久亚洲国产成人精品v| 精华霜和精华液先用哪个| 亚洲婷婷狠狠爱综合网| 国产综合懂色| 国产爱豆传媒在线观看| 美女xxoo啪啪120秒动态图| 99久国产av精品| 国产毛片a区久久久久| 99热网站在线观看| 日韩亚洲欧美综合| 天堂网av新在线| 又粗又爽又猛毛片免费看| 国产高清有码在线观看视频| 精品久久久久久久久久久久久| 不卡一级毛片| 噜噜噜噜噜久久久久久91| 成人鲁丝片一二三区免费| 婷婷亚洲欧美| www.色视频.com| av女优亚洲男人天堂| 伦精品一区二区三区| 午夜a级毛片| 日韩欧美精品免费久久| 联通29元200g的流量卡| 日本与韩国留学比较| av在线播放精品| 国产高清三级在线| 在线观看av片永久免费下载| 我要搜黄色片| 亚洲精品乱码久久久v下载方式| 精品不卡国产一区二区三区| 日韩欧美三级三区| 秋霞在线观看毛片| 不卡视频在线观看欧美| 久久99热这里只有精品18| 欧美色视频一区免费| 日韩成人伦理影院| 久久99精品国语久久久| 久久久久久久午夜电影| 性欧美人与动物交配| 网址你懂的国产日韩在线| 亚洲av.av天堂| 午夜爱爱视频在线播放| 久久久久久久亚洲中文字幕| 欧美区成人在线视频| 亚洲av免费在线观看| 一本久久精品| 色尼玛亚洲综合影院| 精品不卡国产一区二区三区| 爱豆传媒免费全集在线观看| 国产真实乱freesex| 国产午夜精品论理片| 婷婷六月久久综合丁香| 久99久视频精品免费| 少妇熟女aⅴ在线视频| 久久99热这里只有精品18| 99riav亚洲国产免费| 美女被艹到高潮喷水动态| av天堂在线播放| 99国产精品一区二区蜜桃av| 天堂影院成人在线观看| 亚洲一区二区三区色噜噜| 美女国产视频在线观看| 国产淫片久久久久久久久| 色5月婷婷丁香| 日韩人妻高清精品专区| 一区二区三区四区激情视频 | 日韩一本色道免费dvd| 超碰av人人做人人爽久久| 色尼玛亚洲综合影院| 12—13女人毛片做爰片一| 一个人观看的视频www高清免费观看| 久久久久久久亚洲中文字幕| 日韩一本色道免费dvd| 一级二级三级毛片免费看| 成人二区视频| 免费av不卡在线播放| 97在线视频观看| 色哟哟哟哟哟哟| 国产伦理片在线播放av一区 | a级毛片a级免费在线| 日本撒尿小便嘘嘘汇集6| 欧美日韩综合久久久久久| 国产成年人精品一区二区| 秋霞在线观看毛片| 国产精品99久久久久久久久| 大型黄色视频在线免费观看| 日韩三级伦理在线观看| 久久九九热精品免费| av在线亚洲专区| 久久人人精品亚洲av| 一边摸一边抽搐一进一小说| 啦啦啦观看免费观看视频高清| 亚洲精品日韩av片在线观看| 69av精品久久久久久| 一级av片app| 日韩强制内射视频| 伦精品一区二区三区| 国产精品久久久久久av不卡| 五月玫瑰六月丁香| 亚洲精品日韩av片在线观看| 亚洲欧美精品专区久久| 成人三级黄色视频| 精品一区二区三区人妻视频| 高清午夜精品一区二区三区 | 床上黄色一级片| 天美传媒精品一区二区| 欧美zozozo另类| 久久久久九九精品影院| 国产69精品久久久久777片| 国产视频内射| 国产 一区 欧美 日韩| 亚洲,欧美,日韩| 国产高清不卡午夜福利| 欧美+日韩+精品| 精品熟女少妇av免费看| 日韩欧美精品v在线| 色吧在线观看| 黄色视频,在线免费观看| 好男人视频免费观看在线| 一进一出抽搐动态| 欧美另类亚洲清纯唯美| 中文在线观看免费www的网站| 丝袜喷水一区| 久久久久免费精品人妻一区二区| 大香蕉久久网| 久久99精品国语久久久| 欧美日韩精品成人综合77777| 久久久久久久久久黄片| 国产黄色小视频在线观看| 亚洲av免费高清在线观看| 色视频www国产| 国产成人精品一,二区 | 国产精品福利在线免费观看| 久久精品国产清高在天天线| 国产高清视频在线观看网站| 亚洲欧美成人综合另类久久久 | 成人特级黄色片久久久久久久| av天堂在线播放| 哪里可以看免费的av片| 亚洲无线在线观看| 插阴视频在线观看视频| 中文字幕免费在线视频6| 青春草视频在线免费观看| 两个人的视频大全免费| 26uuu在线亚洲综合色| 长腿黑丝高跟| 免费大片18禁| 狂野欧美白嫩少妇大欣赏| 99久久成人亚洲精品观看| 小蜜桃在线观看免费完整版高清| 成年版毛片免费区| 人体艺术视频欧美日本| 欧美日韩在线观看h| 日本免费a在线| 欧美3d第一页| 色5月婷婷丁香| 三级国产精品欧美在线观看| 精品人妻视频免费看| 少妇的逼好多水| 日韩高清综合在线| 亚洲一区高清亚洲精品| 亚洲av成人精品一区久久| 成人午夜高清在线视频| 嘟嘟电影网在线观看| 啦啦啦韩国在线观看视频| 免费av不卡在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 国产白丝娇喘喷水9色精品| 看十八女毛片水多多多| 伦精品一区二区三区| 国产高清不卡午夜福利| 波多野结衣高清无吗| 成人特级av手机在线观看| 性插视频无遮挡在线免费观看| 成人鲁丝片一二三区免费| 日韩欧美在线乱码| 国产爱豆传媒在线观看| 亚洲熟妇中文字幕五十中出| 日韩欧美精品免费久久| 久久久久久久久久久丰满| 大香蕉久久网| 久久亚洲精品不卡| 亚洲国产精品成人综合色| 久久精品久久久久久久性| 亚洲,欧美,日韩| 久久国产乱子免费精品| 少妇裸体淫交视频免费看高清| 欧美色欧美亚洲另类二区| 欧美高清成人免费视频www| 欧美日本亚洲视频在线播放| 一级毛片电影观看 | 日本-黄色视频高清免费观看| 菩萨蛮人人尽说江南好唐韦庄 | 夜夜看夜夜爽夜夜摸| 国产国拍精品亚洲av在线观看| 1024手机看黄色片| 青春草视频在线免费观看| 国产精品综合久久久久久久免费| 九色成人免费人妻av| 一本一本综合久久| av.在线天堂| 日韩欧美三级三区| 2022亚洲国产成人精品| 伦理电影大哥的女人| 天堂av国产一区二区熟女人妻| 蜜桃久久精品国产亚洲av| 亚洲国产精品成人久久小说 | 久久99热6这里只有精品| 成人午夜精彩视频在线观看| 日本黄大片高清| av免费在线看不卡| 18禁裸乳无遮挡免费网站照片| 久久这里有精品视频免费| 午夜福利在线在线| 搡老妇女老女人老熟妇| 99热全是精品| 免费av观看视频| 蜜臀久久99精品久久宅男| 精品久久久久久久久av| 国产精品久久久久久精品电影小说 | 最好的美女福利视频网| 国产av不卡久久| 一个人观看的视频www高清免费观看| 99久国产av精品| 午夜视频国产福利| 男的添女的下面高潮视频| 亚洲av成人av| 日韩欧美精品免费久久| 午夜精品一区二区三区免费看| 人妻夜夜爽99麻豆av| 18禁在线播放成人免费| 欧美丝袜亚洲另类| 熟妇人妻久久中文字幕3abv| 久久精品国产清高在天天线| 午夜激情欧美在线| 伦精品一区二区三区| 少妇被粗大猛烈的视频| 丝袜美腿在线中文| 晚上一个人看的免费电影| www.av在线官网国产| 26uuu在线亚洲综合色| 国产高清视频在线观看网站| 性插视频无遮挡在线免费观看| 国产91av在线免费观看| 婷婷精品国产亚洲av| 日本黄大片高清| 久久久久久久久大av| 搡老妇女老女人老熟妇| 一个人观看的视频www高清免费观看| 国产精品永久免费网站| 国产午夜福利久久久久久| 亚洲18禁久久av| 国产淫片久久久久久久久| 欧美高清成人免费视频www| 久久99蜜桃精品久久| 成人漫画全彩无遮挡| 一级黄片播放器| 国产伦一二天堂av在线观看| 色综合亚洲欧美另类图片| 女的被弄到高潮叫床怎么办| 波多野结衣高清无吗| 国产黄色视频一区二区在线观看 | 国产精品久久久久久久电影| 精品免费久久久久久久清纯| 中文字幕av在线有码专区| 精品久久久久久久久av| 亚洲欧美日韩高清在线视频| 国产精品一及| 高清毛片免费看| 26uuu在线亚洲综合色| 亚洲精品456在线播放app| 亚洲精品日韩在线中文字幕 | 一进一出抽搐gif免费好疼| 给我免费播放毛片高清在线观看| 国产精品国产高清国产av| 两个人视频免费观看高清| 久久鲁丝午夜福利片| 亚洲av一区综合| 欧美高清成人免费视频www| 成人特级av手机在线观看| 麻豆一二三区av精品| 精品午夜福利在线看| 国产亚洲av片在线观看秒播厂 | 青春草国产在线视频 | 青春草亚洲视频在线观看| 国产精品野战在线观看| 国产黄a三级三级三级人| 国产成人精品一,二区 | 丝袜喷水一区| 久久久成人免费电影| 日本三级黄在线观看| 亚洲精品日韩av片在线观看| 免费看日本二区| 欧美bdsm另类| 欧美另类亚洲清纯唯美| 国产高清激情床上av| 国产久久久一区二区三区| 国产蜜桃级精品一区二区三区| 亚洲图色成人| 久久亚洲精品不卡| 亚洲国产精品久久男人天堂| 少妇被粗大猛烈的视频| 校园春色视频在线观看| 精品久久久久久久久亚洲| 麻豆久久精品国产亚洲av| 男人的好看免费观看在线视频| 亚洲电影在线观看av| 久久九九热精品免费| 人妻夜夜爽99麻豆av| 国产精品久久久久久精品电影小说 | 18禁在线播放成人免费| 久久久久免费精品人妻一区二区| 深爱激情五月婷婷| 久久亚洲国产成人精品v| 深爱激情五月婷婷| 国产免费男女视频| 一个人免费在线观看电影| 欧美人与善性xxx| 日本免费a在线| 日韩,欧美,国产一区二区三区 | 亚洲精品影视一区二区三区av| 岛国毛片在线播放| 中国美女看黄片| 全区人妻精品视频| 久久精品国产亚洲av香蕉五月| 干丝袜人妻中文字幕| 国产精品久久久久久亚洲av鲁大| 最近2019中文字幕mv第一页| 亚洲中文字幕一区二区三区有码在线看| 成人性生交大片免费视频hd| 久久久久久九九精品二区国产| 亚洲国产色片| ponron亚洲| 久久精品久久久久久久性| 免费观看a级毛片全部| 亚洲国产精品久久男人天堂| 久久久久久久久大av| 亚洲在线自拍视频| 有码 亚洲区| 尾随美女入室| 国产高清有码在线观看视频| 成人高潮视频无遮挡免费网站| 天天躁日日操中文字幕| 白带黄色成豆腐渣| 18禁在线无遮挡免费观看视频| 你懂的网址亚洲精品在线观看 | 非洲黑人性xxxx精品又粗又长| 人人妻人人澡欧美一区二区| 麻豆av噜噜一区二区三区| 美女大奶头视频| 晚上一个人看的免费电影| 国产高潮美女av| 夜夜爽天天搞| 国产精品人妻久久久久久| 中文精品一卡2卡3卡4更新| 国产精品久久电影中文字幕| 午夜福利在线在线| 久久久久网色| 亚洲av二区三区四区| 久久精品综合一区二区三区| 国内久久婷婷六月综合欲色啪| 波多野结衣高清无吗| 久久久色成人| av福利片在线观看| 欧美日韩乱码在线| 99精品在免费线老司机午夜| 91麻豆精品激情在线观看国产| 干丝袜人妻中文字幕| 久久欧美精品欧美久久欧美|