Attitude analysis of space debris using SLR and light curve data measured with single-photon detector

Attitude is the important parameter for active debris removal and collision avoidance. This paper deduced the spin axis orientation and spin period of the rocket body, CZ-3B R/B (NORAD ID 38253), using the satellite laser ranging and light curve data measured with single-photon detector at Graz station. The epoch method and LC & SLR residuals fitting were combined to determine these values. The derived right ascension angle was around 220°, the declination angle was near 64° and the sidereal period was calculated to be 117.724 s, for 2017-07-03. The results derived from the two distinct methods were mutually validated. Rocket bodies are a major contributor to space debris and this work provides a reference for attitude determination and attitude modelling.     

Spin axis orientation of Ajisai determined from Graz 2 kHz SLR data

The Graz 2 kHz Satellite Laser Ranging (SLR) measurements allow determination of the spin axis orientation of the geodetic satellite Ajisai. The high repetition rate of the laser makes it possible to determine the epoch time when the laser is pointing directly between two corner cube reflector (CCR) rings of the satellite. Identification of many such events during a few (up to 3) consecutive passes allows to state the satellite orientation in the celestial coordinate system. Six years of 2 kHz SLR data (October 2003–October 2009) delivered 331 orientation values which clearly show precession of the axis along a cone centered at 14^h56^m2.8^s in right ascension and 88.512° in declination (J2000.0 celestial reference frame) and with an half-aperture angle θ of 1.405°. The spin axis precesses with a period of 117 days, which is equal to the period of the right ascension of the ascending node of Ajisai’s orbit. We present a model of the axis precession which allows prediction of the satellite orientation – necessary for the envisaged laser time transfer via Ajisai mirrors.     

Intercosmos-19 observations of an additional topside ionization layer: the F3 layer

Spatial properties of an additional ionization layer in the topside ionosphere were investigated using Intercosmos-19 satellite ionograms. The data under analysis were choosen for equinoctial conditions of the high solar activity period (1979–1981). The F3 layer was detected in a narrow longitude sectors (about 60°) between the equatorial anomaly crests. Its intensity has a maximum just above the equator and decreases poleward within ±10° dip. A nighttime F3 layer was observed as well as the daytime events.     

A payload for the study of electric fields and electron density in the equatorial region

A rocket borne payload for simultaneous measurement of the electric field along and perpendicular to the rocket spin axis and the electron density in the medium was developed and flown from Thumba (8° 31′N, 0° 47′S dip) onboard two Centaure rockets for the study of plasma dynamcis in the equatorial E-region. The arrangement of sensors in this payload allows near continuous measurements of some of these parameters to be made.     

Spin parameters of nanosatellite BLITS determined from Graz 2 kHz SLR data

The nanosatellite BLITS (Ball Lens In The Space) is the first object designed as a passive, spherical retroreflector of the Luneburg type, dedicated for Satellite Laser Ranging (SLR). The 2 kHz SLR station Graz measures spin parameters of this satellite, providing information about the rotational dynamics of the body. The measurements obtained during the period from September 26, 2009 to November 24, 2010 show a significant change of the spin configuration. The spin axis was dynamically precessing since the launch and currently is sinus-like behaving between coordinates RA 120°…150°, Dec 30°…60° (J2000 inertial reference frame). The angle between the symmetry axis and the spin axis of BLITS is not constant, but is decreasing since the launch, while its spin period is rather stable with a mean value of 5.613 s (clockwise rotation). The satellite was dynamically changing its attitude during the first three months after deployment; after this time the spin parameters are relatively stable.     

基于多种观测手段的东亚低纬电离层不规则体事件分析

利用海南台站（19.5°N,109.1°E,dip:13.6°N）和磁赤道区的多种地基和天基观测数据,对2011年11月20日观测到的电离层不规则体事件进行了分析.海南台站VHF雷达、电离层闪烁和数字测高仪的综合观测结果表明,当天日落附近发生了强的电离层不规则体事件,主要表现为雷达羽和强闪烁的形态.结合磁赤道区GPS和C/NOFS卫星观测结果进行分析可知,海南台站日落附近出现的雷达羽和强闪烁与南海磁赤道区产生的主等离子体泡存在明显联系.      

Characterization of GPS-TEC over African equatorial ionization anomaly (EIA) region during 2009–2016

