Semi-analytical investigations of high area-to-mass ratio geosynchronous space debris including Earth’s shadowing effects

This paper investigates the long-term perturbations of the orbits of geosynchronous space debris influenced by direct radiation pressure including the Earth’s shadowing effects. For this purpose, we propose an extension of our homemade semi-analytical theory [Valk, S., Lemaître, A., Deleflie, F. Semi-analytical theory of mean orbital motion for geosynchronous space debris under gravitational influence. Adv. Space Res., submitted for publication], based on the method developed by Aksnes [Aksnes, K. Short-period and long-period perturbations of a spherical satellite due to direct solar radiation. Celest. Mech. Dyn. Astron. 13, 89–104, 1976] and generalized into a more convenient non-singular formalism. The perturbations accounting for the direct radiation pressure with the Earth’s shadow are computed on a revolution-by-revolution basis, retaining the original osculating Hamiltonian disturbing function. In this framework, we compute the non-singular mean longitude at shadow entry and shadow exit at every orbital revolution in opposition to classical approaches where the singular eccentric anomalies at shadow entry and shadow exit are computed. This new algorithm is developed using non-singular variables. Consequently, it is particularly suitable for both near-circular and near-equatorial orbits as well as orbits which transit periodically around null eccentricities and null inclinations.The algorithm is tested by means of numerical integrations of the equations, averaged over the short periods, including radiation pressure, J₂, the combined Moon and Sun third body attraction as well as the long-term effects of the 1:1 resonance occurring for geosynchronous objects. As an extension of [Valk, S., Lemaître, A., Anselmo, L. Analytical and semi-analytical investigations of geosynchronous space debris with high area-to-mass ratios influenced by solar radiation pressure. Adv. Space Res., doi:10.1016/j.asr.2007.10.025, 2007b], we especially apply our analysis to space debris with area-to-mass as high as 20 m²/kg. This paper provides numerical and semi-analytical investigations leading to a deep understanding of the long-term evolution of the semi-major axis. Finally, these semi-analytical investigations are compared with accurate numerical integrations of the osculating equations of motion over time scales as high as 25 years.     

Influence of Earth’s shadowing effects on space debris stability

Solar radiation pressure affects the evolution of high area-to-mass geostationary space debris. In this work, we extend the stability study of Valk et al. (2009) by considering the influence of Earth’s shadows on the short- and long-term time evolutions of space debris. To assess the orbits stability, we use the Mean Exponential Growth factor of Nearby Orbits (MEGNO), which is an efficient numerical tool to distinguish between regular and chaotic behaviors. To reliably compute long-term space debris motion, we resort to the Global Symplectic Integrator (GSI) of Libert et al. (2011) which consists in the symplectic integration of both Hamiltonian equations of motion and variational equations. We show how to efficiently compute the MEGNO indicator in a complete symplectic framework, and we also discuss the choice of a symplectic integrator, since propagators adapted to the structure of the Hamiltonian equations of motion are not necessarily suited for the associated variational equations. The performances of our method are illustrated and validated through the study of the Arnold diffusion problem. We then analyze the effects of Earth’s shadows, using the adapted conical and cylindrical Earth’s shadowing models introduced by Hubaux et al. (2012) as the smooth shadow function deriving from these models can be easily included into the variational equations. Our stability study shows that Earth’s shadows greatly affect the global behaviour of space debris orbits by increasing the size of chaotic regions around the geostationary altitude. We also emphasize the differences in the results given by conical or cylindrical Earth’s shadowing models. Finally, such results are compared with a non-symplectic integration scheme.     

One dimensional blood flow in a planetocentric orbit

All life on earth is accustomed to the presence of gravity. When gravity is altered, biological processes can go awry. It is of great importance to ensure safety during a spaceflight. Long term exposure to microgravity can trigger detrimental physiological responses in the human body. Fluid redistribution coupled with fluid loss is one of the effects. In particular, in microgravity blood volume is shifted towards the thorax and head. Sympathetic nervous system-induced vasoconstriction is needed to maintain arterial pressure, while venoconstriction limits venous pooling of blood prevents further reductions in venous return of blood to the heart. In this paper, we modify an existing one dimensional blood flow model with the inclusion of the hydrostatic pressure gradient that further depends on the gravitational field modified by the oblateness and rotation of the Earth. We find that the velocity of the blood flow V_B is inversely proportional to the blood specific volume d, also proportional to the oblateness harmonic coefficient J₂, the angular velocity of the Earth ω_E, and finally proportional to an arbitrary constant c. For c = −0.39073 and ξ_H = −0.5 mmHg, all orbits result to less blood flow velocities than that calculated on the surface of the Earth. From all considered orbits, elliptical polar orbit of eccentricity e = 0.2 exhibit the largest flow velocity V_B = 1.031 m/s, followed by the orbits of inclination i = 45°and 0°. The Earth’s oblateness and its rotation contribute a 0.7% difference to the blood flow velocity.     

