频高图Es层回波的自动识别

提出了一种基于决策树算法的自动识别Es层回波的方法以提高Es识别效率.通过自适应二值化以及中值滤波算法对频高图进行预处理，并提取有效的Es回波区域.利用Es层在虚高轴分布的规律，通过图像投影的方法选择出特征，将这些特征以及时间信息作为决策树算法的输入；通过人工识别方式为每幅频高图标明是否存在Es回波的标签，与决策树算法的输出进行比较.对训练集进行训练、剪枝，获取最接近人工识别结果的决策树，使用测试集对获取的决策树进行测试验证.本文使用云南普洱站（22.7°N，101.5°E）记录的频高图数据作为训练集，分别使用云南普洱站（22.7°N，101.5°E）和四川乐山站（29.5°N，103.7°E）的频高图数据作为测试集对该方法进行测试验证，结果表明该方法对普洱站和乐山站两地Es二跳回波的识别均具有较高的准确率，分别达到84.2%和82.8%.      

Seismogenic ionospheric anomalies associated with the strong Indonesian earthquake occurred on 11 April 2012 (M = 8.5)

Ionospheric perturbations in possible association with a major earthquake (EQ) (M?=?8.5) which occurred in India-Oceania region are investigated by monitoring subionospheric propagation of VLF signals transmitted from the NWC transmitter (F?=?19.8?kHz), Australia to a receiving station at Varanasi (geographic lat. 25.3°N, long 82.99°E), India. The EQ occurred on 11 April 2012 at 08:38:35?h UT (magnitude?≈?8.5, depth?=?10?km, and lat.?=?2.3°N, long.?=?93.0°E). A significant increase of few days before the EQ has been observed by using the VLF nighttime amplitude fluctuation method (fixed frequency transmitter signal). The analysis of total electron contents (TEC) derived from the global positioning system (GPS) at three different stations namely, Hyderabad (latitude 17.38°N, longitude 78.48°E), Singapore (latitude 1.37°N, longitude 103.84°E) and Port Blair (latitude 11.62°N, longitude 92.72°E) due to this EQ has also been presented. Significant perturbation in TEC data (enhancements and depletion) is noted before and after the main shock of the EQ. The possible mechanisms behind these perturbations due to EQ have also been discussed.     

北驼峰区电离层GPS卫星闪烁事件时空特征及对通信的影响

基于子午工程北大深圳站（22.59°N,113.97°E）电离层GPS双频接收机在2011年1月1日至2017年12月31日连续7年的长时间序列闪烁和TEC观测数据,分析不同太阳活条件下华南赤道异常北驼峰区观测到的GPS卫星L波段电离层闪烁事件时空分布特征及其对通信的影响.结果表明:GPS闪烁事件几乎都发生在夜间,且主要发生在春秋分月份;在不同太阳活动条件下,夜间GPS闪烁事件都主要发生在北驼峰区域靠近磁赤道的一侧,且GPS闪烁事件存在明显的东-西侧天区不对称性,即在台站西侧天区发生的闪烁事件明显偏多;在不同太阳活动条件下,弱闪烁事件伴随的TEC耗尽和卫星失锁事件比例相对较低,强闪烁事件则大部分都伴随着TEC耗尽和卫星失锁事件的发生.      

Radio aurora magnetic and streaming aspect sensitivities on 6 simultaneous links at 50 MHz

During August 1981, a 50 MHz c.w. radar system was operated in central Canada to measure auroral scatter amplitudes and Doppler spectra from a scattering region centered near 66° magnetic invariant latitude (L ≈ 6.0). Narrow beams from 3 transmitters, differing in frequency by 1 KHz, were directed to cover a common volume of the ionosphere over a ground location at 56.3°N, 103.5°W. The scattered signals were received on narrow beam antennas at two receiving sites, and recorded in analog form on magnetic tape under the control of an AIM65 microcomputer. The analog tapes were digitized later and FFT-processed to obtain Doppler spectra and amplitudes.The 6 transmission paths were designed to provide several magnetic aspect angles varying by 1.5°-7° from perpendicularity with the earth's field B and two streaming aspect angles differing by ～38°. The objective was to employ controlled geometric factors to study the functional dependency of signal amplitudes and Doppler shifts on magnetic and streaming aspect angles. Several hundred hours of excellent data were obtained in continuous operation during the month of August 1981. Preliminary results will be reported.     

