Numerical modeling of a Global Navigation Satellite System in a general relativistic framework

In this article we model a Global Navigation Satellite System (GNSS) in a Schwarzschild space–time, as a first approximation of the relativistic geometry around the Earth. The closed time-like and scattering light-like geodesics are obtained analytically, describing respectively trajectories of satellites and electromagnetic signals. We implement an algorithm to calculate Schwarzschild coordinates of a GNSS user who receives proper times sent by four satellites, knowing their orbital parameters; the inverse procedure is implemented to check for consistency. The constellation of satellites therefore realizes a geocentric inertial reference system with no a priori realization of a terrestrial reference frame. We perform a simulation of position determination and show that the determination of the four coordinates with a 25–32 digit accuracy takes only around 60 ms. Effects of non-gravitational perturbations on positioning errors are assessed, and methods to reduce them are sketched. In particular, inter-links between satellites could greatly enhance stability and accuracy of the positioning system. Effects of gravitational perturbations are omitted in this paper in order to make a clearer comparison between the relativistic and non-relativistic scheme, but they will be included in subsequent work. We believe that the final algorithm will be a serious alternative to the usual post-Newtonian scheme.     

Performance of GPS-based accelerometry: A simulation experiment

Recent studies have shown that with the availability of high-quality CHAMP and GRACE gravity field models, it is feasible to determine accurate non-gravitational accelerations for low Earth orbiting satellites indirectly from precise GPS satellite-to-satellite observations. Possible applications of this so-called GPS-based accelerometry approach consist of accelerometer calibration and atmospheric density and wind computations. With the growing number of high-quality space-borne GPS receivers, this method could be applied to a large range of satellites. In this paper an extensive simulation study has been carried out, based on real accelerometer data from the GRACE mission, in order to determine the optimal processing strategy and the resulting accuracy of the estimated non-gravitational accelerations. It is shown that the optimal processing strategy consists of a piecewise linear parameterization of the estimated empirical accelerations, together with short 6-h orbit arcs. The GPS-based accelerometry approach makes use of triple-differenced GPS observations and the impact of considering the correlated observation noise was found to be marginal in the presence of other error sources such as GPS ephemeris errors. Using a priori non-gravitational force models improves the recovery of low temporal resolution accelerations, except during huge geomagnetic storms. With this strategy, non-gravitational accelerations can be recovered during high solar activity with an accuracy of better than 10% of the total signal in along-track direction and around 25–40% in cross-track direction, at time resolutions of around 8–20 min. During solar minimum conditions, the relative recovery error will increase to approximately 50% in along-track direction and around 60–70% in cross-track direction, due to the reduced atmospheric drag signal. Unfortunately, GPS-based accelerometry is hardly sensitive in the radial direction.     

Impact of a novel hybrid accelerometer on satellite gravimetry performance

The state-of-the-art electrostatic accelerometers (EA) used for the retrieval of non-gravitational forces acting on a satellite constitute a core component of every dedicated gravity field mission. However, due to their difficult-to-control thermal drift in the low observation frequencies, they are also one of the most limiting factors of the achievable performance of gravity recovery. Recently, a hybrid accelerometer consisting of a regular EA and a novel cold atom interferometer (CAI) that features a time-invariant observation stability and constantly recalibrates the EA has been developed in order to remedy this major drawback. In this paper we aim to assess the value of the hybrid accelerometer for gravity field retrieval in the context of GRACE-type and Bender-type missions by means of numerical closed-loop simulations where possible noise specifications of the novel instrument are considered and different components of the Earth’s gravity field signal are added subsequently. It is shown that the quality of the gravity field solutions is mainly dependent on the CAI’s measurement accuracy. While a low CAI performance (10^?8 to 10^?9?m/s²/Hz^1/2) does not lead to any gains compared to a stand-alone EA, a sufficiently high one (10^?11?m/s²/Hz^1/2) may improve the retrieval performance by over one order of magnitude. We also show that improvements which are limited to low-frequency observations may even propagate into high spherical harmonic degrees. Further, the accelerometer performance seems to play a less prominent role if the overall observation geometry is improved as it is the case for a Bender-type mission. The impact of the accelerometer measurements diminishes further when temporal variations of the gravity field are introduced, pointing out the need for proper de-aliasing techniques. An additional study reveals that the hybrid accelerometer is – contrary to a stand-alone EA – widely unaffected by scale factor instabilities.     

