Non-conservative forces on LAGEOS I and II

The behavior of the empirically determined along-track accelerations from the LAGEOS satellites has been attributed to a variety of physical phenomena. Of these, the models for Yarkovsky thermal drag, anisotropic reflectivity and Yarkovsky-Schach drag are dependent on the spin axis orientation of the satellite. This investigation explores the utilization of these models in an attempt to recover the spin axis orientation history, particularly for the LAGEOS I satellite.     

Drag and energy accommodation coefficients during sunspot maximum

Conditions appropriate to gas-surface interactions on satellite surfaces in orbit have not been successfully duplicated in the laboratory. However, measurements by pressure gauges and mass spectrometers in orbit have revealed enough of the basic physical chemistry that realistic theoretical models of the gas-surface interaction can now be used to calculate physical drag coefficients. The dependence of these drag coefficients on conditions in space can be inferred by comparing the physical drag coefficient of a satellite with a drag coefficient fitted to its observed orbital decay. This study takes advantage of recent data on spheres and attitude stabilized satellites to compare physical drag coefficients with the histories of the orbital decay of several satellites during the recent sunspot maximum. The orbital decay was obtained by fitting, in a least squares sense, the semi-major axis decay inferred from the historical two-line elements acquired by the US Space Surveillance Network. All the principal orbital perturbations were included, namely geopotential harmonics up to the 16th degree and order, third body attraction of the Moon and the Sun, direct solar radiation pressure (with eclipses), and aerodynamic drag, using the Jacchia-Bowman 2006 (JB2006) model to describe the atmospheric density. After adjusting for density model bias, a comparison of the fitted drag coefficient with the physical drag coefficient has yielded values for the energy accommodation coefficient as well as for the physical drag coefficient as a function of altitude during solar maximum conditions. The results are consistent with the altitude and solar cycle variation of atomic oxygen, which is known to be adsorbed on satellite surfaces, affecting both the energy accommodation and angular distribution of the reemitted molecules.     

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.     

Mesospheric ozone densities deduced from solar occultation measurements

Quartz-UV occultation measurements by the satellite Interkosmos-16 have been used to calculate ozone densities at altitudes between 50 and 90 km for the period August to October 1976. Below 65 km densities agree well with the Krueger-Minzner-model. Mesopause densities have been studied in some detail. A certain percentage of the profiles show close correlation with the model of Shimazaki and Laird (with a pronounced minimum below the mesopause) while others fit better to the Park and London model (no minimum). This variability of the ozone density may be caused by different processes in the photo-chemistry of ozone. Two possible causes, the temperature dependent rate coefficients and the odd hydrogen processes are discussed in greater detail.     

Nanosatellite constellation deployment using on-board magnetic torquer interaction with space plasma

One of the advantages that drive nanosatellite development is the potential of multi-point observation through constellation operation. However, constellation deployment of nanosatellites has been a challenge, as thruster operations for orbit maneuver were limited due to mass, volume, and power. Recently, a de-orbiting mechanism using magnetic torquer interaction with space plasma has been introduced, so-called plasma drag. As no additional hardware nor propellant is required, plasma drag has the potential in being used as constellation deployment method. In this research, a novel constellation deployment method using plasma drag is proposed. Orbit decay rate of the satellites in a constellation is controlled using plasma drag in order to achieve a desired phase angle and phase angle rate. A simplified 1D problem is formulated for an elementary analysis of the constellation deployment time. Numerical simulations are further performed for analytical analysis assessment and sensitivity analysis. Analytical analysis and numerical simulation results both agree that the constellation deployment time is proportional to the inverse square root of magnetic moment, the square root of desired phase angle and the square root of satellite mass. CubeSats ranging from 1 to 3?U (1–3?kg nanosatellites) are examined in order to investigate the feasibility of plasma drag constellation on nanosatellite systems. The feasibility analysis results show that plasma drag constellation is feasible on CubeSats, which open up the possibility of CubeSat constellation missions.     

任意偏心率卫星的大气阻力的短周期摄动

通过本文采用的计算方法和CIRA 1972在160～2000公里高度围范的大气密度的分析模型,导出了人造地球卫星大气阻力短周期摄动的一种计算公式,它形式简单,收敛性好,适用于任意偏心率的卫星。     

