李娜，陳金松，丁宗華，許正文，林樂科

（中國電波傳播研究所電波環(huán)境特性及?；夹g(shù)重點實驗室，山東青島266107）

臨近空間大氣擾動風(fēng)場的探測與分析

李娜，陳金松，丁宗華，許正文，林樂科

（中國電波傳播研究所電波環(huán)境特性及?；夹g(shù)重點實驗室，山東青島266107）

目的深入了解掌握臨近空間大氣風(fēng)場的擾動狀態(tài)，給臨近空間飛行器的精確入軌、有效運行提供可靠的環(huán)境信息。方法用昆明站(25.6°N,103.8°E)臨近空間觀測平臺中的MF雷達(dá)開展臨近空間大氣風(fēng)場擾動觀測。基于諧波擬合方法得到昆明地區(qū)2008年8月至2010年8月期間的周日和半日潮汐分量及背景風(fēng)場，分析昆明地區(qū)臨近空間80～100 km之間潮汐波動的變化特征。結(jié)果周日潮汐振幅具有明顯的年變化和半年變化，最大值均出現(xiàn)在10月份和3月份的86～92 km高度范圍內(nèi)，同時第一年度的潮汐幅度明顯強于第二年度，類似準(zhǔn)兩年變化。緯向半日潮汐的變化特征不是非常明顯，而經(jīng)向半日潮汐在88～94 km之間具有年變化，在其他高度處則出現(xiàn)多個峰值。結(jié)論潮汐相位的變化結(jié)果表明潮汐的下行傳播，潮汐振幅的高度變化符合先增大后減小的特征，呈現(xiàn)胞狀結(jié)構(gòu)。

臨近空間；潮汐；MF雷達(dá)

臨近空間是各種高空氣球、高空飛行器、亞軌道飛行器和低軌航天器的駐留區(qū)，臨近空間大氣參數(shù)的擾動和變化過程對臨近空間大飛行器的安全運行具有重要影響。在各種臨近空間環(huán)境參數(shù)中，風(fēng)場是非常重要的一個參數(shù)，它直接影響臨近空間飛行器的精確入軌和可靠運行。同時風(fēng)場中存在豐富的波動現(xiàn)象，比如重力波、潮汐波、行星波等，具有重要的科學(xué)意義。最新研究顯示，由太陽和地磁的周期性變化以及地磁暴引起的高層大氣的周期性擾動（包括風(fēng)場周期性震蕩）會對衛(wèi)星軌道高度參數(shù)產(chǎn)生重要影響[1—5]。

臨近空間風(fēng)場具有顯著的地域性和多時間尺度（包括日變化、半年和年變化等）的變化特征。昆明站位于低緯地區(qū)，靠近青藏高原，其臨近空間大氣的變化特征非常獨特。關(guān)于昆明地區(qū)臨近空間的參數(shù)特性研究可見于文獻(xiàn)[6—8]。

中頻（Medium Frequency,簡稱MF）雷達(dá)可以實現(xiàn)白天60～100 km和晚上80～100 km高度大氣風(fēng)場及電子密度的持續(xù)測量。20世紀(jì)80年代以來，探測技術(shù)和分析方法的發(fā)展，使MF雷達(dá)具有了設(shè)備簡單、運行方便和無人值守的優(yōu)點，因此MF雷達(dá)站點迅速增加[9]。目前，在歐洲、北美、南極、澳大利亞、南非、印度、日本和印度尼西亞等國家和地區(qū)都建有MF雷達(dá)站，這些雷達(dá)為臨近空間風(fēng)場及電子密度研究提供了的重要觀測數(shù)據(jù)。我國于2000年在武漢建立了第一臺MF雷達(dá)[10]，胡雄、張冬婭、張秀川和易帆等人利用該雷達(dá)對武漢中層低熱層的大氣角譜[11]、平均風(fēng)[12]、潮汐變化及相互作用進行了研究[13—14]。由于國內(nèi)在這方面的起步較晚，制約著我國大氣與空間環(huán)境觀測和研究工作的開展。

