我国极区冬季电离层加热实验研究

通过对非相干散射雷达观测数据的处理分析,研究了2008年1月我国在挪威Troms进行的冬季电离层加热实验效应.研究结果表明,电离层临界频率大于泵波频率的O波加热事件扰动效应明显,电子温度存在60%~120%的增强,扰动范围从150 km一直延伸到400 km,电子密度扰动不显著,最大可以观察到12%的密度衰减.受加热影响,离子声波频率有1~2 kHz的增加,离子线谱峰谷比增加,有时伴随有高阶谐振线出现.离子线和等离子线功率存在过冲现象,等离子线的功率剖面存在单峰、双峰和三峰结构,等离子线的功率增强幅度随频率负指数衰减.     

曲靖非相干散射雷达电离层F区日间电子温度变化特征初步分析

利用曲靖非相干散射雷达2017-2018年春夏季观测数据首次分析了电离层日间150~450km电子温度的地方时与高度变化特征及其与电子密度的相关性.发现hmF2及以上的电子温度在日出日落时具有两个峰值,在11∶00-16∶00LT之间变化较小,高度越高午后上升的时间越早;从150km开始迅速增加,在约220km达到最大值,然后开始降低,在约300~350km达到最小值,最后单调上升;200km以下电子温度与电子密度成正相关(主要由热传导控制),200~450km之间存在明显的反相关(光电离过程占主导),电子-离子温度差与电子密度对数之间存在近似线性关系,电子温度逐日变化与光电离因子的变化趋势相似,这种相关性在中午与午后更明显;以上结果与其他非相干散射雷达观测和电离层模型计算结果基本一致,但也存在一些差别,需要结合更多数据深入分析.     

应用经验正交函数估算顶部电离层电子密度剖面

本文基于IRI模型、地面数字测高仪和GNSS TEC数据,提出了一种利用经验正交函数(Empirical Orthogonal Function,简称EOF)估算顶部电离层电子密度剖面的方法,并将其应用于美国Millstone Hill测高仪和GNSS数据以估算顶部电离层电子密度剖面.通过将估算的临界频率、峰值高度、400km以上电子密度分别与测高仪实测临界频率、测高仪实测峰值高度以及非相干散射雷达实测400km以上电子密度作对比以对方法的有效性进行验证.统计结果显示估算临界频率、峰值高度与测高仪实测数据基本一致,400km以上估算电子密度相较于非相干散射雷达实测的绝对误差平均值仅是测高仪推算400km以上电子密度绝对误差平均值的一半左右.所以本文提出的方法可以更加精确地估算顶部电离层电子密度.     

利用GRACE卫星精密微波测距确定星间平均电子密度

本文介绍如何利用GRACE两颗卫星之间K波段双频微波精密测距和轨道数据,得到星间平均电子密度.发展了一种将连续轨道电子密度极小对齐到零的方法,以消除整周模糊度;借助CHAMP卫星朗缪探针测量得到的轨道电子密度基值以及GPS掩星数据计算的等离子体垂直梯度标高,进一步修正了GRACE星间电子密度所固有的偏差;从而得到大约500 km高度上长达近十年的全球电子密度数据.为了检验消除偏差后GRACE星间电子密度数据的可靠性,对比了GRACE卫星过Millstone Hill雷达上空时,非相干散射雷达观测到的大致同时和相近位置的电子密度数据,结果显示,二者之间的线性相关系数为0.97,平均偏差为-7.26%,GRACE星间电子密度总体稍微偏低,偏差的标准差为18.6%.为进一步验证本文方法所得数据的可用价值,利用消除偏差后的电子密度数据,对GRACE卫星与CHAMP卫星在近乎相同的地方时而高度不同的近圆极轨道上飞行的情况下,两颗卫星观测到的电子密度随经度和纬度的全球分布进行了对比分析.多方面的对比检验证明,本文方法得到的几乎连续10年的GRACE高度上全球电子密度数据基本可靠,为电离层气候学与天气学研究提供了宝贵资料.     

顶部电离层和等离子体层电子密度分布——基于GRACE星载GPS信标测量的CT反演

本文利用两颗跟飞的GRACE卫星载GPS信标测量数据和基于差分相对TEC的层析算法,实现了全球范围的顶部电离层和等离子体层(450~5000 km) 层析成像.反演结果表明,利用低轨道卫星载GPS信标测量数据可以有效地重建顶部电离层和等离子体层的全球二维分布图像.对不同地磁活动条件下的天基层析反演结果表明,等离子体层电子密度随纬度的分布是不均匀的;在低纬赤道带,从顶部电离层向上延伸直到等离子体层,以及等离子体层中局地的电离增强云团,经常出现近似垂直于磁力线的电子密度柱状增强结构.     

