短波射线追踪技术中的电离层混合建模方法

在国际参考电离层模型和多层准抛物模型的基础上,提出了一种混合应用两种模型进行电离层建模的新方法.利用射线追踪技术,分别对混合模型和传统国际参考电离层模型下短波射线在电离层中的轨迹进行了仿真,得到了电波群路径.通过与实测电波群路径的对比,结果表明:对中国中纬度地区在电离层混合模型下的射线追踪精度优于传统国际参考电离层模型下的射线追踪精度,同时混合建模方法降低了多层准抛物模型对输入条件的要求,扩展了多层准抛物模型在射线追踪技术中的应用范围.     

电离层Es电波传播模型研究

考虑大气折射的影响,对不同高度Es层条件下大圆路径传播距离与仰角的关系进行计算.基于电离层Es电子密度时空分布特点,建立高阶Es反射模型;运用高阶修正贝塞尔函数表示电离层Es电子密度扰动的自相关函数,综合考虑电离层Es各向异性结构的尺度关系,建立高阶Es散射模型.仿真计算Es层VHF斜入射传播衰减与频率、距离的关系.对于反射模型,Es层厚度越厚,电波衰减越大.对于散射模型,Es层不规则体水平尺度相同,垂直水平漂移的尺度越大,衰减越大;而垂直水平漂移的尺度相同,水平尺度越大,衰减越小.无论反射模型还是散射模型,f₀E_s越高,衰减越小;电波工作频率越高,衰减越大.对比建立的模型与VHF链路的测量结果,明确不同强度的电离层Es反射/散射机制,证明了建立模型的正确性.      

Galileo卫星导航系统广播Nequick模型及其参数拟合

提出了基于IGRF模型的Galileo广播Nequick模型及其参数拟合算法, 解决了Galileo信号仿真中地理场景映射与地磁坐标下的电离层延时修正参数拟合问题. 应用IGRF模型, 可计算出任意给定位置和时间点的地磁参数以及E层、 F₁层、F₂层的电子密度, 从而计算出Galileo电离层修正参数. 仿真结果表明, 该算法拟合的全球电离层延时与IGS提供的实际观测值基本一致, 仿真精度高于一般的经验电离层模型, 实现了Galileo卫星信号的电离层延时修正参数的精确仿真.      

基于一维模型的高频电波加热电离层F层数值研究

基于电离层一维仿真加热模型,详细介绍了模型中电子的动量方程、连续性方程、能量方程和各类参量表达式,利用对角矩阵追赶法数值求解电离层F层加热过程,分析了不同时次电子数密度和电子温度的变化,讨论了不同频率和不同功率电波加热的情形.结果表明:当高频电波加热高电离层时,电子温度迅速增加,并很快趋近于饱和状态;电子密度的变化较为迟缓,但在加热过程中其变化幅度却越来越大;电子密度变化量在沿磁场方向上形成空洞和上下稠团两峰一谷构型;频率越大、功率越高的电波加热时,电子密度的变化也越大,但存在一适值频率的电波对电子温度的影响最小.      

基于最大似然估计的远紫外遥感反演电离层电子密度算法

原子氧135.6 nm夜气辉主要由氧离子O+与电子的辐射复合反应生成,一些星载远紫外遥感观测任务证实135.6 nm夜气辉可用于反演电离层电子密度。针对远紫外临边遥感观测反演电离层电子密度,分析了135.6 nm夜气辉辐射强度与电子密度之间的非线型前向模型,基于离散反演理论设计了从夜间135.6 nm临边观测数据反演电子密度高度分布的反演算法,算法应用最大似然估计通过迭代求解电离层参数的最佳拟合值。通过仿真计算了TIMED卫星上全球紫外成像仪GUVI观测的反演结果,验证了本反演算法的可行性。对GUVI的实际观测数据进行反演,获得了电子密度高度分布。通过与GUVI数据的电离层参数对比分析得出,本文建立的反演模型使N_mF₂被高估,同时使h_mF₂被低估。对于不同的太阳活动强度,N_mF₂和 h_mF₂的系统误差分别在10%和5%以内,能较精确地获得电离层参数。精确获得电离层电子密度信息对于提高空间天气预报及电离层模型的修正具有重要意义。      

