低轨卫星缺失时序数据的模式识别方法

在航天应用中，低轨卫星经常会由于原始数据缺失而影响卫星时序数据模式识别结果，降低准确率。针对该问题提出了一种新型MR-GRU模型，可有效处理缺失时序数据，并获得较好的模式识别准确率。区别于传统模型的补全缺失数据的方法，MR-GRU模型直接在缺失时序数据上运用循环神经网络进行训练，对传统门控循环单元结构进行了改进，增加了两个新变量：掩蔽项和衰减项。掩蔽项作用于输入，衰减项作用于输入和隐层单元输出。MR-GRU模型不仅能够保持时序数据固有的时间特性，还能有效提高模式识别精度。在卫星时序数据上的模式识别试验表明，MR-GRU模型准确率优于传统模型。     

A novel satellite image encryption algorithm based on hyperchaotic systems and Josephus problem

With the rapid growth of the number of Earth observation satellite (EOS) supporting critical applications, it is required to improve the security techniques to protect the sensitive data and images during the transmission between the satellites and the ground stations. This paper introduces a new satellite image encryption algorithm based on the Linear Feedback Shift Register (LFSR) generator, SHA 512 hash function, hyperchaotic systems, and Josephus problem. LFSR generates a matrix that is used to construct the 512-bits value of the hash function. These bits are used to set the initial values and parameters of the proposed encryption algorithm. Firstly, the six dimensions (6-D) hyperchaotic system is divided into three parts, where every two equations are considered as one part. Secondly, the 1-D hyperchaotic logistic-tent system is considered as the controller to select one part. The selected part is used to generate a matrix that is XORed with the original image. Thirdly, the scrambling operation by Josephus sequences is applied to the output of the previous step by scrambling the rows and the columns according to the selected part to produce the pre-encrypted image. Finally, if the number of iterations is less than the required number which is considered as a parameter of the secret key, the previous operations will be repeated in the pre-encrypted image; otherwise, the pre-encrypted image is considered as the final cipher image. Experimental and analyses results show that the proposed algorithm has good performance in terms of high level of security, large enough key-space, tolerance to Single Event Upsets (SEU) as well as low time complexity.     

Establishment of atmospheric weighted mean temperature model in the polar regions

Due to the special geographical location and extreme climate environment, the polar regions (Antarctic and Arctic) have an important impact on global climate change. Atmospheric weighted mean temperature (Tm) is a crucial parameter in the retrieval of precipitable water vapor (PWV) from the zenith wet delay (ZWD) of ground-based Global Navigation Satellite System (GNSS) signal propagation. In this paper, the correlation between weighted mean temperature and surface temperature (Ts) is studied firstly. It is shown that the correlation coefficients between Tm and Ts are 0.93 in the Antarctic and 0.94 in the Arctic. The linear regression Tm model and quadratic function Tm model of the Antarctic and the Arctic are established respectively using the radiosonde profiles of 12 stations in the Antarctic and 58 stations in the Arctic from 2008 to 2015. The accuracies of the linear regression Tm model, the quadratic function Tm model and GPT2w Tm model which is a state-of-the-art global Tm model are verified using the radiosonde profiles from 2016 to 2018 in the Antarctic and Arctic. Root Mean Square (RMS) errors of the linear regression Tm model, the quadratic function Tm model and GPT2w Tm model in the Antarctic are 3.07 K, 2.87 K and 4.32 K respectively, and those in the Arctic are 3.53 K, 3.38 K and 4.82 K, which indicates that the quadratic function Tm model has a higher accuracy compared to linear regression Tm model, and the accuracies of the two regional Tm models are better than that of GPT2w Tm model in the polar regions. In order to better evaluate the accuracy of Tm in the PWV retrieval, the PWV values of radiosondes are used for comparisons as the reference value. The RMS errors of PWV derived from the two Tm models are similar for 1.28 mm in the Antarctic and 1 mm in the Arctic respectively. In addition, the spatial and temporal variation characteristics of Tm are analyzed in the polar regions by spectral analysis of Tm data using fast Fourier transform. The results show that the Tm has obvious seasonality and annual periodicity in the polar regions, and the maximum difference between warm season and cold season is about 63 K. After comparing and analyzing the influences of latitude, longitude and elevation on the Tm in the polar regions, it is found that latitude and elevation have a greater influence on the Tm than the longitude. As the latitude and elevation increase, the Tm decreases, and vice versa in the polar regions.     

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

Identification of seismic activity sources on the subsatellite track by ionospheric plasma disturbances detected with the Sich-2 onboard probes