利用北斗观测数据实时监测中国区域电离层变化

北斗卫星导航信号采用三个频点工作,可以利用伪距双频组差方法解算电离层电子含量,为实时监视中国区域电离层变化提供新的技术手段.中国中低纬度处于电离层赤道异常变化区,在北纬20°±5°区域时常发生较大梯度的电离层变化.利用北斗实时多频伪距和相位观测数据,采用相位平滑伪距方法计算电离层穿刺点电子含量,分析通过北斗系统GEO卫星监测的电离层周日变化特性;采用多面函数方法拟合中国区域1°×1°分辨率的电离层延迟量,每5min绘制一幅中国区域电离层图,观测区域所有电离层穿刺点拟合残差RMS为2.778TECU;分析北斗系统实时监测中国区域电离层异常情况,当发生电离层异常变化时,相邻两天的VTEC（Vertical Total Electronic Content）峰值相差约60TECU.      

Ducted and non-ducted propagation of Omega signals within the plasmasphere

Signals of the transmitter of the Omega navigation system located in the Southern Hemisphere (Reunion, 21.0°S, 55.3°E, L= 1.35) recorded by Interkosmos satellites on passes over the Northern Hemisphere, have shown that these signals can be detected in the conjugate region within a range of L-values between 1.5 – 3.0. Both ducted and non-ducted propagation has been observed. The most frequent propagation is in the whistler mode along the line of forceL= 1.8. The measured signal delays are in accordance with the observed dispersion of whistlers.     

Multi-wavelength coordinated observations of ionospheric irregularity structures from an anomaly crest location during unusual solar minimum of the 24th cycle

The present paper reports coordinated ionospheric irregularity measurements at optical as well as GPS wavelengths. Optical measurements were obtained from Tiny Ionospheric Photometer (TIP) sensors installed onboard the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) satellites. GPS radio signals were obtained from a dual frequency GPS receiver operational at Calcutta (22.58°N, 88.38°E geographic; geomagnetic dip: 32.96°; 13.00°N, 161.63°E geomagnetic) under the SCIntillation Network Decision Aid (SCINDA) program. Calcutta is located near the northern crest of Equatorial Ionization Anomaly (EIA) in the Indian longitude sector. The observations were conducted during the unusually low and prolonged solar minima period of 2008–2010. During this period, four cases of post-sunset GPS scintillation were observed from Calcutta. Among those cases, simultaneous fluctuations in GPS Carrier-to-Noise ratios (C/N_o) and measured radiances from TIP over a common ionospheric volume were observed only on February 2, 2008 and September 25, 2008. Fluctuations observed in measured radiances (maximum 0.95 Rayleigh) from TIP due to ionospheric irregularities were found to correspond well with C/N₀ fluctuations on the GPS links observed from Calcutta, such effects being noted even during late evening hours of 21:00–22:00 LT from locations around 40° magnetic dip. These measurements indicate the existence of electron density irregularities of scale sizes varying over several decades from 135.6 nm to 300–400 m well beyond the northern crest of the EIA in the Indian longitude sector during late evening hours even in the unusually low solar activity conditions.     

Precise orbit determination and precision comparison for FY3C and FY3D with receiver antenna GPS and BDS PCV using spaceborne BDS and GPS observation data

The FY3C and FY3D satellites were equipped with global navigation satellite occultation detector (GNOS) receivers that received both GPS and BDS-2 signals. For further improving precise orbit determination (POD) precisions, we estimated receiver GPS and BDS signal phase center variations (PCV) models with 2° and 5° resolutions and set the different weights for GPS and BDS-2 observations in the combined POD. The BDS-based POD precision using BDS-2 satellite antenna phase center offset (PCO) values from the China Satellite Navigation Office (CSNO) are not as accurate as those obtained from the International GNSS Service (IGS) Multi-GNSS experiments project (MGEX). The estimated receiver GPS and BDS PCV models with 2° and 5° resolutions were estimated from the GPS phase residuals of GPS-based POD and BDS phase residuals of combined POD, respectively. In most cases, the POD precisions using the estimated PCVs with 2° resolution are superior to those with 5° resolution. The precisions of the BDS-based POD and combined POD were both improved by introducing the receiver BDS PCV models. The weighting for GPS and BDS-2 observations can further improve the precision of the combined POD. The tested results of selected weights are better than those with equal weight in the combined POD. The experiment results show that orbital precisions of FY3C are worse than those of FY3D.     

