A simple algorithm for regional mid-latitude ionospheric modelling

A simplified model of the key vertical-incidence ionospheric characteristics such as f₀ F2, M(3000) F2, h'F, f₀F1, and f₀E that are used for prediction of operational parameters of HF telecommunication system in restricted area has been developed. The model is based upon the Fourier analysis of the monthly median values of these characteristics from seven ionospheric stations in Europe. It is shown that only 12 dominant Fourier coefficeints, expressed as function of local time, solar activity and geographic latitude are sufficient to reproduce the main features of the diurnal, seasonal and solar cycle behaviour of the mid-latitude ionosphere under median conditions.     

Fe XVII X-ray lines in solar coronal plasmas

Calculated intensities of the Fe X-ray lines due to transitions 2p⁶ − 2p⁵3d lines (near 15 Å) and 2p⁶ − 2p⁵3s lines (near 17 Å) are compared with measured line intensities in solar and tokamak spectra. For the solar spectra, temperature T_e is obtained from the ratio of the Fe 16.776 Å line to a nearby Fe line. We find excellent agreement for all the major Fe line features in the 15–17 Å region except the Fe 15.015 Å line, the observed flux of which is less than the theoretical by a factor f. We find that f strongly depends on the heliocentric angle θ of the emitting region, being smallest (0.2) when the region is nearest Sun centre, but nearly 1 near the limb. Attributing this to resonance scattering, we are able to deduce the path length and electron density from the observations. Possible application to stellar active regions is given.     

The 2.5 to 5 microns spectrum of comet Halley from the IKS instrument of Vega

Results of the 2.5–5 micron spectroscopic channel of the IKS instrument on Vega are reported and the data reduction process is described. H₂O and CO₂ molecules have been detected with production rates of 10³⁰ s⁻¹ and 1.5 10²⁸ s⁻¹ respectively. Emission features between 3.3 and 3.7 microns are tentatively attributed to CH - bearing compounds - CO is marginally detected with a mixing ratio CO/H₂O 0.2. OH emission and H₂O - ice absorption might also be present in the spectra.     

A Prediction Method for△f0F2min of a Single Station From Interplanetary Parameters Based on Statistical Study

Using 86 CME-interplanetary shock events,the correlation between the peak values of （a） the solar wind parameters(B_z,E_y,P_dyn) and the geomagnetic indices（SYM-H,ASY-H,Kp）, （b） the coupling functions（Borovsky,Akasofu,Newell） and the geomagnetic indices,（c） the solar wind parameters/coupling functions/geomagnetic indices and the ionospheric parameter(Δf₀F_2min), are investigated.The statistical results show that in group（a）,B_{z min} and SYM-H_min have the best correlation,that in group（b）,the best correlation is between the peak values of Akasofu function (A_min) and SYM-H_min,and that in group（c）,the best correlation is between K_pmax andΔf₀F_2min. Based on the statistical results,a method for predicting f₀F₂ of a single station is attempted to be set up.The input is modified B_{z min} and the outputs are SYM-H_min andΔf₀F_2min.Then 25 CME-IPS events that caused geomagnetic storms in 1998 and 2009 are used to check the prediction method. The results show that our method can be used to predict SYM-H_min andΔf₀F_2min.     

The optical to X-ray spectrum of the broad-line ultrasoft AGN RE J2248-511

We present a series of monitoring observations of the ultrasoft broad-line Seyfert galaxy RE J2248-511 with XMM-Newton. Previous X-ray observations showed a transition from a very soft state to a harder state five years later. We find that the ultrasoft X-ray excess has re-emerged, yet there is no change in the hard power-law. Reflection models with a reflection fraction of 15, and Comptonisation models with two components of different temperatures and optical depths (kT₁ = 83 keV, T₁ = 30 eV, τ₁ = 0.8; KT₂ = 3.5 keV, T₂ = 60 eV, τ₂ = 2.8) can be fit to the spectrum, but cannot be constrained. The best representation of the spectrum is a model consisting of two blackbodies (kT₁ = 0.09 ± 0.01 keV, kT₂ = 0.21 ± 0.03 keV) plus a power-law (Γ = 1.8 ± 0.08). We also present simultaneous optical and infrared data showing that the optical spectral slope also changes dramatically on timescales of years. If the optical to X-ray flux comes primarily from a Comptonised accretion disk we obtain estimates for the black hole mass , accretion rate and inclination cos(i)  0.8 of the disk.     

