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

    基于CARISMA地磁臺鏈觀測的IPDP事件的統(tǒng)計分析

    2015-03-16 10:50:09周銘霞袁志剛李海夢王德棟余雄東王貞貞喬錚
    地球物理學報 2015年7期
    關鍵詞:質(zhì)子臺站等離子體

    周銘霞, 袁志剛, 李海夢, 王德棟, 余雄東, 王貞貞, 喬錚

    武漢大學電子信息學院, 武漢 430072

    ?

    基于CARISMA地磁臺鏈觀測的IPDP事件的統(tǒng)計分析

    周銘霞, 袁志剛*, 李海夢, 王德棟, 余雄東, 王貞貞, 喬錚

    武漢大學電子信息學院, 武漢 430072

    亞暴期間磁尾等離子體片離子注入內(nèi)磁層能夠激發(fā)電磁離子回旋(EMIC)波.對應于這種EMIC波,地面磁力儀可觀測到周期逐漸減小的地磁脈動(IPDP).利用GOES衛(wèi)星數(shù)據(jù),地磁指數(shù)和加拿大CARISMA地磁臺站的數(shù)據(jù),我們研究了IPDP事件的產(chǎn)生與亞暴磁尾注入的關系.同時利用CARISMA地磁臺鏈中的MCMU和MSTK兩個臺站,從2005年4月到2014年5月期間的觀測數(shù)據(jù),統(tǒng)計分析了亞暴期間的IPDP事件,研究了IPDP事件的出現(xiàn)率關于季節(jié)和磁地方時的分布特征.我們總共獲得128個兩個臺站同時觀測的IPDP事件.該類事件關于季節(jié)分布的發(fā)生率,冬季最小,為13.28%,春季最大,為32.81%,結果表明IPDP事件關于季節(jié)分布的發(fā)生率受到電離層電導率及亞暴發(fā)生率的影響.兩個臺站同時觀測到的IPDP事件最大出現(xiàn)率出現(xiàn)在15—18 MLT(磁地方時),結果表明IPDP事件主要由亞暴期間產(chǎn)生的能量離子注入內(nèi)磁層,西向漂移遇到等離子體層羽狀結構(Plume)區(qū)的高密度等離子體所激發(fā).

    電磁離子回旋波; 周期逐漸減小的地磁脈動; 亞暴注入; 發(fā)生率

    1 引言

    電磁離子回旋(EMIC)波,由高能離子溫度各向異性(T⊥>T∥)導致的離子回旋不穩(wěn)定性所激發(fā).EMIC波的頻帶范圍為0.1~5 Hz的是Pc1-2波 (Halford et al., 2010).EMIC波在快速離子加熱,冷電子加熱,投擲角散射和沉降導致的環(huán)電流衰減和輻射帶相對論電子快速損失方面起到非常重要的作用 (Jordanova et al., 2001, Fok et al., 2005, Fraser et al., 2012, Su et al., 2010, 2011, Xiao et al., 2011, Yuan et al., 2013, 2014, Zhou et al., 2013, Wang et al., 2014),因此研究EMIC在磁層中的表現(xiàn)特征和發(fā)生率非常重要.在磁暴期間,等離子體層或等離子體羽狀結構(plume)的高密度冷等離子體,降低產(chǎn)生環(huán)電流離子回旋不穩(wěn)定性的門限,從而有利于激發(fā)EMIC波 (Cornwall, 1965, Fraser et al., 2005, Liu et al., 2012, Yuan et al., 2012a, 2012b).等離子體層plume形成于磁暴期間等離子體層黃昏側,而且可從黃昏側向日側磁層頂延伸.不僅等離子體層plume還有等離子體層頂這些高密度梯度區(qū)域都有利于EMIC波的產(chǎn)生與傳播 (Yuan et al., 2010, Yuan et al., 2012a).磁尾等離子體注入的高能質(zhì)子西向漂移遇到等離子體層或等離子體層黃昏隆起時,易激發(fā)EMIC波(Hayakawa et al., 1992).在冷等離子體里沿著地磁場平行傳播的EMIC波的色散關系表明,相互作用的高能質(zhì)子的平行能量在磁赤道有一個最小值,所以磁赤道區(qū)域有利于EMIC波的產(chǎn)生 (Lee and Kwok, 1984, Loto′aniu et al., 2005).磁赤道激發(fā)的EMIC波可沿磁力線傳播到電離層,從而被地磁臺站觀測到Pc1-2波動(Lee et al., 1984, Loto′aniu et al., 2005).

    根據(jù)Pc1-2在功率譜中的形態(tài),Pc1-2可大致分為兩類.第一類是結構型或珠型脈動,這一類事件都具有周期性調(diào)制波幅, 其他大多事件歸為第二類,即非結構型脈動 (Kerttula et al., 2001).周期逐漸減小的地磁脈動(IPDP)就是典型的非結構型的Pc1-2現(xiàn)象.它最顯著的特征在于它的頻率漂移,而且這個特點的出現(xiàn)與亞暴的發(fā)生緊密相關(Saito, 1969, Gendrin, 1971, Gul′elmi, 1974).在功率譜圖中IPDP事件存在著明顯的頻率增長,而且占據(jù)一個比較寬的頻帶.IPDP通常在亞暴開始后的下午或傍晚被觀測到 (Lee and Kwok, 1984).通常 IPDP事件發(fā)生之前,都存在一個AE指數(shù)的增長,并且IPDP波事件的開始時間與AE指數(shù)開始增長的時刻相差0.5~2 h (Hayakawa et al., 1992).以往的研究只有少數(shù)短期觀測的統(tǒng)計分析,研究主要集中在IPDP事件的事例分析.本文選取了MCMU和MSTK兩個臺站從2005年4月到2014年5月這10年間同時觀測到的IPDP事件做了統(tǒng)計分析,研究IPDP事件的發(fā)生率關于磁地方時(MLT)及季節(jié)分布的統(tǒng)計特征.

    2 探測儀器

    Canadian Array for Realtime Investigation of Magnetic Activity(CARISMA)是加拿大地球空間監(jiān)測系統(tǒng)工程(Canadian Geospace Monitoring)的磁力儀監(jiān)測部分.它是一個位于北美大陸,用于測量地磁場的磁力儀設備陣列.我們選取了CARISMA臺鏈中的Fort McMurray(MCMU) 和Ministik Lake(MSTK)兩個臺站的磁通門磁力儀(FGM)的觀測數(shù)據(jù),它們的地理位置分別為(56.65°N,248.79°E)和(53.351°N,247.026°E),L值分別為5.35和4.22.這兩個臺站在磁赤道平面的投影一般處于等離子體層頂附近,非常有利于觀測等離子體層頂附近的EMIC波事件.我們利用FGM的高分辨率(8 Hz)的數(shù)據(jù)進行了IPDP事件分析.

    Geostationary Operational Environmental Satellite(GOES)衛(wèi)星位于地球同步軌道(L=6.6),GOES衛(wèi)星攜帶的高能粒子傳感器(Energetic Particle Sensor, EPS)為我們提供了高能粒子通量數(shù)據(jù)來研究高能離子注入事件.EPS的觀測方向是垂直于GOES的自轉(zhuǎn)軸,且GOES的自轉(zhuǎn)軸與地球的自轉(zhuǎn)軸線大致對齊.EPS探測的高能質(zhì)子通量數(shù)據(jù)為5 min平均值.同時結合地磁指數(shù)AE, AL進行了事例分析.

    3 典型事例分析

    我們根據(jù)GOES衛(wèi)星的質(zhì)子通量數(shù)據(jù),同時利用地磁指數(shù)AE, AL,確定磁尾亞暴期間等離子體片離子注入內(nèi)磁層事件,再利用MCMU和MSTK兩個地磁臺站的FGM的高分辨率(8 Hz)數(shù)據(jù),研究兩個臺站同時觀測到的IPDP事件與亞暴離子注入事件的關系.

    圖1a顯示的是從2010年8月23日到24日時間段GOES 11衛(wèi)星探測的質(zhì)子通量的變化,圖1b顯示的是這兩天的AE, AL指數(shù)的變化.如圖1所示,在2010年8月23日23∶20 UT(14∶26 MLT),質(zhì)子通量有一個快速的增長.同時,AE指數(shù)在之前也一直沒有太大變化,直到2010年8月23日23∶35 UT突然有一個大的增長.為了更好地顯示GOES 11衛(wèi)星測量的質(zhì)子通量的變化與AE, AL指數(shù)變化之間的對應關系,如圖2顯示的是圖1從2010年8月23日23∶00 UT到2010年8月24日04∶00 UT時間段,即圖1中兩條黑色豎線標注的時間段內(nèi)的質(zhì)子通量與AE指數(shù)的變化.由圖2可知在23∶35 UT時AE指數(shù)有一個快速增長.通過研究AE, AL指數(shù)的變化我們可以確定在2010年8月23日23∶35 UT有亞暴發(fā)生,質(zhì)子通量從23∶20 UT開始也有一個大的增長,所以可以確定此時有亞暴離子注入事件產(chǎn)生.

    圖1 2010年8月23日到2010年8月24日質(zhì)子通量(a)和AE,AL指數(shù)(b)的變化 (圖中左邊的黑色豎線是亞暴起始,同時質(zhì)子通量有增長)Fig.1 (a) and (b) the variation of proton flux and AE,AL indices between 23th August 2010 and 24th August 2010, and the left side of black marked line represent the onset of substorm, the same time proton flux has an increase

    地磁分量Bx的實際測量值減去它的滑動平均值之后(滑動窗的長度為128 s)可得出擾動量(ΔBx),圖3是根據(jù)擾動量(ΔBx)經(jīng)過計算得到的功率頻譜圖和原始曲線圖.由圖3可知,MCMU和MSTK兩個臺站同時觀測到的一個Pc1-2事件.該事件發(fā)生于2010年8月24日00∶35 UT (16∶15 MLT),事件持續(xù)時間段為00∶35—01∶10 UT.由于兩個臺站觀測到的Pc1-2事件都明顯存在頻率隨時間的增長,即在頻譜形態(tài)上有一定的斜率,因此我們確定這個Pc1-2事件為一個同時觀測的IPDP事件.MCMU和MSTK兩個臺站同時觀測的IPDP事件在形態(tài)上有些差異,其頻率漂移的斜率在MSTK觀測大于在MCMU觀測,這可能是由于波動的源區(qū)移動及源區(qū)與地磁觀測臺站的相對位置造成的.圖3b功率譜圖中MSTK地磁臺站在2∶00—3∶00 UT也觀測到Pc1-2事件,這也是IPDP事件,但在圖3a中MCMU卻沒有觀測到.這可能是由于IPDP波動產(chǎn)生的源區(qū)距離MCMU地磁臺站太遠,IPDP波動在電離層波導傳輸過程中衰減,導致未被MCMU臺站觀測到.

    圖1—3的研究顯示,在IPDP事件發(fā)生之前,AE指數(shù)快速增長,且GOES 11衛(wèi)星也觀測到質(zhì)子通量快速增長,說明這個IPDP事件是由磁尾亞暴期間等離子體片離子注入內(nèi)磁層后西向漂移,激發(fā)EMIC波產(chǎn)生的.

    圖2 GOES 11的質(zhì)子通量數(shù)據(jù)(a)和AE,AL指數(shù)(b)在2010年8月23日 23∶00 UT 到2010年8月24日04∶00 UT的變化Fig.2 The variation of proton flux from GOES 11 and AE, AL indices between 23∶00 UT 23th August 2010 and 04∶00 UT 24th August 2010

    圖3 MCMU(A)和MSTK(B)在2010年8月24日的同時觀測到的IPDP事件Fig.3 (A) and (B) represent the IPDP wave event simultaneously observed by MCMU and MSTK in 24th August 2010

    4 統(tǒng)計結果

    我們挑選了MCMU和MSTK兩個臺站從2005年4月到2014年5月時間段內(nèi)同時觀測到的Pc1-2事件.當兩個臺站同時觀測的功率譜中存在大于0.1 Hz的頻帶信號并且振幅明顯大于背景平均振幅值,頻率有一個明顯的增長時,就判斷為一個有效的IPDP觀測事件,并且記錄下每個IPDP事件持續(xù)的時間段.將MLT劃分為8個時間段,分別為:0—3 MLT, 3—6 MLT, 6—9 MLT, 9—12 MLT,12—15 MLT, 15—18 MLT, 18—21 MLT, 21—24 MLT. 分別計算IPDP事件落在這8個MLT時間段的發(fā)生率.統(tǒng)計季節(jié)發(fā)生率時,季節(jié)的劃分如表1.

    表1 季節(jié)的劃分Table 1 The definition of the season

    我們所用的FGM數(shù)據(jù)是兩個地磁臺站從2005年4月到2014年5月這10年間的同時觀測數(shù)據(jù),通過統(tǒng)計分析獲得IPDP事件發(fā)生率關于MLT和季節(jié)分布的統(tǒng)計特征.我們總共挑選了128個同時觀測的IPDP事件.圖4a說明兩個臺站的同時觀測到的IPDP事件可以發(fā)生在所有的MLT,但是最大發(fā)生率在15—18 MLT時間段.在不同季節(jié),IPDP事件的發(fā)生率也有所不同.如圖4b所示最小的季節(jié)發(fā)生率在冬季,為13.28%,春夏秋都比較大,其中春季發(fā)生率最大,為32.81%.

    5 討論與結論

    EMIC波產(chǎn)生于磁赤道區(qū)域,由高能離子溫度的各向異性導致的離子回旋不穩(wěn)定性所激發(fā).圖4a顯示IPDP事件可以發(fā)生在所有的磁地方時,但是在15—18 MLT時間段發(fā)生率最大.在不同的季節(jié),IPDP事件的發(fā)生率也有所不同,如圖4b所示,春夏秋都比較大,其中春季最大,達到32.81%,但在冬季發(fā)生率最小,僅為13.28%.圖4a顯示IPDP事件的日發(fā)生率在15—18 MLT時間段呈現(xiàn)峰值,可能是因為離子回旋不穩(wěn)定區(qū)域的西向漂移質(zhì)子遇到高密度的等離子體層plume區(qū)或等離子體層頂都有利于產(chǎn)生EMIC波,等離子體層plume多形成于黃昏側,并向日側延伸(Yuan et al., 2010,Yuan et al., 2012b),所以IPDP事件多出現(xiàn)在黃昏側(Lee and Kwok, 1984, Hayakawa et al., 1992).而且MCMU和MSTK的L值分別為5.35和4.22,是經(jīng)常觀測到等離子體層plume所在的L值范圍,所以在15—18 MLT觀測到的IPDP事件較多,發(fā)生率高.根據(jù)上述圖4b可看出IPDP的季節(jié)分布特征,一是在冬季發(fā)生率最低,僅為13.28%,二是春夏秋發(fā)生率都比較大,其中春季發(fā)生率最高,達到32.81%.這可能與EMIC波的傳播有關系.EMIC波沿地球磁力線從磁層傳播到電離層E層并且向四周傳播.EMIC波在電離層波導中傳播時,電離層的電導率起著非常重要的作用.當缺少光照時,電離程度低,電導率小,不利于Pc1-2傳播 (Kim et al., 2011; 邱奕婷等,2014).CARISMA地磁臺鏈位于北美大陸,冬季缺少光照,因此IPDP波事件在冬季較少,發(fā)生率低.而在夏季由于在北美大陸光照強,電離程度高,電離層電導率大,有利于IPDP波的傳播,所以夏季IPDP事件多,發(fā)生率高.但從圖4b可以看出IPDP事件在春秋季發(fā)生率也比冬季高,而且在春季的發(fā)生率最高,秋季的發(fā)生率分別低于春季和夏季.這可能是因為在春分、秋分時地磁活動增強所致(Saito, 1969).亞暴的強度和地磁活動有一個顯著的年變化,二分點的最大值在春季和秋季之間變化(Mursula et al., 2011).圖5顯示的是從2005年4月到2014年5月時間段內(nèi)春夏秋冬各季節(jié)內(nèi)Kp指數(shù)大于3的數(shù)據(jù)點數(shù)的歸一化出現(xiàn)率.由圖5可知,Kp指數(shù)關于季節(jié)的歸一化出現(xiàn)率是春季高于夏季,夏季高于秋季,秋季高于冬季.因此春季的亞暴強度和地磁活動比秋季強,磁尾等離子體片離子注入事件多,利于IPDP波動的產(chǎn)生,所以春季的發(fā)生率高于秋季.同時可能因為春秋季的電離層電導率一般明顯小于夏季,不利于波動的傳播,且秋季的亞暴強度和地磁活動比春季弱,磁尾等離子體片離子注入事件沒有春季多,所以秋季的發(fā)生率低于夏季.綜上所述說明IPDP事件關于季節(jié)分布的發(fā)生率同時受到電離層電導率和亞暴發(fā)生率的影響.

    圖4 MCMU和MSTK 兩個臺站同時觀測的IPDP事件關于MLT(a)和季節(jié)(b)的分布Fig.4 (a) and (b) represent the occurrence rate distribution of IPDP event simultaneously observed by MCMU and MSTK in MLT and season respectively

    圖5 2005年4月到2014年5月時間段內(nèi)Kp指數(shù) 大于3 的數(shù)據(jù)點數(shù)關于季節(jié)的歸一化出現(xiàn)率Fig.5 From April 2005 until May 2014, normalized rate distribution of Kp>3 in season

    基于GOES衛(wèi)星、加拿大CARISMA地磁臺站的聯(lián)合觀測數(shù)據(jù),本文研究了IPDP事件的產(chǎn)生與亞暴磁尾注入的關系.同時,利用CARISMA地磁臺鏈中的MCMU和MSTK兩個臺站,從2005年4月到2014年5月期間的觀測數(shù)據(jù),統(tǒng)計分析了亞暴期間的IPDP事件,研究了IPDP事件的出現(xiàn)率關于季節(jié)和磁地方時的分布.統(tǒng)計結果表明:IPDP事件關于季節(jié)分布的發(fā)生率,冬季最小,春季最大.IPDP事件關于季節(jié)分布的發(fā)生率受到電離層電導率及亞暴發(fā)生率的影響.兩個臺站同時觀測到的IPDP事件最大出現(xiàn)率出現(xiàn)在15—18 MLT,結果表明IPDP事件主要由亞暴期間產(chǎn)生的高能離子注入內(nèi)磁層,西向漂移遇到等離子體層羽狀結構區(qū)的高密度等離子體所激發(fā).

    致謝 本文所用GOES 11衛(wèi)星數(shù)據(jù)來源于NASA Space Physics Data Facility(http:∥spdf.gsfc.nasa.gov/)網(wǎng)站下載的, 地磁指數(shù)數(shù)據(jù)來源于Geomagnetic Data service (http:∥wdc.kugi.kyoto-u.ac.jp/wdc/Sec3.html).本文所用FGM的數(shù)據(jù)來源于CARISMA.

    Cornwall J M. 1965. Cyclotron instabilities and electromagnetic emission in the ultra low frequency and very low frequency ranges.J.Geophys.Res., 70(1): 61-69, doi: 10.1029/JZ070i001p00061.Fok M C, Ebihara Y, Moore T E, et al. 2005. Geospace storm processes coupling the ring current, radiation belt and plasmasphere. ∥ Burch J, Schulz M, Spence H eds. Inner Magnetosphere Interactions: New Perspectives from Imaging. Geophys. Monogr. Ser., AGU, Washington, D. C., 159, 207-220, doi: 10.1029/159GM16.Fraser B J, Singer H J, Adrian M L, et al. 2005. The relationship between plasma density structure and EMIC waves at geosynchronous orbit. ∥Burch J, Schulz M, Spence H eds. Inner Magnetosphere Interactions: New Perspectives from Imaging, Geophys. Monogr. Ser., AGU, Washington, D. C.,159: 55-70, doi: 10.1029/159GM04.Fraser B J, Morley S K, Grew R S. 2012. Classification of Pc1-2 electromagnetic ion cyclotron waves at geosynchronous orbit. ∥Summers D, Mann I R, Baker D N, et al, eds. Dynamics of the Earth′s Radiation Belts and inner Magnetosphere, AGU, Washington, D.C., doi: 10.1029/2012GM001353

    Gendrin R. 1971. Substorm aspect of magnetic pulsations.SpaceSci.Rev., 11(1): 54-130.

    Gul′elmi A V. 1974. Diagnostics of the magnetosphere and interplanetary medium by means of pulsations.SpaceSci.Rev., 16(3): 331-345.

    Halford A J, Fraser B J, Morley S K. 2010. EMIC wave activity during geomagnetic storm and nonstorm periods: CRRES results.J.Geophys.Res., 115: A12248, doi: 10.1029/2010JA015716.Hayakawa M, Shimakura S, Kobayshi T, et al. 1992. A study of polarization of irregular pulsations of diminishing period and their generation mechanism.Planet.SpaceSci., 40(8): 1081-1091.Jordanova V K, Farrugia C J, Thorne R M, et al. 2001. Modeling ring current proton precipitation by electromagnetic ion cyclotron waves during the May 14-16, 1997, storm.J.Geophys.Res., 106(A1): 7-22, doi: 10.1029/2000JA002008.

    Kerttula R, Mursula K, Pikkarainen T, et al. 2001. Storm-time Pc1 activity at high and middle latitudes.J.Geophys.Res., 106(A4): 6213-6227, doi: 10.1029/2000JA900125.

    Kim H, Lessard M R, Engebretson M J, et al. 2011. Statistical study of Pc1-2 wave propagation characteristics in the high-latitude ionospheric waveguide.J.Geophys.Res., 116: A07227, doi: 10.1029/2010JA016355.

    Lee L C, Kwok Y C. 1984. A mechanism for the IPDP pulsations.J.Geophys.Res., 89(A2): 877-882, doi: 10.1029/JA089iA02p00877.Liu Y H, Fraser B J, Menk F W. 2012. Pc2 EMIC waves generated high off the equator in the dayside outer magnetosphere.Geophys.Res.Lett., 39: L17102, doi: 10.1029/2012GL053082.

    Loto′aniu T M, Fraser B J, Waters C L. 2005. Propagation of electromagnetic ion cyclotron wave energy in the magnetosphere.J.Geophys.Res., 110: A07214, doi: 10.1029/2004JA010816.Mursula K, Tanskanen E, Love J J. 2011. Spring-fall asymmetry of substorm strength, geomagnetic activity and solar wind: Implications for semiannual variation and solar hemispheric asymmetry.Geophys.Res.Lett., 38, L06104, doi: 10.1029/2011GL046751.Qiu Y T, Liu Y H, Zhao H F. 2014. Statistical characteristics of Pc1-2 waves near the cusp latitudes.AdvancesinPolarScience(in Chinese), 26(3): 324-330.

    Saito T. 1969. Geomagnetic pulsations. Space Sci. Rev., 10: 319-412, doi: 10.1007/BF00203620.

    Su Z P, Xiao F L, Zheng H N, et al. 2010. STEERB: A three dimensional code for storm-time evolution of electron radiation belt.J.Geophys.Res., 115: A09208, doi: 10.1029/2009JA015210.Su Z P, Xiao F L, Zheng H N, et al. 2011. CRRES observation and STEERB simulation of the 9 October 1990 electron radiation belt dropout event.Geophys.Res.Lett., 38: L06106, doi: 10.1029/2011GL046873.

    Wang Z Z, Yuan Z G, Li M, et al. 2014. Statistical characteristics of EMIC wave driven relativistic electron precipitation with observations of POES satellites: Revisit.J.Geophys.Res.SpacePhysics, 119(7): 5509-5519, doi: 10.1002/2014JA020082.Xiao F L, Chen L X, He Y H, et al. 2011. Modeling for precipitation loss of ring current protons by electromagnetic ion cyclotron waves.J.Atmos.Sol.Terr.Phys., 73(1): 106-111, doi: 10.1016/j.jastp.2010.01.007.

    Yuan Z G, Deng X H, Lin X, et al. 2010. Link between EMIC waves in a plasmaspheric plume and a detached sub-auroral proton arc with observations of Cluster and IMAGE satellites.Geophys.Res.Lett., 37: L07108, doi: 10.1029/2010GL042711.

    Yuan Z G, Xiong Y, Pang Y, et al. 2012a. Wave-particle interaction in a plasmaspheric plume observed by a Cluster satellite.J.Geophys.Res., 117: A03205, doi: 10.1029/2011JA017152.

    Yuan Z G, Xiong Y, Wang D D, et al. 2012b. Characteristics of precipitating energetic ions/electrons associated with the wave-particle interaction in the plasmaspheric plume.J.Geophys.Res., 117: A08324, doi: 10.1029/2012JA017783.

    Yuan Z G, Li M, Xiong Y, et al. 2013. Simultaneous observations of precipitating radiation belt electrons and ring current ions associated with the plasmaspheric plume.J.Geophys.Res.SpacePhysics, 118(7): 4391-4399, doi: 10.1002/jgra.50432.

    Yuan Z G, Xiong Y, Huang S Y, et al. 2014. Cold electron heating by EMIC waves in the plasmaspheric plume with observations of the Cluster satellite.Geophys.Res.Lett., 41(6): 1830-1837, doi: 10.1002/2014GL059241.

    Zhou Q H, Xiao F L, Yang C, et al. 2013. Observation and modeling of magnetospheric cold electron heating by electromagnetic ion cyclotron waves.J.Geophys.Res.SpacePhysics, 118(11): 6907-6914, doi: 10.1002/2013JA019263.

    附中文參考文獻

    邱奕婷, 劉勇華, 趙浩峰. 2014. 極隙區(qū)緯度Pc1-2波的統(tǒng)計分布特征. 極地研究, 26(3): 324-330.

    (本文編輯 汪海英)

    Statistical study of IPDP events observed by the CARISMA network of magnetometers

    ZHOU Ming-Xia, YUAN Zhi-Gang*, LI Hai-Meng, WANG De-Dong,YU Xiong-Dong, WANG Zhen-Zhen, QIAO Zheng

    SchoolofElectronicInformation,WuhanUniversity,Wuhan430072,China

    During substorms energetic ions injected from the magnetotail plasma sheet can excite electromagnetic ion cyclotron (EMIC) waves in the inner magnetosphere. As a signature of EMIC waves, the interval of pulsations of diminishing periods (IPDP) is observed by ground-based magnetometers. The generating mechanism of IPDP waves is closely associated with substorm injected energetic ions during substorms. In addition, the occurrence rate of IPDP waves in the inner magnetosphere is very important to study the dynamics of the magnetosphere.Firstly, with observations of magnetometers at MCMU and MSTK stations of CARISMA network in Canada and the GOES 11 satellite for the typical IPDP events, the relationship between IPDP and the energetic ions injected from the magnetotail plasma sheet during substorms was revealed. Then based on the observations of magnetometers at MCMU and MSTK stations of CARISMA network between April 2005 and May 2014, we statistically analyzed the occurrence rate distribution of IPDP events in seasons and magnetic local time (MLT).The results for the typical case show that during substorms, energetic ions injected from the magnetotail plasma sheet can excite IPDP waves in the inner magnetosphere. For the statistics, a total of 128 IPDP wave events have been simultaneously observed by both MCMU and MSTK stations. The maximum occurrence rate of the events is between 15 MLT and 18 MLT, corresponding to the region of plasmasphere plumes. The IPDP events have the least occurrence rate in winter (13.28%), the highest occurrence rate is in spring (32.81%). The least occurrence rate in winter should be related to the ionospheric conductivity while the highest in spring is associated with the spring-autumn asymmetry of substorm and geomagnetic activity.The case study and statistics show that IPDP events are mainly generated by energetic protons, which can be injected into the inner magnetosphere during substorms and drift westward and encounter the plasmaspheric plume. In addition, the occurrence rate of IPDP events is affected by the ionospheric conductivity and the occurrence rate of substorms.

    EMIC; IPDP; Substorm injection; Occurrence rate

    10.6038/cjg20150702.

    國家自然科學基金(41174140,41221003,41374168), 高等學校博士點基金(20110141110043), 新世紀優(yōu)秀人才項目(NCET-13-0446)共同資助.

    周銘霞,女,1989年生,碩士研究生,主要從事電離層磁層物理方面的研究.E-mail: mingxiazhou@whu.edu.cn

    *通訊作者 袁志剛, 男,1974年生,教授,主要從事電離層磁層耦合方面的研究.E-mail: y_zgang@vip.163.com

    10.6038/cjg20150702

    P352

    2014-10-13,2015-06-19收修定稿

    周銘霞,袁志剛,李海夢等. 2015. 基于CARISMA地磁臺鏈觀測的IPDP事件的統(tǒng)計分析.地球物理學報,58(7):2235-2241,

    Zhou M X, Yuan Z G, Li H M, et al. 2015. Statistical study of IPDP events observed by the CARISMA network of magnetometers.ChineseJ.Geophys. (in Chinese),58(7):2235-2241,doi:10.6038/cjg20150702.

    猜你喜歡
    質(zhì)子臺站等離子體
    中國科學院野外臺站檔案工作回顧
    氣象基層臺站建設
    西藏科技(2021年12期)2022-01-17 08:46:38
    連續(xù)磁活動對等離子體層演化的影響
    基于低溫等離子體修飾的PET/PVC浮選分離
    質(zhì)子束放療在腫瘤中的研究新進展
    等離子體種子處理技術介紹
    淺談質(zhì)子守恒
    基層臺站綜合觀測業(yè)務管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    “質(zhì)子”號一箭發(fā)雙星
    太空探索(2014年6期)2014-07-10 13:06:11
    MDOS平臺臺站級使用方法及技巧
    国产av码专区亚洲av| 亚洲精品视频女| 考比视频在线观看| 韩国高清视频一区二区三区| 激情五月婷婷亚洲| 菩萨蛮人人尽说江南好唐韦庄| 男女下面插进去视频免费观看 | 精品久久久久久电影网| 亚洲欧美一区二区三区黑人 | a级片在线免费高清观看视频| 天天躁夜夜躁狠狠躁躁| 亚洲四区av| 亚洲,一卡二卡三卡| 777米奇影视久久| 日韩中字成人| 全区人妻精品视频| 亚洲在久久综合| 一级爰片在线观看| 欧美 日韩 精品 国产| 黄色视频在线播放观看不卡| 卡戴珊不雅视频在线播放| 欧美另类一区| 少妇人妻久久综合中文| 超色免费av| 欧美日韩一区二区视频在线观看视频在线| 国产伦理片在线播放av一区| 国产精品久久久av美女十八| 如何舔出高潮| 最后的刺客免费高清国语| 少妇人妻久久综合中文| 九草在线视频观看| 精品一品国产午夜福利视频| 韩国精品一区二区三区 | 男女边摸边吃奶| 全区人妻精品视频| 久久久久国产网址| 久久这里有精品视频免费| av黄色大香蕉| 欧美成人精品欧美一级黄| 最近2019中文字幕mv第一页| 丰满少妇做爰视频| 午夜免费鲁丝| 国产亚洲午夜精品一区二区久久| www.色视频.com| 欧美精品一区二区免费开放| 国产又色又爽无遮挡免| 乱码一卡2卡4卡精品| 18禁动态无遮挡网站| 51国产日韩欧美| 亚洲欧美清纯卡通| 亚洲欧美日韩卡通动漫| 亚洲精品国产色婷婷电影| 少妇的丰满在线观看| 视频区图区小说| 香蕉国产在线看| 91午夜精品亚洲一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 久久久久国产网址| 久久亚洲国产成人精品v| 国产色婷婷99| 亚洲少妇的诱惑av| 欧美日韩亚洲高清精品| 久久人人爽人人片av| 中文乱码字字幕精品一区二区三区| 亚洲高清免费不卡视频| 久久精品久久精品一区二区三区| 久久久精品94久久精品| 亚洲综合色网址| 亚洲一级一片aⅴ在线观看| 免费在线观看黄色视频的| 蜜臀久久99精品久久宅男| 日韩av不卡免费在线播放| 69精品国产乱码久久久| 免费播放大片免费观看视频在线观看| 欧美日韩精品成人综合77777| 中文天堂在线官网| 永久网站在线| 国产av精品麻豆| 日韩熟女老妇一区二区性免费视频| 免费在线观看完整版高清| 久热这里只有精品99| 在线观看一区二区三区激情| 久久精品国产a三级三级三级| 亚洲精品国产av蜜桃| 国产无遮挡羞羞视频在线观看| 99热6这里只有精品| 美女国产高潮福利片在线看| 大香蕉97超碰在线| 丰满少妇做爰视频| 18禁在线无遮挡免费观看视频| 国精品久久久久久国模美| 人人妻人人添人人爽欧美一区卜| 亚洲欧洲精品一区二区精品久久久 | 久久99蜜桃精品久久| 日韩av在线免费看完整版不卡| 成人毛片60女人毛片免费| 午夜福利在线观看免费完整高清在| 欧美最新免费一区二区三区| av有码第一页| 久久鲁丝午夜福利片| 亚洲成av片中文字幕在线观看 | 亚洲成人一二三区av| 国产精品国产三级国产专区5o| 久久精品夜色国产| 九色亚洲精品在线播放| 97人妻天天添夜夜摸| 人妻少妇偷人精品九色| 亚洲国产色片| 丝袜脚勾引网站| av女优亚洲男人天堂| 国产老妇伦熟女老妇高清| 成年人午夜在线观看视频| 精品一区二区三区视频在线| 成年人午夜在线观看视频| 丝袜脚勾引网站| 最近2019中文字幕mv第一页| 欧美xxxx性猛交bbbb| 97精品久久久久久久久久精品| 97精品久久久久久久久久精品| 亚洲精品第二区| 在线免费观看不下载黄p国产| 最近的中文字幕免费完整| 97精品久久久久久久久久精品| 97精品久久久久久久久久精品| 亚洲精品第二区| 久久99热6这里只有精品| 欧美3d第一页| 亚洲国产成人一精品久久久| 成年av动漫网址| 新久久久久国产一级毛片| 91成人精品电影| 夜夜爽夜夜爽视频| 精品一区二区三卡| 一区二区日韩欧美中文字幕 | av网站免费在线观看视频| 国产黄色免费在线视频| 亚洲精品乱码久久久久久按摩| 午夜日本视频在线| 日产精品乱码卡一卡2卡三| 香蕉国产在线看| 日韩伦理黄色片| 日韩制服骚丝袜av| 18禁裸乳无遮挡动漫免费视频| 亚洲av福利一区| 多毛熟女@视频| 国产视频首页在线观看| 天堂8中文在线网| 啦啦啦视频在线资源免费观看| 美女内射精品一级片tv| 国语对白做爰xxxⅹ性视频网站| 2021少妇久久久久久久久久久| 国产精品嫩草影院av在线观看| 精品一区二区三区视频在线| 少妇人妻精品综合一区二区| 国产精品一二三区在线看| 女的被弄到高潮叫床怎么办| 免费黄网站久久成人精品| 亚洲精品一区蜜桃| 国产精品久久久av美女十八| 精品一品国产午夜福利视频| 欧美老熟妇乱子伦牲交| 亚洲在久久综合| 亚洲经典国产精华液单| 免费高清在线观看日韩| 精品熟女少妇av免费看| 国产免费福利视频在线观看| 天天躁夜夜躁狠狠久久av| 超色免费av| 日韩精品有码人妻一区| 最近最新中文字幕免费大全7| 777米奇影视久久| av天堂久久9| 日日爽夜夜爽网站| 婷婷色综合www| 欧美xxxx性猛交bbbb| 久热久热在线精品观看| 亚洲精品一区蜜桃| 亚洲第一av免费看| 国产成人午夜福利电影在线观看| 精品人妻熟女毛片av久久网站| videosex国产| 满18在线观看网站| 丝袜人妻中文字幕| 少妇被粗大猛烈的视频| 七月丁香在线播放| 欧美国产精品一级二级三级| 国产精品女同一区二区软件| 久久 成人 亚洲| 国产精品久久久久成人av| 亚洲五月色婷婷综合| 欧美 日韩 精品 国产| tube8黄色片| 国产av一区二区精品久久| 捣出白浆h1v1| 免费观看av网站的网址| 一级片免费观看大全| 中文字幕人妻熟女乱码| 亚洲av.av天堂| 欧美少妇被猛烈插入视频| 久久久久久久精品精品| 亚洲综合精品二区| 久久人人爽人人片av| 男人舔女人的私密视频| 亚洲国产看品久久| 精品人妻熟女毛片av久久网站| 国产精品久久久久久久电影| 视频在线观看一区二区三区| 欧美日韩精品成人综合77777| 国产欧美另类精品又又久久亚洲欧美| 国产视频首页在线观看| 大片免费播放器 马上看| 水蜜桃什么品种好| 三上悠亚av全集在线观看| 亚洲国产毛片av蜜桃av| 国产毛片在线视频| 国产熟女欧美一区二区| 视频在线观看一区二区三区| 国产精品秋霞免费鲁丝片| 色哟哟·www| 亚洲精品av麻豆狂野| 美女国产高潮福利片在线看| 晚上一个人看的免费电影| 国产日韩欧美在线精品| 久久久久久久久久成人| 中文字幕人妻丝袜制服| 日本与韩国留学比较| 成人国语在线视频| 人人澡人人妻人| 热99国产精品久久久久久7| 欧美3d第一页| 亚洲第一av免费看| 九草在线视频观看| 久久久久久久久久久久大奶| 美女主播在线视频| 国产男人的电影天堂91| 亚洲 欧美一区二区三区| 成年人午夜在线观看视频| 久久久久久人人人人人| 久久狼人影院| 国产精品.久久久| 日韩人妻精品一区2区三区| 搡女人真爽免费视频火全软件| 一级毛片电影观看| videossex国产| 性色av一级| 亚洲国产毛片av蜜桃av| 在线看a的网站| √禁漫天堂资源中文www| 精品少妇内射三级| 狂野欧美激情性bbbbbb| av在线app专区| 在线精品无人区一区二区三| 亚洲国产精品一区二区三区在线| 国产熟女午夜一区二区三区| 国产亚洲午夜精品一区二区久久| 亚洲高清免费不卡视频| 人妻系列 视频| 天天影视国产精品| 大香蕉久久网| 亚洲一区二区三区欧美精品| 一本—道久久a久久精品蜜桃钙片| 欧美日韩视频高清一区二区三区二| 五月天丁香电影| 亚洲图色成人| 欧美老熟妇乱子伦牲交| 亚洲国产精品一区二区三区在线| 亚洲欧洲国产日韩| 久久影院123| 在现免费观看毛片| av线在线观看网站| 亚洲精品美女久久久久99蜜臀 | freevideosex欧美| 99热网站在线观看| 亚洲色图综合在线观看| 韩国av在线不卡| 精品国产一区二区三区四区第35| 亚洲经典国产精华液单| 两性夫妻黄色片 | 国产老妇伦熟女老妇高清| 久久精品国产a三级三级三级| 高清视频免费观看一区二区| 伊人亚洲综合成人网| 国产精品国产三级国产av玫瑰| 黄色毛片三级朝国网站| 啦啦啦啦在线视频资源| 考比视频在线观看| 大话2 男鬼变身卡| 视频区图区小说| 久久久久久久精品精品| 亚洲五月色婷婷综合| 男人添女人高潮全过程视频| 美女脱内裤让男人舔精品视频| 欧美bdsm另类| 久久综合国产亚洲精品| 人妻少妇偷人精品九色| 久久久久久久亚洲中文字幕| 国产男女超爽视频在线观看| av福利片在线| 免费大片18禁| 99热全是精品| 考比视频在线观看| 天天操日日干夜夜撸| 成人国产av品久久久| 国产精品久久久久久av不卡| 欧美亚洲日本最大视频资源| 免费在线观看完整版高清| 亚洲精品自拍成人| 少妇的逼好多水| 免费日韩欧美在线观看| av黄色大香蕉| 亚洲av电影在线进入| 久久久久精品久久久久真实原创| 视频中文字幕在线观看| 青春草国产在线视频| 18+在线观看网站| 亚洲激情五月婷婷啪啪| 有码 亚洲区| 在线 av 中文字幕| 国产一区有黄有色的免费视频| 成人亚洲欧美一区二区av| 亚洲美女视频黄频| 久久韩国三级中文字幕| 久久久久久人人人人人| 不卡视频在线观看欧美| 香蕉国产在线看| 在线观看美女被高潮喷水网站| 久久99热6这里只有精品| kizo精华| 欧美亚洲 丝袜 人妻 在线| 国产视频首页在线观看| 狂野欧美激情性bbbbbb| 一区二区日韩欧美中文字幕 | 免费久久久久久久精品成人欧美视频 | av国产久精品久网站免费入址| 自线自在国产av| 色94色欧美一区二区| 亚洲欧洲国产日韩| av卡一久久| 日韩制服丝袜自拍偷拍| 久久女婷五月综合色啪小说| 精品一区二区三区四区五区乱码 | 97在线视频观看| 黄色视频在线播放观看不卡| 日韩欧美一区视频在线观看| 成人毛片a级毛片在线播放| 有码 亚洲区| 最近2019中文字幕mv第一页| 赤兔流量卡办理| 汤姆久久久久久久影院中文字幕| 国产亚洲最大av| 三级国产精品片| 欧美人与善性xxx| 免费观看a级毛片全部| 亚洲欧美日韩卡通动漫| 自线自在国产av| 久久久久久久久久人人人人人人| 女人精品久久久久毛片| 亚洲欧美清纯卡通| 国产xxxxx性猛交| 母亲3免费完整高清在线观看 | 国产黄色免费在线视频| 久久午夜综合久久蜜桃| 777米奇影视久久| 国产片特级美女逼逼视频| 久久精品熟女亚洲av麻豆精品| 伊人亚洲综合成人网| 欧美精品av麻豆av| 欧美丝袜亚洲另类| 日韩,欧美,国产一区二区三区| 久久久久精品久久久久真实原创| 99久久综合免费| 国产精品麻豆人妻色哟哟久久| 国产成人精品久久久久久| 有码 亚洲区| 99热全是精品| 免费看光身美女| 国产麻豆69| 2021少妇久久久久久久久久久| 成年动漫av网址| 欧美日韩成人在线一区二区| 看免费av毛片| 乱人伦中国视频| 亚洲欧美一区二区三区黑人 | 免费av中文字幕在线| 国产亚洲一区二区精品| 欧美日本中文国产一区发布| 亚洲美女视频黄频| 国内精品宾馆在线| 考比视频在线观看| 亚洲色图 男人天堂 中文字幕 | 国产爽快片一区二区三区| 精品国产露脸久久av麻豆| 午夜激情久久久久久久| 人妻系列 视频| 女人久久www免费人成看片| 国产成人免费观看mmmm| 国产精品久久久久久精品电影小说| 亚洲精华国产精华液的使用体验| 日本色播在线视频| 9191精品国产免费久久| 亚洲三级黄色毛片| 日本与韩国留学比较| 国产精品久久久久久av不卡| 一级爰片在线观看| 国产极品天堂在线| 国产1区2区3区精品| 国产日韩欧美在线精品| 最新的欧美精品一区二区| av在线播放精品| tube8黄色片| 免费高清在线观看视频在线观看| 啦啦啦视频在线资源免费观看| av在线老鸭窝| 男女边摸边吃奶| 高清不卡的av网站| 在线观看一区二区三区激情| 丝袜人妻中文字幕| 国产精品久久久久久av不卡| 欧美97在线视频| 黄片无遮挡物在线观看| 国产精品久久久久久av不卡| 男人添女人高潮全过程视频| 日本av免费视频播放| 国语对白做爰xxxⅹ性视频网站| 赤兔流量卡办理| 国产极品粉嫩免费观看在线| 99久国产av精品国产电影| 国产伦理片在线播放av一区| 美女国产视频在线观看| 久久久久国产精品人妻一区二区| 精品酒店卫生间| 九色成人免费人妻av| 亚洲国产欧美日韩在线播放| 午夜福利在线观看免费完整高清在| 成人18禁高潮啪啪吃奶动态图| 国产在线一区二区三区精| 亚洲精品乱码久久久久久按摩| 美女福利国产在线| 久久精品aⅴ一区二区三区四区 | 色94色欧美一区二区| 精品久久国产蜜桃| 精品99又大又爽又粗少妇毛片| 交换朋友夫妻互换小说| 亚洲,一卡二卡三卡| 男女无遮挡免费网站观看| av免费观看日本| 精品一区二区三区视频在线| 十八禁网站网址无遮挡| 欧美精品国产亚洲| 中文字幕亚洲精品专区| 男的添女的下面高潮视频| 久久 成人 亚洲| 成人午夜精彩视频在线观看| 一区二区av电影网| 亚洲精品第二区| 精品国产一区二区三区久久久樱花| 搡老乐熟女国产| 乱码一卡2卡4卡精品| 成人手机av| 五月天丁香电影| 伦理电影大哥的女人| 国产精品三级大全| 久久精品国产亚洲av天美| 久久久亚洲精品成人影院| 精品国产国语对白av| 亚洲,欧美精品.| xxx大片免费视频| 午夜精品国产一区二区电影| 精品福利永久在线观看| 大片免费播放器 马上看| 日日撸夜夜添| 亚洲精品第二区| 黑丝袜美女国产一区| 亚洲情色 制服丝袜| 日本欧美视频一区| 又大又黄又爽视频免费| 一级黄片播放器| 欧美xxⅹ黑人| 卡戴珊不雅视频在线播放| 亚洲精品一二三| 亚洲av福利一区| 九九爱精品视频在线观看| 蜜臀久久99精品久久宅男| www日本在线高清视频| 90打野战视频偷拍视频| 久久午夜综合久久蜜桃| 久久国产精品男人的天堂亚洲 | 国产精品久久久久久久久免| 伊人亚洲综合成人网| 咕卡用的链子| 亚洲精品美女久久av网站| av一本久久久久| 亚洲欧美色中文字幕在线| 视频区图区小说| 妹子高潮喷水视频| 精品人妻偷拍中文字幕| 日韩免费高清中文字幕av| 黄色配什么色好看| 青春草亚洲视频在线观看| 久久狼人影院| 一级毛片 在线播放| 亚洲精品国产av蜜桃| 大香蕉久久成人网| 久久久久久久久久成人| 2018国产大陆天天弄谢| 日日啪夜夜爽| 最黄视频免费看| 韩国高清视频一区二区三区| 男的添女的下面高潮视频| 亚洲精品国产av成人精品| 女性生殖器流出的白浆| 久久人人爽av亚洲精品天堂| 精品一品国产午夜福利视频| 女人精品久久久久毛片| 精品99又大又爽又粗少妇毛片| 久久99蜜桃精品久久| 日韩 亚洲 欧美在线| 永久网站在线| 最后的刺客免费高清国语| 国产成人精品福利久久| 人人澡人人妻人| 国产成人午夜福利电影在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产日韩欧美亚洲二区| 日韩不卡一区二区三区视频在线| 蜜桃国产av成人99| 久久99精品国语久久久| 国产黄色视频一区二区在线观看| 男女免费视频国产| 人妻人人澡人人爽人人| 国产免费现黄频在线看| 一区二区三区四区激情视频| 高清视频免费观看一区二区| 国产精品一区二区在线不卡| 日韩欧美一区视频在线观看| 啦啦啦在线观看免费高清www| 丝袜喷水一区| 国产成人91sexporn| 日韩制服丝袜自拍偷拍| 最黄视频免费看| 飞空精品影院首页| 亚洲天堂av无毛| 婷婷色av中文字幕| 成人手机av| 激情视频va一区二区三区| 亚洲精品久久久久久婷婷小说| 欧美xxxx性猛交bbbb| a级毛片在线看网站| 亚洲欧美中文字幕日韩二区| 日本黄大片高清| 欧美 亚洲 国产 日韩一| 日本免费在线观看一区| 美国免费a级毛片| 色网站视频免费| 午夜老司机福利剧场| 最近中文字幕2019免费版| 好男人视频免费观看在线| 亚洲av电影在线进入| 热re99久久国产66热| 国产精品一区二区在线观看99| 精品少妇内射三级| 波野结衣二区三区在线| 大片免费播放器 马上看| 亚洲经典国产精华液单| 国产成人a∨麻豆精品| 日韩熟女老妇一区二区性免费视频| 啦啦啦视频在线资源免费观看| 国产欧美亚洲国产| 国产1区2区3区精品| 国产免费现黄频在线看| 好男人视频免费观看在线| 亚洲人成77777在线视频| 午夜91福利影院| 一级毛片黄色毛片免费观看视频| 最近手机中文字幕大全| 亚洲第一区二区三区不卡| 日本欧美国产在线视频| 亚洲国产精品一区二区三区在线| 在线 av 中文字幕| 22中文网久久字幕| 一级黄片播放器| 黄片无遮挡物在线观看| 婷婷成人精品国产| 久久精品久久久久久噜噜老黄| 欧美成人午夜免费资源| 亚洲av.av天堂| 午夜影院在线不卡| 一边亲一边摸免费视频| 欧美变态另类bdsm刘玥| 精品第一国产精品| 精品人妻熟女毛片av久久网站| 国精品久久久久久国模美| 色94色欧美一区二区| 人妻人人澡人人爽人人| 欧美xxⅹ黑人| 国产成人aa在线观看| 2022亚洲国产成人精品| 亚洲成人一二三区av| 又大又黄又爽视频免费| 涩涩av久久男人的天堂| 亚洲精品色激情综合| 久久精品aⅴ一区二区三区四区 | 国产精品一国产av| 中文字幕免费在线视频6| 亚洲国产精品专区欧美| 国产日韩欧美在线精品| 国产精品三级大全| 免费在线观看完整版高清| 亚洲精品av麻豆狂野| 最近中文字幕2019免费版| 国产老妇伦熟女老妇高清| 婷婷色麻豆天堂久久| 国产精品久久久久久久久免| 亚洲国产av新网站| 国产精品麻豆人妻色哟哟久久|