劉勇，陳一定*，劉立波

1 中國(guó)科學(xué)院地質(zhì)與地球物理研究所中國(guó)科學(xué)院地球與行星物理重點(diǎn)實(shí)驗(yàn)室，北京 1000292 中國(guó)科學(xué)院地質(zhì)與地球物理研究所北京空間環(huán)境國(guó)家野外科學(xué)觀測(cè)研究站，北京 1000293 中國(guó)科學(xué)院大學(xué)，北京 100049

?

電離層春秋分不對(duì)稱的地方時(shí)依賴

劉勇1,2,3，陳一定1,2,3*，劉立波1,2

1 中國(guó)科學(xué)院地質(zhì)與地球物理研究所中國(guó)科學(xué)院地球與行星物理重點(diǎn)實(shí)驗(yàn)室，北京 1000292 中國(guó)科學(xué)院地質(zhì)與地球物理研究所北京空間環(huán)境國(guó)家野外科學(xué)觀測(cè)研究站，北京 1000293 中國(guó)科學(xué)院大學(xué)，北京 100049

利用全球203個(gè)電離層測(cè)高儀臺(tái)站的F2層臨界頻率(foF2)和E層臨界頻率(foE)，以及美國(guó)噴氣推進(jìn)實(shí)驗(yàn)室(JPL)提供的電離層總電子含量(TEC)地圖數(shù)據(jù)統(tǒng)計(jì)分析了電離層春秋分(March Equinox and September Equinox，ME and SE)不對(duì)稱的特點(diǎn)．基于電離層參量隨年積日(Day of Year，DoY)和太陽(yáng)活動(dòng)指數(shù)F10.7變化的傅里葉級(jí)數(shù)模型，對(duì)foF2、foE及TEC數(shù)據(jù)分別進(jìn)行最小二乘法擬合，將電離層參量歸算到低太陽(yáng)活動(dòng)(F10.7=80)、中等太陽(yáng)活動(dòng)(F10.7=150)和高太陽(yáng)活動(dòng)(F10.7=200)水平．該方法定量分離了實(shí)際觀測(cè)數(shù)據(jù)中包含的電離層參量隨季節(jié)和太陽(yáng)活動(dòng)的變化，因而得到了更為定量、精確的電離層春秋分不對(duì)稱性特征.分析了不同地方時(shí)(LT)的春秋分不對(duì)稱性指數(shù)(Asymmetry Index，AI)和春秋分差值Δ(=ME-SE)的全球分布特征與太陽(yáng)活動(dòng)依賴性．結(jié)果表明，foE日出時(shí)全球主要表現(xiàn)為9月分點(diǎn)值高于3月分點(diǎn)值，午后春秋分不對(duì)稱性幾乎消失，而日落時(shí)則反轉(zhuǎn)為3月分點(diǎn)值高于9月分點(diǎn)值；foF2日出時(shí)除少數(shù)地區(qū)外也主要表現(xiàn)為9月分點(diǎn)值高于3月分點(diǎn)值，而其他時(shí)段則相反；TEC日出時(shí)低太陽(yáng)活動(dòng)時(shí)的全球及中高太陽(yáng)活動(dòng)時(shí)的低緯地區(qū)表現(xiàn)為9月分點(diǎn)值高于3月分點(diǎn)值，而其他時(shí)段則相反．foE春秋分不對(duì)稱性受太陽(yáng)活動(dòng)影響較弱，而foF2和TEC的春秋分不對(duì)稱隨太陽(yáng)活動(dòng)有明顯的變化，其3月分點(diǎn)值相對(duì)于9月分點(diǎn)值增加．計(jì)算了F2層峰高(hmF2)處對(duì)應(yīng)的氧氮濃度比([O]/[N2]，由大氣模型NRLMSISE-00計(jì)算得到)和hmF2的春秋分不對(duì)稱性，提取了TEC年變化的幅度及相位信息．氧氮濃度比和hmF2的春秋分不對(duì)稱性能夠部分解釋電離層的春秋分不對(duì)稱性，而TEC春秋分不對(duì)稱的全球分布特征可以用TEC年變化的相位的全球分布解釋．

F2層臨界頻率；E層臨界頻率；GPS TEC；電離層；春秋分不對(duì)稱

1 引言

受太陽(yáng)輻射、中性大氣等因素變化的影響，電離層表現(xiàn)出從氣候?qū)W到天氣學(xué)多重時(shí)間尺度的復(fù)雜變化．就氣候?qū)W變化而言，主要包括太陽(yáng)活動(dòng)周變化、季節(jié)變化及周日變化等．對(duì)于電離層季節(jié)變化，如年異常(從全球范圍看，NmF2的年度變化12月比6月大)、半年異常(電子含量在兩分點(diǎn)月份附近出現(xiàn)峰值)、冬季異常(中緯度地區(qū)白天NmF2冬季高于夏季)等，其原因可能與中性大氣成分、熱層中性風(fēng)、日地距離的季節(jié)變化等有關(guān)(Richards,2001;Rishbeth et al.,2000;Buonsanto,1986)．春秋分不對(duì)稱性是電離層季節(jié)變化的一個(gè)重要特征(Balan et al.,1998;Chen et al.,2012)．太陽(yáng)天頂角在春分點(diǎn)與秋分點(diǎn)(本文所討論春分點(diǎn)與秋分點(diǎn)分別指3月分點(diǎn)與9月分點(diǎn))大致相當(dāng)，但即使在太陽(yáng)活動(dòng)水平一致時(shí)，電離層參量有時(shí)在這二分點(diǎn)表現(xiàn)出較強(qiáng)的差異(Liu et al.,2010)，即電離層的春秋分不對(duì)稱性．這種特征可能由中性大氣的季節(jié)變化以及電離層動(dòng)力學(xué)過程的春秋分差異引起(Kawamura et al.,2002;Aruliah et al.,1996)．然而現(xiàn)有的一些電離層經(jīng)驗(yàn)?zāi)Ｊ?，如?guó)際參考電離層(International Reference Ionosphere)(Bilitza et al.,2014)、Fejer的電場(chǎng)模型(Fejer and Scherliess,1997)，都沒有考慮電離層的春秋分不對(duì)稱性的影響．研究該特征對(duì)認(rèn)識(shí)、理解電離層季節(jié)變化以及對(duì)電離層經(jīng)驗(yàn)建模都具有重要意義．