The lower ionosphere response to its disturbances by powerful radio waves

The paper presents data from some campaigns at Sura heating facility in 2011–1016. The experiments on probing of the artificial disturbed region of the lower ionosphere were carried out at two observation sites. One of them was located near Vasil’sursk 1 km from Sura facility (56.1°N; 46.1°E) and the other site was located at the Observatory (55.85°N; 48.8°E) of Kazan State University, 170 km to the East. Investigation of the features of the disturbed region of the lower ionosphere based on its diagnostics by the methods of the vertical sounding and oblique backscattering is the main goal of this paper. Ionosphere disturbance was fulfilled by the effect of the powerful radio wave of the ordinary or extraordinary polarization emitted by transmitters of the Sura facility with effective radiated power ERP = 50–120 MW at the frequency of 4.3, 4.7 and 5.6 MHz. Pumping waves were emitted with period from 30 s to 15 min. The disturbed region of the ionosphere in Vasil’sursk was probed by the vertical sounding technique using the partial reflexion radar at the frequency of 2.95 and 4.7 MHz. For the oblique sounding of the disturbed region the modified ionosonde Cyclon-M, operating at ten frequencies from 2.01 to 6.51 MHz was used at the Observatory site. On many heating sessions simultaneous variations of the probing partial reflection signals in Vasil’sursk and backscattered signals in Observatory were observed at the height at 40–100 km below the reflection height of the pumping wave. These observations were correlated with the pumping periods of the Sura facility. Possible mechanisms of the appearance of the disturbance in the lower ionosphere and its effect on the probing radio waves are discussed.     

Analysis of coastal altimetry in the Mexican Caribbean

A reprocessing of sea-level anomalies (SLA) resulting from X-TRACK coastal altimetry was carried out for the ENVISAT (2002–2010) and TOPEX/POSEIDON-Jason (1992–2019) satellite missions in the coastal area of the Mexican Caribbean. This consisted of applying a tidal correction to coastal altimetry sea level observations. Harmonic analysis of five coastal tide gauge records was performed to estimate the most important tidal components of the area, resulting on M2, N2, O1, S2, K1, MF, and MM. The tidal signal was reconstructed with the seven tidal components using the TPXO9 model. The SLA signals corrected with the seven tidal components were validated with in situ data from coastal tide gauges. The validation showed that the TPXO9 tidal barotropic model (1/30° grid) used to reconstruct the tidal signal with the seven representative tidal components performed better than the FES2012 global model (1/16° grid) that uses 33 tidal components. The reprocessed SLAs showed clear seasonality with significant signals at 4, 6, and 12 months, with the annual signal being the dominant one. In the Mexican Caribbean coastal zone, oceanographic processes with different scales (from coastal to mesoscale) converge, showing their complexity in the different SLA signals observed. The aim of this work is to contribute to the analysis of coastal altimetry data and understanding the sea level variations in the Mexican Caribbean. This work is the first step in the implementation of methodologies that take advantage of coastal satellite altimetry in the Caribbean Sea.     

协同无线通信系统中的分布式自适应功率分配

在系统总功率一定的情况下,通过在协同节点间进行优化的功率分配以改善系统的接收性能.针对单中继协同传输,推导了节点间最优功率分配因子的闭式表达式,给出了中继节点的激活条件,仅当信道的衰落特性满足该条件时,中继节点才发送信号到目的节点,否则进入空闲状态.多节点协同传输场景下,直接根据系统的输出信噪比难以得到节点间最优功率分配因子的闭式表达式.因此,推导了系统输出信噪比的上界,在此基础上,提出了一种新的功率分配方法.将整个功率分配过程分成两步,采用迭代的方法实现节点间的最优功率分配.仿真结果表明,对于单中继节点的情况,所提出的功率分配方案能够获得约1~2 dB的性能增益,随着参与协同传输的节点数的增大,迭代的功率分配可使系统获得显著的性能提升.     