Perspective on space radiation for space flights in 2020–2040

The Sun undergoes several well known periodicities in activity, such as the Schwabe 11 year cycle, the Gleissberg 80–90 year cycle, the Suess 200–210 year cycle and the Halstatt 2200–2300 year cycle. In addition, there is evidence that the 20th century levels of solar activity are unusually high. The years 2020–2040 are expected to coincide with increased activity in human space flight beyond low Earth orbit. The solar cycles and the present level of solar activity are reviewed and their activities during the years 2020–2040 are discussed with a perspective on space radiation and the future program of space flight. It is prudent to prepare for continuing levels of high solar activity as well as for the low levels of the current deep minimum, which has corresponded to high galactic cosmic ray flux.     

Electrostatic accelerometer with bias rejection for gravitation and Solar System physics

Radio tracking of interplanetary probes is an important tool for navigation purposes as well as for testing the laws of physics or exploring planetary environments. The addition of an accelerometer on board a spacecraft provides orbit determination specialists and physicists with an additional observable of great interest: it measures the value of the non-gravitational acceleration acting on the spacecraft, i.e. the departure of the probe from geodesic motion.     

Improved kHz-SLR tracking techniques and orbit quality analysis for LEO missions

Satellite gravity field missions such as CHAMP, GRACE and GOCE are designed as low Earth orbiting spacecraft (LEO) with orbit heights of about 250–500 km. The challenging mission objectives require a very precise knowledge of the satellite orbit position in space. For these missions precise orbit information is typically provided by GPS satellite-to-satellite tracking (SST) observations supported by satellite laser ranging (SLR).     

The importance of restoration of the atmospheric electrical environment in closed Bioregenerative Life Supporting Systems.

Earth based Bioregenerative Life Supporting Systems (BLSS) are subject to 4 main physical factors: gravity, light, temperature and electrical environment. The first 3 are obvious to everyone, the Electrical Environment (EE) is not under the majority of prevailing conditions perceived directly by our senses. The EE is one of the important physical factors directly influencing some plants and in a less obvious way also the majority of plants. There are only two long range forces in nature: the electromagnetic and the gravitational forces. Gravity is very much weaker than the electromagnetic forces FG/FEL=10(-38), where FG is the gravitational Force and FEL are the electromagnetic Forces. The atmospheric electric field prevails all the time over the entire Earth with a mean intensity of 130 V/m. It is therefore a potent factor which may be used by some plants exposed throughout their entire life time to the atmospheric electric field. What effect should the normal atmospheric electric field have on plants? All living plants are good electrical conductors for electrostatic fields. The plants distort the normally vertical field lines, which have to be perpendicular to the plant tissue everywhere in order to avoid the extraction of energy from the field. The meristems concentrate the field lines, thus the electrically charged nutrients are supplied to the growing parts of the plant exposed to the field. This results in electrotropism in some plants. It is very well known that plants do have adaptive capabilities as compared to animals, it is important for their survival, because they cannot run away from trouble. It is found by careful observations of the behaviour of different plants that some plants do respond to the presence of the atmospheric electric field while other plants exposed to the same environment are indifferent to the atmospheric electric field. The plants growing in the Biosphere II were shielded by the metal structure of the Biosphere II. Because these plants which do make use of the atmospheric electric field are also to be found in the Biosphere II and must be deprived of this natural atmospheric electric field and in consequence of some of their natural nutrients. I have experimental evidence that at least some plants do respond to the atmospheric electric field. This effect is the most likely cause of the oxygen depletion and of the carbon dioxide accumulation in Biosphere II. Under the micro gravity encountered in space habitats the restoration of electric fields is even more important for plant growth than it is on the Earth.     

Design and analysis of an extended mission of CE-2: From lunar orbit to Sun–Earth L2 region

