• <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
    老司机亚洲免费影院| 黑人巨大精品欧美一区二区蜜桃 | 亚洲一区二区三区欧美精品| 中国三级夫妇交换| 美女福利国产在线| 欧美日韩亚洲高清精品| 久久精品国产亚洲av涩爱| 80岁老熟妇乱子伦牲交| 高清黄色对白视频在线免费看| 如日韩欧美国产精品一区二区三区| 国产亚洲午夜精品一区二区久久| 熟女人妻精品中文字幕| 伊人久久国产一区二区| 女人久久www免费人成看片| 日本av免费视频播放| 激情五月婷婷亚洲| 久久人人爽人人片av| 成年动漫av网址| 精品熟女少妇av免费看| av片东京热男人的天堂| 丝袜美足系列| 国产亚洲精品第一综合不卡 | 亚洲精品国产av成人精品| 久久久久久久国产电影| 波野结衣二区三区在线| 美女国产视频在线观看| 久久国内精品自在自线图片| 国产又爽黄色视频| 久久久久久久久久久久大奶| 男女国产视频网站| a级毛色黄片| 午夜免费男女啪啪视频观看| 亚洲精品久久久久久婷婷小说| 国产激情久久老熟女| 99久久精品国产国产毛片| 亚洲伊人色综图| 女性生殖器流出的白浆| av免费在线看不卡| 色婷婷久久久亚洲欧美| 久久久精品免费免费高清| 国产福利在线免费观看视频| 九色成人免费人妻av| 国产精品熟女久久久久浪| 日日啪夜夜爽| 国产成人精品久久久久久| 最近最新中文字幕免费大全7| 国产免费一区二区三区四区乱码| 一个人免费看片子| 欧美精品一区二区大全| 在线观看国产h片| 亚洲人与动物交配视频| 精品人妻熟女毛片av久久网站| 国产精品偷伦视频观看了| 国产高清国产精品国产三级| 韩国高清视频一区二区三区| 成年人免费黄色播放视频| 男人操女人黄网站| 99久久综合免费| 韩国精品一区二区三区 | 成年人午夜在线观看视频| 亚洲欧美色中文字幕在线| 精品国产露脸久久av麻豆| 少妇人妻久久综合中文| 在线观看免费高清a一片| 午夜精品国产一区二区电影| 美女脱内裤让男人舔精品视频| 菩萨蛮人人尽说江南好唐韦庄| 大片免费播放器 马上看| 亚洲国产精品国产精品| 插逼视频在线观看| 有码 亚洲区| 国产精品一区www在线观看| 国产成人精品一,二区| 日本-黄色视频高清免费观看| 少妇的丰满在线观看| 亚洲欧美清纯卡通| 久热这里只有精品99| 天天躁夜夜躁狠狠久久av| 国产精品三级大全| 啦啦啦啦在线视频资源| 免费高清在线观看日韩| 又粗又硬又长又爽又黄的视频| 国产片特级美女逼逼视频| 美女中出高潮动态图| 欧美精品一区二区大全| 久久久久久久精品精品| 精品人妻一区二区三区麻豆| 伊人亚洲综合成人网| 免费av中文字幕在线| 亚洲欧美日韩另类电影网站| 伦精品一区二区三区| 日本午夜av视频| av在线播放精品| 天天躁夜夜躁狠狠躁躁| 日韩一区二区视频免费看| 午夜福利网站1000一区二区三区| 少妇人妻久久综合中文| 9191精品国产免费久久| 免费人成在线观看视频色| 母亲3免费完整高清在线观看 | 国产日韩一区二区三区精品不卡| 超碰97精品在线观看| 捣出白浆h1v1| 精品久久久久久电影网| 2021少妇久久久久久久久久久| 欧美97在线视频| 精品人妻偷拍中文字幕| 亚洲 欧美一区二区三区| 男人舔女人的私密视频| 插逼视频在线观看| 久久久久精品久久久久真实原创| 美国免费a级毛片| 天堂俺去俺来也www色官网| 久久狼人影院| 成人黄色视频免费在线看| 日日摸夜夜添夜夜爱| av线在线观看网站| 纵有疾风起免费观看全集完整版| 国产精品一区二区在线不卡| 韩国精品一区二区三区 | 久久久精品区二区三区| 成人亚洲欧美一区二区av| 亚洲成国产人片在线观看| 在现免费观看毛片| 在线精品无人区一区二区三| 桃花免费在线播放| 在线免费观看不下载黄p国产| 久久av网站| 69精品国产乱码久久久| 亚洲婷婷狠狠爱综合网| 亚洲人与动物交配视频| 韩国精品一区二区三区 | 成人漫画全彩无遮挡| 日本av手机在线免费观看| 亚洲欧美一区二区三区国产| 国产精品久久久久久精品古装| 免费少妇av软件| kizo精华| 美女脱内裤让男人舔精品视频| 哪个播放器可以免费观看大片| av片东京热男人的天堂| 亚洲一级一片aⅴ在线观看| 免费大片18禁| 美女xxoo啪啪120秒动态图| 亚洲av电影在线观看一区二区三区| 亚洲成人手机| 热99国产精品久久久久久7| 日韩欧美一区视频在线观看| 一区二区三区精品91| 亚洲欧美日韩卡通动漫| 日韩视频在线欧美| 26uuu在线亚洲综合色| 高清欧美精品videossex| av视频免费观看在线观看| 欧美精品一区二区免费开放| 亚洲欧美清纯卡通| 91精品三级在线观看| 国产成人免费观看mmmm| 亚洲欧美一区二区三区国产| 国产免费现黄频在线看| 亚洲精品国产色婷婷电影| 国产精品久久久久久久久免| 成人亚洲欧美一区二区av| 国产精品免费大片| 九九爱精品视频在线观看| 精品国产露脸久久av麻豆| 男女高潮啪啪啪动态图| 乱人伦中国视频| 永久免费av网站大全| 美女xxoo啪啪120秒动态图| 国产深夜福利视频在线观看| 大话2 男鬼变身卡| 久久热在线av| 午夜影院在线不卡| 街头女战士在线观看网站| 丝袜美足系列| 国产毛片在线视频| 国产欧美另类精品又又久久亚洲欧美| 啦啦啦啦在线视频资源| 亚洲欧美精品自产自拍| av国产精品久久久久影院| 精品亚洲成国产av| 观看av在线不卡| 最近手机中文字幕大全| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久99一区二区三区| 久久影院123| 亚洲图色成人| 国产精品无大码| 99九九在线精品视频| 久久精品aⅴ一区二区三区四区 | 人人妻人人澡人人看| 日本欧美国产在线视频| 久久精品人人爽人人爽视色| 国产精品女同一区二区软件| 国产高清三级在线| 如日韩欧美国产精品一区二区三区| 最后的刺客免费高清国语| 午夜视频国产福利| 女人精品久久久久毛片| 免费看av在线观看网站| 国产福利在线免费观看视频| 亚洲精品av麻豆狂野| 日本欧美视频一区| 视频区图区小说| 亚洲欧美中文字幕日韩二区| 99国产综合亚洲精品| 丰满乱子伦码专区| 天天影视国产精品| 波多野结衣一区麻豆| 久久国内精品自在自线图片| 日本黄色日本黄色录像| 国产精品欧美亚洲77777| 18禁观看日本| 成人手机av| 亚洲综合精品二区| 在线天堂中文资源库| 九色成人免费人妻av| 精品99又大又爽又粗少妇毛片| 欧美精品人与动牲交sv欧美| 国产精品久久久久久av不卡| 亚洲性久久影院| 精品国产乱码久久久久久小说| 毛片一级片免费看久久久久| 国产成人欧美| 日韩欧美一区视频在线观看| 狠狠婷婷综合久久久久久88av| 色婷婷av一区二区三区视频| 精品国产国语对白av| 久久精品国产a三级三级三级| 国产免费现黄频在线看| 精品久久久久久电影网| 性色avwww在线观看| 国产精品麻豆人妻色哟哟久久| 国产成人aa在线观看| 久久这里有精品视频免费| freevideosex欧美| 成人漫画全彩无遮挡| 女人被躁到高潮嗷嗷叫费观| 国产在线免费精品| 日韩不卡一区二区三区视频在线| 美女福利国产在线| 国产深夜福利视频在线观看| 欧美亚洲 丝袜 人妻 在线| 精品久久久精品久久久| 免费大片黄手机在线观看| 国产激情久久老熟女| 老熟女久久久| 国产精品一国产av| 青春草国产在线视频| 免费观看av网站的网址| 美女内射精品一级片tv| 精品人妻熟女毛片av久久网站| 日本猛色少妇xxxxx猛交久久| 中文字幕免费在线视频6| 国产永久视频网站| 99久国产av精品国产电影| 国产欧美亚洲国产| 色视频在线一区二区三区| 日韩精品免费视频一区二区三区 | 激情视频va一区二区三区| 成年美女黄网站色视频大全免费| 大片免费播放器 马上看| 亚洲性久久影院| 亚洲综合色网址| 成年av动漫网址| 黄片播放在线免费| 熟女电影av网| 国产熟女午夜一区二区三区| 久久久久久久精品精品| 丝瓜视频免费看黄片| 精品国产露脸久久av麻豆| 国产亚洲欧美精品永久| 日韩 亚洲 欧美在线| 尾随美女入室| 日韩大片免费观看网站| 老熟女久久久| 日韩不卡一区二区三区视频在线| 中文字幕制服av| 1024视频免费在线观看| 欧美精品亚洲一区二区| 久久99一区二区三区| 99久久综合免费| 欧美+日韩+精品| 一边亲一边摸免费视频| 精品亚洲成a人片在线观看| 欧美精品一区二区免费开放| 在线免费观看不下载黄p国产| 狠狠婷婷综合久久久久久88av| 亚洲精品av麻豆狂野| kizo精华| 中文字幕最新亚洲高清| 欧美精品高潮呻吟av久久| 国产一区二区在线观看日韩| 婷婷色av中文字幕| 蜜臀久久99精品久久宅男| 99久国产av精品国产电影| 久久久国产精品麻豆| 精品第一国产精品| 熟女av电影| 免费久久久久久久精品成人欧美视频 | 最近的中文字幕免费完整| 制服丝袜香蕉在线| 高清欧美精品videossex| 看免费成人av毛片| 精品久久蜜臀av无| 欧美日韩av久久| 国产亚洲精品第一综合不卡 | 久久午夜福利片| 国产日韩欧美在线精品| 免费高清在线观看视频在线观看| 在线观看国产h片| 中文字幕制服av| 一级,二级,三级黄色视频| 18禁在线无遮挡免费观看视频| 精品人妻在线不人妻| 两性夫妻黄色片 | 久久99热6这里只有精品| 人人妻人人爽人人添夜夜欢视频| 国产麻豆69| 中文字幕亚洲精品专区| 中文字幕人妻熟女乱码| 飞空精品影院首页| 欧美bdsm另类| 亚洲美女搞黄在线观看| 看非洲黑人一级黄片| 久久久精品区二区三区| 日本av免费视频播放| 麻豆精品久久久久久蜜桃| 亚洲av欧美aⅴ国产| 国产精品熟女久久久久浪| 全区人妻精品视频| 99久国产av精品国产电影| 欧美日韩一区二区视频在线观看视频在线| 最近手机中文字幕大全| 午夜福利在线观看免费完整高清在| 久久 成人 亚洲| 亚洲在久久综合| 搡老乐熟女国产| 久久毛片免费看一区二区三区| 国产男人的电影天堂91| 亚洲欧美色中文字幕在线| 国产精品国产三级国产av玫瑰| 亚洲图色成人| 女人精品久久久久毛片| 欧美精品高潮呻吟av久久| 最近最新中文字幕免费大全7| av黄色大香蕉| 国产成人免费观看mmmm| 久久久久久久亚洲中文字幕| 久久久国产一区二区| 精品一区二区三卡| 九色成人免费人妻av| 精品视频人人做人人爽| av一本久久久久| 久久精品aⅴ一区二区三区四区 | 高清不卡的av网站| 女人久久www免费人成看片| 黄色毛片三级朝国网站| 久久国产精品大桥未久av| 91精品伊人久久大香线蕉| 女人久久www免费人成看片| 中文字幕av电影在线播放| 国产探花极品一区二区| 国产av国产精品国产| 捣出白浆h1v1| 色视频在线一区二区三区| 99国产精品免费福利视频| 欧美性感艳星| 国产欧美亚洲国产| 亚洲高清免费不卡视频| 美女脱内裤让男人舔精品视频| 国语对白做爰xxxⅹ性视频网站| 欧美最新免费一区二区三区| 欧美日韩国产mv在线观看视频| 免费高清在线观看日韩| 我要看黄色一级片免费的| 在线天堂最新版资源| 我要看黄色一级片免费的| 精品熟女少妇av免费看| 色婷婷久久久亚洲欧美| 久久久亚洲精品成人影院| 97精品久久久久久久久久精品| 18在线观看网站| 亚洲精品第二区| 欧美日韩成人在线一区二区| 亚洲av综合色区一区| 欧美精品一区二区免费开放| 久久久久久久国产电影| 99久久中文字幕三级久久日本| 涩涩av久久男人的天堂| 国产国拍精品亚洲av在线观看| 精品国产一区二区三区四区第35| 免费av中文字幕在线| 18+在线观看网站| 欧美性感艳星| 伦理电影大哥的女人| 91精品国产国语对白视频| 成人免费观看视频高清| 婷婷色综合www| 亚洲av日韩在线播放| 中文乱码字字幕精品一区二区三区| 交换朋友夫妻互换小说| 韩国高清视频一区二区三区| 26uuu在线亚洲综合色| 18在线观看网站| 黄网站色视频无遮挡免费观看| 免费人妻精品一区二区三区视频| 色婷婷久久久亚洲欧美| 欧美最新免费一区二区三区| www.色视频.com| 老司机影院成人| 国产欧美亚洲国产| 日韩精品有码人妻一区| 亚洲精品国产av成人精品| 国产一区二区三区综合在线观看 | h视频一区二区三区| 精品一区二区三区四区五区乱码 | 内地一区二区视频在线| 视频在线观看一区二区三区| 亚洲精品国产av成人精品| 黄色一级大片看看| av一本久久久久| 日本欧美视频一区| 国产成人91sexporn| 飞空精品影院首页| 国产精品偷伦视频观看了| 亚洲人与动物交配视频| 制服丝袜香蕉在线| 一本久久精品| 久久国产精品大桥未久av| 日本av手机在线免费观看| 久久久欧美国产精品| 黑人高潮一二区| 中文字幕亚洲精品专区| 国产黄频视频在线观看| 99热6这里只有精品| 熟女电影av网| 久久 成人 亚洲| 国产成人一区二区在线| 日韩av免费高清视频| 91成人精品电影| 亚洲美女视频黄频| 天堂俺去俺来也www色官网| 久久久久久久久久久久大奶| 日韩欧美精品免费久久| 免费大片黄手机在线观看| 久久热在线av| 少妇被粗大的猛进出69影院 | 中文乱码字字幕精品一区二区三区| 久久国产精品大桥未久av| 2022亚洲国产成人精品| 男女边吃奶边做爰视频| 内地一区二区视频在线| 日本与韩国留学比较| 亚洲人与动物交配视频| 色网站视频免费| 91精品国产国语对白视频| 不卡视频在线观看欧美| 美女xxoo啪啪120秒动态图| 成人午夜精彩视频在线观看| 国产精品人妻久久久影院| 亚洲国产精品国产精品| 免费在线观看黄色视频的| 飞空精品影院首页| av福利片在线| 久久久a久久爽久久v久久| 国产亚洲午夜精品一区二区久久| 精品人妻熟女毛片av久久网站| 久久精品国产综合久久久 | 国产日韩欧美视频二区| 人人妻人人澡人人看| 日本欧美视频一区| 久久人人爽av亚洲精品天堂| 大陆偷拍与自拍| 黄色视频在线播放观看不卡| 五月伊人婷婷丁香| 免费黄频网站在线观看国产| 丰满乱子伦码专区| 性高湖久久久久久久久免费观看| 色婷婷久久久亚洲欧美| 黄片无遮挡物在线观看| 乱码一卡2卡4卡精品| 午夜久久久在线观看| 人人妻人人爽人人添夜夜欢视频| 一区二区三区乱码不卡18| a级毛片在线看网站| 国产精品久久久久久久久免| 桃花免费在线播放| 18+在线观看网站| 少妇猛男粗大的猛烈进出视频| 日本黄大片高清| 一边摸一边做爽爽视频免费| 人人妻人人澡人人看| 国产淫语在线视频| 亚洲色图综合在线观看| 亚洲国产精品一区二区三区在线| 亚洲av综合色区一区| 一级毛片 在线播放| 久久久欧美国产精品| 毛片一级片免费看久久久久| 母亲3免费完整高清在线观看 | 久久久久久久亚洲中文字幕| 国产精品女同一区二区软件| 全区人妻精品视频| 久久久精品区二区三区| 最近最新中文字幕免费大全7| 国产成人午夜福利电影在线观看| 亚洲高清免费不卡视频| 男女边摸边吃奶| 久久久久久人人人人人| 亚洲国产欧美在线一区| 免费看光身美女| 99久国产av精品国产电影| 欧美精品一区二区免费开放| 少妇精品久久久久久久| 天堂8中文在线网| 欧美日韩综合久久久久久| 亚洲精品成人av观看孕妇| 精品一区二区三区四区五区乱码 | xxx大片免费视频| 欧美激情国产日韩精品一区| 国产亚洲精品久久久com| 街头女战士在线观看网站| 一个人免费看片子| 综合色丁香网| 亚洲第一区二区三区不卡| 亚洲精品第二区| 少妇人妻久久综合中文| 欧美老熟妇乱子伦牲交| 日韩欧美精品免费久久| 一区二区三区精品91| 国产亚洲av片在线观看秒播厂| 亚洲av电影在线观看一区二区三区| 国产av精品麻豆| 亚洲五月色婷婷综合| 色婷婷av一区二区三区视频| 熟女av电影| 婷婷色综合www| 两性夫妻黄色片 | 亚洲欧美成人综合另类久久久| 国产淫语在线视频| 久久韩国三级中文字幕| 1024视频免费在线观看| 极品少妇高潮喷水抽搐| 久久综合国产亚洲精品| av.在线天堂| 欧美亚洲日本最大视频资源| 日韩中文字幕视频在线看片| 欧美日本中文国产一区发布| 丝袜脚勾引网站| 精品国产一区二区久久| 亚洲精品一二三| 婷婷色综合www| 欧美日韩视频精品一区| 亚洲国产看品久久| 日日爽夜夜爽网站| 交换朋友夫妻互换小说| 国产深夜福利视频在线观看| 国产精品一区二区在线不卡| 肉色欧美久久久久久久蜜桃| 亚洲国产欧美日韩在线播放| 国产精品三级大全| 黑丝袜美女国产一区| 性色av一级| 成人国语在线视频| 久久精品国产亚洲av涩爱| 曰老女人黄片| 精品一区二区三区四区五区乱码 | 黑人猛操日本美女一级片| 国产男女超爽视频在线观看| 激情视频va一区二区三区| 亚洲av福利一区| 日韩熟女老妇一区二区性免费视频| 日韩 亚洲 欧美在线| 亚洲av男天堂| 精品人妻熟女毛片av久久网站| 丝袜脚勾引网站| www.色视频.com| 亚洲成av片中文字幕在线观看 | 高清在线视频一区二区三区| 9色porny在线观看| 日本黄大片高清| 九九爱精品视频在线观看| 激情五月婷婷亚洲| 免费观看av网站的网址| 纯流量卡能插随身wifi吗| av一本久久久久| 欧美日韩一区二区视频在线观看视频在线| 黄色 视频免费看| 欧美国产精品一级二级三级| 国产av码专区亚洲av| 美女国产高潮福利片在线看| 在线天堂最新版资源| 如何舔出高潮| 伦理电影免费视频| 亚洲成色77777| 欧美丝袜亚洲另类| 国产精品免费大片| 人体艺术视频欧美日本| 欧美 日韩 精品 国产| 男女高潮啪啪啪动态图| 免费观看av网站的网址| 成年av动漫网址| 两个人看的免费小视频| 亚洲av日韩在线播放| 乱人伦中国视频| 天美传媒精品一区二区| 一级,二级,三级黄色视频| tube8黄色片| 午夜福利视频精品| 国产精品国产三级专区第一集| 国产xxxxx性猛交| videossex国产| 欧美人与性动交α欧美精品济南到 |