On the large-scale structure of the tail current as measured by THEMIS

The magnetic field structure and the spatial characteristics of the large-scale currents in the magnetospheric tail were studied during quiet and moderately disturbed geomagnetic conditions in 2009. The magnetic field of the currents other than the tail current was calculated in terms of a paraboloid model of the Earth’s magnetosphere, A2000, and was subtracted from measurements. It was found on the base of obtained tail current magnetic field radial distribution that the inner edge of the tail current sheet is located in the night side magnetosphere, at distances of about 10 R_E and of about 7 R_E during quiet and disturbed periods respectively. During the disturbance of February 14, 2009 (Dst_min ∼ −35 nT), the B_x and the B_z component of the tail current magnetic field near its inner edge were about 60 nT, and −60 nT that means that strong cross-tail current have been developed. The tail current parameters at different time moments during February 14, 2009 have been estimated. Solar wind conditions during this event were consistent with those during moderate magnetic storms with minimum Dst of about −100 nT. However, the magnetospheric current systems (magnetopause and cross-tail currents) were located at larger geocentric distances than typical during the 2009 extremely quiet epoch and did not provide the expected Dst magnitude. Very small disturbance on the Earth’s surface was detected consistent with an “inflated” magnetosphere.     

Symplectic integration of space debris motion considering several Earth’s shadowing models

In this work, we present a symplectic integration scheme to numerically compute space debris motion. Such an integrator is particularly suitable to obtain reliable trajectories of objects lying on high orbits, especially geostationary ones. Indeed, it has already been demonstrated that such objects could stay there for hundreds of years. Our model takes into account the Earth’s gravitational potential, luni-solar and planetary gravitational perturbations and direct solar radiation pressure. Based on the analysis of the energy conservation and on a comparison with a high order non-symplectic integrator, we show that our algorithm allows us to use large time steps and keep accurate results. We also propose an innovative method to model Earth’s shadow crossings by means of a smooth shadow function. In the particular framework of symplectic integration, such a function needs to be included analytically in the equations of motion in order to prevent numerical drifts of the energy. For the sake of completeness, both cylindrical shadows and penumbra transitions models are considered. We show that both models are not equivalent and that big discrepancies actually appear between associated orbits, especially for high area-to-mass ratios.     

Particle acceleration by stochastic fluctuations and dawn-dusk electric field in the Earth’s magnetotail

This work is devoted to investigate the interaction between protons and stochastic time-dependent electromagnetic fields generated by oscillating clouds of finite size, randomly positioned in the x–y  -plane. The geometry of the system is two-dimensional and, beside the time-dependent electromagnetic fluctuations, a steady-state, dawn-dusk electric field, _Ey$E_y$, has been added along the y-direction. The simultaneous presence of the stochastic time-dependent fluctuations and of the constant electric field component in the same system gives rise to two types of acceleration mechanisms operating on test particles: a second order Fermi-like process and a direct acceleration. By performing a parametric study, we extensively study the contribution of the two processes to proton acceleration. The energy values reached by test particles in this simple model are in good agreement with those observed in the Earth’s magnetotail region. Possible applications to the problem of particle acceleration in the terrestrial magnetosphere are widely discussed, and guidelines for future works are drawn.     

Symmetry,reduction and relative equilibria of a rigid body in the J2 problem