Chang’E-2 (CE-2) has firstly successfully achieved the exploring mission from lunar orbit to Sun–Earth L2 region. In this paper, we discuss the design problem of transfer trajectory and at the same time analyze the visible segment of Tracking, Telemetry & Control (TT&C) system for this mission. Firstly, the four-body problem of Sun–Earth–Moon and Spacecraft can be decoupled in two different three-body problems (Sun–Earth + Moon Restricted Three-Body Problems (RTBPs) and Earth–Moon ephemeris model). Then, the transfer trajectory segments in different model are computed, respectively, and patched by Poincaré sections. The full-flight trajectory including transfer trajectory from lunar orbit to Sun–Earth L2 region and target Lissajous orbit is obtained by the differential correction method. Finally, the visibility of TT&C system at the key time is analyzed. Actual execution of CE-2 extended mission shows that the trajectory design of CE-2 mission is feasible.     

GRACE accelerometer calibration by high precision non-gravitational force modeling

The precise modeling and knowledge of non-gravitational forces acting on satellites is of big interest to many scientific tasks and missions. Since 2002, the twin GRACE satellites have measured these forces in a low Earth orbit with highly precise accelerometers, for about 15?years. Besides the significance for the GRACE mission, these measurement data allow the evaluation of modeling approaches and the improvement of force models. Unfortunately, before any scientific usage, the accelerometer measurements need to be calibrated, namely scale factor and bias have to be regularly estimated.In this study we demonstrate an accelerometer calibration approach, solely based on high precision non-gravitational force modeling without any use of empirically or stochastically estimated parameters, using our in-house developed satellite simulation tool XHPS. The aim of this work is twofold, first we use the accelerometer data and the residuals resulting from the calibration to quantitatively analyze and validate different non-gravitational force model approaches. In a second step, we compare the calibration results to three different calibration methods from different authors, based on gravity field recovery, GPS-based precise orbit determination, and based on modeled accelerations.We consider atmospheric drag forces and winds, as well as radiation forces due to solar radiation pressure, albedo, Earth infrared and thermal radiation (TRP) of the satellite itself. For TRP, we investigate different transient temperature calculation approaches for the satellite surfaces with absorbed power from the aforementioned radiation sources. A detailed finite element model of the satellite is utilized for every force, considering orientation, material properties and shadowing conditions for each element.For cross-track and radial direction, which are mainly affected by the radiative forces, our calibration residuals are quite small when drag is not super dominant (1–3?$\text{nm}/{\text{s}}^2$ for total accelerations around $\pm$50?$\text{nm}/{\text{s}}^2$). For these directions the calibration seems to perform better than the other compared methods, where some bigger differences were found. For the drag dominated along-track direction it is vice versa, here our method is not sensitive enough because the difference between modeled and measured drag is bigger (e.g. residuals around 10?$\text{nm}/{\text{s}}^2$ for total accelerations around $\pm$70?$\text{nm}/{\text{s}}^2$ for low solar activity). In along-track direction the orbit determination based methods are more sensitive and produce more reliable results. Results for the complete GRACE mission time span from 2003 to 2017 are shown, covering different seasonal environmental conditions.     

Magnetic torquer induced disturbing signals within GRACE accelerometer data

The GRACE (Gravity Recovery And Climate Experiment) gravity field satellite mission was launched in 2002. Although many investigations have been carried out, not all disturbances and perturbations upon satellite instruments and sensors are resolved yet. In this work the issue of acceleration disturbances onboard of GRACE due to magnetic torquers is investigated and discussed. Each of the GRACE satellites is equipped with a three-axes capacitive accelerometer to measure non-gravitational forces acting on the spacecraft. We used 10 Hz Level 1a raw accelerometer data in order to determine the impact of electric current changes on the accelerometer. After reducing signals which are induced by highly dominating processes in the low frequency range, such as thermospheric drag and solar radiation pressure, which can easily be done by applying a high-pass filter, disturbing signals from onboard instruments such as thruster firing events or heater switch events need to be removed from the previously filtered data. Afterwards the spikes which are induced by the torquers can be very well observed. Spikes vary in amplitude with respect to an increasing or decreasing current used for magnetic torquers, and can be as large as 20 nm/s². Furthermore, we were able to set up a model for the spikes of each scenario with which we were able to compute model spike time series. With these time series the spikes can successfully be removed from the 10 Hz raw accelerometer data. Spectral analysis of the time series reveal that an influence onto gravity field determination due to these effects is very unlikely, but can theoretically not be excluded.     