文中利用昆明站MF雷達(dá)2年的風(fēng)場觀測數(shù)據(jù)，分析了昆明地區(qū)臨近空間大氣風(fēng)場中波動擾動的年和半年尺度的變化特征。

1 MF雷達(dá)工作原理及數(shù)據(jù)處理

1.1 工作原理

60～100 km高度范圍內(nèi)存在幾米到幾百米尺度的電離層不均勻團塊，其界面上存在折射指數(shù)梯度，導(dǎo)致入射電磁波發(fā)生部分反射。該高度范圍內(nèi)還存在其他回波反射和散射機制，如湍流散射和Thompson散射等。MF雷達(dá)主要利用部分反射過程，由于其反射信號強度除與電子密度的梯度成正比外，還與電磁波頻率的平方成反比，因此雷達(dá)工作頻率越低，回波信號越強。當(dāng)頻率太低時，電磁波將在較低高度上發(fā)生全反射，不能夠得到較高高度上的回波，故觀測上通常選擇1.8～2.5 MHz之間的頻率[10]，昆明站MF雷達(dá)選擇的頻率為2.138 MHz。關(guān)于雷達(dá)的主要技術(shù)指標(biāo)參數(shù)：位置為25.6°N,103.8°E，峰值功率64 kW，工作頻率為2.138 MHz，采樣間隔為2 km，觀測模式為FCA，IDI，DAE，DPE，為開始運行時間2007年12月。

該雷達(dá)采用分布式天線系統(tǒng)（Space Antenna，簡稱SA），由4副相互獨立饋電的交叉偶極子天線按等邊三角形布局形成收發(fā)共用天線陣，其中三副天線與另外一副呈120°輻射狀排列，每副交叉偶極子天線的偶極子交叉成90°，布局如圖1所示。天線單元為傳統(tǒng)的半波偶極子，也就是天線陣一共包含8個半波偶極子，可以方便地形成右圓極化和左圓極化模式。雷達(dá)的結(jié)構(gòu)如圖2所示，可以看出，雷達(dá)工作時主要有4部分組成：信號發(fā)射、信號接收、數(shù)據(jù)采集和數(shù)據(jù)分析顯示。

圖1 昆明MF雷達(dá)天線陣布局

圖2 昆明MF雷達(dá)硬件結(jié)構(gòu)

雷達(dá)系統(tǒng)通過兩種不同的方式分別測量大氣風(fēng)場與電子密度。探測風(fēng)場時，雷達(dá)發(fā)射一種偏振波尋常波（O波）或非尋常波（X波），通過分析天線接收到的回波數(shù)據(jù)反演風(fēng)場。昆明MF雷達(dá)反演風(fēng)場時可以采用全相關(guān)分析方法[15—17]或多普勒干涉成像法（Imaging Doppler Interferometry，簡稱IDI），其理論基礎(chǔ)是衍射樣本的空間時間相關(guān)函數(shù)為一組橢球，在不同位置對衍射樣本取樣，計算所得樣本信號的自相關(guān)和互相關(guān)函數(shù)，可以得到橢球參數(shù)，最后通過互相關(guān)函數(shù)的時間置換可以得到視線速度（Apparent Velocity）。由于衍射樣本的運動及不對稱性，需要對自相關(guān)函數(shù)進一步計算才能得到真實速度（True Velocity）[16]。探測電子密度時，雷達(dá)交替發(fā)射O波和X波，由于MF雷達(dá)天線能夠分辨下行波的O波與X波，通過測量兩個特征波的振幅與相位可以得到電子密度等信息[18]。反演電子密度時可以采用微分吸收試驗法或微分相位試驗法（Differential Phase Experiment，簡稱DPE）。昆明MF雷達(dá)目前選擇FCA與DAE交替進行的模式來探測大氣風(fēng)速及電子密度。有關(guān)FCA觀測模式的參數(shù)見表1。

表1 昆明MF雷達(dá)FCA觀測主要參數(shù)

1.2 數(shù)據(jù)預(yù)處理