Auroral signatures of transient processes in the outer magnetosphere

We present an analysis of sporadic and recurrent injections of magnetospheric ions in the midnight auroral oval during substorms and of the associated ionospheric ion outflows. The source of plasma sheet precipitating ions is determined using a simple method, based on the measured relation between the ion inverse velocity and time (l = v × t). This method is applied here to two typical passes of the Interball-Auroral (IA) satellite at distances of 3 R_E above the auroral regions. Substorm related ion injections are shown to be mainly due to time of flight effects. In contrast with particle trajectory computations (Sauvaud et al., 1999), the inverse velocity method does not require magnetic and electric field models and can thus be used systematically for the detection of time of flight dispersed ion structures (TDIS). This allowed us to build a large database of TDIS events and to perform a statistical analysis of their spatial distribution. For the cases presented here the source region of the injected ions is found at radial distances from 18 to 30 R_E near the equatorial magnetosphere. At Interball altitudes ( 3 R_E), ion injections detected at the poleward boundary of the nighside auroral oval are associated with shear Alfvén waves superimposed over large-scale quasi-static current structures. We show that the most poleward TDIS are collocated with a large outflow of ionospheric H⁺ and O⁺ displaying pitch-angle distributions peaked in the pitch-angle range 90°–120°. These ions are thus accelerated perpendicularly to the magnetic field not only in the main auroral acceleration region but also up to at least 3 R_E. The expanding auroral bulge thus constitutes a significant source of H⁺ and O⁺ ions for the mid-tail magnetosphere.     

The GPS measured SITEC caused by the very intense solar flare on July 14, 2000

This work studies the sudden increases in total electron content of the ionosphere caused by the very intense solar flare on July 14, 2000. Total electron content (TEC) data observed from a Global Positioning System (GPS) network are used to calculate the flare-induced TEC increment, δTEC_f, and variation rate, dTEC_f/dt. It is found that both dTEC_f/dt and δTEC_f are closely related with the solar zenith angles. To explain the observation results, we derived a simple relationship between the partial derivative of the flare-induced TEC, ∂TEC_f/∂t, which is a good approximation for dTEC_f/dt, and the solar zenith angle χ, as well as the effective flare radiation flux I_f, according to the well-known Chapman theory of ionization. The derived formula predicted that ∂TEC_f/∂t is proportional to I_f and inverse proportional to Chapman function ch(χ). This theoretical prediction not only explains the correlation of dTEC_f/dt and δTEC_f with χ as shown in our TEC observation, but also gives a way to deduce I_f from TEC observation of GPS network. Thus, the present work shows that GPS observation is a powerful tool in the observation and investigation of solar flare effects on the ionosphere, i.e., the sudden ionospheric disturbances, which is a significant phenomenon of space weather.     

视觉导航下基于H2/H∞的航迹跟踪

为了从理论上研究视觉导航系统的特性,首先在物理上将视觉导航系统划分为传感器(摄像头、激光测距仪)和视觉处理计算机。针对视觉处理计算机,建立视觉位置测量模型。将视觉位置测量模型的不确定性分为参数不确定性和输入干扰2种。然后通过仿真详细分析了视觉位置估计针对不同模型不确定性的敏感程度。最后在引入视觉位置测量模型不确定性的基础上研究H₂/H_∞混合鲁棒控制在航迹跟踪上的应用,并给出了保证航迹控制系统稳定性和跟踪性能情况下的视觉位置测量模型摄动范围,对计算机视觉算法的研究具有一定的指导意义。     

曲靖地区电离层多参数变化特征

