Experimental determination of the main F-layer parameters on a global scale using data from a remote optical probe onboard the satellite Interkosmos-Bulgaria-1300

The optical complex of the satellite Interkosmos-Bulgaria-1300 includes the UV-spectrometer “Photon-1”, and the scanning spectrometer EMO-5 operating in the visible spectrum. Using the intensities of the OI 1356 and 6300 nm emission as measured by these two spectrometers, a theoretical evaluation of the two main F-layer parameters has been made: hmF, the height of the F2-maximum and NmF, the electron density at this height, together with the constant Ao characterizing the Ne(h) distribution in the IRI. The MSIS and CIRA-72 models are used for determining the distribution of the neutral component. If we had used CIRA-72, lower NmF values and higher hmF values had been obtained than with MSIS. The results obtained show that the method is very sensitive to the assumed neutral model, particularly when determining hmF. This work demonstrates the possibility of probing the F-region near its peak on a global scale using comperatively cheap and easily available optical methods.     

Properties of the neutral density and composition in the thermosphere

Measurements of the density and composition of the thermosphere between 150 and 500 km, which were obtained by the S3-1 satellite, have been compared with the Jacchia and MSIS models. The measurements of the densities of O, N₂, N and Ar show some differences from the current models which should be considered during the preparation of the next CIRA model. The Ar measurements are particularly useful in examining the response of the neutral atmosphere to geomagnetic heating. These results are useful in establishing the appropriate lower boundary conditions for modeling of the thermosphere.     

基于轨道预报的热层模式评估

在太阳活动高低年的地磁平静/扰动环境下,利用不同热层大气模式J77,DTM78,MSIS00,JB2008和CHAMP加速度计反演密度,分析有无先验信息条件下的轨道预报误差.结果表明无先验信息的精密轨道预报中,热层模式的性能可能被弹道系数等参数偏差干扰,此时预报误差不能作为模式性能的评价标准.先验信息对轨道预报精度提升非常明显,尤其是地磁扰动期先进热层模式性能得以展现,轨道预报误差为无先验信息情况下的10%～25%.目前热层模式的主要缺陷存在于地磁扰动期.各模式之间的差异是:JB2008模式可以通过线性和单一频率周期项补偿,而J77及DTM78等模式还存在更多频率的误差.本文对不同情况下精密轨道预报的研究结果可为空间碎片碰撞预警等工程实践提供参考.      

Evaluation of recent atmospheric density models

Precision radar tracking data taken in 1985 on two spherical satellites has been used to evaluate atmospheric density models: Jacchia 1971 [CIRA 1972], Jacchia 1977, and MSIS 1983, using Cook's definition of Cd /2/. The satellites have perigee heights of 270 km and 780 km. For each day, using numerical integration, an independent trajectory was computed, including a drag scale factor (S). The difference of S from unity represents the error in the atmospheric model used to compute the drag. At 270 km the three models performed equally well, as S was consistently close to unity. This was true for all hour angles and latitudes. At 780 km the scale factors ranged from 0.1 to 1.8 for all models. However, for the Jacchia 1977 model, the average scale factor was 0.943, and represents the best model at that altitude. Computer timing was also done. For orbit computation, drag models need further improvement, and the Jacchia 1977 model seems the best available today in terms of accuracy and computational efficiency.     

冬季西北地区临近空间气象火箭探测数据分析

在获取冬季西北地区一次临近空间气象火箭探测数据后,将火箭探测温度、密度与MSIS00模式和TIMED/SABER卫星数据进行对比,并将火箭探测风场与HWM07模式和MERRA再分析资料进行对比,分析火箭探测温度误差组成,计算各项温度修正量。结果表明:火箭、卫星、MSIS00模式获取的温度和密度随高度整体变化趋势一致;相对于MSIS00模式,火箭和卫星实测数据能够反映出更多的变化细节,且二者在细节上具有较多一致性。火箭实测风场与MERRA的一致性较好,而与HWM07模式差异较大,在平流层中部火箭探测风场明显强于HWM07模式。相对于HWM07模式和MERRA,火箭探测风场能够体现更多细节,在22 km和45 km附近均探测到较强的风切变。在火箭探测温度的各项修正量中,气动加热、温度滞后、支撑结构热传导及测量电流焦耳效应带来的影响较大,该影响整体上随着高度降低而逐渐减小。分析表明,本次气象火箭获取的探测数据是有效可靠的,但在数据处理方法尤其是温度误差修正等方面还需不断迭代完善。      

Theoretical model of electron temperature and electron density: Comparisons with IRI and MSIS