The J₂ problem is an important problem in celestial mechanics, orbital dynamics and orbital design of spacecraft, as non-spherical mass distribution of the celestial body is taken into account. In this paper, the J₂ problem is generalized to the motion of a rigid body in a J₂ gravitational field. The relative equilibria are studied by using geometric mechanics. A Poisson reduction process is carried out by means of the symmetry. Non-canonical Hamiltonian structure and equations of motion of the reduced system are obtained. The basic geometrical properties of the relative equilibria are given through some analyses on the equilibrium conditions. Then we restrict to the zeroth and second-order approximations of the gravitational potential. Under these approximations, the existence and detailed properties of the relative equilibria are investigated. The orbit–rotation coupling of the rigid body is discussed. It is found that under the second-order approximation, there exists a classical type of relative equilibria except when the rigid body is near the surface of the central body and the central body is very elongated. Another non-classical type of relative equilibria can exist when the central body is elongated enough and has a low average density. The non-classical type of relative equilibria in our paper is distinct from the non-Lagrangian relative equilibria in the spherically-simplified Full Two Body Problem, which cannot exist under the second-order approximation. Our results also extend the previous results on the classical type of relative equilibria in the spherically-simplified Full Two Body Problem by taking into account the oblateness of the primary body. The results on relative equilibria are useful for studies on the motion of many natural satellites, whose motion are close to the relative equilibria.     

Heliotropic dust rings for Earth climate engineering

This paper examines the concept of a Sun-pointing elliptical Earth ring comprised of dust grains to offset global warming. A new family of non-Keplerian periodic orbits, under the effects of solar radiation pressure and the Earth’s J₂ oblateness perturbation, is used to increase the lifetime of the passive cloud of particles and, thus, increase the efficiency of this geoengineering strategy. An analytical model is used to predict the orbit evolution of the dust ring due to solar-radiation pressure and the J₂ effect. The attenuation of the solar radiation can then be calculated from the ring model. In comparison to circular orbits, eccentric orbits yield a more stable environment for small grain sizes and therefore achieve higher efficiencies when the orbit decay of the material is considered. Moreover, the novel orbital dynamics experienced by high area-to-mass ratio objects, influenced by solar radiation pressure and the J₂ effect, ensure the ring will maintain a permanent heliotropic shape, with dust spending the largest portion of time on the Sun facing side of the orbit. It is envisaged that small dust grains can be released from a circular generator orbit with an initial impulse to enter an eccentric orbit with Sun-facing apogee. Finally, a lowest estimate of 1 × 10¹² kg of material is computed as the total mass required to offset the effects of global warming.     

Universal scaling laws for fully-developed magnetic field turbulence near and far upstream of the Earth’s bow shock

We analyze the multifractal scaling of the modulus of the interplanetary magnetic field near and far upstream of the Earth’s bow shock, measured by Cluster and ACE, respectively, from 1 to 3 February 2002. The maximum order of the structure function is carefully estimated for each time series using two different techniques, to ensure the validity of our high-order statistics. The first technique consists of plotting the integrand of the pth order structure function, and the second technique is a quantitative method which relies on the power-law scaling of the extreme events. We compare the scaling exponents computed from the structure functions of magnetic field differences with the predictions obtained by the She–Lévêque model of intermittency in anisotropic magnetohydrodynamic turbulence. Our results show a good agreement between the model and the observations near and far upstream of the Earth’s bow shock, rendering support for the modelling of universal scaling laws based on the Kolmogorov phenomenology in the presence of sheet-like dissipative structures.     

Determination of reconnected flux via remote sensing

Magnetic reconnection is one of the most fundamental processes in the magnetosphere. We present here a simple method to determine the essential parameters of reconnection such as reconnected flux and location of the reconnection site out of single spacecraft data via remote sensing. On the basis of a time-dependent reconnection model, the dependence of the reconnected flux on the magnetic field z-component B_z is shown. The integral of B_z over time is proportional to the reconnected flux and depends on the distance between the reconnection site and the actual position where B_z is measured. This distance can be estimated from analysis of magnetic field B_z data. We apply our method to Cluster measurements in the Earth’s magnetotail.     

Global comparison of the model results of GSM TIP with IRI for summer conditions