選擇2008年8月至2010年8月期間80～100 km高度范圍內(nèi)的風(fēng)場觀測數(shù)據(jù)作為研究對象，其中2009年8月缺失。在數(shù)據(jù)的分析過程中，首先通過線性插值得到每個高度上的小時平均值。其原則是如果每小時內(nèi)的缺失數(shù)據(jù)個數(shù)少于4，那么該小時內(nèi)的數(shù)據(jù)是有效的；如果大于4個，將設(shè)該小時的平均值為NAN，且不參加后續(xù)運算。然后將每個月的數(shù)據(jù)作為一個數(shù)據(jù)序列，如果該數(shù)據(jù)序列中NAN的個數(shù)小于一半，則利用最小二乘法對插值后的數(shù)據(jù)進行諧波擬合，得到背景風(fēng)場、周日潮汐分量和半日潮汐分量。如果數(shù)據(jù)序列中NAN的個數(shù)大于等于一半，則視為潮汐分量分別為NAN。采用的諧波擬合公式為：

式中：y(t)為觀測到經(jīng)向風(fēng)或緯向風(fēng)；y0為背景風(fēng)場；A1與φ1為周日潮汐的振幅和相位；A2與φ2為半日潮汐的振幅和相位。

2 雷達(dá)初步探測結(jié)果

2.1 緯向潮汐分量

圖3給出了昆明地區(qū)2008年8月至2010年8月期間緯向周日潮汐和半日潮汐的振幅和相位隨時間高度的等高變化。其中橫坐標(biāo)表示每年的2月（February）、5月(May)、8月(August)、11月(November)，LT表示地方時，即潮汐振幅最大時對應(yīng)的當(dāng)?shù)貢r。

由圖3可知，周日潮汐分量振幅在5～35 m/s之間變化，半日潮汐分量振幅在1～14 m/s之間變化。相比來說，周日潮汐的年變化和半年變化較為顯著，每年的10月份和3月份潮汐振幅達(dá)到最大值，且2008—2009年的潮汐振幅明顯大于2009—2010年，而半日潮汐振幅則無明顯的年或半年變化特征。半日

圖3 2008—2010期間昆明緯向周日、半日潮汐振幅和相位等高線

振幅也具有2008—2009幅度大于2009—2010年幅度的特征，這可能與低緯地區(qū)的準(zhǔn)兩年變化有關(guān)[19]。在潮汐的垂直結(jié)構(gòu)方面，周日振幅的最大值出現(xiàn)在86～90 km之間，而半日振幅則是在90 km以上。周日潮汐振幅最大值出現(xiàn)的高度可能與風(fēng)剪切有關(guān)，因為對昆明地區(qū)緯向平均風(fēng)變化特征的分析研究中發(fā)現(xiàn)，90 km處緯向風(fēng)由東向風(fēng)急變?yōu)槲飨蝻L(fēng)，通常認(rèn)為風(fēng)剪切現(xiàn)象是重力波造成的[20]。分析結(jié)果與日本的Yamagawa(31.2°N,130.6°E)站的結(jié)果[21]具有明顯的不同。兩個潮汐振幅的共同特點是振幅的高度變化均呈胞狀結(jié)構(gòu)，即先隨高度的增加而增加，之后減小。從數(shù)值上看，周日潮汐大于半日潮汐。昆明站地處我國低緯地區(qū)，海拔2000 m，是周日潮汐較為活躍的區(qū)域，這可能是造成周日潮汐成為該地區(qū)最主要波動形態(tài)的原因。

通常潮汐振幅隨高度的增加而增加是因為大氣密度隨高度降低，所以在中層頂區(qū)域潮汐振幅會變得很大導(dǎo)致其他非線性相互作用的發(fā)生。比如波破碎、潮汐與重力波之間的非線性作用、潮汐與行星波之間的非線性作用等[22—28]，之后振幅則隨高度增加而減小。這是風(fēng)剪切的形成理論。

潮汐相位即潮汐振幅達(dá)到最大值時的時刻。由圖3所示，相位的變化特征進一步驗證了潮汐的年變化和半年變化特征。

2.2 經(jīng)向潮汐分量