An earlier theoretical model (UW-87) accurately predicted the electron temperature in the daytime F-region but suggested N₂ concentrations significantly greater than the predictions of MSIS-86. This discrepancy is resolved when the model is developed to include the effects of vibrationally excited nitrogen molecules and electronically excited oxygen ions on the F-region recombination rate. The revised model (UW-92) continues to predict electron temperatures close to the layer peak with great accuracy but it is now more closely consistent with MSIS. However, the electron temperatures predicted by this model, which are in close agreement with EISCAT observations, are significantly higher than the values predicted by the international Reference Ionosphere.     

Neutral temperatures from Thomson scatter measurements: Comparisons with the CIRA(1972)

Neutral exospheric and lower thermospheric (100–130 km) temperatures from Thomson scatter measurements at Millstone Hill (42°N) are compared with CIRA temperatures with a view towards identifying deficiencies in the CIRA and recommending revisions. CIRA models the observed diurnal mean temperatures (T₀) to within 10% over a wide range of solar conditions (75? F_10.7 ? 250), but consistently underestimates the diurnal temperatures with maximum deviations approaching 50% of observed amplitudes (180–240 K) at solar maximum (1200 K ? T₀ ? 1400 K). The observed semidiurnal amplitudes, which lie in the range of 20K–80K, are always underestimated and frequently by more than 50%. In the lower thermosphere, tidal oscillations of temperature of order 20K–40K occur which are not modelled by CIRA. In addition, an analysis of exospheric temperatures at Millstone Hill during a magnetic disturbance indicates a response within 1–2 hours from storm onset, whereas CIRA assumes a 6.7 hour delay. Although some of these deficiences are addressed by the more recent MSIS model, there exists a sufficient data base to recommend several additional revisions to the CIRA temperatures at this time.     

Semi-annual variation of the thermospheric density

Thermospheric densities were determined from the orbital data of four satellites covering a height range of 266 to 500 km during 1975 to 1982. Using three different thermospheric models (DTM, MSIS, C), the index characterizing the semi-annual density variation was computed and analysed with respect to height and time using a regression by Fourier series. A general formula for the semi-annual density variation was derived. The new formula gives results similar to CIRA 72 at lower heights but displays a greater range of amplitudes at higher altitudes and a shift of extrema with growing height towards later days.     

Comparison between the IRI ion composition and incoherent scatter measurement and theoretical values

Comparisons have been made between the percentage of light ions in the upper ionosphere as predicted by the IRI model and as found in incoherent scatter (ICS) measurements at the stations Millstone Hill, Arecibo and Jicamarca. Major discrepancies are observed in both day and night. The IRI values are always considerably larger than the ICS measurements. Theoretical values are calculated as well, assuming chemical equilibrium and using the MSIS neutral density model /1/. In most cases these theoretical values favour the ICS values; only for the daytime ion composition above Millstone Hill has better agreement with the IRI model been found.     

Recent advances in upper atmospheric structure

A possible quantitative explanation of the semi-annual variation in thermospheric density has been obtained in terms of a semi-annual variation in the computed globally averaged vertical energy carried by propagating tides from the lower and middle atmosphere into the thermosphere. The effect is primarily due to seasonal changes in the distribution of water vapor and in the solar declination angle and Sun-Earth distance. An MSIS-83 empirical model of the thermosphere, representing a revision of the earlier MSIS models, has been prepared. The database used covers a wider range of solar activity than previous models and an improved magnetic storm representation is included. Atomic oxygen profiles in the 100 to 160 km altitude region of the auroral thermosphere have been recalculated from measured quenching of N₂(A³∑_u⁺) using the latest laboratory rates and the results are in good agreement with the mean CIRA 1972 profile. A new empirical model of thermospheric variations with geomagnetic activity has been developed incorporating variations with local magnetic time, latitude dependent terms which can vary with the magnitude of the geomagnetic disturbance, and an altitude dependent expression for the equatorial wave. A new index ML, derived from the AL index, has been developed that appears to have promise to represent the variations of thermospheric species with geomagnetic activity. Satellite measured values of solar UV flux, ground-based observations of CaK plages, sunspot numbers and 10.7 cm solar radio flux have been analyzed for temporal variations. Some differences have been identified and the significance to empirical and theoretical upper atmosphere models is discussed.     

Intra-annual variations of the thermospheric density at 400 km altitude from 1996 to 2006

