NRLMSISE-00大气模型与GRACE和CHAMP卫星大气密度数据的对比分析

利用GRACE(Gravity Recovery And Climate Experiment)和CHAMP(Challenging Mini-Satellite Payload)卫星2002-2008年的大气密度数据与NRLMSISE-00大气模型密度结果进行比较,分析了模型密度误差及其特点.结果显示,NRLMSISE-00大气模型计算的密度值普遍偏大,其相对误差随经纬度变化,在高纬度相对较小;相对误差随地方时变化,在02:00LT和15:00LT左右较大,10:00LT和20:00LT左右较小.通过模型密度相对误差与太阳F_10.7指数的对比分析发现,在太阳活动低年模型相对误差最大,而在太阳活动高年相对误差较小;将模型结果分别与GRACEA/B双星和CHAMP卫星的密度数据进行比较,发现对于轨道高度更高的GRACE卫星轨道,模型相对误差更大;在地磁平静期,相对误差与地磁ap指数(当前3h)相关性不强,但是在大磁暴发生时,误差急剧增大.      

太阳极紫外辐射对热层大气密度的影响

对2001-2021年SOHO卫星的极紫外辐射测量数据,以及CHAMP,GRACE-A和SWARM-C卫星资料推导出的高分辨率大气密度数据进行统计分析,发现大气密度与极紫外测量值的相关系数大于密度与F10.7指数的相关系数,证实极紫外辐射在不同地方时的影响程度存在显著差异,从而驱动大气密度的周日变化。利用三颗卫星的高度差异揭示极紫外辐射对大气密度的加热效应在350～500 km范围随着高度增加而减弱。统计得到极紫外辐射影响在地方时和纬度上的空间差异：对夏季半球的影响大于冬季半球;在白天,对中纬度地区的影响高于赤道和高纬度地区;在夜间,密度对辐射的斜率在夏季半球高纬度地区存在峰值,在冬季半球中纬度存在谷值,模型DTM2000和NRLMSISE00未能准确刻画。为了改进经验模型,提出基于球谐函数的拟合方法,优于主流模型周日效应采用的表达式,对周日效应建模和修正提供有益借鉴。利用昼夜间能量传输和热层大气经向环流机制探讨了统计结果的物理机制。     

Longitudinal variations of the ionospheric trough position

For the first time a comprehensive pattern of the longitudinal effect of the ionospheric trough position was obtained. We present new results with longitudinal variations of the winter trough position as a function of geomagnetic latitude for both hemispheres and conditions of high and low solar activity and all local time hours. We used a large observational data set obtained onboard the Kosmos-900, Interkosmos-19 and CHAMP satellites for quiet geomagnetic conditions. We found that a magnitude of the trough position longitudinal effect averaged for a fixed local time is greater in the daytime (6–8°) than in the nighttime (3–5°). The longitudinal effect magnitude reaches its maximum (16°) in the morning (at 08 LT) in the Southern hemisphere at high solar activity. But on certain days at any solar activity the longitudinal effect magnitude can reach 9–10° even at night. The shape of the longitudinal effect was found to differ significantly in two hemispheres. In the Northern hemisphere the trough is usually closest to the pole in the eastern (American) longitudinal sector, and in the Southern hemisphere the trough is closest in the western (Eurasian) longitudinal sector. The magnitude and shape of the longitudinal effect is also different during low and high solar activity. The Global Self-consistent Model of the Thermosphere, Ionosphere, and Protonosphere (GSM TIP) simulations demonstrate that during low solar activity, the longitudinal variations of the daytime trough position is mainly determined by longitudinal variations of the ionization function, formed due to the longitudinal variations in the solar zenith angle and the atomic oxygen density distribution. The longitudinal variations of the nighttime trough position is formed by the longitudinal variations in ionization of precipitating auroral particles, neutral atmosphere composition, and electric field.     

磁扰和磁静期间海南扩展F起始时间出现率研究