Satellite availability and positioning performance of uncombined precise point positioning using BeiDou-2 and BeiDou-3 multi-frequency signals

To realize the smooth transition from regional BeiDou Navigation Satellite System (BDS-2) to the global one (BDS-3), the integration of BDS-2 and BDS-3 is important for providing continuous, stable and reliable positioning, navigation and timing (PNT) services for global users. This work used 154 globally distributed multi-GNSS (Global Navigation Satellite System) experiment stations spanning 30 days to analyze the satellite availability and positioning performance of uncombined precise point positioning (UC-PPP) under current BDS-2 and BDS-3 constellations. We focused on three issues: the influence of BDS-3 receiver tracking abilities, the positioning performance among different areas, and the benefit of multi-frequency observations. The results show that the elliptical zone caused by poor BDS-2 satellite visibility is disappeared when the evenly distributed BDS-3 medium earth orbit satellites are introduced. When BDS-3 are integrated with BDS-2, the area with the Position Dilution of Precision (PDOP) less than 2 can be expanded to 75° S-75° N and 30° E-150° W. The positioning performance of BDS-3 and BDS-2/BDS-3 UC-PPP are seriously affected by the receiver tracking abilities of BDS-3 signals. When the maximum pseudo-random noise sequences (PRNs) of BDS-3 satellites tracked by stations are within 30 or 37, the positioning accuracy of static UC-PPP can be improved by 22.94% or 8.27% due to the integration of BDS-2 and BDS-3. Besides, the most improvement of BDS-2 and BDS-3 integration is achieved in Asia-Pacific regions, especially for the kinematic UC-PPP or the poor receiver tracking abilities of BDS-3. Similar to the multi-frequency BDS-2 UC-PPP, the benefit of multi-frequency signals for BDS-3 or BDS-2/BDS-3 UC-PPP is also non-vital. The three-dimensional positioning accuracy of BDS-2/BDS-3 multi-frequency UC-PPP in static mode and kinematic mode are 2.24 cm and 5.39 cm, while the corresponding convergence time are 49.62 min and 73.80 min, respectively. Compared with BDS-2, both the positioning accuracy and the convergence time of BDS-2/BDS-3 joint UC-PPP are improved by approximately over 50%, which indicates that BDS-3 has a great potential to provide high-quality PNT services as other global navigation satellite systems.     

瑞利激光雷达反演中层大气温度算法

瑞利激光雷达是探测中层大气温度的重要设备,具有高时空分辨率且不存在探测盲区的优点.为充分利用激光雷达探测到的原始数据,改进了传统的Chanin-Haunchecorne方法,采用由均匀搜索生成温度初值的方式反演大气温度.位于北京延庆（40.3°N,116.2°E）的瑞利激光雷达具有589nm和532nm两个通道,将589nm通道用于计算,532nm通道作为参考.在150~250K区间等间隔选取多个温度初值,利用589nm通道反演60~70km高度范围内的大气温度廓线和大气密度廓线,利用参考密度廓线选取准确的温度初值,反演得到准确的大气温度廓线.将589nm通道和532nm通道反演的温度廓线相比较,发现二者具有较高的一致性.改进后的方法有效利用了信噪比较差的光信号,使589nm通道的温度探测上限从60km提高到70km.      

Observations of whistler-type echoes on signals of a ground VLF transmitter on board the ‘iterkosmos-19’ satellite

Signals of VLF transmitters of the Omega navigation system located in the auroral zone (66.4°N, 13.2°E, L= 5) were recorded by the VLF receiving equipment of the Interkosmos 19 satellite. Signals at frequencies between 10.2 and 13.6 kHz were received in a region above the transmitters, frequently with whistler-type echoes. An analysis of these echoes has shown their predominating occurrence in periods of low geomagnetic activity (Kp<2+). The occurrence region of these phenomena in the outer ionosphere has the dimension of about 1000 km and its position is betweenL= 2.5 and L= 4.4. The delay of echo-signals is practically the same during one satellite pass but its values for different satellite revolutions lie between 2.5 and 3.5 s. The frequency spectrum of these signals can be broadened up to 100 Hz. On the basis of calculations made, it can be shown that the experimental results are generally in accordance with the hypothesis of nonlinear ducting of VLF waves in the magnetosphere.     

