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.     

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

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

Ionospheric response to the storm-time disturbance of 29 May, 2010

The Ionospheric F2-layer peak parameters response to a magnetic storm had been investigated over Ilorin, Nigeria (Lat. 8:53°N, Long. 4.5°E, dip angle, −2.96°), Jicamarca, Peru (11.95°S, 76.87°W, dip angle, 0.8°) and Hermanus, South Africa (34.42°S, 19.22°E, dip angle, −60.77°), using percentage enhancement/depletion values. Our results showed an enhancement in NmF2 at all of these stations. Averagely, pre-noon and post-noon peaks are highest at Ilorin during quiet time. The similar pattern observed for quiet condition between Ilorin and Jicamarca was due to their latitudinal positions. For disturbed NmF2 condition, Jicamarca and Ilorin recorded higher peaks at nighttime than during the daytime for the storms main phase, and the reverse over Hermanus. The nighttime and daytime increases were observed respectively at Ilorin and Hermanus during the recovery period. The hmF2 variation recorded higher enhancement at Jicamarca during the daytime and at Hermanus at nighttime during the main phase. During the recovery phase, the highest enhancement was recorded during the daytime at Jicamarca, and over Hermanus at nighttime. These observations find their explanation in the magnetospheric current, solar wind and E × B drift.     

Estimating the spin axis orientation of the Echostar-2 box-wing geosynchronous satellite

For the first time, the spin axis orientation of an inactive box-wing geosynchronous satellite has been estimated from ground-based optical photometric observations of Echostar-2’s specular reflections. Recent photometric light curves obtained of Echostar-2 over four years suggest that unusually bright and brief specular reflections were occurring twice within an observed spin period. These bright and brief specular reflections suggested two satellite surfaces with surface normals separated by approximately 180°. The geometry between the satellite, the Sun, and the observing location at the time of each of the brightest observed reflections, was used to estimate Echostar-2’s equatorial spin axis orientation coordinates. When considering prograde and retrograde rotation, Echostar-2’s spin axis orientation was estimated to have been located within 30° of either equatorial coordinate pole. Echostar-2’s spin axis was observed to have moved approximately 180° in right ascension, within a time span of six months, suggesting a roughly one year spin axis precession period about the satellite’s angular momentum vector.     

非合作失效航天器绳系拖曳系统研究

空间碎片的不断增加给人类航天活动的开展和在轨资产的安全造成严重威胁。在已经提出的多种空间碎片主动清除方式中,绳系拖曳（tethered space tug,TST）系统因有较好的应用前景而受到广泛关注。部分失效航天器因星上器件损坏且姿态控制系统异常,始终将姿态维持在某一特定指向,针对此类具有典型的非合作特征的大型空间碎片,开展绳系拖曳动力学与控制研究。将拖船和目标均视作刚体,用牛顿法建立了TST系统的动力学模型;根据目标姿态稳定方式分为自旋稳定和三轴稳定两种情况,开展了绳系拖曳动力学分析与控制设计,并考察了系绳在失效航天器表面连接点位置对系统稳定状态的影响。仿真结果表明,拖船和失效航天器会在系绳连接下表现出抗衡特点,最终系统会稳定在不同的均衡状态附近。此研究为安全清除尚有残余姿态控制能力的失效航天器中相关拖曳动力学与控制问题提供了参考。     

北极69°N和78°N空间碎片凝视探测的对比分析

利用北极69°N和78°N两套非相干散射雷达的首次空间碎片联合观测数据进行空间碎片参数（距离、速度、散射截面积、等效直径等）的对比分析,得出以下结论:两部雷达探测的碎片高度均主要分布在500～1100km和1400～1600km区间,但78°N雷达探测的碎片数量较多;空间碎片的径向速度均在-1.5～1.5km…-1区间,其中大部分为负值,说明在此次探测试验中碎片运动方向主要以远离雷达或地球为主;ESR雷达探测的空间碎片射截面积约为10-5～10-2m2,等效直径主要分布在4～10cm,而UHF雷达探测的空间碎片散射截面积约为10-6～10-2m2,等效直径主要分布在2～6cm,说明在同一高度上69°N雷达探测能力更强;经合理设置判据参数后得出重复检测次数,78°N雷达和69°N雷达分别有32次和14次重复检测,两部雷达共有4次重复检测.这些结果为空间碎片检测和建模提供了参考.      