This paper presents the results of the numerical calculations thermosphere/ionosphere parameters which were executed with using of the Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP)and comparison of these results with empirically-based model IRI-2001. Model GSM TIP was developed in West Department of IZMIRAN and solves self-consistently the time-dependent, 3-D coupled equations of the momentum, energy and continuity for neutral particles (O₂, N₂, O), ions (O⁺, H⁺), molecular ions (M⁺) and electrons and largescale eletric field of the dynamo and magnetospheric origin in the range of height from 80 km to 15 Earth’s radii. The empirically derived IRI model describes the E and F regions of the ionosphere in terms of location, time, solar activity and season. Its output provides a global specification not only of N_e but also on the ion and electron temperatures and the ion composition. These two models represent a unique set of capabilities that reflect major differences in along with a substantial approaches of the first-principles model and global database model for the mapping ionosphere parameters. We focus on global distribution of the N_e, T_i, T_e and TEC for the one moment UT and fixed altitudes: 110 km, h_mF2, 300 km and 1000 km. The calculations were executed with using of GSM TIP and IRI models for August 1999, moderate solar activity and quiet geomagnetic conditions. Results present as the global differences between the IRI and GSM TIP models predictions. The discrepancies between model results are discussed.     

Synoptic magnetic field in cycle 23 and in the beginning of the cycle 24

The SOHO/MDI data provide the uniform time series of the synoptic magnetic maps which cover the period of the cycle 23 and the beginning of the cycle 24. It is very interesting period because of the long and deep solar minimum between the cycles 23 and 24. Synoptic structure of the solar magnetic field shows variability during solar cycles. It is known that the magnetic activity contributes to the solar irradiance. The axisymmetrical distribution of the magnetic flux (Fig. 3c) is closely associated with the ‘butterfly’ diagram in the EUV emission (Benevolenskaya et al., 2001). And, also, the magnetic field (B_∥) shows the non-uniform distributions of the solar activity with longitude, so-called ‘active zones’, and ‘coronal holes’ in the mid-latitude. Polar coronal holes are forming after the solar maxima and they persist during the solar minima. SOHO/EIT data in the emission of Fe XII (195 Å) could be a proxy for the coronal holes tracking. The active longitudinal zones or active longitude exist due to the reappearance of the activity and it is clearly seen in the synoptic structure of the solar cycle. On the descending branch of the solar cycle 23 active zones are less pronounced comparing with previous cycles 20, 21 and 22. Moreover, the weak polar magnetic field precedes the long and deep solar minimum. In this paper we have discussed the development of solar cycles 23 and 24 in details.     

Contribution of satellite laser ranging to combined gravity field models

In the framework of satellite-only gravity field modeling, satellite laser ranging (SLR) data is typically exploited to recover long-wavelength features. This contribution provides a detailed discussion of the SLR component of GOCO02S, the latest release of combined models within the GOCO series. Over a period of five years (January 2006 to December 2010), observations to LAGEOS-1, LAGEOS-2, Ajisai, Stella, and Starlette were analyzed. We conducted a series of closed-loop simulations and found that estimating monthly sets of spherical harmonic coefficients beyond degree five leads to exceedingly ill-posed normal equation systems. Therefore, we adopted degree five as the spectral resolution for real data analysis. We compared our monthly coefficient estimates of degree two with SLR and Gravity Recovery and Climate Experiment (GRACE) time series provided by the Center for Space Research (CSR) at Austin, Texas. Significant deviations in C₂₀ were noted between SLR and GRACE; the agreement is better for the non-zonal coefficients. Fitting sinusoids together with a linear trend to our C₂₀ time series yielded a rate of (−1.75 ± 0.6) × 10⁻¹¹/yr; this drift is equivalent to a geoid change from pole to equator of 0.35 ± 0.12 mm/yr or an apparent Greenland mass loss of 178.5 ± 61.2 km³/yr. The mean of all monthly solutions, averaged over the five-year period, served as input for the satellite-only model GOCO02S. The contribution of SLR to the combined gravity field model is highest for C₂₀, and hence is essential for the determination of the Earth’s oblateness.     

Measurements of the radiation quality factor Q at aviation altitudes during solar minimum (2006–2008)