我国自旋稳定通信卫星四种发展方案的比较

以国外自旋稳定通信卫星发展的技术途径为背景,比较了在“东方红二号甲”卫星技术的基础上进一步发展中国自旋稳定通信卫星的不同技术方案;说明方案的选择主要取决于对今后继续发展自旋稳定通信卫星的战略考虑。     

Towards accurate atmospheric mass density determination Using precise positional information of space objects

This paper presents a new method of deriving atmospheric mass densities with a high temporal resolution from precise orbit data of low earth orbiting (LEO) space objects. This method is based on the drag perturbation equation of the semi-major axis of the orbit of LEO space objects which relates the change rate of the semi-major axis to the atmospheric mass density. The effectiveness of the new method is evaluated using the GFZ-ISDC GPS rapid science orbit (RSO) products of the CHAMP satellite over a time period of 3 months. The densities derived using this new method and obtained from accelerometer data are compared and good agreements are achieved. An example of using the derived density to generate good orbit prediction for CHAMP is presented.     

Thermospheric neutral density response to solar forcing

Recent measurements by the Solar EUV (Extreme Ultra Violet) Experiment (SEE) aboard the Thermosphere–Ionosphere–Mesosphere Energetics and Dynamics satellite (TIMED) provide solar EUV spectral irradiance with adequate spectral and temporal resolution, and thus the opportunity to use solar measurements directly in upper atmospheric general circulation models. Thermospheric neutral density is simulated with the NCAR Thermosphere–Ionosphere–Electrodynamic General Circulation Model (TIEGCM) using TIMED/SEE measurements and using the EUVAC solar proxy model. Neutral density is also calculated using the NRLMSISE-00 empirical model. These modeled densities are then compared to density measurements derived from satellite drag data. It is found that using measured solar irradiance in the general circulation model can improve density calculations compared to using the solar proxy model. It is also found that the general circulation model can improve upon the empirical model in simulating geomagnetic storm effects and the solar cycle variation of neutral density.     

Drag coefficient modeling for grace using Direct Simulation Monte Carlo

Drag coefficient is a major source of uncertainty in predicting the orbit of a satellite in low Earth orbit (LEO). Computational methods like the Test Particle Monte Carlo (TPMC) and Direct Simulation Monte Carlo (DSMC) are important tools in accurately computing physical drag coefficients. However, the methods are computationally expensive and cannot be employed real time. Therefore, modeling of the physical drag coefficient is required. This work presents a technique of developing parameterized drag coefficients models using the DSMC method. The technique is validated by developing a model for the Gravity Recovery and Climate Experiment (GRACE) satellite. Results show that drag coefficients computed using the developed model for GRACE agree to within 1% with those computed using DSMC.     

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.     

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.     

Thermosphere density and wind measurements in the equatorial region using a constellation of drag balance nanospacecraft

A mission for in situ thermosphere density and winds measurement is described, based on nanospacecraft equipped with a drag balance instrument (DBI) and a GPS receiver. The mission is based on nanosatellite clusters deployed in three orbital planes. In this study, clusters of 10 nanospacecraft are considered, leading to a mission based on a total of 30 nanospacecraft. The geometry analyzed is a symmetrical one, including an equatorial orbit and two orbits with the same inclination and opposing ascending nodes. The main idea is that, by combining the accurate information on the satellite inertial position and velocity provided by the GPS receiver and the drag acceleration intensity provided by the DBI, due to the orbits’ geometrical configuration, both atmospheric drag and wind can be resolved in a region close to the orbit nodes. Exploiting the Earth oblateness effect, a complete scan of the equatorial regions can be accomplished in the short mission lifetime typical of very low Earth orbit satellites, even in high solar activity peaks, when the expected nanospacecraft lifetime is about 40 days.     

Use of two-line element data for thermosphere neutral density model calibration