We investigate the intra-annual variations of globally averaged thermospheric density at 400 km altitude from 1996 to 2006 by using Artificial Neural Network Method (ANNM). The results indicate that thermospheric density is governed by solar activity, and the absolute error of our model is 13.67%, less than NRLMSISE-00 model. Fourier representation can catch the intra-annual variations more accurately than NRLMSISE-00 model and JB2008 model especially during 2002. We find that the Autumn maximum is slightly greater than Spring maximum during solar minimum, while the reverse is correct during solar maximum. There is a strong linear relation between solar activity and the amplitude of annual/semiannual variations, and the correlation coefficients are 0.9534 and 0.9424, respectively. Moreover, the amplitude ratio of the annual to semiannual variation is about 1.3 averaged, and changes in different years, but it has little relation with solar activity. Besides that, the amplitude of annual variation is larger than semiannual variation during 1996 and 2006 except 1998 and 2000. The relative error of NRLMSISE-00 model is 14.95%, decreasing to 12.49% after revising, and the correlation coefficients between this empirical model and its improved results and the observation are 0.8185 and 0.9210, respectively. Finally, we suggest the revised version of MSIS series of model should use the Fourier representation to express the intra-annual variations.     

A Partially Orthogonal EnKF approach to atmospheric density estimation using orbital debris

A key requirement for accurate trajectory prediction and space situational awareness is knowledge of how non-conservative forces affect space object motion. These forces vary temporally and spatially, and are driven by the underlying behavior of space weather particularly in Low Earth Orbit (LEO). Existing trajectory prediction algorithms adjust space weather models based on calibration satellite observations. However, lack of sufficient data and mismodeling of non-conservative forces cause inaccuracies in space object motion prediction, especially for uncontrolled debris objects. The uncontrolled nature of debris objects makes them particularly sensitive to the variations in space weather. Our research takes advantage of this behavior by utilizing observations of debris objects to infer the space environment parameters influencing their motion.The hypothesis of this research is that it is possible to utilize debris objects as passive, indirect sensors of the space environment. We focus on estimating atmospheric density and its spatial variability to allow for more precise prediction of LEO object motion. The estimated density is parameterized as a grid of values, distributed by latitude and local sidereal time over a spherical shell encompassing Earth at a fixed altitude of 400 km. The position and velocity of each debris object are also estimated. A Partially Orthogonal Ensemble Kalman Filter (POEnKF) is used for assimilation of space object measurements to estimate density.For performance comparison, the scenario characteristics (number of objects, measurement cadence, etc.) are based on a sensor tasking campaign executed for the High Accuracy Satellite Drag Model project. The POEnKF analysis details spatial comparisons between the true and estimated density fields, and quantifies the improved accuracy in debris object motion predictions due to more accurate drag force models from density estimates. It is shown that there is an advantage to utilizing multiple debris objects instead of just one object. Although the work presented here explores the POEnKF performance when using information from only 16 debris objects, the research vision is to utilize information from all routinely observed debris objects. Overall, the filter demonstrates the ability to estimate density to within a threshold of accuracy dependent on measurement/sensor error. In the case of a geomagnetic storm, the filter is able to track the storm and provide more accurate density estimates than would be achieved using a simple exponential atmospheric density model or MSIS Atmospheric Model (when calm conditions are assumed).     

日照边缘区域电离层对耀斑的响应特点研究

利用MSIS模型和背景太阳辐射谱模型,在一定大耀斑辐射谱假设的前提下,计算了耀斑期间日照边缘区域的电子产生率,分析了这一区域电离层电子密度的变化特点.结果表明,大耀斑期间在日照边缘区域,甚至大于太阳天顶角90°的区域都有明显的电子产生率的增加.从不同太阳天顶角处的电子产生率剖面的形态来看,随着天顶角的增加最大电离率减少,但高度增加.计算还显示了在太阳天顶角小于90°的区域内电子产生率的垂直分布有明显的双峰结构,这种结构对应着电离层的E区和F区,但在天顶角大于90°区域,F区的电子产生率要大得多.考虑到离子和电子的复合过程,这一区域的总电子含量的增加主要产生在高F区.      

Ozone reference models for CIRA