相路径与群路径同斜向传输距离的关系

本文从电离层为均匀球面分层的假设, 得出了一组联系斜传播中相路径P、群路径P'和传输距离D的微分关系式, 导出了P～D和P'～D曲线的斜率公式, 并简短地讨论了这些曲线的可能应用.      

非相干雷达探测低电离层的结果分析

位于波多黎各的Arecibo非相干雷达可以获得低电离层电子和离子密度, 利用此非相干雷达数据对中纬度低电离层的运动特征进行研究. 得到了电子密度随时间和高度的变化 情况, 结果显著呈现出周日变化特征, 并分析了电子密度随高度的变化规律. 进一步对数据进行频谱分析, 深入研究低电离层电子密度的周日变化效应. 得到电子密度的高度剖面, 发现从F层底部到E层有明显的等离子体沉降. 低电离层的层结构特征及电子密度变化表明, 在该区域还存在不同程度的等离子体扰动, 由此对低电离层的作用因素 进行分析, 认为大气潮汐或声重波可能对低电离层产生扰动, 即低电离层与大气存在一定程度的耦合作用.      

TIEGCM集合卡尔曼滤波同化模型设计及初步试验

选择参数化的电离层热层理论模型TIEGCM作为背景模型,基于COSMIC掩星观测的电子密度廓线数据,应用集合卡尔曼滤波方法建立全球电离层电子密度同化模型,实现了全球电离层的电子密度同化.同化结果表明,该同化模型能将观测资料有效同化到背景模式中,获得全球三维电离层电子密度.与背景模式相比,同化得到的电子密度相对于观测值的偏差显著下降.对于有同化和无同化参与的试验,N_mF₂的标准偏差分别降低约60%和20%.此外,分组同化与同时同化的结果对比显示,平均偏差改善基本一致,同时同化后的标准偏差在峰值高度以上略有减小.      

电离层Es对HF和VHF传播信号的影响

Es层是存在于电离层中的电子密度非常高的偶发E层（Sporadic E）,其电子密度可达常规E层的100倍.电离层Es能够反射原本穿透F层的VHF低频段（30~150MHz）无线电波,而且对HF高频段（10~30MHz）无线电波传播具有显著影响.运用垂测和斜测观测数据,研究HF频段Es层电波传播特征,得到了不同类型及不同高度Es层的衰减系数.根据f₀E_s的日变化规律,可得HF频段衰减系数的日变化规律,进而分析并得到Es层对短波传播的影响.不存在电离层Es时,通常无法通过电离层实现VHF超远距离通信.为了对VHF链路通过电离层Es的传播衰减进行定量分析,根据EBU多条链路的观测结果,拟合并建立了电离层Es衰减模型.将该模型、ITU模型和观测数据进行对比,发现本文建立的模型准确度更高.利用建立的模型,对电离层Es不同临频f₀E_s条件下接收信号场强和电压随传播距离的变化进行了计算,结果可为VHF链路设计及建立提供参考.      

昆明上空电离层D区电子密度季节变化的首次观测分析

通过分析2008年8月至2009年7月昆明站(25.6°N, 103.8°E) 中频(MF)雷达观测数据, 研究了太阳活动低年电离层D区电子密度的季节变化特性,发现D区电子密度主要呈现半年变化特征, 即在春秋季电子密度较大, 而在夏冬季则较小, 这与国际参考电离层(IRI)预测的年变化趋势不一致, 但与昆明站电离层测高仪的最低回波频率f_min的观测结果相符. 同时比较了D区电子密度半年变化与纬向风半年变化的关系, 发现二者之间保持了非常一致的变化趋势并对这种一致性的内在原因进行了分析.      

Numerical study of traveling ionosphere disturbances with vertical incidence data

