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

    X射線耀斑期間電離層VTEC的時變率

    2016-11-04 02:35:17孫中苗朱永興劉曉剛
    測繪科學與工程 2016年3期
    關(guān)鍵詞:耀斑頂角電離層

    管 斌,孫中苗,朱永興,劉曉剛

    1.信息工程大學地理空間信息學院,河南 鄭州,450052;2.地理信息工程國家重點實驗室,陜西 西安,710054;3.西安測繪研究所,陜西 西安,710054

    ?

    X射線耀斑期間電離層VTEC的時變率

    管斌1,2,3,孫中苗2,3,朱永興2,3,劉曉剛2,3

    1.信息工程大學地理空間信息學院,河南 鄭州,450052;2.地理信息工程國家重點實驗室,陜西 西安,710054;3.西安測繪研究所,陜西 西安,710054

    利用通過GEO衛(wèi)星雙頻載波相位觀測量計算電離層垂直總電子含量時變率(RVTEC)的方法,研究了RVTEC對太陽X射線耀斑的響應。根據(jù)載波相位觀測方程,得到了RVTEC的計算公式。選取3個MGEX測站對2013~2015年間9次X級耀斑的觀測數(shù)據(jù)進行了計算分析,主要結(jié)果包括:①X級耀斑期間RVTEC不一定發(fā)生明顯的增大;②RVTEC與耀斑的級別、太陽天頂角沒有絕對的對應關(guān)系;③所有X3.2級以下耀斑觀測結(jié)果中RVTEC最大不超過0.03TECu/s。

    X射線耀斑;垂直總電子含量變化率;地球靜止軌道;載波相位;電離層

    1 引 言

    地球高層大氣在太陽X射線和極紫外輻射的作用下被加熱、解離,還被部分電離形成了電離層。太陽耀斑爆發(fā)期間向外界輻射多種頻段的電磁輻射,其中的紫外輻射和X射線會造成電離層電子密度的增加,引起多種電離層擾動現(xiàn)象[1],這些現(xiàn)象包括突然頻率偏移、突然相位異常、突然頻率吸收、耀斑的地磁效應及突然總電子含量增加等。

    根據(jù)觀測手段的不同,耀斑主要分為光學耀斑、X射線耀斑等。電離層總電子含量(total electron content, TEC)對X射線耀斑的響應廣受關(guān)注,近年來許多學者通過GPS觀測數(shù)據(jù)對其進行了研究。文獻[2]研究了耀斑位置與TEC突增的關(guān)系,認為TEC增加不僅與耀斑級別有關(guān)系,而且與耀斑日面位置也有關(guān)系。文獻[3]認為耀斑引起的TEC變化率與耀斑的輻射能量成正比,與Chapman函數(shù)(約等于sec(χ),χ為太陽天頂角)成反比,即天頂角越大,TEC變化率越小。文獻[4]認為耀斑期間向日面電離層出現(xiàn)了總電子含量突增事件,并得到“太陽天頂角越大,TEC增幅越小”的結(jié)論。文獻[5]利用中低緯度的GPS跟蹤站在太陽耀斑發(fā)生期間均監(jiān)測到了總電子含量的突增現(xiàn)象,且與X射線輻射通量具有很好的一致性。文獻[6]利用多個測站多條視線方向GPS觀測數(shù)據(jù)的平均值,得到較為明顯的TEC對X級耀斑的響應。文獻[7,8]等也進行了相關(guān)的研究,在此不再贅述。

    上述研究對TEC變化量進行了系統(tǒng)的分析,但對TEC變化率僅描述了現(xiàn)象,沒有定性分析。然而TEC變化率是反映電離層動態(tài)變化的重要參數(shù)之一,同時垂直總電子含量(vertical TEC, VTEC)的時變率(rate of change of VTEC, RVTEC)對X射線耀斑的最大響應值得探究,因此,本文研究了RVTEC對太陽耀斑的響應。由于在太陽耀斑中電離層主要受X射線耀斑的影響,因而本文僅研究RVTEC對X射線耀斑的響應(后文中耀斑均特指X射線太陽耀斑)。

    耀斑期間TEC響應的研究大都通過GPS觀測數(shù)據(jù)來實現(xiàn)[1-8],所得結(jié)果包含了時變部分與空間變化部分。北斗衛(wèi)星導航系統(tǒng)(BeiDou Navigation Satellite System, BDS)中包含了5顆地球靜止軌道(GEO)衛(wèi)星,由于GEO衛(wèi)星的相對靜止特性,利用對GEO的觀測數(shù)據(jù)可以直接得到VTEC的時變率而排除空間變化部分的影響。另外,現(xiàn)有文獻通過GPS觀測數(shù)據(jù)分析TEC的變化率通常采用1 min或30 s采樣間隔的載波相位觀測數(shù)據(jù),然而耀斑的持續(xù)時間由幾分鐘到幾十分鐘不等,利用以上采樣間隔的數(shù)據(jù)無法得到VTEC的小尺度變化細節(jié)。本文采用1s采樣間隔的觀測數(shù)據(jù),可以得到RVTEC更加細致的變化結(jié)果。

    綜上,本文通過GEO衛(wèi)星雙頻載波相位觀測數(shù)據(jù),研究VTEC時變率對太陽X射線耀斑的響應。

    2 計算方法

    文獻[9]給出了通過GPS觀測數(shù)據(jù)計算VTEC在相鄰兩歷元穿刺點上差值的公式:

    ΔVTEC=F[(λ1φ1-λ2φ2)i+1M(zi+1)-(λ1φ1-λ2φ2)iM(zi)+(λ1N1-λ2N2+R1-R2)·(M(zi+1)-M(zi))]

    (1)

    (2)

    其中R是地球平均半徑,α=0.9782,H是電離層單層模型的高度,本文選定H為350km。

    式(1)同樣適用于對BDS衛(wèi)星的觀測。由于GEO的靜地特性,對GEO衛(wèi)星觀測結(jié)果具有特殊之處,即衛(wèi)星天頂角距z在相鄰兩歷元幾乎無變化,有zi+1=zi,故式(1)可以簡化為

    ΔVTEC=F[(λ1φ1-λ2φ2)i+1-(λ1φ1-λ2φ2)i]M(zi)

    (3)

    設(shè)載波相位觀測值的變化量分別為Δφ1=φ1i+1-φ1i,Δφ2=φ2i+1-φ2i,則

    ΔVTEC=F[(λ1Δφ1-λ2Δφ2)]M(zi)

    (4)

    將ΔVTEC除以歷元之間的間隔時間Δt即可得到RVTEC,如式(5):

    (5)

    由式(5),RVTEC的誤差主要來源于載波相位相對變化量的測量誤差,對實測數(shù)據(jù)的誤差分析結(jié)果表明,RVTEC的精度約為4×10-3TECu/s。

    3 數(shù)據(jù)選取

    太陽耀斑按照其X射線峰值的流量可分為A、B、C、M、X五級,所釋放能量依次增大。由于M級以下耀斑對電離層造成的影響不易通過對GNSS衛(wèi)星進行觀測的方式獲取[6],而本文注重于分析RVTEC對耀斑的最大響應,因而本文僅研究X級耀斑的影響。

    采用IGS[12]的MGEX[13]項目數(shù)據(jù)進行計算分析。MGEX項目是多模多頻GNSS實驗項目,其測站具有對多種衛(wèi)星導航系統(tǒng)觀測的能力。該項目從2012年1月開始實施,2012年底才開始陸續(xù)有對BDS信號的觀測,因而本文僅對2013年后的耀斑事件進行分析。設(shè)定30°為對衛(wèi)星觀測的截止仰角,選擇對BDS具有觀測能力并處于能夠有效觀測GEO衛(wèi)星區(qū)域的測站,測站信息見表1。表1中“起始時間”為各測站對BDS具有觀測能力的時間(自2013年始,格式為yymmdd);LT代表當?shù)貢r間,UT代表世界協(xié)調(diào)時。另有兩測站PTGG與SEYG,投入觀測時間較晚,雖然具備對BDS的觀測能力,然而投入運行后的幾次X級耀斑發(fā)生時均不在有效觀測范圍,故表1中未將其列入。

    美國發(fā)射的一系列GOES(Geostationary Operational Environmental Satellite)衛(wèi)星可持續(xù)觀測太陽X射線能量[14],根據(jù)美國空間天氣預報中心(Space Weather Prediction Center, SWPC)發(fā)布的源于GOES衛(wèi)星觀測數(shù)據(jù)的X射線太陽耀斑記錄,2013至2015年間,X級耀斑發(fā)生的次數(shù)分別為12、16、2,耀斑最強的為2014年2月25日發(fā)生的X4.9級耀斑。

    表1測站信息一覽表

    序號測站名起始時間測站坐標30°以上的GEO衛(wèi)星LT1JFNG130101114.491°E,30.515°NGEO01,GEO02,GEO03UT+8h2MAYG14012145.258°E,12.782°SGEO02,GEO05UT+3h3NKLG1403019.672°E,0.354°NGEO05UT+1h

    圖1 2013年至2015年期間X級耀斑發(fā)生情況

    由于耀斑僅對處于向日面地區(qū)的電離層產(chǎn)生影響,TEC增量的大小與耀斑爆發(fā)地方時有直接關(guān)系,中午時段比晚間更大[15],且太陽天頂角越大,TEC變化率越小,TEC增幅越小[3,4],而本文主要研究耀斑對RVTEC的最大影響,因而僅對耀斑發(fā)生期間處于向日面并限定LT處于9:00至15:00之間(后文中稱以上條件為“觀測條件”)的測站進行分析。結(jié)合表1中測站信息以及GOES衛(wèi)星的觀測記錄,可得3年間全部X級耀斑的發(fā)生情況如圖1所示。

    圖1中,“不可有效觀測”代表爆發(fā)期間可用測站均不滿足“觀測條件”的耀斑;“無IGS觀測數(shù)據(jù)”代表爆發(fā)期間有測站滿足“觀測條件”,但由于未知因素影響(可能恰為所發(fā)生耀斑的影響),IGS觀測數(shù)據(jù)質(zhì)量太差而未發(fā)布觀測數(shù)據(jù)的耀斑(如圖中X4.9級耀斑),“觀測數(shù)據(jù)有效”代表爆發(fā)期間有測站滿足“觀測條件”且觀測數(shù)據(jù)有效的耀斑?,F(xiàn)將圖1中“觀測數(shù)據(jù)有效”的X級耀斑信息列于表2,其中時間均為UT時間,格式為hhss,“峰值時間”指X射線通量達到最大值的時間,χ為峰值時間測站觀測太陽的天頂角,其計算方法見文獻[16]。

    表2用于分析的2013~2015年間耀斑信息

    序號發(fā)生日期開始時間結(jié)束時間峰值時間量級觀測測站χ(°)a1305130153023202171.7JFNG29.9b1305140000012001113.2JFNG43.8c1305150125015801481.2JFNG23.4d1311100508051805141.1JFNG50.5e1406101136114411422.2MAYG53.7f1406101136114411422.2NKLG23.3g1406101236130312521.5NKLG31.7h1406110859091009061.0MAYG35.9i1410190417054805031.1JFNG42.8j1410221402145014281.6NKLG51.6

    4 計算與分析

    根據(jù)文獻[3,4],當強烈的X射線耀斑爆發(fā)時,TEC的變化與耀斑的爆發(fā)呈現(xiàn)良好的一致性,即伴隨耀斑X射線的到達,TEC明顯增加,伴隨耀斑的結(jié)束,TEC開始緩慢變化,因而本文僅對耀斑發(fā)生前后一小段時間內(nèi)RVTEC進行分析,對應于表2中各測站的RVTEC響應,如圖2所示。為了將測站對不同GEO衛(wèi)星的觀測結(jié)果呈現(xiàn)于一張圖中,圖2中分別對不同衛(wèi)星的RVTEC加減了相應常數(shù)以示區(qū)分。

    圖2  X級耀斑期間的RVTEC響應

    由表2與圖2知,圖2(a)、(b)、(c)、(g)、(h)、(i)所對應的耀斑發(fā)生期間,測站所觀測的RVTEC并沒有明顯地發(fā)生變化;圖2(d)、(e)、(f)、(j)所對應的耀斑發(fā)生期間有一小段時間,RVTEC發(fā)生了相對明顯的增大。其中,圖2(d)中所觀測到的RVTEC變化是最為明顯的,對三顆GEO衛(wèi)星觀測的RVTEC變化呈現(xiàn)高度一致性,并與耀斑的發(fā)生情況(見表2)相吻合,與之相對應的VTEC相對變化情況如圖3(d)所示。由圖3(d)知,該次耀斑發(fā)生期間JFNG測站附近區(qū)域的VTEC約增加了2TECu。圖3(e)、(f)、(j)為分別對應于圖2(e)、(f)、(j)的VTEC相對變化情況,圖2中其它圖所示耀斑發(fā)生期間VTEC的變化情況不再贅述。

    圖3 耀斑期間VTEC的相對變化

    綜合表2、圖2、圖3可得:

    ①耀斑發(fā)生期間RVTEC并不一定發(fā)生明顯的增大,如圖2(a),雖然X2.8級耀斑發(fā)生期間JFNG測站觀測太陽的天頂角約30°,但RVTEC并未呈現(xiàn)明顯變化。

    ②即使耀斑發(fā)生期間RVTEC出現(xiàn)了明顯的增大,但RVTEC持續(xù)高出正常水平并不貫穿于整個耀斑過程,如圖2(j)所對應的X1.6級耀斑開始于14:02,結(jié)束于14:50并于14:28達到射線通量最大值,但NKLG測站觀測的RVTEC僅在耀斑開始階段有2 min高出正常水平,后回歸到緩慢變化狀態(tài),在圖2(d)、(e)、(f)中具有相同的體現(xiàn)。

    ③耀斑所引起的VTEC的最快變化并不一定發(fā)生在X射線通量最大的時刻,如圖2(j),耀斑于14:28達到射線通量最大值,然而此時RVTEC已經(jīng)恢復到正常變化水平。

    ④耀斑的發(fā)生對RVTEC的影響不可一概而論,相同級別的耀斑對于太陽天頂角相同的地區(qū)所造成的影響可以不同,強度高的耀斑可能并不會引起明顯的RVTEC變化(如圖2(b)所示X3.2級耀斑)。耀斑發(fā)生期間RVTEC與耀斑的級別、太陽的天頂角距沒有絕對的對應關(guān)系,耀斑對RVTEC的影響勢必受空間環(huán)境中其它一系列因素的影響。

    ⑤圖2(e)與圖2(f)分別為MAYG與NKLG兩個測站對同一次耀斑進行觀測計算得到的RVTEC變化序列,耀斑發(fā)生期間分別對應于兩個測站地方時的14:30與12:30,VTEC均呈下降趨勢,期間太陽天頂角分別對53.7°與23.3°,而由圖3(e)與圖3(f)知,VTEC分別最多增大了0.5TECu與0.7TECu。從而,對同一耀斑發(fā)生期間處于向日面不同位置2個測站的觀測結(jié)果比較顯示,即使當VTEC呈下降趨勢時,在耀斑的影響下VTEC也可能階段性增多,且太陽天頂角較小時VTEC增加的幅度較大。

    ⑥耀斑會引起RVTEC的明顯變化,然而概率并不高。綜合所有可用觀測數(shù)據(jù),X3.2級以下耀斑觀測結(jié)果中RVTEC最大不超過0.03TECu/s。

    將以上方法應用于耀斑期間RVTEC的研究,雖然能夠分離開VTEC在空間上變化的影響,直觀地得到耀斑發(fā)生期間VTEC的時變情況,然而也存在一些缺點。主要體現(xiàn)在:①由于目前長期進行BDS系統(tǒng)數(shù)據(jù)觀測并公開數(shù)據(jù)的測站較少,因而直接獲取VTEC時變率的可用觀測數(shù)據(jù)較少;觀測數(shù)據(jù)不足以分析VTEC時變率與耀斑的強度、太陽天頂角之間的關(guān)系。②耀斑發(fā)生期間,向日面測站的觀測數(shù)據(jù)常常因為質(zhì)量較差而不進行發(fā)布,在特大耀斑發(fā)生時,該方法可能因獲取不到數(shù)據(jù)而無法計算RVTEC。由于可用測站有限,電離層穿刺點的分布受限,所得結(jié)果中并不完全包括耀斑發(fā)生期間RVTEC的最大變化情況,相關(guān)的分析有待通過布設(shè)更多的測站進行更全面的觀測來實現(xiàn)。

    5 小 結(jié)

    X射線太陽耀斑發(fā)生時,向日面地球電離層受其影響電子密度會發(fā)生不同程度的增加,而對該問題的研究大都集中于對TEC變化量的分析,沒有關(guān)于TEC變化率響應的定性研究。本文根據(jù)對BDS的GEO衛(wèi)星1s采樣間隔的雙頻載波相位觀測數(shù)據(jù),計算分析了X級X射線耀斑期間RVTEC的變化情況。篩選了3個MGEX項目測站對2013~2015年間的9次耀斑進行分析,得到了RVTEC對X射線耀斑的個例響應特性,所得結(jié)論有益于對電離層變化特性的系統(tǒng)研究。

    [1]Libo Liu, Chen Y. Statistical Analysis of Solar Activity Variations of Total Electron Content Derived at Jet Propulsion Laboratory from GPS Observations[J]. J Geophys Res, 2009,114(A10311).

    [2]Zhang Donghe, Xiao Zuo, Chang Qing. The Correlation of Flare's Location on Solar Disc and the Sudden Increase of Total Electron Content[J]. Chinese Science Bulletin, 2002,47(1):83-85.

    [3]Wan WeiXing, Yuan Hong, Liu Libo, et al. The Sudden Increase in Ionospheric Total Electron Content Caused by the Very Intense Solar Flare on July 14, 2000[J]. SCIENCE IN CHINA (Series A), 2002,45:142-147.

    [4]張東和,孫臻,劉四清等. 大耀斑期間向日面電離層總電子含量的響應個例分析[J]. 空間科學學報, 2004,24(2):96-104.

    [5]劉旭春,易武,李維功等. 利用GPS對太陽耀斑進行監(jiān)測及預報的方法研究[J]. 武漢大學學報·信息科學版, 2008,33(5):484-487.

    [6]黃文耿,古士芬,陳艷紅等. 電離層TEC對小耀斑的響應[J]. 天文學報, 2006,47(2):212-216.

    [7]崔書珍,周金國,彭軍還. 基于GPS的2003年10月28日太陽耀斑的電離層響應研究[J]. 測繪科學, 2009,34(5):55-56.

    [8]劉旭春,張正祿,楊軍. 利用GPS監(jiān)測兩起太陽耀斑事件的結(jié)果分析[J]. 測繪科學, 2007,32(5):11-13.

    [9]Zhang Hongping, Zhu Wenyao, Peng Junhuan, et al. Analysis of the Ionosphere Wave-motion with GPS[J]. Chinese Science Bulletin, 2005,50(12):1373-1381.

    [10]安家春,章迪,杜玉軍等. 極區(qū)電離層梯度的特性分析[J]. 武漢大學學報·信息科學版, 2014,39(1):75-79.

    [11]Denise Dettmering, Marco Limberger, Michael Schmidt. Using DORIS Measurements for Modeling the Vertical Total Electron Content of the Earth's Ionosphere[J]. J Geod, 2014,88:1131-1143.

    [12]Dow J.M., Neilan R.E., Rizos C. The International GNSS Service in a Changing Landscape of Global Navigation Satellite Systems[J]. Journal of Geodesy, 2009,83(3-4):191-198.

    [13]Montenbruck O., Steigenberger P., Khachikyan R., et al. IGS-MGEX : Preparing the Ground for Multi-constellation GNSS Science[J]. Inside GNSS, 2014,9(1):42-49.

    [14]于超,李嘉巍,張效信等. 風云二號C/D衛(wèi)星與GOES衛(wèi)星太陽X射線探測數(shù)據(jù)交叉比對[J]. 地球物理學報, 2012,55(9):2835-2842.

    [15]黃江,林果果,鄧柏昌等. 電離層測高數(shù)據(jù)對太陽耀斑響應的研究[J]. 華北地震科學, 2013,31(4):22-26.

    [16]王存良. 電離層電子濃度計算模型中有關(guān)太陽參數(shù)的計算[J]. 中國電子科學研究院學報, 2013,8(1):86-90.

    Time Variation Rate of Ionospheric Vertical Total Electron Content in X-Ray Flares

    Guan Bin1,2,3, Sun Zhongmiao2,3, Zhu Yongxing2,3, Liu Xiaogang2,3

    1.Institute of Geospatial Information, Information Engineering University, Zhengzhou 450052, China 2. Xi’an Research Institute of Surveying and Mapping, Xi’an 710054, China 3. State Key Laboratory of Geo-information Engineering, Xi’an 710054, China

    Time variation rate of ionospheric vertical total electron content (RVTEC) is calculated using dual-frequency carrier-phase observations of GEO satellites and its response to solar X-ray flare is studied in this paper. The formula for calculating RVTEC is deduced according to carrier-phase observation equation. Three MGEX stations are selected to calculate the observation data of nine times of X-class flares occurred during 2013 to 2015, and the results show that: (1) RVTEC does not always increase obviously during X-class flares; (2) there is no absolute relationship between RVTEC and the level of flare or solar vertex angle; (3) the maximal observation results of RVTEC is 0.03 TECu/s for flares below X3.2 level.

    X-ray flare; time variation rate of vertical total electron content; geostationary orbit; carrier-phase; ionosphere

    2016-02-01。

    國家自然科學基金資助項目(41174017),國家自然科學基金青年科學基金資助項目(41304022)。

    管斌(1988—),男,工程師,主要從事衛(wèi)星測高與衛(wèi)星導航方面的研究。

    P228

    A

    猜你喜歡
    耀斑頂角電離層
    一般三棱鏡最大頂角與折射率的關(guān)系
    一種電離層TEC格點預測模型
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    太陽表面的微型耀斑
    涼亭中的數(shù)學
    “最強”太陽耀斑爆發(fā)會帶來什么?
    頂角為100°的等腰三角形性質(zhì)的應用
    電離層對中高軌SAR影響機理研究
    雷達學報(2017年6期)2017-03-26 07:53:02
    Linux Shell語言在電離層解算中的應用
    太陽耀斑
    一个人免费在线观看的高清视频| 18禁国产床啪视频网站| 亚洲精品一卡2卡三卡4卡5卡| 欧美久久黑人一区二区| x7x7x7水蜜桃| 国产精品av久久久久免费| 精品国产一区二区三区四区第35| 女性被躁到高潮视频| 日本在线视频免费播放| 两性夫妻黄色片| 哪里可以看免费的av片| 日韩成人在线观看一区二区三区| 老司机在亚洲福利影院| 级片在线观看| 欧美色视频一区免费| 久久中文看片网| www日本在线高清视频| 国产亚洲欧美精品永久| 桃红色精品国产亚洲av| 久久人妻福利社区极品人妻图片| 国产精品久久久人人做人人爽| 一区二区三区高清视频在线| 男男h啪啪无遮挡| 亚洲avbb在线观看| 一边摸一边抽搐一进一小说| 欧美黄色片欧美黄色片| 国产三级在线视频| 天天一区二区日本电影三级| 国产精品一区二区免费欧美| av片东京热男人的天堂| 久久午夜亚洲精品久久| 啪啪无遮挡十八禁网站| 久久精品国产亚洲av香蕉五月| 成人特级黄色片久久久久久久| 国产激情欧美一区二区| 国产精品久久电影中文字幕| 一区二区三区精品91| 19禁男女啪啪无遮挡网站| 亚洲欧美日韩无卡精品| 妹子高潮喷水视频| av电影中文网址| 一卡2卡三卡四卡精品乱码亚洲| 禁无遮挡网站| 叶爱在线成人免费视频播放| 欧美激情 高清一区二区三区| 神马国产精品三级电影在线观看 | 日本 欧美在线| 天堂动漫精品| 亚洲自拍偷在线| 亚洲专区中文字幕在线| 婷婷亚洲欧美| 黄色毛片三级朝国网站| 色av中文字幕| 国产麻豆成人av免费视频| 夜夜爽天天搞| 久久精品影院6| 国产在线观看jvid| av中文乱码字幕在线| 国产亚洲av高清不卡| 看黄色毛片网站| 老司机福利观看| 少妇的丰满在线观看| 久久久久国内视频| 欧美日韩乱码在线| 欧美黑人欧美精品刺激| 男人舔女人下体高潮全视频| 999精品在线视频| 日本熟妇午夜| 成人一区二区视频在线观看| www.精华液| 日韩有码中文字幕| 亚洲欧美一区二区三区黑人| АⅤ资源中文在线天堂| 男女床上黄色一级片免费看| 午夜福利18| 久久久久久免费高清国产稀缺| 熟妇人妻久久中文字幕3abv| 精品少妇一区二区三区视频日本电影| 真人做人爱边吃奶动态| 一本精品99久久精品77| 欧美另类亚洲清纯唯美| 亚洲成人久久爱视频| 色精品久久人妻99蜜桃| 99热6这里只有精品| 国产伦人伦偷精品视频| 变态另类丝袜制服| 国产91精品成人一区二区三区| 亚洲国产毛片av蜜桃av| 波多野结衣高清无吗| 黄色成人免费大全| 一本综合久久免费| 一级毛片精品| 99国产精品99久久久久| 两个人看的免费小视频| 一卡2卡三卡四卡精品乱码亚洲| 免费在线观看成人毛片| 亚洲成人久久爱视频| 欧美黑人精品巨大| 亚洲熟女毛片儿| 国产99久久九九免费精品| 色尼玛亚洲综合影院| 色av中文字幕| 国产日本99.免费观看| 国产真实乱freesex| 俺也久久电影网| 99久久无色码亚洲精品果冻| 丝袜美腿诱惑在线| 哪里可以看免费的av片| 国产精品,欧美在线| 国产欧美日韩一区二区三| 波多野结衣高清无吗| 999精品在线视频| 国产精品亚洲一级av第二区| 婷婷精品国产亚洲av| 搞女人的毛片| 亚洲中文日韩欧美视频| 欧美最黄视频在线播放免费| 久久国产亚洲av麻豆专区| 波多野结衣巨乳人妻| 精品欧美国产一区二区三| 国产欧美日韩精品亚洲av| 国产私拍福利视频在线观看| 91麻豆精品激情在线观看国产| 亚洲男人天堂网一区| 一区二区日韩欧美中文字幕| 亚洲成人精品中文字幕电影| 国产av在哪里看| 99国产极品粉嫩在线观看| 国产精品久久久久久亚洲av鲁大| 嫩草影视91久久| 亚洲全国av大片| 神马国产精品三级电影在线观看 | 日本熟妇午夜| 国语自产精品视频在线第100页| 欧洲精品卡2卡3卡4卡5卡区| 国产精品亚洲av一区麻豆| 91麻豆av在线| 亚洲欧美精品综合久久99| 亚洲avbb在线观看| 欧美中文综合在线视频| www.精华液| 欧美国产日韩亚洲一区| 国产高清videossex| 色尼玛亚洲综合影院| 高清在线国产一区| 淫妇啪啪啪对白视频| 看免费av毛片| 午夜日韩欧美国产| 91在线观看av| 91九色精品人成在线观看| 9191精品国产免费久久| 无限看片的www在线观看| 成人精品一区二区免费| 男男h啪啪无遮挡| 亚洲国产日韩欧美精品在线观看 | 91九色精品人成在线观看| 色婷婷久久久亚洲欧美| 99国产精品一区二区三区| 曰老女人黄片| 国产亚洲av高清不卡| 12—13女人毛片做爰片一| 此物有八面人人有两片| 此物有八面人人有两片| 久久久国产欧美日韩av| 亚洲精品一区av在线观看| 日韩精品青青久久久久久| 此物有八面人人有两片| 欧美日本亚洲视频在线播放| 久久久久亚洲av毛片大全| 2021天堂中文幕一二区在线观 | 好男人在线观看高清免费视频 | 久久香蕉精品热| 精品一区二区三区视频在线观看免费| 午夜两性在线视频| 欧美日韩黄片免| 每晚都被弄得嗷嗷叫到高潮| 999久久久国产精品视频| 999久久久国产精品视频| 女人高潮潮喷娇喘18禁视频| 韩国精品一区二区三区| 好男人电影高清在线观看| 精品人妻1区二区| 丁香欧美五月| av欧美777| 日韩欧美一区二区三区在线观看| 欧美午夜高清在线| 亚洲男人天堂网一区| 国产av不卡久久| 久久久水蜜桃国产精品网| 变态另类丝袜制服| 久久久精品欧美日韩精品| 欧美性猛交╳xxx乱大交人| 午夜福利视频1000在线观看| 女性生殖器流出的白浆| 搡老岳熟女国产| 两个人免费观看高清视频| 久久天堂一区二区三区四区| 男女床上黄色一级片免费看| videosex国产| a级毛片在线看网站| 欧美成人一区二区免费高清观看 | 在线十欧美十亚洲十日本专区| 亚洲精品在线美女| 每晚都被弄得嗷嗷叫到高潮| 久久久久亚洲av毛片大全| 成人三级黄色视频| 在线免费观看的www视频| 国产精品自产拍在线观看55亚洲| 午夜亚洲福利在线播放| 欧美日韩福利视频一区二区| 老司机午夜福利在线观看视频| 性色av乱码一区二区三区2| 99国产精品一区二区蜜桃av| 久久天躁狠狠躁夜夜2o2o| 麻豆国产av国片精品| 在线看三级毛片| 法律面前人人平等表现在哪些方面| 级片在线观看| 欧美最黄视频在线播放免费| 精品卡一卡二卡四卡免费| 9191精品国产免费久久| 亚洲色图av天堂| 国产精品亚洲美女久久久| 国产精品一区二区精品视频观看| 麻豆成人av在线观看| 不卡一级毛片| 女人高潮潮喷娇喘18禁视频| 久久亚洲精品不卡| 在线免费观看的www视频| 淫秽高清视频在线观看| 最近最新中文字幕大全免费视频| www.www免费av| 亚洲人成伊人成综合网2020| 欧美中文综合在线视频| 亚洲在线自拍视频| 亚洲av片天天在线观看| 白带黄色成豆腐渣| 亚洲成人久久性| 97人妻精品一区二区三区麻豆 | 久久中文字幕一级| 1024视频免费在线观看| 欧美色视频一区免费| 久久久久国产一级毛片高清牌| 欧美久久黑人一区二区| 在线国产一区二区在线| 欧美日韩瑟瑟在线播放| 日本黄色视频三级网站网址| 韩国精品一区二区三区| 久久国产精品男人的天堂亚洲| 国产私拍福利视频在线观看| 久久久久久国产a免费观看| 老鸭窝网址在线观看| 国产成人一区二区三区免费视频网站| 精品午夜福利视频在线观看一区| 91大片在线观看| 一级毛片精品| 手机成人av网站| 日本免费一区二区三区高清不卡| 午夜a级毛片| 欧美日韩中文字幕国产精品一区二区三区| 国产成人精品无人区| 中文字幕人妻熟女乱码| 欧美成人免费av一区二区三区| 国产精品野战在线观看| 色婷婷久久久亚洲欧美| 身体一侧抽搐| 一级a爱片免费观看的视频| 一二三四社区在线视频社区8| 精品久久久久久久末码| 无人区码免费观看不卡| 色尼玛亚洲综合影院| 精品久久久久久久人妻蜜臀av| 国产伦一二天堂av在线观看| 欧美亚洲日本最大视频资源| 美女免费视频网站| 国产午夜精品久久久久久| 天堂√8在线中文| 久久久久久久精品吃奶| 99久久国产精品久久久| 亚洲av电影不卡..在线观看| 国产精品野战在线观看| 在线观看66精品国产| 精品人妻1区二区| 精品乱码久久久久久99久播| 免费女性裸体啪啪无遮挡网站| 制服诱惑二区| 亚洲av电影不卡..在线观看| 中文字幕高清在线视频| 国产亚洲av高清不卡| 欧美av亚洲av综合av国产av| 十八禁网站免费在线| 中亚洲国语对白在线视频| 很黄的视频免费| 日本一本二区三区精品| 岛国在线观看网站| 国产精品久久久久久亚洲av鲁大| 欧美激情 高清一区二区三区| 少妇裸体淫交视频免费看高清 | 久久99热这里只有精品18| 亚洲av电影在线进入| 一本久久中文字幕| 色av中文字幕| 色老头精品视频在线观看| 男人舔女人的私密视频| 99热只有精品国产| 日日干狠狠操夜夜爽| 午夜激情福利司机影院| 日本一本二区三区精品| 天天一区二区日本电影三级| 欧美+亚洲+日韩+国产| 欧美黄色片欧美黄色片| 夜夜爽天天搞| 女人被狂操c到高潮| 国产精品野战在线观看| 亚洲成人免费电影在线观看| 99re在线观看精品视频| 国产伦人伦偷精品视频| 精品一区二区三区四区五区乱码| 人人澡人人妻人| 日本免费一区二区三区高清不卡| 在线国产一区二区在线| 午夜福利在线在线| 天天添夜夜摸| 成年人黄色毛片网站| 伦理电影免费视频| 女性被躁到高潮视频| 村上凉子中文字幕在线| 亚洲精品国产精品久久久不卡| 亚洲精品久久国产高清桃花| 婷婷丁香在线五月| 午夜日韩欧美国产| 一区二区三区高清视频在线| 亚洲成人国产一区在线观看| 黄片大片在线免费观看| 老司机福利观看| 熟妇人妻久久中文字幕3abv| av天堂在线播放| 女性生殖器流出的白浆| 两性夫妻黄色片| 人人澡人人妻人| 又紧又爽又黄一区二区| 美女午夜性视频免费| 免费观看人在逋| x7x7x7水蜜桃| 中文字幕av电影在线播放| 亚洲真实伦在线观看| 国产精品九九99| 久久香蕉激情| 中文字幕久久专区| 动漫黄色视频在线观看| 黄色女人牲交| 好男人在线观看高清免费视频 | 91成人精品电影| 国产精品自产拍在线观看55亚洲| 美女 人体艺术 gogo| 一区二区三区国产精品乱码| 婷婷精品国产亚洲av在线| 久久久久久久午夜电影| 国产精品美女特级片免费视频播放器 | 免费看日本二区| 人人妻人人看人人澡| 日韩欧美 国产精品| 国产av不卡久久| 午夜精品久久久久久毛片777| 国产精品一区二区免费欧美| 欧美三级亚洲精品| 亚洲精品美女久久av网站| 精品欧美国产一区二区三| 亚洲欧美日韩高清在线视频| 亚洲 欧美 日韩 在线 免费| 色综合亚洲欧美另类图片| 色播亚洲综合网| 看黄色毛片网站| 99国产极品粉嫩在线观看| 露出奶头的视频| 美女 人体艺术 gogo| 精品国产国语对白av| 女人高潮潮喷娇喘18禁视频| 亚洲九九香蕉| 免费看日本二区| 激情在线观看视频在线高清| 男人舔女人下体高潮全视频| 满18在线观看网站| 亚洲av电影不卡..在线观看| 777久久人妻少妇嫩草av网站| 在线永久观看黄色视频| 日本 欧美在线| 伦理电影免费视频| 黄网站色视频无遮挡免费观看| 1024香蕉在线观看| 欧美乱妇无乱码| 色综合婷婷激情| 无遮挡黄片免费观看| 亚洲精品中文字幕一二三四区| 亚洲av电影不卡..在线观看| 久久精品91蜜桃| 亚洲自偷自拍三级| 国产麻豆成人av免费视频| 国产成人一区二区在线| 五月伊人婷婷丁香| 深爱激情五月婷婷| 99热这里只有精品一区| 欧美区成人在线视频| 国产一区亚洲一区在线观看| 可以在线观看的亚洲视频| 日日摸夜夜添夜夜添小说| 人人妻人人澡欧美一区二区| 男女之事视频高清在线观看| 亚洲精品456在线播放app| 国产白丝娇喘喷水9色精品| 国产亚洲91精品色在线| 99久国产av精品| 国产精品乱码一区二三区的特点| 国内揄拍国产精品人妻在线| 一边摸一边抽搐一进一小说| 精品久久久久久久久久免费视频| 精品人妻偷拍中文字幕| 俄罗斯特黄特色一大片| 亚洲av不卡在线观看| 一级黄色大片毛片| 九色成人免费人妻av| 国产日本99.免费观看| 国产精品乱码一区二三区的特点| 97人妻精品一区二区三区麻豆| 亚洲精品影视一区二区三区av| 欧美国产日韩亚洲一区| 国产一级毛片七仙女欲春2| 九九热线精品视视频播放| 国产精品久久久久久久久免| 欧美最新免费一区二区三区| 久久国产乱子免费精品| 成年女人永久免费观看视频| 久久人人爽人人片av| 成人高潮视频无遮挡免费网站| 蜜臀久久99精品久久宅男| 免费av毛片视频| 国产淫片久久久久久久久| 国产乱人偷精品视频| .国产精品久久| 天堂√8在线中文| 51国产日韩欧美| 亚洲最大成人中文| 国产精品免费一区二区三区在线| 国产美女午夜福利| 亚洲专区国产一区二区| 亚州av有码| 国产日本99.免费观看| 网址你懂的国产日韩在线| 亚洲熟妇中文字幕五十中出| 在线免费观看不下载黄p国产| 少妇熟女aⅴ在线视频| 99热网站在线观看| 欧美xxxx性猛交bbbb| 色哟哟哟哟哟哟| 欧美日韩乱码在线| 日本a在线网址| 级片在线观看| 久久欧美精品欧美久久欧美| 岛国在线免费视频观看| 日韩欧美国产在线观看| 大型黄色视频在线免费观看| 日韩欧美一区二区三区在线观看| 亚洲va在线va天堂va国产| 亚洲人成网站在线观看播放| 男女下面进入的视频免费午夜| 最近手机中文字幕大全| 国产亚洲精品av在线| 少妇被粗大猛烈的视频| 春色校园在线视频观看| 在线播放无遮挡| 高清日韩中文字幕在线| 久久99热6这里只有精品| 亚洲成人av在线免费| 综合色av麻豆| 国产av一区在线观看免费| 日韩人妻高清精品专区| 一个人观看的视频www高清免费观看| 国产色婷婷99| 黄片wwwwww| 成人av在线播放网站| 少妇高潮的动态图| 国产成人福利小说| 欧美极品一区二区三区四区| 日韩欧美 国产精品| 69av精品久久久久久| 赤兔流量卡办理| 国产伦精品一区二区三区视频9| 久久精品夜色国产| 三级男女做爰猛烈吃奶摸视频| 亚洲无线在线观看| 午夜久久久久精精品| 久久久午夜欧美精品| 99riav亚洲国产免费| 欧美日韩国产亚洲二区| 97超视频在线观看视频| 免费看日本二区| 国产一级毛片七仙女欲春2| 午夜亚洲福利在线播放| 99热精品在线国产| 精品国产三级普通话版| a级毛色黄片| 国产极品精品免费视频能看的| 一夜夜www| 天堂av国产一区二区熟女人妻| 日韩精品有码人妻一区| 99热6这里只有精品| 高清午夜精品一区二区三区 | 美女黄网站色视频| 91精品国产九色| 22中文网久久字幕| 亚洲va在线va天堂va国产| 欧美激情国产日韩精品一区| 五月玫瑰六月丁香| 国产精品久久电影中文字幕| 亚洲成人精品中文字幕电影| 18禁裸乳无遮挡免费网站照片| 看黄色毛片网站| 亚洲av美国av| 欧美极品一区二区三区四区| 亚洲av中文av极速乱| 国内精品宾馆在线| 亚洲av五月六月丁香网| 深爱激情五月婷婷| 亚洲丝袜综合中文字幕| 国产高清激情床上av| 成人鲁丝片一二三区免费| 天天躁夜夜躁狠狠久久av| 最后的刺客免费高清国语| 伊人久久精品亚洲午夜| or卡值多少钱| 欧美一级a爱片免费观看看| 欧美日韩在线观看h| 亚洲人与动物交配视频| 亚洲第一区二区三区不卡| 免费看美女性在线毛片视频| 精品久久久噜噜| 亚洲欧美成人综合另类久久久 | 悠悠久久av| 午夜免费激情av| 国内精品一区二区在线观看| 亚洲性夜色夜夜综合| 亚洲丝袜综合中文字幕| 亚洲av.av天堂| 久久久久精品国产欧美久久久| 精品人妻视频免费看| 狂野欧美白嫩少妇大欣赏| 午夜福利在线观看免费完整高清在 | 国内精品美女久久久久久| 一区二区三区免费毛片| 国产探花极品一区二区| 亚洲精品国产成人久久av| 精品久久久久久久人妻蜜臀av| 亚洲18禁久久av| 中文字幕人妻熟人妻熟丝袜美| 国产不卡一卡二| 国产69精品久久久久777片| 日韩一本色道免费dvd| 老司机影院成人| 国产高清不卡午夜福利| 欧美激情在线99| 亚洲精华国产精华液的使用体验 | 九九久久精品国产亚洲av麻豆| 极品教师在线视频| 亚洲精品亚洲一区二区| 亚洲欧美中文字幕日韩二区| 亚洲成人av在线免费| 欧美xxxx黑人xx丫x性爽| 婷婷亚洲欧美| 两个人视频免费观看高清| 国产精品99久久久久久久久| 乱人视频在线观看| 日韩欧美在线乱码| 国产黄a三级三级三级人| 国产在线精品亚洲第一网站| 久久精品国产清高在天天线| 内地一区二区视频在线| 色噜噜av男人的天堂激情| 99久国产av精品| 一级a爱片免费观看的视频| 啦啦啦啦在线视频资源| 国产精品福利在线免费观看| a级毛色黄片| 婷婷精品国产亚洲av在线| 国产中年淑女户外野战色| 成人av一区二区三区在线看| 毛片女人毛片| 在线天堂最新版资源| 俺也久久电影网| 日本-黄色视频高清免费观看| 国产亚洲av嫩草精品影院| 欧美区成人在线视频| 成人无遮挡网站| 国产视频一区二区在线看| 高清毛片免费观看视频网站| 国产视频一区二区在线看| 俺也久久电影网| 夜夜看夜夜爽夜夜摸| 俺也久久电影网| 1000部很黄的大片| 国产精品一二三区在线看| 久久国产乱子免费精品| 色哟哟哟哟哟哟| 久久久国产成人免费| 精品一区二区三区人妻视频| 好男人在线观看高清免费视频| 1000部很黄的大片| 熟女人妻精品中文字幕| 欧美+亚洲+日韩+国产| 51国产日韩欧美| 蜜臀久久99精品久久宅男| 国产亚洲欧美98| 91av网一区二区| 九九爱精品视频在线观看| 九九久久精品国产亚洲av麻豆| 国产av一区在线观看免费|