极区电离层加热的数值模拟与实验对比

大功率无线电波能加热电离层等离子体,可以引起电离层电子温度和密度的扰动,实现电离层的人工变态.从电子的连续性方程、动量方程和能量方程出发,我们给出了地面人工大功率无线电波加热电离层的数值模型.通过对方程的数值求解,计算了极地电离层条件下,电子温度、电子密度的加热效应,讨论了泵波参数对加热效应的影响.研究结果表明,电子温度几乎在整个高度上表现为一致性的幅度增强,且在反射高度附近形成温度增强峰很平缓.电子密度在峰上高度附近形成密度谷,谷两侧存在密度增强.加热效应基本随加热功率的增大而增大,随加热频率的增大而减小.使用我国2008年1月在挪威进行的电离层加热实验的电离层参量作为仿真初值,对6个O波加热时刻进行了数值仿真,仿真结果与实验观测基本保持一致.     

磁暴期间极光椭圆与极盖区电离层效应比较研究——F区负暴

联合利用EISCAT和E-Svalbard非相干散射雷达数据,研究l997年5月强磁暴期间向阳侧极盖与极光椭圆区电离层F区负暴.发现在磁暴主相和恢复相初期,极光椭圆和极盖区电离层都在大约l90km高度出现类似F1的峰,F2主峰完全消失,F区电子密度大幅度下降.但离子温度的变化在两个区域很不相同,在极光椭圆区大幅度升高,而在极盖区没有显著变化,反映出引起F区负暴的主要机制在两个区域不尽相同.强对流电场引起大气焦耳加热与离子增温而使O+离子消失的化学反应速率增大所导致的电离损失,对极光椭圆区负暴起主要作用;而输运过程,特别是持续长达数小时的沿场上行离子流,对极盖区负暴起重要作用.磁暴主相期间,当EISCAT雷达位于等离子体对流涡旋转换区下方时,在无焦耳加热与离子摩擦增温的情况下,观测到由顶部电离层O+离子沿场高速外流引起的F区电子密度耗空.     

基于非相干散射雷达和GPS观测研究Millstone Hill地区等离子体层电子含量

本文尝试结合非相干散射雷达和GPS TEC观测数据提取等离子体层总电子含量(PTEC).我们首先描述所用的技术方法,然后具体利用了Millstone Hill台站的观测数据研究该地区上空等离子体层总电子含量(PTEC)的变化情况.我们采用变化标高的Chapman函数对非相干散射雷达测得的电子浓度剖面数据进行拟合,然后通过对剖面积分得到100 km到1000 km高度范围的电离层总电子含量.GPS提供的TEC数据为高度达20200 km的总电子含量,两者之差可近似看成等离子体层的电子含量.本文分别选取太阳活动高年(2000, 2002年)和太阳活动低年(2005,2008年)Millstone Hill台站的静日数据进行研究.结果表明,等离子体层电子含量及其所占GPS TEC的比例具有明显的周日变化.PTEC含量在白天高于夜间,而所占GPS TEC的百分比,夜间明显高于白天.太阳活动高年所选月份等离子体层电子含量在4~14 TECU (1TECU=10¹⁶el/m²) 范围内变化,夜间所占比例可达60%左右.太阳活动低年所选月份等离子体层电子含量在3~7 TECU范围内变化,所占比例夜间最高可达80%左右.我们所得到的结果与前人基于其它观测手段所得结果在变化趋势上一致,在量级上也大致相当.因此,这从一个侧面证明了我们所用方法的可靠性.非相干散射雷达能够探测包括F2层峰值以下及以上高度的电子浓度,利用这一设备所观测得到的资料来推算电离层电子含量将比前人基于电离层垂测仪观测资料进行的推算更具真实性,由此得到的等离子体层电子含量也将更为接近真实情况.     

高纬电离层气候学特征研究--EISCAT雷达观测及与IRI模式的比较