Traditional empirical thermospheric density models are widely used in orbit determination and prediction of low-Earth satellites. Unfortunately, these models often exhibit large density errors of up to around 30% RMS. Density errors translate into orbit errors, adversely affecting applications such as re-entry operations, manoeuvre planning, collision avoidance and precise orbit determination for geodetic missions. The extensive database of two-line element (TLE) orbit data contains a wealth of information on satellite drag, at a sufficiently high spatial and temporal resolution to allow a calibration of existing neutral density models with a latency of one to two days. In our calibration software, new TLE data for selected objects is converted to satellite drag data on a daily basis. The resulting drag data is then used in a daily adjustment of density model calibration parameters, which modify the output of an existing empirical density model with the aim of increasing its accuracy. Two different calibration schemes have been tested using TLE data for about 50 objects during the year 2000. The schemes involve either height-dependent scale factors to the density or corrections to CIRA-72 model temperatures, which affect the density output based on a physical model. Both schemes have been applied with different spherical harmonic expansions of the parameters in latitude and local solar time. Five TLE objects, varying in perigee altitude between 280 and 530 km, were deliberately not used during calibration, in order to provide independent validation. Even with a single daily parameter, the RMS density model error along their tracks can already be reduced from the 30% to the 15% level. Adding additional parameters results in RMS errors lower than 12%.     

Spin motion determination of the Envisat satellite through laser ranging measurements from a single pass measured by a single station

The Satellite Laser Ranging (SLR) technology is used to accurately determine the position of space objects equipped with so-called retro-reflectors or retro-reflector arrays (RRA). This type of measurement allows to measure the range to the spacecraft with high precision, which leads to determination of very accurate orbits for these targets. Non-active spacecraft, which are not attitude controlled any longer, tend to start to spin or tumble under influence of the external and internal torques and forces.If the return signal is measured for a non-spherical non-active rotating object, the signal in the range residuals with respect to the reference orbit is more complex. For rotating objects the return signal shows an oscillating pattern or patterns caused by the RRA moving around the satellite’s centre of mass. This behaviour is projected onto the radial component measured by the SLR.In our work, we demonstrate how the SLR ranging technique from one sensor to a satellite equipped with a RRA can be used to precisely determine its spin motion during one passage. Multiple SLR measurements of one target over time allow to accurately monitor spin motion changes which can be further used for attitude predictions. We show our solutions of the spin motion determined for the non-active ESA satellite Envisat obtained from measurements acquired during years 2013–2015 by the Zimmerwald SLR station, Switzerland. All the necessary parameters are defined for our own so-called point-like model which describes the motion of a point in space around the satellite centre of mass.     

Drag balance Cubesat attitude motion effects on in-situ thermosphere density measurements

The dynamics of Cubesats carrying a drag balance instrument (DBI) for in situ atmosphere density measurements is analyzed. Atmospheric drag force is measured by the displacement of two light plates exposed to the incoming particle flow. This system is well suited for a distributed sensor network in orbit, to get simultaneous in situ local (non orbit averaged) measurements in multiple positions and orbit heights, contributing to the development and validation of global atmosphere models. The implementation of the DBI leads to orbit normal pointing spinning two body system. The use of a spin-magnetic attitude control system is suggested, based only on magnetometer readings, contributing to making the system simple, inexpensive, and reliable. It is shown, by an averaging technique, that this system provides for orbit normal spin axis pointing. The effect of the coupling between the attitude dynamics and the DBI is evaluated, analyzing its frequency content and showing that no frequency components arise, affecting the DBI performance. The analysis is confirmed by Monte Carlo numerical simulation results.     

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.     

由卫星阻力测定的大气密度值与MSIS和CIRA1972模式的部分比对

本文利用我国卫星的阻力资料测定的大气密度值, 分别与MSIS和CIRA1972模式进行比对.比对结果表明:在本文所计算的高度上, MSIS模式比CIRA1972更接近卫星阻力资料测定的密度值.      

Mass loading of planetary magnetospheres by rocky satellites—I: Production of rocky ejecta

Recent hypervelocity studies have been conducted which simulate the collision of interplanetary dust with rocky planetary satellite surfaces. Preliminary flux-mass distributions of micron and submicron ejecta from these hypervelocity impact studies have been determined. Several models of the flux-mass distribution of primary interplanetary dust are used to determine ratios of satellite surface ejecta and primary meteoroid flux-mass distributions. The results are used in a second model to determine the ejecta spatial mass densities near the surface of the satellite.     

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.