利用2002年2月至2007年12月(第23太阳活动周的下降段)近6年的海南DPS-4型电离层测高仪探测数据, 对磁扰和磁静夜晚期间扩展F起始时间出现率进行统计研究. 结合海南电离层观测站所观测到的扩展F类型, 将扩展F区分为频率型、区域型、混合型和强区域型, 分别进行统计分析. 结果表明, 无论磁扰还是磁静夜晚, 混合型扩展F起始时间总出现率最高, 最为活跃, 其次为频率型和强区域型扩展F, 最不活跃的是区域型扩展F; 无论在磁扰还是在磁静夜晚条件下, 混合扩展F起始时间主要围绕在午夜前后, 且磁静时更多地起始于午夜前, 而磁扰时则倾向于延至午夜后, 频率型扩展F在午夜后较高, 而强区域型扩展F则在午夜前较高. 在本次太阳活动下降阶段, 强区域型与区域型扩展F的起始时间出现率逐年与太阳活动呈现一定的相关性. 所得结果有助于分析不同类型扩展F在形态和机制方面的差异.      

Variation in the ionospheric propagating factor M(3000)F2 at Ouagadougou,Burkina Faso

We use hourly monthly median values of propagation factor M(3000)F2 data observed at Ouagadougou Ionospheric Observatory (geographic12.4°N, 1.5°W; 5.9^o dip), Burkina Faso (West Africa) during the years Januar1987–December1988 (average F10.7 < 130 × 10⁻²² W/m²/Hz, representative of low solar flux conditions) and for January 1989–December1990 (average F10.7 ? 130 × 10⁻²² W/m²/Hz, representative of high solar epoch) for magnetically quiet conditions to describe local time, seasonal and solar cycle variations of equatorial ionospheric propagation factor M(3000)F2 in the African region. We show that that seasonal trend between solar maximum and solar minimum curves display simple patterns for all seasons and exhibits reasonable disparity with root mean square error (RMSE) of about 0.31, 0.29 and 0.26 for December solstice, June solstice and equinox, respectively. Variability Σ defined by the percentage ratio of the absolute standard deviation to the mean indicates significant dissimilarity for the two solar flux levels. Solar maximum day (10–14 LT) and night (22–02 LT) values show considerable variations than the solar minimum day and night values. We compare our observations with those of the IRI 2007 to validate the prediction capacity of the empirical model. We find that the IRI model tends to underestimate and overestimate the observed values of M(3000)F2, in particular, during June solstice season. There are large discrepancies, mainly during high solar flux equinox and December solstice between dawn and local midnight. On the other hand, IRI provides a slightly better predictions for M(3000)F2 between 0900 and 1500 LT during equinox low and high solar activity and equinox high sunspot number. Our data are of great importance in the area of short-wave telecommunication and ionospheric modeling.     

The solar cycle variation of plasma parameters in equatorial and mid latitudinal areas during 2005–2010

Based on the ISL data detected by DEMETER satellite, the solar cycle variation in electron density (Ne) and electron temperature (Te) were studied separately in local daytime 10:30 and nighttime 22:30 during 2005–2010 in the 23rd/24th solar cycles. The semi-annual, annual periods and decreasing trend with the descending solar activity were clearly revealed in Ne. At middle and high latitudes, there exhibited phase shift and even reversed annual variation over Southern and Northern hemisphere, and the annual variation amplitudes were asymmetrical at both hemispheres in local daytime. In local nighttime, the annual variations of Ne at south and north hemispheres were symmetrical at same latitudes, but the annual variation amplitudes at different latitudes differed largely, showing obviously zonal features. As for Te, the phase shift in annual variations was not as apparent as Ne with the increase of latitudes at Southern and Northern hemisphere in local daytime. While in local nighttime the reversed annual variations of Te were shown at low latitudinal areas, not at high latitudes as those in Ne. The correlation study on Ne and Te illustrated that, in local daytime, Ne and Te showed strong negative correlation at equator and low latitudes, but during the solar minimum years the correlation between Ne and Te changed to be positive at 25–30° latitudes in March 2009. The correlation coefficient R between Ne and Te also showed semi-annual periodical variations during 2005–2010. While in local nighttime, Ne and Te exhibited relatively weak positive correlation with R being about 0.6 at low latitudes, however no correlation beyond latitudes of 25° was obtained.     