A detailed spin analysis of the San Marco 5 mission

The attitude of the San Marco 5 satellite flown in 1988 has been monitored by several sensors. Thus the history of the spin period is known with a high degree of accuracy. Because of the simple geometry (spherical body with several extremely long antennas) and the good mass balance (used for the accelerometer aboard) of the satellite it was possible to separate different effects of variations of the spin period. The influence of the drag on the spin period has been modeled. The results are used to derive total gas densities along the trajectory and compare them with model densities, in situ measurements from the drag balance instrument and densities derived by orbital drag methods.     

Yarkovsky-Schach effect on space debris motion

The Yarkovsky-Schach effect is a small perturbation affecting Earth satellites and space debris illuminated by the Sun. It was first applied to the orbit of LAGEOS satellites as an explanation of the residuals in orbital elements. In this work, we carry out several numerical integration tests taking into consideration various orbit and rotation parameters, in order to analyse this effect in a broader context. The semi-major axis variations remain small and depend on the spin axis attitude with respect to the Sun. We show that the force amplitude is maximised for orbits inclined with $i\approx$?20–30°. We also observe the influence on other orbital elements, notably on the orbit inclination. However, these effects are clearly observed only on long timescales; in our simulations, we propagated the orbits for 200?y. The Yarkovsky-Schach effect is thus confirmed to have a minuscule magnitude. It should be taken into account in studies requiring high-precision orbit determination, or on expanded timescales.     

On the relationship between clear-sky planetary and surface albedos: a parameterization for simple energy balance climate models

In simple energy balance climate models all physical and dynamical processes are parameterized in terms of the single unknown variable: the surface temperature. To simulate the ice-albedo feedback, the surface albedo is usually assumed to be a function of surface temperature. But to compute the absorbed solar radiation in such models one requires the top-of-the atmosphere albedo: the planetary albedo. In the present study, a simple linear relationship is derived between planetary albedo and surface albedo for the case of clear skies. The relationship is based upon a regression equation derived from simulations and has a standard error of estimate of 0.028. The estimation of planetary albedo from surface albedo is checked by comparing zonally averaged clear-sky planetary albedos, estimated from zonally averaged surface albedos, to satellite determinations of zonally averaged minimum albedos for monthly mean conditions. The minimum albedos are assumed to be representative of the clear-sky planetary albedos. The results show root-mean square differences of 0.05 between the estimated clear-sky planetary albedos and the minimum albedos.More accurate relationships can be obtained if one uses an additional parameter - the solar zenith angle. In this case, the standard errors of estimate are reduced to 0.017 for a zenith angle of 0°, 0.018 for a zenith angle of 60° and 0.021 for a zenith angle of 85°.     

Simulation of long-term rotational dynamics of large space debris: A TOPEX/Poseidon case study

Accurate knowledge of the rotational dynamics of a large space debris is crucial for space situational awareness (SSA), whether it be for accurate orbital predictions needed for satellite conjunction analyses or for the success of an eventual active debris removal mission charged with stabilization, capture and removal of debris from orbit. In this light, the attitude dynamics of an inoperative satellite of great interest to the space debris community, the joint French and American spacecraft TOPEX/Poseidon, is explored. A comparison of simulation results with observations obtained from high-frequency satellite range measurements is made, showing that the spacecraft is currently spinning about its minor principal axis in a stable manner. Predictions of the evolution of its attitude motion to 2030 are presented, emphasizing the uncertainty on those estimates due to internal energy dissipation, which could cause a change of its spin state in the future. The effect of solar radiation pressure and the eddy-current torque are investigated in detail, and insights into some of the satellite’s missing properties are provided. These results are obtained using a novel, open-source, coupled orbit-attitude propagation software, the Debris SPin/Orbit Simulation Environment (D-SPOSE), whose primary goal is the study of the long-term evolution of the attitude dynamics of large space debris.     

Assessing the performance of BDS-3 for multi-GNSS static and kinematic PPP-AR