This study characterizes total electron content (TEC) measured by Global Positioning System (GPS) over African equatorial ionization anomaly (EIA) region for 2009–2016 period during both quiet geomagnetic conditions (Kp?≤?1) and normal conditions (1?>?Kp?≤?4). GPS-TEC data from four equatorial/low-latitude stations, namely, Addis Ababa (ADIS: 9.04°N, 38.77°E, mag. lat: 0.2°N) [Ethiopia]; Yamoussoukro (YKRO: 6.87°N, 5.24°W, mag. lat: 2.6°S) [Ivory Coast]; Malindi (MAL2; 3.00°S, 40.19°E, mag. lat: 12.4°S) [Kenya] and Libreville (NKLG; 0.35°N, 9.67°W, mag. lat: 13.5°S) [Gabon] were used for this study. Interesting features like noontime TEC bite-out, winter anomaly during the ascending and maximum phases of solar cycle 24, diurnal and seasonal variations with solar activity have been observed and investigated in this study. The day-to-day variations exhibited ionospheric TEC asymmetry on an annual scale. TEC observed at equatorial stations (EIA-trough) and EIA-crest reach maximum values between ～1300–1600 LT and ～1300–1600 LT, respectively. About 76% of the high TEC values were recorded in equinoctial months while the June solstice predominantly exhibited low TEC values. Yearly, the estimated TEC values increases or decreases with solar activity, with 2014 having the highest TEC value. Solar activity dependence of TEC within the EIA zone reveals that both F10.7?cm index and EUV flux (24–36?nm) gives a stronger correlation with TEC than Sunspot Number (SSN). A slightly higher degree of dependence is on EUV flux with the mean highest correlation coefficient (R) value of 0.70, 0.83, 0.82 and 0.88 for quiet geomagnetic conditions (Kp?≤?1) at stations ADIS, MAL2, NKLG, and YKRO, respectively. The correlation results for the entire period consequently reveals that SSN and solar flux F10.7?cm index might not be an ideal index as a proxy for EUV flux as well as to measure the variability of TEC strength within the EIA zone. The estimated TEC along the EIA crest (MAL2 and NKLG) exhibited double-hump maximum, as well as post-sunset peaks (night time enhancement of TEC) between ～2100 and 2300 LT. EIA formation was prominent during evening/post-noon hours.     

Spin rate and spin axis orientation of LARES spectrally determined from Satellite Laser Ranging data

Satellite Laser Ranging (SLR) is a powerful technique able to measure spin rate and spin axis orientation of the fully passive, geodetic satellites. This work presents results of the spin determination of LARES – a new satellite for testing General Relativity. 529 SLR passes measured between February 17 and June 9, 2012, were spectrally analyzed. Our results indicate that the initial spin frequency of LARES is f₀ = 86.906 mHz (RMS = 0.539 mHz). A new method for spin axis determination, developed for this analysis, gives orientation of the axis at RA = 12^h22^m48^s (RMS = 49^m), Dec = −70.4° (RMS = 5.2°) (J2000.0 celestial reference frame), and the clockwise (CW) spin direction. The half-life period of the satellite’s spin is 214.924 days and indicates fast slowing down of the spacecraft.     

“嫦娥4号”中继星使命轨道段定轨计算与分析

"嫦娥4号"中继星是"嫦娥4号"探测器实现月球背面着陆与巡视的关键,目前正稳定运行在地-月L2点使命轨道上,该使命轨道为平均周期约14天的南族Halo轨道。因任务的需要,中继星本体系+Z轴需调整指向,处于正对太阳和非正对太阳两种状态。太阳光压在中继星+Z轴对日的情况下会加速卫星的角动量累积,增加卫星卸载喷气频次。基于中继星使命轨道段测控支持条件,采用重叠弧段法对两种状态下的中继星定轨精度进行分析与评估。结果表明,在中继星+Z轴非对日运行状态下,重叠弧段位置误差为1.6 km,速度误差为8 mm/s;在中继星+Z轴对日运行状态下,重叠弧段位置误差为0.6 km,速度误差为3 mm/s,这对中继星的长期运行具有重要参考价值。     

The HESP/R satellite project

The aim of the HESP/R (High Energy Solar Physics/Radiation) satellite project is to obtain data of γ-ray, hard X-ray soft X-ray, EUV, and visible radiation of solar flares at the next solar maximum in order to study physics of flares. The HESP/R will be a spinning satellite of 4–5 rotations per min., and the spin axis is off-set by a small angle (0.5°–1.0°) from the Sun. Total weight will be 400 kg and launched in 1991 with M-3S-II rocket by ISAS.     

Global ionospheric scintillations revealed by GPS radio occultation data with FY3C satellite before midnight during the March 2015 storm

