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

    基于多臺(tái)站OH全天空氣輝成像儀觀測(cè)的中國(guó)中緯地區(qū)重力波傳播特性

    2016-06-30 07:27:20王翠梅李欽增徐寄遙袁韡
    地球物理學(xué)報(bào) 2016年5期
    關(guān)鍵詞:波源傳播

    王翠梅, 李欽增, 徐寄遙, 袁韡

    1 中國(guó)科學(xué)院空間科學(xué)與應(yīng)用研究中心天氣學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100190 2 中國(guó)科學(xué)院大學(xué), 北京 100049

    基于多臺(tái)站OH全天空氣輝成像儀觀測(cè)的中國(guó)中緯地區(qū)重力波傳播特性

    王翠梅1,2, 李欽增1, 徐寄遙1, 袁韡1

    1 中國(guó)科學(xué)院空間科學(xué)與應(yīng)用研究中心天氣學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京100190 2 中國(guó)科學(xué)院大學(xué), 北京100049

    摘要利用位于中國(guó)中緯地區(qū)6個(gè)OH氣輝成像儀2012年1月至2013年12月兩年的觀測(cè)數(shù)據(jù), 我們研究分析了重力波傳播特征.結(jié)果表明重力波的水平波長(zhǎng)、觀測(cè)周期和水平相速度分別主要分布于10~35 km, 4~12 min和30~100 m·s-1范圍. 夏季,重力波主要沿極向方向傳播. 然而, 冬季, 他們有向赤道方向和平行于赤道方向的傳播趨勢(shì). 同時(shí), 我們結(jié)合TRMM衛(wèi)星和ECMWF數(shù)據(jù), 發(fā)現(xiàn)在夏季重力波的北向傳播趨勢(shì)可能主要由觀測(cè)臺(tái)站南方的對(duì)流活動(dòng)導(dǎo)致. 然而, 對(duì)流層頂附近的急流可能在冬季重力波的主要傳播方向方面做出較大貢獻(xiàn). 分析結(jié)果也表明, 低層-中層大氣背景風(fēng)的濾波效應(yīng)僅在夏季與中國(guó)中緯地區(qū)重力波緯向傳播方向各向異性吻合較好.

    關(guān)鍵詞氣輝成像儀; 大氣重力波; 傳播; 波源; 濾波效應(yīng)

    1引言

    大氣重力波在中層-低熱層動(dòng)力學(xué)過(guò)程中起到重要作用(Nakamura et al., 2003).普遍的觀點(diǎn)認(rèn)為,中層-低熱層大氣重力波產(chǎn)生于低層大氣且向上傳播(Taylor and Hapgood, 1988).因此,在重力波向上傳播過(guò)程中,低層-中層大氣風(fēng)場(chǎng)對(duì)重力波的通量和方位角分布起到重要作用(Taylor et al., 1993).

    近幾十年,研究者利用氣輝成像觀測(cè)數(shù)據(jù)對(duì)不同緯度地區(qū)的中層-低熱層高頻重力波傳播特征進(jìn)行了研究.氣輝成像儀具有獲取二維水平結(jié)構(gòu)的能力,被廣泛用于短周期(<1 h) 和小尺度重力波觀測(cè)(水平波長(zhǎng)<100 km).

    中緯地區(qū),研究者(Taylor et al., 1993; Nakamura et al., 1999; Walterscheid et al., 1999; Hecht et al., 2001, 2004; Ejiri et al., 2003; Suzuki et al., 2004; Tang et al., 2005, 2014; Dou et al., 2010; Kim et al., 2010; Li et al, 2011b)利用氣輝成像觀測(cè)數(shù)據(jù),發(fā)現(xiàn)重力波在夏季主要沿東向和極向傳播.然而,很多觀測(cè)報(bào)道了中緯地區(qū)的重力波在冬季主要沿赤道方向傳播(Walterscheid et al., 1999; Ejiri et al., 2003; Suzuki et al., 2004; Tang et al., 2005, 2014; Dou et al., 2010; Kim et al., 2010; Li et al., 2011b).

    低緯地區(qū),Li等(2011a)和Suzuki等(2004)利用OH氣輝成像數(shù)據(jù)分別報(bào)道了Maui(20.7°N,156.3°W)和Australia(12°S)臺(tái)站的重力波在夏季(冬季)趨向于極向(赤道方向)傳播.

    赤道附近,Narayanan和Gurubaran(2013)報(bào)道了Tirunelveli(8.7° N)地區(qū)的重力波在夏季幾乎都向南和西南方向傳播.Medeiros等(2004a)的研究結(jié)果表明,在夏季和冬季,Cariri(7°S,35°W)地區(qū)的重力波分別主要沿東南和東北方向傳播.

    影響重力波傳播方向的因素有很多.其中,波源位置被廣泛認(rèn)為是影響重力波傳播方向的主要因素之一,且對(duì)流活動(dòng)是產(chǎn)生高頻重力波的重要波源之一.模擬(Fovell et al., 1992; Alexander et al., 1995, 2004; Lane et al., 2001; Vadas and Fritts, 2001, 2009; Walterscheid et al., 2001; Song et al., 2003; Vadas et al., 2009)和觀測(cè)(Walterscheid et al., 1999; Hecht et al., 2001; Ejiri et al., 2003; Suzuki et al., 2004; Pautet et al., 2005; Yue et al., 2009)研究給出了對(duì)流活動(dòng)產(chǎn)生重力波的有力證據(jù).另外,有研究表明導(dǎo)管傳播和波源位置共同影響重力波主要傳播方向.Walterscheid等(1999)報(bào)道了Adelaide(35°S)地區(qū)上空的重力波在夏季(冬季)主要沿極向(赤道方向)傳播.他們解釋了產(chǎn)生這種現(xiàn)象是由于重力波陷入導(dǎo)管傳播.Hecht等(2001)通過(guò)分析 Urbana,Illinois(40°N,80°W)臺(tái)站的的氣輝成像觀測(cè)數(shù)據(jù),報(bào)道了重力波在夏季的北向傳播趨勢(shì)是因?yàn)闊釋?dǎo)管和風(fēng)剪切的存在。

    另一個(gè)產(chǎn)生重力波的重要機(jī)制是急流.Shiokawa等(2007)利用FPI干涉儀和557.7 nm氣輝成像儀在Shigaraki(34.8°N,136.1°E)臺(tái)站觀測(cè)到一個(gè)強(qiáng)烈的重力波事件.他們的分析結(jié)果表明,觀測(cè)的重力波產(chǎn)生于日本的東北方向,該處對(duì)流層頂附近存在一個(gè)低壓系統(tǒng)(968 hPa)和急流.

    另外,研究者(Taylor et al., 1993; Stockwell and Lowe, 2001a,2001b; Medeiros et al., 2003,2004a, 2004b, 2005)發(fā)現(xiàn)臨界層濾波效應(yīng)在決定重力波緯向傳播方向各向異性方面起到重要作用.Taylor等(1993)假設(shè)波源位于對(duì)流層且使用CIRA-1986緯向風(fēng)數(shù)據(jù)報(bào)道了Colorado(40.0°N 105.6°W)臺(tái)站的重力波從低層向高層大氣傳播過(guò)程中有可能被臨界層濾波效應(yīng)濾除.Stockwell和Lowe(2001a,2001b)認(rèn)為臨界層的緯向風(fēng)的濾波效應(yīng)解釋了重力波在夏季和冬季分別主要沿東向和西向傳播的現(xiàn)象.Medeiros等(2003, 2004a, 2004b, 2005)使用Taylor等(1993)中的濾波理論和HWM93水平風(fēng)場(chǎng)模型(Hedin et al., 1996)報(bào)道了濾波效應(yīng)決定重力波傳播方向各向異性方面起主要作用.然而,Li等(2011b)利用Xinglong(40.2°N,117.4°E)臺(tái)站的觀測(cè)數(shù)據(jù)對(duì)重力波傳播特性進(jìn)行研究,結(jié)果表明臨界層濾波效應(yīng)僅在夏季與重力波傳播方向吻合較好.他們指出,波源的非均勻分布和導(dǎo)管傳播也同時(shí)影響重力波傳播方向各向異性.

    Medeiros等(2005)通過(guò)分析Atlantic附近的多臺(tái)站觀測(cè)數(shù)據(jù),發(fā)現(xiàn)重力波主要從陸地向海洋傳播.他們認(rèn)為來(lái)自陸地的對(duì)流活動(dòng)在決定沿海地區(qū)重力波傳播方向起重要作用.

    Kim等(2010)利用射線追蹤技術(shù)研究大氣重力波,發(fā)現(xiàn)大部分波動(dòng)來(lái)自于中層高度,他們認(rèn)為這些重力波為二次波.

    本文利用2012年1月至2013年12月兩年位于中國(guó)中緯地區(qū)6個(gè)臺(tái)站的OH氣輝成像儀觀測(cè)數(shù)據(jù),對(duì)中國(guó)中緯地區(qū)重力波傳播特征進(jìn)行了分析.同時(shí),我們結(jié)合TRMM(Tropical Rainfall Measuring Mission)衛(wèi)星觀測(cè)數(shù)據(jù),ECMWF(the European Center for Medium-Range Weather Forecasts)和 MERRA(the Modern Era Retrospective-analysis for Research and Applications)風(fēng)場(chǎng)數(shù)據(jù)重點(diǎn)討論了波源位置和濾波效應(yīng)在重力波傳播方向季節(jié)變化方面的影響.文章結(jié)構(gòu)安排如下: 第2部分給出觀測(cè)數(shù)據(jù)和處理方法.結(jié)果和討論將在第3部分呈現(xiàn).第4部分是全文總結(jié).

    2數(shù)據(jù)和方法

    用于研究重力波傳播特征的數(shù)據(jù)來(lái)自位于中國(guó)中緯地區(qū)6個(gè)臺(tái)站2012年1月至2013年12月的OH氣輝成像儀,它們分別位于: Shuozhou(39.8°N,112.1°E), Xinglong(40.4°N,117.6°E),Donggang(40.0°N,124.0°E),Xinxiang(35.7°N,113.7°E),Linqu(36.2°N,118.7°E)和 Rongcheng(37.3°N,122.5°E).其中,Xinglong(40.4°N,117.6°E)臺(tái)站的OH氣輝成像儀是子午工程(Wang,2010)地基光學(xué)觀測(cè)儀中的一種.每個(gè)OH氣輝成像儀包含溫度穩(wěn)定的濾波輪子,使用Nikon 16mm/2.8D 魚(yú)眼鏡頭,視場(chǎng)角為180°,且由像元數(shù)為1024 pixel×1024 pixel的CCD感光元件組成.整個(gè)觀測(cè)過(guò)程由軟件來(lái)自動(dòng)控制.曝光時(shí)間是1 min.且OH氣輝輻射波段位于715~930 nm,觀測(cè)高度為87±5 km.圖1為OH氣輝成像儀的地理位置示意圖.

    這篇文章中使用的重力波參數(shù)(水平波長(zhǎng)、觀測(cè)周期和水平相速度)是按照如下的步驟提取: 第一步,利用中值濾波法去除原始圖像中星光的污染(Suzuki et al.,2007).第二步,去除Van Rhijin效應(yīng)和大氣消光效應(yīng)(Kubota et al.,2001).第三步,假設(shè)OH氣輝層的高度為87 km,且將氣輝圖像投影到800 km×800 km 的地理坐標(biāo)平面內(nèi)(Baker and Stair,1988).第四步,把連續(xù)的兩幅圖片進(jìn)行相減得到差分(TD)圖像(Tang et al.,2005).最后,我們對(duì)差分圖像進(jìn)行去勢(shì)處理且利用2-D FFT光譜分析方法提取重力波參數(shù)(Garcia et al., 1997; Coble et al, 1998).

    圖1 OH氣輝成像儀地理位置示意圖Fig.1 A map showing locations of OH all-sky airglow imagers

    3結(jié)果和討論

    為了研究重力波的季節(jié)變化特征,我們將全年劃分為四個(gè)季節(jié): 春季(3月—4月)、夏季(5月—8月)、秋季(9月—10月)和冬季(11月—次年2月).利用OH氣輝成像儀對(duì)重力波傳播特征研究的文章(Dou et al.,2010; Ejiri et al., 2003; Li et al., 2011a, 2011b; Medeiros et al., 2004a, 2004b, 2005; Nakamura et al., 2003; Narayanan and Gurubaran, 2013; 王翠梅等, 2014)幾乎呈現(xiàn)出重力波的觀測(cè)周期<1 h,因此在這篇文章中,我們用于研究重力波傳播特征的事件不包含水平相速度趨近于0(觀測(cè)周期遠(yuǎn)大于1 h)和圓心能被看見(jiàn)的環(huán)形重力波(因?yàn)閳A心被看見(jiàn)的環(huán)形重力波,其傳播分布于各個(gè)方向).3.1重力波的傳播特征

    表1列出了各個(gè)觀測(cè)臺(tái)站的重力波觀測(cè)條件以及觀測(cè)的重力波事件數(shù)目.這里,晴空夜晚定義為半小時(shí)及以上的觀測(cè)圖像沒(méi)有云的污染的夜晚.表1表明,夏季的晴空夜晚數(shù)最少.且晴空夜晚數(shù)在Shuozhou,Xinglong 和Donggang(這3個(gè)臺(tái)站在40°N)明顯多于其他3個(gè)臺(tái)站(36°N附近).

    圖2(a—f)分別是2012年1月至2013年12月期間各個(gè)觀測(cè)臺(tái)站的晴空(波動(dòng))小時(shí)數(shù)月分布圖.結(jié)果表明冬季的晴空時(shí)間明顯長(zhǎng)于夏季的,因?yàn)橹袊?guó)中緯地區(qū)夏季多雨.然而,幾乎每個(gè)月的晴空時(shí)間>20 h(除了Xinglong臺(tái)站6月晴空時(shí)間<20 h,因?yàn)樵撛掠^測(cè)天數(shù)較少),所以我們對(duì)重力波傳播特征研究的可信度是比較高的.

    圖3(a—f)反映了各個(gè)觀測(cè)臺(tái)站重力波發(fā)生頻率的月分布情況.這里,重力波發(fā)生頻率定義為產(chǎn)生波動(dòng)的時(shí)間與晴空時(shí)間的比例.結(jié)果表明,Donggang臺(tái)站的重力波發(fā)生頻率在各個(gè)月份幾乎相等.然而,在其他觀測(cè)臺(tái)站,重力波發(fā)生頻率在夏季和冬季高于春季和秋季.這可能由于波源的季節(jié)變化和傳播環(huán)境導(dǎo)致.且后者結(jié)果與幾個(gè)之前的研究結(jié)果相似(Nakamura et al., 1999; Medeiros et al., 2004a, 2004b; Kim et al., 2010; Tang et al., 2014).Nakamura等(1999)和Kim等(2010)分別發(fā)現(xiàn)Shigaraki(34.9°N,136.1°E)和Mt.Bohyun, Korea(36.2°N,128.9°E)地區(qū)的重力波發(fā)生頻率最大值在夏季和冬季.Medeiros等(2004a,2004b)報(bào)道了Cariri(7°S,35°W)和 Brazil(23°S,45°W)臺(tái)站的重力波發(fā)生頻率在夏季和冬季大于春季和秋季.Tang等(2014)的研究結(jié)果表明在CO(40.7°N,104.9°W)觀測(cè)臺(tái)站的重力波發(fā)生頻率表現(xiàn)出明顯的半年變化,且最大值在夏季和冬季,最小值在春季和秋季.

    表1 各個(gè)臺(tái)站的重力波觀測(cè)條件以及觀測(cè)的重力波事件數(shù)

    圖2 晴空時(shí)間(藍(lán)色)和波動(dòng)時(shí)間(紅色)月分布圖(a) Shuozhou; (b) Xinglong; (c) Donggang; (d) Xinxiang; (e) Linqu; (f) Rongcheng.Fig.2 Monthly distributions of clear hours (blue bars) and wave structure hours (red bars) at every observed station

    圖3 重力波發(fā)生頻率月分布圖(a) Shuozhou; (b) Xinglong; (c) Donggang; (d) Xinxiang; (e) Linqu; (f) Rongcheng.Fig.3 The monthly distribution of AGWs occurrence frequency

    圖4(a—f)分別給出每個(gè)觀測(cè)臺(tái)站的重力波參數(shù)分布(其中下標(biāo)1,2和3分別代表水平波長(zhǎng)、觀測(cè)周期和水平相速度).結(jié)果表明水平波長(zhǎng)、觀測(cè)周期和水平相速度分別主要分布于10~35 km,4~12 min和30~100 m·s-1范圍.

    圖5(a—f)分別是每個(gè)觀測(cè)臺(tái)站觀測(cè)的所有重力波事件的傳播方向極向分布.從圖中發(fā)現(xiàn)Rongcheng(37.3°N,122.5°E)和Donggang(40.0°N,124.0°E)臺(tái)站的重力波趨向于東向傳播,可能是因?yàn)檫@兩個(gè)臺(tái)站更接近于沿海且重力波主要來(lái)自大陸對(duì)流.然而,其他四個(gè)觀測(cè)臺(tái)站的重力波傳播方向主要表現(xiàn)為沿經(jīng)向方向,產(chǎn)生這一現(xiàn)象的原因?qū)⒃?.2節(jié)和3.3節(jié)做出討論.

    圖6—9表明重力波傳播方向具有明顯的季節(jié)變化.夏季,重力波在Donggang和Rongcheng趨向于東北方向傳播.在其他4個(gè)觀測(cè)臺(tái)站,幾乎所有的重力波沿北向傳播.冬季,重力波具有強(qiáng)烈的沿赤道(向南)方向和平行于赤道(向西或向東)方向傳播的趨勢(shì).另外,重力波在春季(除了Shuozhou臺(tái)站無(wú)明顯優(yōu)勢(shì)傳播方向)和秋季主要分別沿東北和西北方向傳播.與同緯度的觀測(cè)結(jié)果相比較,幾個(gè)研究者(Wu and Killeen, 1996; Walterscheid et al., 1999; Hecht et al., 2001; Ejiri et al., 2003; Nakamura et al., 2003; Stockwell and Lowe, 2001a; Li et al., 2011a)報(bào)道了重力波在夏季和冬季分別主要沿極向和赤道方向傳播.他們解釋重力波在夏季趨向于極向傳播主要是因?yàn)槌嗟栏浇鼘?duì)流源的影響.然而,中緯和高緯地區(qū)的地形因素和大氣條件可能導(dǎo)致冬季重力波趨向于赤道方向傳播.為了解釋重力波傳播方向在經(jīng)向方向的季節(jié)變化,衛(wèi)星觀測(cè)(Wu and Waters,1996; Nakamura et al.,2003)和環(huán)流模型(Richter et al., 2010)研究結(jié)果表明,夏季多對(duì)流活動(dòng),冬季多急流和鋒系.

    3.2重力波傳播方向影響因素分析

    3.2.1波源

    過(guò)去幾十年,數(shù)值計(jì)算(Horinouchi et al., 2002; Lane et al., 2003; Alexander et al., 2004; Vadas and Fritts, 2009)利用三維模型模擬了對(duì)流活動(dòng)產(chǎn)生重力波,加強(qiáng)了我們對(duì)重力波產(chǎn)生機(jī)制的理解.且隨后的觀測(cè)也驗(yàn)證了重力波產(chǎn)生于對(duì)流活動(dòng).Dewan和Picard(1998)利用天基觀測(cè)數(shù)據(jù),首次報(bào)道了對(duì)流雷暴可以產(chǎn)生重力波.Yue等(2009)利用Yucca Ridge Field Station(40.7°N,104.9°W)臺(tái)站的OH氣輝成像儀觀測(cè)數(shù)據(jù)報(bào)道了強(qiáng)對(duì)流活動(dòng)產(chǎn)生了環(huán)形重力波.

    降雨是反映深對(duì)流活動(dòng)的重要指標(biāo)(Walterscheid et al., 1999; Nakamura et al., 2003).圖10(a—d)是重力波事件發(fā)生前及期間的降雨量季節(jié)分布.結(jié)果表明,觀測(cè)臺(tái)站南方在夏季和春季秋季多降雨,可能在重力波北向傳播趨勢(shì)方面起主要作用.然而,由于觀測(cè)臺(tái)站附近在冬季的降雨較少,因此降雨分布

    圖4 重力波參數(shù)直方圖(a) Shuozhou; (b) Xinglong; (c) Donggang; (d) Xinxiang; (e) Linqu; (f) Rongcheng.(下標(biāo)1,2和3分別代表水平波長(zhǎng)、觀測(cè)周期和水平相速度)Fig.4 Histograms of gravity wave parametersThe subscript of 1,2 and 3 represents horizontal wavelength, observed period and horizontal phase speed, respectively.

    不能反映冬季的重力波傳播特征.這里,降雨數(shù)據(jù)從TRMM衛(wèi)星獲得,其空間分辨率為0.25°latitude×0.25°longitude,觀測(cè)覆蓋的范圍是180°W—180°E,50°S—50°N.

    Buss等(2004)利用ECMWF 數(shù)據(jù)給出Greenland上空對(duì)流層頂附近的急流產(chǎn)生重力波的有力證據(jù).Shiokawa等(2007)利用FPI干涉儀和557.7 nm氣輝成像儀在Shigaraki(34.8°N,136.1°E)臺(tái)站觀測(cè)到一個(gè)與急流有關(guān)的重力波事件.然而,我們通過(guò)ECMWF觀測(cè)數(shù)據(jù)獲得的風(fēng)場(chǎng)矢量圖發(fā)現(xiàn)觀測(cè)臺(tái)站附近上空約300 hPa(~9.2 km)處在冬季經(jīng)常出現(xiàn)兩支急流: 西風(fēng)急流(20°N—45°N)和西北風(fēng)急流(105°E—130°E)(這里沒(méi)有給出冬季平均風(fēng)場(chǎng)矢量圖).這兩支急流可能在影響冬季重力波的主要傳播方向方面起到重要作用.

    3.2.2背景風(fēng)的濾波效應(yīng)

    很多研究表明,低層-中層風(fēng)場(chǎng)的濾波效應(yīng)在決定重力波緯向傳播方向各向異性方面起重要作用(Taylor et al., 1993; Stockwell and Lowe, 2001b; Medeiros et al., 2003).夏季/冬季,重力波傳播方向的東向/西向傳播趨勢(shì)與臨界層西向風(fēng)/東向風(fēng)的濾波效應(yīng)一致(Nakamura et al., 1999).

    為了探究低層-中層大氣濾波效應(yīng)在中層-低熱層觀測(cè)的重力波的緯向傳播方向的影響,圖11(a—d)給出了與各個(gè)季節(jié)在6個(gè)觀測(cè)臺(tái)站的重力波事件相對(duì)應(yīng)的在12UT和18UT的平均緯向風(fēng)剖面(位置: 40°N,120°E).這里,緯向風(fēng)數(shù)據(jù)來(lái)自MERRA,其包含自地面至0.015 hPa 的72個(gè)垂直等壓層,空間分辨率為1/2° latitude×2/3° longitude.圖中表明,在20~70 km高度,夏季盛行西向風(fēng),冬季盛行東向風(fēng).然而風(fēng)的方向在春季和秋季季節(jié)發(fā)生逆轉(zhuǎn).且結(jié)果也表明,低層-中層大氣濾波效應(yīng)與重力波傳播方向僅在夏季吻合較好.正如 Li等(2011a)指出的低層-中層大氣濾波效應(yīng)僅能部分解釋重力波傳播方向各向異性,且波源的非均勻分布和波導(dǎo)傳播對(duì)重力波傳播方向各向異性也起到重要作用.

    圖5 觀測(cè)的所有重力波事件極向分布圖(a) Shuozhou; (b) Xinglong; (c) Donggang; (d) Xinxiang; (e) Linqu; (f) Rongcheng.Fig.5 Polar histograms of propagation directions of all gravity wave events

    圖6 春季重力波傳播方向極向圖(a) Shuozhou; (b) Xinglong; (c) Donggang; (d) Xinxiang; (e) Linqu; (f) Rongcheng.Fig.6 Polar histograms of propagation directions in spring

    圖7 夏季重力波傳播方向極向圖(a) Shuozhou; (b) Xinglong; (c) Donggang; (d) Xinxiang; (e) Linqu; (f) Rongcheng.Fig.7 Polar histograms of propagation directions in summer

    圖8 秋季重力波傳播方向極向圖(a) Shuozhou; (b) Xinglong; (c) Donggang; (d) Xinxiang; (e) Linqu; (f) Rongcheng.Fig.8 Polar histograms of propagation directions in autumn

    圖9 冬季重力波傳播方向極向圖(a) Shuozhou; (b) Xinglong; (c) Donggang; (d) Xinxiang; (e) Linqu; (f) Rongcheng.Fig.9 Polar histograms of propagation directions of all gravity wave events in winter

    圖10 降雨季節(jié)分布(a) 春季; (b) 夏季; (c) 秋季; (d) 冬季.Fig.10 The seasonal distribution of precipitation(a) Spring; (b) Summer; (c) Autumn; (d) Winter.

    圖11 緯向風(fēng)剖面圖(a) 春季; (b) 夏季; (c) 秋季; (d) 冬季.Fig.11 The profile of zonal wind velocity

    4結(jié)論

    本文基于分布在中國(guó)中緯地區(qū)2012年1月至2013年12月6個(gè)臺(tái)站的OH氣輝成像儀觀測(cè)數(shù)據(jù),研究分析了重力波的傳播特征.同時(shí),我們結(jié)合TRMM衛(wèi)星觀測(cè)、ECMWF和MERRA風(fēng)場(chǎng)數(shù)據(jù)討論了影響重力波傳播方向各向異性的因素.主要結(jié)果總結(jié)如下:

    (1) 重力波發(fā)生頻率在夏季和冬季高于春季和秋季季節(jié).

    (2) 水平波長(zhǎng)、觀測(cè)周期和水平相速度分別主要分布于10~35 km,4~12 min和30~100 m·s-1范圍.

    (3) 重力波傳播方向表現(xiàn)出明顯的季節(jié)變化.重力波在夏季主要沿極向方向傳播.然而,在冬季主要沿赤道方向和平行于赤道方向傳播.

    (4) 重力波在夏季和春季秋季的北向傳播趨勢(shì)可能主要由觀測(cè)臺(tái)站南方對(duì)流活動(dòng)導(dǎo)致.然而,對(duì)流層附近急流可能在決定冬季重力波主要傳播方向方面起到重要作用.

    波源的非均勻分布、臨界層的濾波效應(yīng)和導(dǎo)管傳播均在重力波傳播方向各向異性方面起重要作用.在未來(lái)的工作中,我們將結(jié)合風(fēng)場(chǎng)觀測(cè)和溫度數(shù)據(jù)對(duì)產(chǎn)生重力波傳播方向各向異性的因素進(jìn)行更深入的分析.

    致謝我們非常感謝ECMWF和MERRA提供水平風(fēng)場(chǎng)數(shù)據(jù).我們也感謝TRMM團(tuán)隊(duì)提供降雨數(shù)據(jù).同時(shí),我們也非常感謝中國(guó)子午工程提供數(shù)據(jù).

    References

    Alexander M J, Holton J R, Durran D R. 1995. The gravity wave response above deep convection in a squall line simulation.J.Atmos.Sci., 52(12): 2212-2226.

    Alexander M J, May P T, Beres J H. 2004. Gravity waves generated by convection in the Darwin area during the Darwin Area Wave Experiment.J.Geophys.Res., 109: D20S04, doi: 10.1029/2004JD004729.

    Baker D J, Stair A T. 1988. Rocket measurements of the altitude distributions of the hydroxyl airglow.Phys.Scripta, 37(4): 611-622.

    Buss S, Hertzog A, Hostettler C, et al. 2004. Analysis of a jet stream induced gravity wave associated with an observed ice cloud over Greenland.Atmos.Chem.Phys., 4: 1183-1200.

    Coble M R, Papen G C, Gardner C S. 1998. Computing two-dimensional unambiguous horizontal wavenumber spectra from OH airglow images.IEEET.Geosci.RemoteSensing, 36(2): 368-382.

    Dewan E M, Picard R H. 1998. Mesospferic bores.J.Geophys.Res., 103(D6): 6295-6305.

    Dou X K, Li T, Tang Y H, et al. 2010. Variability of gravity wave occurrence frequency and propagation direction in the upper mesosphere observed by the OH imager in Northern Colorado.J.Atmos.Sol.Terr.Phys., 72(5-6): 457-462.

    Ejiri M K, Shiokawa K, Ogawa T, et al. 2003. Statistical study of short-period gravity waves in OH and OI nightglow images at two separated sites.J.Geophys.Res., 108(D21): 4679, doi: 10.1029/2002JD002795.

    Fovell R, Durran D, Holton J R. 1992. Numerical simulations of convectively generated stratospheric gravity waves.J.Atmos.Sci., 49(16): 1427-1442.

    Garcia F J, Taylor M J, Kelley M C. 1997. Two-dimensional spectral analysis of mesospheric airglow image data.Appl.Optics, 36(29): 7374-7385.

    Hecht J H, Walterscheid R L, Hickey M P, et al. 2001. Climatology and modeling of quasi-monochromatic atmospheric gravity waves observed over Urbana Illinois.J.Geophys.Res., 106(D6): 5181-5195, doi: 10.1029/2000JD900722.

    Hecht J H, Kovalam S, May P T, et al. 2004. Airglow imager observations of atmospheric gravity waves at Alice Springs and Adelaide, Australia during the Darwin Area Wave Experiment (DAWEX).J.Geophys.Res., 109: D20S05, doi: 10.1029/2004JD004697.

    Hedin A E, Fleming E L, Monson A H, et al. 1996. Empirical wind model for the upper, middle and lower atmosphere.JournalofAtmosphericandTerrestrialPhysics, 58(13): 1421-1447.

    Horinouchi T, Nakamura T, Kosaka J I. 2002. Convectively generated mesoscale gravity waves simulated throughout the middle atmosphere.GeophysRes.Lett., 29(21): 3-1-3-4, doi: 10.1029/2002GL016069. Kim Y H, Lee C S, Chung J K, et al. 2010. Seasonal variations of mesospheric gravity waves observed with an airglow all-sky camera at Mt. Bohyun, Korea (36° N).J.Astron.SpaceSci., 27(3): 181-188, doi: 10.5140/JASS.2010.27.3.181. Kubota M, Fukunishi H, Okano S. 2001. Characteristics of medium- and large-scale TIDs over Japan derived from OI 630-nm nightglow observation.Earth,PlanetsSpace, 53(7): 741-751.

    Lane T P, Reeder M J, Clark T L. 2001. Numerical modeling of gravity wave generation by deep tropical convection.J.Atmos.Sci., 58(10): 1249-1274.

    Lane T P, Sharman R D, Clark T L, et al. 2003. An investigation of turbulence generation mechanisms above deep convection.J.Atmos.Sci., 60(10): 1297-1321, doi: 10.1175/1520-0469.

    Li Q, Xu J, Yue J, et al. 2011b. Statistical characteristics of gravity wave activities observed by an OH airglow imager at Xinglong, in northern China.Ann.Geophys., 29(8): 1401-1410.

    Li Z H, Liu A Z, Lu X, et al. 2011a. Gravity wave characteristics from OH airglow imager over Maui.J.Geophy.Res., 116: D22115, doi: 10.1029/2011JD015870.

    Medeiros A F, Taylor M J, Takahashi H, et al. 2003. An investigation of gravity wave activity in the low-latitude upper mesosphere: propagation direction and wind filtering.J.Geophys.Res., 108(D14): 4411. Medeiros A F, Buriti R A, Machado E A, et al. 2004a. Comparison of gravity wave activity observed by airglow imaging at two different latitudes in Brazil.J.Atmos.Sol.-Terr.Phys., 66(6-9): 647-654. Medeiros A F, Takahashi H, Batista P P, et al. 2004b. Observations of atmospheric gravity waves using airglow all-sky CCD imager at Cachoeira Paulista, Brazil (23°S, 45°W).GeofisicaInternational, 43(1): 29-39. Medeiros A F, Takahashi H, Buriti R A, et al. 2005. Atmospheric gravity wave propagation direction observed by airglow imaging in the South American sector.J.Atmos.Sol.-Terr.Phys., 67(17-18): 1767-1773. Nakamura T, Higashikawa A, Tsuda T, et al. 1999. Seasonal variations of gravity wave structures in OH airglow with a CCD imager at Shigaraki.EarthPlanetsSpace, 51(7-8): 897-906.

    Nakamura T, Aono T, Tsuda T, et al. 2003. Mesospheric gravity waves over a tropical convective region observed by OH airglow imaging in Indonesia.Geophys.ResLett., 30(17): 1882, doi: 10.1029/2003 GL017619.

    Narayanan V L, Gurubaran S. 2013. Statistical characteristics of high frequency gravity waves observed by OH airglow imaging from Tirunelveli (8.7°N).J.Atmos.Sol.-Terr.Phys., 92: 43-50.

    Pautet P D, Taylor M J, Liu A Z, et al. 2005. Climatology of short-period gravity waves observed over northern Australia during the Darwin area wave experiment (DAWEX) and their dominant source regions.J.Geophys.Res., 110: D03S90.

    Richter J H, Sassi F, Garcia R R. 2010. Toward a physically based gravity wave source parameterization in a general circulation model.J.Atmos.Sci., 67(1): 136-156.

    Shiokawa K, Suzuki S, Otsuka Y, et al. 2007. An intense gravity wave near the mesopause region observed by a Fabry-Perot interferometer and an airglow imager.J.Geophys.Res., 112: D07106, doi: 10.1029/2006JD007385.

    Song I S, Chun H Y, Lane T P. 2003. Generation mechanisms of convectively forced internal gravity waves and their propagation to the stratosphere.J.Atmos.Sci., 60(16): 1960-1980.

    Stockwell R G, Lowe R P. 2001a. Airglow imaging of gravity waves: 1. Results from a small network of OH nightglow scanning imagers.J.Geophys.Res., 106(D15): 17185-17203, doi: 10.1029/2001JD900035. Stockwell R G, Lowe R P. 2001b. Airglow imaging of gravity waves: 2. Critical layer filtering.J.Geophys.Res., 106(D15): 17205-17220, doi: 10.1029/2001JD900036.

    Suzuki S, Shiokawa K, Otsuka Y, et al. 2004. Statistical characteristics of gravity waves observed by an all-sky imager at Darwin, Australia.J.Geophys.Res., 109: D20S07, doi: 10.1029/2003JD004336. Suzuki S, Shiokawa K, Otsuka Y, et al. 2007. Gravity wave momentum flux in the upper mesosphere derived from OH airglow imaging measurements.EarthPlanetsSpace, 59(5): 421-428.

    Tang J, Swenson G R, Liu A Z, et al. 2005. Observational investigations of gravity wave momentum flux with spectroscopic imaging.J.Geophys.Res., 110: D09S09, doi: 10.1029/2004JD004778.

    Tang Y H, Dou X K, Li T, et al. 2014. Gravity wave characteristics in the mesopause region revealed from OH airglow imager observations over Northern Colorado.J.Geophys.Res., 119(1): 630-645, doi: 10.1002/2013JA018955. Taylor M J, Ryan E H, Tuan T F, et al. 1993. Evidence of preferential directions for gravity wave propagation due to wind filtering in the middle atmosphere.J.Geophys.Res., 98(A4): 6047-6057, doi: 10.1029/92JA02604.

    Taylor M J, Hapgood M A. 1988. Identification of a thunderstorm as a source of short period gravity waves in the upper atmospheric nightglow emissions.Planet.SpaceSci., 36(10): 975-985.

    Vadas S L, Fritts D C. 2001. Gravity wave radiation and mean responses to local body forces in the atmosphere.J.Atmos.Sci., 58(16): 2249-2279.

    Vadas S L, Fritts D C. 2009. Reconstruction of the gravity wave field from convective plumes via ray tracing.Ann.Geophys., 27(1): 147-177.

    Vadas S L, Yue J, She C Y, et al. 2009. A model study of the effects of winds on concentric rings of gravity waves from a convective plume near Fort Collins on 11 May 2004.J.Geophys.

    Walterscheid R L, Schubert G, Brinkman D G. 2001. Small-scale gravity waves in the upper mesosphere and lower thermosphere generated by deep tropical convection.J.Geophys.Res., 106(D23): 31825-31832, doi: 10.1029/2000JD000131.

    Wang C. 2010. New chains of space weather monitoring stations in China.SpaceWeather, 8: S08001, doi: 10.1029/2010SW000603.Wang C M, Li Q Z, Xu J Y, et al. 2014. Statistical characteristics analysis of atmospheric gravity waves with OH all-sky airglow imagers at low-latitude region of China.ChineseJ.Geophys. (in Chinese), 57(11): 3659-3667, doi: 10.6038/cjg20141120.

    Wu D L, Waters J W. 1996. Satellite observations of atmospheric variances: A possible indication of gravity waves.Geophys.Res.Lett., 23(24): 3631-3634.

    Wu Q, Killeen T L. 1996. Seasonal dependence of mesospheric gravity waves (<100 km) at Peach Mountain Observatory, Michigan.Geophys.Res.Lett., 23(17): 2211-2214.

    Yue J, Vadas S L, She C Y, et al. 2009. Concentric gravity waves in the mesosphere generated by deep convective plumes in the lower atmosphere near Fort Collins, Colorado.J.Geophys.Res., 114: D06104, doi: 10.1029/2008JD011244.

    附中文參考文獻(xiàn)

    王翠梅, 李欽增, 徐寄遙等. 2014. 基于OH全天空氣輝成像儀觀測(cè)的中國(guó)低緯地區(qū)的重力波傳播統(tǒng)計(jì)特征. 地球物理學(xué)報(bào), 57(11): 3659-3667, doi: 10.6038/cjg20141120.

    (本文編輯胡素芳)

    Gravity wave characteristics from multi-stations observation with OH all-sky airglow imagers over mid-latitude regions of China

    WANG Cui-Mei1,2, LI Qin-Zeng1, XU Ji-Yao1, YUAN Wei1

    1StateKeyLaboratoryofSpaceWeather,CenterforSpaceScienceandAppliedResearch,

    ChineseAcademyofSciences,Beijing100190,China2UniversityofChineseAcademyofSciences,Beijing100049,China

    AbstractGravity waves were observed by 6 OH airglow all-sky imagers at mid-latitude regions of China from January 2012 to December 2013. The result shows that the horizontal wavelengths, observed periods, and horizontal phase speeds were typically in the range of 10~35 km, 4~12 min, and 30~100 m/s, respectively. Gravity waves prefer to polarward propagation in summer. However, they tend to equatorward and parallel to the equator propagation during winter. Through TRMM satellite and ECMWF data sets, we find that the convective activity at the south of observation stations might mainly explain the northward propagation preference in summer. However, the jet streams near tropopause might play dominant role in determining preferred propagation direction of gravity waves during winter. Simultaneously, the result also suggests that the filtering effect of background wind from the low and middle atmospheric is matched well with propagation direction of AGWs just in summer in the mid-latitude of China in MLT region.

    KeywordsAirglow imager; Atmospheric gravity waves; Propagation; Wave source; Filtering effect

    基金項(xiàng)目國(guó)家自然科學(xué)基金(41229001,41331069,41404121)和國(guó)家重點(diǎn)基礎(chǔ)研究專(zhuān)項(xiàng)經(jīng)費(fèi)(2011CB811405)聯(lián)合資助.

    作者簡(jiǎn)介王翠梅,女,1989年生,中國(guó)科學(xué)院大學(xué)博士研究生,研究方向是中高層大氣. E-mail:cmwang@spaceweather.ac.cn

    doi:10.6038/cjg20160502 中圖分類(lèi)號(hào)P421 Res., 114: D06103, 10.1029/2008JD010753. Walterscheid R L, Hecht J H, Vincent R, et al. 1999. Analysis and interpretation of airglow and radar observations of quasi-monochromatic gravity waves in the upper mesosphere and lower thermosphere over Adelaide, Australia (35°S, 138°E). J. Atmos. Sol.-Terr. Phys., 61(6): 461-478, 10.1016/S1364-6826(99)00002-4.

    收稿日期2015-11-13,2016-02-25收修定稿

    王翠梅, 李欽增, 徐寄遙等. 2016. 基于多臺(tái)站OH全天空氣輝成像儀觀測(cè)的中國(guó)中緯地區(qū)重力波傳播特性.地球物理學(xué)報(bào),59(5):1566-1577,doi:10.6038/cjg20160502.

    Wang C M, Li Q Z, Xu J Y, et al. 2016. Gravity wave characteristics from multi-stations observation with OH all-sky airglow imagers over mid-latitude regions of China.ChineseJ.Geophys. (in Chinese),59(5):1566-1577,doi:10.6038/cjg20160502.

    猜你喜歡
    波源傳播
    多普勒效應(yīng)之新討論
    雙波源同向交替入射下大鼠的電磁暴露仿真
    中國(guó)文學(xué)作品外譯策略研究
    新媒體背景下湖湘文化的傳播效果評(píng)價(jià)研究
    淺論呂劇藝術(shù)的傳承與傳播
    戲劇之家(2016年19期)2016-10-31 17:25:42
    當(dāng)代傳播視野下的昆曲現(xiàn)象
    戲劇之家(2016年19期)2016-10-31 17:11:16
    新媒體環(huán)境下網(wǎng)絡(luò)輿情傳播
    當(dāng)前紙媒微信公眾號(hào)運(yùn)營(yíng)的突出問(wèn)題與策略建議
    基于運(yùn)動(dòng)學(xué)基礎(chǔ)知識(shí)的多普勒效應(yīng)公式推導(dǎo)
    低阻抗沖擊電壓標(biāo)準(zhǔn)波源
    大片电影免费在线观看免费| 精品久久久久久久久亚洲| 亚洲av电影在线观看一区二区三区 | av在线蜜桃| 日产精品乱码卡一卡2卡三| 欧美日韩视频精品一区| 日韩三级伦理在线观看| av在线播放精品| 久久久色成人| 久久99热6这里只有精品| 欧美3d第一页| 国产日韩欧美亚洲二区| 国产亚洲一区二区精品| 18禁在线播放成人免费| 亚洲精品视频女| 亚洲精品国产av成人精品| 日韩av免费高清视频| 国产精品久久久久久精品电影| 美女主播在线视频| 久久99精品国语久久久| av在线蜜桃| 麻豆久久精品国产亚洲av| 黄色一级大片看看| 中国美白少妇内射xxxbb| 久久精品国产亚洲网站| 亚洲婷婷狠狠爱综合网| 一级片'在线观看视频| 日本黄大片高清| 一级爰片在线观看| 亚洲av免费高清在线观看| 欧美少妇被猛烈插入视频| 高清av免费在线| 亚洲伊人久久精品综合| 赤兔流量卡办理| videossex国产| 国产一区二区三区av在线| 午夜视频国产福利| 亚洲国产成人一精品久久久| 欧美少妇被猛烈插入视频| 日韩av不卡免费在线播放| 99热这里只有精品一区| 男女国产视频网站| 男女边摸边吃奶| 直男gayav资源| 午夜福利视频精品| 久久热精品热| 国产成人a区在线观看| 欧美zozozo另类| 久久久久久九九精品二区国产| 亚洲国产精品成人久久小说| 欧美三级亚洲精品| 91狼人影院| 欧美97在线视频| av国产精品久久久久影院| 亚洲,一卡二卡三卡| 亚洲精品国产av成人精品| 久久精品久久久久久久性| 国产精品麻豆人妻色哟哟久久| 亚洲国产日韩一区二区| 免费观看a级毛片全部| 欧美精品国产亚洲| 国产精品久久久久久av不卡| 内地一区二区视频在线| 亚洲最大成人中文| eeuss影院久久| 一本色道久久久久久精品综合| 纵有疾风起免费观看全集完整版| 婷婷色综合大香蕉| 青青草视频在线视频观看| 亚洲欧美日韩另类电影网站 | 人妻夜夜爽99麻豆av| 欧美日韩一区二区视频在线观看视频在线 | 伊人久久精品亚洲午夜| 国产午夜精品一二区理论片| 国产永久视频网站| 欧美国产精品一级二级三级 | 欧美 日韩 精品 国产| 极品教师在线视频| 日韩免费高清中文字幕av| 国产午夜精品久久久久久一区二区三区| 日韩欧美一区视频在线观看 | 久久久久久久午夜电影| 国产国拍精品亚洲av在线观看| 日韩伦理黄色片| 国产精品国产三级国产专区5o| 嫩草影院新地址| 午夜免费鲁丝| 在线观看美女被高潮喷水网站| 91久久精品国产一区二区三区| 婷婷色综合大香蕉| 人人妻人人澡人人爽人人夜夜| 高清av免费在线| 高清日韩中文字幕在线| 欧美三级亚洲精品| 国产精品99久久久久久久久| 亚洲av成人精品一区久久| 少妇的逼水好多| 国产白丝娇喘喷水9色精品| 成人美女网站在线观看视频| 国内精品美女久久久久久| 精品一区二区免费观看| 99久国产av精品国产电影| 极品教师在线视频| 久久久久久久久久久免费av| 免费看a级黄色片| 有码 亚洲区| 国模一区二区三区四区视频| 欧美成人午夜免费资源| 免费看日本二区| 中国三级夫妇交换| 人人妻人人澡人人爽人人夜夜| 人妻系列 视频| 午夜老司机福利剧场| 国产精品久久久久久精品电影| 草草在线视频免费看| 亚洲婷婷狠狠爱综合网| 毛片一级片免费看久久久久| 久久久久久久精品精品| 国产视频内射| 欧美变态另类bdsm刘玥| 一区二区三区精品91| 99久久精品热视频| 91精品伊人久久大香线蕉| 国产精品一二三区在线看| 一级毛片我不卡| 国产成人福利小说| 久久久精品94久久精品| 亚洲内射少妇av| 亚洲欧洲日产国产| 精品一区在线观看国产| 亚洲怡红院男人天堂| 婷婷色综合大香蕉| 中文天堂在线官网| 午夜福利网站1000一区二区三区| 97超碰精品成人国产| 好男人在线观看高清免费视频| 久久久久网色| 一区二区三区乱码不卡18| 人人妻人人澡人人爽人人夜夜| 中文字幕免费在线视频6| 嫩草影院新地址| 日韩av在线免费看完整版不卡| 在线播放无遮挡| 亚洲精品乱码久久久v下载方式| 免费黄网站久久成人精品| 麻豆精品久久久久久蜜桃| 亚洲最大成人手机在线| 好男人视频免费观看在线| 国产成人精品婷婷| 大陆偷拍与自拍| 大陆偷拍与自拍| 国产成人精品福利久久| 精品久久久久久久人妻蜜臀av| 国产久久久一区二区三区| 国内精品宾馆在线| 亚洲,一卡二卡三卡| 久久精品综合一区二区三区| 国产永久视频网站| 色综合色国产| 久久久a久久爽久久v久久| 精品人妻一区二区三区麻豆| 免费观看的影片在线观看| 69人妻影院| av在线天堂中文字幕| 王馨瑶露胸无遮挡在线观看| 亚洲精品影视一区二区三区av| 国模一区二区三区四区视频| 最近手机中文字幕大全| 18禁在线播放成人免费| 免费高清在线观看视频在线观看| 国产老妇伦熟女老妇高清| 亚洲,一卡二卡三卡| 亚洲自拍偷在线| 99久久精品一区二区三区| 国产爽快片一区二区三区| 下体分泌物呈黄色| 91精品伊人久久大香线蕉| 亚洲成人中文字幕在线播放| 国产白丝娇喘喷水9色精品| av线在线观看网站| 免费黄色在线免费观看| 秋霞在线观看毛片| 国产成人精品婷婷| 亚洲最大成人中文| 亚洲四区av| 亚洲成人久久爱视频| av国产久精品久网站免费入址| 色视频在线一区二区三区| 色哟哟·www| 91久久精品国产一区二区三区| 精品酒店卫生间| 亚洲欧美精品自产自拍| 只有这里有精品99| 亚洲伊人久久精品综合| 久久久久久久久大av| 欧美xxⅹ黑人| 国产午夜精品一二区理论片| 麻豆成人av视频| www.av在线官网国产| 久久久午夜欧美精品| 国产成人精品婷婷| 亚洲欧洲国产日韩| 日韩制服骚丝袜av| 亚洲av免费在线观看| 你懂的网址亚洲精品在线观看| 国产精品久久久久久久电影| 国产精品福利在线免费观看| 精品国产乱码久久久久久小说| 看十八女毛片水多多多| 欧美日韩综合久久久久久| 国产高清不卡午夜福利| h日本视频在线播放| 成人亚洲精品一区在线观看 | 亚洲自拍偷在线| av在线观看视频网站免费| 国产国拍精品亚洲av在线观看| 久久97久久精品| 国产一区二区亚洲精品在线观看| 国产成年人精品一区二区| 亚洲内射少妇av| 亚洲美女视频黄频| 免费播放大片免费观看视频在线观看| 亚洲精品日韩av片在线观看| 交换朋友夫妻互换小说| 麻豆久久精品国产亚洲av| 精华霜和精华液先用哪个| 久久久成人免费电影| 高清日韩中文字幕在线| 国产伦精品一区二区三区视频9| 99热国产这里只有精品6| 一边亲一边摸免费视频| 我的女老师完整版在线观看| 日本-黄色视频高清免费观看| 麻豆精品久久久久久蜜桃| av线在线观看网站| 在线播放无遮挡| 男女那种视频在线观看| 内地一区二区视频在线| 国产色婷婷99| 亚洲人成网站在线播| 嘟嘟电影网在线观看| 国产欧美日韩一区二区三区在线 | 久久精品久久久久久久性| 王馨瑶露胸无遮挡在线观看| 美女国产视频在线观看| 交换朋友夫妻互换小说| 99re6热这里在线精品视频| 日韩在线高清观看一区二区三区| 精品久久久久久久人妻蜜臀av| 午夜老司机福利剧场| 热99国产精品久久久久久7| 国产老妇伦熟女老妇高清| 免费大片18禁| 欧美另类一区| 欧美日韩精品成人综合77777| 久久99精品国语久久久| 国产精品成人在线| 精品国产三级普通话版| 中文字幕av成人在线电影| 自拍偷自拍亚洲精品老妇| 成人毛片a级毛片在线播放| 大又大粗又爽又黄少妇毛片口| 久久精品久久久久久久性| 亚洲成人一二三区av| 亚洲国产欧美人成| 最近最新中文字幕免费大全7| 亚洲国产高清在线一区二区三| 大陆偷拍与自拍| 免费人成在线观看视频色| 国产精品麻豆人妻色哟哟久久| 免费看日本二区| 午夜老司机福利剧场| 青春草亚洲视频在线观看| 国产成人精品婷婷| 国产男人的电影天堂91| 日韩中字成人| 成人毛片a级毛片在线播放| 免费看光身美女| 成人综合一区亚洲| 国产精品久久久久久久久免| 精品人妻熟女av久视频| 日日摸夜夜添夜夜爱| 亚洲内射少妇av| 97人妻精品一区二区三区麻豆| 日日啪夜夜撸| h日本视频在线播放| 国产黄片美女视频| 黄色一级大片看看| 91精品一卡2卡3卡4卡| 大又大粗又爽又黄少妇毛片口| 亚洲三级黄色毛片| 69av精品久久久久久| 一级片'在线观看视频| 日韩大片免费观看网站| 国产高清三级在线| 精品人妻视频免费看| 久久久久久九九精品二区国产| 在线a可以看的网站| 久久久久国产网址| 亚洲国产欧美人成| 边亲边吃奶的免费视频| 欧美老熟妇乱子伦牲交| 亚洲va在线va天堂va国产| 亚洲,一卡二卡三卡| 日本av手机在线免费观看| 我要看日韩黄色一级片| 欧美精品国产亚洲| 欧美bdsm另类| 美女主播在线视频| 亚洲一区二区三区欧美精品 | 精品久久久精品久久久| 国产成人精品婷婷| 七月丁香在线播放| 欧美精品国产亚洲| 久久99热这里只有精品18| 久久久久久久久久人人人人人人| 亚洲成人久久爱视频| 18禁在线播放成人免费| 国产在线一区二区三区精| 男插女下体视频免费在线播放| 国内揄拍国产精品人妻在线| 街头女战士在线观看网站| 精品人妻视频免费看| 高清午夜精品一区二区三区| 黄色日韩在线| 日本黄大片高清| 亚洲欧美一区二区三区黑人 | 久久99热6这里只有精品| 97热精品久久久久久| 啦啦啦中文免费视频观看日本| 日本与韩国留学比较| 日韩人妻高清精品专区| 日韩成人伦理影院| 亚洲丝袜综合中文字幕| 国产又色又爽无遮挡免| 国产成人精品一,二区| 国产精品不卡视频一区二区| 夜夜爽夜夜爽视频| 国产一区有黄有色的免费视频| 亚洲av免费高清在线观看| 亚洲精品第二区| 日本欧美国产在线视频| 欧美一级a爱片免费观看看| 九九在线视频观看精品| 成人亚洲欧美一区二区av| 免费av不卡在线播放| 九九在线视频观看精品| 丝袜喷水一区| 99久久人妻综合| 麻豆成人午夜福利视频| 极品少妇高潮喷水抽搐| 少妇裸体淫交视频免费看高清| 亚洲av福利一区| 成人漫画全彩无遮挡| 男男h啪啪无遮挡| 免费看日本二区| 精品人妻一区二区三区麻豆| videos熟女内射| 国产精品不卡视频一区二区| 亚洲av二区三区四区| 搞女人的毛片| 成年女人在线观看亚洲视频 | 午夜激情福利司机影院| 久久99热6这里只有精品| 国产精品蜜桃在线观看| 美女被艹到高潮喷水动态| 日韩人妻高清精品专区| 爱豆传媒免费全集在线观看| 国产亚洲最大av| 久久人人爽av亚洲精品天堂 | av在线app专区| 最近中文字幕高清免费大全6| 又爽又黄无遮挡网站| 亚洲国产欧美在线一区| 狂野欧美激情性xxxx在线观看| 亚洲欧洲国产日韩| 日韩一本色道免费dvd| 91久久精品国产一区二区三区| 日韩制服骚丝袜av| 亚洲国产欧美在线一区| 国精品久久久久久国模美| 插阴视频在线观看视频| 久久精品国产亚洲av天美| 在线精品无人区一区二区三 | 久久久久性生活片| 午夜精品一区二区三区免费看| 91午夜精品亚洲一区二区三区| 青春草亚洲视频在线观看| 美女视频免费永久观看网站| 女人久久www免费人成看片| 精品一区二区免费观看| 午夜亚洲福利在线播放| 国产精品伦人一区二区| 大香蕉久久网| 亚洲激情五月婷婷啪啪| 亚洲精品,欧美精品| 亚洲av福利一区| 免费观看无遮挡的男女| 啦啦啦啦在线视频资源| 免费av观看视频| 久久99蜜桃精品久久| 日本爱情动作片www.在线观看| 成年人午夜在线观看视频| 成人亚洲精品av一区二区| 免费人成在线观看视频色| 国内精品美女久久久久久| 九色成人免费人妻av| 国产免费视频播放在线视频| 日本黄色片子视频| 亚洲av中文字字幕乱码综合| 免费看a级黄色片| 国产高清国产精品国产三级 | 一级毛片黄色毛片免费观看视频| 国产伦精品一区二区三区视频9| 国产精品嫩草影院av在线观看| 国产一区有黄有色的免费视频| 免费观看在线日韩| 亚洲av.av天堂| 日本欧美国产在线视频| 成人美女网站在线观看视频| 在线看a的网站| 亚洲av免费在线观看| 久久久久精品久久久久真实原创| 亚洲成人久久爱视频| 久久久精品免费免费高清| 只有这里有精品99| 国产高清有码在线观看视频| 久久久精品免费免费高清| 久久99蜜桃精品久久| 一个人观看的视频www高清免费观看| 亚洲精品成人av观看孕妇| 亚洲国产高清在线一区二区三| a级毛片免费高清观看在线播放| 日韩一区二区三区影片| 少妇高潮的动态图| 国产高清有码在线观看视频| 十八禁网站网址无遮挡 | 成人欧美大片| 嫩草影院精品99| 午夜免费观看性视频| 99九九线精品视频在线观看视频| 人人妻人人看人人澡| 偷拍熟女少妇极品色| 国产精品国产av在线观看| 午夜亚洲福利在线播放| 麻豆成人av视频| 91精品伊人久久大香线蕉| 国产欧美亚洲国产| 亚洲欧美清纯卡通| 欧美激情久久久久久爽电影| 丝袜美腿在线中文| 久久久久国产网址| 欧美成人一区二区免费高清观看| av天堂中文字幕网| 国产亚洲av嫩草精品影院| 寂寞人妻少妇视频99o| 国产中年淑女户外野战色| 久久精品国产亚洲网站| 午夜免费男女啪啪视频观看| 成人国产麻豆网| 欧美xxxx性猛交bbbb| 另类亚洲欧美激情| 久热久热在线精品观看| 在线天堂最新版资源| 神马国产精品三级电影在线观看| 日韩欧美精品免费久久| 久久亚洲国产成人精品v| 成人国产麻豆网| 乱系列少妇在线播放| 亚洲精品一二三| 最新中文字幕久久久久| 最后的刺客免费高清国语| 欧美人与善性xxx| a级毛片免费高清观看在线播放| 国产片特级美女逼逼视频| 天堂中文最新版在线下载 | 亚洲最大成人av| 国产综合懂色| 精品久久久久久久久av| 秋霞伦理黄片| 亚洲,一卡二卡三卡| 久久鲁丝午夜福利片| 亚洲国产最新在线播放| 精品国产露脸久久av麻豆| 欧美精品国产亚洲| 人妻 亚洲 视频| 色视频在线一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 美女国产视频在线观看| 18禁在线无遮挡免费观看视频| 人人妻人人澡人人爽人人夜夜| av女优亚洲男人天堂| 久久99蜜桃精品久久| 美女视频免费永久观看网站| 最新中文字幕久久久久| 国产欧美日韩一区二区三区在线 | 欧美xxxx黑人xx丫x性爽| 午夜免费观看性视频| 搡老乐熟女国产| av在线老鸭窝| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 高清视频免费观看一区二区| videossex国产| 免费黄网站久久成人精品| 身体一侧抽搐| 丰满少妇做爰视频| 成人亚洲欧美一区二区av| 欧美成人精品欧美一级黄| 极品教师在线视频| 全区人妻精品视频| 一区二区三区乱码不卡18| 伊人久久精品亚洲午夜| 中文欧美无线码| 69av精品久久久久久| 美女被艹到高潮喷水动态| 免费av毛片视频| 尤物成人国产欧美一区二区三区| 久久99蜜桃精品久久| 久久6这里有精品| 小蜜桃在线观看免费完整版高清| 3wmmmm亚洲av在线观看| 国产精品蜜桃在线观看| 边亲边吃奶的免费视频| 亚洲欧洲日产国产| 高清在线视频一区二区三区| 国产男女内射视频| 18+在线观看网站| 亚洲精品影视一区二区三区av| 国产高清不卡午夜福利| 韩国av在线不卡| 国产黄色免费在线视频| 又爽又黄a免费视频| 亚洲av成人精品一二三区| 久久久精品免费免费高清| 我要看日韩黄色一级片| 国产精品久久久久久久电影| 免费观看的影片在线观看| 日韩欧美精品v在线| 国产精品久久久久久精品古装| 亚洲av免费在线观看| 联通29元200g的流量卡| 卡戴珊不雅视频在线播放| 国产精品人妻久久久久久| 在线观看免费高清a一片| 欧美日本视频| 听说在线观看完整版免费高清| 国产av码专区亚洲av| a级毛片免费高清观看在线播放| av福利片在线观看| 69人妻影院| 亚洲电影在线观看av| 日韩亚洲欧美综合| 日韩伦理黄色片| 日日啪夜夜爽| 久久综合国产亚洲精品| 久久久久国产网址| 51国产日韩欧美| 久久久亚洲精品成人影院| 国产久久久一区二区三区| 成人国产av品久久久| 国产精品成人在线| 在线观看人妻少妇| 成人国产av品久久久| 免费黄频网站在线观看国产| av专区在线播放| 精品国产乱码久久久久久小说| 久久人人爽人人片av| 亚洲精品久久久久久婷婷小说| 丰满少妇做爰视频| 少妇高潮的动态图| 亚洲婷婷狠狠爱综合网| 欧美97在线视频| 国产真实伦视频高清在线观看| 久久精品久久久久久久性| 久久久午夜欧美精品| 久久精品久久久久久久性| 亚洲第一区二区三区不卡| 丝袜美腿在线中文| 午夜免费男女啪啪视频观看| 嫩草影院精品99| 夜夜爽夜夜爽视频| av免费观看日本| 99久久精品国产国产毛片| 国产探花极品一区二区| 午夜激情久久久久久久| 国产69精品久久久久777片| 国产黄频视频在线观看| 日韩亚洲欧美综合| 久热这里只有精品99| 久久久久久久久大av| 乱系列少妇在线播放| 久久精品熟女亚洲av麻豆精品| 少妇被粗大猛烈的视频| 国产亚洲av片在线观看秒播厂| 欧美日韩在线观看h| 在线观看免费高清a一片| 亚洲精品一二三| 中文精品一卡2卡3卡4更新| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美精品v在线| 国产精品99久久99久久久不卡 | 国国产精品蜜臀av免费| 18禁动态无遮挡网站| 免费大片黄手机在线观看| 久久久久久伊人网av| 一级av片app| 亚洲在久久综合| 久久影院123| 国产成人精品一,二区| 国产亚洲5aaaaa淫片| 在现免费观看毛片| 身体一侧抽搐| 2021少妇久久久久久久久久久| 人妻夜夜爽99麻豆av| 欧美xxxx性猛交bbbb|