昆明地區(qū)2008年8月至2010年8月期間經(jīng)向周日潮汐和半日潮汐的振幅和相位隨時間高度的等高變化如圖4所示。首先經(jīng)向周日潮汐振幅和半日潮汐振幅的變化范圍分別是1～40 m/s和1～20 m/s，稍大于緯向分量振幅幅度。經(jīng)向周日潮汐與緯向表現(xiàn)出一致的變化特點，即同時具有年變化和半年變化，最大值出現(xiàn)在10月份和3月份的86～90 km之間。半日潮汐與緯向不同，經(jīng)向半日潮汐具有非常明顯的年變化，即最大值出現(xiàn)在每年12—4月份的88～94 km之間。由圖4可知，2008—2009年經(jīng)向潮汐分量幅度大于2009—2010年經(jīng)向潮汐分量幅度，這與潮汐波動的準(zhǔn)兩年變化有關(guān)，更多的特征驗證需要更多的數(shù)據(jù)支持。從高度區(qū)域上來看，經(jīng)向半日潮汐在94 km以上活動更頻繁，比如經(jīng)向半日潮汐振幅分別在9，11，2，5月出現(xiàn)峰值，同時在84 km以下也出現(xiàn)類似特征。這可能與中層頂區(qū)域波動相互之間的非線性作用有關(guān)。經(jīng)向潮汐振幅的高度分布同樣經(jīng)歷了隨高度先增大后減小的變化特征。

圖4 2008—2010年期間昆明經(jīng)向周日、半日潮汐振幅和相位等高線

3 結(jié)論

大氣潮汐是全球大氣中常見的一種波動現(xiàn)象，周期為一天或其整數(shù)分之一，一般可向上傳播至熱層底部（電離層E區(qū)）甚至更高，在中間層-低熱層大氣水平和垂直風(fēng)場中都具有很強的幅度，是影響中性大氣-電離層耦合過程的重要因素，也是影響臨近空間飛行器和航天器精確飛行的重要因素。潮汐波的變化在時間和空間上均具有多尺度的特征。從時間尺度來看，潮汐既具有長周期的變化，諸如季節(jié)變化、年變化、半年變化、準(zhǔn)兩年變化甚至更長周期的規(guī)則與異常變化，又存在短周期變化，諸如日變化。從空間尺度來看，潮汐既具有局域地區(qū)變化，又具有隨緯度或經(jīng)度的全球尺度變化。文中研究工作主要是基于昆明站MF雷達(dá)2008—2010年期間風(fēng)場觀測數(shù)據(jù)分析昆明地區(qū)臨近空間大氣風(fēng)場中潮汐波動的變化特征。

1）緯向（東西向）周日潮汐和經(jīng)向（南北向）周日潮汐均具有年變化和半年變化，最大值均出現(xiàn)在10月份和3月份的86～92 km高度范圍內(nèi)。2008—2009年潮汐幅度大于2009—2010年潮汐幅度，這與潮汐的準(zhǔn)兩年變化有關(guān)。以往的研究證明，準(zhǔn)兩年波動在低緯地區(qū)起主導(dǎo)作用[20—22]。潮汐的垂直結(jié)構(gòu)都經(jīng)歷了先隨高度增加而增加，達(dá)到最大值后隨高度減小。

2）緯向半日潮汐沒有表現(xiàn)出明顯的特征，但經(jīng)向半日潮汐的一個重要特征是88～94 km之間為年變化，即12月—4月期間的潮汐幅度要大于其他月份。94 km以上高度半日潮汐振幅出現(xiàn)多個峰值，比如在9，11，2，5月份經(jīng)向半日潮汐振幅明顯大于其他月份，而在84 km以下也出現(xiàn)類似特征。

文中僅給出了昆明地區(qū)兩年的風(fēng)場潮汐變化，嚴(yán)格上得出昆明地區(qū)潮汐具有年變化特征和半年變化特征有待于更多數(shù)據(jù)的支持。潮汐的準(zhǔn)兩年變化是最近研究的熱點，文中指出昆明地區(qū)大氣潮汐具有類準(zhǔn)兩年變化，下一步是對更多的數(shù)據(jù)進行處理分析。

