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

    東亞地區(qū)夏季Es的分布特性

    2016-12-21 02:06:13杜鵬車海琴滿莉
    電波科學學報 2016年5期
    關(guān)鍵詞:中值山川電離層

    杜鵬 車海琴 滿莉

    (中國電波傳播研究所,青島 266107)

    ?

    東亞地區(qū)夏季Es的分布特性

    杜鵬 車海琴 滿莉

    (中國電波傳播研究所,青島 266107)

    東亞地區(qū)是全球Es出現(xiàn)的峰值區(qū)域.以往關(guān)于東亞異常的研究和討論具有一定的局限性,或者由于缺乏中國區(qū)域的資料而論據(jù)不足,或者缺乏中國、日本Es中心的詳細對比.文中利用中國及周邊地區(qū)23個垂測站二十余年的Es數(shù)據(jù)進行空間插值分析,給出東亞地區(qū)夏季Es出現(xiàn)率的空間分布等值線,并深入分析東亞異常中心區(qū)域的時間變化特征.研究表明,中國重慶、廣州或者日本山川、國分寺、秋田是夏季Es的高發(fā)區(qū)域,東亞異常中心始終位于其中某個地區(qū),具體位置隨著統(tǒng)計指標改變;高發(fā)區(qū)域內(nèi)的Es強度呈現(xiàn)出年變化和日變化,中國高發(fā)區(qū)域Es強度的晝夜差別比日本高發(fā)區(qū)域更為顯著,但后者具有更為復雜的年變化和晝夜變化特征,這對揭示中緯Es形成原因具有一定的意義.

    突發(fā)E層(Es);東亞異常;出現(xiàn)率;高發(fā)區(qū)域

    DOI 10.13443/j.cjrs.2016032801

    引 言

    突發(fā)E層(Sporadic E,ES)是電離層內(nèi)較常見的一種隨機異?,F(xiàn)象,能夠顯著改變無線電波的傳輸路徑,對短波通信、衛(wèi)星導航、雷達偵測等業(yè)務(wù)造成影響.根據(jù)火箭[1-3]、非相干散射雷達[4-5]、中高層大氣雷達和激光雷達[6-7]等觀察結(jié)果:Es是電離增強的不均勻薄層,具有極高的電子密度和梯度,一般位于90~120 km高度范圍內(nèi),呈現(xiàn)出層狀、片狀或塊狀空間結(jié)構(gòu),水平尺度可達幾十甚至上百千米,但垂直方向厚度約1~3 km.

    Es一直是電離層研究的熱點和難點,人們利用多種探測手段對Es形態(tài)結(jié)構(gòu)進行了詳細研究,并對其形成機制給出一定的理論解釋.通過全球電離層垂測站記錄的統(tǒng)計分析,Smith[8]早在1957年基于國際地球物理年(International Geophysical Year,IGY)期間的Es觀測數(shù)據(jù),就發(fā)現(xiàn)了電離層東亞異常和夏季異常,但迄今為止這些特性還不能完全得到解釋.Whitehea[9-10]發(fā)展并完善了風剪切作用壓縮金屬離子形成中緯Es的理論,并指出Es與地磁活動、流星雨缺乏關(guān)聯(lián),雖然后者是Es層金屬離子的主要來源.激光雷達觀測進一步驗證了突發(fā)金屬原子層(如突發(fā)納層和鉀層)與Es之間存在一定聯(lián)系.Mathews[11]試圖結(jié)合風剪切、外部電場和潮汐風系統(tǒng)來解釋中緯Es,并指出對流電場和等離子體不穩(wěn)定性分別在極區(qū)Es和赤道Es形成過程中起主導作用.在突發(fā)E層散射(SESCAT)試驗中,Voiculescu和Haldoupis等人[12]首先發(fā)現(xiàn)行星波對甚高頻(Very High Frequency,VHF)雷達后向散射回波和中緯Es出現(xiàn)率的作用.Pancheva和Haldoupis等人[13-14]給出了行星波與Es相關(guān)性的直接證據(jù),并進一步通過分析中間層和低熱層(Mesosphere and Lower Thermosphere,MLT)風場的潮汐和半日潮汐幅度擾動的周期,得出行星波通過調(diào)制潮汐間接影響Es的觀點.Zuo and Wan[15]通過研究武漢地區(qū)Es出現(xiàn)率的準6日、12小時和24小時擾動與行星波之間的相關(guān)性,為行星波通過調(diào)制潮汐間接作用于Es的觀點提供了新證據(jù).Kagan和Kelley通過激光雷達測量首次利用光學方法獲得可視化中緯Es二維水平結(jié)構(gòu)[16-17],并具體討論了生成Es不均勻體的不穩(wěn)定性機制[18].

    Es地域形態(tài)統(tǒng)計特征很早就引起了學者的關(guān)注.北半球中緯地區(qū)是Es高發(fā)區(qū),國際研究早已給出東亞異常的定論,即東亞中低緯地區(qū)夏季Es的出現(xiàn)率特別高,Es強度也特別強,遠遠高于全球其它同緯度地區(qū).遺憾的是,歷史上Smith得出這一結(jié)論時,全球垂測站數(shù)量較少,尤其缺乏中國地區(qū)Es觀測數(shù)據(jù)的支持.隨著衛(wèi)星技術(shù)的發(fā)展,利用衛(wèi)星信標探測Es空間分布的方法逐漸成熟,如2005年Wu和Ao等人[19]利用GPS-CHAMP(CHAMP 衛(wèi)星)任務(wù)期間對Es進行掩星探測,給出全球范圍內(nèi)的Es形態(tài)分布以及時變特征,指出Es對磁傾角有明顯的依賴性,其研究成果有效地豐富了全球Es形態(tài)特征.具體的,在北半球磁緯度為30°至60°處夏季(6-8月)中緯度Es出現(xiàn)最為頻繁,除了東經(jīng)90°至130°較大范圍(即東亞異常區(qū)域)外,在磁緯度為30°至60°、東經(jīng)175°至185°附近,以及南半球夏季(10-12月)磁緯度-55°至-35°、東經(jīng)115°至125°附近和磁緯-60°至-55°,東經(jīng)175°至185°附近也存在著Es出現(xiàn)率非常高的區(qū)域.

    衛(wèi)星信標探測技術(shù)在一度程度上拓展了Es研究手段,但也有自身缺陷,例如缺少全天候探測能力,探測點分布不均.相比較而言,全球范圍的電離層垂測站數(shù)量越來越多,全時段觀測數(shù)據(jù)的積累時間越來越長.雖然Es截止頻率(foEs)與垂測設(shè)備的靈敏度有關(guān),但Smith根據(jù)垂測數(shù)據(jù)分析給出結(jié)論后,又通過VHF斜向探測資料加以驗證.因此,利用電離層垂測數(shù)據(jù)開展Es形態(tài)學研究仍然具有正確性和明顯優(yōu)勢,如Pietrella和Bianchi[20]利用32年觀測數(shù)據(jù)對羅馬Es出現(xiàn)率進行了研究,指出四季中日出前和日出后Es出現(xiàn)率較低,夏季白天的Es出現(xiàn)率最大.

    我國位于東亞異常區(qū)域,對Es形態(tài)學研究而言是需要重點關(guān)注的區(qū)域.1990年李玉輝[21]利用國內(nèi)1964-1968年9個電離層垂測站的歷史數(shù)據(jù),首次給出我國夏季(5-8月)foEs>7 MHz的Es出現(xiàn)率分布,對國際無線電咨詢委員(International Radio Consultative Committee,CCIR)報告中全球Es出現(xiàn)率分布進行了修正.隨后,李鈞[22]對東亞異常中心提出了異議,認為其中心位于我國的蘭州和重慶,而非日本東南部.譚輝[23]對重慶和蘭州foEs>5 MHz的Es出現(xiàn)率進行了詳細對比,認為兩地Es顯著的差異性除了緯度差別外還與蘭州位于地磁Sq電流焦點有關(guān).左小敏[24]利用武漢測量數(shù)據(jù)探討了中緯Es出現(xiàn)率與潮汐和行星波之間的關(guān)系,趙海生和許正文等人[25]利用中國和日本觀測臺站的數(shù)據(jù)對Es的強度特性和空時分布特點進行了探討,獲得了十分有意義的結(jié)果.

    盡管如此,上述相關(guān)研究或者只針對中國地區(qū)的數(shù)據(jù),或者只關(guān)注中國臺站與日本臺站數(shù)據(jù)的對比,并未涵蓋完整的中國周邊地區(qū).本文利用東亞地區(qū)多個臺站二十余年的歷史觀測數(shù)據(jù)進行分析,獲得東亞地區(qū)夏季Es出現(xiàn)率等值分布圖,并對東亞異常中心位置及中心區(qū)域的Es時間變化特性進行深入探討,所述結(jié)果有效地完善了東亞尤其中國地區(qū)的Es出現(xiàn)率的分布特性,對揭示中緯Es的形成機制也有一定的積極意義.

    1 觀測數(shù)據(jù)資料及分析方法

    利用測高儀從地面垂直向電離層發(fā)射無線電信號是最傳統(tǒng)、使用最廣泛的電離層探測方法,目前全世界已建成兩百多個電離層觀測站,并實現(xiàn)大部分觀測數(shù)據(jù)(以1 h為間隔)的國際交換.經(jīng)過多年發(fā)展,我國境內(nèi)已形成了布局合理的電離層垂直探測網(wǎng).以往觀測站主要分布在中東部,西部地區(qū)觀測站點較少,雖然近年來加大了西部觀測站點建設(shè),但這些臺站積累的數(shù)據(jù)仍然有限.為保證Es統(tǒng)計樣本的豐富度,這里仍選用國內(nèi)歷史數(shù)據(jù)最豐富的臺站,并從美國地球物理中心公開網(wǎng)站下載我國周邊國家地區(qū)主要觀測站的Es數(shù)據(jù).利用這些數(shù)據(jù)作為中國觀測數(shù)據(jù)的外圍邊界,提高空間插值分析的精度和可信度,由此獲得東亞地區(qū)Es出現(xiàn)率的空間分布.

    考慮樣本的時間重合度,選擇1968-1988年夏季(5-8月份)烏魯木齊、滿洲里、長春、北京、蘭州、重慶、武漢、廣州、海口、拉薩、秋田、國分寺、沖繩、稚內(nèi)、山川、阿拉木圖、卡拉干達、新卡扎林斯克、阿什哈巴德、伊爾庫茨克、新西伯利亞、斯維爾德洛夫斯克、托木斯克共23個臺站的foEs數(shù)據(jù).各臺站位置如表1所示.

    表1 1968-1988年夏季Es觀測臺站及其位置

    描述Es的最重要參量為foEs,衡量Es強度的重要標準是一段時間內(nèi)Es的出現(xiàn)率和中值.在Es空間形態(tài)學統(tǒng)計中,慣常作法[8,20-21]是考察一段時間內(nèi)foEs大于某一特定頻率fT的出現(xiàn)率(即大于fT的有效數(shù)據(jù)除以總統(tǒng)計樣本,一般選取fT為5、7 MHz),記為P(foEs>fT).

    文中取P(foEs>5 MHz)和P(foEs>7 MHz)作為Es出現(xiàn)率的統(tǒng)計指標,采用克里金插值方法[26-27]對東亞地區(qū)的foEs數(shù)據(jù)進行統(tǒng)計分析,獲得Es出現(xiàn)率的空間分布圖,然后詳細分析高發(fā)區(qū)域內(nèi)Es強度時間變化特性.

    2 東亞地區(qū)Es出現(xiàn)率的空間分布

    電離層Es存在顯著的晝夜日變化,這里分別考察1968-1988年夏季全天、白天(07:00-18:00)和夜晚(19:00-06:00)的Es出現(xiàn)率的空間分布特征,圖1~圖3給出東亞地區(qū)相應(yīng)時段的空間等值線分布圖.

    從圖1~圖3可以看出,日本境內(nèi)的山川、國分寺和秋田以及我國西南的重慶均為夏季Es高發(fā)地區(qū),Es出現(xiàn)率明顯高于東亞其它相鄰緯度地區(qū).這幾個臺站均位于磁緯34°至55°、東經(jīng)100°至140°范圍內(nèi),與Wu和Ao等人[19]利用掩星觀測的結(jié)果相吻合.在不同時段下(全天、白天和夜晚),東亞異常區(qū)域的中心會發(fā)生偏移,夏季白天的Es出現(xiàn)率明顯高于夜間.其中,P(foEs>5 MHz)的白天最大值為58.7%,出現(xiàn)在重慶,夜間最大值為27.9%,出現(xiàn)在國分寺;P(foEs>7 MHz)的白天最大值為26.4%,出現(xiàn)在山川,夜間最大值為11.6%,出現(xiàn)在秋田.

    (a) 全天P(foEs>5 MHz) (b) 全天P(foEs>7 MHz)圖1 1968-1988年5-8月,東亞全天Es出現(xiàn)率分布

    (a) 白天P(foEs>5 MHz) (b) 白天P(foEs>7 MHz)圖2 1968-1988年5-8月,東亞白天Es出現(xiàn)率分布

    (a) 夜晚P(foEs>5 MHz) (b) 夜晚P(foEs>7 MHz)圖3 1968-1988年5-8月,東亞夜晚Es出現(xiàn)率分布

    表2 歷年夏季東亞異常區(qū)域中心位置及其Es出現(xiàn)率

    表2給出每年夏季全天統(tǒng)計條件下,東亞異常中心及相應(yīng)的Es出現(xiàn)率.對于全天結(jié)果而言,考慮foEs>5 MHz的Es出現(xiàn)率時,21年間東亞異常中心2/3的年份出現(xiàn)在日本境內(nèi),其中國分寺7次、山川4次、秋田3次,1/3的年份出現(xiàn)在中國境內(nèi),其中重慶5次、廣州2次;考慮foEs>7 MHz的Es出現(xiàn)率時,東亞異常中心位于日本境內(nèi)的概率高達80%,其中山川最為頻繁共出現(xiàn)10次.

    如果考慮電離層Es的日夜變化,具體統(tǒng)計值可能略有變化,這里不再贅述.

    3 東亞異常中心Es強度的時間變化

    由表2可知,東亞異常中心出現(xiàn)在我國重慶、廣州或者日本山川、國分寺、秋田一帶.為進一步研究東亞Es高發(fā)地區(qū)的Es變化規(guī)律,并便于中、日相鄰緯度臺站數(shù)據(jù)的對比,圖4給出有關(guān)臺站夏季月中值的日變化平均曲線.就月份而言,夏季Es最強出現(xiàn)在6月,其次是7月,5月和8月的Es稍弱.就一天中不同時段而言,09:00-15:00期間Es強度較高,最大值出現(xiàn)在當?shù)貢r10:00-11:00左右;22:00-06:00期間Es強度較弱,最小值出現(xiàn)在日出前(當?shù)貢r05:00附近).

    不同臺站的Es強度的具體日變化略有不同,相對日本山川、國分寺、秋田站而言,中國重慶、廣州站Es強度的晝夜變化較大,差異性也更為顯著,但日本三站的Es日變化更為復雜.具體的,在Es較強的09:00-15:00期間,處于相鄰緯度的重慶和山川的foEs中值非常接近,強于國分寺、秋田和廣州;在Es較弱的22:00-06:00期間,重慶和廣州的foEs中值低于日本的山川和國分寺;在17:00-21:00期間,重慶和廣州的Es強度隨時間一直減弱,但日本山川、國分寺、秋田的Es強度會先上升后下降,在某一時間出現(xiàn)較為明顯的亞峰,呈現(xiàn)出一定的準半日變化特性.

    這里以緯度相鄰的重慶和山川為例,對比分析Es強度的日變化和半日變化特性.為了避免Es突發(fā)性引起的數(shù)據(jù)缺失,選取1968年7月和1977年6月的foEs小時測量序列,分別作傅里葉變化和Morlet連續(xù)小波變換[14,24].限于篇幅,這里只給出1968年7月的頻譜圖和時間-周期圖,如圖5所示.可以看到,重慶和山川的日變化(24小時周期)分量最顯著,重慶站的半日變化(12小時周期)分量要稍弱于山川站.1977年6月份所述兩站的半日變化振幅較大,相對日變化來說較為明顯;高發(fā)地區(qū)其它臺站夏季月份foEs時間序列的頻譜也有類似特征,不再展開討論.

    (a) 5月份 (b) 6月份

    (a) 7月份 (d) 8月份圖4 Es高發(fā)地區(qū)foEs月中值的日變化

    (a) 頻譜圖(重慶) (b) 時間周期圖(重慶)

    (c) 頻譜圖(山川) (d) 時間周期圖(山川)圖5 1968年7月重慶和山川foEs頻譜和時間-周期圖

    為詳細討論東亞異常中心Es強度的逐年變化特性,圖6給出1968-1988年夏季,重慶、山川、蘭州、國分寺、廣州和秋田6個臺站夏季foEs中值的時-年變化特性.其中,重慶與山川、蘭州與國分寺的緯度接近,除蘭州站缺乏1973年夏季數(shù)據(jù)外,其余臺站在此期間的數(shù)據(jù)都是完整的.

    在前面第2節(jié)中提到,東亞地區(qū)夏季白天(07:00-18:00)的Es出現(xiàn)率明顯高于夜間(19:00-06:00)時段.從圖6中東亞Es高發(fā)區(qū)域foEs中值的時-年變化特性圖中還可以看到,不同時段的Es出現(xiàn)率也呈現(xiàn)出不同的特點,表3給出1968-1988年夏季白天和夜間不同時段Es高發(fā)地區(qū)P(foEs>5 MHz)的出現(xiàn)率.其中上午時段為08:00-12:00,下午時段為13:00-17:00,午前時段為10:00-12:00,午后時段為13:00-15:00,19:00-次日00:00為日落后至子夜時段即前半夜,01:00-05:00為子夜至日出前時段即后半夜.

    就各臺站自身而言,午前時段的Es出現(xiàn)率稍高于午后時段,前半夜的Es出現(xiàn)率明顯高于后半夜.如表3和圖6所示,就不同臺站對比而言,重慶和廣州午后時段的Es強度要略高于日本三站,但日本三站前半夜時段的Es強度高于中國三站,這進一步體現(xiàn)了重慶、廣州地區(qū)Es強度顯著的晝夜差異性.

    表3 夏季Es高發(fā)區(qū)不同時段P(foEs>5 MHz)出現(xiàn)率(%)

    從圖6還可以看到,就Es強度整體最強并且緯度相鄰的重慶和山川來說,1968-1975年09:00-15:00期間重慶foEs中值大大高于山川,1976-1988年期間剛好相反.與各臺站自身相比,1975年之前中國重慶、廣州地區(qū)白天時段的foEs中值要顯著高于1975年之后的結(jié)果,這一時期日本山川、國分寺、秋田地區(qū)白天的foEs中值要低于后期的結(jié)果.

    值得注意的是,電離層變化受到太陽活動的控制,上述期間的歷史foEs觀測數(shù)據(jù)完整、采樣間隔保持不變,1975年之前重慶和廣州foEs中值顯著高于1975年之后的結(jié)果很可能與太陽活動周期有密切關(guān)系.在太陽活動準11年周期、Hale周期(準22年)以及Gleissberg周期等共同作用下,第19周出現(xiàn)罕見的太陽活動(1957年為19周太陽活動高年,太陽黑子數(shù)達到有記錄以來的歷史峰值),這在一定程度上影響著foEs的長期變化趨勢,這與文獻[25]中重慶和廣州站foEs月中值出現(xiàn)長期下降趨勢吻合.考察Es高發(fā)地區(qū)1968-1988年foEs24小時統(tǒng)計中值和全天foEs>5 MHz出現(xiàn)率的全年變化,可以看到中日六站的Es強度均呈現(xiàn)出與太陽活動一致的準11年周期.對Es強度的太陽周期變化和長期趨勢變化將在今后深入探討.

    圖6 Es高發(fā)地區(qū)夏季foEs統(tǒng)計中值的時-年變化

    4 結(jié)論和討論

    利用中國及周邊地區(qū)23個觀測臺站21年5-8月的歷史觀測數(shù)據(jù),研究了東亞地區(qū)夏季Es出現(xiàn)率的空間分布,分析了東亞異常中心區(qū)域Es的時間分布,可以得出如下結(jié)論:

    1) 東亞地區(qū)作為夏季Es的高發(fā)區(qū)域,Es出現(xiàn)率分布在時間和空間上呈現(xiàn)出復雜的不均勻特性,其峰值區(qū)域位于磁緯34°至55°、東經(jīng)100°至140°范圍內(nèi),可能出現(xiàn)在我國重慶、廣州或者日本山川、國分寺、秋田中的某個地區(qū),具體的出現(xiàn)位置和出現(xiàn)率峰值與測量年份、統(tǒng)計時段和統(tǒng)計指標有關(guān),其中位于磁緯度45°至50°之間的山川、重慶和國分寺是Es峰值出現(xiàn)最為頻繁的地區(qū).

    若改變某個統(tǒng)計參數(shù),Es出現(xiàn)率值會相應(yīng)改變,峰值位置也可能會在上述地區(qū)中變化.如1968-1988年夏季foEs>5 MHz的全天最大出現(xiàn)率為39.9%,出現(xiàn)在國分寺,白天最大出現(xiàn)率為58.7%,出現(xiàn)在重慶;foEs>7 MHz的全天和白天最大出現(xiàn)率分別為18.9%和26.4%,均出現(xiàn)在山川.以每年全天24小時foEs>5 MHz的出現(xiàn)率作為統(tǒng)計指標時,東亞異常中心有1/3的年份出現(xiàn)在中國境內(nèi),其中重慶出現(xiàn)5次;2/3的年份出現(xiàn)在日本境內(nèi),其中國分寺和山川共出現(xiàn)11次.以全天24小時foEs>7 MHz的出現(xiàn)率作為統(tǒng)計指標時,東亞異常中心位于日本境內(nèi)的概率高達80%,僅山川就出現(xiàn)10次.

    2) 東亞Es高發(fā)地區(qū)夏季Es強度呈現(xiàn)出明顯的晝夜日變化特性和一定的半日變化特征.

    其中,夜間Es強度顯著低于白天,白天09:00-15:00為Es較強時段,夜間22:00-06:00為Es較弱時段,其中午前時段(10:00-12:00)的Es出現(xiàn)率稍高于午后時段(13:00-15:00),前半夜的Es出現(xiàn)率高于后半夜.

    中國高發(fā)地區(qū)(重慶、廣州)夏季Es日變化比日本高發(fā)地區(qū)(山川、國分寺和秋田)顯著,但日本高發(fā)地區(qū)的夏季Es半日變化特征相對來說更為明顯.

    3) 中、日高發(fā)地區(qū)夏季Es強度具有不同的逐年變化.1975年前Es強度較強時段,重慶Es中值遠高于同時期日本臺站以及1975年后重慶的觀測結(jié)果;1975年以后這一差距縮小,甚至在某些年份,重慶、廣州Es中值低于該時段日本臺站的觀測結(jié)果.

    這里給出的Es出現(xiàn)概率統(tǒng)計結(jié)果體現(xiàn)了Es空間分布形態(tài)的豐富多樣性,不僅在一定程度上完善了東亞尤其中國地區(qū)的Es出現(xiàn)率的分布特性,也有助于理解中緯Es的形成機制.目前,風剪切作用壓縮金屬離子形成中緯Es的物理機制已經(jīng)被大多數(shù)學者所接受,該理論預(yù)測壓縮機制在地磁場水平分量較大的地方起作用.對foEs>5 MHz全天出現(xiàn)概率與國際地磁參考場(International Geomagnetic Reference Field,IGRF)模型給出的地磁場水平分量作相關(guān)分析,可以發(fā)現(xiàn)所述21年夏季期間二者間的相關(guān)系數(shù)在0.66~0.88之間,具有較明顯的正相關(guān)性.

    日變化和半日變化是Es準周期性結(jié)構(gòu)的一種表現(xiàn)形式,也是Es研究的熱點之一.引言中已經(jīng)指出,已經(jīng)有大量研究通過單個不同地區(qū)Es和同期風場觀測數(shù)據(jù)之間的比較,分析Es數(shù)據(jù)中的12小時和24小時潮汐分量,給出大氣潮汐風場引起Es半日變化和日變化的證據(jù).通過多個地區(qū)的GPS掩星觀測資料,發(fā)現(xiàn)夏季中緯度Es確實具有較強的日變化和半日變化特性.本文初步給出了夏季Es高發(fā)地區(qū)不同臺站Es的日變化和半日變化差異性,可以對Es形成原因的研究提供某些支持.此外,了解東亞Es地區(qū)的Es長周期性變化特征有利于建立Es與其控制及影響因素之間的關(guān)聯(lián),這將是今后進一步研究的重點.

    致謝:感謝中國電波傳播研究所、中國科學院地質(zhì)與地球物理研究所以及美國地球物理中心提供電離層Es觀測數(shù)據(jù).

    [1]ISTOMIN V G.Absolute concentration of ion components of the earth’s atmosphere at altitudes between 100 and 200 km[J].Planetary and space science, 1963, 11(2):169-172.

    [2]BERNHARDT P, et al.Multi-layer structure of mid-latitude sporadic-E observed during the SEEK-2 campaign[J].Annales geophysicae, 2005, 23(7):2347-2355.

    [3]OYAMA K I, HIBINO K, et al.Energetics and structure of the lower E region associated with sporadic E layer[J].Annales geophysicae, 2008, 26(9):2929-2936.

    [4]SWARTZ W E, Collins S C, et al.First observations of an F-region turbulent upwelling coincident with severe E-region plasma and neutral perturbations[J].Journal of atmospheric and solar-terrestrial physics, 2002:1545-1556.

    [5]DAMTIE B, et al.High resolution observations of sporadic-E layers within the polar cap ionosphere using a new incoherent scatter radar experiment[J].Annales geophysicae, 2002,20(9):1429-1438.

    [6]KANE T, GRIME B, et al.Joint observations of sodium enhancements and field-aligned ionospheric irregularities[J].Geophysical research letters, 2001, 28(28):1375-1378.

    [7]LARSEN M F, HYSELL D L, et al.Imaging coherent scatter radar, incoherent scatter radar, and optical observations of quasiperiodic structures associated with sporadic E layers[J].Journal of geophysical research atmospheres, 2007, 112(A6):23-27.

    [8]SMITH J.Ionospheric sporadic E[M]//International series of monographs on electromagnetic waves:Volume 2[M].Britain:Pregamon Press, 1962:3-12.

    [9]WHITEHEAD J D.The formation of the sporadic-E layer in the temperate zones[J].Journal of atmospheric and terrestrial physics, 1960, 188(188):567.

    [10]WHITEHEAD J D.Recent work on mid-latittude and equatourial sporadic E[J].Journal of atmospheric and terrestrial physics, 1989, 51(5):401-424.

    [11]MATHEWS J D.Sporadic E:current views and recent progress[J].Journal of atmospheric and solar-terrestrial physics, 1998, 60(4):413-435.

    [12]VOICULESCU M, HALDOUPIS C,SCHLEGEL K.Evidence for planetary wave effects on midlatitude backscatter and sporadic E layer occurrence[J].Geophysical research letters, 1999, 26(8):1105-1108.

    [13]HALDOUPIS C, PANCHEVA D.Planetary waves and midlatitude sporadic E layers:Strong experimental evidence for a close relationship[J].Journal of geophysical research atmospheres, 2002, 107(A6).DOI:10.1029/2001JA000212.

    [14]PANCHEVA D, HALDOUPIS C, et al.Evidence for a role for modulation atmospheric tides in the dependence of sporadic E layers on planetary waves[J].Journal of geophysical research atmospheres, 2003, 108(A5).

    [15]ZUO X M, WAN W X.Planetary wave oscillations in sporadic E layer occurrence at Wuhan[J].Earth planets and space, 2008, 60(6):647-652.

    [16]KAGAN L M, KELLEY M C.A wind-driven gradient drift mechanism for mid-latitude E-region ionospheric irregularities[J].Geophysical research letters, 1998, 25(22):4141-4144.

    [17]KAGAN L M, KELLEY M C.A thermal mechanism for generation of small-scale irregularities in the ionospheric E region[J].Journal of geophysical research space physics, 2000, 105(A3):5291-5303.

    [18]KAGAN L M, KELLEY M C, et al.The structure of electromagnetic wav-induced 557.7-nm emission associated with a sporadic-E event over Arecibo[J].Physical review letters, 2000, 85(1):218-221.

    [19]Wu D L, Ao C O, et al.Sporadic E morphology from GPS-CHAMP radio occultation[J].Journal of geoph-ysical research, 2005, 110:A01306.DOI:10.1029/2004JA010701.

    [20]PIETRELLA M, BIANCHI C.Occurrence of sporadic-E layer over the ionospheric station of Rome:Analysis of data for thirty-two years[J].Advances in space research, 2009, 44(1):72-81.

    [21]李玉輝.中國夏季Es出現(xiàn)概率的分布(foEs>7 MHz) [J].空間科學學報, 1990, 10(2):156-158.

    LI Y H.Distribution of Es occurrence probability in the summer over China (foEs>7 MHz)[J].Chinese journal of space science, 1990, 10(2):156-158.(in Chinese)

    [22]李鈞.電離層我國地區(qū)特性研究的進展[C]//第四屆全國電波傳播學術(shù)會議論文集(上).武漢:1991:3-12.

    LI J.The development of the research on the characteristics of the ionosphere in China[C]// Proceedings of the forth national conference on radio propagation.Wuhan:1991:3-12.

    [23]譚輝.重慶與蘭州電離層Es出現(xiàn)率的比較研究[J].空間科學學報, 2000, 20(4):373-379.

    TAN H.Comparison and investigation of Es occurrence of the ionosphere over Chongqing and Lanzhou[J].Chinese journal of space science, 2000, 20(4):373-379.(in Chinese)

    [24]左小敏.中緯電離層突發(fā)E層中的波動信號及長期變化特性[D].武漢:中國科學院研究生院,2008.

    ZUO X M.Wave signatures and long-term trend seen through mid-latitude sporadic E-layers [D].Wuhan:University of Chinese Academy of Sciences, 2008.(in Chinese)

    [25]趙海生, 許正文, 吳振森, 等.東亞地區(qū)Es特性研究[J].空間科學學報, 2014, 24(1):38-45.

    ZHAO H S, XU Z W, WU Z S, et al.Sporadic E morphology of East Asia[J].Chinese journal of space science, 2014, 24(1):38-45.(in Chinese)

    [26]王世凱, 焦培南, 柳文.改進的Kriging 技術(shù)實時重構(gòu)區(qū)域電離層foF2的分布[J].電波科學學報, 2006, 21(2):166-171.

    WANG S K, JIAO P N, LIU W.Improved Kriging technique of ionospheric parameter foF2 instantaneous mapping[J].Chinese journal of radio science, 2006, 21(2):166-171.(in Chinese)

    [27]陳歡歡, 李星, 丁文秀.Surfer8.0等值線繪制中的十二種插值方法[J].工程地球物理學報, 2007, 4(1):52.

    CHEN H H, LI X, DING W X.Twelve kinds of gridding methods of surfer 8.0 in isoline drawing[J].Chinese journal of engineering geophysics, 2007, 4(1):52.(in Chinese)

    杜鵬 (1975-),男,山東人,1997年畢業(yè)于桂林電子工業(yè)學院無線電工程系,2015年結(jié)業(yè)于西安電子科技大學電子工程學院工程碩士班,現(xiàn)為中國電波傳播研究所高級工程師.主要從事電離層的觀測與研究,數(shù)據(jù)信號的測試、分析處理和數(shù)據(jù)庫建立等工作.

    車海琴 (1982-),女,湖北人,2005年畢業(yè)于南開大學信息技術(shù)科學學院,2008年畢業(yè)于中國電波傳播研究所無線電物理專業(yè),現(xiàn)為中國電波傳播研究所工程師.主要從事電離層物理、等離子體技術(shù)等工作.

    滿莉 (1984-),女,山東人,2005年畢業(yè)于大連交通大學,2010年畢業(yè)于中國海洋大學工程學院,現(xiàn)為中國電波傳播研究所工程師.主要從事電離層數(shù)據(jù)的分析與研究,數(shù)據(jù)庫建立等工作.

    Distribution characteristics of the summertime sporadic E in East Asia

    DU Peng CHE Haiqin MAN Li

    (ChinaResearchInstituteofRadiowavePropagation,Qingdao266107,China)

    The occurrence of global sporadic E reaches its peak in East Asia, known as the East Asia Anomaly.The relevant investigations and discussions are previously limited by either observational data unavailable in China or short of detailed comparisons of results in the popular areas of China with that in Japan, both of which are central areas of global Es occurrence probabilities.Based on spatial interpolation of Es data in more than twenty years from twenty-three vertical sounding stations in China and the surrounding countries,the results are drawn in the contour lines of spatial distribution of summer Es occurrence in East Asia, following a thorough analysis on characteristics of time variations in the center of the East Asia abnormal region.It is shown that Chongqing and Guangzhou in China together with Yamagawa, Kokubunji and Akita in Japan have remarkable Es occurrence all the time in summer, some of which is the potential peak of the East Asia abnormal region.The peak’s specific location varies with different statistical index.The Es intensities in the abnormal region center are characterized with annual and diurnal variations, which are more complicated in the popular areas of Japan than in that of China, while the diurnal differences of the Es strength are more profound in the latter region, which may be of great significance to understand the cause of mid-latitude Es.

    sporadic E layer (Es);the East Asia Anomaly;occurrence probability;popular area

    杜鵬, 車海琴, 滿莉.東亞地區(qū)夏季Es的分布特性[J].電波科學學報,2016,31(5):969-977.

    10.13443/j.cjrs.2016032801

    DU P, CHE H Q, MAN L.Distribution characteristics of the summertime sporadic E in East Asia [J].Chinese journal of radio science,2016,31(5):969-977.(in Chinese).DOI:10.13443/j.cjrs.2016032801

    2016-03-28

    電波環(huán)境特性及?;夹g(shù)國家重點實驗室基金

    P352

    A

    1005-0388(2016)05-0969-09

    聯(lián)系人:滿莉 E-mail:manlionly@qq.com

    猜你喜歡
    中值山川電離層
    一種電離層TEC格點預(yù)測模型
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    山川依舊 風味不改
    退耕還林還草 再造秀美山川
    綠色中國(2019年19期)2019-11-26 07:13:18
    退耕還林還草 綠了山川 富了百姓
    綠色中國(2019年17期)2019-11-26 07:04:40
    山川
    寶藏(2018年1期)2018-04-18 07:39:48
    Lagrange中值定理的巧妙應(yīng)用
    微分中值定理教法研討
    電離層對中高軌SAR影響機理研究
    雷達學報(2017年6期)2017-03-26 07:53:02
    后中值波電流脈沖MIG焊工藝
    老司机在亚洲福利影院| 黄色一级大片看看| 丝袜人妻中文字幕| 国产高清国产精品国产三级| 老司机在亚洲福利影院| 巨乳人妻的诱惑在线观看| 伦理电影大哥的女人| 十八禁人妻一区二区| 一本—道久久a久久精品蜜桃钙片| 高清不卡的av网站| 天堂俺去俺来也www色官网| 亚洲精品久久成人aⅴ小说| 精品国产一区二区三区久久久樱花| 欧美日韩视频精品一区| 9色porny在线观看| 如何舔出高潮| 日韩,欧美,国产一区二区三区| 国产精品国产三级专区第一集| av又黄又爽大尺度在线免费看| 国产免费一区二区三区四区乱码| 97在线人人人人妻| 狠狠婷婷综合久久久久久88av| 操出白浆在线播放| av卡一久久| 欧美成人精品欧美一级黄| 1024香蕉在线观看| 成人黄色视频免费在线看| 国产亚洲午夜精品一区二区久久| 欧美日韩精品网址| 久久久精品免费免费高清| 欧美日韩综合久久久久久| 精品国产乱码久久久久久小说| 久久久精品94久久精品| 国产精品国产三级国产专区5o| 久久婷婷青草| 电影成人av| 欧美国产精品一级二级三级| 老司机在亚洲福利影院| 考比视频在线观看| 如日韩欧美国产精品一区二区三区| av一本久久久久| 丝袜脚勾引网站| 国产色婷婷99| 人成视频在线观看免费观看| 欧美少妇被猛烈插入视频| 美女国产高潮福利片在线看| 欧美黑人精品巨大| av在线观看视频网站免费| 伦理电影免费视频| 黄网站色视频无遮挡免费观看| 赤兔流量卡办理| av网站免费在线观看视频| 亚洲国产欧美在线一区| 99久国产av精品国产电影| 亚洲精品一区蜜桃| 国产成人精品久久久久久| 成人国产av品久久久| 黄片小视频在线播放| 亚洲激情五月婷婷啪啪| 亚洲欧洲精品一区二区精品久久久 | 亚洲av成人不卡在线观看播放网 | 一区二区日韩欧美中文字幕| 久久99精品国语久久久| 色视频在线一区二区三区| 蜜桃国产av成人99| 午夜福利在线免费观看网站| 国产精品人妻久久久影院| 亚洲精品国产av成人精品| 国产精品麻豆人妻色哟哟久久| 人体艺术视频欧美日本| 欧美日韩亚洲国产一区二区在线观看 | 久久韩国三级中文字幕| 国产在线一区二区三区精| 人人妻人人澡人人爽人人夜夜| 亚洲天堂av无毛| 色视频在线一区二区三区| 如何舔出高潮| av.在线天堂| 少妇的丰满在线观看| 欧美 亚洲 国产 日韩一| 欧美av亚洲av综合av国产av | 午夜激情av网站| 久久国产精品男人的天堂亚洲| 51午夜福利影视在线观看| 国产97色在线日韩免费| 免费日韩欧美在线观看| 秋霞伦理黄片| 亚洲熟女毛片儿| 80岁老熟妇乱子伦牲交| 国产一区二区三区av在线| 天堂8中文在线网| 国产精品一二三区在线看| 午夜免费男女啪啪视频观看| 国产熟女午夜一区二区三区| 午夜91福利影院| 丁香六月天网| 欧美激情高清一区二区三区 | 亚洲国产看品久久| 国产精品 国内视频| 丰满乱子伦码专区| 婷婷色av中文字幕| 人人妻人人澡人人爽人人夜夜| 一级毛片 在线播放| 好男人视频免费观看在线| 老熟女久久久| 亚洲七黄色美女视频| 国产亚洲欧美精品永久| 建设人人有责人人尽责人人享有的| 人妻一区二区av| 亚洲一卡2卡3卡4卡5卡精品中文| 精品久久久精品久久久| 亚洲视频免费观看视频| 街头女战士在线观看网站| 日日摸夜夜添夜夜爱| 丝瓜视频免费看黄片| 大陆偷拍与自拍| 亚洲精品乱久久久久久| 热re99久久精品国产66热6| 精品少妇黑人巨大在线播放| 亚洲av男天堂| 国产日韩欧美视频二区| 一本一本久久a久久精品综合妖精| 一区在线观看完整版| 日日摸夜夜添夜夜爱| 男女午夜视频在线观看| 色精品久久人妻99蜜桃| 国产99久久九九免费精品| 免费黄网站久久成人精品| 天天添夜夜摸| www.精华液| 久久人人爽人人片av| 尾随美女入室| 日日撸夜夜添| 一级a爱视频在线免费观看| 亚洲欧美成人精品一区二区| 日本av免费视频播放| 中文字幕av电影在线播放| 日日啪夜夜爽| 亚洲精品视频女| 亚洲情色 制服丝袜| 久久久久久人人人人人| 精品免费久久久久久久清纯 | 水蜜桃什么品种好| 九草在线视频观看| 又粗又硬又长又爽又黄的视频| 国产精品熟女久久久久浪| 考比视频在线观看| 99热国产这里只有精品6| 黄色怎么调成土黄色| 成年av动漫网址| 老司机影院毛片| 日本wwww免费看| 高清欧美精品videossex| h视频一区二区三区| 一区二区三区精品91| 校园人妻丝袜中文字幕| 日韩精品免费视频一区二区三区| 只有这里有精品99| 亚洲国产最新在线播放| 精品国产一区二区三区四区第35| 街头女战士在线观看网站| 午夜日韩欧美国产| 久久久精品免费免费高清| 男人操女人黄网站| 不卡av一区二区三区| 两个人免费观看高清视频| 免费高清在线观看视频在线观看| 天天躁日日躁夜夜躁夜夜| av视频免费观看在线观看| 久久性视频一级片| 亚洲四区av| 国产又色又爽无遮挡免| 99热网站在线观看| 丰满饥渴人妻一区二区三| 日韩熟女老妇一区二区性免费视频| 各种免费的搞黄视频| 亚洲欧美日韩另类电影网站| 大陆偷拍与自拍| 久久久久精品国产欧美久久久 | 黄色毛片三级朝国网站| 日韩 欧美 亚洲 中文字幕| 久久99精品国语久久久| 在现免费观看毛片| videosex国产| 青青草视频在线视频观看| 伊人久久国产一区二区| 大香蕉久久成人网| 精品卡一卡二卡四卡免费| 亚洲精品久久成人aⅴ小说| 婷婷色综合大香蕉| av不卡在线播放| 宅男免费午夜| 亚洲欧美成人综合另类久久久| 国产精品熟女久久久久浪| 亚洲男人天堂网一区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美精品综合一区二区三区| 成人亚洲精品一区在线观看| 精品一区二区三区av网在线观看 | 国产熟女午夜一区二区三区| a 毛片基地| 欧美日本中文国产一区发布| 18禁裸乳无遮挡动漫免费视频| 免费黄色在线免费观看| 欧美日韩综合久久久久久| 国产高清国产精品国产三级| 日韩 欧美 亚洲 中文字幕| 69精品国产乱码久久久| 最近2019中文字幕mv第一页| 国产av精品麻豆| 欧美中文综合在线视频| av在线app专区| 亚洲激情五月婷婷啪啪| 热99国产精品久久久久久7| 国产伦理片在线播放av一区| 成人三级做爰电影| 亚洲av综合色区一区| 在线天堂最新版资源| 国产一区二区 视频在线| 久久久国产欧美日韩av| 久久99精品国语久久久| 久久97久久精品| 丝瓜视频免费看黄片| 欧美激情高清一区二区三区 | av视频免费观看在线观看| 国产麻豆69| 青春草国产在线视频| 叶爱在线成人免费视频播放| 久久人人爽av亚洲精品天堂| 日韩电影二区| 国产成人精品久久久久久| 亚洲精品久久午夜乱码| 少妇被粗大的猛进出69影院| 桃花免费在线播放| 中文字幕亚洲精品专区| 欧美亚洲 丝袜 人妻 在线| 高清欧美精品videossex| 亚洲在久久综合| 国产男人的电影天堂91| 亚洲av在线观看美女高潮| 婷婷色麻豆天堂久久| 国产免费一区二区三区四区乱码| 啦啦啦 在线观看视频| 国产精品人妻久久久影院| 在线亚洲精品国产二区图片欧美| 国产精品 国内视频| 亚洲国产最新在线播放| 性高湖久久久久久久久免费观看| 一级片'在线观看视频| 欧美日韩亚洲高清精品| 超色免费av| 国产精品久久久久成人av| 国产在线视频一区二区| 搡老乐熟女国产| 在线观看人妻少妇| 少妇人妻久久综合中文| 国产精品蜜桃在线观看| 超色免费av| 亚洲精品日本国产第一区| 国产一区二区激情短视频 | 最近最新中文字幕免费大全7| 亚洲欧美一区二区三区久久| 国产成人系列免费观看| xxx大片免费视频| 女性被躁到高潮视频| 国产精品麻豆人妻色哟哟久久| 母亲3免费完整高清在线观看| 国产片内射在线| 亚洲精品国产av蜜桃| 在线天堂最新版资源| 熟女少妇亚洲综合色aaa.| 一边摸一边抽搐一进一出视频| 国产精品久久久人人做人人爽| 日本色播在线视频| 观看美女的网站| 蜜桃国产av成人99| 99国产综合亚洲精品| 国产精品免费大片| 国产爽快片一区二区三区| 国产欧美亚洲国产| 午夜日本视频在线| 国产精品女同一区二区软件| 丝袜喷水一区| 日韩一本色道免费dvd| 啦啦啦中文免费视频观看日本| 亚洲精品一区蜜桃| 女性生殖器流出的白浆| 欧美精品一区二区免费开放| 在线观看国产h片| 曰老女人黄片| 日韩中文字幕欧美一区二区 | h视频一区二区三区| 香蕉丝袜av| 可以免费在线观看a视频的电影网站 | 麻豆精品久久久久久蜜桃| 秋霞伦理黄片| 久久99一区二区三区| 18禁裸乳无遮挡动漫免费视频| 五月天丁香电影| 久久精品熟女亚洲av麻豆精品| 成人毛片60女人毛片免费| 在线观看国产h片| 国产无遮挡羞羞视频在线观看| 午夜福利网站1000一区二区三区| 一本大道久久a久久精品| 国产成人精品福利久久| 丁香六月欧美| 性高湖久久久久久久久免费观看| 91精品国产国语对白视频| 久久av网站| 亚洲成国产人片在线观看| 老熟女久久久| 精品久久久久久电影网| 晚上一个人看的免费电影| 高清黄色对白视频在线免费看| 观看美女的网站| 18禁国产床啪视频网站| 国产精品 欧美亚洲| 久久国产精品男人的天堂亚洲| av片东京热男人的天堂| 亚洲av在线观看美女高潮| 操出白浆在线播放| 自线自在国产av| 欧美人与性动交α欧美精品济南到| 亚洲五月色婷婷综合| 免费女性裸体啪啪无遮挡网站| 美女福利国产在线| 国产精品一区二区在线观看99| 一区二区三区四区激情视频| 五月开心婷婷网| 久久影院123| 久久久欧美国产精品| 午夜激情久久久久久久| 精品一品国产午夜福利视频| 国产亚洲欧美精品永久| 欧美黑人精品巨大| 精品久久蜜臀av无| 老汉色∧v一级毛片| 欧美激情高清一区二区三区 | 制服丝袜香蕉在线| 国产av国产精品国产| 久久精品熟女亚洲av麻豆精品| 女人精品久久久久毛片| 在线观看一区二区三区激情| 国产精品亚洲av一区麻豆 | 国产精品 欧美亚洲| 日韩精品免费视频一区二区三区| 2018国产大陆天天弄谢| 国产精品香港三级国产av潘金莲 | 中文乱码字字幕精品一区二区三区| 国产精品二区激情视频| 久久久久久久久免费视频了| 日韩,欧美,国产一区二区三区| 日韩成人av中文字幕在线观看| 一区福利在线观看| 国产亚洲午夜精品一区二区久久| h视频一区二区三区| 国产激情久久老熟女| 成年动漫av网址| 美女扒开内裤让男人捅视频| 国产又爽黄色视频| 多毛熟女@视频| 美女中出高潮动态图| 午夜激情久久久久久久| 欧美变态另类bdsm刘玥| 黑人欧美特级aaaaaa片| 久久亚洲国产成人精品v| 制服人妻中文乱码| 男女床上黄色一级片免费看| xxxhd国产人妻xxx| 成人亚洲精品一区在线观看| 精品卡一卡二卡四卡免费| 黑丝袜美女国产一区| 亚洲国产最新在线播放| 亚洲欧美精品自产自拍| 一边摸一边抽搐一进一出视频| 久久99一区二区三区| 美女国产高潮福利片在线看| 视频在线观看一区二区三区| 国产精品国产三级国产专区5o| av在线观看视频网站免费| 亚洲欧美激情在线| 制服诱惑二区| 自线自在国产av| 日韩一卡2卡3卡4卡2021年| 男女午夜视频在线观看| 一级片'在线观看视频| 在线观看免费视频网站a站| 欧美日韩亚洲国产一区二区在线观看 | 日日啪夜夜爽| 自拍欧美九色日韩亚洲蝌蚪91| 男人爽女人下面视频在线观看| 成人毛片60女人毛片免费| a 毛片基地| 亚洲精品日韩在线中文字幕| 韩国高清视频一区二区三区| 街头女战士在线观看网站| 老司机靠b影院| av免费观看日本| 18禁观看日本| 性少妇av在线| 午夜影院在线不卡| 丰满乱子伦码专区| 不卡视频在线观看欧美| 国产成人a∨麻豆精品| 国产一区二区三区综合在线观看| 热re99久久国产66热| 亚洲av电影在线进入| 亚洲国产最新在线播放| 国产女主播在线喷水免费视频网站| 国产精品一区二区精品视频观看| 亚洲美女黄色视频免费看| 女的被弄到高潮叫床怎么办| 精品一区二区三区四区五区乱码 | 亚洲人成网站在线观看播放| 天美传媒精品一区二区| 国产精品.久久久| 丝袜美足系列| 免费观看av网站的网址| 一级a爱视频在线免费观看| 一区二区三区精品91| 欧美成人精品欧美一级黄| 国产精品秋霞免费鲁丝片| 女人高潮潮喷娇喘18禁视频| 熟女少妇亚洲综合色aaa.| 国产午夜精品一二区理论片| 看免费av毛片| 叶爱在线成人免费视频播放| 黄色视频在线播放观看不卡| 欧美精品一区二区免费开放| 久久精品国产a三级三级三级| 亚洲精品乱久久久久久| 国产午夜精品一二区理论片| 亚洲色图综合在线观看| 欧美日韩综合久久久久久| www日本在线高清视频| 欧美精品人与动牲交sv欧美| 免费看av在线观看网站| 国产一区二区在线观看av| 欧美精品高潮呻吟av久久| 国产野战对白在线观看| 日本vs欧美在线观看视频| 天天影视国产精品| 久久久精品94久久精品| 亚洲国产日韩一区二区| 咕卡用的链子| 交换朋友夫妻互换小说| 亚洲av综合色区一区| 久久久精品国产亚洲av高清涩受| 欧美精品高潮呻吟av久久| 一级毛片电影观看| 欧美激情极品国产一区二区三区| 高清欧美精品videossex| a级毛片在线看网站| 亚洲精品中文字幕在线视频| 精品久久久久久电影网| 久久婷婷青草| xxxhd国产人妻xxx| 欧美日韩一区二区视频在线观看视频在线| 黄色视频在线播放观看不卡| 日本午夜av视频| 亚洲图色成人| 成人三级做爰电影| 在线精品无人区一区二区三| 女人爽到高潮嗷嗷叫在线视频| 午夜福利网站1000一区二区三区| 亚洲av综合色区一区| 免费观看a级毛片全部| 中文字幕制服av| 亚洲欧美成人综合另类久久久| 国产日韩欧美在线精品| 国产一区二区 视频在线| 国产免费现黄频在线看| 2021少妇久久久久久久久久久| 亚洲一区二区三区欧美精品| 欧美亚洲日本最大视频资源| 18禁裸乳无遮挡动漫免费视频| 精品人妻在线不人妻| 人人澡人人妻人| 久久鲁丝午夜福利片| 99久久99久久久精品蜜桃| 亚洲国产精品一区二区三区在线| 热re99久久精品国产66热6| 菩萨蛮人人尽说江南好唐韦庄| 久久婷婷青草| 一区二区三区四区激情视频| 大话2 男鬼变身卡| 国产精品久久久久久精品古装| 日日撸夜夜添| 1024香蕉在线观看| 精品国产国语对白av| 国产成人精品无人区| 另类精品久久| 各种免费的搞黄视频| 熟女av电影| 好男人视频免费观看在线| 美女大奶头黄色视频| 欧美日韩亚洲高清精品| 成人三级做爰电影| 亚洲国产日韩一区二区| 亚洲精品日本国产第一区| 亚洲美女黄色视频免费看| 男女午夜视频在线观看| 成人18禁高潮啪啪吃奶动态图| 久久女婷五月综合色啪小说| 一二三四在线观看免费中文在| av免费观看日本| 国产日韩欧美在线精品| 校园人妻丝袜中文字幕| 久久精品国产亚洲av高清一级| 亚洲三区欧美一区| 少妇人妻久久综合中文| 亚洲精华国产精华液的使用体验| 男男h啪啪无遮挡| 国产成人精品久久二区二区91 | 亚洲精品aⅴ在线观看| 丰满饥渴人妻一区二区三| 国产成人免费观看mmmm| 亚洲,欧美精品.| 亚洲国产欧美网| 亚洲精品视频女| 丝瓜视频免费看黄片| 欧美国产精品一级二级三级| 国产免费现黄频在线看| 少妇人妻精品综合一区二区| 成人国产麻豆网| 成年动漫av网址| 日日摸夜夜添夜夜爱| 交换朋友夫妻互换小说| 熟妇人妻不卡中文字幕| 黑人猛操日本美女一级片| 色94色欧美一区二区| 狂野欧美激情性bbbbbb| 精品久久久久久电影网| 亚洲欧美成人精品一区二区| 日韩免费高清中文字幕av| 91国产中文字幕| 男女午夜视频在线观看| 人人澡人人妻人| 精品免费久久久久久久清纯 | 欧美最新免费一区二区三区| 嫩草影院入口| 久久久久久久久久久久大奶| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧洲精品一区二区精品久久久 | 亚洲av成人精品一二三区| 成年动漫av网址| 男男h啪啪无遮挡| 午夜日本视频在线| 精品亚洲成国产av| 黑人欧美特级aaaaaa片| 亚洲欧美一区二区三区黑人| 婷婷成人精品国产| 亚洲av福利一区| 国产探花极品一区二区| 国产精品嫩草影院av在线观看| 午夜av观看不卡| 午夜影院在线不卡| 亚洲色图综合在线观看| 久久免费观看电影| 涩涩av久久男人的天堂| 人人妻人人添人人爽欧美一区卜| 国产成人精品久久二区二区91 | 91aial.com中文字幕在线观看| 午夜影院在线不卡| 日日爽夜夜爽网站| 亚洲激情五月婷婷啪啪| 在线观看人妻少妇| 久久99热这里只频精品6学生| 久久精品国产亚洲av高清一级| 久久久久久人人人人人| 伦理电影免费视频| 免费观看人在逋| 超碰97精品在线观看| 大片电影免费在线观看免费| 啦啦啦 在线观看视频| 亚洲欧美精品综合一区二区三区| 免费不卡黄色视频| 精品少妇黑人巨大在线播放| a 毛片基地| 99热网站在线观看| 亚洲精品美女久久久久99蜜臀 | 亚洲视频免费观看视频| 免费在线观看视频国产中文字幕亚洲 | 久久精品国产亚洲av涩爱| 下体分泌物呈黄色| 国产成人91sexporn| 国产精品二区激情视频| 欧美日韩av久久| 操美女的视频在线观看| 精品人妻熟女毛片av久久网站| 黄色视频在线播放观看不卡| av片东京热男人的天堂| 久久久久精品人妻al黑| 最近的中文字幕免费完整| 国产精品亚洲av一区麻豆 | 亚洲一区中文字幕在线| 菩萨蛮人人尽说江南好唐韦庄| 深夜精品福利| 午夜激情av网站| 十分钟在线观看高清视频www| 午夜福利在线免费观看网站| 激情视频va一区二区三区| 亚洲国产欧美在线一区| 午夜福利免费观看在线| 亚洲熟女毛片儿| av网站在线播放免费| 久久国产精品大桥未久av| 999精品在线视频| 男女高潮啪啪啪动态图| 51午夜福利影视在线观看| 嫩草影视91久久| 成人漫画全彩无遮挡| 欧美精品人与动牲交sv欧美| 涩涩av久久男人的天堂|