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

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

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

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

基于MARSIS的火星低纬度电离层统计特性研究

在对欧空局火星快车探测器搭载的MARSIS雷达的浅表层探测数据进行校准过程中,获得了火星电离层的总电子含量(total electron content,TEC)观测数据。利用该数据,计算火星低纬度地区电离层的峰值电子密度和标高;并对其进行统计分析发现,在低纬度地区,火星冬季电离层的标高和峰值电子密度均较夏季高,即冬季电离层较夏季更显著,且春季电离层的电子密度梯度最大。     

三频信标高精度TEC测量新方法

电子总含量(TEC)是电离层探测的主要参量之一, 作为层析(CT)的输入参量, TEC测量精度直接影响电离层CT成像的结果. 过去主要采用双频信标测量TEC, 由于相位积分常数的求解、系统硬件延迟等误差, 使得TEC测量结果不能满足电离层CT高精度重建成像的要求. 三频相干信标技术的出现, 使得电离层天基测量技术有了新的发展. 提出了基于三频信标的传播时延-相位联合测量反演TEC的方法, 融合三频信标在电子密度随机起伏探测和相位积分常数计算两方面的优势, 进一步提高了TEC的测量精度. 模拟结果显示利用此方法的三频信标TEC测量结果提高了电离层CT的精度.      

基于小波分解与重构方法研究电离层偶发E层

电离层不规则体对卫星导航、通信、雷达系统等有重要影响.通过数值模拟及与实测数据的对比,论证基于小波分解与重构方法实现利用掩星数据反演电离层不规则体的可行性.以电离层偶发E层为例,利用国际参考电离层(IRI)模拟背景电离层电子密度分布,利用掩星探测的水平电子密度总含量δht反演不规则体信息,并与模拟数据进行比较.对200...     

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

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

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

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

利用COSMIC掩星数据监测电离层的异常变化

分析了COSMIC掩星数据反演电子密度的方法,利用实例研究反演方法的特点,并采用ISR非相干散射雷达获取的电子密度数据进行验证,进而反演了长三角区域SHAO（IGS）站上空在日全食和太阳风暴期间的电子密度廓线图. 通过与平静日期间电离层电子密度进行比较,发现日全食及太阳风暴导致电离层发生的异常变化,从而提出COSMIC掩星数据反演电子密度在监测电离层变化时所具有的优势.      

电离层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链路设计及建立提供参考.      

武昌低电离层LF观测的季节变化

通过分析武昌低电离层LF观测数据,得到了中低纬度低电离层电子密度剖面的季节变化特征及其与太阳天顶角周年变化的关系,给出LF相位周年振荡幅度的年际变化曲线．     

化学物质释放人工改变电离层

考虑中性气体在电离层高度的扩散过程和相应的电离层离子化学过程，研究了利用主动化学物质释放来改变电离层的方法，理论计算了H2O和SF6两种气体释放后电离层随时间的响应过程，结果表明，在电离层高度上气体的扩散过程非常迅速，电离层F区的电子密度有很大程度的减少，而扩散慢且化学反应快的气体对电离层的影响更大，就更加有利于电离层洞的形成。     

电离层H2O释放扰动效应的数值模拟研究

为充分研究化学物质在电离层释放的扰动效应和后期发展效果,基于化学物质在电离层的扩散模型、化学反应和电离层扩展F的控制模型,通过电离层H₂O的释放,研究电子e,H₂O,O+和H₂O+共4种粒子的分布状态,分析点源、多源和线源释放对电离层的扰动效果,比较不同高度、不同量和不同时间释放的影响结果,模拟夜间释放后期所激发的扩展F发展差异.结果表明,H₂O在电离层释放后,能有效耗散背景电子形成空洞,O+和H₂O+数密度呈椭圆形分布;点源、多源和运动目标线源等不同释放方式对电离层的扰动效果不同,证实了人工影响一定形态和区域电离层的可能性;H₂O释放扰动幅度,低层大于高层,白天强于夜晚,释放量越多扰动越突出;夜间化学释放能激发扩展F,并且释放量越多,激发效果越好.      

Photochemical response of neutral and ionized middle atmosphere composition to the strong solar proton event of October 1989

One of the strongest solar proton events (SPE) occurred in October 1989. Its forcing of the middle atmosphere chemistry including ionized components in the D-region is examined. The ionization rate, and ozone, NO and OH density temporal and spatial (vertical) deviations induced by the SPE, calculated by a 1-D time-dependent photochemical model separately for daytime and nighttime (not shown here), are used in a 1-D model of the lower ionosphere to calculate the response of ionized components to combined forcing by ionization rate and neutral chemical composition disturbances. The radio wave absorption caused by electron density disturbances after the SPE is calculated and compared with observations. The computed ozone values are compared with observations, as well.     

Longitudinal distribution of thermospheric densities and ionization rates in the upper atmosphere of Mars at mid-latitude using MGS ACC data