Global observations of S4 amplitude scintillation index by the GPS Occultation Sounder (GNOS) on FengYun-3 C (FY3C) satellite reveal global dynamic patterns of a strong pre-midnight scintillations in F-region of the ionosphere during the St. Patrick’s Day geomagnetic super storm of 17–19 March 2015. The observed strong scintillations mainly occurred in the low latitudes, caused by equatorial plasma bubbles. During the main storm phase (March 17), the scintillations were first triggered in the New Zealand sector near 160°E longitudes, extending beyond 40°S dip latitude. They were also enhanced in the Indian sector, but significantly suppressed in East Asia near 120°E longitude and in Africa around 30°E longitude. During the initial recovery phase (March 18–19), the global scintillations were seldom observed in GNOS data. During the later recovery phase (after March 19), the scintillations recovered to the pre-storm level in Indian, African, and American sectors, but not in East Asian and any of Pacific sectors. These results closely correlate with observations of the density depletion structures by the Communication/Navigation Outage Forecasting System (C/NOFS) satellite, and ground-based instruments. Such consistency indicates reliability of our scintillation sensing approach even in a case-by-case comparison study. The prompt penetration electric field and disturbance dynamo electric field are suggested as the main factors that control the enhancement and inhibition of the scintillations during the storm, respectively.     

利用近地点磁场探测数据确定卫星自旋轴参数

阐述了利用近地点磁场探测数据确定卫星自旋轴参数的理论方法和实施步骤,并说明了这种研究对卫星运行和科学探测的重要性.特别强调了需要注意的基本条件,即卫星必须自旋稳定且近地点不很高(1000 km以下).这种方法关键的步骤是,根据卫星轨道数据定出模型磁场数值,比较近地点星载磁强计探测数据和近地点地磁模型数值确定卫星自旋轴的指向.通过对TC-1和TC-2卫星姿态的具体计算,对确定精度和应用效果进行了分析和比较.结果表明,在实际的卫星应用过程中此方法和措施非常有效,在科学分析和将来的卫星运行工程中具有重要的应用意义.      

一种基于连续星图观测的自旋姿态估计方法

传统的利用地球敏感器和太阳敏感器作为测量仪器的自旋卫星姿态确定方法存在系统误差和安装误差等,从而导致自旋姿态确定误差较大的问题,文章提出了一种利用星敏感器获取的连续星图估计卫星自旋姿态参数的新方法。该方法以卫星的自旋轴和旋转角速度作为状态变量,通过星敏感器连续跟踪拍摄的恒星的成像位置作为观测量,利用无迹卡尔曼滤波估计出卫星的自旋姿态参数。仿真结果表明,在星敏感器的精度为3″时,该方法的自旋轴估计精度为0.3448″,自旋角速度估计精度为10^-4(°)/s数量级。     

TC-2卫星探测伸杆在技术状态更改设计中的仿真分析

讨论了TC-2卫星探测伸杆在技术状态更改中的结构和力学问题, 并对TC-2卫星 探测伸杆锁紧支座进行了模态分析、结构静强度计算和解锁机构展开的仿真分 析. 结果表明, 伸杆锁紧支座在更换火工切割器后, 支座频率变化不大, 不会 对整星的基频产生影响, 改进后的锁紧支座结构能够承受动载荷, 最大应力和 应变均在材料的许用范围之内, 锁紧支座在打开过程中不会与任何零件发生干 涉, 满足卫星对伸杆部件的要求.      

Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.     

Characterization of the low latitude plasma density irregularities observed using C/NOFS and SCINDA data

Complex electrodynamic processes over the low latitude region often result in post sunset plasma density irregularities which degrade satellite communication and navigation. In order to forecast the density irregularities, their occurrence time, duration and location need to be quantified. Data from the Communication/Navigation Outage Forecasting System (C/NOFS) satellite was used to characterize the low latitude ion density irregularities from 2011 to 2013. This was supported by ground based data from the SCIntillation Network Decision Aid (SCINDA) receivers at Makerere (Geographic coordinate 32.6°E, 0.3°N, and dip latitude ?9.3°N) and Nairobi (Geographic coordinate 36.8°E, ?1.3°N, and dip latitude ?10.8°N). The results show that irregularities in ion density have a daily pattern with peaks from 20:00 to 24:00 Local Time (LT). Scintillation activity at L band and VHF over East Africa peaked in 2011 and 2012 from 20:00 to 24:00 LT, though in many cases scintillation at VHF persisted longer than that at L band. A longitudinal pattern in ion density irregularity occurrence was observed with peaks over 135–180°E and 270–300°E. The likelihood of ion density irregularity occurrence decreased with increasing altitude. Analysis of C/NOFS zonal ion drift velocities showed that the largest nighttime and daytime drifts were in 270–300°E and 300–330°E longitude regions respectively. Zonal irregularity drift velocities over East Africa were for the first time estimated from L-band scintillation indices. The results show that the velocity of plasma density irregularities in 2011 and 2012 varied daily, and hourly in the range of 50–150 m s^?1. The zonal drift velocity estimates from the L-band scintillation indices had good positive correlation with the zonal drift velocities derived from VHF receivers by the spaced receiver technique.     