Parallel,adaptive, multi-object trajectory integrator for space simulation applications

Computer simulation is a very helpful approach for improving results from space born experiments. Initial-value problems (IVPs) can be applied for modeling dynamics of different objects – artificial Earth satellites, charged particles in magnetic and electric fields, charged or non-charged dust particles, space debris. An ordinary differential equations systems (ODESs) integrator based on applying different order embedded Runge–Kutta–Fehlberg methods is developed. These methods enable evaluation of the local error. Instead of step-size control based on local error evaluation, an optimal integration method is selected. Integration while meeting the required local error proceeds with constant-sized steps. This optimal scheme selection reduces the amount of calculation needed for solving the IVPs. In addition, for an implementation on a multi core processor and parallelization based on threads application, we describe how to solve multiple systems of IVPs efficiently in parallel.     

A Partially Orthogonal EnKF approach to atmospheric density estimation using orbital debris

A key requirement for accurate trajectory prediction and space situational awareness is knowledge of how non-conservative forces affect space object motion. These forces vary temporally and spatially, and are driven by the underlying behavior of space weather particularly in Low Earth Orbit (LEO). Existing trajectory prediction algorithms adjust space weather models based on calibration satellite observations. However, lack of sufficient data and mismodeling of non-conservative forces cause inaccuracies in space object motion prediction, especially for uncontrolled debris objects. The uncontrolled nature of debris objects makes them particularly sensitive to the variations in space weather. Our research takes advantage of this behavior by utilizing observations of debris objects to infer the space environment parameters influencing their motion.The hypothesis of this research is that it is possible to utilize debris objects as passive, indirect sensors of the space environment. We focus on estimating atmospheric density and its spatial variability to allow for more precise prediction of LEO object motion. The estimated density is parameterized as a grid of values, distributed by latitude and local sidereal time over a spherical shell encompassing Earth at a fixed altitude of 400 km. The position and velocity of each debris object are also estimated. A Partially Orthogonal Ensemble Kalman Filter (POEnKF) is used for assimilation of space object measurements to estimate density.For performance comparison, the scenario characteristics (number of objects, measurement cadence, etc.) are based on a sensor tasking campaign executed for the High Accuracy Satellite Drag Model project. The POEnKF analysis details spatial comparisons between the true and estimated density fields, and quantifies the improved accuracy in debris object motion predictions due to more accurate drag force models from density estimates. It is shown that there is an advantage to utilizing multiple debris objects instead of just one object. Although the work presented here explores the POEnKF performance when using information from only 16 debris objects, the research vision is to utilize information from all routinely observed debris objects. Overall, the filter demonstrates the ability to estimate density to within a threshold of accuracy dependent on measurement/sensor error. In the case of a geomagnetic storm, the filter is able to track the storm and provide more accurate density estimates than would be achieved using a simple exponential atmospheric density model or MSIS Atmospheric Model (when calm conditions are assumed).     

Instability of the present LEO satellite populations

Several studies conducted during 1991–2001 demonstrated, with some assumed launch rates, the future unintended growth potential of the Earth satellite population, resulting from random, accidental collisions among resident space objects. In some low Earth orbit (LEO) altitude regimes where the number density of satellites is above a critical spatial density, the production rate of new breakup debris due to collisions would exceed the loss of objects due to orbital decay.     

Might astronauts one day be treated like return samples?

The next time humans set foot on the Moon or another planet, will we treat the crew like we would a sample return mission when they come back to Earth? This may seem a surprising or even provocative question, but it is one we need to address. The hurdles and hazards of sending humans to Mars – for example, the technology constraints and physiological and psychological challenges – are many; but let us not forget the need to protect populations and environments from the risk of contamination [United Nations, treaty on principles governing the activities of states in the exploration and use of outer space, including the Moon and other celestial bodies (the “Outer Space Treaty”) referenced 610 UNTS 205 - resolution 2222(XXI) of December 1966].     

R3D-B2 – Measurement of ionizing and solar radiation in open space in the BIOPAN 5 facility outside the FOTON M2 satellite