As of 2021, a total of four different GNSS constellations – namely, GPS, GLONASS, Galileo, and BDS-3 – can be used with Full Operational Capability (FOC). In this work, the contribution of BDS-3 FOC to GPS + GLONASS + Galileo (GRE) PPP-AR is investigated, considering the three different cut-off angles (7°, 30°, and 45°) and different lengths of static observation sessions (24-, 12-, 6-, 3-, 1-, 0.5-, 0.25-hour). The data of 31 IGS-MGEX stations is processed with GRE PPP-AR and GREC3 (GPS + GLONASS (using float mode) + Galileo + BDS-3) PPP-AR modes. The results showed that BDS-3 degraded the horizontal (except for 24-h sessions) and vertical accuracy of static GRE PPP-AR solutions regardless of the elevation cutoff angle and observation time. The kinematic results showed that BDS-3 significantly contributed to the accuracy of GRE kinematic PPP-AR for 30° and 45° cutoff angles. The convergence time analysis showed that BDS-3 only contributes to GRE kinematic PPP-AR for the vertical component.     

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.     

22 Years of AJISAI spin period determination from standard SLR and kHz SLR data

Satellite Laser Ranging (SLR) is a powerful and efficient technique to measure spin parameters of satellites equipped with corner cube reflectors. We obtained spin period determination of the satellite AJISAI from SLR data only: 17246 pass-by-pass estimates from standard 1–15 Hz SLR data (14/Aug/1986–30/Dec/2008) and 1444 pass-by-pass estimates (9/Oct/2003–30/Dec/2008) from data of the first 2 kHz SLR system from Graz, Austria. A continuous history of the slowing down of AJISAI spin is derived from frequency analysis, and corrected for the apparent effects. The apparent corrections, elaborated here, allowed very accurate determination of AJISAI initial spin period: 1.4855 ± 0.0007 [s]. The paper identifies also non-gravitational effects as a source of the periodical changes in the rate of slowing down of the satellite.     

Deorbiter CubeSat mission design

This paper presents the mission design for a CubeSat-based active debris removal approach intended for transferring sizable debris objects from low-Earth orbit to a deorbit altitude of 100 km. The mission consists of a mothership spacecraft that carries and deploys several debris-removing nanosatellites, called Deorbiter CubeSats. Each Deorbiter is designed based on the utilization of an eight-unit CubeSat form factor and commercially-available components with significant flight heritage. The mothership spacecraft delivers Deorbiter CubeSats to the vicinity of a predetermined target debris, through performing a long-range rendezvous maneuver. Through a formation flying maneuver, the mothership then performs in-situ measurements of debris shape and orbital state. Upon release from the mothership, each Deorbiter CubeSat proceeds to performing a rendezvous and attachment maneuver with a debris object. Once attached to the debris, the CubeSat performs a detumbling maneuver, by which the residual angular momentum of the CubeSat-debris system is dumped using Deorbiter’s onboard reaction wheels. After stabilizing the attitude motion of the combined Deorbiter-debris system, the CubeSat proceeds to performing a deorbiting maneuver, i.e., reducing system’s altitude so much so that the bodies disintegrate and burn up due to atmospheric drag, typically at around 100 km above the Earth surface. The attitude and orbital maneuvers that are planned for the mission are described, both for the mothership and Deorbiter CubeSat. The performance of each spacecraft during their operations is investigated, using the actual performance specifications of the onboard components. The viability of the proposed debris removal approach is discussed in light of the results.     

Velocity pointing error analysis for symmetric spinning thrusting Cubesat

To solve the problem of thrust vector misalignment from the Cubesat center of mass during orbital maneuver, spin-stabilized method is applied to eliminate velocity pointing error. Spinning thrusting Cubesat model involves the effects of mass variation and jet damping is established. Analytical solutions for the angular velocity, nutation angle, Euler angle, and inertial velocity with nonzero initial conditions are derived. Simulations show that the analytical solutions closely match numerical simulations. Based on the analytical solutions, the velocity pointing error influencing factors is analyzed. The results show that the velocity pointing error caused by initial transverse velocity, nutation angle and transverse disturbance torque can be reduced by raising the spin rate, but the initial Euler angle need to be limited. Also, the spinning thrusting maneuver can allow for a lower spin rate by increasing the axial moment of inertia.     

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.     

A deformation model of flexible,HAMR objects for accurate propagation under perturbations and the self-shadowing effects

