田玉芳, 呂達(dá)仁
1 中國(guó)科學(xué)院大氣物理研究所中層大氣與全球環(huán)境探測(cè)重點(diǎn)實(shí)驗(yàn)室(LAGEO), 北京 100029 2 中國(guó)科學(xué)院大學(xué), 北京 100049
?
北京MST雷達(dá)探測(cè)中間層-低熱層觀測(cè)結(jié)果初步分析
田玉芳1,2, 呂達(dá)仁1*
1 中國(guó)科學(xué)院大氣物理研究所中層大氣與全球環(huán)境探測(cè)重點(diǎn)實(shí)驗(yàn)室(LAGEO), 北京100029 2 中國(guó)科學(xué)院大學(xué), 北京100049
摘要北京MST雷達(dá)是子午工程建設(shè)的國(guó)內(nèi)僅有的兩部MST雷達(dá)之一,為研究其在中間層-低熱層MLT區(qū)域的探測(cè)能力以及數(shù)據(jù)可靠性,本文應(yīng)用北京MST雷達(dá)2012、2013兩年高模式數(shù)據(jù),從數(shù)據(jù)獲取率、與廊坊流星雷達(dá)測(cè)風(fēng)對(duì)比以及風(fēng)場(chǎng)時(shí)空分布特征三個(gè)方面進(jìn)行初步分析.結(jié)果是:(1) 數(shù)據(jù)獲取率日變化特征為:白天65~100 km均可獲取數(shù)據(jù),數(shù)據(jù)獲取率的高值區(qū)主要集中在70~80 km,最大值可達(dá)80%;夜間主要集中在80~100 km,數(shù)據(jù)獲取率在30%及以下.表明該MST雷達(dá)白天可以探測(cè)到電離層D層和E層低層,夜間D層消失,只探測(cè)到E層低層.季節(jié)變化特征為:夏季白天可獲取數(shù)據(jù)的時(shí)間和高度區(qū)間都比較大,春季次之,冬季最小.夏季白天以及日落后1 h內(nèi)可探測(cè)到120 km.(2) 對(duì)北京MST雷達(dá)與廊坊流星雷達(dá)2012年5月份、80~100 km高度區(qū)間測(cè)量的水平風(fēng)進(jìn)行對(duì)比分析,二者測(cè)風(fēng)結(jié)果在時(shí)空分布上有很好的一致性,表明MST雷達(dá)探測(cè)數(shù)據(jù)是可靠的.(3) 2012年和2013年相應(yīng)月份平均的緯向風(fēng)、經(jīng)向風(fēng)時(shí)空分布特征有較高的一致性,并與HWM07模式結(jié)果也基本一致.上述初步分析結(jié)果表明,北京MST雷達(dá)對(duì)中間層-低熱層60~120 km高度區(qū)域已具備較強(qiáng)的探測(cè)能力,所得結(jié)果將可用于MLT過(guò)程揭示與驅(qū)動(dòng)因子研究,并可與該高度上其他探測(cè)手段作綜合研究.
關(guān)鍵詞北京MST雷達(dá); 中間層-低熱層; 探測(cè)能力; 數(shù)據(jù)可靠性; 水平風(fēng)場(chǎng)時(shí)空分布特征
Two years of data (January 2012 to December 2013) in high mode were used. From the aspect of data acquisition rate, comparison of wind measurements with the nearby Langfang meteor radar and spatial and temporal characteristics of wind, preliminary analyses were made.
Analysis results indicate that: (1) Data acquisition rate shows diurnal variation characteristics: during daytime the MST radar can obtain data in 65~100 km. High-ratio areas of data acquisition rate were mainly concentrated in the 70~80 km with a maximum of up to 80%.At nighttime available data altitude range focused on 80~100 km with data acquisition rate at 30% and below. This indicates that the MST radar can detect ionospheric D layer and lower E layer during daytime, at night time D layer disappears, only lower E layer detected. Seasonal variation characteristics of data acquisition rate are that the time and height range of available data were large in summer daytime, followed by spring, and winter was minimum. During summer daytime and until one hour after sunset, the MST radar can detect up to 120 km. (2) Horizontal winds in 80~100 km measured by Beijing MST radar and Langfang meteor radar in May 2012 were compared. The results showed good agreement in spatial and temporal distribution, indicating that the data of MST radar observations is reliable. (3) The spatial and temporal distribution of averaged zonal wind and meridional wind of corresponding month in year 2012 and 2013 had high consistency, and were also fairly consistent with HWM07 model results.
The preliminary analysis results above show that the Beijing MST radar has strong detection ability in the MLT region (60~120 km), and the observed data can be furtherly used to reveal the processes and driving factors of the MLT region.
1引言
中間層-低熱層(MLT)區(qū)域是中性大氣與電離層大氣之間的過(guò)渡區(qū)域,也是中層大氣的一個(gè)重要區(qū)域,從小尺度的湍流到大尺度的波動(dòng)以及它們的耗散在中層大氣的能量和動(dòng)量變化方面都發(fā)揮著重要的作用(Kumar et al.,2014).然而由于探測(cè)能力的限制,對(duì)該區(qū)域的理解和認(rèn)識(shí)仍是遠(yuǎn)遠(yuǎn)不夠的.流星雷達(dá)、MF雷達(dá)、探空火箭以及衛(wèi)星觀測(cè)等均可用于該區(qū)域的觀測(cè),但是時(shí)空分辨率以及探測(cè)高度范圍方面不盡如人意(Kumar et al.,2008).自從1974年世界上第一部MST雷達(dá)問(wèn)世(Woodman and Guillen,1974),為MLT區(qū)域大氣探測(cè)提供了強(qiáng)有力的手段.
MST(Mesosphere Stratosphere Troposphere)雷達(dá)可以探測(cè)對(duì)流層、平流層以及中間層大氣.它利用大氣折射率不規(guī)則體對(duì)雷達(dá)電磁波的散射回波,來(lái)研究大氣的動(dòng)力學(xué)特征——風(fēng)場(chǎng)、波動(dòng)、湍流、大氣穩(wěn)定性以及與大氣折射指數(shù)相關(guān)的濕度、溫度、電子密度及其梯度等信息.
大氣折射率的不均勻引起對(duì)電磁波的散射或反射,在不同高度的大氣層有不同的機(jī)制.對(duì)流層及以下湍流散射起主要作用,有效散射體空間尺度為入射波長(zhǎng)的一半;對(duì)流層頂以及平流層用米波或者分米波雷達(dá)探測(cè)時(shí)存在鏡面反射;而在中間層-低熱層大氣中離子、電子熱運(yùn)動(dòng)形成的電子密度隨機(jī)起伏而產(chǎn)生的熱散射起主要作用(Gage and Balsley,1980).
大氣湍流散射是晴空大氣探測(cè)雷達(dá)回波信號(hào)的主要成因.根據(jù)大氣湍流理論,湍流可以看成多尺度渦旋的疊加,尺度可以小到毫米量級(jí),大到千米量級(jí).大氣中主要的湍渦尺度隨高度是變化的,隨著高度增加,湍渦尺度逐漸增大,小尺度渦旋逐漸減小(Fukao, 2007).因此,MST雷達(dá)工作在較長(zhǎng)的波長(zhǎng),通常為幾米,頻率在50 MHz左右,來(lái)保證有足夠的靈敏度探測(cè)到高層大氣.
晴空大氣探測(cè)雷達(dá)的回波信號(hào)非常微弱,而且隨著高度增加迅速減小.即使是功能強(qiáng)大的MST雷達(dá),在上平流層和下中間層區(qū)域是很難探測(cè)到回波的,因?yàn)樵谶@一高度區(qū)間內(nèi)來(lái)自湍流散射的折射率不規(guī)則體對(duì)雷達(dá)電磁波的散射隨著大氣密度變稀薄而大大減小,而且這一高度也沒(méi)有自由電子產(chǎn)生熱散射使得回波加強(qiáng)(Van Zandt,2000; Woodman and Guillen,1974).因此MST雷達(dá)的有效探測(cè)高度范圍一般在2~25 km,以及60~95 km左右.
表1 國(guó)際已有部分MST雷達(dá)基本信息
MST雷達(dá)技術(shù)從1974年開(kāi)始不斷發(fā)展,已有多個(gè)國(guó)家已經(jīng)布設(shè)并投入使用.表1中列出了國(guó)際上一些MST雷達(dá)的基本信息.
1980s中期,在中國(guó)科學(xué)院支持下大氣物理研究所開(kāi)始獨(dú)立研制大型相控陣天線VHF Doppler雷達(dá).該雷達(dá)建于河北香河觀測(cè)站,工作頻率為78 MHz.自1993年開(kāi)始陸續(xù)開(kāi)展觀測(cè)試驗(yàn),獲取了水平風(fēng)、垂直速度、湍流強(qiáng)度等多種類(lèi)型的觀測(cè)結(jié)果(李衛(wèi)和呂達(dá)仁,1997;Li and Lü,1998).這是國(guó)內(nèi)第一臺(tái)獨(dú)立研制并具有多功能的中層大氣研究雷達(dá)(呂達(dá)仁,陳洪濱,2003).
近年來(lái)我國(guó)MST雷達(dá)得到了較大的發(fā)展.子午工程建設(shè)2部MST雷達(dá),即北京和武漢MST雷達(dá).武漢MST雷達(dá)(114°8′8″E, 29°31′58″N)于2011年3月完成試運(yùn)行(Zhao et al.,2013);北京MST雷達(dá),坐落在河北香河(116°59′24.00″E, 39°45′14.40″N),由中國(guó)科學(xué)院大氣物理研究所負(fù)責(zé)主持建設(shè),2011年5月完成系統(tǒng)連調(diào),2011年7月底完成單位驗(yàn)收.
北京MST雷達(dá)可以獲取到其上空中間層-低熱層(MLT)區(qū)域?qū)氋F的三維風(fēng)場(chǎng)、湍流信息等實(shí)測(cè)資料.在資料應(yīng)用之前,有必要對(duì)MST雷達(dá)在該區(qū)域的探測(cè)能力以及數(shù)據(jù)可靠性進(jìn)行分析.比如將MST雷達(dá)與其他觀測(cè)手段做比較,部分已有的研究包括:Ratnam等(2001)將Indian Gadanki MST雷達(dá)與MF雷達(dá)探測(cè)65~80 km的水平風(fēng)做比較,發(fā)現(xiàn)二者的緯向風(fēng)分量有好的一致性,經(jīng)向風(fēng)分量一致性差一些;而與CIRA-86模式比較結(jié)果顯示緯向風(fēng)分量除秋分外,其他季節(jié)模式值均大于MST雷達(dá)探測(cè)值.Kumar等(2008)應(yīng)用1995—2006較長(zhǎng)時(shí)間段內(nèi)Indian Gadanki MST雷達(dá)資料與該時(shí)間段內(nèi)的MF雷達(dá)、火箭探空、模式(HWM93)以及衛(wèi)星觀測(cè)(HRDI/UARS)的65~85 km高度范圍內(nèi)的水平風(fēng)做比較,發(fā)現(xiàn)幾種觀測(cè)手段在趨勢(shì)上有好的一致性,但在數(shù)值上有較大的差異.其中火箭探空與MST雷達(dá)探測(cè)結(jié)果有好的一致性,HRDI觀測(cè)值大于MST雷達(dá)且經(jīng)向風(fēng)比緯向風(fēng)差值更大些,MF雷達(dá)觀測(cè)值比MST雷達(dá)小,HWM93模式的風(fēng)速也比MST雷達(dá)測(cè)值小.Rao等(2014) 比較了Sri Venkateswara University (SVU) 流星雷達(dá)、MF雷達(dá)以及Gadanki MST雷達(dá)在70~98 km探測(cè)水平風(fēng)情況,發(fā)現(xiàn)MST雷達(dá)與流星雷達(dá)探測(cè)經(jīng)向風(fēng)和緯向風(fēng)的一致性較差,而MST雷達(dá)與MF雷達(dá)探測(cè)的緯向風(fēng)分量一致性較好.
國(guó)外應(yīng)用MST雷達(dá)探測(cè)資料對(duì)中間層-低熱層(MLT)區(qū)域已有較多的研究.Selvaraj等(2014)應(yīng)用Gadanki雷達(dá)對(duì)中間層回波的散射截面,以及湍流參量包括湍流能量耗散率、湍流內(nèi)尺度、外尺度等進(jìn)行計(jì)算分析,并與之前的相關(guān)研究結(jié)果做對(duì)比.Lehmacher等(2009)利用Jicamarca雷達(dá)以及Belova等(2012)應(yīng)用India Gadanki MST雷達(dá)對(duì)低緯地區(qū)中間層雷達(dá)截面以及回波特征進(jìn)行研究.Fukao等(1994)應(yīng)用MU radar,Guo等(2007)應(yīng)用Jicamarca Radar,Rao等(2001)和Chandra等(2008,2012)應(yīng)用Gadankir adar對(duì)中間層區(qū)域的湍流特征進(jìn)行了分析.Kubo等(1997),Kamala等(2003),Kumar等(2007)研究了中間層回波特征.Ratnam等(2008)應(yīng)用MST雷達(dá)與衛(wèi)星以及火箭探空資料對(duì)中間層半年振蕩(MSAO)和準(zhǔn)兩年振蕩(MQBO)進(jìn)行了研究.此外,還有中間層垂直速度(Eswaraiah et al.,2011, 2012),風(fēng)場(chǎng)特征(Kumar et al.,2008)以及KHI不穩(wěn)定(Lehmacher et al.,2007)的研究.北京MST雷達(dá)對(duì)中間層低熱層(MLT)區(qū)域探測(cè)的應(yīng)用與研究,將極大促進(jìn)對(duì)該地大氣層MLT區(qū)域的理解和認(rèn)識(shí).
本文主要從數(shù)據(jù)獲取率、與廊坊流星雷達(dá)測(cè)風(fēng)對(duì)比以及探測(cè)的風(fēng)場(chǎng)特征三個(gè)方面對(duì)北京MST雷達(dá)在MLT區(qū)域的探測(cè)能力、數(shù)據(jù)可靠性以及風(fēng)場(chǎng)觀測(cè)結(jié)果進(jìn)行初步分析.為今后應(yīng)用北京MST雷達(dá)開(kāi)展MLT區(qū)域研究奠定基礎(chǔ).
2北京MST雷達(dá)簡(jiǎn)介
北京MST雷達(dá)采用全固態(tài)數(shù)字陣列脈沖多普勒體制,由24×24付三單元八木天線組成矩形有源相控天線陣系統(tǒng).工作頻率:50±1 MHz;探測(cè)范圍:低模式為3.5~10 km,中模式為11~25 km,高模式為60~90 km及以上;徑向速度分辨率:≤0.2 m·s-1;時(shí)間分辨率≤30 min;最大測(cè)量徑向速度≥35 m·s-1,風(fēng)向?yàn)?~360°.本文應(yīng)用MST雷達(dá)高模式數(shù)據(jù),表2中給出具體的雷達(dá)高模式技術(shù)參數(shù).
3MST雷達(dá)探測(cè)中間層-低熱層(MLT)區(qū)域數(shù)據(jù)獲取率分析
數(shù)據(jù)獲取率是指在一段時(shí)間的探測(cè)中,探測(cè)數(shù)據(jù)通過(guò)質(zhì)量控制的次數(shù)與總探測(cè)次數(shù)的百分比,它是表征雷達(dá)探測(cè)能力的參數(shù)之一(胡明寶等,2008).
表2 北京MST雷達(dá)高模式技術(shù)參數(shù)
本研究中采用三倍標(biāo)準(zhǔn)差法對(duì)數(shù)據(jù)異常點(diǎn)進(jìn)行判別與剔除,有探測(cè)數(shù)據(jù)且通過(guò)異常點(diǎn)檢驗(yàn)的次數(shù)確定為探測(cè)數(shù)據(jù)通過(guò)質(zhì)量控制的次數(shù).
選取2012、2013年兩年北京MST雷達(dá)高模式數(shù)據(jù),統(tǒng)計(jì)兩年對(duì)應(yīng)月份的平均數(shù)據(jù)獲取率,繪制了各月份的平均數(shù)據(jù)獲取率時(shí)間-高度剖面圖,如圖1所示.
圖1中疊加的粉色直線分別是平均的日出和日落時(shí)間.可以看出數(shù)據(jù)獲取率有明顯的日變化特征:白天從65 km到100 km均可獲取數(shù)據(jù),數(shù)據(jù)獲取率的高值區(qū)主要集中在70~80 km高度區(qū)間內(nèi),最大值可達(dá)80%;夜間可獲取到數(shù)據(jù)的區(qū)域主要集中在80~100 km高度區(qū)間內(nèi),且數(shù)據(jù)獲取率在30%以下.
此外,還可以看出數(shù)據(jù)獲取率時(shí)間-高度剖面分布有季節(jié)變化特征,因此繪制MST雷達(dá)季節(jié)平均的數(shù)據(jù)獲取率時(shí)間-高度剖面,如圖2所示.
由圖2可見(jiàn)夏季白天可獲取數(shù)據(jù)的時(shí)間和高度區(qū)間均較其他季節(jié)大,春季次之,冬季最小;四個(gè)季節(jié)白天數(shù)據(jù)獲取率最大值均在70~80 km范圍內(nèi),正午前后達(dá)到最大.春季夜間在日落后到凌晨,90~100 km處數(shù)據(jù)獲取率較大,夏季白天以及日落后1 h可探測(cè)到120 km.
為進(jìn)一步分析白天和夜間數(shù)據(jù)獲取率隨高度的變化特征,將各季節(jié)白天、夜間的數(shù)據(jù)獲取率分別做平均,繪制各季節(jié)白天和夜間平均數(shù)據(jù)獲取率隨高度變化圖,如圖3所示,圖4疊加了標(biāo)準(zhǔn)差.
圖1 2012—2013年1—12月平均數(shù)據(jù)獲取率時(shí)間-高度剖面Fig.1 Time-height plot of monthly average data acquisition rate of 2012 and 2013
圖2 2012—2013年季節(jié)平均的數(shù)據(jù)獲取率時(shí)間-高度剖面Fig.2 Time-height plot of seasonal mean data acquisition rate of 2012 and 2013
圖3 各季節(jié)白天、夜間平均數(shù)據(jù)獲取率隨高度變化Fig.3 The profile of seasonal average data acquisition rate for daytime and nighttime
圖4 各季節(jié)白天、夜間平均數(shù)據(jù)獲取率隨高度變化疊加標(biāo)準(zhǔn)差圖Fig.4 The profile of seasonal average data acquisition rate for daytime and nighttime superimposed standard deviation
由圖3可以看出,各季節(jié)白天的平均數(shù)據(jù)獲取率隨高度變化均有2個(gè)峰值,第一峰值在70~80 km高度區(qū)間,且夏季的峰值最大,約為41%,其次是秋季38%,春季36%,冬季最小28%;第二峰值在90~100 km高度區(qū)間,夏季的峰值最大28%,春季次之19%,秋季和冬季最小15%.值得注意的是夏季白天低至60 km處數(shù)據(jù)獲取率仍達(dá)到5%,而其他季節(jié)均接近0.
各季節(jié)夜間在80~100 km高度區(qū)間可獲取數(shù)據(jù),基本都在92 km附近平均數(shù)據(jù)獲取率達(dá)到最大值,春季最大,約為19%,夏季次之15%,冬季11%,秋季10%.
圖4可以更加清楚地看到各季節(jié)白天平均數(shù)據(jù)獲取率隨高度變化的雙峰結(jié)構(gòu),以及夜間的單峰特征.各季節(jié)在白天數(shù)據(jù)獲取率最大值高度附近平均數(shù)據(jù)獲取率標(biāo)準(zhǔn)差大,夜間平均數(shù)據(jù)獲取率標(biāo)準(zhǔn)差較小.
綜合以上分析,MST雷達(dá)探測(cè)中間層-低熱層數(shù)據(jù)獲取率的變化特征與MST雷達(dá)的探測(cè)原理以及電離層自身的晝夜,季節(jié)變化密切相關(guān).因?yàn)樵?0 km以上,MST發(fā)射出電磁波的散射機(jī)制主要以熱散射為主,也就是電離層中電子密度不均勻形成的不規(guī)則散射體使得電磁波產(chǎn)生后向散射.
60~90 km是電離層的D層,主要的電離源是太陽(yáng)X射線,強(qiáng)烈地受太陽(yáng)直接輻射控制,該層的電子密度隨太陽(yáng)高度的變化而迅速變化,并隨太陽(yáng)高度變化呈現(xiàn)出較大的日變化特征.自由電子主要出現(xiàn)在白天,夜間基本消失,所以夜間在這個(gè)高度區(qū)間內(nèi)基本無(wú)雷達(dá)回波.90~130 km是電離層E層,自由電子白天夜間均存在,夜間電子密度減小但并不消失.所以夜間MST雷達(dá)在90 km附近可以獲取數(shù)據(jù),白天MST雷達(dá)在65~100 km均可獲取數(shù)據(jù),呈現(xiàn)出數(shù)據(jù)獲取率的日變化特征.
由于電離層D層、E層除日變化外,還有季節(jié)變化,因此MST雷達(dá)的數(shù)據(jù)獲取率也相應(yīng)地呈現(xiàn)出季節(jié)變化特征.
4北京MST雷達(dá)與廊坊流星雷達(dá)高空測(cè)風(fēng)對(duì)比分析
MST雷達(dá)是獲取大氣層各高度風(fēng)場(chǎng)的重要探測(cè)手段.北京MST雷達(dá)同樣可獲取風(fēng)場(chǎng)及其變化特征.為進(jìn)一步確定北京MST雷達(dá)探測(cè)的可靠性,將北京MST雷達(dá)測(cè)得的MLT區(qū)域經(jīng)向、緯向風(fēng)數(shù)據(jù)與附近廊坊流星雷達(dá)的測(cè)值作對(duì)比分析.廊坊流星雷達(dá)位于(39.4°N, 116.7°E),與北京MST雷達(dá)直線距離45 km左右.流星雷達(dá)工作頻率為35 MHz,可探測(cè)的高度范圍在70~110 km,時(shí)間分辨率為1 h,垂直分辨率2 km,每小時(shí)的數(shù)據(jù)是一天中有效觀測(cè)不少于18 h的平均值(Xiao et al.,2014).MST雷達(dá)高模式垂直分辨率1.12 km或1.13 km,時(shí)間分辨率小于等于30 min,本研究中采用的是時(shí)間分辨率30 min數(shù)據(jù).根據(jù)兩部雷達(dá)共同觀測(cè)的時(shí)間,采用2012年5月份數(shù)據(jù)進(jìn)行對(duì)比分析.
2012年5月流星雷達(dá)和MST雷達(dá)風(fēng)場(chǎng)數(shù)據(jù)的獲取率時(shí)空分布如圖5所示.
由圖5可見(jiàn),流星雷達(dá)在80~96 km的數(shù)據(jù)獲取率基本達(dá)到100%,而MST雷達(dá)的數(shù)據(jù)獲取率整體偏低,白天在70~80 km高度區(qū)間數(shù)據(jù)獲取率較高,在40%~70%之間.
流星雷達(dá)是利用流星(電離)余跡對(duì)無(wú)線電波的散射來(lái)實(shí)現(xiàn)對(duì)高空風(fēng)場(chǎng)探測(cè)的,流星現(xiàn)象白天夜間都會(huì)發(fā)生.每小時(shí)一次的風(fēng)場(chǎng)數(shù)據(jù)是時(shí)空上的平均值,因此流星雷達(dá)的風(fēng)場(chǎng)數(shù)據(jù)在時(shí)空分布上都比較連續(xù).而MST雷達(dá)是根據(jù)電離層中電子密度不均勻形成的不規(guī)則散射體使得電磁波產(chǎn)生后向散射從而獲取信息,電離層中電子密度分布的特征以及MST雷達(dá)瞬時(shí)的觀測(cè)決定了MST雷達(dá)獲取的數(shù)據(jù)比較離散而且數(shù)據(jù)獲取率較低.但另一方面,MST雷達(dá)在可測(cè)高度上時(shí)間與垂直分辨率較高.
80 km以下以及100 km以上流星雷達(dá)觀測(cè)到的流星數(shù)明顯減少,對(duì)測(cè)風(fēng)結(jié)果的準(zhǔn)確性與可靠性會(huì)有影響.流星雷達(dá)可用來(lái)連續(xù)監(jiān)測(cè)MLT區(qū)域的風(fēng)場(chǎng),但是他們的探測(cè)結(jié)果在80 km以上才是可靠的(Ratnam et al.,2001;Kumar et al.,2008).因此將MST雷達(dá)與流星雷達(dá)測(cè)風(fēng)對(duì)比的高度選為80~100 km.首先初步整體對(duì)比流星雷達(dá)與MST雷達(dá)高空風(fēng)場(chǎng)探測(cè)的情況,將二者2012年5月份的經(jīng)向和緯向風(fēng)速數(shù)據(jù)各自做平均,繪制時(shí)間-高度分布圖,如圖6所示.
由圖6可見(jiàn),80~100 km高度區(qū)間內(nèi),流星雷達(dá)與MST雷達(dá)經(jīng)向風(fēng)和緯向風(fēng)的時(shí)空分布特征基本一致.對(duì)于緯向風(fēng)來(lái)說(shuō),90 km以上基本為西風(fēng);流星雷達(dá)經(jīng)向風(fēng)隨時(shí)間推移南風(fēng)、北風(fēng)交替出現(xiàn),MST雷達(dá)經(jīng)向風(fēng)雖分布離散,但是仍表現(xiàn)出同樣的特征.
整體對(duì)比之后,進(jìn)一步將流星雷達(dá)數(shù)據(jù)按MST雷達(dá)觀測(cè)高度插值,同時(shí)將MST雷達(dá)每小時(shí)2次的觀測(cè)數(shù)據(jù)做平均,作為MST雷達(dá)每小時(shí)的數(shù)據(jù),將高度插值后的流星雷達(dá)廓線數(shù)據(jù)與MST雷達(dá)每小時(shí)廓線數(shù)據(jù)進(jìn)行對(duì)比分析,共有504組對(duì)比廓線.
首先繪制流星雷達(dá)與MST雷達(dá)經(jīng)向風(fēng)、緯向風(fēng)散點(diǎn)圖,并做最小二乘擬合,結(jié)果如圖7所示.
圖7中散點(diǎn)為不同高度上所有可對(duì)比點(diǎn),對(duì)比數(shù)據(jù)共有1849對(duì).緯向風(fēng)最小二乘擬合的斜率為0.93,截距為-7.7;相關(guān)系數(shù)為0.71,平均差值為-8.72 m·s-1,差值的標(biāo)準(zhǔn)差為26.94 m·s-1.經(jīng)向風(fēng)最小二乘擬合的斜率為0.87,截距為-0.5;相關(guān)系數(shù)為0.59,平均差值為0.89 m·s-1,差值的標(biāo)準(zhǔn)差為25.10 m·s-1.
圖5 流星雷達(dá)和MST雷達(dá)2012年5月份風(fēng)場(chǎng)數(shù)據(jù)的獲取率時(shí)空分布Fig.5 The time-height plot of wind dada acquisition rate for meteor radar and MST radar in May 2012
圖6 廊坊流星雷達(dá)與北京MST雷達(dá)2012年5月平均經(jīng)向風(fēng)與緯向風(fēng)時(shí)間-高度分布Fig.6 The time-height plot of monthly mean zonal wind and meridional wind for Langfang meteor radar and Beijing MST radar in May 2012
為進(jìn)一步了解MST雷達(dá)與流星雷達(dá)在各高度上風(fēng)場(chǎng)的對(duì)比情況,繪制了二者平均經(jīng)向風(fēng)、緯向風(fēng)隨高度變化,以及二者探測(cè)經(jīng)向風(fēng)、緯向風(fēng)平均差值隨高度變化圖,如圖8所示.
由圖8可以看出:緯向風(fēng)在80~87 km高度區(qū)間內(nèi),MST雷達(dá)與流星雷達(dá)探測(cè)的平均緯向風(fēng)均為東風(fēng),且MST雷達(dá)探測(cè)的風(fēng)速大于流星雷達(dá);87~100 km二者探測(cè)的平均緯向風(fēng)均為西風(fēng),且MST雷達(dá)探測(cè)的風(fēng)速小于流星雷達(dá).MST雷達(dá)與流星雷達(dá)探測(cè)緯向風(fēng)的差值在80~98 km區(qū)間內(nèi)基本保持在-10 m·s-1.
圖7 流星雷達(dá)與MST雷達(dá)2012年5月經(jīng)向風(fēng)、緯向風(fēng)散點(diǎn)圖以及最小二乘擬合結(jié)果,其中N為所有高度上總的對(duì)比點(diǎn)數(shù),r為相關(guān)系數(shù),d為MST雷達(dá)與流星雷達(dá)測(cè)風(fēng)的平均差值,std為二者差值的標(biāo)準(zhǔn)差 Fig.7 The scatter diagram and the least squares fitting results of zonal and meridional wind for Langfang meteor radar and Beijing MST radar in May 2012.N indicates the total number of comparison data pairs at all heights, r is correlation coefficient, d represents the mean difference between wind measurement of MST radar and meteor radar, std is the standard deviation of d
圖8 MST雷達(dá)與流星雷達(dá)探測(cè)的平均緯向、經(jīng)向風(fēng)以及二者的平均緯向、經(jīng)向風(fēng)差值隨高度分布Fig.8 The profile of mean zonal wind, mean meridional wind, mean zonal wind difference and the mean meridional wind difference for Langfang meteor radar and Beijing MST radar
經(jīng)向風(fēng)在82~100 km高度區(qū)間內(nèi),MST雷達(dá)與流星雷達(dá)探測(cè)的平均經(jīng)向風(fēng)均為北風(fēng),而且基本保持在-10 m·s-1;MST雷達(dá)與流星雷達(dá)探測(cè)經(jīng)向風(fēng)的差值在82~98 km區(qū)間內(nèi)最小,基本保持在±5 m·s-1以內(nèi).
綜合以上分析,可知MST雷達(dá)與流星雷達(dá)在80~100 km測(cè)風(fēng)結(jié)果在時(shí)空分布上有高度的一致性.二者探測(cè)的風(fēng)速數(shù)值有差異,這與兩部雷達(dá)處于不同位置,探測(cè)到的是不同經(jīng)緯度的高空風(fēng)場(chǎng)有關(guān).但是不同高度上二者探測(cè)的經(jīng)向風(fēng)、緯向風(fēng)差值基本保持不變,說(shuō)明MST雷達(dá)與流星雷達(dá)探測(cè)的一致性以及各自的系統(tǒng)穩(wěn)定性.同時(shí),以上結(jié)果也表明MST雷達(dá)探測(cè)得到的數(shù)據(jù)是可靠的.
此外,MST雷達(dá)白天在65~80 km數(shù)據(jù)獲取率明顯高于80~100 km,可對(duì)電離層D層中性風(fēng)場(chǎng),湍流信息等進(jìn)行有效而高時(shí)空分辨率的探測(cè).
5北京MST雷達(dá)探測(cè)的水平風(fēng)月平均結(jié)果初步分析
在已獲知MST雷達(dá)在中間層-低熱層MLT區(qū)域的探測(cè)能力以及初步了解其獲取數(shù)據(jù)可靠的基礎(chǔ)上,本文分析了雷達(dá)探測(cè)的水平風(fēng)的時(shí)空分布特征.將MST雷達(dá)2012年、2013年以及水平風(fēng)模式HWM07月平均的經(jīng)向、緯向風(fēng)數(shù)據(jù)做初步對(duì)比.繪制各月平均的緯向風(fēng)、經(jīng)向風(fēng)時(shí)空剖面圖,分別見(jiàn)圖9和圖10.HWM07模式應(yīng)用50多年的衛(wèi)星資料,火箭探空, FPI,非相干散射雷達(dá),中頻雷達(dá),測(cè)風(fēng)測(cè)溫激光雷達(dá)以及低層大氣的NCEP數(shù)據(jù)等,為地面到外逸層(0~500 km)范圍內(nèi)水平風(fēng)場(chǎng)提供了一個(gè)統(tǒng)計(jì)表示(Drob et al.,2008).
圖9 2012年(a)、2013年(b)、HWM07(c)緯向風(fēng)月平均時(shí)間-高度剖面Fig.9 The time-height plot of monthly mean zonal wind for 2012 (a), 2013 (b), HWM07 (c)
圖10 2012年(a)、2013年(b)、HWM07(c)經(jīng)向風(fēng)月平均時(shí)間-高度剖面Fig.10 The time-height plot of monthly meridional wind for 2012(a), 2013(b), HWM07(c)
由圖9可以看到,2012和2013年MST雷達(dá)各月平均的緯向風(fēng)時(shí)空分布特征有較高的一致性,而且與HWM07模式結(jié)果也基本一致.具體的時(shí)空分布特征為:1—3月60~90 km基本都為西風(fēng),90~100 km午后為東風(fēng);4月表現(xiàn)出過(guò)渡特征,70~90 km為東風(fēng);5—7月60~85 km為東風(fēng),85~100 km為西風(fēng),且風(fēng)速均較大;8月也表現(xiàn)出過(guò)渡特征,雖然風(fēng)向特征與5—7月一致,但是風(fēng)速明顯減小,尤其是60~80 km;9、10月在90~100 km一天中東風(fēng)西風(fēng)交替出現(xiàn),60~80 km為西風(fēng),且10月風(fēng)速大于9月;11、12月60~90 km為西風(fēng),90~100 km午后為東風(fēng).
由圖10可以看到,經(jīng)向風(fēng)月平均時(shí)空分布特征比較復(fù)雜,但是2012和2013年MST雷達(dá)各月平均的經(jīng)向風(fēng)時(shí)空分布特征基本一致,而且與HWM07模式結(jié)果也基本一致.具體的時(shí)空分布特征可以參見(jiàn)圖中所示.圖9和10的結(jié)果進(jìn)一步表明了北京MST雷達(dá)探測(cè)中間層-低熱層(MLT)區(qū)域數(shù)據(jù)的可靠性.
6結(jié)論
利用子午工程北京MST雷達(dá)2012、2013年兩年的高模式觀測(cè)數(shù)據(jù),進(jìn)行了數(shù)據(jù)質(zhì)量檢驗(yàn),在此基礎(chǔ)上做了以下初步分析,其主要結(jié)論如下:
(1) 數(shù)據(jù)獲取率的日變化與季節(jié)變化特征.分析了中間層-低熱層(MLT)60~120 km區(qū)域數(shù)據(jù)獲取率,可以看到明顯日變化特征:白天從65 km到100 km均可獲取數(shù)據(jù),數(shù)據(jù)獲取率的高值區(qū)主要集中在70~80 km高度區(qū)間內(nèi),最大值可達(dá)80%;夜間可獲取到數(shù)據(jù)的區(qū)域主要集中在80~100 km高度區(qū)間內(nèi),數(shù)據(jù)獲取率在30%以下.結(jié)果表明,該MST雷達(dá)白天可以探測(cè)到電離層D層和E層低層,夜間D層消失,只探測(cè)到E層低層.此外,數(shù)據(jù)獲取率還有明顯的季節(jié)變化特征:夏季白天可獲取數(shù)據(jù)的時(shí)間和高度區(qū)間都比較大,春季次之,冬季最小.夏季白天以及日落后1h可探測(cè)到120 km.
(2) MST雷達(dá)與流星雷達(dá)測(cè)風(fēng)對(duì)比分析.對(duì)北京MST雷達(dá)與廊坊流星雷達(dá)2012年5月份,80~100 km高度區(qū)間測(cè)量的水平風(fēng)進(jìn)行對(duì)比分析,測(cè)風(fēng)結(jié)果在時(shí)空分布上有較高的一致性.二者探測(cè)的風(fēng)速數(shù)值有差異,這與兩部雷達(dá)相距45 km,探測(cè)到的是不同位置的高空風(fēng)場(chǎng)有關(guān),但是不同高度上二者探測(cè)的經(jīng)向風(fēng)、緯向風(fēng)差值基本保持不變.以上均表明MST雷達(dá)探測(cè)得到的數(shù)據(jù)是可靠的.
(3) MLT區(qū)域60~120 km MST雷達(dá)觀測(cè)水平風(fēng)場(chǎng)的月平均時(shí)空分布特征.2012年和2013年兩年相應(yīng)月份緯向風(fēng)、經(jīng)向風(fēng)的時(shí)空分布特征有高的一致性,此外與HWM07模式對(duì)比也基本一致,進(jìn)一步表明MST雷達(dá)探測(cè)數(shù)據(jù)的可靠性.
上述初步結(jié)果表明,北京MST雷達(dá)對(duì)中間層-低熱層60~120 km高度區(qū)域已具備好的探測(cè)能力,而且探測(cè)資料是可靠的.下一步將分析觀測(cè)資料所揭示的該區(qū)域的動(dòng)力學(xué)特征與MLT結(jié)構(gòu)變化,闡明太陽(yáng)輻射與活動(dòng)性以及大氣層環(huán)流的特征與作用.
致謝本項(xiàng)成果使用國(guó)家重大科技基礎(chǔ)設(shè)施子午工程科學(xué)數(shù)據(jù).北京MST雷達(dá)是子午工程建設(shè)的一項(xiàng)重大設(shè)備.中國(guó)科學(xué)院大氣研究所香河觀測(cè)站MST雷達(dá)觀測(cè)組負(fù)責(zé)長(zhǎng)期觀測(cè).王勇、章文星、段樹(shù)三位正研級(jí)高工承擔(dān)了雷達(dá)研發(fā)和改進(jìn)的技術(shù).廊坊流星雷達(dá)2012年5月份的數(shù)據(jù)是中國(guó)科學(xué)院空間中心提供的,在此一并表示感謝.田玉芳感謝凌超在HWM07模式方面的相關(guān)幫助.
References
Balsley B B, Ecklund W L, Carter D A, et al. 1980. The MST radar at Poker-Flat, Alaska.RadioScience, 15(2): 213-223, doi: 10.1029/RS015i002p00213. Belova E, Kirkwood S, Rao T N, et al. 2012. Spectral characteristics and scatter cross-section of low latitude mesospheric echoes measured by the Indian MST radar at Gadanki.AnnalesGeophysicae, 30(6): 983-990, doi: 10.5194/angeo-30-983-2012. Chandra H, Sinha H S S, Das U, et al. 2008. First mesospheric turbulence study using coordinated rocket and MST radar measurements over Indian low latitude region.AnnalesGeophysicae, 26(9): 2725-2738. Chandra H, Sinha H S S, Patra A K, et al. 2012. Low-latitude mesospheric turbulence investigated using coordinated MST radar and rocket-borne observations from India.JournalofGeophysicalResearch-Atmospheres, 117(D22), doi: 10.1029/2011jd016946.
Czechowsky P, Schmidt G, Rüster R. 1984. The mobile SOUSY Doppler radar: technical design and first results.RadioScience, 19(1): 441-450, doi: 10.1029/RS019i001p00441.
Drob D P, Emmert J T, Crowley G, et al. 2008. An empirical model of the Earth′s horizontal wind fields: HWM07.J.Geophys.Res., 113: A12304, doi: 10.1029/2008JA013668. Eswaraiah S, Venkat Ratnam M, Krishna Murthy B V, et al. 2011. Low-latitude mesospheric vertical winds observed using VHF radar.JournalofGeophysicalResearch:Atmospheres, 116: D22117, doi: 10.1029/2011jd016385.
Eswaraiah S, Ratnam M V, Portnyagin Y I, et al. 2012. Comparison of mesospheric mean vertical winds between MST radar and empirical model over a tropical station.CurrentScience, 102(9): 1298-1303.
Fukao S, Yamanaka M D, Ao N, et al. 1994. Seasonal variability of vertical eddy diffusivity in the middle atmosphere. 1. Three-year observations by the middle and upper atmosphere radar.JournalofGeophysicalResearch-Atmospheres, 99(D9): 18973-18987, doi: 10.1029/94jd00911.
Fukao S. 2007. Recent advances in atmospheric radar study.J.Meteor.Soc.Jpn., 85B: 215-239, doi: 10.2151/jmsj.85B.215. Gage K S, Balsley B B. 1980. On the scattering and reflection mechanisms contributing to clear air radar echoes from the troposphere, stratosphere, and mesophere.RadioScience, 15(2): 243-257, doi: 10.1029/RS015i002p00243.
Guo L, Lehmacher G A, Kudeki E, et al. 2007. Turbulent kinetic energy dissipation rates and eddy diffusivities in the tropical mesosphere using Jicamarca radar data.AdvancesinSpaceResearch, 40(6): 744-750, doi: 10.1016/j.asr.2007.05.068.
Hocking W K. 1997. Recent advances in radar instrumentation and techniques for studies of the mesosphere, stratosphere, and troposphere.RadioScience, 32(6): 2241-2270, doi: 10.1029/97rs02781.
Hooper D, Thomas L. 1995. Aspect sensitivity of VHF scatterers in the troposphere and stratosphere from comparisons of powers in off-vertical beams.JournalofAtmosphericandTerrestrialPhysics, 57(6): 655-663, doi: 10.1016/0021-9169(94)00056-T. Hu M B, Zheng G G, Xiao W J. 2008. The statistical analysis of wind profiler radar data availability.ModernRadar(in Chinese), 30(10): 14-16.Kamala S, Rao D N, Chakravarty S C, et al. 2003. Vertical structure of mesospheric echoes from the Indian MST radar.JournalofAtmosphericandSolar-TerrestrialPhysics, 65(1): 71-83, doi: 10.1016/s1364-6826(02)00281-x.
Kato S, Ogawa T, Tsuda T, et al. 1984. The middle and upper-atmosphere radar: First results using a partial system.RadioScience, 19(6): 1475-1484, doi: 10.1029/RS019i006p01475.
Kubo K, Sugiyama T, Nakamura T, et al. 1997. Seasonal and interannual variability of mesospheric echoes observed with the middle and upper atmosphere radar during 1986-1995.Geophys.Res.Lett., 24(10): 1211-1214, doi: 10.1029/97gl01063.
Kumar G K, Ratnam M V, Patra A K, et al. 2007. Climatology of low-latitude mesospheric echo characteristics observed by Indian mesosphere, stratosphere, and troposphere radar.JournalofGeophysicalResearch, 112(D6): doi: 10.1029/2006jd007609.
Kumar G K, Ratnam M V, Patra A K, et al. 2008. Low-latitude mesospheric mean winds observed by Gadanki mesosphere-stratosphere-troposphere (MST) radar and comparison with rocket, High Resolution Doppler Imager (HRDI), and MF radar measurements and HWM93.JournalofGeophysicalResearch-Atmospheres, 113(D19): doi: 10.1029/2008jd009862. Kumar G K, Kumar K K, Singer W, et al. 2014. Mesosphere and lower thermosphere zonal wind variations over low latitudes: Relation to local stratospheric zonal winds and global circulation anomalies.JournalofGeophysicalResearch-Atmospheres, 119(10): 5913-5927, doi: 10.1002/2014jd021610.
Lehmacher G A, Guo L, Kudeki E, et al. 2007. High-resolution observations of mesospheric layers with the Jicamarca VHF radar.AdvancesinSpaceResearch, 40(6): 734-743, doi: 10.1016/j.asr.2007.05.059.
Lehmacher G A, Kudeki E, Akgiray A, et al. 2009. Radar cross sections for mesospheric echoes at Jicamarca.AnnalesGeophysicae, 27(7): 2675-2684.Li W, Lü D R. 1998. Preliminary results of atmospheric turbulence observation with Xianghe MST/ST radar.ActaMeteorologicaSinica, 12(4): 394-399.
Lü D R, Chen H B. 2003. Advances in middle atmosphere physics research.ChineseJournalofAtmosphericSciences(in Chinese), 27(4): 750-769.Rao D N, Ratnam M V, Rao T N, et al. 2001. Seasonal variation of vertical eddy diffusivity in the troposphere, lower stratosphere and mesosphere over a tropical station.AnnalesGeophysicae, 19(8): 975-984.
Rao P B, Jain A R, Kishore P, et al. 1995. Indian MST radar. 1. System description and sample vector wind measurements in st mode.RadioScience, 30(4): 1125-1138, doi: 10.1029/95rs00787. Rao S V B, Eswaraiah S, Ratnam M V, et al. 2014. Advanced meteor radar installed at Tirupati: System details and comparison with different radars.JournalofGeophysicalResearch-Atmospheres, 119(21): 11893-11904, doi: 10.1002/2014jd021781.Ratnam M V, Rao D N, Rao T N, et al. 2001. Mean winds observed with Indian MST radar over tropical mesosphere and comparison
with various techniques.AnnalesGeophysicae, 19(8): 1027-1038.
Ratnam M V, Kumar G K, Murthy B V K, et al. 2008. Long-term variability of the low latitude mesospheric SAO and QBO and their relation with stratospheric QBO.Geophys.Res.Lett., 35(21): doi: 10.1029/2008gl035390.
R?ttger J, Rastogi P K, Woodman R F. 1979. High-resolution VHF radar observations of turbulence structures in the mesosphere.Geophys.Res.Lett., 6(7): 617-620, doi: 10.1029/GL006i007p00617.R?ttger J, Liu C H, Chao J K, et al. 1990. The Chung-Li VHF radar: Technical layout and a summary of initial results.RadioScience, 25(4): 487-502, doi: 10.1029/RS025i004p00487.
Sato T, Tsuda T, Kato S, et al. 1985. High-resolution MST observations of turbulence by using the MU radar.RadioScience, 20(6): 1452-1460, doi: 10.1029/RS020i006p01452.
Selvaraj D, Patra A K, Chandra H, et al. 2014. Scattering cross section of mesospheric echoes and turbulence parameters from Gadanki radar observations.JournalofAtmosphericandSolar-TerrestrialPhysics, 119: 162-172, doi: 10.1016/j.jastp.2014.08.004.
Van Zandt T E. 2000. A brief history of the development of wind-profiling or MST radars.AnnalesGeophysicae, 18(7): 740-749. Woodman R F, Guillen A. 1974. Radar observations of winds and turbulence in the stratosphere and mesosphere.JournaloftheAtmosphericSciences, 31(2): 493-505, doi: 10.1175/1520-0469(1974)031〈0493:roowat〉2.0.co;2.
Xiao C Y, Hu X, Xu Q C. 2014. Response of mesosphere and lower thermosphere wind over mid-latitude to the 2013 major sudden stratospheric warming event.∥ Proceedings of the 2014 XXXIth URSI General Assembly and Scientific Symposium (URSI GASS). Beijing: IEEE, 1-4, doi: 10.1109/URSIGASS.2014.6929843.
Zhao Z Y, Zhou C, Qing H Y, et al. 2013. Wuhan Atmosphere Radio Exploration (WARE) radar: System design and online winds measurements.RadioScience, 48(3): 326-333, doi: 10.1002/rds.20040.
附中文參考文獻(xiàn)
胡明寶, 鄭國(guó)光, 肖文建. 2008. 風(fēng)廓線雷達(dá)數(shù)據(jù)獲取率的統(tǒng)計(jì)分析. 現(xiàn)代雷達(dá), 30(10): 14-16.
李衛(wèi), 呂達(dá)仁. 1997. 香河MST雷達(dá)大氣湍流探測(cè)初步結(jié)果.∥ 地球環(huán)境和氣候變化探測(cè)與過(guò)程研究. 北京:氣象出版社.
呂達(dá)仁, 陳洪濱. 2003. 平流層和中層大氣研究的進(jìn)展. 大氣科學(xué), 27(4): 750-769.
(本文編輯胡素芳)
基金項(xiàng)目國(guó)家重大科研儀器設(shè)備研制專(zhuān)項(xiàng)“多波段多大氣成分主被動(dòng)綜合探測(cè)系統(tǒng)”(41127901)(NSFC)資助.
作者簡(jiǎn)介田玉芳,女,1989年生,博士生,主要從事晴空大氣探測(cè)雷達(dá)探測(cè)原理與應(yīng)用研究. E-mail:tianyufang@mail.iap.ac.cn *通信作者呂達(dá)仁,男,1940年生,研究員,長(zhǎng)期從事大氣科學(xué)、日地科學(xué)與空間對(duì)地觀測(cè)基礎(chǔ)與重大應(yīng)用研究. E-mail:ludr@mail.iap.ac.cn
doi:10.6038/cjg20160204 中圖分類(lèi)號(hào)P412
收稿日期2015-09-06,2015-12-14收修定稿
Preliminary analysis of Beijing MST radar observation results in the mesosphere-lower thermosphere
TIAN Yu-Fang1,2, Lü Da-Ren1*
1KeyLaboratoryofMiddleAtmosphereandGlobalEnvironmentObservation(LAGEO),InstituteofAtmosphericPhysics,ChineseAcademyofSciences,Beijing100029,China2UniversityofChineseAcademyofSciences,Beijing100049,China
AbstractBeijing MST radar is one of the Chinese Meridian Project constructions of only two domestic MST radars. It was built by Institute of Atmospheric Physics, Chinese Academy of Sciences, located at Xianghe field observatory,Hebei province (39.754°N,116.962°E). System integrated test was finished in May 2011 and completed acceptance in the end of July 2011. After 1-month joint tests from 21 September 2011 to 20 October 2011, Beijing MST radar operated continuously to observe real-time winds and turbulence in the altitude range of 3.5~10 km (low mode), 11~25 km (middle mode) and 60~90 km (even above) (high mode). Beijing MST radar operates at 50 MHz with a peak power of 172kW and 104 m2 of antenna area. The purpose of this paper is to investigate the detection capability and data reliability of Beijing MST radar in the mesosphere-lower thermosphere (MLT) region.
KeywordsBeijing MST radar; Mesosphere-lower thermosphere; Detection capability; Data reliability; Spatial and temporal distribution of horizontal wind
田玉芳, 呂達(dá)仁. 2016. 北京MST雷達(dá)探測(cè)中間層-低熱層觀測(cè)結(jié)果初步分析.地球物理學(xué)報(bào),59(2):440-452,doi:10.6038/cjg20160204.
Tian Y F, Lü D R. 2016. Preliminary analysis of Beijing MST radar observation results in the mesosphere-lower thermosphere.ChineseJ.Geophys. (in Chinese),59(2):440-452,doi:10.6038/cjg20160204.