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

    極隙區(qū)緯度Pc1-2波的統(tǒng)計(jì)分布特征

    2014-03-08 01:25:17邱奕婷劉勇華趙浩峰
    極地研究 2014年3期
    關(guān)鍵詞:磁層中山站戴維斯

    邱奕婷 劉勇華 趙浩峰

    (1南京信息工程大學(xué),江蘇 南京210044;2中國(guó)極地研究中心,上海200136)

    0 引言

    相較于中低緯度,在高緯地區(qū)尤其是在極隙區(qū)附近,超低頻(ULF)波,也稱地磁脈動(dòng)或地磁微脈動(dòng),非?;钴S。這些波動(dòng)事件的發(fā)生往往與當(dāng)時(shí)的磁層活動(dòng)情況,如極區(qū)磁暴、亞暴、極光、磁層內(nèi)帶電粒子的能量通量變化,以及太陽(yáng)風(fēng)特性和太陽(yáng)活動(dòng)等一系列相關(guān)現(xiàn)象有密切關(guān)系。研究極區(qū)超低頻波現(xiàn)象可以獲得極區(qū)空間電磁環(huán)境的大量信息,對(duì)進(jìn)一步研究磁層物理和極地空間環(huán)境有重要價(jià)值。

    超低頻波的一個(gè)重要分支是Pc1-2波(0.1—5 Hz)。這是一種電磁離子回旋波,由離子回旋不穩(wěn)定性激發(fā),其頻率接近離子回旋頻率[1-2]。在地面,尤其是在高緯極隙區(qū)觀測(cè)到的Pc1-2波具有豐富形態(tài),根據(jù)是否出現(xiàn)重復(fù)的珠型結(jié)構(gòu)可將Pc1-2波分為結(jié)構(gòu)型和非結(jié)構(gòu)型兩種[3-5]。有研究顯示,結(jié)構(gòu)型Pc1-2波常在中低緯出現(xiàn),而非結(jié)構(gòu)型Pc1-2波則主要出現(xiàn)在高緯地區(qū)[3,6-9]。

    盡管在地球磁層的閉合磁力線上觀測(cè)到許多Pc1-2波,它們由從磁尾向日對(duì)流的環(huán)電流或等離子體片離子驅(qū)動(dòng),但一些高緯觀測(cè)臺(tái)站的觀測(cè)結(jié)果顯示 Pc1-2波或與極隙區(qū)有關(guān)[4,10-13]。Bolshakova等[14]利用在63°—85°地磁緯度范圍的 6個(gè)地面臺(tái)站的數(shù)據(jù)進(jìn)行統(tǒng)計(jì),結(jié)果表明當(dāng)極隙區(qū)的赤道向邊界非常接近觀測(cè)點(diǎn)的時(shí)候,Pc1-2波的發(fā)生率有最大值。此時(shí)的場(chǎng)線是穿過(guò)磁鞘和極隙區(qū)的,于是他們推測(cè)Pc1-2波在日側(cè)的磁鞘區(qū)域產(chǎn)生。

    雖然國(guó)外一些學(xué)者對(duì)Pc1-2波的季節(jié)、磁地方時(shí)、頻率和振幅變化已經(jīng)做過(guò)相關(guān)研究,如Kuwashima等統(tǒng)計(jì)了1977—1979年南極低緯和高緯兩個(gè)臺(tái)站Memambetsu和Syowa的Pc1波事件,發(fā)現(xiàn)在高緯臺(tái)站Pc1波集中發(fā)生在日側(cè),且在磁中午(~13 MLT)向后1 h左右有最大值,它最常被觀測(cè)到的時(shí)候?yàn)榇呵锛荆琒yowa站Pc1波的平均頻率峰值在0.3—0.4 Hz左右,但是文中并沒(méi)有涉及振幅隨磁地方時(shí)和季節(jié)的變化情況。Anderson等對(duì)1990年南極點(diǎn)站(South Pole)和麥克默多站 (McMurdo)以及1988年麥克默多站和Sondrestromfjord的磁流體(HM)哨聲(Pc1-2波的一種)數(shù)據(jù)進(jìn)行了研究,結(jié)果顯示這3個(gè)高緯臺(tái)站的磁流體哨聲也明顯地發(fā)生在日側(cè)靠近磁中午的時(shí)候,不過(guò)由于數(shù)據(jù)的不完整或缺失,他們從可用的數(shù)據(jù)中統(tǒng)計(jì)發(fā)現(xiàn)磁流體哨聲發(fā)生率在2、3月份有最大值,同時(shí)他們也并沒(méi)有研究Pc1-2波的頻率和振幅變化情況[15-16]。

    中國(guó)南極中山站(69.37°S,76.38°E)的建立為我們研究高緯地區(qū)地磁與高空大氣物理創(chuàng)造了極好的條件,它與澳大利亞南極戴維斯站(68.58°S,77.97°E)的直線距離約120 km,同位于磁層極隙區(qū)緯度,兩站均安裝了相同的感應(yīng)式磁力計(jì),可組成短基線研究空間ULF波,數(shù)據(jù)記錄在磁盤(pán),精度可達(dá)1 ms[17]。

    目前,國(guó)內(nèi)關(guān)于中山站Pc1-2波的傳播特性的相關(guān)研究,尤其是頻率和振幅方面的研究非常少。這里我們利用動(dòng)態(tài)譜和互譜分析方法,分析南極中山站和戴維斯站觀測(cè)的感應(yīng)式磁力計(jì)數(shù)據(jù),研究Pc1-2波的特性包括其頻次、振幅和頻率等波參數(shù)的時(shí)空變化特征,可以進(jìn)一步深入了解極隙區(qū)緯度Pc1-2波的傳播特性。

    1 數(shù)據(jù)及分析方法

    我們選取南極中山站和戴維斯站2004年3、6、9和12月共4個(gè)月觀測(cè)的感應(yīng)式磁力計(jì)數(shù)據(jù)(數(shù)據(jù)分別來(lái)自中國(guó)第20次南極科學(xué)考察隊(duì)和澳大利亞南極考察隊(duì))分析Pc1-2波現(xiàn)象。文中數(shù)據(jù)處理包括數(shù)據(jù)的格式轉(zhuǎn)化和譜分析。兩站的感應(yīng)式磁力計(jì)輸出為二進(jìn)制文件,包括南北(H)和東西(D)兩個(gè)分量,可通過(guò)軟件編程將其轉(zhuǎn)換為ASCII碼文件。完成格式轉(zhuǎn)換后,采用256點(diǎn)漢寧窗函數(shù),對(duì)數(shù)據(jù)進(jìn)行動(dòng)態(tài)譜分析來(lái)選擇Pc1-2波事件(圖1),當(dāng)動(dòng)態(tài)譜圖中0.1—5 Hz頻帶信號(hào)明顯增強(qiáng)且持續(xù)30 min以上時(shí),就初步判斷為有效Pc1-2波事件發(fā)生了,且所選事件為兩站同時(shí)發(fā)生的事件。事實(shí)上,對(duì)比兩站的Pc1-2波事件之后我們發(fā)現(xiàn)Pc1-2波在兩站都能被同時(shí)觀測(cè)到。然后對(duì)相應(yīng)的時(shí)間段數(shù)據(jù)進(jìn)行互譜分析,將頻率設(shè)定在0.1—5 Hz,當(dāng)最大振幅大于3倍的背景平均振幅值,且最大振幅大于0.01 nT2Hz時(shí),就獲得一個(gè)有效的波事件。

    圖1 中山站(a)和戴維斯站(b)2004年3月9日10:45—14:45 UT的Pc1-2波動(dòng)態(tài)譜:H分量(上幀)和D分量(下幀)Fig.1.The dynamic spectrum of Pc1-2 waves during 10:45—14:45 UT,9 March,2004 at Zhongshan Station(a)and Davis Station(b):H component(upper panel)and D component(lower panel)

    圖1是2004年3月9日10:45—14:45 UT南極中山站和戴維斯站觀測(cè)的Pc1-2波的動(dòng)態(tài)譜。從圖中可以看出在該事件中波的頻率范圍在0.1—0.35 Hz之間,屬于Pc1-2波,且該事件持續(xù)了近4小時(shí)左右。觀察所有波事件的動(dòng)態(tài)譜(文中以3月9日的波事件為例)可發(fā)現(xiàn)圖中動(dòng)態(tài)譜并沒(méi)有表現(xiàn)出“珠型”結(jié)構(gòu),也印證了在高緯地區(qū)出現(xiàn)的主要為非結(jié)構(gòu)型Pc1-2波。

    圖2為以 2004年 3月9日10:45—10:51 UT為例的兩站波事件原始數(shù)據(jù)記錄和相應(yīng)的經(jīng)濾波處理后的波形。從圖中也可看出,兩站的數(shù)據(jù)波形非常相似,濾波后波形周期在4—6 s左右,處在Pc1-2波的頻段范圍(0.2—10 s)。

    圖2 中山站(a)和戴維斯站(b)感應(yīng)式磁力計(jì)數(shù)據(jù):H分量原始記錄(第1幀,從上到下)、濾波后的H分量(第2幀)、D分量原始記錄(第3幀)和濾波后的D分量(第4幀)Fig.2.The data from the inductionmagnetometer at Zhongshan Station(a)and Davis Station(b):The original data(the first panel)and the filtered data(the second panel)on H component,the original data(the third panel)and the filtered data(the forth panel)on D component

    我們選擇互譜振幅最大值對(duì)應(yīng)頻率為在兩站傳播的波動(dòng)的最相關(guān)頻率即中心頻率,對(duì)應(yīng)的互譜相位為該頻率波動(dòng)在兩站之間傳播的相位差[18]。

    最后,統(tǒng)計(jì)事件的頻率和振幅的季節(jié)變化和磁地方時(shí)變化特征,將各天相同時(shí)段內(nèi)的波事件的中心頻率和最大互譜振幅求平均值并按磁地方時(shí)排列,就得到了在兩站傳播的Pc1-2波的頻率和振幅日變化。將每個(gè)月內(nèi)波事件的中心頻率和振幅求平均值并按季節(jié)順序排列則得到了Pc1-2波的頻率和振幅季節(jié)變化。

    2 結(jié)果和分析

    對(duì)南極中山站和戴維斯站2004年3、6、9和12月的感應(yīng)式磁力計(jì)數(shù)據(jù)做上述處理后,共篩選獲得2 932個(gè)Pc1-2波事件,并獲得Pc1-2波的發(fā)生頻次、中心頻率和振幅的季節(jié)和磁地方時(shí)變化特征。

    2.1 Pc1-2波頻次的季節(jié)和磁地方時(shí)變化

    圖3(a)為這4個(gè)月的Pc1-2波事件的發(fā)生頻次統(tǒng)計(jì),從圖中明顯可以看出在中山站-戴維斯站地區(qū)Pc1-2波發(fā)生頻次表現(xiàn)出顯著的季節(jié)變化。其中,在3月和9月觀測(cè)的波事件較多,分別觀測(cè)到1 506和768個(gè)事件,各占51.4%和26.2%;12月次之,為 546個(gè),占 18.6%;6月最少,為 112個(gè),占3.8%。

    圖3(b)給出了這4個(gè)月的Pc1-2波事件的磁地方時(shí)變化統(tǒng)計(jì)。在南極中山站和戴維斯站,磁地方時(shí)MLT≈UT+2 h。從圖中可以看出,在中山站和戴維斯站,Pc1-2波一般分布在0400—1500 UT這個(gè)范圍,且主要分布在午前(0800—1000 UT)靠近磁中午附近,有部分事件分布在晨側(cè) (0400—0700 UT)和午后 (1100—1500 UT)。

    2.2 Pc1-2波頻率的季節(jié)和磁地方時(shí)變化

    圖4為Pc1-2波的平均中心頻率的季節(jié)變化和磁地方時(shí)變化分布。從這些結(jié)果可以看出,Pc1-2波在6月(南極冬季)的平均中心頻率最大,且在午后靠近磁中午附近的平均中心頻率比午前的大。

    Y1的模型回歸方程為:Y1=-21405.4+2460.99x1+1021.11x2+227.592x3-193.636x12-7.40594x22-0.868309x32-157.577x1x2+5.18425x1x3-3.37819x2x3。

    2.3 Pc1-2波振幅的季節(jié)和磁地方時(shí)變化

    圖5為Pc1-2波的振幅季節(jié)變化和磁地方時(shí)變化分布。從這些結(jié)果可以看出,Pc1-2波的平均振幅在6月份最小,在10 MLT(0800 UT)出現(xiàn)最大值。

    圖3 2004年P(guān)c1-2波事件發(fā)生頻次的季節(jié)變化(左幀)和磁地方時(shí)變化(右?guī)〧ig.3.The variation of the Pc1-2 waves occurrencewith respect to the season(leftpanel)and themagnetic local time(rightpanel)in 2004

    圖4 Pc1-2波的平均中心頻率季節(jié)變化(左幀)和磁地方時(shí)變化(右?guī)┓植糉ig.4.The season(the left panel)and themagnetic local time(the right panel)variation of the average center frequency of the Pc1-2 waves

    圖5 Pc1-2波的振幅季節(jié)變化(左幀)和磁地方時(shí)變化(右?guī)┓植糉ig.5.The season(the left panel)and themagnetic local time(the right panel)variation of the amplitude of the Pc1-2 waves

    3 討論

    我們選取南極中山站和戴維斯站在2004年3、6、9和12月同時(shí)觀測(cè)的感應(yīng)式磁力計(jì)數(shù)據(jù),分別考察極隙區(qū)緯度Pc1-2波在南極秋、冬、春、夏四季的活動(dòng)特征,獲得了南極中山站-戴維斯站地區(qū)Pc1-2波發(fā)生頻次、波頻率和振幅的季節(jié)和地方時(shí)變化。

    圖3左幀示出的Pc1-2波的季節(jié)分布特征,有兩個(gè)明顯的特點(diǎn),一是南極冬季(6月)Pc1-2波的發(fā)生率很低,僅為3.8%,二是春秋季發(fā)生率高,分別達(dá)51.4%和26.2%,兩者相加達(dá)77.6%。這些結(jié)果與Kuwashima等[18]的結(jié)論一致,他們通過(guò)地面觀測(cè)也得出Pc1-2波事件在高緯地區(qū)春秋季出現(xiàn)頻率較高,認(rèn)為出現(xiàn)這種現(xiàn)象的原因可能與電離層電導(dǎo)率的變化有關(guān),6月份南極缺少光照,電離程度低,電離層電導(dǎo)率小,不利于Pc1-2波的傳播[15],故Pc1-2波事件在6月份較少。Natsuo Sato等[19]對(duì)非結(jié)構(gòu)型Pc1-2波在極光帶地磁共軛臺(tái)站的季節(jié)和日變化發(fā)生率特性也進(jìn)行過(guò)研究,結(jié)果顯示在地面臺(tái)站觀測(cè)到的非結(jié)構(gòu)型Pc1-2波強(qiáng)烈地受到電離層頂日光效應(yīng)的控制。在南北半球,光照可能會(huì)引起波從磁層通過(guò)電離層向地面?zhèn)鞑サ牟粚?duì)稱性,在日光照射的半球,由于電離層電子密度的增加,會(huì)產(chǎn)生沿地磁場(chǎng)線的電子密度梯度,且在夏季電子密度梯度增加,從而為Pc1-2波的產(chǎn)生和傳播創(chuàng)造條件。但是,這只能解釋南極冬季比夏季觀測(cè)的Pc1-2波少,而不能解釋南極春秋季具有最高的Pc1-2波發(fā)生率,因?yàn)橐话銇?lái)說(shuō)春秋季電離層的電導(dǎo)率會(huì)明顯小于夏季。Pc1-2波主要發(fā)生在春秋季,這與地磁擾動(dòng)在春秋季明顯增加的特點(diǎn)一致,說(shuō)明Pc1-2波的發(fā)生與磁層磁場(chǎng)擾動(dòng)程度具有一定的相關(guān)性。一般認(rèn)為,地磁擾動(dòng)水平在春秋季增強(qiáng),是由于地磁偶極軸對(duì)黃道面傾斜,磁層磁場(chǎng)與太陽(yáng)風(fēng)磁場(chǎng)的夾角隨地球公轉(zhuǎn)而變化,結(jié)果在春秋季節(jié)更加容易觸發(fā)太陽(yáng)風(fēng)磁層之間的磁重聯(lián),導(dǎo)致磁層磁場(chǎng)擾動(dòng)增加,此即所謂Russell-McPherron效應(yīng)[20]。而地磁活動(dòng)的半年變化與地球磁偶極軸傾角 (dipole tilt angle)有關(guān),Ables等[21-22]對(duì) 2008年斯科特站 (Scott Base)和凱西站(Casey)的Pc1-2波數(shù)據(jù)進(jìn)行了研究,他們的圖1所示結(jié)果也說(shuō)明了Pc1-2波在春秋季有最大發(fā)生率,且出現(xiàn)的時(shí)間在由地球磁偶極軸傾角引起的行星際磁場(chǎng)(IMF)的最大地磁效應(yīng)附近,也就是Russell-McPherron效應(yīng)的表現(xiàn),地球磁偶極軸傾角在春秋季比在冬夏季要小,這再一次說(shuō)明了日側(cè)磁重聯(lián)對(duì)Pc1-2波產(chǎn)生機(jī)制的重要性。Feygin和Kerttula等[23-24]對(duì) Pc1-2波的研究中也指出,Pc1波的動(dòng)態(tài)變化很可能是由與地磁活動(dòng)的發(fā)展相關(guān)的波動(dòng)源區(qū)位置的徑向發(fā)展造成的。

    另一方面,根據(jù)圖3右?guī)窘Y(jié)果,Pc1-2波發(fā)生率在磁中午附近呈現(xiàn)峰值,這可能是由于在磁中午兩站正處于軟電子沉降區(qū),電離層的電導(dǎo)率大,有利于Pc1-2波的傳播,從而使Pc1-2波在靠近磁中午附近發(fā)生較多[25]。我們注意到,在中山站地區(qū)觀測(cè)的Pc1-2波,主要發(fā)生在日側(cè),在夜側(cè)基本消失。這可能與我們選取波事件的要求有關(guān),即波動(dòng)持續(xù)30 min及以上,從而使得我們獲得的Pc1-2波事件非常典型。這可能排除了夜側(cè)與磁暴、亞暴活動(dòng)有關(guān)的Pc1-2波事件。這些長(zhǎng)時(shí)間持續(xù)發(fā)生的Pc1-2波主要發(fā)生在日側(cè),尤其是磁中午附近,表明其源區(qū)在磁中午附近。中山站地處極隙區(qū)緯度,L~14,靠近磁層頂。在白天,由于太陽(yáng)風(fēng)的壓縮,外磁層沿磁力線的磁場(chǎng)強(qiáng)度分布已經(jīng)不同于內(nèi)磁層,其極小值已經(jīng)遠(yuǎn)離赤道,在南北半球高緯區(qū)域形成了兩個(gè)磁場(chǎng)極小值區(qū)域。該結(jié)構(gòu)可以捕獲能量離子,從而形成Pc1-2波的激發(fā)源區(qū)[26],其中心位置就在磁中午附近。同時(shí),日側(cè)磁層容易受到太陽(yáng)風(fēng)的壓縮,這會(huì)增加被壓縮區(qū)域等離子體的溫度各項(xiàng)異性,從而增加 Pc1-2波發(fā)生的機(jī)會(huì)[27]。

    圖4和圖5顯示,中山站-戴維斯站地區(qū)Pc1-2波的平均頻率在冬季(6月)最大,而其平均振幅在冬季最小,這又是由于受到了電離層電導(dǎo)率的影響。如上所述,這些Pc1-2波起源于磁層尤其是磁中午附近。當(dāng)它們向下傳播到電離層時(shí),將發(fā)生模式轉(zhuǎn)換并繼續(xù)在電離層峰值電子濃度附近形成的波導(dǎo)中水平傳播[28-29]。受電離層波導(dǎo)結(jié)構(gòu)的影響,Pc1-2波在電離層中傳播時(shí)將受到最低截止頻率限制,即波頻率低于一定值時(shí),Pc1-2波不能在電離層中傳播[30]。而且,研究表明,該最低截止頻率與電離層電子濃度成反比[31]。在南極中山站-戴維斯站地區(qū),冬季處于極夜,缺少太陽(yáng)光照射,電離層電子濃度低,電導(dǎo)率低,因而Pc1-2波的截止頻率高,所以觀測(cè)到的Pc1-2波平均頻率高。

    Kim等[32]對(duì)高緯電離層波導(dǎo)中Pc1-2波的傳播特性進(jìn)行了統(tǒng)計(jì)研究,他們指出Pc1-2波的頻率變化在很大程度上依賴于電離層電導(dǎo)率,且波的能量衰減隨著頻率的增加而增加,因?yàn)椴ǖ哪芰克p跟碰撞頻率有關(guān)。這也就解釋了在冬季Pc1-2波頻率增加(圖 4)而振幅減小(圖 5)。同時(shí),Yahnin等[33]在北半球Sodankyla的地面觀測(cè)中也指出Pc1波的頻率在夏季有最小值,這與本文圖4的結(jié)論也是一致的,即在南半球中山站和戴維斯站Pc1-2波的頻率在冬季有最大值。這是因?yàn)橄鄬?duì)于較低頻率,較高頻率的波其波長(zhǎng)變小,從而會(huì)更易受到波導(dǎo)的限制,因此波的傳播對(duì)波導(dǎo)邊界的尺度更敏感,而波導(dǎo)的物理尺度是受電離層的日照條件影響的。因此,冬季受日照影響小,Pc1-2波傳播的頻率要高于夏季的Pc1-2波頻率。

    4 結(jié)論

    本文利用南極中山站和戴維斯站的感應(yīng)式磁力計(jì)測(cè)得的2004年3、6、9、12月的數(shù)據(jù),對(duì)極隙區(qū)緯度Pc1-2波事件進(jìn)行了統(tǒng)計(jì)分析,經(jīng)篩選,共得到2 932個(gè)事件。主要結(jié)果如下:

    (1)在中山站和戴維斯站,Pc1-2波事件出現(xiàn)在3、9月較多,分別有1 506和768個(gè),各占51.4%和26.2%;12月次之,為 546個(gè),占 18.6%;6月最少,為112個(gè),占3.8%;

    (2)Pc1-2波主要出現(xiàn)在午前(0800—1000 UT)靠近磁中午附近,有部分事件出現(xiàn)在晨側(cè)和午后;

    (3)Pc1-2波在6月(南極冬季)的平均中心頻率最大,且在午后靠近磁中午附近的平均中心頻率比午前的大;

    (4)Pc1-2波的平均振幅在6月份最小,在10 MLT(0800 UT)出現(xiàn)最大值。

    這些結(jié)果表明,極隙區(qū)緯度的Pc1-2波在很大程度上受極區(qū)季節(jié)變化和電離層電導(dǎo)率的影響。

    1 Cornwall JM.Cyclotron instabilities and electromagnetic emission in the ultra-low frequency and very low frequency ranges.Journal of Geophysical Research,1965,70(1):61—69.

    2 Fraser B J,Nguyen T S.Is the plasma pause a preferred source region of electromagnetic ion cyclotron waves in themagnetosphere.Journal of Atmospheric and Solar-Terrestrial Physics,2001,63(11):1225—1247.

    3 Fukunishi H,Toya T,Koike K,etal.Classification of hydromagnetic emissions based on frequency-time spectra.Journal of Geophysical Research,1981,86(A11):9029—9039.

    4 Menk FW,F(xiàn)raser B J,Hansen H J,etal.Identification of themagnetospheric cusp and cleftusing Pc 1-2 ULF pulsations.Journal of Atmospheric and Terrestrial Physics,1992,54(7-8):1021—1042.

    5 Menk FW,F(xiàn)raser B J,Hansen H J,et al.Multistation observations of Pc1-2 ULF pulsations in the vicinity of the polar cusp.Journal of Geomagnetism and Geoelectricity,1993,45(10):1159—1173.

    6 Fraser-Smith A C.Some statistics on Pc1 geomagneticmicropulsation occurrence atmiddle latitudes:inverse relation with sunspot cycle and semi-annual period.Journal of Geophysical Research,1970,75(25):4735—4745.

    7 Saito T.Geomagnetic pulsations.Space Science Reviews,1969,10(3):319—412.

    8 Gendrin R E,Troitskaya V A.Preliminary results of amicropulsation experimentat conjugate points.Radio Science,1965,69(8):1107—1116.

    9 Mursula K,Kangas J,Pikkarainen T.Properties of structured and unstructured Pc1 pulsations at high latitudes:variation over the 21st solar cycle.SolarWind Sources of Magnetospheric Ultra-Low-Frequency Waves,Geophysical Monograph Series,1994,81:409—415.

    10 Kangas J,Guglielmi A,Pokhotelov O.Morphology and physics of short-period magnetic pulsations-a review.Space Science Reviews,1998,83(3-4):435—512.

    11 Fraser B J,Loto’aniu TM,Singer H J.Electromagnetic ioncyclotron waves in themagnetosphere//Takahashi K,Chi P J,Denton R E,et al.Magnetospheric ULFWaves:Synthesis and New Directions.Geophysical Monograph Series.Washington:AGU,2006,169:195—212.

    12 Dyrud L P,Engebretson M J,Posch JL,et al.Ground observations and possible source regions of two types of Pc 1-2 micropulsations at very high latitudes.Journal of Geophysical Research,1997,102(A12):27011—27027.

    13 Engebretson M J,Moen J,Posch JL,etal.Searching for ULF signatures of the cusp:observations from search coilmagnetometers and auroral imagers in Svalbard.Journal of Geophysical Research,2009,114(A6):A06217.

    14 Bolshakova O V,Troitskaya V A,Ivanov K G.High-latitude Pc1-2 geomagnetic pulsations and their connection with location of the dayside polar cusp.Planetary and Space Science,1980,28(1):1—7.

    15 Kuwashima M,Toya T,Kawamura M,etal.Comparative study ofmagnetic Pc1 pulsations between low latitudes and high latitudes:statistical study.National Institute Polar Research Memoirs,1981,18:101—117.

    16 Anderson H M,Engebretson M J,Arnoldy R L,etal.Statistical study of hydromagnetic chorus events at very high latitudes.Journal of Geophysical Research,1995,100(A3):3681—3692.

    17 Fraser B J,McNabb PW,Menk FW,et al.A personal computer inductionmagnetometer system for recording geomagnetic pulsation.ANARE Re-search Notes,1991.

    18 劉勇華,劉瑞源,楊少峰,等.極隙區(qū)Pc3頻段脈動(dòng)的傳播特性.極地研究,2000,12(2):121—128.

    19 Sato N,Saemundsson T.Unstructured Pc1-2 pulsations observed at geomagnetically conjugate stations in the auroral zone.Journal of Geomagnetism and Geoelectricity,1990,42(5):653—662.

    20 Russell C T,McPherron R L.Semiannual variation of geomagnetic activity.Journal of Geophysical Research,1973,78(1):92—108.

    21 Nowada M,Shue JH,Russell C T.Effects of dipole tiltangle on geomagnetic activity.Planetary and Space Science,2009,57(11):1254—1259.

    22 Ables S,Sciffer M,F(xiàn)raser B J,etal.A study of polar cap Pc1-2:source locations and wave propagation//Proceedings of the9th Australian Space Science Conference.2010:105—114.

    23 Feygin F Z,Kleimenova N G,Pokhotelov O A,et al.Nonstationary pearl pulsations as a signature ofmagnetospheric disturbances.Annales Geophysicae,2000,18(5):517—522.

    24 Kerttula R,Mursula K,Pikkarainen T,et al.Effect ofmagnetic storm intensity on Pc1 activity at high and mid-latitudes.Journal of Atmospheric and Solar-Terrestrial Physics,2001,63(5):503—511.

    25 劉順林,賀龍松,劉瑞源.南極中山站冬季電離層的平均特性.極地研究,1997,9(3):192—197.

    26 Liu Y H,F(xiàn)raser B J,Menk FW.Pc2 EMIC waves generated high off the equator in the dayside outermagnetosphere.Geophysical Research Letters,2012,39(17):L17102,doi:10.1029/2012GL053082.

    27 Anderson B J,Hamilton D C.Electromagnetic ion cyclotron waves stimulated by modestmagnetospheric compressions.Journal of Geophysical Research,1993,98(A7):11369—11382.

    28 Tepley L,Landshoff R K.Waveguide theory for ionospheric propagation of hydromagnetic emissions.Journal of Geophysical Research,1966,71(5):1499—1504.

    29 Lysak R L.Propagation of Alfvén waves through the ionosphere:dependence on ionospheric parameters.Journal of Geophysical Research,1999,104(A5):10017—10030.

    30 Greifinger C,Greifinger P S.Theory of hydromagnetic propagation in the ionospheric waveguide.Journal of Geophysical Research,1968,73(23):7473—7490.

    31 Manchester R N.An Investigation of the Propagation of Hydromagnetic Emissions away from Their Generation Latitude to Both High and Low Latitudes.Australia:University of Newcastle,NSW,1968:116.

    32 Kim H,Lessard M R,Engebretson M J,etal.Statistical study of Pc1-2 wave propagation characteristics in the high-latitude ionospheric waveguide.Journal of Geophysical Research,2011,116(A7):A07227,doi:10.1029/2010JA016355.

    33 Yahnin A G,Yahnina T A,Demekhov A G,et al.Seasonal variations in the frequency of Pc1 geomagnetic pulsations and variations in the latitude of localized proton precipitations.Geomagnetizm and Aeronomiya,2004,44(3):282—287.

    猜你喜歡
    磁層中山站戴維斯
    南極中山站和長(zhǎng)城站日照時(shí)數(shù)的氣候變化特征及其影響因素
    基于對(duì)抗神經(jīng)網(wǎng)絡(luò)有限角度的磁層邊界CT重構(gòu)技術(shù)
    磁鞘快速流引起的磁層頂凹陷事件
    等離子體鞘層效應(yīng)對(duì)磁層探測(cè)電場(chǎng)儀設(shè)計(jì)的影響
    進(jìn)一步放大細(xì)膩通透的優(yōu)勢(shì) Davis Acoustics(戴維斯聲學(xué)) MODEL M/MODEL S/BASSON 88
    愛(ài)的就是這股Hi-Fi味 Davis Acoustics(戴維斯聲學(xué))Balthus 70
    地球磁層亞暴統(tǒng)計(jì)分析
    南極洲——中山站
    “雪龍”號(hào)搭載南極科考隊(duì)抵達(dá)中山站
    今日華人(2018年2期)2018-01-16 18:23:33
    太空藝術(shù)家唐·戴維斯
    太空探索(2016年2期)2016-07-12 09:57:27
    美女国产高潮福利片在线看| 999久久久国产精品视频| 欧美精品高潮呻吟av久久| 欧美人与善性xxx| 国产熟女欧美一区二区| 又大又黄又爽视频免费| 色综合欧美亚洲国产小说| 亚洲av电影在线观看一区二区三区| 无限看片的www在线观看| 国产精品蜜桃在线观看| 熟女av电影| 国产熟女欧美一区二区| 操美女的视频在线观看| 9191精品国产免费久久| 天堂俺去俺来也www色官网| 精品一品国产午夜福利视频| 精品福利永久在线观看| 亚洲成国产人片在线观看| 校园人妻丝袜中文字幕| 亚洲av综合色区一区| 国产av一区二区精品久久| 久久精品人人爽人人爽视色| 日韩精品有码人妻一区| 亚洲精品av麻豆狂野| 久久精品久久久久久久性| 国产av一区二区精品久久| 亚洲国产精品一区二区三区在线| 一级毛片 在线播放| 久久韩国三级中文字幕| 乱人伦中国视频| 啦啦啦 在线观看视频| 只有这里有精品99| 成人影院久久| 亚洲精品成人av观看孕妇| 男人添女人高潮全过程视频| 国产福利在线免费观看视频| av网站在线播放免费| 婷婷色综合www| 亚洲av成人不卡在线观看播放网 | 亚洲国产欧美在线一区| 丁香六月欧美| 狠狠精品人妻久久久久久综合| 亚洲国产精品999| 中文天堂在线官网| 69精品国产乱码久久久| 一区二区三区四区激情视频| 亚洲成人一二三区av| 国产不卡av网站在线观看| 天天操日日干夜夜撸| 操美女的视频在线观看| av不卡在线播放| 伊人久久国产一区二区| 女性生殖器流出的白浆| av福利片在线| 亚洲国产欧美网| 狠狠精品人妻久久久久久综合| 在线精品无人区一区二区三| 亚洲天堂av无毛| 欧美日韩视频高清一区二区三区二| 99热全是精品| 人妻一区二区av| 18在线观看网站| 亚洲国产中文字幕在线视频| 免费黄色在线免费观看| 亚洲视频免费观看视频| 91国产中文字幕| 日本av手机在线免费观看| 日韩中文字幕欧美一区二区 | 男女免费视频国产| 可以免费在线观看a视频的电影网站 | 国产精品蜜桃在线观看| 久久国产精品大桥未久av| 久久精品久久久久久噜噜老黄| 啦啦啦中文免费视频观看日本| 丰满饥渴人妻一区二区三| 国产一卡二卡三卡精品 | 欧美国产精品va在线观看不卡| 丁香六月天网| 制服丝袜香蕉在线| xxx大片免费视频| 亚洲自偷自拍图片 自拍| 男女午夜视频在线观看| 国产福利在线免费观看视频| av天堂久久9| 99热全是精品| 亚洲中文av在线| 国产极品粉嫩免费观看在线| 亚洲欧美一区二区三区久久| 丁香六月天网| 久久久亚洲精品成人影院| 夫妻午夜视频| 亚洲精品中文字幕在线视频| videosex国产| 中文字幕色久视频| 老司机亚洲免费影院| 亚洲欧美色中文字幕在线| 天堂中文最新版在线下载| 在线观看三级黄色| 97精品久久久久久久久久精品| 两个人免费观看高清视频| 亚洲av欧美aⅴ国产| 国产精品免费大片| 亚洲欧美一区二区三区国产| 黄片小视频在线播放| 黄片无遮挡物在线观看| 我的亚洲天堂| 男人舔女人的私密视频| 一级毛片黄色毛片免费观看视频| 亚洲精品国产一区二区精华液| 丝袜脚勾引网站| av视频免费观看在线观看| 日韩av在线免费看完整版不卡| 免费av中文字幕在线| 国产日韩欧美视频二区| 亚洲精品美女久久久久99蜜臀 | 丰满迷人的少妇在线观看| 免费人妻精品一区二区三区视频| 日本av手机在线免费观看| 中文字幕亚洲精品专区| 欧美激情高清一区二区三区 | 尾随美女入室| 精品国产乱码久久久久久小说| 国产1区2区3区精品| 熟妇人妻不卡中文字幕| 中文字幕人妻丝袜制服| 亚洲欧洲精品一区二区精品久久久 | 少妇精品久久久久久久| 亚洲精品自拍成人| 久久久久久久精品精品| 女的被弄到高潮叫床怎么办| 色网站视频免费| 成人国语在线视频| 亚洲欧美日韩另类电影网站| 精品少妇黑人巨大在线播放| 如日韩欧美国产精品一区二区三区| 青青草视频在线视频观看| 亚洲国产欧美日韩在线播放| 叶爱在线成人免费视频播放| 欧美少妇被猛烈插入视频| 天堂俺去俺来也www色官网| 成年人免费黄色播放视频| 视频区图区小说| 日韩免费高清中文字幕av| 国产免费又黄又爽又色| 国产精品一国产av| 亚洲欧美清纯卡通| 久久天堂一区二区三区四区| 久久久精品国产亚洲av高清涩受| 免费日韩欧美在线观看| 午夜福利,免费看| 99热国产这里只有精品6| 日韩欧美精品免费久久| 国产午夜精品一二区理论片| 久久精品久久久久久噜噜老黄| 亚洲综合精品二区| kizo精华| 国产成人免费观看mmmm| 一边摸一边做爽爽视频免费| 亚洲av日韩精品久久久久久密 | 国产激情久久老熟女| 一边摸一边做爽爽视频免费| 成人亚洲欧美一区二区av| 日日撸夜夜添| 亚洲三区欧美一区| 叶爱在线成人免费视频播放| 丰满迷人的少妇在线观看| 男人舔女人的私密视频| 欧美精品人与动牲交sv欧美| 国产99久久九九免费精品| av线在线观看网站| 亚洲精品自拍成人| 一级毛片我不卡| 久久精品熟女亚洲av麻豆精品| 一区二区三区精品91| 男男h啪啪无遮挡| 美女午夜性视频免费| 亚洲精品日韩在线中文字幕| 亚洲成人手机| 午夜免费男女啪啪视频观看| 极品人妻少妇av视频| 亚洲欧美中文字幕日韩二区| 欧美激情极品国产一区二区三区| 久久久久视频综合| 国产成人欧美| 纵有疾风起免费观看全集完整版| av在线播放精品| 国产成人系列免费观看| 日本av免费视频播放| 欧美日韩国产mv在线观看视频| 午夜免费观看性视频| 国产精品无大码| 亚洲欧美一区二区三区久久| 日韩一区二区视频免费看| 中文欧美无线码| 亚洲婷婷狠狠爱综合网| av视频免费观看在线观看| 97精品久久久久久久久久精品| 亚洲,欧美精品.| 国产极品粉嫩免费观看在线| 在线免费观看不下载黄p国产| 五月开心婷婷网| 亚洲av福利一区| 建设人人有责人人尽责人人享有的| 成人亚洲精品一区在线观看| 国产精品一区二区在线不卡| 亚洲伊人色综图| 欧美精品一区二区大全| 国产激情久久老熟女| 欧美 日韩 精品 国产| 欧美激情高清一区二区三区 | 成人亚洲精品一区在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲av在线观看美女高潮| 国产女主播在线喷水免费视频网站| 成人漫画全彩无遮挡| 伊人久久国产一区二区| 汤姆久久久久久久影院中文字幕| 午夜福利影视在线免费观看| 日韩,欧美,国产一区二区三区| 一本色道久久久久久精品综合| 一区二区三区激情视频| 男女免费视频国产| 色播在线永久视频| 熟女av电影| av片东京热男人的天堂| 麻豆乱淫一区二区| 香蕉丝袜av| 9色porny在线观看| 无遮挡黄片免费观看| 欧美亚洲 丝袜 人妻 在线| 丰满少妇做爰视频| 捣出白浆h1v1| 久久久久人妻精品一区果冻| 无遮挡黄片免费观看| 只有这里有精品99| 日本欧美国产在线视频| 亚洲av电影在线进入| 伦理电影免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 操出白浆在线播放| 一边摸一边做爽爽视频免费| 纵有疾风起免费观看全集完整版| 精品久久久精品久久久| 精品酒店卫生间| www.精华液| 天天躁夜夜躁狠狠躁躁| 成人手机av| 午夜久久久在线观看| 成年av动漫网址| 亚洲av中文av极速乱| 一区二区三区乱码不卡18| 久久亚洲国产成人精品v| 女人精品久久久久毛片| 久久人人97超碰香蕉20202| 美女大奶头黄色视频| 丝袜美足系列| 日韩,欧美,国产一区二区三区| 2018国产大陆天天弄谢| 涩涩av久久男人的天堂| 在线观看免费午夜福利视频| 黄网站色视频无遮挡免费观看| 亚洲第一区二区三区不卡| 久久久久久久精品精品| 大香蕉久久网| 国产片内射在线| videos熟女内射| 国产精品成人在线| 女人爽到高潮嗷嗷叫在线视频| 人妻一区二区av| 亚洲av福利一区| 久久人人97超碰香蕉20202| 午夜日本视频在线| a级毛片黄视频| 国产男女超爽视频在线观看| 国精品久久久久久国模美| 少妇被粗大的猛进出69影院| 欧美日韩亚洲国产一区二区在线观看 | 综合色丁香网| 亚洲成人国产一区在线观看 | 在现免费观看毛片| 国产精品av久久久久免费| 9色porny在线观看| 免费看不卡的av| 亚洲国产欧美在线一区| 成年av动漫网址| 免费高清在线观看视频在线观看| 九色亚洲精品在线播放| 亚洲av日韩在线播放| 精品少妇内射三级| 看免费成人av毛片| 欧美日韩视频高清一区二区三区二| 青春草视频在线免费观看| 亚洲国产av新网站| 熟妇人妻不卡中文字幕| 中文字幕人妻熟女乱码| 超碰97精品在线观看| 国产女主播在线喷水免费视频网站| 色播在线永久视频| 久久久久久久久久久免费av| 日日爽夜夜爽网站| 精品人妻熟女毛片av久久网站| 日韩中文字幕视频在线看片| 女人高潮潮喷娇喘18禁视频| 亚洲av成人精品一二三区| 亚洲成人一二三区av| 男女边摸边吃奶| 又黄又粗又硬又大视频| 国产高清不卡午夜福利| 亚洲婷婷狠狠爱综合网| 国产免费现黄频在线看| 国精品久久久久久国模美| a级片在线免费高清观看视频| 亚洲国产欧美在线一区| 欧美av亚洲av综合av国产av | av欧美777| 亚洲一区高清亚洲精品| 亚洲少妇的诱惑av| 亚洲人成伊人成综合网2020| 韩国精品一区二区三区| 女警被强在线播放| 在线观看www视频免费| 亚洲国产欧美一区二区综合| 精品福利观看| av网站免费在线观看视频| 一本综合久久免费| 日韩欧美免费精品| 99久久99久久久精品蜜桃| 啦啦啦 在线观看视频| 两个人看的免费小视频| 亚洲,欧美精品.| 日本在线视频免费播放| 精品人妻在线不人妻| 91麻豆av在线| 亚洲国产看品久久| 亚洲欧美精品综合久久99| 亚洲精品一区av在线观看| 操出白浆在线播放| 精品国产亚洲在线| 国产亚洲av嫩草精品影院| 侵犯人妻中文字幕一二三四区| 亚洲av成人不卡在线观看播放网| 精品无人区乱码1区二区| 一进一出抽搐gif免费好疼| 久久中文字幕人妻熟女| 搡老熟女国产l中国老女人| videosex国产| 国产成人精品久久二区二区免费| 亚洲av电影在线进入| 看片在线看免费视频| 99热只有精品国产| 精品第一国产精品| 亚洲av美国av| 欧美一级a爱片免费观看看 | 男女做爰动态图高潮gif福利片 | 午夜福利18| 免费在线观看黄色视频的| 国产精品一区二区精品视频观看| 国产精品 国内视频| 久久精品亚洲熟妇少妇任你| 国产高清有码在线观看视频 | 性欧美人与动物交配| 国产av在哪里看| 精品国产一区二区久久| 欧美+亚洲+日韩+国产| 可以在线观看的亚洲视频| 亚洲成人久久性| 禁无遮挡网站| 国产精品1区2区在线观看.| cao死你这个sao货| 国产精品1区2区在线观看.| 啪啪无遮挡十八禁网站| 久热爱精品视频在线9| 欧美在线一区亚洲| 9热在线视频观看99| 久久久久久久精品吃奶| 美女国产高潮福利片在线看| 亚洲精华国产精华精| 老司机在亚洲福利影院| 免费看十八禁软件| 动漫黄色视频在线观看| 国产熟女午夜一区二区三区| 午夜福利高清视频| 免费在线观看亚洲国产| 久久人人精品亚洲av| 国产成人一区二区三区免费视频网站| 日本五十路高清| 亚洲人成伊人成综合网2020| 久久人人精品亚洲av| 亚洲免费av在线视频| 午夜影院日韩av| 最新在线观看一区二区三区| 欧美在线黄色| 国产成人啪精品午夜网站| 久久亚洲精品不卡| 丰满人妻熟妇乱又伦精品不卡| 少妇被粗大的猛进出69影院| 亚洲情色 制服丝袜| 精品乱码久久久久久99久播| 欧美一级a爱片免费观看看 | 亚洲,欧美精品.| 国产精品久久电影中文字幕| 日韩欧美一区二区三区在线观看| 亚洲天堂国产精品一区在线| 久久精品国产清高在天天线| 午夜免费鲁丝| 琪琪午夜伦伦电影理论片6080| 免费一级毛片在线播放高清视频 | 亚洲午夜精品一区,二区,三区| 69av精品久久久久久| 老司机在亚洲福利影院| 嫁个100分男人电影在线观看| 亚洲 欧美 日韩 在线 免费| 亚洲一卡2卡3卡4卡5卡精品中文| 伊人久久大香线蕉亚洲五| 国产精品秋霞免费鲁丝片| 国产成人欧美| aaaaa片日本免费| 人人妻人人澡欧美一区二区 | 国产91精品成人一区二区三区| 久久人妻av系列| 国产精品二区激情视频| 国产激情欧美一区二区| 黄色视频,在线免费观看| 又黄又粗又硬又大视频| 老鸭窝网址在线观看| 九色亚洲精品在线播放| 欧美乱妇无乱码| 12—13女人毛片做爰片一| 黄色片一级片一级黄色片| 亚洲精品在线观看二区| 自拍欧美九色日韩亚洲蝌蚪91| 欧美在线黄色| 国内久久婷婷六月综合欲色啪| 久久中文看片网| 一级毛片女人18水好多| 亚洲人成77777在线视频| 咕卡用的链子| 欧美在线一区亚洲| 日韩中文字幕欧美一区二区| 别揉我奶头~嗯~啊~动态视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲成国产人片在线观看| av欧美777| 国产aⅴ精品一区二区三区波| 国产99白浆流出| 国产精品久久久人人做人人爽| 亚洲欧洲精品一区二区精品久久久| 男女床上黄色一级片免费看| 男男h啪啪无遮挡| 欧美激情高清一区二区三区| 九色亚洲精品在线播放| 好男人在线观看高清免费视频 | 久久精品人人爽人人爽视色| 一本综合久久免费| av超薄肉色丝袜交足视频| 亚洲中文字幕日韩| 国产熟女午夜一区二区三区| 国产亚洲精品第一综合不卡| 午夜成年电影在线免费观看| 免费高清视频大片| 国产熟女xx| 久久久久久久精品吃奶| 亚洲 欧美 日韩 在线 免费| 久久久精品国产亚洲av高清涩受| 美女午夜性视频免费| 级片在线观看| 女生性感内裤真人,穿戴方法视频| 91九色精品人成在线观看| 99国产精品一区二区蜜桃av| 亚洲欧美激情在线| 操出白浆在线播放| 久久国产精品男人的天堂亚洲| 国产高清videossex| 亚洲熟妇中文字幕五十中出| 精品久久久精品久久久| 中亚洲国语对白在线视频| 亚洲精品国产区一区二| 免费在线观看亚洲国产| 亚洲片人在线观看| 看免费av毛片| 国产亚洲精品久久久久5区| 777久久人妻少妇嫩草av网站| 成年人黄色毛片网站| 亚洲伊人色综图| 成熟少妇高潮喷水视频| 日本精品一区二区三区蜜桃| 欧美日本视频| 1024视频免费在线观看| 日韩免费av在线播放| 欧美中文日本在线观看视频| 日韩欧美一区二区三区在线观看| 亚洲精品av麻豆狂野| 免费女性裸体啪啪无遮挡网站| 亚洲精品一卡2卡三卡4卡5卡| 欧美色视频一区免费| av中文乱码字幕在线| 欧美另类亚洲清纯唯美| 国产欧美日韩精品亚洲av| av电影中文网址| 中亚洲国语对白在线视频| 一区二区三区高清视频在线| 最新在线观看一区二区三区| 1024香蕉在线观看| 热re99久久国产66热| 久久久久久久精品吃奶| 麻豆久久精品国产亚洲av| av片东京热男人的天堂| 丰满的人妻完整版| 精品高清国产在线一区| 一个人观看的视频www高清免费观看 | 欧美老熟妇乱子伦牲交| 亚洲中文日韩欧美视频| 日本 欧美在线| 欧美一区二区精品小视频在线| 少妇的丰满在线观看| 国产免费av片在线观看野外av| 窝窝影院91人妻| 亚洲欧美一区二区三区黑人| 亚洲中文字幕一区二区三区有码在线看 | 又紧又爽又黄一区二区| 国产成人一区二区三区免费视频网站| 亚洲片人在线观看| 黄色毛片三级朝国网站| 亚洲欧美一区二区三区黑人| 国产一级毛片七仙女欲春2 | 少妇熟女aⅴ在线视频| 青草久久国产| 天天躁狠狠躁夜夜躁狠狠躁| 91老司机精品| 成在线人永久免费视频| 亚洲午夜精品一区,二区,三区| 麻豆国产av国片精品| 国产精品久久久av美女十八| 黑人欧美特级aaaaaa片| 叶爱在线成人免费视频播放| 亚洲三区欧美一区| 十分钟在线观看高清视频www| 国产91精品成人一区二区三区| 色综合婷婷激情| 日本 欧美在线| 好男人电影高清在线观看| 欧美中文综合在线视频| 亚洲无线在线观看| 又紧又爽又黄一区二区| 欧美色欧美亚洲另类二区 | 免费一级毛片在线播放高清视频 | 色哟哟哟哟哟哟| 亚洲成国产人片在线观看| 亚洲男人天堂网一区| 中文字幕精品免费在线观看视频| 桃红色精品国产亚洲av| 亚洲美女黄片视频| 人人澡人人妻人| 免费女性裸体啪啪无遮挡网站| 妹子高潮喷水视频| 中文字幕色久视频| 久久精品影院6| 色老头精品视频在线观看| 欧美日韩一级在线毛片| 午夜a级毛片| 一级毛片高清免费大全| 午夜福利在线观看吧| av天堂在线播放| 国产日韩一区二区三区精品不卡| 国内毛片毛片毛片毛片毛片| 夜夜看夜夜爽夜夜摸| 亚洲成国产人片在线观看| 国产精品一区二区三区四区久久 | 激情视频va一区二区三区| 欧美黄色淫秽网站| 日本精品一区二区三区蜜桃| 女性被躁到高潮视频| 国产精品美女特级片免费视频播放器 | 欧美日韩福利视频一区二区| 在线观看一区二区三区| 真人做人爱边吃奶动态| 这个男人来自地球电影免费观看| 91成人精品电影| 午夜免费鲁丝| 色播亚洲综合网| 波多野结衣高清无吗| 日韩欧美国产在线观看| 精品国产乱码久久久久久男人| 不卡一级毛片| 精品久久久久久成人av| 成人亚洲精品一区在线观看| 久久香蕉精品热| 久久精品国产99精品国产亚洲性色 | 久久性视频一级片| 免费少妇av软件| 久久国产亚洲av麻豆专区| 1024视频免费在线观看| 亚洲自拍偷在线| 免费观看精品视频网站| 国产精品 国内视频| 日本 欧美在线| 涩涩av久久男人的天堂| 日韩大尺度精品在线看网址 | 老汉色∧v一级毛片| 亚洲男人的天堂狠狠| 露出奶头的视频| 亚洲电影在线观看av| 人妻久久中文字幕网| 国产精品爽爽va在线观看网站 | 亚洲精品在线观看二区| 国产精品一区二区免费欧美| 日本a在线网址| 亚洲 国产 在线| 最近最新中文字幕大全电影3 | 在线观看免费日韩欧美大片| 国产一区在线观看成人免费| 亚洲中文日韩欧美视频| 国产色视频综合| 国产精品香港三级国产av潘金莲| 黑人操中国人逼视频|