In radiation protection, the Q-factor has been defined to describe the biological effectiveness of the energy deposition or absorbed dose to humans in the mixed radiation fields at aviation altitudes. This particular radiation field is generated by the interactions of primary cosmic particles with the atoms of the constituents of the Earth’s atmosphere. Thus the intensity, characterized by the ambient dose equivalent rate H∗(10), depends on the flight altitude and the energy spectra of the particles, mainly protons and alpha particles, impinging on the atmosphere. These charged cosmic projectiles are deflected both by the interplanetary and the Earth’s magnetic field such that the corresponding energy spectra are modulated by these fields. The solar minimum is a time period of particular interest since the interplanetary magnetic field is weakest within the 11-year solar cycle and the dose rates at aviation altitudes reach their maximum due to the reduced shielding of galactic cosmic radiation. For this reason, the German Aerospace Center (DLR) performed repeated dosimetric on-board measurements in cooperation with several German airlines during the past solar minimum from March 2006 to August 2008. The Q-factors measured with a TEPC range from 1.98 at the equator to 2.60 in the polar region.     

Swarm – An Earth Observation Mission investigating Geospace

The Swarm mission was selected as the 5th mission in ESA’s Earth Explorer Programme in 2004. This mission aims at measuring the Earth’s magnetic field with unprecedented accuracy. This will be done by a constellation of three satellites, where two will fly at lower altitude, measuring the gradient of the magnetic field, and one satellite will fly at higher altitude. The measured magnetic field is the sum of many contributions including both magnetic fields and currents in the Earth’s interior and electrical currents in Geospace. In order to separate all these sources electric field and plasma measurements will also be made to complement the primary magnetic field measurements. Together these will allow the deduction of information on a series of solid earth processes responsible for the creation of the fields measured. The completeness of the measurements on each satellite and the constellation aspect, however, implies simultaneous observations of a unique set of important electrodynamical parameters crucial for the understanding of the physical processes in Geospace, which are an important part of the objectives of the International Living With a Star Programme, ILWS. In this paper an overview of the Swarm science objectives, the mission concept, the scientific instrumentation, and the expected contribution to the ILWS programme will be summarized.     

Thermal neutrons’ response to the GLEs

Thermal neutrons’ flux near the Earth’s crust is very sensitive regarding different processes and phenomena both in the near-Earth space and in the Earth’s crust by reason of the dual nature of the thermal neutron flux. Its first source is associated with high-energy particles of cosmic rays penetrating into the Earth’s atmosphere and interacting with its elements. The second source originates from the radioactive gases contained in the Earth’s crust. So the contributions of these two sources are specified by the Earth’s crust conditions and its movements, on one hand and variations of high-energy particles flux near the Earth.     

Semi-analytical theory of mean orbital motion for geosynchronous space debris under gravitational influence

This paper provides a hamiltonian formulation of the equations of motion of an artificial satellite or space debris orbiting the geostationary ring. This theory of order 1 has been formulated using canonical and non-singular elements for eccentricity and inclination. The analysis is based on an expansion in powers of the eccentricity and of the inclination. The theory accounts for the influence of the Earth gravity field expanded in spherical harmonics, paying a particular attention to the resonance occurring for geosynchronous objects. The luni-solar perturbations are also taken into account. We present the resonant motion and its main characteristics: equilibria, stability, fundamental frequencies and width of the resonant area by comparison with a basic analytical model. Finally, we show some results concerning the long term dynamics of a typical space debris under the influence of the gravitational field of the Earth and the luni-solar interactions.     

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).     

Evolution of periodic orbits near the Lagrangian point L2

A study of the evolution of the periodic and the quasi-periodic orbits near the Lagrangian point L₂, which is located to the right of the smaller primary on the line joining the primaries and whose distance from the more massive primary is greater than the distance between the primaries, in the framework of restricted three-body problem for the Sun–Jupiter, Earth–Moon (relatively large mass ratio) and Saturn–Titan (relatively small mass ratio) systems is made. Two families of periodic orbits around the smaller primary are identified using the Poincaré surface of section method – family I (initially elliptical, gradually becomes egg-shaped with the increase in the Jacobi constant C and elongated towards the more massive primary) and family II (initially egg-shaped orbits elongated towards L₂ and gradually becomes elliptical with the increase in C). The family I in the Sun–Jupiter and Saturn–Titan systems contains two separatrix caused by third-order and fourth-order resonances, while the Earth–Moon system has only one separatrix which is caused by third-order resonances. Also in the Sun–Jupiter and the Saturn–Titan systems, family I merge with family II, around Jacobian constant 3.0393 and 3.0163, respectively, while in the Earth–Moon system, family II evolves separately from two different branches. The two branches merge at C = 3.184515. In the Earth–Moon system, the family II contains a separatrix due to third-order resonances which is absent in the other two systems.