Using a cylindrical Langmuir probe and the authors’ proprietary two-channel pressure transducer, ionospheric plasma parameter distributions along the orbit of the Sich-2 satellite (Ukraine, 2011–2012) were measured. This paper is concerned with identifying the space–time location of ionospheric plasma disturbance sources, including the epicenters of actual earthquakes (before or during the satellite flyover) and incipient earthquakes on the subsatellite track, from the measured distributions of the electron density and temperature and the neutral particle temperature along the satellite orbit. To do this, the measured ionospheric plasma parameter distributions are connected to the coordinates on the subsatellite track.It is shown that local disturbances in the electron density and temperature and neutral particle temperature distributions in the satellite orbit in the ionosphere may serve as indicators of seismic activity on the subsatellite track. The epicenters of incipient earthquakes may be set off from other plasma parameter disturbance sources associated with seismic activity using information provided by special monitoring and survey centers that monitor the current seismic situation.     

Semi-annual variations of chlorophyll concentration on the Eastern Tropical Pacific coast of Mexico

An analysis of the main remotly sensed oceanographic variables was conducted to characterize the dominant drivers that modulate the spatial-temporal variability of coastal phytoplankton biomass in the northern limit of the Eastern Tropical Pacific. The phytoplankton biomass was analyzed using monthly average satellite chlorophyll (Chlo) concentration data from MODIS sensors for the period 2003–2017. Oceanographic conditions were analyzed by considering (i) Sea Surface Temperature (SST) high-resolution data from the GHRSST project, (ii) wind stress calculated with data from the zonal and meridional components of the CCMP product, (iii) climatological precipitation, and (iv) climatological river flow. Chlo time series and spatial variability were analyzed using Hovmöller diagrams and Empirical Orthogonal Functions (EOF), respectively. A strong semi-annual signal in Chlo concentration along the coast was observed: the first peak occurs in winter-spring (5.0 mg·m⁻³) and a second one in summer-autumn (6.5 mg·m⁻³). A high year-round average of Chlo concentration (3.0–15.0 mg·m⁻³) was maintained by a 10 km wide along-shore coast, with maximum values spatially associated with river mouths, of which Santiago river registered the highest Chlo values (20.0 mg·m⁻³). Surface oceanographic conditions showed a marked annual cycle with warmer conditions (30–31 °C) from July to October and colder ones (23–24 °C) from December to April; during the cold period, seasonal wind stress stimulated coastal upwelling. The EOF showed that 70% of Chlo concentration variability was controlled by the semi-annual pattern, which responded to coastal upwelling conditions during cold period, while rivers outflows influenced high Chlo concentration during the warm period. These results highlighted the importance for land-ocean interface to sustain coastal ecosystems' biological production and the major role of watersheds as sources of nutrients to maintain high biological production during warm periods in transitional tropical-subtropical zones.     

星载硬质吸波材料真空功率耐受性能验证

星载吸波材料是复杂空间环境条件下天线、微波部组件大功率使用中满足隔离度要求的核心部件。本文研究了星载吸波材料真空功率耐受性能和吸波材料原材料制备工艺关键要素之间的关系,文中首先介绍了星载吸波材料电磁波吸收机理,其电磁参数直接关系到电磁波吸收性能优异与否;接着给出了影响电磁参数稳定性的原材料制备关键要素,以及成型材料的机加工艺特点;随后构建了一套星载吸波材料功率耐受性能的验证平台,开展功率试验;对两种工艺固化方法制备的原材料,分别制作了波导型吸收负载试验件进行试验验证。结果表明：（1）未进行高温预处理的吸波材料,残存未固化的小分子,会导致在高温真空工况下可凝挥发物析出增多,与外导体镀银层发生氧化反应,进而使负载组件的驻波变化率较大;（2）经过高温预处理后的吸波材料,在高温下的真空质损和可凝挥发物均得到了有效控制,其电磁参数也趋于稳定,负载组件的驻波变化率试验前后差异不大。因此,吸波材料原材料工艺制备过程中高温预处理属于关键要素;该工艺固化方法的有效实施将有助于星载吸波材料的应用,提高航天器在轨服役可靠性和安全性。     

ADS数据的处理和显示

基于航空卫星通信的飞行监视,即自动相关监视（ADS）,是航空卫星移动业务（AMSS）通信应用新领域,是建立在航空卫星通信基础上的一种新的飞行监视技术。该文介绍了基本ADS报告信息电文的组成及各部分的含义,着重阐述将所收到的ADS数据进行处理,并以直观易懂的表格及图象形式显示出来,便于确切地知道飞机的四维位置（纬度、经度、高度和时间）。     

卫星沿轨道运动对地面覆盖带外沿轨迹的确定

根据单参数平面曲线族包络线求解原理,给出了卫星沿轨道运动对地面覆盖带外沿轨迹的解,井讨论了地球旋转对地面覆盖带的影响。