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

    暴時(shí)場(chǎng)向電流及磁層—電離層耦合過程

    2022-05-14 05:58:40
    關(guān)鍵詞:磁層磁暴電離層

    王 慧

    武漢大學(xué)電子信息學(xué)院 空間物理系,武漢 430072

    0 引言

    場(chǎng)向電流(field-aligned currents,FACs)是高緯電離層與磁層之間沿著磁力線流動(dòng)的電流,是磁層質(zhì)量、能量、動(dòng)量傳輸?shù)诫婋x層的重要途徑,在磁層—電離層耦合中扮演著重要角色.1908 年,挪威物理學(xué)家Kristian Birkeland 在研究北極圈極光時(shí)最早發(fā)現(xiàn)FACs,當(dāng)極光出現(xiàn)時(shí),磁力儀的指針會(huì)改變方向.Birkeland 認(rèn)為這意味著在上空有電流流動(dòng).他通過著名的Terella 球?qū)嶒?yàn),解釋說:“我們之所以能看到極光,是因?yàn)樗鼈兪莵碜蕴?yáng)的等離子體”.他理論推測(cè),極光是由電流沿著磁力線從外太空流向高層大氣形成的.是否真實(shí)存在場(chǎng)向電流的爭(zhēng)論持續(xù)了50 多年.隨著首顆人造衛(wèi)星Sputnik 發(fā)射升空,人類進(jìn)入太空時(shí)代后,1966 年FACs 終于被衛(wèi)星(1963 38C)磁場(chǎng)觀測(cè)所證實(shí)[雖然Zmuda等(1966)最初的解釋是磁流體波].60 多年來科學(xué)家最終證實(shí)了Birkeland 的預(yù)言.在1967 年的Birkeland 會(huì)議上,學(xué)者們一致提議將場(chǎng)向電流稱為“Birkeland 電流”,得到了國(guó)際地磁和大氣聯(lián)合會(huì)的認(rèn)可.Iijima 和Potemra(1978)根據(jù)Triad 衛(wèi)星磁場(chǎng)觀測(cè)數(shù)據(jù)得到FACs 的統(tǒng)計(jì)模式,此時(shí)距離FACs 提出已有70 年.如今,科學(xué)家們深信,地球空間許多重要的物理現(xiàn)象都與場(chǎng)向電流相關(guān).

    此后,隨著低軌衛(wèi)星的不斷發(fā)射升空,如S3-2、Dynamics Explorer-1(DE-1)、Magsat、ISEE 1 號(hào)和2 號(hào)、Viking、DE 2、FAST、?rsted、Magsat、DMSP、CHAMP、AMPER、ST-5、Swarm 等,空間磁場(chǎng)資料日益豐富.科學(xué)家們使用這些豐富的衛(wèi)星本地測(cè)量數(shù)據(jù),開展了大量的高緯FACs 研究,包括統(tǒng)計(jì)學(xué)分析了太陽(yáng)風(fēng)、行星際磁場(chǎng)(interplanetary magnetic field,IMF)對(duì)FACs 的影響,地磁擾動(dòng)期間(磁暴和亞暴期間)FACs 的動(dòng)力學(xué)演變特征,以及不同時(shí)間尺度FACs 的分布特征(例如Smiddy et al.,1980;Burch et al.,1983;Zanetti et al.,1983;Kelly et al.,1986;Ohtani et al.,1995;Weimer,2001;Eriksson et al.,2002;Papitashvili et al.,2002;Ohtani et al.,2005;Wang et al.,2005;Rother et al.,2007;Gjerloev et al.,2011;Ritter et al.,2013;Anderson et al.,2014;Le et al.,2016;Ohtani et al.,2018;王慧等,2020;Wang and Lühr,2021).這些細(xì)致深入的研究為我們理解場(chǎng)向電流的時(shí)空分布和產(chǎn)生機(jī)制提供了有價(jià)值的信息,具有重要的理論意義.

    本文總結(jié)了近年來國(guó)內(nèi)外關(guān)于場(chǎng)向電流的研究進(jìn)展,其中第1 節(jié)為場(chǎng)向電流的產(chǎn)生及其主要類型;第2 節(jié)為磁暴期間場(chǎng)向電流的變化特征;第3 節(jié)為亞暴期間場(chǎng)向電流的空間構(gòu)型;第4 節(jié)討論了場(chǎng)向電流的開放性問題,以供今后進(jìn)一步研究參考.

    1 場(chǎng)向電流的產(chǎn)生及其類型

    通常認(rèn)為,場(chǎng)向電流起源于磁層,來自太陽(yáng)風(fēng)與磁層之間的相互作用.在晨昏地方時(shí)扇區(qū)存在兩片明顯的大尺度FACs(如圖1 所示),靠極側(cè)的電流元為1 區(qū)(Region 1,R1)FACs,在晨側(cè)流入電離層,在昏側(cè)流出電離層.靠赤道側(cè)的電流元為2 區(qū)(Region 2,R2)FACs,流向與R1 FACs 相反(Iijima and Potemra,1976).R1 FACs 起源于磁層頂和磁尾瓣區(qū)域,R2 FACs 與極光電激流指數(shù)的變化密切相關(guān),源于內(nèi)磁層的環(huán)電流區(qū)域(Sato and Iijima,1979).Imajo 等(2018)利用磁層Arase 衛(wèi)星和電離層AMPERE 衛(wèi)星進(jìn)行協(xié)同觀測(cè),從觀測(cè)上證實(shí)了R2 場(chǎng)向電流來自內(nèi)磁層環(huán)電流區(qū).子夜前地方時(shí)扇區(qū)通常有三片電流系,電流從中心流出電離層,從兩側(cè)流入電離層,這與Harang 間斷不連續(xù)性有關(guān).地磁擾動(dòng)期,場(chǎng)向電流強(qiáng)度增強(qiáng),并向低緯側(cè)移動(dòng),電流區(qū)域擴(kuò)寬.亞暴期間,磁尾越尾電流片坍塌,電流通過R1 場(chǎng)向電流形成閉合回路,在晨側(cè)流入電離層,在昏側(cè)流出電離層,而場(chǎng)向電流在電離層通過西向極光電激流形成閉合通道.

    圖1 行星際磁場(chǎng)IMF By>0 時(shí)南半球場(chǎng)向電流的平均位置分布圖.黑色代表流入電離層的電流,白色代表流出電離層的電流(修改自Wang et al.,2005).MLAT為地磁緯度,MLT 為地磁地方時(shí).R0(1,2)表示區(qū)域0(1,2)場(chǎng)向電流Fig.1 Average location of upward (white) and downward(black) FACs peaks in the Southern Hemisphere for IMF By >0 orientation (modified from Wang et al.,2005)

    Wang 等(2005)利用CHAMP 衛(wèi)星磁場(chǎng)數(shù)據(jù),對(duì)南半球場(chǎng)向電流特征進(jìn)行全面、系統(tǒng)的統(tǒng)計(jì)學(xué)研究,重點(diǎn)考察了磁層源區(qū)的電流源和電壓源特性,發(fā)現(xiàn)場(chǎng)向電流峰值密度和Pedersen 電導(dǎo)率之間存在顯著的線性相關(guān)性(如圖2 所示).場(chǎng)向電流和電導(dǎo)率的相關(guān)系數(shù)不高(只有0.53),這可能是因?yàn)槲闹惺褂昧松闲泻拖滦蟹逯祱?chǎng)向電流的平均值,而不是特定的區(qū)域0、區(qū)域1 或區(qū)域2 場(chǎng)向電流;另外如果電導(dǎo)率經(jīng)驗(yàn)?zāi)P偷恼`差較大,也會(huì)影響相關(guān)分析結(jié)果.電離層電導(dǎo)率不僅影響電流強(qiáng)度,還影響其位置,在無陽(yáng)光照射時(shí),正午的電流向低緯移動(dòng)了2o左右.如果磁層動(dòng)力學(xué)過程扮演了電壓源的角色,那么場(chǎng)向電流應(yīng)該較強(qiáng)烈依賴于電離層電導(dǎo)率,另一方面,如果磁層為電流源,則磁層電場(chǎng)和電離層電導(dǎo)率之一將改變來適應(yīng)電流的變化.觀測(cè)表明,白天場(chǎng)向電流受電離層光照條件(即電導(dǎo)率)的影響,使該部分類似于電壓源.但是在夜晚時(shí)段,場(chǎng)向電流密度不受電離層光照條件的影響,使其類似于電流源.

    圖2 日側(cè)場(chǎng)向電流的平均值隨電離層Pedersen 電導(dǎo)率的散點(diǎn)分布圖(修改自Wang et al.,2005)Fig.2 Correlation analysis of the average density of fieldaligned current in the daytime and the solar irradiation induced ionospheric Pedersen conductivity (modified from Wang et al.,2005)

    已有研究表明,R1 FACs 密度高于R2 FACs,如Le 等(2010)發(fā)現(xiàn),在平靜或中等地磁擾動(dòng)情況下,晨昏側(cè)R1 和R2 FACs 始終處于非平衡狀態(tài),R1 FACs 密度一般大于R2 FACs.Wang 等(2005)發(fā)現(xiàn)晨昏兩側(cè)R1 FACs 對(duì)重聯(lián)電場(chǎng)的響應(yīng)要強(qiáng)于R2.Christiansen 等(2002)發(fā)現(xiàn)平靜期(擾動(dòng)期)R1 電流密度比R2 大約40%(約50%).Coxon 等(2014)利用AMPERE 衛(wèi)星,發(fā)現(xiàn)強(qiáng)重聯(lián)電場(chǎng)條件下R1 FACs 是R2 FACs 的1.15 倍,與跨極蓋電流回路以及亞暴電流楔有關(guān).王慧等(2020)利用CHAMP 衛(wèi)星數(shù)據(jù)首次研究了R2 電流可能大于R1電流的分布特征及其與行星際磁場(chǎng)的關(guān)系.結(jié)果表明,這類事件均發(fā)生在行星際磁場(chǎng)北向旋轉(zhuǎn)和亞暴恢復(fù)相期間(如圖3 所示).當(dāng)行星際磁場(chǎng)北向旋轉(zhuǎn)時(shí),極蓋電勢(shì)和R1 電流迅速下降,因?yàn)樗鼈冎苯优c磁層頂相連.由于R2 電流起源于內(nèi)磁層,對(duì)行星際磁場(chǎng)的響應(yīng)比R1 FACs 慢幾十分鐘,而R2電流在亞暴恢復(fù)階段可能處于發(fā)展階段.因此,R1電流很快減弱,但是R2 電流減弱得更為緩慢,所以會(huì)出現(xiàn)R2 電流高于R1 電流的情形,這與兩種電流的產(chǎn)生機(jī)制有關(guān).

    圖3 當(dāng)R2 場(chǎng)向電流高于R1 場(chǎng)向電流時(shí),行星際磁場(chǎng)Bz 分量和極光電激流指數(shù)AL 指數(shù)的時(shí)間變化圖.(a)晨側(cè);(b)昏側(cè).0 時(shí)刻表示CHAMP 衛(wèi)星觀測(cè)到R2>R1 FACs 的時(shí)刻(修改自王慧等,2020)Fig.3 When R2 field-aligned currents are higher than R1 field-aligned currents,the variations of IMF Bz and auroral current index AL.(a) In the dawn sector;(b) In the dusk sector.Time 0 represents the time when the CHAMP satellite observed R2 >R1 FACs event (modified from Wang et al.,2020)

    當(dāng)行星際磁場(chǎng)為北向時(shí),NBZ(northward IMFBz)FACs 出現(xiàn)在R1 FACs 的靠極側(cè),且流向與R1 電流相反(Iijima and Potemra,1976).當(dāng)IMFBy>0 時(shí),北半球昏側(cè)NBZ FACs 區(qū)域變窄,晨側(cè)NBZ FACs 區(qū)域擴(kuò)大,IMFBy<0 時(shí)晨昏變化相反,南 半 球 與 北 半 球 變 化 相 反(Iijima and Shibaji,1987).如果IMFBy 分量占主導(dǎo),正午時(shí)段電離層極尖區(qū)將出現(xiàn)兩片極性相反的FACs(Wilhjelm et al.,1978),如圖1 所示,稱作DPY(disturbance polar Y)FACs(Clauer and Friis-Christensen,1988).當(dāng)IMFBy>0(<0)時(shí),北半球靠極側(cè)的DPY FACs 流入(流出)電離層,而較低緯度的DPY FACs 極性相反.DPY 靠極側(cè)的電流元也被稱作極幔、極隙或區(qū)域0(Region 0,R0)FACs(Iijima and Potemra,1982;Erlandson et al.,1988),命名的差異反映了人們對(duì)DPY 電流的磁層源區(qū)的看法存在較多爭(zhēng)議.人們關(guān)于NBZ 和DPY FACs 之間的聯(lián)系有不同看法.通常認(rèn)為NBZ 和DPY FACs 是兩個(gè)獨(dú)立的電流類型,它們之間并無內(nèi)在聯(lián)系(例如Iijima and Shibaji,1987).然而,一部分學(xué)者注意到兩者之間密切的聯(lián)系,當(dāng)IMF 從北轉(zhuǎn)為昏向時(shí),NBZ FACs 可以演化為DPY FACs,人們進(jìn)一步指出IMF 方向不同會(huì)引起磁重聯(lián)位形的變化,從而導(dǎo)致這兩種電流元相互轉(zhuǎn)化(例如Fedder et al.,1997;Vennerstrom et al.,2002;Wang et al.,2008).Wang 等(2008)利 用 Space weather modeling framework 模型研究了冬至季節(jié)NBZ 和DPY 場(chǎng)向電流的轉(zhuǎn)換關(guān)系,考察了地球偶極傾角的影響.研究發(fā)現(xiàn),當(dāng)IMF 由北向逐漸轉(zhuǎn)為昏向時(shí),北半球和南半球NBZ 電流片都發(fā)生了逆時(shí)針方向轉(zhuǎn)動(dòng),而春分時(shí)節(jié)順時(shí)針方向轉(zhuǎn)動(dòng),南半球NBZ 電流最后變成DPY 電流,但北半球的晨側(cè)NBZ 電流變?yōu)镈PY R0 電流,而昏側(cè)NBZ 電流則轉(zhuǎn)到夜側(cè),并未形成DPY R1 電流(如圖4 所示).

    圖4 當(dāng)行星際磁場(chǎng)從北向轉(zhuǎn)為昏向時(shí),SWMF 模擬北半球電離層場(chǎng)向電流的空間分布(修改自Wang et al.,2008).UT 為世界時(shí)Fig.4 Field-aligned currents (μA/m2) in the Northern (winter) Hemisphere under varying northward and duskward IMF (modified from Wang et al.,2008)

    除上述電流元(即R1、R2、DPY、NBZ 電流)外,人們還發(fā)現(xiàn)了一種與IMFBx 分量相關(guān)的場(chǎng)向電流.徑向行星際磁場(chǎng)通常產(chǎn)生于行星際日冕物質(zhì)拋射的尾部,此時(shí)IMF 平行或反平行于太陽(yáng)風(fēng)速度(Neugebauer and Goldstein,1997).一般認(rèn)為,當(dāng)行星際磁場(chǎng)為徑向時(shí),磁層—電離層處于平靜狀態(tài),以往的模擬和觀測(cè)研究都忽視了IMFBx 的影響.但是,許多研究表明,當(dāng)行星際磁場(chǎng)為徑向時(shí),弓形激波、磁鞘和磁層頂區(qū)都會(huì)產(chǎn)生劇烈的擾動(dòng)(Verigin et al.,2001;Merka et al.,2003;Lin and Wang,2005;Hietala et al.,2009;Omidi et al.,2009;Shue et al.,2009;Amata et al.,2011;Hietala et al.,2012).Tang 等(2013)模擬結(jié)果表明徑向行星際磁場(chǎng)期間磁層頂會(huì)出現(xiàn)開爾文—亥姆霍茲不穩(wěn)定性.Wang 等(2014)利用CHAMP 衛(wèi)星觀測(cè)數(shù)據(jù)揭示,在徑向行星際磁場(chǎng)中,日側(cè)極隙區(qū)有較強(qiáng)的FACs,流向與DPY 電流相反,而在子夜附近并不出現(xiàn)增強(qiáng)的電流(如圖5 所示).研究表明,當(dāng)IMF 呈徑向分布時(shí),可與地球磁尾瓣開放磁力線重聯(lián),形成向日對(duì)流和場(chǎng)向電流元.Shi 等(2020)利用AMPERE 衛(wèi)星觀測(cè)數(shù)據(jù),也發(fā)現(xiàn)在IMFBx 分量占主導(dǎo)期間,極尖區(qū)場(chǎng)向電流有所增強(qiáng).

    圖5 徑向IMF 期間,CHAMP 衛(wèi)星觀測(cè)的北半球場(chǎng)向電流隨世界時(shí)和地磁緯度的分布.(a)子夜地方時(shí)扇區(qū);(b)正午地方時(shí)扇區(qū).正值代表留出電離層的場(chǎng)向電流,負(fù)值代表流入電離層的場(chǎng)向電流(修改自Wang et al.,2014).MLT 為地磁地方時(shí)Fig.5 Latitudinal versus universal time variations of FACs observed by CHAMP in the nighttime (a) and daytime (b) in the Northern Hemisphere.Positive values represent upward FACs and negative values represent downward FACs (modified from Wang et al.,2014)

    以上各類場(chǎng)向電流均起源于磁層,而科學(xué)家們利用數(shù)值模型提出了另一類產(chǎn)生于電離層電導(dǎo)率非均勻分布的場(chǎng)向電流.這種場(chǎng)向電流位于電離層晨昏線處,稱作界限電流(terminator current,Zhu et al.,2014).在極區(qū)日側(cè)有光照區(qū)域,電離層電導(dǎo)率較高,極光橢圓區(qū)域電導(dǎo)率也比較高,其他區(qū)域電導(dǎo)率較低,因此電離層電導(dǎo)率有明顯的空間梯度,一般從夜側(cè)指向白天側(cè).根據(jù)UTAH 大學(xué)的數(shù)據(jù)同化模型(ionospheric dynamics and electrodynamics with data assimilation,IDED-DA),Zhu 等(2014)發(fā)現(xiàn)在晨昏邊界上有一片沿電導(dǎo)率終端分布的場(chǎng)向電流,流出電離層,這種電流起源于電離層電導(dǎo)率的空間梯度.研究發(fā)現(xiàn),這種界限電流在冬季流出電離層,夏季流入電離層,主要出現(xiàn)在地磁平靜期,磁暴期間幾乎無法分辨,冬季電流強(qiáng)度與太陽(yáng)活動(dòng)成正比,夏季較高的太陽(yáng)活動(dòng)導(dǎo)致電流片變寬.

    2 磁暴期間場(chǎng)向電流的分布特征

    磁暴期間,整個(gè)地球空間會(huì)發(fā)生一系列擾動(dòng).近十多年的研究表明,較強(qiáng)而持續(xù)的南向IMFBz使得巨大的太陽(yáng)風(fēng)能量得以通過磁重聯(lián)機(jī)制進(jìn)入地球空間,從而引發(fā)磁暴.隨著太陽(yáng)風(fēng)能量的注入,場(chǎng)向電流會(huì)發(fā)生巨大改變,電流中心向赤道側(cè)移動(dòng),電流強(qiáng)度顯著增強(qiáng)(例如Wang et al.,2006).但是,不同磁暴事件期間場(chǎng)向電流的位置和強(qiáng)度的變化規(guī)律及其與太陽(yáng)風(fēng)行星際磁場(chǎng)的聯(lián)系存在較大差別.Anderson 等(2002)利用Iridium 衛(wèi)星研究了兩次磁暴期間(1999 年9 月22 日至23 日和10 月21日至22 日)場(chǎng)向電流的空間分布,發(fā)現(xiàn)R1 和R2 FACs 的位置與Kp 成反比.Liou 等(2005)研究了2002 年4 月17 日磁暴事件,發(fā)現(xiàn)場(chǎng)向電流在磁暴初期向低緯移動(dòng)了5°,在主相期間向低緯移動(dòng)了5°~6°.Wang 等(2006)調(diào)查了2003 年兩次大磁暴期間南北半球的FACs 強(qiáng)度及其緯度位置的演變,發(fā)現(xiàn)磁暴各個(gè)階段日夜側(cè)的變化并不相同,向陽(yáng)側(cè)的低緯擴(kuò)展直接受南向IMFBz 控制;而在背陽(yáng)側(cè),這一擴(kuò)展及恢復(fù),比行星際參數(shù)變化滯后約3 個(gè)小時(shí),但與表征磁層環(huán)電流Dst 指數(shù)幾乎同步變化.磁暴主相期間,背陽(yáng)側(cè)強(qiáng)場(chǎng)向電流緯度分布范圍大大擴(kuò)展,可達(dá)25°以上,并出現(xiàn)多達(dá)10片以上的多電流片結(jié)構(gòu)(如圖6所示).Zou 等(2014)研究了6 個(gè)磁暴,發(fā)現(xiàn)場(chǎng)向電流的位置與SYM-H 的變化一致,在磁暴期間向低緯移動(dòng)了5°~10°.Xiong 等(2014)研究了場(chǎng)向電流所在極光橢圓的位置與重聯(lián)電場(chǎng)、AL 和Dst 指數(shù)的關(guān)系.Lyon 等(2016)研究了2013 年St.Patrick 磁暴,發(fā)現(xiàn)南向IMF 期間R1 和R2 電流向低緯移動(dòng)1°~2°,在另一個(gè)南向IMF 期間則向低緯移動(dòng)2°~3°,且出現(xiàn)小尺度多片電流元結(jié)構(gòu).Le 等(2016)對(duì)2015 年3 月17 日磁暴進(jìn)行了研究,結(jié)果表明,持續(xù)時(shí)間較長(zhǎng)的南向IMF 導(dǎo)致場(chǎng)向電流向低緯度移動(dòng)8°,而電流的最低緯度出現(xiàn)在最強(qiáng)南向IMF 期間.Lukianova(2020)利用Swarm 衛(wèi)星數(shù)據(jù)研究了2017 年9 月份磁暴事件,發(fā)現(xiàn)電流的赤道邊界受控于SYM-H 指數(shù),夜晚電流密度變化較大,尤其是亞暴發(fā)生期間,昏側(cè)R2 電流密度強(qiáng)于晨側(cè)電流,可能與部分環(huán)電流有關(guān).

    圖6 2003 年10 月29 日至31 日大磁暴期間向陽(yáng)側(cè)場(chǎng)向電流隨磁暴時(shí)間和磁緯度的變化圖.(a)北半球;(b)南半球(修改自Wang et al.,2006)Fig.6 Latitudinal and storm time variation of the daytime FACs in comparison with the IMF Bz during the storm on October 29-31,2003.(a) The northern hemisphere;(b) The southern hemisphere (modified from Wang et al.,2006)

    南向IMF 期間,場(chǎng)向電流向低緯移動(dòng),主要是由于日側(cè)磁重聯(lián)引起極蓋區(qū)磁通量增加,極蓋面積增大.亞暴的發(fā)生通常與南向IMF 的北向轉(zhuǎn)變有關(guān),因此在亞暴發(fā)生之前場(chǎng)向電流也會(huì)向低緯移動(dòng).Coxon 等(2014)采用時(shí)序疊加方法,研究發(fā)現(xiàn),在亞暴發(fā)生前2 小時(shí)后場(chǎng)向電流的平均位置向低緯移動(dòng)了1°~2°,15 min 后達(dá)到最低緯度,亞暴發(fā)生后,極光橢圓開始收縮,電流向極側(cè)移動(dòng).主要是因?yàn)槿諅?cè)磁重聯(lián)增加引起極蓋面積增大,而磁尾磁重聯(lián)增加導(dǎo)致極蓋面積減少.磁暴過程中,場(chǎng)向電流的強(qiáng)度也將增大,通常伴隨著電流向赤道方向移動(dòng),當(dāng)?shù)竭_(dá)最低緯度時(shí)達(dá)到峰值強(qiáng)度(Coxon et al.,2014).一些研究表明,R1 和R2 場(chǎng)向電流的強(qiáng)度受控于日—夜側(cè)磁重聯(lián)過程(Milan et al.,2015).磁暴期間,場(chǎng)向電流的強(qiáng)度呈南北半球不對(duì)稱性,北半球的電流強(qiáng)度大于南半球,這可能與地磁場(chǎng)的構(gòu)型有關(guān).在磁暴主相期間,場(chǎng)向電流的強(qiáng)度可發(fā)生飽和現(xiàn)象(Anderson and Korth,2007),這是因?yàn)镽1 場(chǎng)向電流的增加,減弱了日側(cè)磁場(chǎng),從而降低磁重聯(lián)率.Wang 等(2006)發(fā)現(xiàn)當(dāng)南向IMFBz 小于-30 nT 時(shí),場(chǎng)向電流朝赤道的擴(kuò)展也會(huì)出現(xiàn)非線性飽和,且冬季半球飽和緯度高于夏季半球;磁暴期間場(chǎng)向電流密度比磁靜期要高約5 倍,南(夏季)半球比北(冬季)半球高約2.5 倍;進(jìn)一步發(fā)現(xiàn)沿緯度積分的場(chǎng)向電流密度主要受太陽(yáng)風(fēng)動(dòng)力學(xué)壓強(qiáng)而不是行星際磁場(chǎng)的控制.Lukianova(2020)也報(bào)道了磁暴期間場(chǎng)向電流類似的飽和效應(yīng).

    磁暴期間,場(chǎng)向電流向赤道側(cè)移動(dòng)的幅度,通常呈現(xiàn)晨昏不對(duì)稱性.Anderson 等(2005)統(tǒng)計(jì)學(xué)發(fā)現(xiàn)磁暴主相期間,昏側(cè)電流比晨側(cè)位于更低緯度,相差2.4°,可能與IMFBy 分量有關(guān),也可能與部分環(huán)電流有關(guān),也可能與昏側(cè)亞極光區(qū)極化流有關(guān).Green 等(2009)發(fā)現(xiàn)這一現(xiàn)象在夏季比冬季更明顯.場(chǎng)向電流的平均位置可以代表極光橢圓的位置,因此可用來表征日—夜側(cè)磁場(chǎng)重聯(lián)率的相對(duì)強(qiáng)弱(Clausen et al.,2016).磁暴期間,小尺度場(chǎng)向電流能引起電離層電子密度不規(guī)則結(jié)構(gòu)(Liu et al.,2021).上行場(chǎng)向電流中的高能電子沉降可以產(chǎn)生極光,但昏側(cè)的R1 場(chǎng)向電流是否與極光有關(guān),還有待進(jìn)一步研究(Carter et al.,2016).

    3 亞暴期間場(chǎng)向電流的分布特征

    亞暴期間,磁尾中性片電流會(huì)發(fā)生中斷,通過場(chǎng)向電流和電離層電流形成回路,子夜后向下流入電離層,子夜前向上流出電離層.場(chǎng)向電流與電離層夜側(cè)的西向電激流相連,構(gòu)成亞暴電流楔(McPherron et al.,1973).這個(gè)單亞暴電流楔(substorm current wedge,SCW)模型提出40 多年后得到了更新,在R1 電流靠地球一側(cè)加入R2 電流,稱為雙SCW 模型(Sergeev et al.,2011).磁尾高速等離子體流將磁通量從重聯(lián)區(qū)傳輸?shù)絻?nèi)磁層,改變子夜附近的磁場(chǎng),到了近地空間,向東西方向運(yùn)動(dòng),導(dǎo)致磁通量管出現(xiàn)扭曲,從而在遠(yuǎn)離地區(qū)一側(cè)形成R1 場(chǎng)向電流,在靠地球一側(cè)形成R2 場(chǎng)向電流.在磁赤道平面,R2 電流可以通過徑向電流與R1 電流形成回路,也可以通過部分環(huán)電流與R2 電流自身形成回路.Wang 等(2006)統(tǒng)計(jì)了亞暴初始位置附近上行和下行場(chǎng)向電流與亞暴初始的相對(duì)位置分布,結(jié)果表明,最接近亞暴初始位置的為上行 FACs,在靠極側(cè)和靠赤道側(cè)為向下的電流.然而,以往的研究對(duì)于下行FACs 的位置以及上行和下行FACs 之間的間距究竟有多大一直沒有統(tǒng)一的結(jié)論,不同的亞暴事件結(jié)果也并不相同.

    人們關(guān)于單、雙SCW 模型仍存在一定爭(zhēng)議.Murphy 等(2013)使用AMPERE 數(shù)據(jù),對(duì)三個(gè)亞暴事件進(jìn)行了研究,其中R1 和R2 電流在其中兩個(gè)亞暴事件中均有增強(qiáng),但在第三個(gè)亞暴事件中沒有明顯變化.Clausen 等(2013)研究了場(chǎng)向電流隨亞暴的變化特征,發(fā)現(xiàn)白天R1 電流雖然在亞暴發(fā)生之前占主導(dǎo)地位,但在亞暴發(fā)生后,夜側(cè)R1 電流占主導(dǎo)地位,R2 電流沒有增強(qiáng),因此支持單SCW理論.Coxon 等(2014)使用時(shí)序疊加法,對(duì)3 年的AMPERE 數(shù)據(jù)進(jìn)行分析,發(fā)現(xiàn)亞暴初始之后10 min的R1 和R2 電流開始向低緯擴(kuò)展,均顯著增強(qiáng),兩者的比值也有所增加,因此支持雙SCW 理論.

    亞暴期間實(shí)際的FACs 結(jié)構(gòu)可能比等效SCW模型更為復(fù)雜.Murphy 等(2013)發(fā)現(xiàn)周期性的上行和下行FACs 出現(xiàn)在子夜附近約8 個(gè)小時(shí)的地方時(shí)扇區(qū)內(nèi),當(dāng)沿緯度對(duì)FACs 進(jìn)行平均后,電流結(jié)構(gòu)與SCW 比較一致,且上行的FACs 與離散極光密切相關(guān).Coxon 等(2017)的時(shí)序疊加分析表明,亞暴初始附近的電流在亞暴發(fā)生之前減弱,這與亞暴初始之前極光強(qiáng)度減弱一致,表明兩者的緊密聯(lián)系.Forsyth 等(2014)利用Cluster 數(shù)據(jù)發(fā)現(xiàn)SCW是由很多小尺度楔狀流組成的,楔狀流與磁尾多個(gè)向地快速離子流有關(guān).Liu 等(2015)利用THEMIS衛(wèi)星統(tǒng)計(jì)研究了亞暴期間小尺度楔狀流,發(fā)現(xiàn)黎明時(shí)段朝向地球的電流比遠(yuǎn)離地球的FACs 多,因此凈電流朝向地球;黃昏時(shí)段情形正好相反,這些小尺度楔狀流結(jié)合形成大尺度的SCW.Nakamura 等(2016)發(fā)現(xiàn)等離子體片邊界層內(nèi)存在小尺度場(chǎng)向電流,尺度只有幾十千米,遠(yuǎn)低于離子尺寸,因此主要源于磁尾磁重聯(lián)之后的向地電子流,而不是離子流.目前還沒有明確亞暴期間大尺度SCW 和小尺度楔狀流究竟是哪一種占優(yōu)勢(shì),而且無法定量理解小尺度楔狀流對(duì)SCW 的貢獻(xiàn)有多大.Ohtani 和Gjerloev(2020)發(fā)現(xiàn)中緯度地磁擾動(dòng)與AL 指數(shù)有很好的相關(guān)性,因而認(rèn)為大尺度SCW 應(yīng)為全球性的分布,小尺度楔狀流不是主導(dǎo)因素.

    在磁暴期間,亞暴電流楔的中心可向晨側(cè)移動(dòng),即兩片R1 場(chǎng)向電流都出現(xiàn)在清晨時(shí)段,這與越尾電流的崩塌位置向晨側(cè)偏移有關(guān).Ohtani 等(2018)發(fā)現(xiàn)磁暴期間晨側(cè)出現(xiàn)單SCW 現(xiàn)象,冬季較容易發(fā)生,這主要是因?yàn)槎緲O蓋區(qū)電導(dǎo)率偏低,導(dǎo)致晨側(cè)下行凈場(chǎng)向電流增強(qiáng).另一個(gè)影響因素是IMFBy,當(dāng)IMFBy<0 時(shí)晨側(cè)SCW 更容易發(fā)生,這是由于在晨向IMF 情形下,正午下行的凈場(chǎng)向電流的中心向晨側(cè)偏移,導(dǎo)致晨側(cè)R1 電流增強(qiáng).

    另有研究表明,夜側(cè)存在兩個(gè)SCW 系統(tǒng),分別發(fā)生在子夜前和子夜后.子夜前的西向電激流相對(duì)于子夜后的電激流處于更高緯度(Gjerloev and Hoffman,2014).Wang 和Lühr(2021)統(tǒng)計(jì)研究了亞暴期間極光電激流和兩側(cè)場(chǎng)向電流的相對(duì)位置分布,得到夜側(cè)極光電激流的位置的最可幾分布圖(如圖7 所示).從圖中可以看出子夜后電流比子夜前處于更低的位置,從這點(diǎn)上看,似乎支持Gjerloev 和Hoffman(2014)的觀點(diǎn).

    圖7 夜側(cè)極光電激流的事件數(shù)隨地磁緯度的分布.Orbit-1 表示亞暴發(fā)生之前的軌道,Orbit0-2 表示亞暴發(fā)生之后三個(gè)軌道.不同顏色代表不同地方時(shí)扇區(qū).每個(gè)地方時(shí)扇區(qū)的平均值已經(jīng)給出(修改自Wang et al.,2021)Fig.7 The event number of the peak auroral electrojet as a function of MLAT at four subsequent orbits in six different MLT sections in the nighttime.The average MLat in each MLT is shown in the subfigures (modified from Wang and Lühr,2021)

    通常認(rèn)為SCW 電流通過西向極光電激流在電離層中形成閉合回路.然而由于電離層電導(dǎo)率的快速變化,電離層中會(huì)出現(xiàn)多片場(chǎng)向電流.除西向電激流外,還能形成子午向電流.也就是說,除了向上流出和向下流入電離層的R1 FACs 之外,在極光橢圓靠極側(cè)出現(xiàn)了另外的與R1 電流流向相反的R0電流.這樣,向上的R1 場(chǎng)向電流,除了通過西向電激流和向下的R1 電流形成遠(yuǎn)程閉合回路之外,還和極光橢圓靠極側(cè)邊界的R0 電流在子午方向上形成一個(gè)當(dāng)?shù)亻]合回路.Wang 和Lühr(2021)利用CHAMP 衛(wèi)星近10 年的極光電激流和場(chǎng)向電流的觀測(cè)數(shù)據(jù),統(tǒng)計(jì)學(xué)研究了亞暴時(shí)期兩者的空間相對(duì)分布,發(fā)現(xiàn)遠(yuǎn)程西向電激流與經(jīng)向電流的比值在亞暴期間增加了250%,表明亞暴期間遠(yuǎn)程通道更為有效.

    4 未來展望

    本文對(duì)極區(qū)場(chǎng)向電流的類型和地磁擾動(dòng)期間場(chǎng)向電流的時(shí)空變化特征進(jìn)行了系統(tǒng)總結(jié)和回顧.雖然已經(jīng)有數(shù)十年的研究歷史,且在電流的分布、起源和時(shí)空演化等方面已經(jīng)取得了重要的研究進(jìn)展,但是由于觀測(cè)手段的限制,場(chǎng)向電流目前仍然存在著許多尚未解決的問題.下面列舉幾個(gè)方面的問題,為希望開展此方面工作的讀者提供一些參考.

    (1)場(chǎng)向電流的南北半球和晨昏地方時(shí)不對(duì)稱性問題:磁暴期間,場(chǎng)向電流的強(qiáng)度呈南北半球不對(duì)稱性,昏側(cè)比晨側(cè)位于更低地磁緯度,其中具體物理機(jī)制并不清楚.

    (2)場(chǎng)向電流的小尺度波狀結(jié)構(gòu)的形成機(jī)理:正午時(shí)段極尖—極隙區(qū)小尺度場(chǎng)向電流的時(shí)空分布精度及其磁層起源并不清晰,是否與磁層頂?shù)亩啻未胖芈?lián)過程有關(guān)?是否與磁層電磁波動(dòng)過程有關(guān)?

    (3)亞暴電流楔中的小尺度電流楔問題:亞暴電流楔中存在小尺度電流,其產(chǎn)生機(jī)制與時(shí)空分布和演化并不清晰.是否與磁尾高速流的快速變化有關(guān)?是否與磁尾的電磁波動(dòng)過程有關(guān)?

    這些科學(xué)問題需要結(jié)合更豐富的觀測(cè)手段和更先進(jìn)的模擬技術(shù),進(jìn)行更加細(xì)致、深入的研究.建議利用Cluster、THEMIS、Van Allen Probes、MMS等磁層衛(wèi)星結(jié)合Swarm、AMPERE 等低軌道電離層衛(wèi)星的協(xié)同觀測(cè),利用地球空間磁層—電離層—熱層三維數(shù)值模型,研究不同尺度的場(chǎng)向電流的分布規(guī)律和形成機(jī)制,有助于人們了解太陽(yáng)風(fēng)能量向近地空間的傳播、轉(zhuǎn)化和耗散過程.

    致謝

    感謝德國(guó)CHAMP 和Swarm 衛(wèi)星團(tuán)隊(duì)提供的磁場(chǎng)數(shù)據(jù)、美國(guó)密歇根大學(xué)提供SWMF 模型、得克薩斯大學(xué)空間科學(xué)中心提供DMSP數(shù)據(jù),以及NASA/GSFC OMNIWeb(https://omniweb.gsfc.nasa.gov)提供太陽(yáng)風(fēng)和行星際磁場(chǎng)的數(shù)據(jù).

    猜你喜歡
    磁層磁暴電離層
    小諸葛漫“話”天氣之極光與地磁暴
    航天器的“隱形殺手”:地磁暴
    軍事文摘(2023年14期)2023-07-28 08:39:58
    《從太陽(yáng)風(fēng)暴到行星磁暴》??á螅┘恼Z(yǔ)
    一種電離層TEC格點(diǎn)預(yù)測(cè)模型
    基于對(duì)抗神經(jīng)網(wǎng)絡(luò)有限角度的磁層邊界CT重構(gòu)技術(shù)
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    磁鞘快速流引起的磁層頂凹陷事件
    等離子體鞘層效應(yīng)對(duì)磁層探測(cè)電場(chǎng)儀設(shè)計(jì)的影響
    地球磁層亞暴統(tǒng)計(jì)分析
    磁暴恢復(fù)相后期超低頻波觀測(cè)分析
    日韩欧美精品免费久久| av一本久久久久| 亚洲精品国产一区二区精华液| 亚洲激情五月婷婷啪啪| 90打野战视频偷拍视频| 午夜福利一区二区在线看| 亚洲国产成人一精品久久久| 午夜日韩欧美国产| 国产高清国产精品国产三级| 多毛熟女@视频| 色综合欧美亚洲国产小说| 亚洲精品aⅴ在线观看| 综合色丁香网| 在线观看三级黄色| 一级毛片 在线播放| 狂野欧美激情性xxxx| 色吧在线观看| 99热国产这里只有精品6| 叶爱在线成人免费视频播放| 国产成人一区二区在线| 日韩欧美一区视频在线观看| 日本欧美视频一区| 天天躁夜夜躁狠狠躁躁| 欧美成人午夜精品| 亚洲婷婷狠狠爱综合网| 人人澡人人妻人| 久久综合国产亚洲精品| 男女高潮啪啪啪动态图| 操出白浆在线播放| 看非洲黑人一级黄片| 亚洲av中文av极速乱| 国产精品一国产av| 亚洲美女视频黄频| 亚洲欧洲国产日韩| 亚洲美女黄色视频免费看| 精品福利永久在线观看| 这个男人来自地球电影免费观看 | 中国三级夫妇交换| 高清av免费在线| 别揉我奶头~嗯~啊~动态视频 | 九色亚洲精品在线播放| 老司机影院成人| 亚洲 欧美一区二区三区| 久久国产亚洲av麻豆专区| 亚洲 欧美一区二区三区| 久久韩国三级中文字幕| 欧美日韩精品网址| 欧美日韩av久久| 国产精品99久久99久久久不卡 | 亚洲国产欧美一区二区综合| 在线观看免费日韩欧美大片| 免费观看a级毛片全部| 欧美日韩一级在线毛片| 日韩一卡2卡3卡4卡2021年| xxx大片免费视频| 久久午夜综合久久蜜桃| 热re99久久国产66热| 老司机深夜福利视频在线观看 | 亚洲少妇的诱惑av| 国产成人欧美| 国产在线一区二区三区精| 晚上一个人看的免费电影| 国产在线一区二区三区精| 乱人伦中国视频| 国产精品麻豆人妻色哟哟久久| 午夜影院在线不卡| 精品一品国产午夜福利视频| 国产精品二区激情视频| 国产免费现黄频在线看| 五月开心婷婷网| 免费女性裸体啪啪无遮挡网站| 欧美日韩亚洲综合一区二区三区_| 午夜福利网站1000一区二区三区| 久久久久久久久久久久大奶| 如日韩欧美国产精品一区二区三区| 色婷婷av一区二区三区视频| 中文字幕色久视频| 日韩av免费高清视频| 午夜久久久在线观看| 男女国产视频网站| 久久热在线av| 一级毛片我不卡| 老司机亚洲免费影院| 亚洲五月色婷婷综合| 国产精品蜜桃在线观看| 婷婷成人精品国产| 午夜激情av网站| 成年人午夜在线观看视频| 一二三四在线观看免费中文在| 日韩 亚洲 欧美在线| 久久久久国产精品人妻一区二区| 最近2019中文字幕mv第一页| 久久久久久久精品精品| av在线播放精品| av国产久精品久网站免费入址| 黄色视频在线播放观看不卡| 国产无遮挡羞羞视频在线观看| 国产一级毛片在线| 女的被弄到高潮叫床怎么办| 婷婷成人精品国产| 欧美亚洲日本最大视频资源| 日日啪夜夜爽| 亚洲第一av免费看| 国产不卡av网站在线观看| 美女高潮到喷水免费观看| 久久久久精品人妻al黑| 国产成人精品福利久久| 色94色欧美一区二区| 欧美日本中文国产一区发布| 国产av精品麻豆| 99re6热这里在线精品视频| 久久婷婷青草| 精品一品国产午夜福利视频| 秋霞伦理黄片| 国产不卡av网站在线观看| 九草在线视频观看| 国产av精品麻豆| 伊人亚洲综合成人网| 国产精品成人在线| 亚洲成国产人片在线观看| 男女国产视频网站| 欧美日韩视频高清一区二区三区二| 久久青草综合色| 一区二区三区乱码不卡18| 国产在视频线精品| 97人妻天天添夜夜摸| 2018国产大陆天天弄谢| 久久久精品94久久精品| 一级毛片 在线播放| 午夜久久久在线观看| 亚洲精品久久午夜乱码| 观看美女的网站| 校园人妻丝袜中文字幕| 少妇的丰满在线观看| 国产精品无大码| 久久免费观看电影| 欧美日韩视频精品一区| 天堂中文最新版在线下载| 久久国产精品大桥未久av| 亚洲欧美精品自产自拍| 女人久久www免费人成看片| 国产高清不卡午夜福利| 夫妻性生交免费视频一级片| 麻豆av在线久日| 少妇人妻久久综合中文| 97人妻天天添夜夜摸| 一级爰片在线观看| 欧美黑人欧美精品刺激| 日韩制服丝袜自拍偷拍| 国产精品二区激情视频| 国产精品一国产av| svipshipincom国产片| 青青草视频在线视频观看| 男人爽女人下面视频在线观看| 啦啦啦在线免费观看视频4| 侵犯人妻中文字幕一二三四区| 精品视频人人做人人爽| 亚洲国产欧美一区二区综合| 一边亲一边摸免费视频| av女优亚洲男人天堂| h视频一区二区三区| 永久免费av网站大全| 男女国产视频网站| 亚洲精品第二区| 国产片特级美女逼逼视频| 晚上一个人看的免费电影| 久久久久精品人妻al黑| 国产熟女午夜一区二区三区| 国产精品 欧美亚洲| 亚洲精品第二区| 1024香蕉在线观看| 亚洲激情五月婷婷啪啪| 国产精品 欧美亚洲| 国产日韩欧美视频二区| 国产成人一区二区在线| 久久性视频一级片| 久久人人爽av亚洲精品天堂| 日韩 亚洲 欧美在线| 考比视频在线观看| 日韩av在线免费看完整版不卡| 国产精品免费大片| 久久97久久精品| 18禁动态无遮挡网站| 国产精品二区激情视频| 亚洲av中文av极速乱| 中文精品一卡2卡3卡4更新| 国产精品三级大全| 久久韩国三级中文字幕| 欧美日韩av久久| 99久久人妻综合| 国产精品香港三级国产av潘金莲 | 国产成人午夜福利电影在线观看| 久久精品人人爽人人爽视色| 欧美亚洲 丝袜 人妻 在线| √禁漫天堂资源中文www| av女优亚洲男人天堂| 日韩精品免费视频一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 日本av免费视频播放| 成人国产麻豆网| www.精华液| 国产在视频线精品| 18禁动态无遮挡网站| 国产精品.久久久| 午夜福利网站1000一区二区三区| 精品午夜福利在线看| 王馨瑶露胸无遮挡在线观看| 韩国av在线不卡| 在线观看www视频免费| 极品少妇高潮喷水抽搐| 深夜精品福利| 尾随美女入室| 伊人久久国产一区二区| 国产成人91sexporn| 欧美精品一区二区免费开放| 亚洲精品美女久久av网站| 波野结衣二区三区在线| 久久精品人人爽人人爽视色| 日韩熟女老妇一区二区性免费视频| a级片在线免费高清观看视频| 亚洲国产日韩一区二区| 赤兔流量卡办理| 1024香蕉在线观看| 九草在线视频观看| 丝袜在线中文字幕| 中文字幕亚洲精品专区| 久久久久久人妻| 亚洲人成网站在线观看播放| 久久鲁丝午夜福利片| 人成视频在线观看免费观看| 欧美精品一区二区大全| 亚洲人成77777在线视频| 亚洲,欧美,日韩| 99精国产麻豆久久婷婷| 精品酒店卫生间| 亚洲精品久久久久久婷婷小说| 久久99一区二区三区| 只有这里有精品99| 国产欧美日韩综合在线一区二区| 天天操日日干夜夜撸| 中文天堂在线官网| 亚洲色图 男人天堂 中文字幕| e午夜精品久久久久久久| 一区二区av电影网| av片东京热男人的天堂| 少妇的丰满在线观看| www.精华液| 操美女的视频在线观看| 亚洲国产av影院在线观看| 香蕉国产在线看| 1024视频免费在线观看| 国产有黄有色有爽视频| 亚洲美女黄色视频免费看| 亚洲精品aⅴ在线观看| 国产精品成人在线| 国产黄色免费在线视频| 日本午夜av视频| 桃花免费在线播放| 亚洲精品视频女| 国产精品偷伦视频观看了| 69精品国产乱码久久久| av网站在线播放免费| 黄片播放在线免费| 一区二区三区四区激情视频| 午夜福利影视在线免费观看| 中文字幕制服av| 久久鲁丝午夜福利片| 成年美女黄网站色视频大全免费| 日本wwww免费看| 亚洲一级一片aⅴ在线观看| 国产免费福利视频在线观看| 欧美黑人欧美精品刺激| 亚洲av日韩在线播放| 久久久久精品久久久久真实原创| 国产日韩一区二区三区精品不卡| 午夜激情av网站| 精品人妻一区二区三区麻豆| 亚洲精品久久成人aⅴ小说| 久久午夜综合久久蜜桃| 国产极品天堂在线| 99热网站在线观看| av在线播放精品| 午夜日韩欧美国产| 另类精品久久| 另类精品久久| 天天躁夜夜躁狠狠躁躁| 赤兔流量卡办理| 欧美日韩一区二区视频在线观看视频在线| 男女之事视频高清在线观看 | 国产欧美亚洲国产| 国产黄色免费在线视频| 91aial.com中文字幕在线观看| av免费观看日本| 亚洲伊人久久精品综合| 亚洲欧洲日产国产| 99久久99久久久精品蜜桃| 如何舔出高潮| 国产免费福利视频在线观看| 国产日韩一区二区三区精品不卡| 亚洲图色成人| 中文字幕另类日韩欧美亚洲嫩草| 国产精品亚洲av一区麻豆 | 黄色怎么调成土黄色| 91成人精品电影| 欧美精品一区二区大全| 国产1区2区3区精品| 老司机靠b影院| 国产成人精品久久二区二区91 | av免费观看日本| 久久精品熟女亚洲av麻豆精品| 亚洲成人国产一区在线观看 | 成年av动漫网址| 午夜激情av网站| 卡戴珊不雅视频在线播放| 999精品在线视频| 久久久精品免费免费高清| 欧美人与善性xxx| 国产精品 国内视频| 欧美人与性动交α欧美精品济南到| 亚洲成国产人片在线观看| 9191精品国产免费久久| 成年人午夜在线观看视频| 亚洲精品国产av蜜桃| a级片在线免费高清观看视频| 亚洲av日韩精品久久久久久密 | 亚洲精品国产色婷婷电影| 2021少妇久久久久久久久久久| 一二三四在线观看免费中文在| 超碰成人久久| 国产在视频线精品| 国产av国产精品国产| 亚洲精品国产区一区二| 韩国精品一区二区三区| 午夜免费男女啪啪视频观看| 精品久久久久久电影网| 97精品久久久久久久久久精品| 男女午夜视频在线观看| 午夜福利网站1000一区二区三区| 久久99热这里只频精品6学生| 97在线人人人人妻| 99精品久久久久人妻精品| 亚洲精品国产区一区二| 老司机影院毛片| 波野结衣二区三区在线| 国产亚洲欧美精品永久| 黄色视频在线播放观看不卡| kizo精华| 各种免费的搞黄视频| 精品亚洲成国产av| 交换朋友夫妻互换小说| 咕卡用的链子| 国产午夜精品一二区理论片| 丰满乱子伦码专区| 少妇 在线观看| 波野结衣二区三区在线| 欧美国产精品一级二级三级| 欧美日韩综合久久久久久| 69精品国产乱码久久久| 亚洲欧美一区二区三区国产| 岛国毛片在线播放| 9色porny在线观看| 一本一本久久a久久精品综合妖精| 精品国产乱码久久久久久小说| 黄网站色视频无遮挡免费观看| 一级毛片黄色毛片免费观看视频| 日本av手机在线免费观看| 亚洲精品中文字幕在线视频| 天天添夜夜摸| 观看美女的网站| 亚洲av综合色区一区| 亚洲精品日韩在线中文字幕| 91精品伊人久久大香线蕉| 男人添女人高潮全过程视频| av在线观看视频网站免费| 国产日韩一区二区三区精品不卡| 99re6热这里在线精品视频| 大码成人一级视频| 国产精品一二三区在线看| 国产黄频视频在线观看| 欧美日韩国产mv在线观看视频| 亚洲成av片中文字幕在线观看| 一本色道久久久久久精品综合| 视频在线观看一区二区三区| 国产99久久九九免费精品| 免费看不卡的av| 亚洲美女黄色视频免费看| 丰满饥渴人妻一区二区三| 久久精品人人爽人人爽视色| 久久热在线av| 日韩欧美一区视频在线观看| 日本vs欧美在线观看视频| 久久久久久久国产电影| 美女高潮到喷水免费观看| 国产精品一区二区在线不卡| 亚洲国产最新在线播放| 成年动漫av网址| 亚洲av男天堂| 亚洲精品中文字幕在线视频| 超色免费av| 国产精品久久久久久精品古装| 天天操日日干夜夜撸| 天天躁日日躁夜夜躁夜夜| 肉色欧美久久久久久久蜜桃| 男男h啪啪无遮挡| 国产男人的电影天堂91| 老汉色∧v一级毛片| 欧美xxⅹ黑人| 亚洲国产欧美在线一区| 亚洲专区中文字幕在线 | 宅男免费午夜| 中文天堂在线官网| 亚洲成人手机| 久久精品熟女亚洲av麻豆精品| 欧美日韩视频精品一区| av在线播放精品| 宅男免费午夜| 国产激情久久老熟女| 日韩电影二区| 中文字幕最新亚洲高清| 一区二区三区精品91| 国产成人免费无遮挡视频| 久久天堂一区二区三区四区| 一级片'在线观看视频| 亚洲国产欧美日韩在线播放| 日韩伦理黄色片| 少妇人妻精品综合一区二区| 91精品三级在线观看| av一本久久久久| 深夜精品福利| 亚洲av福利一区| 老司机亚洲免费影院| 一区二区日韩欧美中文字幕| 久久99一区二区三区| 欧美国产精品一级二级三级| 国产精品久久久久成人av| 免费黄网站久久成人精品| 女人久久www免费人成看片| e午夜精品久久久久久久| 黄色视频在线播放观看不卡| 久久毛片免费看一区二区三区| 在线观看免费日韩欧美大片| 精品一区二区三区av网在线观看 | 电影成人av| 久久久久精品国产欧美久久久 | av网站在线播放免费| 狂野欧美激情性bbbbbb| 久久99一区二区三区| 黄片播放在线免费| 亚洲婷婷狠狠爱综合网| 人妻 亚洲 视频| 18在线观看网站| av片东京热男人的天堂| 成人影院久久| 三上悠亚av全集在线观看| 亚洲成色77777| 少妇人妻 视频| 婷婷色综合大香蕉| 亚洲一级一片aⅴ在线观看| av电影中文网址| 欧美亚洲日本最大视频资源| 欧美黄色片欧美黄色片| 国产精品无大码| av网站在线播放免费| 色视频在线一区二区三区| 在线观看一区二区三区激情| 国产黄色免费在线视频| 高清不卡的av网站| 一个人免费看片子| 欧美日韩成人在线一区二区| 王馨瑶露胸无遮挡在线观看| 亚洲欧美色中文字幕在线| 男女下面插进去视频免费观看| www.熟女人妻精品国产| 美女福利国产在线| 嫩草影视91久久| 欧美日本中文国产一区发布| 性少妇av在线| 精品少妇一区二区三区视频日本电影 | 免费观看av网站的网址| 熟女少妇亚洲综合色aaa.| 色婷婷久久久亚洲欧美| 老司机影院成人| 久久综合国产亚洲精品| 欧美 日韩 精品 国产| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美精品综合一区二区三区| 日本午夜av视频| 亚洲av综合色区一区| 制服诱惑二区| 男人舔女人的私密视频| 操出白浆在线播放| 国产黄色视频一区二区在线观看| 日韩熟女老妇一区二区性免费视频| 久久久久久久久免费视频了| 嫩草影视91久久| 大话2 男鬼变身卡| 亚洲综合色网址| 侵犯人妻中文字幕一二三四区| 中文字幕人妻丝袜一区二区 | 久久ye,这里只有精品| 伊人久久大香线蕉亚洲五| 亚洲第一区二区三区不卡| 国产精品麻豆人妻色哟哟久久| 操出白浆在线播放| 日韩 欧美 亚洲 中文字幕| 最新在线观看一区二区三区 | 美女中出高潮动态图| 国产日韩欧美亚洲二区| 国产熟女午夜一区二区三区| 一级黄片播放器| 另类精品久久| 国产激情久久老熟女| 亚洲国产最新在线播放| 国产午夜精品一二区理论片| 亚洲国产av新网站| 国产精品免费视频内射| 国产av精品麻豆| 制服人妻中文乱码| 久久人人97超碰香蕉20202| 一本一本久久a久久精品综合妖精| 操出白浆在线播放| 婷婷色综合www| 日韩 欧美 亚洲 中文字幕| 一级,二级,三级黄色视频| 啦啦啦中文免费视频观看日本| 韩国精品一区二区三区| 精品国产乱码久久久久久男人| 黄片无遮挡物在线观看| 免费女性裸体啪啪无遮挡网站| 日韩中文字幕视频在线看片| 下体分泌物呈黄色| 午夜免费男女啪啪视频观看| 久久这里只有精品19| 亚洲av成人不卡在线观看播放网 | 久久精品国产综合久久久| 午夜福利,免费看| 中国三级夫妇交换| 热99国产精品久久久久久7| 国产免费又黄又爽又色| 日本av免费视频播放| 在线 av 中文字幕| 日本黄色日本黄色录像| 少妇人妻 视频| 一本久久精品| 成人免费观看视频高清| 国产人伦9x9x在线观看| 久久毛片免费看一区二区三区| 久久这里只有精品19| 嫩草影视91久久| 汤姆久久久久久久影院中文字幕| 成人漫画全彩无遮挡| 久热这里只有精品99| 色婷婷av一区二区三区视频| 亚洲欧美中文字幕日韩二区| 性高湖久久久久久久久免费观看| 精品一区在线观看国产| 无遮挡黄片免费观看| 国产国语露脸激情在线看| 水蜜桃什么品种好| 欧美xxⅹ黑人| 在线 av 中文字幕| 久久精品人人爽人人爽视色| 久久这里只有精品19| 99久久人妻综合| 街头女战士在线观看网站| av国产久精品久网站免费入址| 精品少妇黑人巨大在线播放| 最近最新中文字幕免费大全7| 久久久久久久大尺度免费视频| 亚洲熟女精品中文字幕| 国产老妇伦熟女老妇高清| 亚洲在久久综合| 欧美日韩视频精品一区| av国产久精品久网站免费入址| 国产亚洲精品第一综合不卡| 日韩大码丰满熟妇| 国产免费一区二区三区四区乱码| 国产欧美亚洲国产| 校园人妻丝袜中文字幕| 午夜免费观看性视频| 别揉我奶头~嗯~啊~动态视频 | 日日爽夜夜爽网站| 午夜福利一区二区在线看| 桃花免费在线播放| 欧美亚洲日本最大视频资源| 日韩av不卡免费在线播放| 国语对白做爰xxxⅹ性视频网站| 不卡视频在线观看欧美| 大码成人一级视频| 热99久久久久精品小说推荐| 欧美日韩一区二区视频在线观看视频在线| 精品少妇一区二区三区视频日本电影 | 亚洲熟女精品中文字幕| 飞空精品影院首页| www.熟女人妻精品国产| 精品国产露脸久久av麻豆| 人妻人人澡人人爽人人| 日韩,欧美,国产一区二区三区| 观看美女的网站| 国产xxxxx性猛交| 国产97色在线日韩免费| 老熟女久久久| 别揉我奶头~嗯~啊~动态视频 | av福利片在线| 80岁老熟妇乱子伦牲交| 人人澡人人妻人| 久久久久久久久久久久大奶| 亚洲欧美一区二区三区国产| 亚洲一级一片aⅴ在线观看| 少妇被粗大的猛进出69影院| 咕卡用的链子| 老鸭窝网址在线观看| 97在线人人人人妻| 亚洲欧洲日产国产| 天天添夜夜摸|