利用电波环境观测网曲靖站电离层垂直探测仪数据（2008－2018）,分析了该站电离层多参数（包括f_min,f₀E_s,h'E_s,f₀E,f₀F₁,f₀F₂,h'F,MUF(3000)F₂等）随太阳活动、季节、地方时的变化特征。结果表明,f_min,f₀E_s,h'E_s随太阳活动变化不明显,f₀E,f₀F₁,f₀F₂,h'F,MUF(3000)F₂与太阳活动变化具有相关性,f₀F₂经常发生日落增强现象,同时各自具有独立性;与拉萨、乌鲁木齐地区f₀F₂,h'F,f₀E_s变化特征对比,发现地方时与月份变化趋势相同,但是极值出现的位置有所不同,这可能与各地区的地理纬度和地形等因素有关;h'F夜间大于白天,在日出日落时段有突然上升现象,冬季的h'F一般都小于其他季节;MUF(3000)F₂与f₀F₂的变化特征相似,白天值高于夜晚值,春秋分季高于夏冬季,太阳活动高年日落后MUF(3000)F₂一直维持较高值并持续到03:00－04:00 LT。      

On the adjust of minor constituents density in middle atmosphere during solar proton events

Profiles of O₃ partial pressure and of other minor atmospheric constituents (NO, NO₂, HCL, HF and H₂O), observed in the middle atmosphere during Solar proton events (20.04.1998; 05.04.2000), were analysed. Conclusions were drawn that under SCR impact a short-term O₃ partial pressure increase and destruction of some freon constituents took place.     

The dark matter halo structure of “fossil” groups candidates

We investigated properties of four isolated giant elliptical galaxies with extended X-ray halo using ASCA data. The derived size of X-ray halo, X-ray luminosity, and gravitational mass of the dark halo are unusually large those of X-ray halo of a single galaxy, but are typical for X-ray halos of groups and poor clusters of galaxies. The measured temperatures and abundances of the X-ray halo gas in these galaxies are also similar to those of the groups and poor clusters. Based on these results we identified these galaxies as “isolated X-ray overluminous elliptical galaxy” (IOLEG). The radial profiles of dark halo in these objects were derived from X-ray data. It is found that some are similar to those of compact groups while others are the same as those of normal ellipticals. The dark halos of lOLEGs are thus indistinguishable from those of groups (and poor clusters), which appears to be consistent with a widely believed idea that lOLEGs are a product of dynamical evolution of a compact group. However, mass-to-light ratios of IOLEGs (M₂₀₀/L_B  100–1000) are far greater than those of Hickson compact groups M₂₀₀/L_B  40–60). Since it is hard to consider that total optical luminosity of a compact group decreases by an order of magnitude in the course of dynamical evolution, such difference in the observed mass-to-light ratio between IOLEGs and Hickson compact groups strongly suggests that most IOLEGs have not evolved from compact groups which are observed at present.     

Intensity profiles of the multiple scattering light near the cometary nucleus

The multiple scattering of solar radiation in the cometary atmosphere is treated with the method of successive scattering. Referring to in situ measurements of comet Halley about the size and spatial distributions of dust, the optical thickness τ₁ of dust has been estimated, i.e. τ₁=0.03 at wavelength λ=0.62μm in a quiet time, but τ₁=0.3 when the outbursts/jets occur. In the derivation of τ₁, optical properties of dust including a mixing ratio of absorbing to silicate grains, are determined based on the polarimetry of P/Halley at λ=0.62μm observed during the phase angles over Nov. 1985 to May 1986 at the Dodaira Station of Tokyo Astronomical Observatory.

It is found that a temporary enhancement of τ₁ leads an increase of the upward reflected intensity when the surface albedo A of the nucleus is less than 0.04, but the reverse is true when A>0.04. On the other hand, the intensity of the downward radiation at the surface of the nucleus always decreases as an increase of τ₁.     

Neutral upper atmospheres of Venus and Mars

Numerous measurements of the neutral upper atmosphere above 100 km have been made from spacecraft over Venus and over Mars. The Venus exospheric temperatures are unexpectedly low (less than 300°K near noon and less than 130°K near midnight). These very low temperatures may be partially caused by collisional excitation of CO₂ vibrational states by atomic oxygen and partially by eddy cooling. The Venus atmosphere is unexpectedly insensitive to solar EUV variability. On the other hand, the Martian dayside exospheric temperature varies from 150°K to 400°K over the 11-year solar cycle, where CO₂ 15-μm cooling may be less effective because of lower atomic oxygen mixing ratios. On Venus, temperature increases with altitude on the dayside (thermosphere), but decreases with altitude from 100 to 150 km on the nightside (cryosphere). However, dayside Martian temperatures near solar minimum for maximum planet-sun distance and low solar activity are essentially isothermal from 40 km to 200 km. During high solar activity, the thermospheric temperatures of Mars sharply increase. The Venus neutral upper atmosphere contains CO₂, O, CO, C, N₂, N, He, H, D and hot nonthermal H, O, C, and N, while the dayside Mars neutral upper atmosphere contains CO₂, O, O₂, CO, C, N₂, He, H, and Ar. There is evidence on Venus for inhibited day-to-night transport as well as superrotation of the upper atmosphere. Both atmospheres have substantial wave activity. Various theoretical models used to interpret the planetary atmospheric data are discussed.     