A broad climatology of very high latitude substorms

Magnetic data from a newly commissioned Indian Antarctic station Bharati (corrected geomagnetic (CGM) coordinates 74.7°S, 97.2°E) and closely-spaced IMAGE chain observatories (∼100° magnetic meridian in Northern hemisphere) has been analyzed to study the climatology of substorms which were localized poleward of the standard auroral oval. We considered four austral summers (year 2007–2010) when data from Bharati was available. Several very high latitude substorms were observed in this duration when the solar activity remained unexpectedly low for a long time. Various features of very high latitude substorms, e.g., local time dependence, interplanetary state, hemispherical asymmetry and their nightside low latitude signatures are examined. Events studied here, suggested the following properties of substorms occurring at very high latitudes: (1) maximum occurrence was observed near magnetic midnight (21:00–02:00 MLT). (2) In contradiction to earlier reports, many substorms were observed even during negative IMF Bz condition. In addition, majority of substorms occurred during low or moderate solar wind streams. (3) Magnetic signatures were often pronounced in the winter hemisphere. (4) Even if widely used standard AE indices fail to monitor very high latitude substorms, their low latitude signatures are often evident.     

The asymmetrical features in electron density during extreme solar minimum

The variations of plasma density in topside ionosphere during 23rd/24th solar cycle minimum attract more attentions in recently years. In this analysis, we use the data of electron density (Ne) from DEMETER (Detection of Electromagnetic Emissions Transmitted from Earthquake Regions) satellite at the altitude of 660–710 km to investigate the solstitial and equinoctial asymmetry under geomagnetic coordinate system at LT (local time) 1030 and 2230 during 2005–2010, especially in solar minimum years of 2008–2009. The results reveal that ΔNe (December–June) is always positive over Southern Hemisphere and negative over northern part whatever at LT 1030 or 2230, only at 0–10°N the winter anomaly occurs with ΔNe (December–June) > 0, and its amplitude becomes smaller with the declining of solar flux from 2005 to 2009. The ΔNe between September and March is completely negative during 2005–2008, but in 2009, it turns to be positive at latitudes of 20°S–40°N at LT 1030 and 10°S–20°N at LT 2230. Furthermore, the solstitial and equinoctial asymmetry index (AI) are calculated and studied respectively, which all depends on local time, latitude and longitude. The notable differences occur at higher latitudes in solar minimum year of 2009 with those in 2005–2008. The equinoctial AI at LT 2230 is quite consistent with the variational trend of solar flux with the lowest absolute AI occurring in 2009, the extreme solar minimum, but the solstitial AI exhibits abnormal enhancement during 2008 and 2009 with bigger AI than those in 2005–2007. Compared with the neutral compositions at 500 km altitude, it illustrates that [O/N₂] and [O] play some roles in daytime and nighttime asymmetry of Ne at topside ionosphere.     

第23太阳活动周热层大气密度对太阳辐射指数F10.7响应的模拟研究

利用NCAR-TIEGCM计算了第23太阳活动周期间（1996—2008年）400km高度上的大气密度,并统计分析大气密度对太阳辐射指数FF_10.7的响应.结果表明,在第23太阳活动周内,大气密度的变化趋势与太阳辐射指数FF_10.7的变化趋势基本一致,但是大气密度在不同年份、不同月份对太阳辐射指数FF_10.7的响应存在差异.第23太阳活动周内太阳辐射极大值和极小值之比大于4,而大气密度的极大值与极小值之比则大于10.太阳辐射低年的年内大气密度变化不到2倍,而太阳辐射高年的年内大气密度变化可达2倍甚至3倍.大气密度与FF_10.7指数在北半球高纬的相关系数比南半球高纬的相关系数大.在低纬地区,太阳辐射高年大气密度与FF_10.7指数的相关系数比低年的大.不同纬度上,大气密度与太阳辐射指数FF_10.7的27天变化值之间的相关系数都大于其与81天变化值之间的相关系数.      

Climatology of global,hemispheric and regional electron content variations during the solar cycles 23 and 24