利用1988～1999年欧洲非相干散射EISCAT（European Incoherent Scatter）雷达观测数据,对不同太阳活动周相、不同季节的极光椭圆区电离层F区电子密度进行统计分析,研究其气候学特征,并与IRI 2001模式比较.EISCAT观测到的电子密度显示出显著的太阳活动高年“冬季异常”和太阳活动低年半年变化等现象.EISCAT实测电子密度随时间和高度的平均二维分布和500 km高度以下总电子含量TEC,从总体来看与IRI 2001模式预测结果符合较好.但高年在TEC达到最大值前后,IRI 2001模式预测的电子密度高度剖面与EISCAT观测结果有显著差别：F2峰以上IRI 2001模式预测的电子密度过大,造成TEC明显高于雷达观测值.另外,在太阳活动下降相,EISCAT观测显示出明显的半年周期季节变化特征,但IRI 2001模式未能预测出此下降相季节变化.     

电离层诸参量对大气重力波的响应与重力波参数反演

本文利用欧洲非相干散射雷达数据,分析研究了电离层不同等离子体参量对大气重力波的响应之间的关系.应用这种关系,发展了一种在垂直于地磁场的电场可以忽略等简化假设下,由电子密度和离子沿场速度的同时测量数据,反演求解较高（约250km以上）F区中引起TID的重力波传播参数的方法.用此法对一典型TID事件进行分析计算,所得结果与全波解数值研究结果很好符合.     

Electron velocity distribution function in a plasma with temperature gradient and in the presence of suprathermal electrons: application to incoherent-scatter plasma lines

The plasma dispersion function and the reduced velocity distribution function are calculated numerically for any arbitrary velocity distribution function with cylindrical symmetry along the magnetic field. The electron velocity distribution is separated into two distributions representing the distribution of the ambient electrons and the suprathermal electrons. The velocity distribution function of the ambient electrons is modelled by a near-Maxwellian distribution function in presence of a temperature gradient and a potential electric field. The velocity distribution function of the suprathermal electrons is derived from a numerical model of the angular energy flux spectrum obtained by solving the transport equation of electrons. The numerical method used to calculate the plasma dispersion function and the reduced velocity distribution is described. The numerical code is used with simulated data to evaluate the Doppler frequency asymmetry between the up- and downshifted plasma lines of the incoherent-scatter plasma lines at different wave vectors. It is shown that the observed Doppler asymmetry is more dependent on deviation from the Maxwellian through the thermal part for high-frequency radars, while for low-frequency radars the Doppler asymmetry depends more on the presence of a suprathermal population. It is also seen that the full evaluation of the plasma dispersion function gives larger Doppler asymmetry than the heat flow approximation for Langmuir waves with phase velocity about three to six times the mean thermal velocity. For such waves the moment expansion of the dispersion function is not fully valid and the full calculation of the dispersion function is needed.     

等离子体湍流对电子的加速

本文讨论了等离子体湍流对电子加速的两种模型:（1）假定在空间中存在一个空间均匀的等离子体湍流区,当具有一定初始分布的电子束通过此湍流区时,研究湍流场对电子束的加速过程;（2）在某一封闭的区域中,存在着具有一定初始分布和空间均匀的等离子体,当某种类型的等离子体波突然传入此等离子体区,然后考察此区中电子的加速过程。在这两种模型中,可能存在着某种电子消失机制。假定湍谱是幂指数形式,我们给出了不同类型湍流扩散系数的普遍形式。利用较简单的数学方法,求解了包括消失过程的一维准线性动力学方程,对于给定的初始分布,得出了分布函数的解析解,并给出了平均能量时间关系的表达式。另外,对于特定的湍谱指数,解出了当平行电场和湍流同时存在时的分布函数。最后,对所得结果进行了数值分析和讨论。     

考虑电子吸附效应的低电离层加热研究

基于低电离层自洽加热模型，综合考虑了低电离层中电子的复合效应及典型吸附效应，本文数值仿真了大功率高频无线电波持续加热低电离层所产生的电子温度、电子密度的扰动，并且首次模拟分析了由于电子温度扰动造成的加热电波自吸收效应.结果表明：电子吸收大功率加热电波能量导致了电子温度的增加，同时改变了电离层的吸收指数，引起了加热电波的自吸收效应.加热电波的自吸收效应对低电离层较高区域的电子温度扰动有重要的抑制作用.因此，随着加热频率的减小或有效辐射功率的增大，低电离层较低区域的电子温度增量明显增大而在高度100 km以上区域的电子温度增量始终较小.另一方面，随着电子温度的增加，电子的复合系数减小而电子的吸附系数增加，导致了电子密度在低电离层中较高区域出现正扰动而在较低区域出现负扰动.当饱和电子温度较大时，继续减小加热频率或增大有效辐射功率对电子密度扰动所造成的改变较小，尤其当电子温度超出复合系数和吸附系数的温度敏感区间.此外，电子温度与电子密度的饱和时间相差较大，电子温度的饱和时间为微秒量级而电子密度的饱和时间为秒量级.     