Interplanetary shocks and geomagnetic activity during solar maximum (2000) and solar minimum (1995–1996)

Plasma and magnetic field parameter variations through fast forward interplanetary shocks were correlated with the peak geomagnetic activity index Dst in a period from 0 to 3 days after the shock, during solar maximum (2000) and solar minimum (1995–1996). Solar wind speed (V) and total magnetic field (B_t) were the parameters with higher correlations with peak Dst index. The correlation coefficients were higher during solar minimum (r² = 56% for V and 39% for B_t) than during solar maximum (r² = 15% for V and 12% for B_t). A statistical distribution of geomagnetic activity levels following interplanetary shocks was obtained. It was observed that during solar maximum, 36% and 28% of interplanetary shocks were followed by intense (Dst  −100 nT) and moderate (−50  Dst < −100 nT) geomagnetic activity, whereas during solar minimum 13% and 33% of the shocks were followed by intense and moderate geomagnetic activity. It can be concluded that the upstream/downstream variations of V and B_t through the shocks were the parameters better correlated with geomagnetic activity level, and during solar maximum a higher relative number of interplanetary shocks can be followed by intense geomagnetic activity than during solar minimum. One can extrapolate, for forecasting goals, that during a whole solar cycle a shock has a probability of around 50% to be followed by intense/moderate geomagnetic activity.     

基于拉曼努金和的非均匀多载波调制系统

针对高速移动通信,多普勒效应和多径效应导致信道非均匀,提出一种基于拉曼努金和的频谱可调的非均匀调制多载波系统.首先,根据拉曼努金和的正交性和周期性,推导了拉曼努金傅里叶正反变换的完全重建条件,进而建立基于拉曼努金和的多载波调制系统（RFMT,Ramanujan Fourier Multi-Tone system）.由于拉曼努金和的非均匀频谱分布性质及频率共振性质,RFMT在不同载波通道的误比特率不同,可实现对数据的不均等保护.在非均匀信道下RFMT具有优于OFDM的抗多径能力,仿真验证了RFMT的抗多径有效性.在使用迫零均衡算法时,10-5误比特率下,根据信道情况设计的非均匀多载波系统RFMT对E_b/N₀的要求可以比OFDM低4 dB.      

国际参考电离层模型输出参数峰值高度性能

电离层峰值高度H_mF₂是描述电离层形态的重要参数之一,国际参考电离层模型IRI-2016中融入了大量电离层测高仪和无线电掩星探测数据,用以提升H_mF₂的预测精度.本文利用太阳活动低年（2007—2010年）气象、电离层和气候卫星联合观测系统COSMIC探测数据描述全球范围内COSMIC H_mF₂的三维形态变化,对比分析了IRI-2016与IRI-2012模型的预测结果,同时分析了IRI-2016模型输出H_mF₂的性能.结果表明,IRI模型在中高纬度地区的输出结果高于COSMIC反演结果,而赤道及低纬地区则大都偏低.与IRI-2012模型相比,IRI-2016模型的输出结果在夜间至凌晨时段呈现较为明显的纬向梯度变化且大部分区域输出值偏高,但在白天时段赤道附近区域的输出值大都偏低.上述结果为电离层IRI模型的完善提供了一定参考.      

Towards validation of SCIAMACHY lunar occultation NO2 vertical profiles