文中主要目的是描繪昆明地區(qū)臨近空間環(huán)境擾動風(fēng)場的變化特征，為融合多手段的探測數(shù)據(jù)、結(jié)合動力學(xué)理論等建立臨近空間環(huán)境擾動風(fēng)場模型，進一步為臨近空間武器裝備的高效運行等提供理論保障。

[1]XU J,SMITH A K,LIU H,et al.Seasonal and Quasibiennial Variations in the Migrating Diurnal Tide Observed by Thfermosphere,Ionosphere,Mesosphere, Energetics and Dynamics(TIMED)[J].J Geophys Res, 2010,114(D13):267—275.

[2]LIU X,XU J,ZHANG S,et al.Thermospheric Planetary Wave-type Oscillations Observed by FPIs over Xinglong and Millstone Hill[J].J Geophys Res Space Physics, 2015,119(8):6891—6901.

[3]YI Fan.Short-term Variability and Temporary Structures of Tides and Mean Wind in the Polar Summer Mesospher[J].J Atmos Solar Terr Phys,2001,63: 749—757.

[4]ZHAO L,CHEN J,DING Z,et al.First Observations of Tidal Oscilations by an MF Radar over Kunming (25.6°N,103.8°E)[J].Journal of Atmospheric and Solar-Terrestrial Physics,2012,s78–79:44—52.

[5]XU J,WANG W,LEI J,et al.The Effect of Periodic Variations of Thermospheric Density on CHAMP and GRACE Orbits[J].J Geophys Res,2011,116(A2):1—11.

[6]李娜,陳金松,趙蕾.基于MF雷達(dá)觀測的D區(qū)日食效應(yīng)的研究[J].地球物理學(xué)報,2011,54(4):942—949.

[7]LI Na,CHEN Jin-song,DING Zong-hua,et al.The Mean Winds Observed by an MF Radar over Kunming (25.6°N,103.8°E)[J].Journal of Atmosphere Solarterrestrical Physics,2015(122):58—65.

[8]YI Wen,CHEN Jin-song,Ma Chun-bo,et al.Observation of Upper Atmospheric Temperature by Kunming All-sky Meteor Radar[J].Chinese Journal of Geophysics,2014, 57(5):750—760.

[9]HOLDSWORTH D A,REID I M.A Simple Model of Atmospheric Radar Backscatter Description and Application to the Full Correlation Analysis of Spaced Antenna Data[J].Radio Sci,1995,30(4):1263—1280

[10]李鳳琴,胡雄,張訓(xùn)械,等.武漢中層大氣中頻雷達(dá)及其初步探測結(jié)果[J].空間科學(xué)學(xué)報,2002,22(1):65—71.

[11]胡雄,曾楨,張冬婭,等.武漢中層?低熱層大氣角譜中頻雷達(dá)觀測[J].空間科學(xué)學(xué)報,2003,23(4): 256—261.

[12]張冬婭,胡雄,張訓(xùn)械,等.北緯30°N中間層和低熱層大氣平均風(fēng)中頻雷達(dá)觀測[J].空間科學(xué)學(xué)報,2005, 25(4):267—272.

[13]張紹東,易帆,胡雄.武漢上空(30°N,114°E)潮汐及其相互作用的MF雷達(dá)觀測[J].空間科學(xué)學(xué)報,2003, 23(6):430—435.

[14]張秀川,易帆.武漢上空中層頂區(qū)域潮汐的MF雷達(dá)觀測[J].空間科學(xué)學(xué)報,2007,27(1):29—34.

[15]BRIGGS B H.The Analysis of Spaced Sensor Records by Correlation Techniques[J].MAP Handbook,1984,13: 166—186.

[16]ANDREW D M.VHF Boundary Layer Radar and Rass [D].Adelaide:The University of Adelaide,2001.

[17]BRIGGS B H,PHILLIPS G J,SHINN D H.The Analysis of Observations on Spaced Receivers of the Fading of Radio Signals[J].Proc Phys Soc,1950,63B:106—121.

[18]GARDNER F F,PAWSEY J L.Study of the Ionospheric D-region Using Partial Reflection[J].J Atmos Terr Phys, 1953,53(3):321—344.