Coordinated Akebono and EISCAT observations of suprathermal ion outflows in the nightside inverted-V region

We have investigated ion outflows observed by the Akebono satellite and the EISCAT radar in the nightside auroral region on February 16, 1993. The Akebono satellite at about 7000 km altitude observed the region of suprathermal ion outflows and inverted-V type electron precipitation alternately with a horizontal separation of 70–150 km at the ionospheric level. These two regions corresponded to the upward and downward field-aligned current region, respectively, and intense ELF waves were observed in the ion outflow region. From the EISCAT VHF radar observation (Common Program 7 mode), it has been suggested that the ion outflow region and the enhanced electron temperature region were aligned along geomagnetic field lines with vertical and horizontal separations of 200–400 and 70–80 km, respectively and these two regions convected equatorward across the EISCAT radar at Tromsø site. Based on these results, we propose a model for this ion outflow as follows. In the nightside auroral region, downward FAC regions exist near the edge of the inverted-V type electron precipitation regions. ELF waves are excited probably by a plasma instability due to the upward thermal electron beam carrying the downward FACs, and these ELF waves cause transverse ion heating at the top of the ionosphere. The produced ion conics contribute significantly to ion outflow.     

Observations of flux rope – associated particle bursts with GEOTAIL in the distant tail

Geotail energetic particle, magnetic field data and plasma observations (EPIC, MGF and CPI experiments) have been examined for a number of energetic particle bursts in the distant tail (120Re < |X_GSM| < 130 Re), associated with moving magnetic field structures, following substorm onsets. The features obtained from this data analysis are consistent with the distant magnetotail dynamics determined first by ISEE3 observations and explained in terms of the neutral line model. At the onset of the bursts, before plasma sheet entrance, energetic electrons appear as a field-aligned beam flowing in the tailward direction, followed by anisotro-pic ions. Within the flux rope region, suprathermal ions exhibit a convective anisotropy, which allows determination of the plasma flow velocity, assuming that the anisotropy arises from the Compton-Getting effect. The velocities thus determined in the plasma sheet are estimated to be 200–650 km/s, and compare favourably with the velocities derived from the CPI electron and proton experiment. The estimated length of magnetic field structures varies between 28 and 56 Re and depends on the strength of the westward electrojet intensification. Finally, the three structures reported here show clear magnetic field signatures of flux rope topology. The existence of a strong magnetic field aligned approximately along the Y-axis and centred on the north-to-south excursion of the field, and the bipolar signature in both By and/or Bz components, is consistent with the existence of closed field lines extending from Earth and wrapping around the core of the flux rope structure.     

Transport of thermal plasma above the auroral ionosphere in the presence of electrostatic ion-cyclotron turbulence

The electron component of intensive electric currents flowing along the geomagnetic field lines excites turbulence in the thermal magnetospheric plasma. The protons are then scattered by the excited electromagnetic waves, and as a result the plasma is stable. As the electron and ion temperatures of the background plasma are approximately equal each other, here electrostatic ion-cyclotron (EIC) turbulence is considered. In the nonisothermal plasma the ion-acoustic turbulence may occur additionally. The anomalous resistivity of the plasma causes large-scale differences of the electrostatic potential along the magnetic field lines. The presence of these differences provides heating and acceleration of the thermal and energetic auroral plasma. The investigation of the energy and momentum balance of the plasma and waves in the turbulent region is performed numerically, taking the magnetospheric convection and thermal conductivity of the plasma into account. As shown for the quasi-steady state, EIC turbulence may provide differences of the electric potential of δ V ≈ 1–10 kV at altitudes of 500 < h < 10 000 km above the Earth’s surface. In the turbulent region, the temperatures of the electrons and protons increase only a few times in comparison with the background values.     

Enhanced incoherent scatter plasma lines