High frequency ionosphere vertical sounding, as an important and representative application for detecting the ionosphere and studying the characteristics of radio propagation, can be utilized to monitor the ionosphere continuously variation and to acquire the ionosphere asymmetrical features of diverse scale above the ionosphere vertical sounding stations. This is a first article on real time application of numerical methods to obtain the parameters of traveling ionosphere disturbances (TIDs) using vertical incident ionograms. In this paper, the distribution of ionosphere electron density with TIDs is constructed using a background ionosphere model superimposed a perturbation theory model. The background ionosphere electron density is modelled by the inversion of vertical incident ionograms which are observed before the appearance of the disturbance. Based on the fourth order Adams-Bashforth-Moulton (the so-called ABM) predictor corrector method, instead of Runge-Kutta method, the fast digital ray tracing method is established. According to process of the disturbed trace simulation and parameters inversion, the characteristic parameters of ionosphere disturbance at different detection time can be obtained in real time. The numerical analysis of TIDs is then captured completely.     

Estimation of ground range on the sweep frequency backscatter leading edge

The high frequency management system with backscatter radar supplies the real time ionosphere channel conditions to high frequency users, which leads to the demand for the ground range between the radar location and the scatters on the distant ground. The ionosphere electron density profile is usually inversed to obtain the ground range. An inversion algorithm, with which the ground range on the leading edge of the backscatter ionograms can be obtained without electron density, is presented in this paper. The ray path geometry of the backscatter sounding and the change in the group path on the leading edge with operating frequency are used to derive the ground range. Synthesized backscatter ionogram and experimental backscatter ionograms are processed to validate the algorithm. The results indicate that the algorithm is usable for high frequency management system.     

Location of layered structures in the ionosphere and atmosphere by use of GPS occultation data

A method is introduced to locate the layered structures in the atmosphere and ionosphere based on simultaneous observations of radio wave intensity and phase variations in trans-ionospheric satellite-to-satellite links. The method determines location of a tangent point on the trans-ionospheric ray trajectory where gradient of refractivity is perpendicular to the ray trajectory and influence of a layered structure on radio wave parameters is maximal. An estimate of the location of a layer can be obtained from a combination of the phase and intensity variations. This new technique was applied to measurements provided during FORMOSAT-3 and CHAMP radio occultation (RO) missions. For the considered RO events the location of the inclined plasma layer in the lower ionosphere is found and the electron density distribution is retrieved. The method is checked by measuring the location of the tangent point on the ray trajectory in the neutral gas in the atmosphere. The results showed a fairly good agreement.     

Physical and model interpretation of HF radio propagation on the St. Petersburg–Longyearbyen (Svalbard) path

HF radio wave observations have been carried out with an oblique ionospheric sounding (OIS) method on the radio path from St. Petersburg to Longyearbyen (Svalbard), and experimental ionograms were obtained for December 2001. These ionograms have been analysed to investigate the impact of the main ionospheric trough (MIT) and magnetic disturbances on the signals on this path. The observations during weakly disturbed (K_р = 2) magnetic conditions on 14–15 December 2001 were compared with predictions from ray-tracing through a numerical model of the ionosphere. The ray-tracing computer program synthesizes the OIS ionograms by means of the “shooting method”. This method calculates trajectories of HF radio waves for different values of elevation angle and transmission frequency. There was a variety of calculated trajectories, from which we choose those which reach the receiver, and the selected paths provide a synthesis of the oblique ionograms. To simulate HF radio wave propagation, we apply a three-dimensional distribution of the electron density calculated with the mathematical model of the high-latitude ionosphere developed in the Polar Geophysical Institute (PGI). These numerical simulations permit us to interpret specific peculiarities of the OIS data such as abnormal propagation modes, increased delays of signals, enhanced MOF (maximum observed frequency) values etc. New results of the study are summarised as follows. (1) An unusual feature of the propagation along the path is the change of propagation mechanism during substorms on entering a path midpoint (or 1-hop reflection point) to the MIT. (2) Even weak substorms, having the distinguished intensities, lead to the appearance of different types of irregularities observed by the CUTLASS radar and therefore to the different propagation modes and F2MOF values. (3) The PGI model of the ionosphere was first used for ray-tracing at high latitudes. The model results are basically in a good qualitative agreement with experimental observations. This model provides the satisfactory agreement between the calculated and experimental F2MOF values while not correctly representing the fine structure of the experimental OIS ionograms at night. An agreement between the calculated and experimental data is better for day and evening hours than at night.     

