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

    利用轉移熵研究引起磁暴擾動的太陽風參數(shù)重要性排序*

    2022-06-20 09:03:54于佳斌佟繼周方少峰胡曉彥
    空間科學學報 2022年3期
    關鍵詞:磁層磁暴太陽活動

    于佳斌 佟繼周 方少峰 胡曉彥

    1(中國科學院國家空間科學中心 北京 100190)

    2(中國科學院大學 北京 100049)

    0 引言

    太陽風–磁層相互作用形成了一個多尺度耦合的復雜系統(tǒng),太陽風是引起地磁暴的最主要驅動源,例如當行星際磁場南向時,行星際磁場可與地球磁場發(fā)生磁場重聯(lián),較易引發(fā)磁暴。因而,理解太陽風變化對地磁場的影響機制,是預報地磁暴的關鍵之一。

    有關太陽風與地磁擾動之間的相關性研究,新方法不斷出現(xiàn),這些方法主要基于太陽風和地磁參數(shù)的解析關系、相關系數(shù)或預報模型,線性回歸、統(tǒng)計相關等方法,已被證實是認識磁暴時地磁變化的有效方法,例如Liu 等[1]研究了太陽活動上升年的地磁Kp指數(shù)與行星際磁場南向分量Bz、太陽風速度vsw、太陽風溫度Tsw與太陽風數(shù)密度Dsw的相關性,相關系數(shù)分別為0.66、0.03、0.58 和–0.45,發(fā)現(xiàn)有一類磁暴與Bz地磁場重聯(lián)無關。Cane 等[2]統(tǒng)計結果表明行星際磁場南向分量Bz與 磁暴擾動指數(shù)最小值Dstmin的相關系數(shù)為0.74;Wu[3]統(tǒng)計分析了1995-1998 年Wind 數(shù)據(jù),證實vswBzmin與Dstmin的相關系數(shù)為0.87,略高于行星際磁場南向分量最小值Bzmin與磁暴擾動指數(shù)最小值Dstmin的相關系數(shù)0.81。Zhang 等[4]發(fā)現(xiàn)行星際電場E與Dst指數(shù)有很好的相關性,并且在中、強等地磁活動期間,存在顯著的突變特征曲線,相對于vsw,,vswB2和 能量耦合函數(shù)ε=,行星際電場E的突變特征曲線更易識別。Zhao 等[5]的研究表明在諸多參數(shù)中,就單一因素來說,行星際電場E對磁暴強度影響最大,行星際磁場南向分量Bz對磁暴強度影響次之;Khabarova等[6]和Ahmed 等[7]也得到了相似的結果。除此之外,Iyemori 等[8]首次成功地應用線性預測濾波方法,基于太陽風參數(shù)預測了地磁活動,Ji 等[9]和Rast?tter等[10]用經(jīng)驗和物理模型、Ahmed 等[7]用神經(jīng)網(wǎng)絡對磁暴進行預報研究。然而,地磁系統(tǒng)對太陽風變化的響應是非線性的,線性統(tǒng)計不足以表達其完備的相關性,同時在使用各種預報模型時,不同的輸入?yún)?shù)對磁暴的影響和重要性仍不明確,關于哪個太陽風參數(shù)對磁暴的發(fā)生更重要,還未有統(tǒng)一的認識。

    近年來,基于信息論的新方法[11]被提出并成功地應用于揭示地球磁層–電離層系統(tǒng)響應太陽風變化的復雜動力學特征。Michelis 等[12]以地磁活動指數(shù)AL和Sym-H分別定量化亞暴和磁暴過程,進行了雙變量轉移熵分析,結果表明亞暴和磁暴之間的信息流動方向取決于全球地磁活動水平。Wing 等[13]利用信息論進行太陽風–地球輻射帶系統(tǒng)研究,利用互信息、條件互信息和轉移熵對能量范圍為1.8~3.5 MeV的地球同步軌道高能電子通量的太陽風驅動因素進行了研究,得出了基于信息論分析的重要性排序。轉移熵可以提供強大的無模型統(tǒng)計量,在太陽風–地球輻射帶系統(tǒng)研究已經(jīng)得到了驗證,但目前還未在太陽風–磁層整體耦合系統(tǒng)中得到應用。本文利用轉移熵研究太陽風參數(shù)與地磁指數(shù)之間的非線性關系,應用轉移熵來衡量磁暴時的太陽風驅動因素并進行量化,對不同太陽風參數(shù)的重要性進行研究。

    1 數(shù)據(jù)與方法

    1.1 數(shù)據(jù)來源

    用于研究的太陽風及地磁活動指數(shù)的數(shù)據(jù)均來源于美國國家航空航天局(NASA)的OMNI*http://omniweb.gsfc.nasa.gov/公開數(shù)據(jù),具體包括分鐘分辨率數(shù)據(jù)(HRO)和小時分辨率數(shù)據(jù)(LRO),時間跨度從1996 年到2018 年。OMNI 整合了來自ACE、WIND、IMP8 等航天器的太陽風磁場和等離子體數(shù)據(jù)集,并根據(jù)實際的太陽風速度時移到地球的弓激波鼻尖處。此外,本研究還使用了世界地磁數(shù)據(jù)中心(WDC)提供的地磁Sym-H指數(shù)數(shù)據(jù)、比利時SILSO 中心提供的太陽黑子數(shù)據(jù)、和加拿大自然資源部的太陽10.7 cm 射電流量F10.7數(shù)據(jù)。

    本研究利用美國空間天氣研究中心(SWRC)提供的空間天氣數(shù)據(jù)庫(DONKI)**https://kauai.ccmc.gsfc.nasa.gov/DONKI/中提取的磁暴事件列表,共計93 個事例,主要分布在第24 太陽活動周期間的2010 年到2018 年。

    1.2 相關系數(shù)與轉移熵

    1.2.1 相關系數(shù)算法

    相關系數(shù)又稱皮爾遜積矩相關系數(shù),適用于度量兩個變量x與y之間的線性相關性,其值介于1 與–1之間,即

    直接推導可知

    時移相關系數(shù)[13]是皮爾遜積矩相關系數(shù)度量向兩個時間序列x(t)和y(t+τ)的拓展,即

    其中τ為時間延遲。當相關函數(shù)有多個峰或沒有明顯的不對稱時,這種分析的結果可能不是特別清楚。此外,相關系數(shù)是兩個變量之間線性關系的一個度量,不能用于描述非線性關系。

    1.2.2 轉移熵算法

    度量兩個變量x和y之間依賴關系,更通用的方法是通過考慮

    其中,P(x,y)為x和y的聯(lián)合概率。如果式(2)等號成立,則表示變量x與y相互獨立;相比式(1),除線性關系外,式(2)可以表達更多的非線性關系。

    信息論中互信息[14]表示一個隨機變量中包含的關于另一個隨機變量的信息量,代表兩個隨機變量的依賴程度。

    對于離散的概率分布,互信息的定義如下:

    從式(3)和(4)可以看出,互信息M(x,y)對式(2)進行了定量描述,相比于相關系數(shù),更好地表達了兩個變量x與y之間的依賴關系。進一步,如果x與y服從聯(lián)合高斯概率分布,互信息具有如下表達式[14]:

    互信息的取值范圍為[ 0,∞),經(jīng)過轉化Λ(x,y)的取值范圍位于 [ 0,1]。特別地,當Λ(x,y)=0 時,互信息M(x,y)=0,從其定義可知p(x,y)=p(x)p(y),即x與y相互獨立。

    Λ(x,y)是一個同時包含線性和非線性關系的度量,而互相關系數(shù)r(x,y)和Λ(x,y)之間的差異

    可用于揭示變量之間除線性關系之外更復雜的關聯(lián)關系。

    當已知z∈{z1,z2,...,zn}≡N3時,衡量x與y的互信息即條件互信息[13]有

    轉移熵是一個量化兩個系統(tǒng)相關程度或混沌程度的參數(shù),可用于處理有限長度的信號數(shù)據(jù)。轉移熵定義由Schreiber[15]于2000 年提出,同時考慮源序列和目標序列的數(shù)據(jù)長度,通用的轉移熵定義為(默認單位為nats)

    其 中,yp(t)=[y(t),y(t ??),...,y(t ?k?)],k+1為歷史時間序列的長度,?為最小的取樣間隔,τ為時間延遲。轉移熵可以被視為一種特殊形式的條件互信息[13],即

    轉移熵可以衡量系統(tǒng)中已知目標歷史序列yp的情況下,有多少信息從輸入序列X傳遞到目標序列Y的下一狀態(tài)。轉移熵與相關系數(shù)的不同在于方向性,即TX→Y TY →X。

    轉移熵算法容易受到有限大小樣本效應的影響,也容易受到數(shù)值效應的影響,因此建立一個閾值來界定轉移熵計算結果是否顯著是必要的。為了檢驗計算結果的統(tǒng)計學意義,將源序列x隨機打亂重新抽樣,這一步驟保證了置換時間序列具有與原始序列相同的均值、方差、自相關函數(shù),但破壞了非線性關系。這個過程的目的是消除兩個序列數(shù)據(jù)x與y之間的所有潛在關聯(lián)關系,因此計算得到的Tsur(X)→Y理論值為零。但是在有限時間序列中,由于有限樣本效應Tsur(X)→Y很少為零,故而可以得到轉移熵TX→Y高于Tsur(X)→Y即可認為是顯著的。通過100 次隨機的零假設抽樣計算Tsur(X)→Y,若根據(jù)置換序列計算出的新Tsur(X)→Y值的95%小于原來TX→Y,則可以認為原始序列計算得來的TX→Y顯著。

    1.2.3 算法比較

    運用時移相關系數(shù)、轉移熵分析行星際磁場南向分量Bz與地磁指數(shù)Sym-H之間的相關性。使用1996-2018 年小時精度數(shù)據(jù)對時移相關系數(shù)和轉移熵進行逐年計算。針對每一年的數(shù)據(jù),將時移相關系數(shù)rx→y(τ)=r[x(t),y(t+τ)]的最大時間延遲設置為120 h,以步長1 h 逐步計算,即每年共120 個數(shù)據(jù)點,然后取相關系數(shù)絕對值最大的數(shù)據(jù)點作為衡量該年Bz和地磁指數(shù)Sym-H之間的相關系數(shù)計算結果,有

    同理,計算轉移熵時采用同樣的時間延遲設置和最大值選取方法,由此可得圖1。圖1 中1996 年行星際磁場南向分量與地磁指數(shù)的最大時移相關系數(shù)r(Bz,Sym-H)=0.39;在太陽活動高年2000 年r(Bz,Sym-H)=0.35,同極小年的差異率為11%。另外,1996 年和2000 年行星際磁場南向分量與地磁Sym-H指數(shù)之間的轉移熵T(Bz →Sym-H)分別為0.18 nats、0.26 nats,差異率為31%。相關系數(shù)在太陽活動高低年的差異性相對較小,轉移熵在太陽活動高低年表現(xiàn)出明顯的差異。與傳統(tǒng)相關系數(shù)的對比表明,轉移熵可以有效表征太陽風–磁層關系的相關性。

    如圖1 所示,轉移熵的變化趨勢表現(xiàn)出類似太陽活動周變化的雙峰結構。太陽活動水平在一定程度上影響著空間天氣,太陽黑子數(shù)(Nss)是描述太陽活動水平的重要指標,計算得到太陽黑子和行星際磁場南向分量與地磁Sym-H指數(shù)之間的轉移熵TBz→Sym-H=T(Bz →Sym-H)的相關系數(shù)

    圖1 行星際磁場南向分量Bz 與地磁Sym-H 指數(shù)的轉移熵和相關系數(shù)Fig.1 Transfer entropy and correlation coefficient between I MF Bz and Sym-H

    太陽黑子和行星際磁場南向分量與地磁Sym-H指數(shù)之間的轉移熵TBz→Sym-H的轉移熵

    太陽黑子和行星際磁場南向分量與地磁Sym-H指數(shù)之間的相關系數(shù)rBz,Sym-H=r(Bz,Sym-H)的相關系數(shù):

    對比結果得知,

    以太陽活動周期的角度看,各太陽風參數(shù)向地磁Sym-H指數(shù)信息傳遞的強弱與太陽活動水平的周期性變化一致。另外,太陽10.7 cm 射電流量(F10.7)是綜合衡量太陽色球、過渡區(qū)和日冕極紫外輻射強度的一種常用指數(shù),以F10.7作為衡量太陽活動水平的指標可得

    以F10.7刻畫太陽活動強度,得到與太陽黑子數(shù)相同的結果。因此,轉移熵可以更好地定量化太陽活動與太陽風–磁層系統(tǒng)的相關性。

    2 太陽活動水平對太陽風–磁層相互作用的影響

    以往太陽風–磁層耦合關系多參數(shù)影響因素的重要性研究,往往僅集中在個別太陽風參數(shù),例如行星際磁場南向分量(Bz)或行星際電場(E),而與其他參數(shù),例如行星際磁場大?。˙)、太陽風速度(vsw)、太陽風數(shù)密度(Dsw)、太陽風溫度(Tsw)、太陽風動壓(Psw)等的綜合研究較少涉及。Schwenn 等[16]發(fā)現(xiàn),地磁擾動在大小和方向上都與行星際磁場波動密切相關。長期南向的行星際磁場與地球磁場之間的相互連接允許太陽風能量傳輸?shù)降厍虼艑覽17]。包括Kissinger等[18]在內(nèi)的多項研究表明了太陽風速度在磁暴產(chǎn)生中的作用,事實上持續(xù)和增強的太陽風速度以及南向或北向的行星際磁場分量通常與行星際物質拋射有關,其被認為是磁暴的重要原因。另外,增強的太陽風密度也是一個經(jīng)常影響磁暴強度的重要參數(shù)。太陽風密度的增大會導致日側磁層頂?shù)膲嚎s,從而導致磁層頂流、場向流和橫向流的增加。Ahmed 等[7]的研究表明太陽風密度與其他參數(shù)協(xié)同具有更好的相關性,而太陽風溫度則表現(xiàn)較差。Xie 等[19]多篇論文描述了高速太陽風動壓與磁暴的關系。Mcpherron等[20]的研究表明磁暴時的環(huán)電流積聚原則上可以直接由太陽風電場驅動。Kane[21]表示磁暴期間太陽風速度與行星際磁場南向分量vswBz的共同作用效果明顯優(yōu)于vsw;Wang 等[22]研究表明,對于磁暴的形成,的重要性強于行星際磁場南向分量增強的持續(xù)時間 ?t。

    標識地磁活動水平的典型地磁指數(shù),例如Dst指數(shù),表征了全球對稱赤道環(huán)電流的強度,是四個低緯天文臺的地球磁場水平分量的小時平均值[7],而Sym-H為1 min 分辨率的Dst指數(shù)。

    2.1 太陽風參數(shù)與地磁指數(shù)的轉移熵

    實驗使用1996-2018 年小時精度的太陽風數(shù)據(jù),主要分析參數(shù)包括行星際磁場(B)、行星際磁場南向分量(Bz)、太陽風速度(vsw)、太陽風等離子體數(shù)密度(Dsw)、太陽風溫度(Tsw)、太陽風動壓(Psw)、行星際電場(E),以及小時精度的地磁Sym-H指數(shù)數(shù)據(jù)。

    圖2 展示了利用轉移熵方法分析各太陽風參數(shù)與地磁Sym-H指數(shù)的相關性結果,各太陽風參數(shù)向地磁指數(shù)的轉移熵的相對值表現(xiàn)出一定的穩(wěn)定性,即:1996-2018 年中,行星際電場與地磁指數(shù)之間的信息傳遞TE→Sym-H均最高,隨后是和,最低是和位于中部且相互差異較小,排名在各年中略有不同。各太陽風參數(shù)與地磁Sym-H指數(shù)轉移熵的變化趨勢表現(xiàn)類似于圖1 中行星際磁場南向分量與地磁Sym-H指數(shù)的轉移熵相似的特性:(1)在太陽活動高年,Tsw→Sym-H更為顯著;(2)對于任一太陽風參數(shù),Tsw→Sym-H在第24 太陽活動周的峰值比23 太陽活動周要低。

    圖2 各太陽風參數(shù)對地磁Sym-H 指數(shù)的轉移熵Fig.2 Transfer entropy of multiple solar wind parameters to Sym-H

    2.2 轉移熵與太陽活動水平相關性分析

    如圖2 所示,除行星際磁場南向分量外,行星際磁場、太陽風速度、太陽風等離子體數(shù)密度、行星際電場及地磁Sym-H指數(shù)轉移熵表現(xiàn)出與太陽活動周類似的雙峰結構。同時將Tsw→Sym-H與太陽活動指數(shù)太陽黑子數(shù)、F10.7進行相關性分析,結果列于表1。

    由 表1 可 知,r(Nss,Tsw→Sym-H)位 于[0.61,0.83]之間、r(F10.7,Tsw→Sym-H)位于[0.58,0.82]之間,均表現(xiàn)出較強的相關性。同時獲得T(Nss,Tsw→Sym-H)位于[0.15,0.36]之間、T(F10.7,Tsw→Sym-H)位于[0.12,0.33]之間,即太陽活動指數(shù)與Tsw→Sym-H存在明顯的信息傳遞。由此可知,以小時精度的數(shù)據(jù)分析,行星際磁場南向分量與太陽活動之間的強相關并非偶然,其他各太陽風參數(shù)與太陽活動水平也存在較強的關聯(lián),佐證了Liu 等[23]對于地球近地空間環(huán)境變化對太陽活動的依賴性。

    表1 T sw→Sym-H同太陽活動指數(shù)的相關性Table 1 Correlation betweenTsw→Sym-H and solar radiation

    3 引起磁暴擾動的太陽風參數(shù)重要性排序

    3.1 數(shù)據(jù)選擇

    磁暴是指整個地球磁層發(fā)生的持續(xù)十幾個小時到幾十個小時的一種劇烈地磁擾動,中低緯度地磁臺站水平分量的顯著減小為磁暴的主要特征。Wu 等[24]的研究表明,太陽擾動日變化在主相期間最強,在恢復相期間逐漸減弱。在磁暴初相和主相期間,一般會出現(xiàn)太陽風–磁層發(fā)電機效應,將導致磁層對流電場快速滲透到中低緯度地區(qū);高緯焦耳加熱產(chǎn)生全球熱層擾動場,主要出現(xiàn)在恢復相期間[25]。磁暴主相期間,低緯水平磁場擾動(?H)在各磁地方時(MLT)均明顯下降并且出現(xiàn)晨昏不對稱分布;磁暴恢復相期間,行星際電場(E)為負值時,?H較主相期間弱。已有的研究表明:主相階段主要源于環(huán)電流的增強,而影響恢復相的主要為電荷損失機制,包括電荷交換、庫侖散射和波粒相互作用等[26]。因此,由于不同的物理成因,應分別研究磁暴主相和恢復相,目前本文將研究范圍限定于磁暴主相階段。

    一般根據(jù)Dst峰值的不同,將地磁暴分為小磁暴(?50 nT

    3.2 重要性排序方法

    以往研究結果表明,地磁指數(shù)對太陽風變化的響應時間普遍小于3 h[11,28],本研究將響應時間范圍設定在0 至500 min,與Stumpo 等[11]研究磁層電離層對太陽風的響應時間范圍一致。以往的轉移熵往往集中在單個時間序列,而磁暴期間的重要性系數(shù)為磁暴集合(見表2)上的太陽風參數(shù)對地磁Sym-H指數(shù)的轉移熵平均值,故而定義

    表2 2010-2018 年磁暴事件Table 2 Geomagnetic storm events from 2010 to 2018

    其中〈·〉代表在表2 磁暴集合上的平均值。

    相應地由轉移熵定義的重要性系數(shù)為

    以及定義信噪比(S/N),顯著性檢驗指標

    轉移熵算法通過置換時間序列[15]來創(chuàng)建零假設下的概率分布來計算噪聲,噪聲的均值和標準偏差是由100 次隨機的零假設抽樣計算Tsur(X)→Sym-H平均而來。

    3.3 重要性排序

    對行星際磁場南向分量(Bz)和地磁Sym-H指數(shù)之間的轉移熵進行計算,分析兩者之間的信息流動。首先對每個磁暴事件進行獨立的轉移熵分析,而后按照式(12)以相同的時間延遲τ計算多個磁暴事件轉移熵(τ)的平均值(τ)。由圖3 可以看出,轉移熵在τ=60 min時達到最大值Tpeak= 0.244 nats、Γ=0.196 nats,信噪比δS/N=4.15和顯著性αsig=36.83σ。轉移熵明顯大于噪聲,即信息傳遞顯著。而基本都處于背景噪聲的內(nèi)部或邊緣,不顯示任何明顯的時間延遲特征,信息傳遞不顯著。

    圖3 行星際磁層南向分量 Bz和 地磁Sym-H指數(shù)的轉移熵Fig.3 Transfer entropy of B z and Sym-H

    經(jīng)比較,上述結果與Runge 等[28]的結論基本一致,即行星際磁場南向分量(Bz)擾動是導致行星際條件變化的重要部分;Stumpo 等[11]的研究也表明,在時間延遲τ=60 min 時行星際磁場南向分量(Bz)與Sym-H指數(shù)之間有較強的信息傳遞。但是已有研究并未涉及太陽風溫度(Tsw)、行星際電場(E)等太陽風參數(shù)的綜合比較。進一步利用轉移熵分析磁暴期間其他太陽風參數(shù),包括行星際磁場(B)、太陽風速度(vsw)、太陽風等離子體數(shù)密度(Dsw)、太陽風溫度(Tsw)、太陽風動壓(Psw)和行星際電場(E),對Sym-H指數(shù)的影響,結果如圖4 所示。

    由圖4 可知,各太陽風參數(shù)對地磁Sym-H指數(shù)的轉移熵均高于噪聲;而則處于噪聲邊緣,缺乏顯著性。最大值出現(xiàn)的時間普遍小于3 h,反應出磁層對太陽風變異性的響應是突變式的激增,表現(xiàn)為短時間內(nèi)一個強烈的信息傳遞。表3 根據(jù)各太陽風參數(shù)到地磁Sym-H指數(shù)的轉移熵,得出了各太陽風參數(shù)的重要性系數(shù)序列。

    由表3 可知,以轉移熵構建的重要性系數(shù)為指標,排在第一位的是行星際電場(E),ΓE=0.200 nats。第二位是行星際磁場南向分量(Bz),ΓBz=0.196 nats,略低于ΓE。然后排名依次是太陽風速度(vsw)、太陽風溫度(Tsw)、太陽風密度(Dsw)、行星際磁場(B)。最低的是太陽風動壓(Psw),ΓPsw=0.135 nats。

    值得注意的是圖4(f),行星際電場(E)對地磁Sym-H指數(shù)的轉移熵在τ=60 min時達到最大值Tpeak= 0.252 nats、Γ=0.200 nats,信噪比δS/N= 3.92 和顯著性αsig=37.43σ,也有明顯的信息傳遞,甚至比更強烈。

    表3 中排在第三位的是太陽風速度(vsw),結合圖4(d)可知在τ=140 min 時達到最大值Tpeak= 0.210 nats、Γ=0.159 nats,信噪比δS/N=3.15和顯著性αsig=29.92σ。行星際磁場南向分量和磁層的重聯(lián)是引發(fā)磁暴等空間天氣事件的重要因素。Zhang 等[4]進一步得出,行星際磁場南向分量與太陽風速度的聯(lián)合作用又體現(xiàn)于行星際電場的時變曲線,統(tǒng)一表現(xiàn)為行星際電場,太陽風速度的變化影響著磁層對流電場,從而影響環(huán)電流的形成快慢與強度。因此,從信息傳遞的角度分析,行星際電場擾動是磁暴發(fā)生的主要驅動因素。以往太陽風溫度(Tsw)、太陽風密度(Dsw)、行星際磁場(B)、太陽風動壓(Psw)等綜合研究較少,從轉移熵的角度出發(fā),本文給出了這些太陽風參數(shù)在引發(fā)磁暴時的相對重要性。

    圖4 行星際磁場(B )、太陽風速度(vsw)、太陽風等離子體數(shù)密度(D sw )、太陽風溫度(Tsw)、太陽風動壓(Psw)、行星際電場(E)對地磁Sym-H 指數(shù)的轉移熵Fig.4 Transfer entropy of interplanetary magnetic field (B),solar wind speed (vsw),solar wind plasma number density (D sw) ,solar wind temperature (Tsw) solar wind dynamic pressure (Psw),the interplanetary electric field (E) and Sym-H

    表3 太陽風參數(shù)對地磁Sym-H 指數(shù)重要性系數(shù)排序Table 3 Importance ranking of solar wind parameters to geomagnetic index Sym-H

    以相關系數(shù)為基礎進行排序,由表3 中的|r|列可知,前三位仍是行星際電場(E)、行星際磁場南向分量(Bz)和太陽風速度(vsw),但是行星際磁場(B)的相對重要性提升至第四位,太陽風溫度(Tsw)下降至第六位。根據(jù)1.2 節(jié)的討論可知,相關系數(shù)僅能挖掘變量之間的線性關系,而從概率角度出發(fā)的轉移熵可以挖掘變量之間更一般的非線性依賴關系,轉移熵和相關系數(shù)排序的差異,表明了太陽風參數(shù)與地磁Sym-H指數(shù)之間關系的復雜性。進一步通過式(7)計算DM來對差異進行定量刻畫,DM >0 揭示太陽風–磁層之間線性關系之外更復雜的關聯(lián)性。由此可知,相對于相關系數(shù),轉移熵方法可以更多地挖掘參數(shù)間的非線性關系,故而導致了重要性排序的差異。

    綜上所述,行星際電場(E)的變化是引發(fā)地磁暴的主要驅動因素,這與Zhao 等[5]的研究結果基本一致。行星際磁場南向分量(Bz)與地磁Sym-H指數(shù)之間的重要性系數(shù)略低于ΓE,但是明顯高于其他太陽風參數(shù)的Γsw。太陽風動壓(Psw)與地磁Sym-H指數(shù)之間的重要性系數(shù)最低;其他各太陽風參數(shù)在引發(fā)磁暴時的相對關聯(lián)程度由表3 給出。

    4 結論

    應用轉移熵研究了地磁擾動的太陽風驅動因素,可以得出如下結論。

    (1)從太陽活動周期的角度分析,各太陽風參數(shù)向地磁Sym-H指數(shù)信息傳遞的強弱與太陽活動水平的周期性變化一致。

    (2)磁暴期間各太陽風參數(shù)向地磁Sym-H指數(shù)信息傳遞,最相關的變量是行星際電場(E),行星際磁場南向分量(Bz)次之,然后排名依次是太陽風速度(vsw)、太陽風溫度(Tsw)、太陽風等離子體數(shù)密度(Dsw)、行星際磁場(B),最弱的是太陽風動壓(Psw),各參數(shù)信息傳遞強弱排名見表3。

    (3)行星際電場與地磁指數(shù)之間的信息傳遞TE→Sym-H在τmax=60 min 處達到峰值,比TBz→Sym-H略高,佐證了Zhao 等[5]的研究。

    在研究太陽風–磁層系統(tǒng)響應太陽風條件變化時,運用轉移熵可以推斷太陽風–磁層的非線性關系。本文將轉移熵方法推廣至磁暴事件集合上的平均轉移熵,定義了磁暴期間太陽風參數(shù)與地磁Sym-H指數(shù)的重要性系數(shù),并綜合考慮了7 種太陽風參數(shù),以重要性系數(shù)為依據(jù)獲得了引發(fā)磁暴擾動的太陽風參數(shù)重要性排序。轉移熵定量描述太陽風與磁層參數(shù)之間的非線性信息傳遞,為探索磁層對行星際太陽風擾動的響應過程提供一條數(shù)據(jù)驅動的研究途徑,為模型研究提供參數(shù)選擇的依據(jù),由此可加速建模過程。未來轉移熵在空間天氣領域可以通過更多事件統(tǒng)計研究開展廣泛應用。

    致謝數(shù)據(jù)分析環(huán)境由國家科技資源共享服務平臺之一的國家空間科學數(shù)據(jù)中心(http://www.nssdc.ac.cn)提供。

    猜你喜歡
    磁層磁暴太陽活動
    小諸葛漫“話”天氣之極光與地磁暴
    科學24小時(2024年5期)2024-05-17 12:59:56
    航天器的“隱形殺手”:地磁暴
    軍事文摘(2023年14期)2023-07-28 08:39:58
    《從太陽風暴到行星磁暴》專刊(Ⅲ)寄語
    第24太陽活動周中國地區(qū)電離層閃爍統(tǒng)計特性研究
    基于對抗神經(jīng)網(wǎng)絡有限角度的磁層邊界CT重構技術
    第23和24太陽活動周高緯地磁感應電流分布特性
    磁鞘快速流引起的磁層頂凹陷事件
    等離子體鞘層效應對磁層探測電場儀設計的影響
    地球磁層亞暴統(tǒng)計分析
    磁暴恢復相后期超低頻波觀測分析
    夜夜骑夜夜射夜夜干| 丝袜在线中文字幕| 亚洲av男天堂| av网站免费在线观看视频| www.av在线官网国产| 亚洲av成人精品一区久久| 天堂俺去俺来也www色官网| 亚洲av福利一区| 成人免费观看视频高清| 少妇精品久久久久久久| 国产日韩欧美视频二区| 亚洲国产色片| 久久 成人 亚洲| 亚洲国产精品一区三区| 夫妻午夜视频| 少妇的逼好多水| 少妇人妻久久综合中文| 国产一区二区在线观看日韩| 欧美日韩在线观看h| 精品熟女少妇av免费看| 日韩,欧美,国产一区二区三区| 成人综合一区亚洲| av在线app专区| 亚洲国产成人一精品久久久| 少妇人妻精品综合一区二区| 在线 av 中文字幕| 国产在视频线精品| 国产精品久久久久久av不卡| 免费观看性生交大片5| 久久久久人妻精品一区果冻| 精品久久久久久久久av| av又黄又爽大尺度在线免费看| 精品国产一区二区久久| 亚洲美女黄色视频免费看| 亚洲av国产av综合av卡| 寂寞人妻少妇视频99o| 欧美日韩综合久久久久久| 亚洲在久久综合| 午夜视频国产福利| 久久久久久久亚洲中文字幕| 涩涩av久久男人的天堂| 欧美成人午夜免费资源| 97超视频在线观看视频| 国产 一区精品| 久久久午夜欧美精品| 国内少妇人妻偷人精品xxx网站| 亚洲成色77777| 国产欧美亚洲国产| 亚洲性久久影院| 另类亚洲欧美激情| 国产在线男女| 亚洲三级黄色毛片| 国产亚洲av片在线观看秒播厂| 亚洲成人一二三区av| 亚洲精品日本国产第一区| 国语对白做爰xxxⅹ性视频网站| 亚洲高清免费不卡视频| 精品卡一卡二卡四卡免费| 精品少妇内射三级| 在线观看免费高清a一片| 91精品国产国语对白视频| 中文字幕av电影在线播放| 久久青草综合色| 中文字幕av电影在线播放| 国产高清有码在线观看视频| 天堂8中文在线网| 午夜免费鲁丝| av天堂久久9| 18禁在线播放成人免费| 一级片'在线观看视频| 亚洲欧洲日产国产| 777米奇影视久久| 国产精品国产三级国产av玫瑰| 日日撸夜夜添| 日本vs欧美在线观看视频 | 成人特级av手机在线观看| 成人影院久久| 激情五月婷婷亚洲| 黑人猛操日本美女一级片| 黑人高潮一二区| 亚洲精品色激情综合| 大码成人一级视频| 最近最新中文字幕免费大全7| 国产欧美另类精品又又久久亚洲欧美| 久久久久视频综合| 新久久久久国产一级毛片| 男人狂女人下面高潮的视频| 欧美精品亚洲一区二区| 亚洲人与动物交配视频| 在线免费观看不下载黄p国产| 老司机影院成人| 国产片特级美女逼逼视频| 日韩 亚洲 欧美在线| 最近的中文字幕免费完整| 国产成人aa在线观看| 久久久久网色| 狂野欧美白嫩少妇大欣赏| 三级国产精品片| 美女视频免费永久观看网站| 色视频www国产| 免费看日本二区| 中文字幕免费在线视频6| 老熟女久久久| 在线观看一区二区三区激情| 青春草视频在线免费观看| 99热这里只有是精品50| 丝瓜视频免费看黄片| 香蕉精品网在线| 香蕉精品网在线| 亚洲图色成人| 又爽又黄a免费视频| 亚洲国产精品国产精品| 丁香六月天网| 免费黄频网站在线观看国产| 夜夜爽夜夜爽视频| 国产精品人妻久久久久久| 国产在线免费精品| 黑丝袜美女国产一区| 午夜福利在线观看免费完整高清在| 免费观看av网站的网址| 22中文网久久字幕| 久久久午夜欧美精品| 国产一区二区三区av在线| 精品久久久久久久久亚洲| 在线观看三级黄色| 美女脱内裤让男人舔精品视频| 全区人妻精品视频| av网站免费在线观看视频| 日本av免费视频播放| 日韩制服骚丝袜av| 大码成人一级视频| 成人影院久久| 免费观看无遮挡的男女| 美女内射精品一级片tv| 插阴视频在线观看视频| 精品一区二区三区视频在线| 乱系列少妇在线播放| 国产高清有码在线观看视频| 亚洲一级一片aⅴ在线观看| 天堂中文最新版在线下载| 久热这里只有精品99| 国产视频内射| 精品少妇内射三级| 国产高清国产精品国产三级| 日韩不卡一区二区三区视频在线| 国产女主播在线喷水免费视频网站| 日韩一区二区视频免费看| 国产日韩欧美视频二区| 国产黄片视频在线免费观看| 99久久中文字幕三级久久日本| av网站免费在线观看视频| av视频免费观看在线观看| 欧美日韩亚洲高清精品| 亚洲欧美清纯卡通| 亚洲欧洲精品一区二区精品久久久 | 国产精品成人在线| 亚洲怡红院男人天堂| 爱豆传媒免费全集在线观看| 国产精品三级大全| 久久99热这里只频精品6学生| 亚洲综合精品二区| 欧美xxxx性猛交bbbb| 国产69精品久久久久777片| 欧美成人午夜免费资源| 一级爰片在线观看| 两个人免费观看高清视频 | 啦啦啦在线观看免费高清www| 夫妻午夜视频| 偷拍熟女少妇极品色| 亚洲婷婷狠狠爱综合网| 麻豆成人午夜福利视频| 国产69精品久久久久777片| 老熟女久久久| 久热这里只有精品99| 伦理电影大哥的女人| 夜夜看夜夜爽夜夜摸| 国产极品天堂在线| 日韩精品有码人妻一区| 一级毛片 在线播放| a级毛片免费高清观看在线播放| 一区二区三区四区激情视频| 色婷婷久久久亚洲欧美| 男人狂女人下面高潮的视频| 欧美一级a爱片免费观看看| 国产伦理片在线播放av一区| 大陆偷拍与自拍| 亚洲一区二区三区欧美精品| 搡女人真爽免费视频火全软件| 啦啦啦在线观看免费高清www| 偷拍熟女少妇极品色| 久久6这里有精品| 国产精品熟女久久久久浪| 男人狂女人下面高潮的视频| 久久精品国产亚洲av天美| 欧美精品高潮呻吟av久久| 国产伦在线观看视频一区| 波野结衣二区三区在线| 亚洲欧美精品自产自拍| 国产一区二区三区av在线| 国产伦精品一区二区三区视频9| 极品人妻少妇av视频| 日韩不卡一区二区三区视频在线| 国产一区二区在线观看av| 亚洲av男天堂| 大码成人一级视频| a级毛片免费高清观看在线播放| 久久久久人妻精品一区果冻| 久久青草综合色| 久久国产亚洲av麻豆专区| 精品少妇内射三级| 国产精品嫩草影院av在线观看| 久久久久网色| 另类精品久久| 亚洲av不卡在线观看| 在线观看一区二区三区激情| 午夜视频国产福利| 免费看光身美女| 精品国产一区二区三区久久久樱花| 最近中文字幕2019免费版| 成人亚洲欧美一区二区av| 一级毛片黄色毛片免费观看视频| 日韩av不卡免费在线播放| 大片电影免费在线观看免费| 免费人妻精品一区二区三区视频| 日本黄色日本黄色录像| 国产高清三级在线| 狂野欧美激情性bbbbbb| 日本av手机在线免费观看| a级毛片免费高清观看在线播放| 一区在线观看完整版| 夜夜看夜夜爽夜夜摸| 国内少妇人妻偷人精品xxx网站| 日本-黄色视频高清免费观看| 亚洲国产精品999| 免费黄网站久久成人精品| 黑人巨大精品欧美一区二区蜜桃 | 我的女老师完整版在线观看| √禁漫天堂资源中文www| 夜夜骑夜夜射夜夜干| 欧美丝袜亚洲另类| 国产在线一区二区三区精| 欧美变态另类bdsm刘玥| 久久99一区二区三区| 欧美日韩在线观看h| 五月玫瑰六月丁香| 亚洲精品久久午夜乱码| 久久国内精品自在自线图片| 少妇粗大呻吟视频| 日韩 欧美 亚洲 中文字幕| 老司机影院毛片| 久久久久久久久久久久大奶| 黄色片一级片一级黄色片| 男女高潮啪啪啪动态图| 一边摸一边抽搐一进一出视频| 纯流量卡能插随身wifi吗| 午夜影院在线不卡| 亚洲精品美女久久av网站| 桃红色精品国产亚洲av| 欧美日韩福利视频一区二区| 国产激情久久老熟女| 精品高清国产在线一区| 国产亚洲一区二区精品| 国产亚洲精品第一综合不卡| 日本黄色日本黄色录像| 最近最新中文字幕大全免费视频| 亚洲欧美成人综合另类久久久| 亚洲精品国产色婷婷电影| 夜夜骑夜夜射夜夜干| 欧美日韩亚洲高清精品| 国产精品久久久人人做人人爽| 欧美激情高清一区二区三区| 制服诱惑二区| 男女国产视频网站| 亚洲av成人一区二区三| 麻豆乱淫一区二区| 久久亚洲国产成人精品v| 麻豆av在线久日| 国产欧美日韩综合在线一区二区| 精品第一国产精品| 日韩免费高清中文字幕av| 婷婷成人精品国产| 日韩 欧美 亚洲 中文字幕| 久久国产精品大桥未久av| 国产一区二区在线观看av| 黑人欧美特级aaaaaa片| 下体分泌物呈黄色| 国产亚洲精品第一综合不卡| av有码第一页| 一区二区av电影网| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产精品成人久久小说| 涩涩av久久男人的天堂| 成人亚洲精品一区在线观看| 国产欧美日韩一区二区精品| 亚洲欧洲日产国产| 两性午夜刺激爽爽歪歪视频在线观看 | 天堂中文最新版在线下载| 精品国产乱码久久久久久小说| 91大片在线观看| 一二三四在线观看免费中文在| 久久毛片免费看一区二区三区| 性色av乱码一区二区三区2| 国产av精品麻豆| 嫁个100分男人电影在线观看| 精品国产一区二区久久| 啦啦啦视频在线资源免费观看| 久久久精品区二区三区| 亚洲精品久久成人aⅴ小说| 永久免费av网站大全| 极品人妻少妇av视频| 99久久综合免费| 成人av一区二区三区在线看 | 精品少妇久久久久久888优播| 国产高清视频在线播放一区 | 国产亚洲午夜精品一区二区久久| 99久久综合免费| 国产亚洲一区二区精品| 在线观看免费视频网站a站| 日本a在线网址| 欧美黄色淫秽网站| 天天躁日日躁夜夜躁夜夜| 伊人久久大香线蕉亚洲五| 国产免费现黄频在线看| 亚洲人成电影免费在线| 啦啦啦中文免费视频观看日本| 麻豆av在线久日| 熟女少妇亚洲综合色aaa.| 免费在线观看视频国产中文字幕亚洲 | h视频一区二区三区| 国产免费福利视频在线观看| av网站在线播放免费| 一级黄色大片毛片| 美女高潮喷水抽搐中文字幕| 十八禁人妻一区二区| 国产男女内射视频| 高清在线国产一区| 精品少妇一区二区三区视频日本电影| h视频一区二区三区| 亚洲欧美一区二区三区黑人| 国产激情久久老熟女| 亚洲午夜精品一区,二区,三区| 国产精品香港三级国产av潘金莲| 视频在线观看一区二区三区| 高清在线国产一区| 一区二区三区激情视频| 日韩,欧美,国产一区二区三区| 亚洲免费av在线视频| 欧美日韩av久久| 国产精品二区激情视频| 国产精品成人在线| 91成年电影在线观看| 操美女的视频在线观看| 精品卡一卡二卡四卡免费| 欧美亚洲日本最大视频资源| 女性生殖器流出的白浆| 欧美日韩一级在线毛片| 欧美日韩国产mv在线观看视频| 丝袜喷水一区| 99国产精品一区二区蜜桃av | av视频免费观看在线观看| 一级片免费观看大全| 亚洲国产欧美一区二区综合| 一本色道久久久久久精品综合| 高清欧美精品videossex| 99国产极品粉嫩在线观看| 国内毛片毛片毛片毛片毛片| 91麻豆av在线| 激情视频va一区二区三区| 欧美黑人欧美精品刺激| 国产一区二区三区av在线| 黄色a级毛片大全视频| 麻豆乱淫一区二区| 黄片小视频在线播放| 亚洲色图 男人天堂 中文字幕| 午夜日韩欧美国产| 亚洲第一av免费看| av超薄肉色丝袜交足视频| 久久青草综合色| 国产精品麻豆人妻色哟哟久久| 中国美女看黄片| 久久热在线av| 老汉色∧v一级毛片| 国产亚洲欧美在线一区二区| 一本色道久久久久久精品综合| 欧美性长视频在线观看| 久久中文看片网| 欧美日韩福利视频一区二区| 国产亚洲精品一区二区www | 日韩大片免费观看网站| 99国产精品一区二区三区| av视频免费观看在线观看| 在线观看免费日韩欧美大片| 亚洲精品久久午夜乱码| 亚洲全国av大片| 满18在线观看网站| 淫妇啪啪啪对白视频 | av有码第一页| 桃花免费在线播放| 黑人巨大精品欧美一区二区mp4| 亚洲久久久国产精品| 老司机午夜福利在线观看视频 | 中国国产av一级| 999精品在线视频| 国产亚洲欧美在线一区二区| 女人精品久久久久毛片| 麻豆乱淫一区二区| 成人亚洲精品一区在线观看| 大码成人一级视频| 久久久精品区二区三区| 视频区欧美日本亚洲| 可以免费在线观看a视频的电影网站| 国产成+人综合+亚洲专区| 80岁老熟妇乱子伦牲交| 国产欧美日韩综合在线一区二区| www日本在线高清视频| 欧美午夜高清在线| 99久久精品国产亚洲精品| 亚洲一区中文字幕在线| 国产91精品成人一区二区三区 | 视频在线观看一区二区三区| 青青草视频在线视频观看| 国产亚洲av片在线观看秒播厂| av网站免费在线观看视频| 人人妻人人添人人爽欧美一区卜| 极品少妇高潮喷水抽搐| 中文欧美无线码| 欧美激情久久久久久爽电影 | 岛国在线观看网站| 久久av网站| 久久久久久亚洲精品国产蜜桃av| 亚洲中文字幕日韩| 日本a在线网址| 亚洲国产精品一区三区| 日本五十路高清| 国产黄色免费在线视频| 91精品伊人久久大香线蕉| 久久国产精品男人的天堂亚洲| 91国产中文字幕| 色婷婷久久久亚洲欧美| 日韩,欧美,国产一区二区三区| 成人国语在线视频| 麻豆乱淫一区二区| 一个人免费在线观看的高清视频 | 久久av网站| 香蕉国产在线看| 午夜福利影视在线免费观看| kizo精华| e午夜精品久久久久久久| 亚洲精品成人av观看孕妇| 五月开心婷婷网| 精品一区在线观看国产| 亚洲成人免费av在线播放| 午夜精品久久久久久毛片777| 国产精品99久久99久久久不卡| 啦啦啦 在线观看视频| 美女扒开内裤让男人捅视频| 女人被躁到高潮嗷嗷叫费观| 高清黄色对白视频在线免费看| 久久久久久久精品精品| 欧美在线黄色| 久久精品久久久久久噜噜老黄| 国产精品一二三区在线看| 亚洲九九香蕉| 久久99热这里只频精品6学生| h视频一区二区三区| 少妇人妻久久综合中文| 久久人人爽人人片av| 免费观看av网站的网址| 精品国内亚洲2022精品成人 | 麻豆av在线久日| 超碰成人久久| 国产一区二区激情短视频 | 久久亚洲精品不卡| 免费不卡黄色视频| 成人av一区二区三区在线看 | 在线十欧美十亚洲十日本专区| 美女高潮到喷水免费观看| 欧美精品av麻豆av| 美女午夜性视频免费| 国产欧美日韩综合在线一区二区| av免费在线观看网站| 欧美日韩亚洲综合一区二区三区_| 我的亚洲天堂| 人妻 亚洲 视频| 桃花免费在线播放| 亚洲国产成人一精品久久久| 免费看十八禁软件| 亚洲专区字幕在线| 亚洲人成77777在线视频| 黑人欧美特级aaaaaa片| 国产又爽黄色视频| 国产精品秋霞免费鲁丝片| 欧美日韩黄片免| 视频在线观看一区二区三区| 日韩免费高清中文字幕av| 国产一区二区激情短视频 | 80岁老熟妇乱子伦牲交| 婷婷丁香在线五月| av线在线观看网站| 制服诱惑二区| 久热爱精品视频在线9| 狠狠婷婷综合久久久久久88av| 淫妇啪啪啪对白视频 | 国产免费福利视频在线观看| netflix在线观看网站| av有码第一页| 精品国内亚洲2022精品成人 | 波多野结衣一区麻豆| 99热全是精品| 久久人妻福利社区极品人妻图片| 丁香六月欧美| 亚洲成av片中文字幕在线观看| 欧美大码av| 久久国产精品男人的天堂亚洲| 亚洲国产av新网站| 中文字幕制服av| 亚洲五月婷婷丁香| 色老头精品视频在线观看| 国产有黄有色有爽视频| 热re99久久国产66热| 女人久久www免费人成看片| 久久午夜综合久久蜜桃| 日韩电影二区| 久久久久视频综合| 久久精品国产a三级三级三级| 亚洲人成电影免费在线| 超碰97精品在线观看| 嫁个100分男人电影在线观看| 久热爱精品视频在线9| 成年美女黄网站色视频大全免费| 在线观看舔阴道视频| 国产成人精品在线电影| 黄色怎么调成土黄色| 欧美午夜高清在线| √禁漫天堂资源中文www| 成年动漫av网址| 大片电影免费在线观看免费| 免费少妇av软件| 亚洲综合色网址| 国产精品麻豆人妻色哟哟久久| 少妇人妻久久综合中文| 三上悠亚av全集在线观看| 国产黄频视频在线观看| 一级黄色大片毛片| 伊人亚洲综合成人网| 97人妻天天添夜夜摸| 久久久久久久久久久久大奶| 老司机在亚洲福利影院| 亚洲一码二码三码区别大吗| av电影中文网址| 大陆偷拍与自拍| 亚洲中文日韩欧美视频| 捣出白浆h1v1| av国产精品久久久久影院| 亚洲精品美女久久av网站| 免费在线观看日本一区| 亚洲精品粉嫩美女一区| 啦啦啦视频在线资源免费观看| 日韩,欧美,国产一区二区三区| 天堂8中文在线网| 人人妻人人爽人人添夜夜欢视频| 精品第一国产精品| 亚洲国产精品999| 国产一区二区 视频在线| 欧美黄色片欧美黄色片| 国产一区二区 视频在线| 亚洲久久久国产精品| 人人澡人人妻人| 深夜精品福利| 精品国产乱子伦一区二区三区 | 欧美大码av| 久久亚洲国产成人精品v| 黑人欧美特级aaaaaa片| 一本一本久久a久久精品综合妖精| 黄色a级毛片大全视频| 啦啦啦 在线观看视频| 亚洲中文日韩欧美视频| 最近中文字幕2019免费版| 精品久久久精品久久久| 久久人妻福利社区极品人妻图片| 女人爽到高潮嗷嗷叫在线视频| 黄色视频,在线免费观看| 天天躁日日躁夜夜躁夜夜| 国产有黄有色有爽视频| 高清欧美精品videossex| 亚洲精品中文字幕一二三四区 | 欧美乱码精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 夜夜骑夜夜射夜夜干| 少妇裸体淫交视频免费看高清 | 大陆偷拍与自拍| 亚洲精品国产av成人精品| 国内毛片毛片毛片毛片毛片| 亚洲精品国产av成人精品| 美女脱内裤让男人舔精品视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品中文字幕一二三四区 | 一边摸一边抽搐一进一出视频| 他把我摸到了高潮在线观看 | 大片电影免费在线观看免费| 欧美在线黄色| 国产亚洲午夜精品一区二区久久| 亚洲视频免费观看视频| 天堂8中文在线网| www.av在线官网国产| 精品久久久精品久久久| 欧美在线一区亚洲| 又紧又爽又黄一区二区| 后天国语完整版免费观看| 亚洲av成人不卡在线观看播放网 | 国产精品久久久久久精品古装| 中亚洲国语对白在线视频| 男人爽女人下面视频在线观看| 中国美女看黄片| 老司机靠b影院|