在電離層季節(jié)變化期間，太陽(yáng)輻射也在顯著改變(除太陽(yáng)活動(dòng)低年)，引起電離層變化．從觀測(cè)數(shù)據(jù)中通常難以準(zhǔn)確地將這兩種因素導(dǎo)致的電離層變化完全定量分離開，從而得到電離層春秋分不對(duì)稱特征的更為定量的信息，有時(shí)甚至?xí)鹫`導(dǎo)．例如，Unnikrishnan等(Unnikrishnan et al.,2002)分析Palehua地區(qū)TEC的春秋分不對(duì)稱特征，得出TEC春秋分不對(duì)稱在太陽(yáng)活動(dòng)高年與太陽(yáng)活動(dòng)低年截然相反(1981太陽(yáng)活動(dòng)高年9月份TEC高，而1984太陽(yáng)活動(dòng)低年3月份TEC高)的結(jié)論．從太陽(yáng)輻射變化的角度分析原因，可以發(fā)現(xiàn)，1981年9月份附近的太陽(yáng)輻射水平明顯強(qiáng)于3月份附近的太陽(yáng)輻射水平，從而導(dǎo)致TEC值在9月份高于3月份．由此可見，對(duì)于電離層季節(jié)變化研究，分離出太陽(yáng)輻射變化的影響十分必要．在本文中，為了得到不同太陽(yáng)活動(dòng)條件下的電離層季節(jié)變化，以分析其春秋分不對(duì)稱性，我們通過建立電離層參量(foF2、foE、TEC)隨太陽(yáng)活動(dòng)指數(shù)F10.7和年積日DoY變化的模型，剝離觀測(cè)到的電離層季節(jié)變化中的太陽(yáng)輻射變化的影響；利用該模型計(jì)算低太陽(yáng)活動(dòng)(F10.7=80)、中等太陽(yáng)活動(dòng)(F10.7=150)和高太陽(yáng)活動(dòng)(F10.7=200)情況下電離層參量的季節(jié)變化，得到其春秋分不對(duì)稱性指數(shù)(AIfoF2、AIfoE、AITEC)和春秋分差值(ΔfoF2、ΔfoE、ΔTEC)，從而分析不同地方時(shí)的春秋分不對(duì)稱指數(shù)和春秋分差值在全球的分布特點(diǎn)及其隨太陽(yáng)活動(dòng)的變化特征，最后簡(jiǎn)要探討出現(xiàn)這些特征的原因．

2 數(shù)據(jù)及分析方法

所用的全球203個(gè)電離層測(cè)高儀臺(tái)站foF2及foE的每小時(shí)值數(shù)據(jù)下載自SPIDR網(wǎng)站(http:∥spidr.ngdc.noaa.gov/spidr)．為得到可靠的foF2及foE數(shù)據(jù)，選擇經(jīng)手動(dòng)標(biāo)定及編輯過的數(shù)據(jù)(Pezzopane and Scotto,2007)．?dāng)?shù)據(jù)空間分布和時(shí)間分布如圖1所示，數(shù)據(jù)時(shí)間范圍依不同的臺(tái)站而有所不同，大致為1957年至2014年，或其間的若干年．TEC數(shù)據(jù)由美國(guó)噴氣推進(jìn)實(shí)驗(yàn)室(JPL)提供(下載自ftp:∥cddis.gsfc.nasa.gov/pub/gps/products/ionex/)，時(shí)間分辨率為每?jī)蒛T時(shí)一張全球TEC圖，空間分辨率為2.5°(Latitude)×5°(Longitude)．為了研究TEC春秋分不對(duì)稱性對(duì)地方時(shí)的變化規(guī)律，將TEC數(shù)據(jù)按地方時(shí)線性插值到地方時(shí)網(wǎng)格．太陽(yáng)活動(dòng)指數(shù)F10.7由加拿大彭蒂克頓無(wú)線電觀測(cè)站(Pentictin Radio Observatory)提供，每天一個(gè)數(shù)據(jù)．大氣模型NRLMSISE-00是一個(gè)中性大氣的經(jīng)驗(yàn)?zāi)Ｐ?，提供從地表到逃逸層的中性組分、密度和溫度(Picone et al.,2002)，本文中所用到的氧氮濃度比數(shù)據(jù)由該模型計(jì)算得到．考慮到地磁活動(dòng)的影響，同時(shí)保留較多的數(shù)據(jù)，將Ap>27的數(shù)據(jù)去除．

圖1 測(cè)高儀數(shù)據(jù)分布(a)臺(tái)站位置分布；(b)數(shù)據(jù)時(shí)間覆蓋.Fig.1 Spatial and temporal distributions of ionosonde data(a) Lonosonde locations;(b) Temporal coverage of the data.

眾多研究表明，可以用二次多項(xiàng)式函數(shù)較好地描述foF2隨F10.7的變化關(guān)系(Liu et al.,2004,2006；Chen et al.,2008)．因此建立foF2與F10.7的關(guān)系如下：

(1)

其中A，B，C為系數(shù)．而foF2的季節(jié)變化成分通常主要包括年變化分量與半年變化分量，可以用傅里葉級(jí)數(shù)表示如下(Liu et al.,2009;Yang et al.,2012)：

(2)

其中c0為年平均值，DoY為年積日，ε為略去殘差．將年變化分量中的兩項(xiàng)與半年變化分量中的兩項(xiàng)分別合并，得到

(3)

(4)

其中，A1、A2分別為年變化和半年變化的振幅，φ1、φ2分別為年變化和變年變化的相位．

由此，我們對(duì)每個(gè)站每個(gè)地方時(shí)的foF2建立如下模型(Liu et al.,2009)：

(5)

其中Ai,Bi,Ci,Di,Ei,Fi(i=0，1，2)為系數(shù)．基于該模型，對(duì)觀測(cè)的foF2值進(jìn)行最小二乘擬合，確定模型系數(shù)，進(jìn)而通過模型計(jì)算不同太陽(yáng)活動(dòng)水平(不同F(xiàn)10.7值)條件下foF2隨年積日的變化，分析春秋分不對(duì)稱特征．

以日本Akita站(臺(tái)站代碼：AK539，位置：39.7°N，140.1°E)地方時(shí)14∶00為例，說明本文中的數(shù)據(jù)分析方法．對(duì)該臺(tái)站地方時(shí)為14∶00的foF2數(shù)據(jù)建立式(5)所示模型．將觀測(cè)數(shù)據(jù)foF2及其相應(yīng)DoY、F10.7代入模型，通過最小二乘擬合確定模型系數(shù)Ai,Bi,Ci,Di,Ei,Fi(i=0，1，2)．?dāng)M合的結(jié)果與觀測(cè)值的比較如圖2所示．從圖2中可以看出觀測(cè)值和擬合值比較一致，可以清晰地看到AK539站地方時(shí)14∶00的半年異常特征．模型值與觀測(cè)值的相對(duì)誤差(RE)及其期望(μ)、均方根誤差(RMSE)按下式計(jì)算：

(6)

(7)

(8)

其中foF2obs為觀測(cè)值，foF2mod為模型值，REi為每個(gè)觀測(cè)點(diǎn)與相應(yīng)模型值的相對(duì)誤差，RMSE為RE的標(biāo)準(zhǔn)差．為保證模型擬合的可靠性，我們將均方根誤差控制在30%以下，否則認(rèn)為該組數(shù)據(jù)擬合得不好，將該組數(shù)據(jù)剔除． AK539站地方時(shí)14∶00擬合的相對(duì)誤差分布如圖3所示，其擬合相對(duì)誤差的期望μ為0.0486%，均方根誤差為12.2191%．圖4 為所有臺(tái)站各地方時(shí)foF2擬合的μ及RMSE分布．去除RMSE>30%的數(shù)據(jù)后，所有數(shù)據(jù)擬合的μ的平均值為0.0207%，RMSE平均值為15.55%．同理，foE和TEC分別去除RMSE>20%和RMSE>35%的數(shù)據(jù)后，所有數(shù)據(jù)擬合的μ的平均值分別為-0.0937%、0.1465%，RMSE平均值分別為6.41%、18.87%．

定義春秋分不對(duì)稱性指數(shù)AIfoF2為foF2在3月分點(diǎn)與9月分點(diǎn)的差值與foF29月分點(diǎn)的值之比，再乘以100%，即foF2兩分點(diǎn)值的相對(duì)差：

(9)

foF2春秋分絕對(duì)差ΔfoF2為

(10)

圖2 AK539站地方時(shí)14∶00 foF2的擬合結(jié)果與觀測(cè)數(shù)據(jù)的比較(a) 原始數(shù)據(jù)與擬合數(shù)據(jù)的時(shí)間序列；(b) 擬合foF2隨DoY與F10.7的變化；(c) 觀測(cè)的foF2隨DoY與F10.7的變化.Fig.2 Comparison between fitted and observed foF2 at 14∶00 LT over Akita ionosonde station(a) Temporal series of observed and fitted foF2;(b) Fitted foF2 as a function of DoY and F10.7; (c) Scatter plots of observed foF2 versus DoY and F10.7.

同理，foE、TEC的春秋分不對(duì)稱性指數(shù)AIfoE、AITEC及其春秋分絕對(duì)差值ΔfoE、ΔTEC分別為

(11)

(12)

(13)

ΔTEC=TECME-TECSE.

(14)

由此，定量得到了Akita站地方時(shí)為14∶00在不同太陽(yáng)活動(dòng)水平下的春秋分不對(duì)稱性特征．由foF2不對(duì)稱性指數(shù)的定義可以看出,若AIfoF2的值為正則表示foF23月分點(diǎn)值大于9月分點(diǎn)值，若AIfoF2的值為負(fù)則表示foF23月分點(diǎn)值小于9月分點(diǎn)值，若AIfoF2的值為零，則表示foF23月分點(diǎn)值與9月分點(diǎn)值相等；AIfoF2的絕對(duì)值越大，則春秋分不對(duì)稱性越強(qiáng)．如圖5所示，當(dāng)F10.7分別為80、150、200時(shí)，AIfoF2分別為9.13%、8.99%、9.53%，ΔfoF2分別為0.69 MHz、1.04 MHz、1.26 MHz，即Akita站14時(shí)的foF2有3月分點(diǎn)值大于9月分點(diǎn)值的春秋分不對(duì)稱性．

圖3 AK539站地方時(shí)14∶00 foF2擬合的相對(duì)誤差分布Fig.3 Error distribution of fitted foF2 with respect to observed foF2 at AK539 at 14∶00 LT

圖4 (a)所有測(cè)高儀站foF2擬合相對(duì)誤差的期望值；(b) 所有測(cè)高儀站foF2擬合的均方根誤差分布Fig.4 (a) Averages of the relative error of fitted foF2 with respect to observed foF2 over all stations; (b) RMSE of fitted foF2 with respect to observed foF2 over all stations

圖5 AK539站14∶00不同太陽(yáng)活動(dòng)條件foF2春秋分不對(duì)稱性Fig.5 Equinoctial asymmetry of Akita foF2 at 14∶00 LT

3 結(jié)果

3.1 foF2春秋分不對(duì)稱性

圖6揭示了不同太陽(yáng)活動(dòng)水平，不同地方時(shí)foF2春秋分不對(duì)稱性指數(shù)AIfoF2及其春秋分差值ΔfoF2的分布情況．圖中三條藍(lán)綠色的線分別代表地磁傾角為20°、0°、-20°，紅色圓點(diǎn)代表該臺(tái)站處AIfoF2>0(ΔfoF2>0)，即foF23月分點(diǎn)值高于9月分點(diǎn)值；藍(lán)色圓點(diǎn)代表該臺(tái)站處AIfoF2<0(ΔfoF2<0)，即foF23月分點(diǎn)值低于9月分點(diǎn)值．從圖6中可以看出，日出時(shí)全球除少數(shù)臺(tái)站外foF2春秋分不對(duì)稱主要表現(xiàn)為9月分點(diǎn)值高于3月分點(diǎn)值，且隨著太陽(yáng)活動(dòng)增強(qiáng)，這種差異變得更弱．低太陽(yáng)活動(dòng)(F10.7=80)、中等太陽(yáng)活動(dòng)(F10.7=150)和高太陽(yáng)活動(dòng)(F10.7=200)，AIfoF2平均值分別為-5.02%，-3.53%，-2.54%，表明隨著太陽(yáng)活動(dòng)增強(qiáng)，日出時(shí)段foF2春秋分相對(duì)差異有所減弱；而ΔfoF2平均值分別為-0.25 MHz、-0.27 MHz、-0.24 MHz，表明日出時(shí)foF2春秋分絕對(duì)差受太陽(yáng)活動(dòng)影響較?。谖绾?、日落和午夜時(shí)，foF2的春秋分不對(duì)稱全球主要表現(xiàn)為3月分點(diǎn)值高于9月分點(diǎn)值，且隨著太陽(yáng)活動(dòng)增強(qiáng)，特別是從低太陽(yáng)活動(dòng)到中等太陽(yáng)活動(dòng)，foF23月分點(diǎn)值與9月分點(diǎn)值之相對(duì)差與絕對(duì)差均增大．午后和日落時(shí)的AIfoF2高于午夜的AIfoF2，特別是在低太陽(yáng)活動(dòng)條件下；同時(shí)，午后和日落的ΔfoF2在各太陽(yáng)活動(dòng)水平條件下均高于午夜的ΔfoF2．另外，圖6b中的午后時(shí)段，在120°E附近的ΔfoF2剖面表明在該經(jīng)度附近低太陽(yáng)活動(dòng)時(shí)foF2春秋分不對(duì)稱主要出現(xiàn)在低緯赤道異常峰地區(qū)，而高太陽(yáng)活動(dòng)時(shí)中緯的foF2春秋分不對(duì)稱也很顯著，這與Chen 等人的研究結(jié)果一致(Chen et al.,2012)．

3.2 foE春秋分不對(duì)稱性

圖7所示的是foE春秋分不對(duì)稱性指數(shù)AIfoE及其春秋分絕對(duì)差值ΔfoE的分布情況．由于電離層E層峰夜間基本消失，因此本文不考慮foE夜間的春秋分不對(duì)稱性．從圖7a中可以看出，AIfoE與AIfoF2情況類似．日出時(shí)，foE春秋分不對(duì)稱性主要表現(xiàn)為9月分點(diǎn)值高于3月分點(diǎn)值．在低太陽(yáng)活動(dòng)(F10.7=80)、中等太陽(yáng)活動(dòng)(F10.7=150)和高太陽(yáng)活動(dòng)(F10.7=200)情況下，AIfoE平均值分別為-4.84%、-5.69%、-5.74%，ΔfoE 平均值分別為-0.11 MHz、-0.14 MHz、-0.15 MHz，即日出時(shí)foE春秋分不對(duì)稱性隨著太陽(yáng)活動(dòng)增強(qiáng)而略有增大；E層春秋分絕對(duì)差值表現(xiàn)出的不對(duì)稱程度弱于F層，但其相對(duì)差值表現(xiàn)出的不對(duì)稱性卻強(qiáng)于F層．午后和日落時(shí)，foE春秋分不對(duì)稱性主要表現(xiàn)為3月分點(diǎn)值高于9月分點(diǎn)值,且日落時(shí)不對(duì)稱性明顯強(qiáng)于午后時(shí)不對(duì)稱性．圖7b展示了ΔfoE的分布情況，其結(jié)果與AIfoE分布類似．

圖6 foF2春秋分不對(duì)稱的分布．每個(gè)子圖左下角括號(hào)里為AIfoF2或ΔfoF2均值 (a) 春秋分不對(duì)稱性指數(shù)AIfoF2的分布;(b) 春秋分差值ΔfoF2的分布Fig.6 Scatter plots of (a) equinoctial asymmetry index and (b) equinoctial difference of foF2 versus geographic coordinates

圖7 foE春秋分不對(duì)稱的分布．每個(gè)子圖左下角括號(hào)里為AIfoE或ΔfoE均值． (a) 春秋分不對(duì)稱性指數(shù)AIfoE的分布；(b) 春秋分差值ΔfoE的分布.Fig.7 Scatter plots of (a) equinoctial asymmetry index and (b) equinoctial difference of foE versus geographic coordinates

圖8 foE春秋分不對(duì)稱性指數(shù)與foF2春秋分不對(duì)稱性指數(shù)的地方時(shí)變化圖中實(shí)線為各地方時(shí)所有臺(tái)站的春秋分不對(duì)稱指數(shù)的均值，誤差棒為其標(biāo)準(zhǔn)偏差.Fig.8 Local time variations of the equinoctial asymmetry indices of (left) foE and (right) foF2The solid lines are the averages of the equinoctial asymmetry indices of all ionosonde stations at each local time, and the error bars indicate the standard deviations of the indices.

圖9 TEC春秋分不對(duì)稱的分布(a) 春秋分不對(duì)稱性指數(shù)AITEC的分布；(b) 春秋分差值ΔTEC的分布.Fig.9 Longitudinal and latitudinal variations of (a) Equinoctial asymmetry index,(b) Equinoctial difference of TEC

圖8所示為AIfoE與AIfoF2隨地方時(shí)變化，其中實(shí)線分別代表各地方時(shí)AIfoE、AIfoF2的平均值，誤差棒為所有臺(tái)站AIfoE、AIfoF2在各地方時(shí)的標(biāo)準(zhǔn)差．從圖中可以看出，就平均值而言，AIfoE日出時(shí)為負(fù)值，在午后變得相對(duì)較小，到日落時(shí)變?yōu)橄鄬?duì)較大的正值．即早晨時(shí)段foE 9月分點(diǎn)值大于3月分點(diǎn)值，而下午時(shí)段則變?yōu)?月分點(diǎn)值大于9月分點(diǎn)值．太陽(yáng)活動(dòng)水平對(duì)AIfoE的地方時(shí)變化特征的影響并不顯著．而AIfoF2隨地方時(shí)變化的主要特點(diǎn)是，平均而言，日出時(shí)偏負(fù)值(尤其在低太陽(yáng)活動(dòng)時(shí))，而其余時(shí)段主要為正值．也就是說，日出時(shí)foF2主要表現(xiàn)為9月分點(diǎn)值高于3月分點(diǎn)值，而其余時(shí)段則相反．低太陽(yáng)活動(dòng)條件下foF2夜間的春秋分不對(duì)稱性弱．隨著太陽(yáng)活動(dòng)增強(qiáng)，日出時(shí)段foF2的春秋分不對(duì)稱性有減弱趨勢(shì)，而其他地方時(shí)段foF2的春秋分不對(duì)稱性有增強(qiáng)趨勢(shì)．

3.3 TEC春秋分不對(duì)稱性

圖9所示為TEC春秋分不對(duì)稱性分布．圖中黑色的線代表地磁傾角為0°．從圖9a中可以看出，日出時(shí)段，TEC在低太陽(yáng)活動(dòng)(F10.7=80)條件下全球主要表現(xiàn)為9月分點(diǎn)值高于3月分點(diǎn)值，這與foF2的春秋分不對(duì)稱性及foE的春秋分不對(duì)稱性類似，但在高太陽(yáng)活動(dòng)下，僅在低緯地區(qū)出現(xiàn)9月分點(diǎn)值高于3月分點(diǎn)值的春秋分不對(duì)稱性．而在其他時(shí)段，低太陽(yáng)活動(dòng)時(shí)，在午后和日落時(shí)段低緯地區(qū)主要表現(xiàn)為3月分點(diǎn)值高于9月分點(diǎn)值，在午夜時(shí)段傾角赤道附近有明顯的3月分點(diǎn)值高于9月分點(diǎn)值的不對(duì)稱性；中高太陽(yáng)活動(dòng)時(shí)，全球都主要表現(xiàn)為3月分點(diǎn)值高于9月分點(diǎn)值，且春秋分不對(duì)稱性程度較大(大于foE和foF2的春秋分不對(duì)稱性程度)．圖9b中ΔTEC展示了與AITEC類似的結(jié)果．Liu等(Liu et al.,2010)研究了GPS TEC在低太陽(yáng)活動(dòng)白天午后時(shí)春秋分不對(duì)稱性的全球分布特征，結(jié)果表明GPS TEC 3月分點(diǎn)值高于9月分點(diǎn)值的春秋分不對(duì)稱性主要出現(xiàn)在低緯地區(qū)，通過圖9中 F10.7=80，LT=1400的ΔTEC分布與其對(duì)比，結(jié)果符合得較好．

4 討論

熱層成分及電離層動(dòng)力學(xué)過程的春秋分差異可能是引起電離層春秋分不對(duì)稱的重要原因．Balan等(Balan et al.,1998)通過分析非相干散射雷達(dá)的數(shù)據(jù)表明，日本Shigaraki(34.85°N，136.10°E)地區(qū)300 km高度白天向極風(fēng)的速度春分日比秋分日低了20 m·s-1．因此，電離層在春分日相較秋分日有所抬升，春分日電離層有較低的復(fù)合損失率，秋分日電離層有較高的復(fù)合損失率．但由于缺少全球的中性風(fēng)觀測(cè)數(shù)據(jù)，而中性風(fēng)影響峰高的分布，因此本文嘗試分析氧氮濃度比及hmF2的春秋分不對(duì)稱性．foE春秋分不對(duì)稱性的原因需要進(jìn)一步研究．

對(duì)電離層F2層而言，電子產(chǎn)生率與氧原子濃度有關(guān)，而損失率與分子成分(氮分子與氧分子，二者對(duì)F2層復(fù)合損失的貢獻(xiàn)基本相當(dāng))的濃度有關(guān)．因此，[O]/[N2]是衡量F2層峰區(qū)光化平衡電子密度的重要指標(biāo)．由于原子成分與分子成分標(biāo)高的差異，靜壓平衡狀態(tài)下熱層[O]/[N2]隨高度增加而顯著減小．因而，在太陽(yáng)輻射與熱層狀態(tài)不變的情況下，F(xiàn)2層峰高h(yuǎn)mF2越大對(duì)應(yīng)的foF2就越高．有研究表明，hmF2與foF2有較好的正相關(guān)性(Danilov and Konstantinova,2013)．由于熱層中性風(fēng)等因素在兩個(gè)分點(diǎn)季節(jié)可能不一致(Balan et al.,1998)，全球F2層峰高可能也有春秋分不對(duì)稱性，從而影響到F2

層電子密度的春秋分不對(duì)稱性．因此，考察hmF2在兩分點(diǎn)的值．根據(jù)求取峰高的經(jīng)驗(yàn)公式(Dudeney,1983;Marin et al.,2001;McNamara,2008)：

hmF2=-176+1490

(15)

求得各臺(tái)站的hmF2數(shù)據(jù)，根據(jù)上述方法對(duì)hmF2的春秋分不對(duì)稱性進(jìn)行分析．如圖10所示，從統(tǒng)計(jì)上平均而言，hmF2除低太陽(yáng)活動(dòng)時(shí)的早晨與日落時(shí)段有3月分點(diǎn)值略低于9月分點(diǎn)值的情況之外，其他條件下AIhmF2平均值為正(AIhmF2代表hmF2春秋分不對(duì)稱性指數(shù))，即hmF23月分點(diǎn)值高于9月分點(diǎn)值．即hmF2的春秋分不對(duì)稱性與白天foF23月分點(diǎn)值高于9月分點(diǎn)值的不對(duì)稱性基本一致，但日出時(shí)段hmF2的春秋分不對(duì)稱性與foF2的春秋分不對(duì)稱性不一致．

熱層成分的季節(jié)變化是導(dǎo)致電離層一些季節(jié)變化特征的重要原因．研究表明，電離層季節(jié)變化中的半年異常和冬季異常就與氧氮濃度比[O]/[N2]的季節(jié)變化有關(guān)(Rishbeth et al.,2000;Yu et al.,2004)．在其他條件不變情況下，高的[O]/[N2]對(duì)應(yīng)F2層高的電子密度．熱層成分的季節(jié)變化也可能是引起電離層電子密度春秋分不對(duì)稱的原因．同時(shí)考慮F2層峰高的變化和熱層成分的變化，對(duì)hmF2處對(duì)應(yīng)的氧氮濃度比[O]/[N2]的春秋分不對(duì)稱性進(jìn)行分析．其中，氧原子濃度[O] 和氮?dú)夥肿訚舛萚N2]的數(shù)據(jù)由NRLMSISE-00模型(Picone et al.,2002)計(jì)算得到．如圖10所示，從統(tǒng)計(jì)上看，幾乎在低、中、高各太陽(yáng)活動(dòng)水平條件下從日出到日落的所有地方時(shí)段(低太陽(yáng)活動(dòng)時(shí)的日落除外)，AIrON2平均值均大于零(其中AIrON2代表氧氮濃度比[O]/[N2]春秋分不對(duì)稱性指數(shù)),即F2層峰高度的氧氮濃度比[O]/[N2]表現(xiàn)為3月分點(diǎn)值高于9月分點(diǎn)值．因此，可以部分解釋除日出時(shí)段外的foF2春秋分不對(duì)稱性(3月分點(diǎn)值高于9月分點(diǎn)值)．但對(duì)于日出時(shí)段foF2為什么會(huì)出現(xiàn)9月分點(diǎn)值高于3月分點(diǎn)值，還需要進(jìn)一步的研究．

電離層最主要的季節(jié)變化分量包括年變化與半年變化分量．從數(shù)學(xué)上而言，半年變化分量不論其相位如何，都不會(huì)引起電離層的春秋分不對(duì)稱性；而年變化分量，只要其幅度較大，其相位的改變會(huì)顯著影響到電離層的春秋分不對(duì)稱性．因此，本文從數(shù)學(xué)上考察TEC年變化的幅值與相位，以試圖理解TEC春秋分不對(duì)稱性的全球分布特點(diǎn)．圖11所展示的是低太陽(yáng)活動(dòng)地方時(shí)14∶00的TEC年變化幅值與相位分布圖．其中Aannual、φannual分別為年變化的幅值和相位．φannual代表了TEC年分量(即周期為12個(gè)月余弦變化)的峰值出現(xiàn)的時(shí)間．可以看出，TEC年變化的幅值和相位與圖9中的低太陽(yáng)活動(dòng)地方時(shí)14∶00的TEC春秋分不對(duì)稱性的全球分布特征符合得較好．例如，低太陽(yáng)活動(dòng)14∶00時(shí)低緯地區(qū)主要表現(xiàn)為3月分點(diǎn)值高于9月分點(diǎn)值，是因?yàn)槠淠曜兓至康南辔晃挥?月份，且其振幅在低緯較大．

圖10 hmF2及hmF2高度對(duì)應(yīng)的氧氮濃度比[O]/[N2]的春秋分不對(duì)稱性的地方時(shí)變化Fig.10 Local time variations of the equinoctial asymmetries of (left) hmF2 and (right) the ratio of [O]/[N2] at hmF2

圖11 低太陽(yáng)活動(dòng)地方時(shí)14∶00 TEC年變化的幅值(Aannual)與相位(φannual)圖Fig.11 Longitudinal and latitudinal variations of the amplitude and phase of TEC annual component under 14∶00 LT,solar minimum conditions

5 結(jié)論

本文通過建立電離層三個(gè)參量foF2、foE、TEC隨年積日與太陽(yáng)活動(dòng)的變化模型，研究了電離層的春秋分不對(duì)稱性．結(jié)果表明，電離層的春秋分不對(duì)稱特征有明顯的地方時(shí)依賴性．主要結(jié)果如下：

(1) 對(duì)于foF2，日出時(shí)全球除少數(shù)測(cè)高儀臺(tái)站外，主要表現(xiàn)為9月分點(diǎn)值高于3月分點(diǎn)值，其他時(shí)段則主要表現(xiàn)為3月分點(diǎn)值高于9月分點(diǎn)值，隨著太陽(yáng)活動(dòng)增強(qiáng)，foF2的春秋分不對(duì)稱性指數(shù)增大．

(2) 對(duì)于foE的春秋分不對(duì)稱性，日出和日落表現(xiàn)出截然相反的特點(diǎn)．日出時(shí)主要是9月分點(diǎn)值高于3月分點(diǎn)值，而下午及日落時(shí)段則變?yōu)?月分點(diǎn)值高于9月分點(diǎn)值，午后的春秋分不對(duì)稱性較弱．foE的春秋分不對(duì)稱性受太陽(yáng)活動(dòng)的影響較?。?/p>

(3) 對(duì)于TEC，日出時(shí)低緯地區(qū)主要為9月分點(diǎn)值高于3月分點(diǎn)值，而中高緯地區(qū)在低太活動(dòng)時(shí)表現(xiàn)為9月分點(diǎn)值高于3月分點(diǎn)值，中高太陽(yáng)活動(dòng)時(shí)3月分點(diǎn)值高于9月分點(diǎn)值；其他時(shí)段，除低太陽(yáng)活動(dòng)時(shí)3月分點(diǎn)值高于9月分點(diǎn)值主要出現(xiàn)在低緯地區(qū)外，中高太陽(yáng)活動(dòng)時(shí)全球主要為3月分點(diǎn)值高于9月分點(diǎn)值．

hmF2及氧氮濃度比的春秋分差異可以部分解釋F2層的春秋分不對(duì)稱的特征．此外，TEC年變化分量的相位和幅值分布可以較好地解釋TEC的春秋分不對(duì)稱的分布特征，但其機(jī)制需要更進(jìn)一步的研究．

致謝 測(cè)高儀觀測(cè)數(shù)據(jù)和TEC數(shù)據(jù)分別由SPIDR(http:∥spidr.ngdc.noaa.gov/)和ftp:∥cddis.gsfc.nasa.gov/pub/gps/products/ionex/網(wǎng)站提供；兩位審稿專家對(duì)本文提出建設(shè)性修改意見．在此一并表示感謝．

Aruliah A L,Farmer A D,Rees D,et al.1996.The seasonal behavior of high-latitude thermospheric winds and ion velocities observed over one solar cycle.J.Geophys.Res.,101(A7):15701-15711.

Bailey G J,Su Y Z,Oyama K I.2000.Yearly variations in the low-latitude topside ionosphere.Ann.Geophys.,18(7):789-798,doi:10.1007/s00585-000-0789-0.

Balan N,Otsuka Y,Fukao S.1997.New aspects in the annual variation of the ionosphere observed by the MU radar.Geophys.Res.Lett.,24(18):2287-2290.

Balan N,Otsuka Y,Bailey G J,et al.1998.Equinoctial asymmetries in the ionosphere and thermosphere observed by the MU radar.J.Geophys.Res.,103(A5):9481-9495.

Bardhan A,Aggarwal M,Sharma D K,et al.2014.Equinoctial asymmetry in low latitude ionosphere as observed by SROSS-C2 satellite.J.Atmos.Sol.-Terr.Phys.,117:101-109.

Bilitza D,Altadill D,Zhang Y L,et al.2014.The International Reference Ionosphere 2012—a model of international collaboration.J.Space Weather Space Clim.,4:A07.

Buonsanto M J.1986.Possible effects of the changing Earth-Sun distance on the upper atmosphere.South Pacific J.Nat.Sci.,8:58-65.

Chen Y D,Liu L B,Le H J.2008.Solar activity variations of nighttime ionospheric peak electron density.J.Geophys.Res.,113:A11306,doi:10.1029/2008JA013114.

Chen Y D,Liu L B,Wan W X,et al.2012.Equinoctial asymmetry in solar activity variations of NmF2and TEC.Ann.Geophys.,30:613-622,doi:10.5194/angeo-30-613-2012.

da Rosa A V,Waldman H,Bendito J,et al.1973.Response of the ionospheric electron content to fluctuations in solar activity.J.Atmos.Terr.Phys.,35(8):1429-1442.

Danilov A D,Konstantinova A V.2013.Relation between changes in foF2and hmF2within various time intervals.Geomagn.Aeron.,53(5):629-634,doi:10.1134/S0016793213050058.

Dudeney J R.1983.The accuracy of simple methods for determining the height of the maximum electron concentration of the F2-layer from scaled ionospheric characteristics.J.Atmos.Terr.Phys.,45(8-9):629-640.

Fejer B G,Scherliess L.1997.Empirical models of storm time equatorial zonal electric fields.J.Geophys.Res.,102(A11):24047-24056.

French A G.1966.Seasonal variation of some F-region parameters at sun-spot minimum.J.Atmos.Terr.Phys.,28(1):9-16.

Fukao S,Sato T,Kimura I,et al.1979.Seasonal mean structure of the night-time F2region over Arecibo.J.Atmos.Terr.Phys.,41(12):1205-1221,doi:10.1016/0021-9169(79)90024-2.

Fukao S,Oliver W L,Onishi Y,et al.1991.F-region seasonal behavior as measured by the MU radar.J.Atmos.Terr.Phys.,53(6-7):599-618.

Ganesh K S.1965.Solar cycle and seasonal variations in the equatorial ionospheric F-region.J.Atmos.Terr.Phys.,27(5):645-655,doi:10.1016/0021-9169(65)90132-7.

Jakowski N,Bettac H D,Lazo B,et al.1981.Seasonal variations of the columnar electron content of the ionosphere observed in Havana from July 1974 to April 1975.J.Atmos.Terr.Phys.,43(1):7-11.

Kawamura S,Balan N,Otsuka Y,et al.2002.Annual and semiannual variations of the midlatitude ionosphere under low solar activity.J.Geophys.Res.,107:SIA 8-1-SIA 8-10,doi:10.1029/2001JA000267.

Liu H X,Lühr H,Watanabe S.2007.Climatology of the equatorial thermospheric mass density anomaly.J.Geophys.Res.,112:A05305,doi:10.1029/2006JA012199.

Liu L B,Wan W X,Ning B Q.2004.Statistical modeling of ionospheric foF2over Wuhan.Radio Sci.,39:RS2013,doi:10.1029/2003RS003005.

Liu L B,Wan W X,Ning B Q,et al.2006.Solar activity variations of the ionospheric peak electron density.J.Geophys.Res.,111:A08304,doi:10.1029/2006JA011598.

Liu L B,Zhao B Q,Wan W X,et al.2009.Seasonal variations of the ionospheric electron densities retrieved from Constellation Observing System for Meteorology,Ionosphere,and Climate mission radio occultation measurements.J.Geophys.Res.,114:A02302,doi:10.1029/2008JA013819.

Liu L B,He M S,Yue X A,et al.2010.Ionosphere around equinoxes during low solar activity.J.Geophys.Res.,115:A09307,doi:10.1029/2010JA015318.

Ma R P,Xu J Y,Liao H Z.2003.The features and a possible mechanism of semiannual variation in the peak electron density of the low latitude F2 layer.J.Atmos.Sol.Terr.Phys.,65(1):47-57.

Marin D,Mikhailov A V,de la Morena B A,et al.2001.Long-term hmF2trends in the Eurasian longitudinal sector from the ground-based ionosonde observations.Ann.Geophys.,19(7):761-772.

McNamara L F.2008.Accuracy of models of hmF2used for long-term trend analyses.Radio Sci.,43:RS2002,doi:10.1029/2007RS003740.

Millward G H,Rishbeth H,Fuller-Rowell T J,et al.1996.Ionospheric F2layer seasonal and semiannual variations.J.Geophys.Res.,101(A3):5149-5156.

Natali M P,Meza A.2010.Annual and semiannual VTEC effects at low solar activity based on GPS observations at different geomagnetic latitudes.J.Geophys.Res.,115:D18106,doi:10.1029/2010JD014267.

Pavlov A V,Pavlova N M.2013.Variations in statistical parameters of the NmF2equinoctial asymmetry with latitude and solar activity near noon.Adv.Space Res.,51(11):2018-2034.

Pezzopane M,Scotto C.2007.Automatic scaling of critical frequency foF2and MUF(3000)F2:A comparison between Autoscala and ARTIST 4.5 on Rome data.Radio Sci.,42:RS4003,doi:10.1029/2006RS003581.

Picone J M,Hedin A E,Drob D P,et al.2002.NRLMSISE-00 empirical model of the atmosphere:Statistical comparisons and scientific issues.J.Geophys.Res.,107(A12):SIA 15-1-SIA 15-16,doi:10.1029/2002JA009430.

Rao P V S R,Ram P S,Jayachandran P T,et al.1996.Seasonal variation in ionospheric electron content and irregularities over Waltair—A comparison with SLIM model.Adv.Space Res.,18(6):259-262.

Richards P G.2001.Seasonal and solar cycle variations of the ionospheric peak electron density:Comparison of measurement and models.J.Geophys.Res.,106(A7):12803-12819.

Rishbeth H,Müller-Wodarg I C F,Zou L,et al.2000.Annual and semiannual variations in the ionospheric F2-layer:II.Physical discussion.Ann.Geophys.,18(8):945-956,doi:10.1007/s00585-000-0945-6.

Roble R G,Salah J E,Emery B A.1977.The seasonal variation of the diurnal thermospheric winds over Millstone Hill during solar cycle maximum.J.Atmos.Terr.Phys.,39(4):503-511.

Schmelovsky K H.1963.Diurnal and seasonal variations of electron density distribution in the upper F-region.J.Atmos.Terr.Phys.,25(5):241-247,doi:10.1016/0021-9169(63)90020-5.

Titheridge J E.1973.The electron content of the southern mid-latitude ionosphere,1965-1971.J.Atmos.Terr.Phys.,35(5):981-1001,doi:10.1016/0021-9169(73)90077-9.

Tsai H F,Liu Y J,Tsai W H,et al.2001.Seasonal variations of the ionospheric total electron content in Asian equatorial anomaly regions.J.Geophys.Res.,106(A12):30363-30369.

Unnikrishnan K,Nair R B,Venugopal C.2002.Harmonic analysis and an empirical model for TEC over Palehua.J.Atmos.Sol.-Terr.Phys.,64(17):1833-1840,doi:10.1016/S1364-6826(02)00187-6.

Yang J,Liu L B,Chen Y D,et al.2012.Does the equatorial ionospheric peak electron density really record the lowest during the recent deep solar minimum?.Chinese J.Geophys.,55(5):457-465.

Yu T,Wan W X,Liu L B,et al.2004.Global scale annual and semi-annual variations of daytime NmF2in the high solar activity years.J.Atmos.Sol.Terr.Phys.,66(18):1691-1701,doi:10.1016/j.jastp.2003.09.018.

Zhang S R,Holt J M.2007.Ionospheric climatology and variability from long-term and multiple incoherent scatter radar observations:Climatology in eastern American sector.J.Geophys.Res.,112:A06328,doi:10.1029/2006JA012206.

(本文編輯 胡素芳)

Local time dependence of ionospheric equinoctial asymmetry

LIU Yong1,2,3,CHEN Yi-Ding1,2,3*,LIU Li-Bo1,2

1 Key Laboratory of Earth and Planetary Physics,Institute of Geology and Geophysics, Chinese Academy of Sciences,Beijing 100029,China2 Beijing National Observatory of Space Environment,Institute of Geology and Geophysics, Chinese Academy of Sciences,Beijing 100029,China3 University of the Chinese Academy of Sciences,Beijing 100049,China

In this paper,the critical frequency of the F2-layer (foF2) and the critical frequency over the E-layer (foE) of global 203 ionosonde stations and Jet Propulsion Laboratory (JPL) total electron content (TEC) map data were used to investigate the characteristics of ionospheric equinoctial (March Equinox and September Equinox，ME and SE) asymmetry.By fitting the three parameters foF2,foE and TEC in terms of a function of Day of Year (DoY) and the solar activity index F10.7,the three parameters were normalized to the low solar activity level of F10.7=80,the moderate solar activity level of F10.7=150,and the high solar activity level of F10.7=200.Thus,by this method,the actually observed seasonal and solar activity variations of ionospheric parameters,which are interlaced with each other,can be quantitatively separated to acquire more accurate characteristics of ionospheric equinoctial asymmetry.The equinoctial asymmetry index (AI) and the equinoctial difference (Δ=ME-SE) of the three ionospheric parameters at different local times were analyzed to get the global features and the solar activity dependence of equinoctial asymmetry.In general,foE is mainly higher at SE than at ME over the global ionosonde stations at sunrise,and the equinoctial asymmetry of foE almost disappears in the afternoon,while foE is higher at ME than at SE at sunset.foF2is higher at SE than at ME over the most of the global ionosonde stations at sunrise,while it is higher at ME than at SE at other local times.With respect to TEC,the equinoctial asymmetry of higher TEC at SE than at ME exists over the global at low solar activity level but mainly occurs at low latitudes at moderate and high solar activity levels at sunrise,while TEC is mainly higher at ME than at SE at other local times.The dependence of the equinoctial asymmetry of foE is weak,while for foF2and TEC,solar activity dependence of the equinoctial asymmetry is evident.The equinoctial asymmetries of the oxygen-nitrogen ratio ([O]/[N2],computed from the atmospheric model NRLMSISE-00) at the peak height of the F2-layer (hmF2) and hmF2were investigated;they can be used to partially explain the equinoctial asymmetry of foF2.The amplitude and the phase of TEC annual variation were calculated.The phase of TEC annual variation can be used to account for the global distribution of TEC equinoctial asymmetry.

foF2;foE;GPS TEC;Ionosphere;Equinoctial asymmetry

劉勇，陳一定，劉立波.2016.電離層春秋分不對(duì)稱的地方時(shí)依賴.地球物理學(xué)報(bào)，59(11):3941-3954,

10.6038/cjg20161101.

Liu Y,Chen Y D,Liu L B.2016.Local time dependence of ionospheric equinoctial asymmetry.Chinese J.Geophys.(in Chinese),59(11):3941-3954,doi:10.6038/cjg20161101.

國(guó)家自然科學(xué)基金(41274161，41231065，41321003)和973項(xiàng)目課題(2012CB825604)共同資助．

劉勇，男，1989年生，漢族，四川內(nèi)江人，碩士生在讀，主要研究方向?yàn)殡婋x層物理．E-mail:455509433@qq.com

*通迅作者 陳一定，男，1981年生，漢族，安徽懷寧人，博士，副研究員，主要從事電離層物理研究．E-mail:chenyd@mail.iggcas.ac.cn

10.6038/cjg20161101

P352

2015-12-17，2016-04-24收修定稿