• <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)波源
    18+在线观看网站| 如日韩欧美国产精品一区二区三区 | 国语对白做爰xxxⅹ性视频网站| 欧美bdsm另类| 人人妻人人爽人人添夜夜欢视频 | 欧美国产精品一级二级三级 | 欧美区成人在线视频| 国产黄片视频在线免费观看| 美女脱内裤让男人舔精品视频| 国产深夜福利视频在线观看| 精品少妇内射三级| 国产一区二区在线观看日韩| 国产永久视频网站| 成人综合一区亚洲| 亚洲美女视频黄频| 欧美变态另类bdsm刘玥| 国产精品久久久久久精品电影小说| 国产男人的电影天堂91| 99视频精品全部免费 在线| 久久精品国产亚洲av天美| 国产免费一区二区三区四区乱码| 国产男女内射视频| freevideosex欧美| 久久影院123| 国产av精品麻豆| 99热这里只有精品一区| 午夜老司机福利剧场| 国产欧美另类精品又又久久亚洲欧美| 日本-黄色视频高清免费观看| 久久国产精品男人的天堂亚洲 | 观看美女的网站| 成年美女黄网站色视频大全免费 | 亚洲va在线va天堂va国产| 久久精品国产自在天天线| 免费av不卡在线播放| 久久精品国产亚洲av涩爱| 极品少妇高潮喷水抽搐| 日韩亚洲欧美综合| 偷拍熟女少妇极品色| 桃花免费在线播放| 五月伊人婷婷丁香| 人人妻人人澡人人爽人人夜夜| 男人和女人高潮做爰伦理| 亚洲精品,欧美精品| 九草在线视频观看| 嘟嘟电影网在线观看| 性色av一级| 亚洲精品国产av蜜桃| 成人美女网站在线观看视频| 国产成人午夜福利电影在线观看| 欧美最新免费一区二区三区| 日韩精品免费视频一区二区三区 | 色哟哟·www| 亚洲精品日韩在线中文字幕| a级片在线免费高清观看视频| 国产日韩欧美在线精品| 亚洲av不卡在线观看| 最近中文字幕高清免费大全6| 亚洲精品色激情综合| 国产伦精品一区二区三区视频9| 在线天堂最新版资源| 精品一区二区三卡| 如何舔出高潮| 精品久久久久久久久av| 黄片无遮挡物在线观看| av福利片在线观看| 国产在线一区二区三区精| 精品一区二区三区视频在线| 亚洲精品久久午夜乱码| av免费在线看不卡| 国产黄频视频在线观看| 国产在线免费精品| 久久久欧美国产精品| 热re99久久精品国产66热6| 美女福利国产在线| 亚洲精品,欧美精品| 久久av网站| 国产乱人偷精品视频| 人妻少妇偷人精品九色| 国模一区二区三区四区视频| 男女边吃奶边做爰视频| 国产真实伦视频高清在线观看| 亚洲第一区二区三区不卡| 欧美变态另类bdsm刘玥| 午夜免费男女啪啪视频观看| 婷婷色av中文字幕| 下体分泌物呈黄色| 久久久久国产精品人妻一区二区| 免费看av在线观看网站| 中文字幕亚洲精品专区| 韩国av在线不卡| 美女cb高潮喷水在线观看| 午夜免费鲁丝| 免费av中文字幕在线| 女人精品久久久久毛片| 欧美日韩视频精品一区| 国产精品一区二区性色av| 少妇高潮的动态图| 国产伦理片在线播放av一区| 国产日韩欧美在线精品| 国产欧美日韩一区二区三区在线 | 美女国产视频在线观看| 久久热精品热| 亚洲中文av在线| 国产精品人妻久久久久久| 在线观看一区二区三区激情| 狠狠精品人妻久久久久久综合| 在线观看国产h片| av免费观看日本| 国内少妇人妻偷人精品xxx网站| 日韩在线高清观看一区二区三区| 久久精品夜色国产| 精品久久久久久久久av| 国产av精品麻豆| 久久精品国产亚洲网站| 免费久久久久久久精品成人欧美视频 | 欧美成人午夜免费资源| 午夜激情福利司机影院| 国产片特级美女逼逼视频| 丰满迷人的少妇在线观看| 国产精品免费大片| 精品少妇久久久久久888优播| 麻豆成人午夜福利视频| 欧美成人精品欧美一级黄| 最近2019中文字幕mv第一页| 精品久久久久久久久av| 成人无遮挡网站| 欧美成人精品欧美一级黄| 日韩电影二区| 不卡视频在线观看欧美| 有码 亚洲区| 91久久精品电影网| 三级国产精品欧美在线观看| 午夜91福利影院| 人体艺术视频欧美日本| 国产欧美日韩一区二区三区在线 | 成人毛片a级毛片在线播放| 精品人妻熟女av久视频| 一本大道久久a久久精品| 亚洲综合色惰| 亚洲内射少妇av| 各种免费的搞黄视频| 日本爱情动作片www.在线观看| h日本视频在线播放| 日韩视频在线欧美| 中文字幕免费在线视频6| 久久99热6这里只有精品| 狂野欧美激情性bbbbbb| 国产精品嫩草影院av在线观看| 99九九在线精品视频 | 曰老女人黄片| 中文字幕精品免费在线观看视频 | 日本欧美视频一区| 极品教师在线视频| 蜜桃久久精品国产亚洲av| 老女人水多毛片| 日韩在线高清观看一区二区三区| 大香蕉久久网| 夫妻性生交免费视频一级片| 老司机影院毛片| 久久久久久久久久人人人人人人| a级片在线免费高清观看视频| 熟女av电影| 91精品一卡2卡3卡4卡| a级一级毛片免费在线观看| 黄色欧美视频在线观看| 午夜影院在线不卡| 91成人精品电影| 极品少妇高潮喷水抽搐| 亚洲av成人精品一区久久| 十八禁高潮呻吟视频 | 黄色日韩在线| 国产精品99久久久久久久久| 国产精品不卡视频一区二区| 男女无遮挡免费网站观看| 涩涩av久久男人的天堂| 久久国产精品男人的天堂亚洲 | 国产精品.久久久| 日韩,欧美,国产一区二区三区| 伊人久久精品亚洲午夜| 成年人免费黄色播放视频 | 成人国产麻豆网| 我的女老师完整版在线观看| 一级毛片 在线播放| 18禁裸乳无遮挡动漫免费视频| 国产亚洲91精品色在线| 精品久久久久久久久av| 国产亚洲91精品色在线| 亚洲精品色激情综合| 成人国产麻豆网| 街头女战士在线观看网站| 久久99一区二区三区| 国产有黄有色有爽视频| 大香蕉久久网| 少妇高潮的动态图| a级片在线免费高清观看视频| 欧美变态另类bdsm刘玥| 赤兔流量卡办理| 丁香六月天网| 亚洲精品一区蜜桃| 日本黄色片子视频| 在线观看www视频免费| 少妇熟女欧美另类| 日本欧美视频一区| 男女无遮挡免费网站观看| 免费人成在线观看视频色| 大又大粗又爽又黄少妇毛片口| 日本欧美视频一区| 两个人免费观看高清视频 | 久久久久久久久久久久大奶| 国产女主播在线喷水免费视频网站| 人妻夜夜爽99麻豆av| 少妇被粗大猛烈的视频| 国产男人的电影天堂91| 免费av中文字幕在线| 亚洲精品色激情综合| 欧美bdsm另类| 最新中文字幕久久久久| 九九久久精品国产亚洲av麻豆| .国产精品久久| 久久久久久久久久人人人人人人| 你懂的网址亚洲精品在线观看| 男女免费视频国产| 国产黄色免费在线视频| 又大又黄又爽视频免费| 少妇人妻 视频| 美女主播在线视频| 99久久精品热视频| 少妇人妻久久综合中文| 精品一区二区免费观看| 黑人高潮一二区| 九色成人免费人妻av| 青春草亚洲视频在线观看| 最近中文字幕高清免费大全6| 好男人视频免费观看在线| 国产淫片久久久久久久久| 最近的中文字幕免费完整| 韩国高清视频一区二区三区| 午夜激情福利司机影院| 麻豆乱淫一区二区| 18禁动态无遮挡网站| 亚洲国产精品国产精品| av在线老鸭窝| 亚洲三级黄色毛片| 午夜视频国产福利| 久久久精品免费免费高清| 人人妻人人看人人澡| 日韩在线高清观看一区二区三区| 在线观看一区二区三区激情| 国产色爽女视频免费观看| 精品人妻一区二区三区麻豆| 久久久久久人妻| 麻豆精品久久久久久蜜桃| 午夜精品国产一区二区电影| 国产日韩欧美亚洲二区| 中文字幕久久专区| 亚洲精品一二三| 另类亚洲欧美激情| 亚洲久久久国产精品| videos熟女内射| av在线app专区| 男女无遮挡免费网站观看| 中国美白少妇内射xxxbb| 国产伦在线观看视频一区| 国产高清不卡午夜福利| 亚洲欧美精品自产自拍| 91午夜精品亚洲一区二区三区| 国产精品无大码| 精品国产乱码久久久久久小说| 少妇人妻 视频| 各种免费的搞黄视频| 少妇的逼好多水| 欧美日韩视频高清一区二区三区二| 久久综合国产亚洲精品| 王馨瑶露胸无遮挡在线观看| 一本一本综合久久| 国产亚洲一区二区精品| 久久精品国产a三级三级三级| 亚洲av成人精品一区久久| 精品视频人人做人人爽| 国产伦在线观看视频一区| 色5月婷婷丁香| 国产精品99久久久久久久久| 色94色欧美一区二区| 伊人亚洲综合成人网| 国产精品福利在线免费观看| 草草在线视频免费看| 啦啦啦在线观看免费高清www| 日韩成人av中文字幕在线观看| 在线观看免费高清a一片| 曰老女人黄片| 久久久久精品性色| 丝袜在线中文字幕| 极品教师在线视频| 老女人水多毛片| 午夜视频国产福利| 精品酒店卫生间| 精品人妻偷拍中文字幕| 欧美日韩精品成人综合77777| 一级毛片 在线播放| 又粗又硬又长又爽又黄的视频| 天堂8中文在线网| 日日摸夜夜添夜夜爱| 久久 成人 亚洲| 久久久久国产精品人妻一区二区| 99热这里只有是精品在线观看| 国产精品一二三区在线看| 久久久久久久久久久丰满| 午夜视频国产福利| 老司机影院毛片| av播播在线观看一区| 国产精品国产三级国产av玫瑰| 国产高清国产精品国产三级| 人妻人人澡人人爽人人| 欧美丝袜亚洲另类| 少妇精品久久久久久久| tube8黄色片| 久久久欧美国产精品| 午夜精品国产一区二区电影| 日本91视频免费播放| 欧美激情国产日韩精品一区| 免费观看av网站的网址| 91aial.com中文字幕在线观看| 亚洲国产av新网站| 久久久精品94久久精品| 少妇丰满av| 国产男人的电影天堂91| 18禁在线无遮挡免费观看视频| 麻豆精品久久久久久蜜桃| 少妇精品久久久久久久| 国产探花极品一区二区| 插逼视频在线观看| 观看av在线不卡| 成人特级av手机在线观看| 一本大道久久a久久精品| 中文字幕亚洲精品专区| 亚洲精华国产精华液的使用体验| 国产高清不卡午夜福利| 如日韩欧美国产精品一区二区三区 | 街头女战士在线观看网站| 色婷婷av一区二区三区视频| 多毛熟女@视频| 热re99久久国产66热| 最近2019中文字幕mv第一页| 国产色爽女视频免费观看| 五月开心婷婷网| 只有这里有精品99| 国产国拍精品亚洲av在线观看| 校园人妻丝袜中文字幕| 亚洲欧美一区二区三区国产| 国产精品麻豆人妻色哟哟久久| 色94色欧美一区二区| 街头女战士在线观看网站| 国产精品人妻久久久影院| a级毛片免费高清观看在线播放| 好男人视频免费观看在线| 天堂俺去俺来也www色官网| 精品亚洲成国产av| 精品酒店卫生间| 国产午夜精品久久久久久一区二区三区| 黑人高潮一二区| 好男人视频免费观看在线| 国产精品久久久久久久久免| 国产又色又爽无遮挡免| 亚洲欧美日韩东京热| 亚洲精品乱码久久久久久按摩| 日韩精品有码人妻一区| 久久久久人妻精品一区果冻| 男人和女人高潮做爰伦理| 亚洲美女搞黄在线观看| 高清视频免费观看一区二区| 久久免费观看电影| 日韩av免费高清视频| 日韩欧美精品免费久久| 特大巨黑吊av在线直播| 色5月婷婷丁香| 中文天堂在线官网| 亚洲一级一片aⅴ在线观看| 一区二区三区四区激情视频| 狂野欧美白嫩少妇大欣赏| 日本欧美视频一区| 老司机影院毛片| 亚洲欧洲国产日韩| av在线播放精品| 丰满饥渴人妻一区二区三| 人妻 亚洲 视频| 大话2 男鬼变身卡| 菩萨蛮人人尽说江南好唐韦庄| 啦啦啦中文免费视频观看日本| 国产69精品久久久久777片| av专区在线播放| 夜夜看夜夜爽夜夜摸| 26uuu在线亚洲综合色| 偷拍熟女少妇极品色| 97精品久久久久久久久久精品| 日本欧美视频一区| av.在线天堂| 免费不卡的大黄色大毛片视频在线观看| 18禁在线无遮挡免费观看视频| 久久久久久人妻| 波野结衣二区三区在线| 少妇丰满av| 黄色怎么调成土黄色| 亚洲精品成人av观看孕妇| tube8黄色片| 五月开心婷婷网| 亚洲欧洲国产日韩| 毛片一级片免费看久久久久| 久久久久久久国产电影| 黄色配什么色好看| 欧美bdsm另类| av国产精品久久久久影院| 久久99热6这里只有精品| 久久av网站| 日本av免费视频播放| 久久综合国产亚洲精品| 人妻系列 视频| 亚洲精品一二三| 国产精品国产三级专区第一集| 一级毛片电影观看| a 毛片基地| 欧美97在线视频| 欧美精品一区二区大全| 啦啦啦在线观看免费高清www| 中文字幕av电影在线播放| 亚洲不卡免费看| 国产高清不卡午夜福利| 亚洲综合色惰| 午夜福利影视在线免费观看| 乱码一卡2卡4卡精品| av网站免费在线观看视频| 成人国产av品久久久| 亚洲国产成人一精品久久久| 国产一区有黄有色的免费视频| 国产一区二区三区综合在线观看 | 美女福利国产在线| 婷婷色av中文字幕| 内射极品少妇av片p| 男女边摸边吃奶| 国产成人午夜福利电影在线观看| 欧美日韩综合久久久久久| 色婷婷av一区二区三区视频| 成人无遮挡网站| 亚洲av在线观看美女高潮| 欧美日韩一区二区视频在线观看视频在线| 国产欧美另类精品又又久久亚洲欧美| 国产在视频线精品| 精品视频人人做人人爽| 美女大奶头黄色视频| 成人特级av手机在线观看| 一级毛片我不卡| 精品亚洲乱码少妇综合久久| 91午夜精品亚洲一区二区三区| 在线观看美女被高潮喷水网站| 在线亚洲精品国产二区图片欧美 | 日韩视频在线欧美| 亚洲精品视频女| 欧美性感艳星| 国产欧美日韩一区二区三区在线 | 十八禁网站网址无遮挡 | 性色av一级| av福利片在线观看| 大码成人一级视频| 午夜激情福利司机影院| 亚洲性久久影院| 草草在线视频免费看| 国产精品久久久久成人av| a级一级毛片免费在线观看| 欧美精品一区二区免费开放| 精品少妇内射三级| 日本欧美视频一区| 国产黄频视频在线观看| 大片免费播放器 马上看| 精品久久久噜噜| 亚洲精品亚洲一区二区| 午夜免费鲁丝| 在线观看人妻少妇| 男女边摸边吃奶| av国产精品久久久久影院| 国产在线一区二区三区精| 寂寞人妻少妇视频99o| 一级片'在线观看视频| 亚洲精品成人av观看孕妇| 黄色一级大片看看| 亚洲综合色惰| 亚洲欧美一区二区三区国产| 亚洲国产日韩一区二区| 99久久中文字幕三级久久日本| 狂野欧美白嫩少妇大欣赏| 国产探花极品一区二区| 久久久久精品久久久久真实原创| 中文字幕精品免费在线观看视频 | 午夜av观看不卡| 99久国产av精品国产电影| 日韩一区二区视频免费看| 亚洲熟女精品中文字幕| 亚洲va在线va天堂va国产| 大又大粗又爽又黄少妇毛片口| freevideosex欧美| 黄色毛片三级朝国网站 | 欧美bdsm另类| 久久国内精品自在自线图片| av视频免费观看在线观看| 观看av在线不卡| 亚洲欧美清纯卡通| 人妻 亚洲 视频| 日本色播在线视频| 午夜91福利影院| 国产av一区二区精品久久| 久久国产乱子免费精品| 男男h啪啪无遮挡| 日本91视频免费播放| 日日啪夜夜撸| 一区二区三区免费毛片| 亚洲精品一二三| 十分钟在线观看高清视频www | 午夜影院在线不卡| 少妇人妻一区二区三区视频| 最近手机中文字幕大全| 日韩电影二区| 亚洲av成人精品一二三区| 最近的中文字幕免费完整| 欧美日韩综合久久久久久| 久久午夜综合久久蜜桃| 91成人精品电影| 亚洲真实伦在线观看| 卡戴珊不雅视频在线播放| 午夜视频国产福利| 菩萨蛮人人尽说江南好唐韦庄| 国产中年淑女户外野战色| av.在线天堂| 久久6这里有精品| 欧美变态另类bdsm刘玥| 日日撸夜夜添| 亚洲国产成人一精品久久久| 亚洲精品亚洲一区二区| 日本欧美国产在线视频| 日韩熟女老妇一区二区性免费视频| 午夜福利视频精品| 精品少妇黑人巨大在线播放| 成人漫画全彩无遮挡| 91成人精品电影| 亚洲欧美成人综合另类久久久| 麻豆乱淫一区二区| 亚洲av日韩在线播放| 美女内射精品一级片tv| 午夜福利,免费看| 国产精品三级大全| 日本欧美视频一区| 久久久久精品久久久久真实原创| 在线免费观看不下载黄p国产| 中文字幕人妻丝袜制服| av卡一久久| 亚洲性久久影院| 一级a做视频免费观看| 国产亚洲av片在线观看秒播厂| 午夜av观看不卡| 久久 成人 亚洲| 免费看不卡的av| 亚洲自偷自拍三级| 国产日韩欧美视频二区| 插逼视频在线观看| 男男h啪啪无遮挡| 在线精品无人区一区二区三| 激情五月婷婷亚洲| 午夜免费观看性视频| 自拍欧美九色日韩亚洲蝌蚪91 | 男女免费视频国产| 好男人视频免费观看在线| 成人综合一区亚洲| 国产精品秋霞免费鲁丝片| 亚洲国产最新在线播放| 一级毛片aaaaaa免费看小| 精品一区二区免费观看| 国产精品欧美亚洲77777| 久久久a久久爽久久v久久| 日韩不卡一区二区三区视频在线| 亚洲国产欧美日韩在线播放 | 精品久久久精品久久久| 在现免费观看毛片| 国产欧美亚洲国产| 噜噜噜噜噜久久久久久91| 99热这里只有是精品50| 十分钟在线观看高清视频www | 日本猛色少妇xxxxx猛交久久| 久久精品国产a三级三级三级| 91久久精品电影网| 亚洲成人手机| 一个人看视频在线观看www免费| 九九久久精品国产亚洲av麻豆| 中文乱码字字幕精品一区二区三区| 亚洲欧美日韩另类电影网站| kizo精华| 亚洲自偷自拍三级| 亚洲av不卡在线观看| 中文欧美无线码| 久久久久网色| 亚洲电影在线观看av| 午夜免费鲁丝| h视频一区二区三区| 精品午夜福利在线看| 建设人人有责人人尽责人人享有的| 精品国产一区二区三区久久久樱花| 免费av中文字幕在线| 丝袜喷水一区| 日本av免费视频播放| 精品卡一卡二卡四卡免费| 18+在线观看网站| 最近最新中文字幕免费大全7| 纯流量卡能插随身wifi吗| 久久久亚洲精品成人影院| 美女内射精品一级片tv| 亚洲成色77777| 成人18禁高潮啪啪吃奶动态图 | 国产熟女午夜一区二区三区 |