The Accelerometer Experiment (ACC) onboard Mars Global Surveyor (MGS) measured 1600 density profiles in the upper atmosphere of Mars during aerobraking. These measurements reveal large-scale and small-scale structure in the thermosphere of Mars. Here, the measurements of mass density for 115 orbits (#P0670–P0789) from November 1 to 30, 1998, under spring equinox and medium solar activity conditions (average F_10.7 ∼ 137) during phase 2 of the aerobraking in the thermosphere of Mars at different altitudes and longitudes are presented for northern mid-latitude (17–42°N) in the dayside atmosphere using ACC onboard MGS. From these mass densities, the neutral densities of different gases are derived from their mixing ratios. Using these neutral densities, the longitudinal distribution of photoionization rates and photoelectron impact ionization rates are calculated at wavelength range 1–102.57 nm due to EUV and soft X-ray radiation under photochemical controlled region using Analytical Yield Spectrum approach (AYS). These conditions are appropriate for MGS Phase 2 aerobraking period from which the accelerometer data is used. Under the photochemical equilibrium condition, the electron density near the peak varies as the square root of the total peak ionization rate. Using this fact, an attempt is being made to estimate the mean primary and secondary peak electron density by averaging the longitudinal variations of total peak ionization rates in the northern mid-latitude (17–42°N) ionosphere of Mars, as there is no radio science measurement at this latitude region by MGS.     

Characterization of ionospheric amplitude scintillations using wavelet entropy detrended GNSS data

The extensive monitoring networks of Global Navigation Satellite System (GNSS) ionospheric scintillation have been established to continuously log observation data. Further, the amplitude scintillation index and the phase scintillation index, which are derived from scintillation observations, are anticipated to accommodate the accuracy requirement of both the user level and the monitoring station level. However, raw scintillation observations essentially measure superposed waveform impairments of GNSS signals propagating through ionosphere and troposphere. It implies that fluctuations of raw scintillation observations are caused by multiple factors from the entire radio propagation environment. Hence, it is crucial to characterize ionospheric scintillations from GNSS observation data. And the characterization is implemented through extracting fluctuations of raw observations merely induced by ionospheric scintillations. Designed to address this problem by means of Fourier filtering detrending, the present work investigates the influence of varying detrending cutoff frequencies on wavelet statistical energy and wavelet entropy distributions of scintillation data. It consequently derives criteria on the optimum detrending cutoff frequency for three types of raw amplitude scintillation data, which are classified by their wavelet energy distributions. Results of the present work verify that detrending with specific optimum cutoff frequencies rather than the fixed and universally applicable one renders the validity and credibility of characterizing ionospheric scintillations as the part of GNSS observation fluctuations purely induced by ionosphere electron density irregularities whose scale sizes are comparable with or smaller than the Fresnel scale.     

由高频无线电波反射回波参数反演电离层运动和结构的高度剖面

本文从广义电波射线微分方程组出发,导出了时变、不均匀和各向异性电离层介质中高频无线电波反射回波的可测参数与电离层的运动和结构参数之间的积分方程组,并给出了利用现代数字测高仪中的多普勒图、定向图、频高图和测距图,反演电离层的结构和运动随高度的分布剖面的计算方法.根据这种反演方法,单站探测电离层运动和结构剖面时有多种可测参量的组合方式.文中提出了一种最佳的参量组合方案,为探测设备的改进和观测记录的分析提供了理论依据.      

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.     

The D-region ion composition and electron density response to strong solar proton events (model simulations)

Fluxes of energetic solar protons penetrate deep into the Earth’s polar cap middle atmosphere. Interacting with molecules of the air they cause additional dissociation and ionization, and the formed NO_x, OH_y and ions enter chemical and ion-molecular reactions. Induced changes of the ionospheric D-layer are modeled by a 1D model of lower ionosphere with chemistry, using neutral species concentrations calculated by a 1D photochemical time-dependent model. Changes of the electron and ion densities, and the most important ionospheric parameters are calculated after SPE with the onset on July 14, 2000 and the results are compared with our results obtained previously for the October 19, 1989 SPE. It is shown that not only electron density increases after SPE, but also the amount of clusters. It is found that the magnitude of the ionospheric response depends on season.     

Regional 4-D modeling of the ionospheric electron density from satellite data and IRI

Accurate knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying ionosphere physics. During the last decade Global Navigation Satellite Systems (GNSS) have become a promising tool for monitoring ionospheric parameters such as the total electron content (TEC). In this contribution we present a four-dimensional (4-D) model of the electron density consisting of a given reference part, i.e., the International Reference Ionosphere (IRI), and an unknown correction term expanded in terms of multi-dimensional base functions. The corresponding series coefficients are calculable from the satellite measurements by applying parameter estimation procedures. Since satellite data are usually sampled between GPS satellites and ground stations, finer structures of the electron density are modelable just in regions with a sufficient number of ground stations. The proposed method is applied to simulated geometry-free GPS phase measurements. The procedure can be used, for example, to study the equatorial anomaly.