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

    利用COSMIC低軌衛(wèi)星對GPS信號的頂部TEC觀測資料研究等離子體層電子含量的變化特征

    2016-07-29 07:56:59張滿蓮劉立波萬衛(wèi)星寧百齊
    地球物理學報 2016年1期

    張滿蓮, 劉立波, 萬衛(wèi)星, 寧百齊

    中國科學院地質(zhì)與地球物理研究所,中國科學院地球與行星物理重點實驗室, 北京 100029

    ?

    利用COSMIC低軌衛(wèi)星對GPS信號的頂部TEC觀測資料研究等離子體層電子含量的變化特征

    張滿蓮, 劉立波, 萬衛(wèi)星, 寧百齊

    中國科學院地質(zhì)與地球物理研究所,中國科學院地球與行星物理重點實驗室, 北京100029

    摘要本文嘗試利用COSMIC低軌衛(wèi)星對GPS信號的頂部TEC觀測資料研究等離子體層電子含量(簡稱PEC)的變化規(guī)律.首先介紹從低軌衛(wèi)星對GPS的頂部TEC觀測資料提取等離子體層垂直電子含量的方法,然后利用該方法提取2008年全年的PEC數(shù)據(jù),進而研究了2008年這一太陽活動低年PEC隨地磁緯度(MLAT)、磁地方時(MLT)以及不同季節(jié)的變化規(guī)律.此外,還利用提取的120°E和300°E經(jīng)度鏈上的數(shù)據(jù)對比研究了PEC的經(jīng)度變化情形.研究結(jié)果表明:(1)PEC主要集中分布在磁赤道±45°之間的一個繞地球的環(huán)帶狀區(qū)域中;(2)PEC表現(xiàn)出以下的晝夜變化規(guī)律特征:白天時段之值高于夜間,約在12—16MLT之間達到最高峰值,而最小PEC值出現(xiàn)在日出前大約4—5MLT左右的時段;(3)相比其他季節(jié)月份而言,PEC在北半球夏季月份(5—8月)具有最小值;(4)PEC存在明顯的經(jīng)度變化,不同經(jīng)度鏈上的PEC存在不同的季節(jié)變化特征.

    關鍵詞等離子體層; 電離層; 總電子含量(TEC); 等離子體層電子含量

    1引言

    電波信號受到電離成分的作用會發(fā)生折射而產(chǎn)生延遲效應,從而對衛(wèi)星定位、導航授時以及遙感、遙測等產(chǎn)生額外的誤差.這種誤差與電波傳播路徑上單位截面柱體里的積分電子含量(通常稱為總電子含量TEC)有關(Davíes, 1990).全球?qū)Ш蕉ㄎ幌到y(tǒng)(GPS)的電波信號穿過地球大氣上空的電離層和等離子體層.盡管等離子體層電子密度比電離層電子密度要稀薄很多,但是由于其覆蓋的高度范圍很大,等離子體層的電子含量仍然占有相當?shù)姆至?特別是在夜間和太陽活動平靜時期,等離子體層電子含量相對總電子含量TEC而言占有相當大的比例,白天一般約為10%左右,而夜間則可高達60%(Yizengaw et al., 2008).因此,不僅電離層對信號傳播有影響,等離子體層對信號傳播的影響也不可忽視.在許多實際應用中(例如用GPS單頻接收機進行定位時),常常要用到電離層修正模型來消除由電離成分對電波的延遲效應而產(chǎn)生的定位誤差.但是目前一般的電離層修正模型,例如國際參考電離層(IRI)模型 (Bilitza, 2001; Bilitza and Reinisch, 2008)、NeQuick模型(Leitinger et al., 2005; Nava et al., 2008),通常只包括電離層的電子含量,而沒有將等離子體層的含量考慮在內(nèi).近年來國際參考電離層工作組開始倡導與鼓勵人們開展等離子體層電子含量模型的研究,以期將IRI模型擴展到等離子體層范圍(Bilitza et al., 2011).因此,通過對等離子體層電子含量變化特征的研究,可為建立或完善等離子體層電子含量模型奠定基礎,為將國際參考電離層模型擴展到等離子體層做出貢獻,為導航定位等實際應用中的誤差修正提供一定的參考.另一方面,電離層與等離子體層有很強的耦合關系,通過對等離子體層電子含量變化規(guī)律的研究,能夠加深理解電離層與等離子體層之間所發(fā)生的各種動力學耦合過程的作用機制.

    前人對等離子體層電子含量的提取技術(shù),一般可歸納為以下幾種方法:(1)由兩顆不同高度衛(wèi)星(如GPS與另一低高度衛(wèi)星)觀測的TEC數(shù)據(jù),由兩者之差得到等離子體層電子含量 (Lunt et al., 1999; Chen and Yao, 2015).這種方法的缺點是兩顆衛(wèi)星所在的高度范圍不一定剛好合適于等離子體層的高度范圍.(2)利用地面垂直測高儀的觀測數(shù)據(jù)來估算電離層電子含量ITEC,結(jié)合對GPS觀測的TEC,由兩者之差得到等離子體層電子含量 (Belehaki et al., 2004; Zhang et al., 2006).這種方法中,由于只有F2層峰高以下的底部電離層電子密度數(shù)據(jù)是真實的觀測結(jié)果,而F2層峰高以上的頂部電離層電子密度是通過假設F2層峰高以上的電子密度符合Chapman剖面分布而外推估算得到 (Reinisch and Huang, 2001).因而由垂直測高儀觀測資料估算得到的電離層ITEC可能會有較大的誤差,從而導致最后求得的等離子體層電子含量誤差較大.(3)利用非相干散射雷達(ISR)觀測資料與地面GPS接收機觀測資料相結(jié)合 (種小燕等,2013).非相干散射雷達能夠探測包括F2層峰高以下底部電離層及峰高以上的頂部電離層高度的電子密度.利用這一設備所觀測得到的資料來推算電離層電子含量將比基于電離層垂測儀觀測資料進行的推算更具真實性,從而可得到精度較高的等離子體層電子含量.然而,非相干散射雷達全球覆蓋率極為有限,且觀測時間不連續(xù),這為研究等離子體層電子含量的全球分布特征與建模帶來局限.此外,上述這些方法還由于涉及到不同的觀測設備或觀測手段而帶來的觀測誤差導致所得到的結(jié)果有較大的不確定性.COSMIC低軌衛(wèi)星(LEO)的精確定軌(POD)接收機對GPS信號的頂部測量為我們提供了大量的等離子體層電子含量信息的數(shù)據(jù).COSMIC低軌衛(wèi)星對GPS信號的這種頂部TEC觀測資料具有良好的全球覆蓋性與觀測時間的連續(xù)性,因此特別適合用于研究等離子體層電子含量的全球分布特征與不同時間尺度的變化性,從而為等離子體層電子含量的全球建模提供數(shù)據(jù)基礎.本文將利用COSMIC低軌衛(wèi)星對GPS衛(wèi)星信號的頂部TEC觀測資料來提取等離子體層電子含量,主要技術(shù)涉及由LEO-GPS衛(wèi)星連線傳播路徑上的總電子含量演算得到等離子體層垂直電子含量(PEC,Plasmaspheric Electron Content).這種方法的優(yōu)點在于涉及的觀測資料單一一致,所帶來的觀測及反演所積累的誤差小,因而能夠得到精度較高的等離子體層電子含量.然后將利用提取的等離子體層電子含量數(shù)據(jù)研究等離子體層電子含量隨地磁緯度(MLAT)、磁地方時(MLT)、季節(jié)的變化規(guī)律,同時也抽取兩個經(jīng)度鏈上(120°E和300°E)的數(shù)據(jù)對比研究PEC的經(jīng)度變化性.本文第2節(jié)將簡要介紹原始數(shù)據(jù)來源及提取PEC數(shù)據(jù)采用的方法,第3節(jié)給出研究結(jié)果及對結(jié)果的解釋,第4節(jié)為結(jié)論.

    2原始數(shù)據(jù)來源及提取PEC數(shù)據(jù)的方法簡介

    本文所用原始數(shù)據(jù)來源于COSMIC/FORMSAT-3(The Constellation Observing System for Meteorology, Ionosphere and Climate/Formosa Satellite 3) 星座的低軌衛(wèi)星(LEO)上的精確定軌(POD)接收機對GPS信號的頂部TEC測量數(shù)據(jù).COSMIC/FORMOSAT-3是中國臺灣和美國于2006年4月15日聯(lián)合發(fā)射的由6顆低軌(LEO)微衛(wèi)星組成的一個星座.每顆微衛(wèi)星攜帶一個GPS掩星實驗(GOX)載荷,主要設計用于對地球大氣和電離層進行無線電掩星測量.衛(wèi)星軌道周期約100 min,傾角72°.COSMIC 觀測數(shù)據(jù)由以下網(wǎng)址提供:(1)FS-3/COSMIC臺灣分析中心(TACC, http:∥tacc.cwb.gov.tw/en/);(2)COSMIC 數(shù)據(jù)分析和存檔中心(CDAAC, http:∥www.cosmic.ucar.edu/cdaac/).本文用到的數(shù)據(jù)主要是COSMIC星座LEO衛(wèi)星上用于精確定軌(POD)的衛(wèi)星天線對GPS衛(wèi)星接收信號而計算得到的LEO-GPS沿線的TEC數(shù)據(jù),即COSMIC產(chǎn)品數(shù)據(jù)目錄level1b下的podTEC數(shù)據(jù).該數(shù)據(jù)產(chǎn)品已經(jīng)是LEO-GPS沿線的絕對TEC值,已經(jīng)去除了接收機本身的硬件誤差 (Yue et al., 2011).本文用到的是2008年的全年數(shù)據(jù).

    COSMIC星座LEO衛(wèi)星飛行的高度大約在800 km高度左右,在太陽活動低年,這一高度可近似看成是電離層與等離子體層的過渡高度,因此LEO衛(wèi)星高度到GPS衛(wèi)星高度之間的垂直高度積分電子含量可近似看成等離子體層電子含量PEC.由于podTEC數(shù)據(jù)給出的是LEO-GPS沿線路徑上的斜路徑積分電子含量,為了求取垂直方向(即天頂方向)的PEC,必須選取適當?shù)霓D(zhuǎn)換因子f(ε)把斜路徑方向的podTEC轉(zhuǎn)換成垂直方向的PEC:

    (1)

    本文采用Foelsche和Kirchengast(2002)以及Yue等(2011)所介紹和引用的轉(zhuǎn)換因子f(ε),其表達式為:

    (2)

    (3)

    (4)

    其中ε為LEO-GPS衛(wèi)星連線方向的抬升角,Re為地球半徑,Hppt是LEO-GPS衛(wèi)星之間電波射線在等離子體層中的穿刺點高度,Horb是低軌衛(wèi)星所在的高度.Hppt可選為幾百到幾千公里.Yue等(2011)選取Hppt值為3000 km.我們試過不同的Hppt值,雖然用不同的Hppt值計算得到的PEC值會有些微差別,但并不影響本文研究所得到的主要結(jié)論.因此,本文按Yue等(2011)選取3000 km 為Hppt值.圖1是2008年1月27日全天當中LEO-GPS衛(wèi)星連線電波信號在Hppt=3000 km處的穿刺點坐標分布.從圖中可見衛(wèi)星信號穿刺點非常密集地覆蓋了全球所有空間,因而這一數(shù)據(jù)源很適宜于研究等離子體層電子含量的全球分布特征.

    3等離子體層電子含量(PEC)變化特征

    3.1等離子體層電子含量(PEC)隨地磁緯度的變化和磁地方時的晝夜變化

    為了研究PEC隨地磁緯度(MLAT)和磁地方時(MLT)的變化規(guī)律,我們選取太陽活動低年2008年的全年數(shù)據(jù),采用前述轉(zhuǎn)換方法演算出天頂方向的PEC,然后對全年數(shù)據(jù)求取平均,得到如圖2所示的2008年PEC隨MLAT和MLT的平均分布狀況.從圖2中可看出,PEC主要集中分布在磁赤道±45°之間的一個繞地球的環(huán)帶狀區(qū)域中.這與通常的理論預期和前人利用其他手段的觀測結(jié)果是一致的,由此表明我們采用的提取等離子體層電子含量的方法是可靠的.在晝夜變化規(guī)律方面,從圖2中可看出PEC呈現(xiàn)出以下的變化特征:白天時段之值高于夜間時段,PEC約在12—16MLT之間的某個時刻達到最高峰值(7~8 TECu左右),而最小PEC值大約出現(xiàn)在日出前(約4—5MLT左右)的時段.等離子體層電子含量的這一晝夜變化特征與等離子體層和電離層之間存在強烈的電動力學耦合作用有關.我們知道,等離子體層的離子來源與白天電離層離子沿磁力線的上行流動有關.等離子體層是電離層離子的一個“蓄水池”(Singh et al., 2011).白天電離層離子會沿著磁力線向上漂移與擴散填充到等離子體層中,而晚上等離子體層中的離子會返回到低高度范圍維持夜間的F區(qū)電離層.這一耦合作用過程產(chǎn)生的效果便導致了我們所觀測到的等離子體層電子含量PEC的上述晝夜變化特征.

    圖1 2008年1月27日一天內(nèi)COSMIC低軌衛(wèi)星與GPS衛(wèi)星之間連線在3000 km高度處穿刺點的地理分布情況

    圖2 2008年等離子體層電子含量PEC隨磁地方時和地磁緯度的平均分布狀況

    3.2等離子體層電子含量(PEC)的季節(jié)變化

    為研究PEC隨不同季節(jié)的變化規(guī)律,我們將2008年的PEC數(shù)據(jù)按不同季節(jié)劃分,求取了四個不同季節(jié)的PEC的平均變化情形,所得結(jié)果如圖3所示.圖中左上、右上、左下及右下分別對應的是2008年北半球春(3、4月份)、夏(5、6、7、8月份)、秋(9、10月份)、冬(1、2、11、12月份)四個季節(jié)的PEC的季節(jié)平均.從圖3中可看出,四個不同的季節(jié)PEC隨MLT及MLAT的變化形態(tài)與3.1節(jié)中的結(jié)果完全一致:等離子體層電子含量主要集中分布在磁赤道±45°之間的一個繞地球的環(huán)帶狀區(qū)域中;PEC白天之值高于夜間,在約12—16MLT之間達到最高峰值,而在日出前的時段(約4—5MLT左右)PEC達到最小值.然而,PEC也表現(xiàn)出以下的季節(jié)變化特征:北半球夏季PEC白天的峰值低于其他所有季節(jié)的白天峰值,而在北半球春季PEC白天的峰值在四個季節(jié)中則是最高的,秋、冬季之值則介于這二者之間,且冬季之值又高于秋季.Chen和Yao (2015)利用GPS-TEC和COSMIC LEO衛(wèi)星掩星觀測資料研究2011年3、6、9、12月份PEC的季節(jié)變化.他們所研究的2011年處于太陽活動的上升年,本文所用數(shù)據(jù)采用的是處在太陽活動低年的2008年,因此很顯然其得到的PEC數(shù)值要明顯高于本文得到的值(請參見Chen and Yao (2015) 中的圖3).但其顯示的季節(jié)變化特征與我們上述描述的季節(jié)變化特征基本一致.這些變化特征基本上與電離層的總電子含量ITEC的變化特征表現(xiàn)是一致的.考慮前面3.1節(jié)中所述的原因(電離層與等離子體層之間存在很強的耦合關系),PEC的這些季節(jié)變化特征都能得到合理的解釋.

    3.3等離子體層電子含量(PEC)的經(jīng)度變化

    總而言之,將物化探技術(shù)應用到礦產(chǎn)勘查工作以及地質(zhì)效果考察之中,相關技術(shù)工作者需要合理、科學的應用該技術(shù),這樣才能有效勘查礦產(chǎn)資源,有效提高礦產(chǎn)勘查工作的效率。

    為研究PEC是否存在經(jīng)度變化,我們特別提取120°E和300°E經(jīng)度鏈上的數(shù)據(jù),分別計算了2008年北半球春、夏、秋、冬四個季節(jié)中這兩個經(jīng)度鏈上的PEC的變化情形,結(jié)果如圖4—5中所示.其中圖4顯示的是120°E經(jīng)度鏈上的PEC變化情形,圖5顯示的是300°E經(jīng)度鏈上的PEC變化情形.從圖4—5中可看出,不同經(jīng)度鏈上的PEC具有明顯的不同的變化規(guī)律.對于北半球春、夏、秋這三個季節(jié)120°E經(jīng)度鏈上的PEC的白天峰值明顯高于300°E經(jīng)度鏈上的峰值.而對于北半球冬季則是300°E經(jīng)度鏈上的值大于120°E經(jīng)度鏈上的值.此外,這兩個經(jīng)度鏈上的PEC還表現(xiàn)出不同的季節(jié)變化規(guī)律:120°E經(jīng)度鏈上的PEC,白天峰值在春秋季稍高于冬夏季,亦即表現(xiàn)為以半年變化為主的特征;而300°E經(jīng)度鏈上的PEC則表現(xiàn)出非常明顯的以年變化為主的季節(jié)變化特征:其PEC白天峰值在北半球冬季最大、夏季最小,而春、秋季之值介于前兩者之間.特別是比較北半球冬季和夏季的情形時,會發(fā)現(xiàn)300°E經(jīng)度鏈上的PEC值差別非常大,形成很鮮明的對比.而120°E經(jīng)度鏈上的PEC值在這兩個季節(jié)的差別不是那么明顯.對于PEC在這兩個經(jīng)度鏈上表現(xiàn)出的不同季節(jié)變化特征,我們給出如下的解釋:電離層與等離子體層之間存在強烈耦合,等離子體層的離子來源與白天電離層離子沿磁力線的上行流動有關.因此PEC這種不同經(jīng)度上的不同季節(jié)變化特征應該與等離子體層下面的電離層有關.電離層中最關鍵的參量之一是F2層峰值電子濃度NmF2(或臨界頻率foF2).我們知道,NmF2(電離層的總電子含量ITEC也是)在全球不同區(qū)域會表現(xiàn)出各種不同的季節(jié)變化上的異?,F(xiàn)象,包括(1)半年異常 (semiannual anomaly):NmF2在春秋季月份達到最大值;(2)冬季異常(又叫季節(jié)異常, seasonal anomaly):NmF2在冬季月份達到最大值;(3)年度異常(也叫非季節(jié)異常, non-seasonal anomaly):全球平均而言NmF2在12月份之值高于6月份之值.其中年度異常(3)可理解為北半球的冬季異常強于南半球的結(jié)果.Rishbeth (1998)對電離層這些季節(jié)變化異?,F(xiàn)象給出非常詳盡的解釋(有興趣的讀者可仔細參閱該文章),其主要原因可歸納為:(1)太陽天頂角隨季節(jié)的變化導致太陽輻射對中性大氣的電離產(chǎn)生率的變化;(2)熱層大氣環(huán)流引起的背景中性大氣層氧氮濃度比[O/N2]的季節(jié)變化導致離子電離率和復合率的變化;(3)熱層大氣環(huán)流引起的風場對電離層的抬升與降低作用導致離子電離率和復合率的變化,該作用的大小與地磁場位形在全球的分布有關.這三個因素在不同的經(jīng)度、緯度其影響程度不同,就導致全球不同區(qū)域的電離層表現(xiàn)出不同的季節(jié)變化特征.Torr M R和Torr D G(1973)以及Rishbeth(1998)給出了不同太陽活動水平下全球不同區(qū)域foF2的主導季節(jié)變化特征.從Rishbeth(1998)文章中的圖5可以看出,在太陽活動低年,120°E經(jīng)度鏈附近的中、低緯度上電離層是以半年變化為主導的;而300°E經(jīng)度鏈附近的中低緯度,北半球以冬季異常(12月份至點季最大,6月份至點季最小)為主導特征,南半球則以正常的夏季(12月份至點季)最大、冬季(6月份至點季)最小為主導特征.電離層在這兩個經(jīng)度鏈扇區(qū)的季節(jié)變化特征,通過與等離子體層的動力學/電動力學耦合反映到等離子體層中,從而使PEC表現(xiàn)出相似的季節(jié)變化特征:120°E經(jīng)度鏈上的PEC表現(xiàn)為在兩分點達最大值的半年變化特征,而300°E經(jīng)度鏈上的PEC表現(xiàn)為在12月份至點季達最大值、在6月份至點季達最小值的年變化特征.

    圖3 2008年等離子體層電子含量PEC的季節(jié)變化,其中左上、右上、左下、右下圖分別對應的是北半球春(3—4月份)、夏(5—8月份)、秋(9—10月份)、冬(1—2、11—12月份)四個季節(jié)的PEC變化情形

    圖4 120°E經(jīng)度鏈上等離子體層電子含量PEC在2008年北半球春(3 —4月份)、夏(5 —8月份)、秋(9 —10月份)、冬(1—2、11—12月份)四個季節(jié)里的變化情形

    圖5 300°E經(jīng)度鏈上等離子體層電子含量PEC在2008年北半球春(3—4月份)、夏(5—8月份)、秋(9—10月份)、冬(1—2、11—12月份)四個季節(jié)里的變化情形

    4結(jié)論

    本文初步嘗試利用COSMIC低軌衛(wèi)星對GPS信號的頂部TEC觀測資料產(chǎn)品數(shù)據(jù)podTEC提取等離子體層電子含量PEC.并利用所提取的PEC數(shù)據(jù)研究了2008年太陽活動低年PEC隨磁地方時(MLT)、地磁緯度(MLAT)和不同季節(jié)的變化規(guī)律.同時對比研究了120°E和300°E經(jīng)度鏈上PEC的不同變化特征.總結(jié)本文前面的結(jié)果,可以得出以下結(jié)論:(1)等離子體層電子含量主要集中分布在磁赤道±45°之間的一個繞地球的環(huán)帶狀區(qū)域中;(2)PEC表現(xiàn)的晝夜變化規(guī)律為白天之值高于夜間,約在12—16MLT之間達到最高峰值(量級大約為7~8TECu左右),而最小PEC值大約出現(xiàn)在日出前(約4—5MLT)的時段;(3)相比其他季節(jié)月份而言,等離子體層電子含量PEC在北半球夏季月份(即5—8月)具有最小值;(4)等離子體層電子含量PEC存在明顯的經(jīng)度變化,不同經(jīng)度鏈上的PEC存在不同的季節(jié)變化特征:120°E經(jīng)度鏈上的PEC表現(xiàn)為在兩分點達最大值的半年變化特征,而300°E經(jīng)度鏈上的PEC表現(xiàn)為在12月份至點季達最大值、在6月份至點季達最小值的年變化特征.

    致謝本文研究中所用到的COSMIC低軌衛(wèi)星對GPS的觀測數(shù)據(jù)來源于COSMIC數(shù)據(jù)分析和存檔中心(CDAAC, http:∥www.cosmic.ucar.edu/cdaac/).

    References

    Belehaki A, Jakowski N, Reinisch B W. 2004. Plasmaspheric electron content derived from GPS TEC and digisonde ionograms.Adv.SpaceRes., 33(6): 833-837. Bilitza D. 2001. International reference ionosphere 2000.RadioSci., 36(2): 261-275.

    Bilitza D, Reinisch B W. 2008. International reference ionosphere 2007: improvements and new parameters.Adv.SpaceRes., 42(4): 599-609, doi: 10.1016/j.asr.2007.07.048. Bilitza D, McKinnell L A, Reinisch B, et al. 2011. The international reference ionosphere today and in the future.J.Geod., 85(12): 909-920, doi: 10.1007/s00190-010-0427-x.

    Chen P, Yao Y B. 2015. Research on global plasmaspheric electron content by using LEO occultation and GPS data.Adv.SpaceRes., 55(9): 2248-2255.

    Chong X Y, Zhang M L, Zhang S R, et al. 2013. An investigation on plasmaspheric electron content derived from ISR and GPS observations at Millstone Hill.ChineseJ.Geophys. (in Chinese), 56(3): 738-745, doi: 10.6038/cjg20130303.Davíes K. 1990. Ionospheric Radio (IEEE Electromagnetic Waves Series, Vol. 31). London: Peter Peregrinus Ltd.

    Foelsche U, Kirchengast G. 2002. A simple “geometric” mapping function for the hydrostatic delay at radio frequencies and assessment of its performance.Geophys.Res.Lett., 29(10): 111-1-111-4, doi: 10.1029/2001GL013744.

    Leitinger R, Zhang M L, Radicella S M. 2005. An improved bottomside for the ionospheric electron density model NeQuick.AnnalsofGeophysics, 48(3): 525-534.

    Lunt N, Kersley L, Bishop G J, et al. 1999. The contribution of the protonosphere to GPS total electron content: experimental measurements.RadioSci., 34(5): 1273-1280.

    Nava B, Co?sson P, Radicella S M. 2008. A new version of the NeQuick ionosphere electron density model.J.Atmos.Solar-Terr.Phys., 70(15): 1856-1862.

    Reinisch B W, Huang X Q. 2001. Deducing topside profiles and total electron content from bottomside ionograms.Adv.SpaceRes., 27(1): 23-30.

    Rishbeth H. 1998. How the thermospheric circulation affects the ionospheric F2-layer.J.Atmos.Solar-Terr.Phys., 60(14): 1385-1402.

    Singh A K, Singh R P, Siingh D. 2011. State studies of Earth′s plasmasphere: A review.PlanetaryandSpaceScience, 59(9): 810-834. Torr M R, Torr D G. 1973. The seasonal behaviour of the F2-layer of the ionosphere.J.Atmos.Solar-Terr.Phys., 35(12): 2237-2251. Yizengaw E, Moldwin M B, Galvan D, et al. 2008. Global plasmaspheric TEC and its relative contribution to GPS TEC.J.Atmos.Solar-Terr.Phys., 70(11-12): 1541-1548. Yue X, Schreiner W S, Hunt D C, et al. 2011. Quantitative evaluation of the low Earth orbit satellite based slant total electron content determination.SpaceWeather, 9: S09001, doi: 10.1029/2011SW000687.

    Zhang M L, Radicella S M, Shi J, et al. 2006. Comparison among IRI, GPS-IGS and ionogram-derived total electron contents.Adv.SpaceRes., 37(5): 972-977.

    附中文參考文獻

    種小燕, 張滿蓮, 張順榮等. 2013. 基于非相干散射雷達和GPS觀測研究Millstone Hill地區(qū)等離子體層電子含量. 地球物理學報, 56(3): 738-745, doi: 10.6038/cjg20130303.

    (本文編輯何燕)

    基金項目國家自然科學基金(41274163)、國家重點基礎研究發(fā)展計劃(973計劃) (2012CB825604)、中國科學院重點部署項目(KZZD-EW-01-03)、國家自然科學基金創(chuàng)新群體項目(41321003)以及國家自然科學基金重點項目(41231065)資助.

    作者簡介張滿蓮,女,1962年生,副研究員,主要從事電離層等離子體層變化特性及經(jīng)驗建模研究. E-mail: zhangml@mail.iggcas.ac.cn

    doi:10.6038/cjg20160101 中圖分類號P352

    收稿日期2015-08-02,2015-09-21收修定稿

    Variation of the plasmaspheric electron content derived from the podTEC observations of COSMIC LEO satellites to GPS signals

    ZHANG Man-Lian, LIU Li-Bo, WAN Wei-Xing, NING Bai-Qi

    KeyLaboratoryofEarthandPlanetaryPhysics,InstituteofGeologyandGeophysics,ChineseAcademyofSciences,Beijing100029,China

    AbstractWe tried to study the variations of the plasmaspheric electron content (PEC) using the PEC data derived from the podTEC observation of the COSMIC low Earth orbit (LEO) satellite to the GPS satellite signals. We first give a brief introduce to the method we used to convert the slant podTEC to the vertical PEC. Then we used the converted PEC data of the year 2008 to study the variations of PEC with the geomagnetic latitude (MLAT), magnetic local time (MLT) and with four different seasons. Besides, we made a study on the longitudinal variation of PEC using the extracted PEC data from two different longitudes (120°E and 300°E). Our study showed that: (1) The distribution of PEC is mainly confined to a region within ±45° of the magnetic equator of the Earth; (2) PEC shows a well-defined diurnal variation pattern with higher values during daytime hours than during nighttime hours. PEC reaches its peak value at the hour around 12—16MLT, whereas it reaches its minimum value at around 4—5MLT. (3) PEC has a lowest value in the June solstice season (May—August) compared with other seasons. (4) PEC shows an evident longitudinal variation and it has different seasonal variations for different longitudes.

    KeywordsPlasmasphere; Ionosphere; Total Electron Content (TEC); Plasmaspheric electron content

    張滿蓮, 劉立波, 萬衛(wèi)星等. 2016. 利用COSMIC低軌衛(wèi)星對GPS信號的頂部TEC觀測資料研究等離子體層電子含量的變化特征.地球物理學報,59(1):1-7,doi:10.6038/cjg20160101.

    Zhang M L, Liu L B, Wan W X, et al. 2016. Variation of the plasmaspheric electron content derived from the podTEC observations of COSMIC LEO satellites to GPS signals.ChineseJ.Geophys. (in Chinese),59(1):1-7,doi:10.6038/cjg20160101.

    久久久久国产一级毛片高清牌| a在线观看视频网站| 精品熟女少妇八av免费久了| 国产片内射在线| 欧美日韩瑟瑟在线播放| 亚洲自拍偷在线| 久久久水蜜桃国产精品网| 麻豆av在线久日| 可以免费在线观看a视频的电影网站| 国产成人精品久久二区二区免费| 大型黄色视频在线免费观看| 淫秽高清视频在线观看| av片东京热男人的天堂| 亚洲专区字幕在线| 丁香欧美五月| 69av精品久久久久久| 国产成人一区二区三区免费视频网站| 国产v大片淫在线免费观看| 我的亚洲天堂| 国产伦一二天堂av在线观看| 大香蕉久久成人网| 丝袜人妻中文字幕| 亚洲精品久久国产高清桃花| 嫁个100分男人电影在线观看| 男女之事视频高清在线观看| 色播亚洲综合网| 免费观看人在逋| 日本a在线网址| 99国产精品一区二区蜜桃av| 老鸭窝网址在线观看| 好男人在线观看高清免费视频 | 久久亚洲真实| 国产精品综合久久久久久久免费| 91字幕亚洲| 91成人精品电影| 一区福利在线观看| 午夜a级毛片| 国产欧美日韩精品亚洲av| 黄色视频,在线免费观看| 国产精品久久久久久精品电影 | tocl精华| 99久久国产精品久久久| 午夜老司机福利片| 欧美黄色淫秽网站| 成人亚洲精品av一区二区| 久久热在线av| 成人欧美大片| 日日摸夜夜添夜夜添小说| 国产伦一二天堂av在线观看| 色老头精品视频在线观看| 亚洲成人久久性| 亚洲片人在线观看| 国产精品 国内视频| 法律面前人人平等表现在哪些方面| 成人国产综合亚洲| 宅男免费午夜| 好男人在线观看高清免费视频 | 欧美性猛交╳xxx乱大交人| 国产成人一区二区三区免费视频网站| 欧美 亚洲 国产 日韩一| 午夜影院日韩av| 亚洲性夜色夜夜综合| 国内精品久久久久久久电影| 欧美激情久久久久久爽电影| 国产私拍福利视频在线观看| 别揉我奶头~嗯~啊~动态视频| 夜夜躁狠狠躁天天躁| 精品国产超薄肉色丝袜足j| 亚洲第一青青草原| 国产极品粉嫩免费观看在线| 欧美日本视频| 精品国产一区二区三区四区第35| 大香蕉久久成人网| 我的亚洲天堂| 久久久久久人人人人人| 香蕉丝袜av| 亚洲一区高清亚洲精品| 女性被躁到高潮视频| 精华霜和精华液先用哪个| 满18在线观看网站| 久久伊人香网站| 亚洲五月色婷婷综合| 午夜免费观看网址| 亚洲一区二区三区不卡视频| 在线观看一区二区三区| 成年人黄色毛片网站| 久久青草综合色| 日韩免费av在线播放| 午夜福利免费观看在线| 国产亚洲欧美精品永久| 黄网站色视频无遮挡免费观看| 亚洲精品色激情综合| 麻豆成人午夜福利视频| 亚洲av成人一区二区三| 一a级毛片在线观看| 一个人观看的视频www高清免费观看 | 久久久久久久久中文| 国产男靠女视频免费网站| 欧美日韩亚洲综合一区二区三区_| 亚洲免费av在线视频| 国产免费男女视频| 少妇被粗大的猛进出69影院| 一级作爱视频免费观看| 中文资源天堂在线| 男女做爰动态图高潮gif福利片| 最近最新免费中文字幕在线| 啦啦啦观看免费观看视频高清| 午夜福利高清视频| 国产黄色小视频在线观看| 欧美成人免费av一区二区三区| 成人国产综合亚洲| АⅤ资源中文在线天堂| 亚洲免费av在线视频| 久久久久国产精品人妻aⅴ院| 男人舔女人下体高潮全视频| 男人操女人黄网站| 国产亚洲精品一区二区www| 18禁黄网站禁片免费观看直播| av超薄肉色丝袜交足视频| 九色国产91popny在线| 一区二区三区高清视频在线| 男女之事视频高清在线观看| 久久久久久久久免费视频了| 最好的美女福利视频网| 久久精品91无色码中文字幕| 啦啦啦韩国在线观看视频| 久久久久亚洲av毛片大全| 少妇裸体淫交视频免费看高清 | 精品国产乱子伦一区二区三区| 日本免费一区二区三区高清不卡| 中文字幕另类日韩欧美亚洲嫩草| 日本三级黄在线观看| 国产黄片美女视频| 亚洲精品美女久久av网站| 麻豆国产av国片精品| 一级毛片高清免费大全| 中文字幕人成人乱码亚洲影| 国产欧美日韩一区二区精品| 不卡av一区二区三区| 国产av一区二区精品久久| 国产色视频综合| 真人做人爱边吃奶动态| 国产男靠女视频免费网站| 亚洲av美国av| 国产成人一区二区三区免费视频网站| 日日摸夜夜添夜夜添小说| 999久久久精品免费观看国产| 久久久久久大精品| 国产精品久久电影中文字幕| 成人av一区二区三区在线看| 日日干狠狠操夜夜爽| 日韩欧美国产在线观看| 亚洲精品在线观看二区| 十八禁人妻一区二区| 中文字幕av电影在线播放| 亚洲精品一区av在线观看| 亚洲精品中文字幕一二三四区| 国产人伦9x9x在线观看| 欧美又色又爽又黄视频| 丝袜人妻中文字幕| 91在线观看av| 99热只有精品国产| 国产成人欧美在线观看| 久久国产精品影院| 黑丝袜美女国产一区| 特大巨黑吊av在线直播 | 国产av又大| 中文字幕高清在线视频| 精品国产美女av久久久久小说| 精品欧美国产一区二区三| 99久久99久久久精品蜜桃| 久久精品夜夜夜夜夜久久蜜豆 | 首页视频小说图片口味搜索| 国产一级毛片七仙女欲春2 | 不卡一级毛片| av视频在线观看入口| 精品久久久久久久久久免费视频| 美女免费视频网站| 亚洲成人久久爱视频| 岛国在线观看网站| 国产亚洲欧美在线一区二区| 在线播放国产精品三级| 国产野战对白在线观看| 成人手机av| 午夜日韩欧美国产| 国产精品免费一区二区三区在线| 狂野欧美激情性xxxx| 中文亚洲av片在线观看爽| 亚洲av日韩精品久久久久久密| 欧美成人免费av一区二区三区| 国产午夜福利久久久久久| 国内毛片毛片毛片毛片毛片| 怎么达到女性高潮| 欧美av亚洲av综合av国产av| 嫁个100分男人电影在线观看| 成人精品一区二区免费| 国产精品香港三级国产av潘金莲| 99热6这里只有精品| 欧美不卡视频在线免费观看 | 黄色片一级片一级黄色片| 欧美日韩瑟瑟在线播放| 欧美日韩中文字幕国产精品一区二区三区| 中文在线观看免费www的网站 | 国产亚洲精品久久久久5区| 亚洲天堂国产精品一区在线| 欧美在线黄色| av片东京热男人的天堂| 国产成人欧美| 亚洲五月天丁香| 亚洲在线自拍视频| 亚洲美女黄片视频| 操出白浆在线播放| 亚洲熟女毛片儿| a级毛片a级免费在线| √禁漫天堂资源中文www| 国产高清视频在线播放一区| 久久天躁狠狠躁夜夜2o2o| 成人精品一区二区免费| 久久婷婷人人爽人人干人人爱| 后天国语完整版免费观看| 亚洲va日本ⅴa欧美va伊人久久| 波多野结衣高清作品| 少妇的丰满在线观看| 国产99白浆流出| 亚洲精品中文字幕一二三四区| 国产v大片淫在线免费观看| 精品国产超薄肉色丝袜足j| 男女之事视频高清在线观看| 一级毛片高清免费大全| 亚洲久久久国产精品| 成熟少妇高潮喷水视频| 女人高潮潮喷娇喘18禁视频| 麻豆成人午夜福利视频| 午夜激情av网站| 久久久久久大精品| 免费在线观看日本一区| 久久精品国产清高在天天线| 成熟少妇高潮喷水视频| 韩国精品一区二区三区| 亚洲精华国产精华精| 50天的宝宝边吃奶边哭怎么回事| 久久99热这里只有精品18| 国内精品久久久久久久电影| 午夜视频精品福利| 国产精品 国内视频| 欧美日本视频| 国产成年人精品一区二区| 可以免费在线观看a视频的电影网站| 悠悠久久av| 欧美黑人巨大hd| 国产成人精品久久二区二区免费| 欧美黄色片欧美黄色片| 亚洲国产精品久久男人天堂| 美女午夜性视频免费| 亚洲成人久久爱视频| 日本成人三级电影网站| 夜夜爽天天搞| www日本在线高清视频| 亚洲av日韩精品久久久久久密| 天堂影院成人在线观看| 国产又色又爽无遮挡免费看| 成年免费大片在线观看| 欧美又色又爽又黄视频| 后天国语完整版免费观看| 欧美日韩精品网址| 亚洲一码二码三码区别大吗| 美国免费a级毛片| 19禁男女啪啪无遮挡网站| 国产av一区在线观看免费| 精品久久久久久,| 别揉我奶头~嗯~啊~动态视频| 黄色成人免费大全| av视频在线观看入口| 国产亚洲精品综合一区在线观看 | 后天国语完整版免费观看| 国产黄片美女视频| 免费观看人在逋| 丁香六月欧美| av欧美777| 亚洲色图 男人天堂 中文字幕| 99久久99久久久精品蜜桃| 国产黄片美女视频| 欧美乱码精品一区二区三区| 国产精品久久久av美女十八| 日韩高清综合在线| 天堂影院成人在线观看| 最好的美女福利视频网| 一区二区三区精品91| 欧美另类亚洲清纯唯美| 亚洲中文字幕一区二区三区有码在线看 | 亚洲成国产人片在线观看| 女同久久另类99精品国产91| 欧美日韩中文字幕国产精品一区二区三区| 熟女少妇亚洲综合色aaa.| 欧美成人午夜精品| 在线天堂中文资源库| 中出人妻视频一区二区| 国产99白浆流出| 两个人视频免费观看高清| 天堂影院成人在线观看| 神马国产精品三级电影在线观看 | 搡老妇女老女人老熟妇| 国产又黄又爽又无遮挡在线| 天天躁夜夜躁狠狠躁躁| 99久久无色码亚洲精品果冻| 成人一区二区视频在线观看| 伊人久久大香线蕉亚洲五| 听说在线观看完整版免费高清| 欧美日韩亚洲国产一区二区在线观看| 国产精品1区2区在线观看.| 黄片小视频在线播放| av福利片在线| 国产亚洲精品综合一区在线观看 | 国产精品,欧美在线| tocl精华| 亚洲av成人一区二区三| 欧美av亚洲av综合av国产av| 视频在线观看一区二区三区| 精品一区二区三区av网在线观看| 国产高清视频在线播放一区| 黄色视频,在线免费观看| 亚洲avbb在线观看| 亚洲精品国产区一区二| 亚洲精品粉嫩美女一区| 国产激情欧美一区二区| 国产一区二区在线av高清观看| 国产在线精品亚洲第一网站| 深夜精品福利| 黄色丝袜av网址大全| avwww免费| 亚洲第一欧美日韩一区二区三区| www日本黄色视频网| 国产亚洲欧美98| 亚洲av成人不卡在线观看播放网| 老司机午夜十八禁免费视频| 国产男靠女视频免费网站| 三级毛片av免费| 亚洲国产欧美日韩在线播放| 国产欧美日韩一区二区精品| 欧美一区二区精品小视频在线| 后天国语完整版免费观看| 精品国产超薄肉色丝袜足j| 十八禁人妻一区二区| 亚洲专区国产一区二区| 18禁裸乳无遮挡免费网站照片 | 亚洲一码二码三码区别大吗| 岛国在线观看网站| 免费人成视频x8x8入口观看| 国产精品免费视频内射| 久久久久九九精品影院| 亚洲精品在线美女| 好男人在线观看高清免费视频 | 首页视频小说图片口味搜索| 美女免费视频网站| 亚洲国产中文字幕在线视频| 久久99热这里只有精品18| 国产精品 欧美亚洲| 久久精品aⅴ一区二区三区四区| 校园春色视频在线观看| av欧美777| 国产精品久久电影中文字幕| 他把我摸到了高潮在线观看| 黄片播放在线免费| 波多野结衣高清作品| 50天的宝宝边吃奶边哭怎么回事| 国产午夜精品久久久久久| 脱女人内裤的视频| 成人免费观看视频高清| 搡老熟女国产l中国老女人| 精品久久久久久久毛片微露脸| 看黄色毛片网站| 国产亚洲欧美在线一区二区| 韩国av一区二区三区四区| 99国产综合亚洲精品| 在线观看舔阴道视频| 成在线人永久免费视频| 亚洲最大成人中文| a级毛片在线看网站| 成年版毛片免费区| 18禁裸乳无遮挡免费网站照片 | 国产精品久久久久久人妻精品电影| 欧美日韩瑟瑟在线播放| 精品久久久久久久末码| 黄网站色视频无遮挡免费观看| 午夜福利视频1000在线观看| 欧美亚洲日本最大视频资源| 亚洲成人久久爱视频| 97碰自拍视频| 国产日本99.免费观看| 2021天堂中文幕一二区在线观 | 黄片播放在线免费| 亚洲熟妇熟女久久| 精品国产国语对白av| 亚洲熟妇中文字幕五十中出| 亚洲精品av麻豆狂野| videosex国产| 无限看片的www在线观看| 少妇裸体淫交视频免费看高清 | 久久婷婷成人综合色麻豆| 国产精华一区二区三区| 悠悠久久av| 国产麻豆成人av免费视频| 午夜福利在线在线| 精品日产1卡2卡| 日本 av在线| 欧美国产日韩亚洲一区| 久久久久久久久久黄片| 黄色毛片三级朝国网站| 精品国产乱码久久久久久男人| 日本熟妇午夜| 国产精品久久视频播放| 999精品在线视频| bbb黄色大片| 天天躁夜夜躁狠狠躁躁| 桃色一区二区三区在线观看| avwww免费| 中文字幕最新亚洲高清| 日本精品一区二区三区蜜桃| 午夜日韩欧美国产| 久久精品国产亚洲av高清一级| 白带黄色成豆腐渣| 亚洲五月婷婷丁香| 成人国产一区最新在线观看| 国产单亲对白刺激| 国产精品久久久久久精品电影 | 免费观看人在逋| 成人国产综合亚洲| 国产精华一区二区三区| www.自偷自拍.com| 免费看日本二区| av在线播放免费不卡| www.熟女人妻精品国产| 亚洲熟妇中文字幕五十中出| netflix在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 黄色成人免费大全| 精品欧美国产一区二区三| www.自偷自拍.com| 老司机福利观看| 午夜福利一区二区在线看| 黄色毛片三级朝国网站| 亚洲第一av免费看| 精品欧美一区二区三区在线| 久久精品国产亚洲av香蕉五月| 中文在线观看免费www的网站 | 免费在线观看黄色视频的| 国产男靠女视频免费网站| 亚洲第一欧美日韩一区二区三区| 久久精品国产清高在天天线| 亚洲av中文字字幕乱码综合 | 精品少妇一区二区三区视频日本电影| 午夜精品久久久久久毛片777| 日韩精品青青久久久久久| 十分钟在线观看高清视频www| 国产一区在线观看成人免费| 成人18禁高潮啪啪吃奶动态图| 天堂√8在线中文| 亚洲国产欧美网| 久久精品国产亚洲av高清一级| 两性午夜刺激爽爽歪歪视频在线观看 | 国产av一区在线观看免费| 黄网站色视频无遮挡免费观看| 一个人观看的视频www高清免费观看 | 亚洲男人天堂网一区| 人人妻人人澡人人看| 亚洲成人久久性| 精品欧美国产一区二区三| 夜夜看夜夜爽夜夜摸| 1024视频免费在线观看| 亚洲色图av天堂| 亚洲第一电影网av| 老司机靠b影院| 亚洲成av人片免费观看| 国产亚洲av高清不卡| 国产精品久久久久久人妻精品电影| 搡老熟女国产l中国老女人| videosex国产| 欧美日本视频| 1024视频免费在线观看| 亚洲专区中文字幕在线| 校园春色视频在线观看| 久久精品影院6| 一本综合久久免费| 国产成人影院久久av| 日本免费一区二区三区高清不卡| 欧美成人免费av一区二区三区| 黄片小视频在线播放| 看黄色毛片网站| 大型av网站在线播放| 夜夜看夜夜爽夜夜摸| 亚洲人成伊人成综合网2020| 久久亚洲精品不卡| 一级毛片女人18水好多| 精品免费久久久久久久清纯| 亚洲电影在线观看av| 免费在线观看成人毛片| 国产主播在线观看一区二区| 免费在线观看日本一区| 亚洲精品中文字幕一二三四区| 国产亚洲精品久久久久久毛片| 免费在线观看黄色视频的| 巨乳人妻的诱惑在线观看| 欧美大码av| 亚洲人成77777在线视频| 日韩精品中文字幕看吧| 麻豆国产av国片精品| 国产成人系列免费观看| 精品久久久久久久久久免费视频| 久久人妻福利社区极品人妻图片| 又紧又爽又黄一区二区| 日本一区二区免费在线视频| avwww免费| 啦啦啦韩国在线观看视频| 亚洲av第一区精品v没综合| 在线观看www视频免费| 亚洲欧美激情综合另类| 久久精品国产亚洲av香蕉五月| 中文字幕av电影在线播放| 亚洲一码二码三码区别大吗| 一区二区三区高清视频在线| 搡老熟女国产l中国老女人| 亚洲国产欧美一区二区综合| www.自偷自拍.com| 欧美国产日韩亚洲一区| 国产精品爽爽va在线观看网站 | 国产精品爽爽va在线观看网站 | 日日干狠狠操夜夜爽| 亚洲中文字幕一区二区三区有码在线看 | 琪琪午夜伦伦电影理论片6080| 岛国在线观看网站| 中文字幕人妻熟女乱码| av欧美777| 丁香六月欧美| 一级片免费观看大全| 一进一出抽搐gif免费好疼| 亚洲av电影不卡..在线观看| 19禁男女啪啪无遮挡网站| 成人国产综合亚洲| 亚洲成人久久爱视频| 国产主播在线观看一区二区| 看黄色毛片网站| 麻豆av在线久日| 午夜免费成人在线视频| 久久热在线av| 亚洲欧美一区二区三区黑人| 一区二区三区高清视频在线| 午夜久久久在线观看| 亚洲中文av在线| 免费无遮挡裸体视频| 黄色片一级片一级黄色片| 亚洲精品久久成人aⅴ小说| 麻豆一二三区av精品| 精品久久久久久久末码| 婷婷精品国产亚洲av在线| 欧美大码av| 91av网站免费观看| 老汉色∧v一级毛片| 婷婷丁香在线五月| 国产一区二区激情短视频| 黄色成人免费大全| 99久久精品国产亚洲精品| 国产私拍福利视频在线观看| 中文字幕精品免费在线观看视频| 1024手机看黄色片| 男女下面进入的视频免费午夜 | 亚洲av成人av| 白带黄色成豆腐渣| 非洲黑人性xxxx精品又粗又长| 在线观看www视频免费| 亚洲第一欧美日韩一区二区三区| 国产黄a三级三级三级人| av免费在线观看网站| 免费在线观看成人毛片| 久久久久久久久久黄片| 757午夜福利合集在线观看| 一级作爱视频免费观看| 男女做爰动态图高潮gif福利片| 亚洲国产欧美网| 国产精华一区二区三区| 国产一区在线观看成人免费| 成人精品一区二区免费| 一本久久中文字幕| 18禁黄网站禁片免费观看直播| 日日摸夜夜添夜夜添小说| 1024视频免费在线观看| 久久国产乱子伦精品免费另类| 午夜免费鲁丝| 大型av网站在线播放| www.自偷自拍.com| xxx96com| 亚洲一码二码三码区别大吗| 此物有八面人人有两片| 久久午夜综合久久蜜桃| 香蕉丝袜av| 午夜视频精品福利| 国产高清videossex| 国产黄色小视频在线观看| 日韩免费av在线播放| 中文字幕最新亚洲高清| а√天堂www在线а√下载| 亚洲 国产 在线| 精品第一国产精品| 亚洲av熟女| 成人国产一区最新在线观看| 九色国产91popny在线| 69av精品久久久久久| 嫩草影视91久久| 日本一本二区三区精品| 亚洲av熟女| av中文乱码字幕在线| 啦啦啦观看免费观看视频高清| 最近最新免费中文字幕在线| 男女午夜视频在线观看| 可以在线观看毛片的网站| 免费看十八禁软件|