[19]VERYARD R G,EBDON R A.Fluctuations in Tropical Stratospheric Winds[J].Meteorol,1961,125:90.

[20]SHARMA A K,ROKADE M V,KONDALA R,et al. Comparative Study of MLT Mean Winds Using MF Radars Located at 16.8N and 8.7N[J].J Earth Syst Sci, 2010,119(4):461—470.

[21]NAMBOOTHIRI S P,KISHORE P,IGARASHI K,et al.MF Radar Observations of Mean Winds over Yamagama (31.2°N,130.6°E)and Wakkanai(45.4°N,141.7°E) [J].JAtmos Sol-Terr Phys,2000,62:1177—1187.

[22]FRITTS D C,ISLER J R.Mean Motions and Tidal and Two-day Wave Structure and Variability in the Mesosphere and Lower Thermosphere over Hawaii[J].J Atmos Sci,1994,51:2145—2164.

[23]LIEBERMAN R S,HAYS P B.An Estimate of the Momentum Deposition in the Lower Thermosphere by the Observed Diurnal Tide[J].J Atmos Sci,1994,51: 3094—3105.

[24]NAKAMURA T,FRITTS D C,ISLER J R,et al.Shortperiod Fluctuations of the Diurnal tide Observed with Low-latitude MF and Meteor Radars during CADRE: Evidence for Gravity Wave/Tidal Interactions[J].J Geophys Res,1997,102(26):26225—26238.

[25]TEITELBAUM H,VIAL F.Momentum Transfer to the Thermosphere by Atmospheric Tides[J].J Geophys Res, 1981,86:9693—9697.

[26]WU Q,ORTLAND D A,KILLEEN T L,et al.Global Distribution and Inter-annual Variations of Mesospheric and Lower Thermospheric Neutral Wind Diurnal Tide:1. Migrating Tide[J].J Geophys Res,2008,113(6):211—215.

[27]WU Q,ORTLAND D A,KILLEEN T L,et al.Global Distribution and Inter-annual Variations of Mesospheric and Lower Thermospheric Neutral Wind Diurnal Tide:2. Nonmigrating Tide[J].J Geophys Res,2008,113(5): 237—245.

[28]Xu J,Smith A K,Liu H L,et al.Season and Quasibiennial Variations in the Migrating Diurnal Tide Observed by Thermosphere,Ionosphere,Mesosphere, Energetics and Dynamics(TIMED)[J].J Geophys Res, 2009,114(D13):267—275.

Wind-oscillation Measurement and Study in Near Space

L I Na,CHEN Jin-song,DING Zong-hua,X U Zheng-wen,LIN Le-ke
(China Research Institute of Radio-wave Propagation,National Key Laboratory of Electromagnetic Environment, Qingdao,266107,China)

ObjectiveTo understand variations of oscillation wind and provide exact environmental status for precision entering and running of aerocraft in near space.MethodsThe experiments on wind field turbulence were performed at Kunming station(25.6°N,103.8°E)with Medium frequency radar(MF radar)in the near space observation platform. Based on data on diurnal and semidiurnal tide component and background wind field from Kunming from August 2008 to August 2010 with harmonic fitting method,characteristics of tide fluctuation in 80～100 km of near space in Kunming were analyzed.ResultsThe variations periods of diurnal tide behaved annually and semiannually,although the amplitude of tidal in the first year was larger than that in the second year,which may be thought as quasi-biennial oscillation.The maximum values always occurred in October and March at 86～92 km.But the zonal semidiurnal tidal component has no this feature,at least not obviously.The case became complicated in meridional semidiurnal tidal that annual variation only dominated at 88～94 km region;and several peak values occurred at other height.ConclusionChanges of tide phase show the downward transmission of tide.The height change of tide meets the feature of increasing with height and then decreasing.It shows a cystiform structure.

near space,tidal,MF radar

10.7643/issn.1672-9242.2017.07.007

TN011

A

1672-9242(2017)07-0035-06

2017-03-13；

2017-04-06

國家自然科學(xué)青年基金(41604129)

李娜（1981—），女，河南新鄉(xiāng)人，工程師，主要研究方向為臨近空間大氣環(huán)境探測。