Adaptation of the bottomside electron density profile of the IRI model to data of topside radiosounding

A method is proposed for reconstructing the electron density profiles N(h) of the IRI model from ionograms of topside satellite sounding of the ionosphere. An ionograms feature is the presence of traces of signal reflection from the Earth's surface. The profile reconstruction is carried out in two stages. At the first stage, the N(h) –profile is calculated from the lower boundary of the ionosphere to the satellite height (total profile) by the method presented in this paper using the ionogram. In this case, the monotonic profile of the topside ionosphere is calculated by the classical method. The profile of the inner ionosphere is represented by analytical functions, the parameters of which are calculated by optimization methods using traces of signal reflection, both from the topside ionosphere and from the Earth. At the second stage, the profile calculated from the ionogram is used to obtain the key parameters: the height of the maximum hmF2 of the F2 layer, the critical frequency foF2, the values of B0 and B1, which determine the profile shape in the F region in the IRI model. The input of key parameters, time of observation, and coordinates of sounding into the IRI model allows obtaining the IRI-profile corrected to real experimental conditions. The results of using the data of the ISIS-2 satellite show that the profiles calculated from the ionograms and the IRI profiles corrected from them are close to each other in the inner ionosphere and can differ significantly in the topside ionosphere. This indicates the possibility of obtaining a profile in the inner ionosphere close to the real distribution, which can significantly expand the information database useful for the IRTAM (IRI Realmax Assimilative Modeling) model. The calculated profiles can be used independently for local ionospheric research.     

南京地区斜向加热电离层效应数值模拟

与垂直加热相比,斜向加热电离层具有灵活性高、影响范围广和便于实际操作等优势.在非偏区考虑欧姆吸收,基于电子能量方程和连续性方程构建电波斜向加热低电离层的物理模型,并利用IRI-2007和NRLMSISE-00经验模型提供的背景参数对南京地区斜向加热低电离层进行数值模拟,对比不同加热条件下电子温度和电子密度的扰动情况.研究结果表明:电波加热效果随入射仰角和有效辐射功率的增大而增大;电子温度和电子密度增幅随电波频率增大而减小;X波模比O波模造成的电子温度扰动幅度和电子密度扰动幅度更大,同时X波模比O波模能更快地使电子温度和电子密度达到稳定状态;一定范围内较小仰角、较低频率、较大有效辐射功率的电波能使电子密度更快达到稳定,后两者还能加快电子温度达到稳定的过程;电子温度达到稳定所需时间随入射电波仰角呈单峰变化,仰角为62°时达到最大.      

Modelling bottom and topside electron density and TEC with profile data from topside ionograms

The paper describes the technique that has been implemented to model the electron density distribution above and below the F2 peak making use of only the profiles obtained from the INTERCOSMOS-19 topside ionograms. Each single profile from the satellite height to the ionosphere peak has been fitted by a semi-Epstein layer function of the type used in the DGR model with shape factor variable with altitude. The topside above the satellite height has been extrapolated to match given values of plasmaspheric electron densities to obtain the full topside profile. The bottomside electron density has been calculated by using the maximum electron density and its altitude estimated from the topside ionogram as input for a modified version of the DGR derived profiler that uses model values for the foF1 and foE layers of the ionosphere. Total electron content has also been calculated. Longitudinal cross sections of vertical profiles from latitudes 50° N to 50° S latitude are shown for low and high geomagnetic activity. These cross sections indicate the equatorial anomaly effect and the changes of the shape of low latitude topside ionosphere during geomagnetic active periods. These results and the potentiality of the technique are discussed.     

用经验模式参数研究电离层发电机效应

本文利用MSIS-86和IRI-86模式的基本参数,求得热层风系和电导率的三维分布;继而从发电机理论得到了电离层电位及层电流密度在北半球的分布和变化。本研究将大量数据统计平均的大气模式参数与热层、电离层理论研究联系了起来。