Effects of X2-class solar flare events on ionospheric GPS-TEC and radio waves over Brazilian sector

In this investigation, we present and discuss the effects of 6 X2-class solar flare events in the ionospheric F region over Brazilian sector that occurred during 2013 to 2015. For this investigation, we present the vertical total electron content (VTEC) observations from nearly 120 Global Positioning System (GPS) receivers all over the Brazilian sector for each event. Also, ionospheric sounding observations obtained in São José dos Campos (23.2°S, 45.9°W, dip latitude 17.6°S; hereafter referred to as SJC), under the southern crest of the equatorial ionospheric anomaly (EIA), Brazil, are presented. The observations show that the greatest TEC impact occurs with the EUV fluxes increases lasting for more than one hour and when the solar active region is located close to the solar disc center. We present a detailed study of the efficiency of the EUV flux with wavelengths ranging from 0.1 to 190?nm for the F region ionization. The largest increase of ΔTEC occurs below the magnetic equator line, covering mainly the central, northeast, southeast and south regions, which includes the equatorial ionospheric anomaly (EIA) region. The ionograms show partial or total fade out in the echoes traces observed causing blackouts of radio signals of up to 60?min, which can have serious consequences to technological systems of public and private agencies around Brazilian sector. This study can help to better understand the effects of solar flares in the ionospheric F region.     

Temporal evolution of the EIA along 95°E as obtained from GNSS TEC measurements and SAMI3 model

The total electron content (TEC) derived from GNSS measurements at a trans-hemispheric meridional chain of ground stations around 95°E longitude are used to study the quiet time inter-hemispheric structure and dynamics of the equatorial ionization anomaly (EIA) during the period March 2015 to February 2016. The stations are Dibrugarh (27.5°N, 95°E, 43° dip), Kohima (25.6°N, 94.1°E, 39° dip), Aizawl (23.7°N, 92.8°E, 36° dip), Port Blair (11.63°N, 92.71°E, 9° dip) and Cocos Islands (12.2°S, 96.8°E, 43° dip). The observation shows that the northern crest of the EIA lies in the south of 23°N (Aizawl) in all seasons but recedes further south towards the equator during December solstice. The largest poleward expansion of the northern (southern) EIA is observed in the March equinox (December solstice). The equinoctial and hemispherical asymmetry of TEC is noted. The winter anomaly is observed in the northern hemisphere but not in the southern hemisphere. The highest midday TEC over any station is observed in the March equinox. The TEC in southern summer (December solstice) is significantly higher than that in the northern summer (June solstice). The observed northern EIA contracts equatorward in the postsunset period of solstice but the southern EIA persists late into the midnight in the December solstice. The asymmetry may be attributed to the different geographic location of the magnetically conjugate stations. The SAMI3 simulations broadly capture the EIA structure and the inter-hemispheric asymmetry during solstices. The difference between observations and the SAMI3 is higher in March equinox and December solstice. The higher E?×?B vertical drift in the 90–100°E sector and the large geographic-geomagnetic offset in observing stations may have contributed to the observed differences.     

Mean winds of the mesosphere and lower thermosphere (60–110 km): A global distribution from radar systems (MF,Meteor, VHF)

During the last decade a large number of radars (～12) have been developed, which have produced substantial quantities of tidally-corrected mean winds data. The distribution of the radars is not global, but many areas are well covered: the Americas with Poker Flat (65°N), Saskatoon (52°N), Durham (43°N), Atlanta (34°N), Puerto Rico (18°N); Europe with Kiruna (68°), Garchy (47°N) and Monpazier (44°N); and Oceania with Christchurch (44°S), Adelaide (35°S), Townsville (20°S), and Kyoto (35°N). Zonal and meridional wind height-time cross-sections from $\text{60}\text{80}$ km (MF/Meteor Radar) to ～110 km have been prepared for the last 5–6 years. They are compared with cross-sections from CIRA-72 for zonal winds, and Groves (1969) for meridional winds.It is shown that while CIRA-72 is still a useful model for many purposes, significant differences exist between it and the new radar data. The latter demonstrate important seasonal, latitudinal, longitudinal and hemispheric variations. The new meridional cross-sections are of great value. The common features with Groves (1969) are the equatorward cells in summer near 85 km; however their strength (～10 ms^?1) and size are less. Systematic and somewhat different variations emerge at higher (?52°N) and middle (35–44°) latitudes.     