We present the results of study on the variations of ionospheric total electron content (TEC) by using global, hemispheric, and regional electron contents computed from the global ionospheric maps (GIMs) for the period from 1999 to 2020. For a low and moderate solar activity, the global and regional electron contents vary linearly with solar 10.7 cm radio flux and EUV flux. While a saturation effect in the electron content verses EUV and F10.7 is found during the high solar activity periods at all regions, the maximum effect is observed at low-latitudes followed by high and mid-latitudes region. The extent of saturation effect is more pronounced for F10.7 as compared to EUV. A wavelet transform is applied to global and hemispheric electron contents to examine the relative strength of different variations. The semi-annual variations dominate in the northern hemisphere, whereas annual variations dominate in the southern counterpart. The amplitude of annual variations in southern hemisphere is found to be higher than northern counterpart at all latitudes. This asymmetry in the amplitude of annual variation is maximum at low-latitudes, followed by mid and high-latitudes, respectively. The semi-annual variations are in-phase in both hemisphere and follow the solar cycle. The northern hemisphere depicts relatively large amplitude of semi-annual variations and exhibit the maximum effect at high-latitudes.     

Cross-hemisphere comparison of mid-latitude ionospheric variability during 1996–2009: Juliusruh vs. Hobart

Analysis of a long-time series of hourly median characteristics of the ionospheric plasma at two mid-latitude locations in the Northern and Southern hemisphere, Juliusruh (54.6N; 13.4E) and Hobart (42.9S; 147.3E), reveals patterns of their synchronous and independent variability. We studied timelines of GPS vTEC, ionogram-derived F2-layer peak electron density NmF2, ionospheric equivalent slab thickness τ, and their ratios at two locations during the complete 23rd solar cycle and its following period of the extremely low solar activity in 2008–2009. This study has also involved the comparative analysis of the observed data versus the model predictions by IRI-2012. During the high solar activity in 2000–2002, seasonal variations show a complicated cross-hemisphere behavior influenced by the winter and semi-annual anomalies, with the largest noon-time values of TEC and NmF2 observed around equinoxes. Strength of the winter anomaly in NmF2 was significantly greater at Juliusruh in comparison with Hobart. The winter anomaly in GPS vTEC values was much weaker than in NmF2 for the Northern hemisphere mid-latitudes and was entirely absent at the Southern hemisphere. Cross-hemisphere analysis of the equivalent slab thickness shows its clear seasonal dependence for all levels of solar activity: the day-time maximum τ_max is observed during local summer, whereas the day-time minimum τ_min is observed during local winter. The night-time values of τ were higher compared to the day-time values during the winter and equinox seasons. Comparative model-data study shows rather good IRI performance of the day-time NmF2 for mid-latitudes of both hemispheres and rather noticeable overestimations for the mid-night NmF2 values during high solar activity. Analysis of IRI vTEC demonstrates the model limitations, related with the absence of the plasmaspheric part, and actual demand in a reliable and standard ionosphere–plasmasphere model for analysis of GPS vTEC.     

The equatorial ionization anomaly at the topside F region of the ionosphere along 75°E

Electron density measured by the Indian satellite SROSS C2 at the altitude of ∼500 km in the 75°E longitude sector for the ascending half of the solar cycle 22 from 1995 to 1999 are used to study the position and density of the equatorial ionization anomaly (EIA). Results show that the latitudinal position and peak electron density of the EIA crest and crest to trough ratios of the anomaly during the 10:00–14:00 LT period vary with season and from one year to another. Both EIA crest position and density are found to be asymmetric about the magnetic equator and the asymmetry depends on season as well as the year of observation, i.e., solar activity. The latitudinal position of the crest of the EIA and the crest density bears good positive correlation with F_10.7 and the strength of the equatorial electrojet (EEJ).     

第24太阳活动周中国地区电离层闪烁统计特性研究