Detailed model calculations of auroral secondary and photoelectron distributions for varying conditions have been used to calculate the theoretical enhancement of incoherent scatter plasma lines. These calculations are compared with EISCAT UHF radar measurements of enhanced plasma lines from both the E and F regions, and published EISCAT VHP radar measurements. The agreement between the calculated and observed plasma line enhancements is good. The enhancement from the superthermal distribution can explain even the very strong enhancements observed in the auroral E region during aurora, as previously shown by Kirk-wood et al. The model calculations are used to predict the range of conditions when enhanced plasma lines will be seen with the existing high-latitude incoherent scatter radars, including the new EISCAT Svalbard radar. It is found that the detailed structure, i.e. the gradients in the suprathermal distribution, are most important for the plasma line enhancement. The level of superthermal flux affects the enhancement only in the region of low phase energy where the number of thermal electrons is comparable to the number of suprathermal electrons and in the region of high phase energy where the suprathermal fluxes fall to such low levels that their effect becomes small compared to the collision term. To facilitate the use of the predictions for the different radars, the expected signal-to-noise ratios (SNRs) for typical plasma line enhancements have been calculated. It is found that the high-frequency radars (Søndre Strømfjord, EISCAT UHF) should observe the highest SNR, but only for rather high plasma frequencies. The VHP radars (EISCAT VHP and Svalbard) will detect enhanced plasma lines over a wider range of frequencies, but with lower SNR.     

Effect of suprathermal electrons on the intensity and Doppler frequency of electron plasma lines

In an incoherent scattering radar experiment, the spectral measurement of the so-called up- and downshifted electron plasma lines provides information about their intensity and their Doppler frequency. These two spectral lines correspond, in the backscatter geometry, to two Langmuir waves travelling towards and away from the radar. In the daytime ionosphere, the presence of a small percentage of photoelectrons produced by the solar EUV of the total electron population can excite or damp these Langmuir waves above the thermal equilibrium, resulting in an enhancement of the intensity of the lines above the thermal level. The presence of photo-electrons also modifies the dielectric response function of the plasma from the Maxwellian and thus influences the Doppler frequency of the plasma lines. In this paper, we present a high time-resolution plasma-line data set collected on the EISCAT VHF radar. The analysed data are compared with a model that includes the effect of a suprathermal electron population calculated by a transport code. By comparing the intensity of the analysed plasma lines data to our model, we show that two sharp peaks in the electron suprathermal distribution in the energy range 20–30 eV causes an increased Landau damping around 24.25 eV and 26.25 eV. We have identified these two sharp peaks as the effect of the photoionisation of N2 and O by the intense flux of monochromatic HeII radiation of wavelength 30.378 nm (40.812 eV) created in the chromospheric network and coronal holes. Furthermore, we see that what would have been interpreted as a mean Doppler drift veloCity for a Maxwellian plasma is actually a shift of the Doppler frequency of the plasma lines due to suprathermal electrons.     

行星际扰动和地磁活动对GEO相对论电子影响

利用1988—2010年小时平均的GOES卫星数据,对地球同步轨道(GEO)相对论电子变化进行了统计分析,研究了相对论电子通量(Fe)增强事件的发展过程,探讨了利于相对论电子通量增强的太阳风和地磁活动条件.主要结论如下:(1)GEO相对论电子通量即使是峰值,也具有明显的地方时特性,最大电子通量出现在磁正午时.午/夜电子通量比率随着太阳风速度(Vsw)增加而增大;在Dst-50nT时相对论电子具有规则的地方时变化.在太阳活动下降相,电子通量与各参数的相关性较好,与其相关性最好的Vsw、Kp指数以及三次根号下的太阳风密度(N)分别出现在电子通量前39~57h、57~80h和12~24h.(2)强(日平均电子通量峰值Femax≥104 pfu)相对论电子事件,在距离太阳活动谷年前两年左右和春秋分期间发生率最高,较弱(104Femax≥103 pfu)事件无此特点;大部分强相对论电子事件中,电子通量在磁暴主相开始增加,较弱事件中则在恢复相开始回升.(3)太阳风密度变化对相对论电子事件的发展具有重要指示作用.电子通量在太阳风密度极大值后0~1天达到极小值,太阳风密度极小值后0~2天达到极大值.(4)90%以上相对论电子事件是在磁暴及高速太阳风的条件下发生的,与其伴随的行星际参数和地磁活动指数极值满足以下条件:Vswmax516km/s,Dstmin-31nT,Nmin2.8cm-3,Nmax14.1cm-3,Bzmin-2.9nT,AEmax698nT.(5)磁暴过程中,Dstmin后日平均电子通量大于103 pfu的发生概率为53%左右,强/弱相对论电子事件占总数比例分别为36%/64%左右,磁暴强度对其无影响.磁暴过程中的Vsw、N和AE指数大小对于能否引起相对论电子增强起着指示作用.