重力梯度力矩作用下近地卫星自旋运动规律分析

为研究近地卫星自旋运动规律,建立了近地卫星在受摄动影响的轨道上运行并受重力梯度力矩作用下的姿态运动模型,推导了自旋角速率满足一定条件下自旋运动的进动角、章动角、自旋角的解析解,对重力梯度作用下的自旋姿态运动规律进行了仿真分析,并用仿真计算结果验证了解析解的正确性。在轨道面缓慢进动情况下,当卫星绕最大主惯量轴自旋时,给出了自旋角速率取值范围表达式,在该取值范围内卫星自旋运动能够跟随轨道面一起进动,自旋轴以恒定的平均角速率进动,章动角在小范围内波动。建立的自旋姿态运动模型和分析结论可用于近地卫星姿态失控后的姿态确定和预测、在轨姿态设计及在轨备份等。     

A comparison of seasonal variations of sea level in the southern Baltic Sea from altimetry and tide gauge data

In this paper, seasonal sea level variations have been determined at five locations in the Baltic Sea from satellite altimetry for the period 1993–2015. The results were compared to tide gauge water level data. Annual and semi-annual amplitudes have been investigated for both sea level anomalies and tide gauge water level. It was found that the two independent observations of sea level variations along the Polish coast are in good agreement both in terms of their annual and semi-annual amplitudes and their annual and semi-annual phases. The annual cycles in the sea level variations measured by altimetry and tide gauge reach maximum values at approximately the same month (November/December).Moreover, this article shows the differences between the annual and semi-annual amplitudes and phases in the sea level anomalies and water level data within the same time frame. The difference in the annual amplitudes between the satellite altimetry and the tide gauge results is between 0.33?cm and 1.53?cm. The maximum differences in the annual cycle of the sea level changes were found at the Swinoujscie station. The correlations between the original series and the calculated curves were determined, and the relationship between the amplitudes and the phases were investigated. The correlation between the annual variations observed from the two independent observation techniques is 0.92.To analyse the dynamics of the change in sea level, the linear trend was estimated from the satellite altimetry and tide gauge time series both in the original time series of the data and in the time series in which seasonal variations were removed. In addition, we calculated the estimated errors of regression and how many years’ worth of data are needed to obtain an accuracy of 0.1?mm per year. The estimated errors of regression showed that to get an accuracy of 0.1?mm per year, we need 100?years of data.     

Spectral filter for signal identification in the kHz SLR measurements of the fast spinning satellite Ajisai

The high repetition rate satellite laser ranging (SLR) measurements to the fast spinning satellites contain a frequency signal caused by the rotational motion of the corner cube reflector (CCR) array. The spectral filter, developed here, is based on the Lomb algorithm, and is tested with the simulated and the observed high repetition rate SLR data of the geodetic satellite Ajisai (spin period ∼2 s). The filter allows for the noise elimination from the SLR data, and for identification of the returns from the single CCRs of the array – even for the low return rates. Applying the spectral filter to the simulated SLR data increases the S/N ratio by a factor 40–45% for all return rates. Filtering out the noise from the observed data strengthens the frequency signal by factor of ∼25 for the low return rates, which significantly helps to determine the spin phase of the satellite. The spectral filter is applied to the Graz SLR data and the spin rates of Ajisai are determined by two different methods: the frequency analysis and the phase determination of the spinning retroreflector array.The analysis of more than 8 years of the Graz SLR measurements indicates an exponential spin rate trend: f = 0.67034 exp(−0.0148542 Y) [Hz], RMS = 0.085 mHz, where Y is the year since launch. The highly accurate spin rate information demonstrates periodic changes related to the precession of the orbital plane of Ajisai, as it determines the amount of energy received by the satellite from the Sun. The rate of deceleration of Ajisai is not constant: the half life period of the satellite’s spin oscillates around 46.7 years with an amplitude of about 5 years.