On the occurrence and strength of multi-frequency multi-GNSS Ionospheric Scintillations in Indian sector during declining phase of solar cycle 24

This study presents unique perspectives of occurrence and strength of low latitude ionospheric scintillations on multiple signals of Global Navigation Satellite System (GNSS) and its frequency dependence using continuous observation records of 780 nights. A robust comparative analysis is performed using scintillation index, S4 and its variation during pre-midnight and post-midnight duration from a GNSS receiver located at Waltair (17.7°N, 83.3°E), India, covering period from July 2014 to August 2016. The results, generally exhibit the impact of declining phase of solar cycle 24 on occurrence and strength of scintillations, which, however, is evidently different over different frequencies transmitted from different GNSS systems. A deeper quantitative analysis uniquely reveals that apart from the solar cycle and seasonal effects, the number of visible satellites of a selected GNSS markedly affect the occurrence and also the strength. Processing scheme of adopting 6 hourly time windows of pre-midnight and post-midnight brought a novel result that the strength and occurrence of strong scintillations decrease with declining solar activity during pre-midnight hours but remarkably increase for moderate and weak scintillations during post-midnight. The physical processes that dominate the post-midnight equatorial ionosphere are invoked to explain such variations that are special during declining solar activity. Finally, inter-GNSS signal analysis in terms of the effect of strong, moderate and weak scintillations is presented with due consideration of number of satellite passes affected and frequency dependence of mean S4. The quantitative results of this study emphasize for the first time effect of low latitude scintillation on GNSS signals in Indian zone under changing background solar and seasonal conditions.     

Ionospheric perturbations due to earthquakes as determined from VLF and GPS-TEC data analysis at Agra,India

Employing SoftPAL receiver, amplitude variations of VLF transmitter signals NWC (19.8?kHz) and NPM (21.4?kHz) are analyzed at Agra station in India (Geograph. lat. 27.2°N, long. 78°E) ±15?days from five major earthquakes of magnitude M?=?6.9–8.5 occurred in Indian subcontinent during the years 2011–2013. We apply nighttime fluctuation (NF) method and show that in almost all cases the trend decreases and dispersion and NF increase on the same days corresponding to each earthquake about 11–15?days prior to the main shock. Assuming that the ionospheric perturbations are caused by atmospheric gravity waves (AGW), we also calculate AGW modulation index for each case and find its values increased on the days amplitude fluctuations take place. Its value is decreased in one case only where the perturbations may be attributed to penetration of seismogenic electric field.In order to support the above results we also present GPS-TEC data analyzed by us corresponding to three of the above earthquakes. We study the TEC anomalies (unusual enhancements) and find that in one case the precursory period is almost the same as that found in NF method.     

HXMT卫星空间环境监测器及初步观测结果

硬X射线调制望远镜（HXMT）卫星是中国首个专门进行天文探测的空间科学实验卫星,运行于高度约550km、倾角约43°的低地球轨道.星载空间环境监测器为星上科学任务开展提供背景辐射实测资料.该监测器采用固体探测器望远镜系统和扇形阵列全新组合设计,可获取轨道空间高能质子和高能电子能谱、方向综合动态结果,给出更为全面的粒子辐射分布图像.初步探测结果显示,卫星运行轨道遭遇的带电粒子辐射集中分布在经度80°W-20°E,纬度0°-40°S的南大西洋异常区,粒子辐射在该区域表现出不同程度的方向差异分布,高能电子方向差异分布显著强于高能质子.2017年9月空间环境扰动期间,爆发的太阳质子事件并未对该轨道粒子辐射产生影响,而地磁活动导致该轨道穿越经度120°W-60°E,纬度40°-43°N的北美上空和经度60°-120°E,纬度43°-40°S的澳大利亚西南区域时遭遇增强粒子辐射影响,增强的粒子辐射表现出极强的方向分布.