利用第24太阳活动周中国多个地区GNSS电离层闪烁监测站数据,统计分析中国中低纬地区电离层闪烁特性.结果显示:电离层闪烁主要发生在春秋分及夜间20:00—02:00LT时段;在28°N以南地区,纬度越低电离层闪烁强度和发生概率越高;电离层闪烁发生概率与太阳活动呈正相关,太阳活动上升年电离层闪烁发生概率高于下降年;不同强度地磁活动条件下,电离层闪烁均可能发生,且与地磁活动强度整体呈负相关.通过研究电离层闪烁统计特性,可以为电离层闪烁机理的深入研究、预报及工程应用提供参考.      

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.     

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.     

武汉地区电离层TEC和NmF2及板厚的季节变化

通过利用武汉电离层观测站(114.4°E,30.6°N)1980-1990年对E8T-Ⅱ卫星信标的法拉第旋转测量的TEC(电子浓度总含量)数据,以及由测高仪测量的1980-1990年间的f0F2(F2层临界频率)数据,研究了武汉地区TEC,NmF2(最大电子浓度)和板厚的季节变化,同时比较了IRI和武汉单站模式在预测NmF2季节性方面的有效性.武汉单站模式在预测NmF2季节性变化方面优于IRI模式.      

Studies of electron and ion temperatures at 500 km altitude during sunrise using Indian SROSS C2 satellite

Solar dependence of electron and ion temperatures (T_e and T_i) in the ionosphere is studied using RPA data onboard SROSS C2 at an altitude of ∼500 km and 77°E longitude during early morning hours (04:00–07:00 LT) for three solar activities: solar minimum, moderate and maximum during winter, summer and equinox months in 10°S–20°N geomagnetic latitude. In winter the morning overshoot phenomenon is observed around 06:00 LT (T_e enhances to ∼4000 K) during low-solar activity and to T_e ∼ 3800 K, during higher solar activity. In summer, it is observed around 05:30 LT, but the rate of T_e enhancement is higher during moderate solar activity (∼2700 K/hr) than the low-solar activity (∼1700 K/hr). During equinox, this phenomenon is delayed and is observed around 06:00 LT (∼4200 K) during all three activities.     

海口站与Huancayo站扩展F的差异

本文统计分析了我国海口站太阳黑子高年和低年的扩展F的资料,得出了该站扩展F的出现率随地方时、季节和太阳活动周期的某些变化规律并与Huancayo站扩展F的资料进行了比较,指出了两站扩展F的差异。      

基于FY-3C掩星数据偶发E层的研究

利用FY-3C极轨卫星提供的2014年6月至2015年5月的GPS无线电掩星数据,统计分析了全球范围内抽样频率为50 Hz的C/A码SNR扰动情况,进而对偶发E层进行了研究.结果表明:偶发E层在夏季半球中纬地区的扰动强度远远大于冬季半球同一纬度地区的扰动强度,偶发E层在纬度40°附近扰动明显增强;在E层100 km高度附近,Es层在10:00 LT和22:00 LT达到峰值;Es层在夏季半球的出现率明显高于冬季半球;FY-3C卫星的掩星观测结果与COSMIC系统的观测结果较一致,可以利用FY-3C卫星的掩星数据研究电离层偶发E层等的变化.      

嫦娥四号着陆区织女陨坑电场环境数值模拟

由于太阳紫外、X射线和太阳风粒子作用,近月表形成尘埃等离子体环境,而月表陨坑地形使得这种尘埃等离子体环境更为复杂.本文以位于嫦娥四号着陆区的织女陨坑为研究对象,基于高程数据构造了该陨坑的三维模型.根据太阳-月球关系和陨坑地理坐标信息,计算了陨坑白天任意时刻的有效太阳辐照度分布,探讨了不同时刻陨坑内外的光照面积占比,得到陨坑随地方时而发生的遮蔽效应特征.同时,基于月表充电方程计算了织女陨坑在不同地方时条件下的平衡表面电势、德拜鞘高度和电场强度分布,发现陨坑自身遮蔽效应对坑内电场环境影响十分明显.以坑底中心为例,讨论了地方时和纬度对类织女陨坑的平衡表面电势、德拜鞘高度及电场强度的影响,结果表明三者变化特征均以正午时刻及赤道为界呈对称分布,越接近12:00LT或者越接近赤道,坑底中心的平衡表面电势和电场强度越高,德拜鞘高度越低.