Vertical profiles of stratospheric nitrogen dioxide (NO₂) have been retrieved from moderate resolution lunar occultation transmission spectra measured by Scanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) on board the European Environmental Satellite (ENVISAT). These measurements were taken over the high southern latitude of 50°–90° during the period of 2003–2005. To assess the accuracy of the retrieved NO₂ profiles, the SCIAMACHY nighttime NO₂ profiles were compared with NO₂ profiles retrieved from sunrise solar occultation spectra measured by the Halogen Occultation Experiment (HALOE) and the Stratospheric Aerosol and Gas Experiments II (SAGE II) using a photochemical correction model. The validation results show good agreement of SCIAMACHY lunar occultation NO₂ with scaled HALOE and SAGE II profiles. The relative mean differences (rmd) with scaled HALOE profiles are within −13% to +5% and standard deviations (rms) of the relative differences are within 3–19% between 25 and 38 km. The rmd and rms with scaled SAGE II NO₂ profiles are in the range of −9 to +7 and 10–17% respectively between 22 and 39 km.     

Cometary atmosphere (gas and dust)

There is important progress now in the identifications and measurements of primary (parent) molecules in the inner coma of Comet Halley. H₂O, CO₂ and CO are definitely in the list, CH and some complicate organic molecules are suspected. Gas production rate for water vapor is Q_H2O 10³⁰ s⁻¹. The bulk of data doesn't contradict to the Whipple model of nucleus (with clathrate modification). Pronounced spatial structure of gaseous flow in the coma was observed, but in general measured properties of neutral gas in the coma of Comet Halley are not very different from predicted. Situation for dust is different. In situ dust measurements show that size spectrum and optical properties of particles in coma are substantively declining from predicted on the base of groundbased photometry. However there are discrepancies between Vega and Giotto dust counter data. Dust in the inner coma didn't prevent the succesful imaging of nucleus by TV on Vega 1 and 2.     

SHAKING: Adjusted spherical harmonics adding KrigING method for near real-time ionospheric modeling with multi-GNSS observations

Precise positioning based on Global Navigation Satellite System (GNSS) technique requires high accuracy ionospheric total electron content (TEC) correction models to account for the ionospheric path delay errors. We present an adjusted Spherical Harmonics Adding KrigING method (SHAKING) approach for regional ionospheric vertical TEC (VTEC) modeling in real time. In the proposed SHAKING method, the VTEC information over the sparse observation data area is extrapolated by the Adjusted Spherical Harmonic (ASH) function, and the boundary distortion in regional VTEC modeling is corrected by the stochastic VTEC estimated using Kriging interpolation. Using real-time GPS, GLONASS and BDS-2/3 data streams of the Crust Movement Observation Network of China (CMONOC), the SHAKING-based regional ionospheric VTEC maps are re-constructed over China and its boundary regions. Compared to GNSS VTECs derived from the independent stations, the quality of SHAKING solution improves by 13–31% and 6–33% with respect to the ASH-only solution during high and low geomagnetic periods, respectively. Compared to the inverse distance weighting (IDW) generated result, significant quality improved of SHAKING-based VTEC maps is also observed, especially over the edge areas with an improvement of 60–80%. Overall, the proposed SHAKING method exhibits notable advantage over the existing regional VTEC modeling techniques, which can be used for regional TEC modeling and associated high-precision positioning applications.     

Proposal of new models of the bottom-side B0 and B1 parameters for IRI

The time series of hourly electron density profiles N(h) obtained from 27 ionosonde stations distributed world-wide have been used to obtain N(h) average profiles on a monthly basis and to extract the expected bottom-side parameters that define the IRI profile under quiet conditions. The time series embrace the time interval from 1998 to 2006, which practically contains the entire solar cycle 23. The Spherical Harmonic Analysis (SHA) has been used as an analytical technique for modeling globally the B0 and B1 parameters as general functions on a spherical surface. Due to the irregular longitudinal distribution of the stations over the globe, it has been assumed that the ionosphere remains approximately constant in form for a given day under quiet conditions for a particular coordinate system. Since the Earth rotates under a Sun-fixed system, the time differences have been considered to be equivalent to longitude differences. The time dependence has been represented by a two-degree Fourier expansion to model the annual and semiannual variations and the year-by-year analyses of the B0 and B1 have furnished nine sets of spherical harmonic coefficients for each parameter. The spatial–temporal yearly coefficients have been further expressed as linear functions of Rz12 to model the solar cycle dependence. The resultant analytical model provides a tool to predict B0 and B1 at any location distributed among the used range of latitudes (70°N–50°S) and at any time that improves the fit to the observed data with respect to IRI prediction.