A new type of space debris in near geosynchronous orbit (GEO) was recently discovered and later identified as exhibiting unique characteristics associated with high area-to-mass ratio (HAMR) objects, such as high rotation rates and high reflection properties. Observations have shown that this debris type is very sensitive to environmental disturbances, particularly solar radiation pressure, due to the fact that its motion depends on the actual effective area, orientation of that effective area, reflection properties and the area-to-mass ratio of the object is not stable over time. Previous investigations have modelled this type of debris as rigid bodies (constant area-to-mass ratios) or discrete deformed body; however, these simplifications will lead to inaccurate long term orbital predictions. This paper proposes a simple yet reliable model of a thin, deformable membrane based on multibody dynamics. The membrane is modelled as a series of flat plates, connected through joints, representing the flexibility of the membrane itself. The mass of the membrane, albeit low, is taken into account through lump masses at the joints. The attitude and orbital motion of this flexible membrane model is then propagated near GEO to predict its orbital evolution under the perturbations of solar radiation pressure, Earth’s gravity field (J₂), third body gravitational fields (the Sun and Moon) and self-shadowing. These results are then compared to those obtained for two rigid body models (cannonball and flat rigid plate). In addition, Monte Carlo simulations of the flexible model by varying initial attitude and deformation angle (different shape) are investigated and compared with the two rigid models (cannonball and flat rigid plate) over a period of 100?days. The numerical results demonstrate that cannonball and rigid flat plate are not appropriate to capture the true dynamical evolution of these objects, at the cost of increased computational time.     

ETALON spin period determination from kHz SLR data

Using kHz Satellite LASER Ranging (SLR) data of the SLR station Graz only, we determined the spin periods of the two ETALON satellites – launched into high orbits of about 20,000 km – and their spin period increase during 3 years. The determined spin period values and spin period increase rates at 2004-01-01 are: T_ET1 = 63 s + 0.484 s/year, and T_ET2 = 65.5 s + 0.401 s/year.     

Modeling of effect of polarization on UV sky radiance during twilight

The spatial distribution of the vector of the Stokes parameters characterizing the radiance intensity and the radiance polarization describes the radiation field in the atmosphere. A simplified treatment of light as the scalar has only restricted application. A few studies compared previously results of the vector and scalar radiative transfer models and showed that scalar models are in error by up to 10% for many cases. Though several observational conditions were exploited, an effect of polarization on modeling of UV radiance has not been investigated yet for twilight. The paper presents a preliminary study of modeled UV radiance during twilight taking into account polarization. The intensity and the degree of linear polarization of the scattered UV radiance for two cases of the ground-based observations are discussed. In the first case, radiation incoming from the zenith for the solar zenith angles (SZA) from 90° to 98° is under investigation. Radiation in the solar principal plane for the beginning of twilight (SZA = 90.1°) was calculated in the second case. The study showed that the UV radiation field in the twilight atmosphere can be handled correctly only using the vector theory. The errors of scalar radiative transfer strongly depend on wavelength, line of an observation and solar position. The revealed distortion of the zenith radiance caused by using of the scalar approximation reaches maximum of 15% at 340 nm for the solar zenith angle (SZA) equal to 98°. The shorter wavelengths have the smaller errors, about 5% at 305 nm for SZA = 98°, due to the larger part of the single scattered radiance. The error of the scalar modeling may be as large as −17% for radiance incoming from the horizon for SZA = 90.1°. Scalar radiative transfer models underestimate the integral intensity in the principal plane up to 3–4% ± 0.5% at SZA = 90.1° for wavelengths from 320 to 340 nm. This should be taken into account in problems of radiative budget estimation and remote sensing of the atmosphere exploiting the twilight period.     

航天器自主交会对接的视觉相对导航方法

针对航天器自主交会对接在最终逼近阶段的相对导航问题,研究了基于特征光标的航天器视觉相对导航方法。首先建立了小角度假设条件下的相对导航模型。给出了航天器相对位置和采用欧拉角描述的相对姿态的近似解析算法。其次利用相对姿态的解析解作为观测信息,采用降阶的无迹卡尔曼滤波方法对目标航天器的姿态角和惯性角速度进行了估计,间接得到了相对姿态和相对姿态角速度的估计值。数值仿真表明提出的解析算法能够有效获得相对位置和相对姿态测量信息,通过降阶滤波可明显提高相对姿态确定精度。