The ozone reference model which is to be incorporated in the COSPAR International Reference Atmosphere (CIRA) is described and compared with other measurements of the Earth's ozone distribution. Ozone vertical structure models from approximately 25 to 90 km are provided combining data from five recent satellite experiments (Nimbus-7 LIMS, Nimbus-7 SBUV, AE-2 SAGE, Solar Mesosphere Explorer (SME) UVS, and SME IR). The results include the latest improvements in the SBUV algorithms using the most recent estimates of ozone cross sections. Also, the latest refinements in SME algorithms are incorporated. These algorithm improvements have improved agreement between the satellite data sets. Standard deviations are provided of monthly zonal means, and an estimate of the interannual variability is given. The models based on satallite data compare well with the Krueger and Minzner mid-latitude model incorporated into the U. S. Standard Atmosphere which is based on rocket and balloon measurements. Other comparisons are shown with Umkehr and more recent balloon data. Models are also provided of total columnar ozone reflecting recent improved estimates of ozone cross section. Information is provided on semiannual and annual variations. Other systematic variations including estimates of diurnal variations in the mesosphere will be included in the CIRA document.     

Comparison of measured and predicted F2 peak altitude

Different models for the F2 peak altitude hmF2 are compared with mean values determined from incoherent scatter measurements of the radar stations Jicamarca/Peru and Arecibo/Puerto Rico at all local times. Our investigation shows that the daytime hmF2 is fairly well represented by the most recent models, whereas at nighttime, higher peak altitudes are measured than predicted by either model. The observed diurnal structure is slightly misinterpreted by the models above Jicamarca and is strongly disagreeing above Arecibo.     

Modelling of ionospheric temperature profiles

The amount of measured temperature data accumulated in recent years allows and asks for improvement and refinement of the rather crude temperature models employed in the International Reference Ionosphere 1979. By combining the mission-oriented models by BRACE, THEIS /2/ for the AE-C satellites and by SPENNER, PLUGGE /1/ for the AEROS-A satellite, a much better diurnal and latitudinal reliability can be obtained. It is also suggested that IRI should have the option to make use of the strong anti-correlation between electron temperature and density in cases where actual measured densities are available. For daytime condition, incorporation of a density dependent model into IRI can significantly enhance the prediction quality of IRI in the altitude range 300 to 600 km.Furtheron, the solar activity dependence of the electron temperature and of the density-temperature-relation are investigated by comparing the low solar activity models with more recent AE-C and -DE data.     

Variability study of ionospheric total electron content at crest of equatorial anomaly in China from 1997 to 2007

The variation of TEC data at Wuhan station (geographic coordinate: 30.5°N, 114.4°E; geomagnetic coordinate: 19.2°N, 183.8°E) at crest of equatorial anomaly in China from January 1997 to December 2007 were analyzed. Variability with solar activity, annual, semiannual, diurnal and seasonal variation were also analyzed. The MSIS00 model and ISR model were used to analyze the possible mechanisms of the variabilities found in the results. The TEC data in 1997 and 2001 deduced from another crest station Xiamen (geographic coordinate: 24.4°N, 118.1°E; geomagnetic coordinate: 13.2°N, 187.4°E) were used to contrast. Analysis results show that long-term variations of TEC at Xiamen station are mainly controlled by the variations of solar activities. Typical diurnal variation behaves as a minimum of the TEC in the pre-dawn hours around 05:00–06:00LT and a maximum on the afternoon hours around 13:00–15:00LT. Some features like the semiannual anomaly and winter anomaly in TEC have been reported. The anomaly may be the result of common action of the electric field over the magnetic equatorial and the [O/N₂] at the crest station.     

变比冲连续小推力轨迹优化方法综述

电推进是目前常用于深空探测任务的推进方式。电推进发动机的比冲是变化的,所以变比冲电推进模型更加符合工程实际,对于变比冲连续小推力轨迹优化方法的研究具有工程应用前景。常用的变比冲电推进简化模型有三种：比冲在区间内变化、比冲随距离变化和功率分档变化,它们都对应于真实发动机系统工作的原理和方式。关于此三种模型下的轨迹优化方法,常用的有直接法和间接法,但是由于存在方法的不足和模型的特殊性,会在直接法和间接法的基础上增加改进的策略。本文综述变比冲电推进的简化模型和适用于变比冲连续小推力轨迹优化设计的方法。     

Cosmic-ray modulation and the heliospheric magnetic field

The heliospheric magnetic field plays a key role in any model for the modulation of cosmic rays. It enters into all diffusion coefficients, and its magnitude, spatial gradient and direction determine drifts patterns of cosmic rays in the heliosphere. While the first axisymmetric model of E.N. Parker proved quite successful to explain in situ measurements in the ecliptic plane, new insight into the origin and the nature of the field, especially at high heliographic latitudes, has led to the development of complex fully-three-dimensional, time-dependent models. In this review, we discuss a selection of models for the heliospheric magnetic field, and discuss how some of the more recent Fisk-type models affect the modulation of cosmic rays.