Solar and space radiation have been monitored using the R3D-B2 radiation risks radiometer-dosimeter on board a recent space flight on the Russian satellite Foton M2 within the ESA Biopan 5 facility mounted on the outside of the satellite exposed to space conditions. The solar radiation has been assayed in four wavelength bands (UV-C, 170–280 nm, UV-B, 280–315 nm), UV-A (315–400 nm) and PAR (photosynthetic active radiation, 400–700 nm). The data show an increasing tumbling rotation of the satellite during the mission. The photodiodes do not show a cosine response to the incident light which has been corrected. After calibration of the signals using the extraterrestrial spectrum, doses have been calculated for each orbit, for each day and for the total mission as basic data for the biological material which has been exposed in parallel in the Biopan facility. Cosmic ionizing radiation has been monitored and separated in 256 deposited energy spectra, which were further used for determination of the absorbed dose rate and flux. Basic data tables were prepared to be used by other Biopan 5 experiments. The paper summarizes the results for the Earth radiation environment at the altitude (262–304 km) of the Foton M2 spacecraft. Comparisons with the predictions of NASA Earth radiation environment experimental models AE-8 and AP-8, and the PSB97 model are also presented, which calculate the fluxes of ionizing radiation from a simulation. AP-8 is a model for trapped radiation.     

Application of China–Brazil Earth resources satellite in China

The launch and successful operation of Chinese–Brazil Earth resources satellite (CBERS-1) in China has accelerated the application of space technology in China. These applications include agriculture, forestry, water conservation, land resources, city planning, environment protection and natural hazards monitoring and so on. The result of these applications provides a scientific basis for government decision making and has created great economic and social benefits in Chinese national economy construction. In this paper we present examples and provide auxiliary documentation of additional applications of the data from Earth resource monitoring.     

Recovery of the Earth’s gravity field from formation-flying satellites: Temporal aliasing issues

Temporal and mean gravity field models derived from the twin-satellite, leader–follower mission GRACE have provided a new type of information for Earth sciences. In this contribution, we study the potential of various alternative satellite formations for gravity field determination in the post-GRACE era in a simulation environment. In particular, the effects of spherical harmonic truncation and of temporal aliasing in the processing of gravity products from such future formations are investigated.     

Dynamical cartography of Earth satellite orbits

We have carried out a numerical investigation of the coupled gravitational and non-gravitational perturbations acting on Earth satellite orbits in an extensive grid, covering the whole circumterrestrial space, using an appropriately modified version of the SWIFT symplectic integrator, which is suitable for long-term ($～$120?years) integrations of the non-averaged equations of motion. Hence, we characterize the long-term dynamics and the phase-space structure of the Earth-orbiter environment, starting from low altitudes ($～$400?km) and going up to the GEO region and beyond. This investigation was done in the framework of the EC-funded “ReDSHIFT” project, with the purpose of enabling the definition of passive debris removal strategies, based on the use of physical mechanisms inherent in the complex dynamics of the problem (i.e., resonances). Accordingly, the complicated interactions among resonances, generated by different perturbing forces (i.e., lunisolar gravity, solar radiation pressure, tesseral harmonics in the geopotential) are accurately depicted in our results, where we can identify the regions of phase space where the motion is regular and long-term stable and regions for which eccentricity growth and even instability due to chaotic behavior can emerge. The results are presented in an “atlas” of dynamical stability maps for different orbital zones, with a particular focus on the (drag-free) range of semimajor axes, where the perturbing effects of the Earth’s oblateness and lunisolar gravity are of comparable order. In some regions, the overlapping of the predominant lunisolar secular and semi-secular resonances furnish a number of interesting disposal hatches at moderate to low eccentricity orbits. All computations were repeated for an increased area-to-mass ratio, simulating the case of a satellite equipped with an on-board, area-augmenting device. We find that this would generally promote the deorbiting process, particularly at the transition region between LEO and MEO. Although direct reentry from very low eccentricities is very unlikely in most cases of interest, we find that a modest “delta-v” ($\mathrm{\Delta }V$) budget would be enough for satellites to be steered into a relatively short-lived resonance and achieve reentry into the Earth’s atmosphere within reasonable timescales ($～$50?years).