Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
The Dynamical Origin of the Graviton Mass in the Non-Linear Theory of Massive Gravity
Next Article in Special Issue
A Theory of Inertia Based on Mach’s Principle
Previous Article in Journal / Special Issue
Gravitational Radiation, Vorticity And Super–Energy: A Conspicuous Threesome
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A HERO for General Relativity

Ministero dell’Istruzione, dell’Università e della Ricerca (M.I.U.R.)-Istruzione, Viale Unità di Italia 68, 70125 Bari (BA), Italy
Universe 2019, 5(7), 165; https://doi.org/10.3390/universe5070165
Submission received: 14 June 2019 / Revised: 29 June 2019 / Accepted: 1 July 2019 / Published: 5 July 2019
(This article belongs to the Special Issue Rotation Effects in Relativity)

Abstract

:
HERO (Highly Eccentric Relativity Orbiter) is a space-based mission concept aimed to perform several tests of post-Newtonian gravity around the Earth with a preferably drag-free spacecraft moving along a highly elliptical path fixed in its plane undergoing a relatively fast secular precession. We considered two possible scenarios—a fast, 4-h orbit with high perigee height of 1047 km and a slow, 21-h path with a low perigee height of 642 km . HERO may detect, for the first time, the post-Newtonian orbital effects induced by the mass quadrupole moment J 2 of the Earth which, among other things, affects the semimajor axis a via a secular trend of ≃4–12 cm yr 1 , depending on the orbital configuration. Recently, the secular decay of the semimajor axis of the passive satellite LARES was measured with an error as little as 0 . 7 cm yr 1 . Also the post-Newtonian spin dipole (Lense-Thirring) and mass monopole (Schwarzschild) effects could be tested to a high accuracy depending on the level of compensation of the non-gravitational perturbations, not treated here. Moreover, the large eccentricity of the orbit would allow one to constrain several long-range modified models of gravity and accurately measure the gravitational red-shift as well. Each of the six Keplerian orbital elements could be individually monitored to extract the G J 2 / c 2 signature, or they could be suitably combined in order to disentangle the post-Newtonian effect(s) of interest from the competing mismodeled Newtonian secular precessions induced by the zonal harmonic multipoles J of the geopotential. In the latter case, the systematic uncertainty due to the current formal errors σ J of a recent global Earth’s gravity field model are better than 1 % for all the post-Newtonian effects considered, with a peak of 10 7 for the Schwarzschild-like shifts. Instead, the gravitomagnetic spin octupole precessions are too small to be detectable.

1. Introduction

The (slow) motion of a test particle moving in spacetime (weakly) deformed by the mass-energy content of an isolated, axially symmetric rotating body of mass M, angular momentum S , polar and equatorial radii R p , R e , ellipticity ε = 1 R p 2 / R e 2 , dimensionless quadrupole mass moment J 2 exhibits several post-Newtonian (pN) features. Some of them have never been put to the test so far because of their smallness; they are the gravitoelectric and gravitomagnetic effects associated with the asphericity of the central body induced by its mass quadrupole and spin octupole moments, respectively [1,2,3,4,5,6].
Instead, the pN orbital effects which have been extensively tested so far in several terrestrial and astronomical scenarios are the gravitoelectric and gravitomagnetic precessions due to the mass monopole and spin dipole moments, respectively. The former is responsible for the time-honored, previously anomalous perihelion precession of Mercury [7], whose explanation by Reference [8] was the first empirical success of his newly born theory of gravitation. It was later repeatedly measured with radar measurements of Mercury itself [9,10], of other inner planets [11,12], and of the asteroid Icarus [13,14] as well. Binary pulsars [15] and Earth’s artificial satellites [16,17] have also been used. The latter is the so-called Lense-Thirring effect [18,19,20,21,22], which is currently under scrutiny in the Earth’s surrounding with the geodetic satellites of the LAGEOS family; see, for example, Renzetti [23], and Lucchesi et al. [24], and references therein. Another gravitomagnetic effect—the Pugh-Schiff rates of change of orbiting gyroscopes [25,26]—was successfully tested in the field of the Earth with the dedicated Gravity Probe B (GP-B) spaceborne mission a few years ago [27,28] to a 19 % accuracy level, despite the originally expected error level of the order of 1 % [29].
By assuming the validity of general relativity, its mass quadrupole and spin octupole accelerations are, to the first post-Newtonian (pN) order,
A pN M J 2 = μ J 2 R e 2 c 2 r 4 3 2 5 ξ 2 1 r ^ 2 ξ S ^ v 2 4 μ r 6 5 ξ 2 1 v r 2 ξ v S v 2 μ r 3 ξ 2 1 r ^ ,
A pN S ε = 3 G S R e 2 ε 2 7 c 2 r 5 v × 5 ξ 7 ξ 2 3 r ^ + 3 1 5 ξ 2 S ^ ,
where G is the Newtonian constant of gravitation, μ G M is the gravitational parameter of the primary, c is the speed of light in vacuum, S ^ is the unit vector of the rotational axis,
ξ S ^ · r ^
is the cosine of the angle between the directions of the body’s angular momentum and the orbiter’s position vector r , and
v r v · r ^ ,
v S v · S ^
are the components of the particle’s velocity v along the radial direction and the primary’s spin respectively. The averaged rates of change of the semimajor axis a, the eccentricity e, the inclination I, the longitude of the ascending node Ω and the argument of pericenter ω induced by Equations (1) and (2) were analytically calculated for a general orientation of S ^ in space by Iorio [30,31]; previous derivations of the gravitoelectric mass quadrupole effects in the particular case of an equatorial coordinate system with its reference z axis aligned with S ^ can be found in Brumberg [1], Soffel [3], Will [4], Soffel et al. [6].
In this paper, we will preliminarily explore the perspectives of measuring, for the first time, some consequences of Equations (1) and (2) by suitably designing a dedicated drag-free satellite-based mission around the Earth encompassing a highly eccentric geocentric orbit exploiting the frozen perigee configuration; we provisionally name it as “HERO” (Highly Eccentric Relativity Orbiter). For some embryonal thoughts about the possibility of using an Earth’s spacecraft to measure the pN gravitoelectric effects proportional to G J 2 / c 2 , see Iorio [30,32]; for deeper investigations concerning a possible probe around Jupiter to measure them and the pN gravitomagnetic signature proportional to G S ε 2 / c 2 , see Iorio [31,32]. About the propagation of the electromagnetic waves in the deformed spacetime of an oblate body and the perspectives of measuring the resulting deflection due to Jupiter with astrometric techniques, see, e.g., Abbas, Bucciarelli & Lattanzi [33], Crosta & Mignard [34], Kopeikin & Makarov [35], Le Poncin-Lafitte & Teyssandier [36], and references therein. Other more or less similar mission concepts can be found, for example, in Angélil et al. [37], Schärer et al. [38,39]. We will show that the size of the secular rate of a predicted by Equation (1) falls within the recently reached experimental accuracy in measuring phenomenologically such a kind of an effect with the existing passive LARES satellite [24]. Be that as it may, we will show that, as a by-product, other general relativistic features of motion could also be measured with high accuracy, at least as far as the systematic error due to the current formal level of mismodeling in the competing classical precessions due to the zonal harmonic coefficients J , = 2 , 3 , 4 , of the multipolar expansion of the Earth’s gravity potential is concerned. To this aim, it is crucial to assess the level of possible cancellation of the non-gravitational perturbations by the drag-free technology. Its evaluation is outside the scope of the present paper. However, we will look in detail at the atmospheric drag, which is one of the major disturbing non-conservative accelerations inducing relevant competing signatures, especially on a. For the sake of simplicity, we will assume a spherical shape for a passive spacecraft.
The high eccentricity1 of the suggested orbit of HERO would allow also for accurate tests of the gravitational redshift provided that the spacecraft is endowed with accurate atomic clocks; see Table A2 and Table A6 for the expected sizes of it. For recent tests of such an effect performed with the H-maser clocks carried onboard the satellites GSAT0201 (5-Doresa) and GSAT0202 (6-Milena) of the Galileo constellation by exploiting their fortuitous rather high eccentricity due to their erroneous orbital injection, see Delva et al. [40], Herrmann et al. [41]. It is interesting to recall that the formerly proposed Space-Time Explorer and QUantum Equivalence Space Test (STE-QUEST) space mission [42] was pre-selected by the European Space Agency (ESA) for the Cosmic Vision M3 launch opportunity; although it was not finally selected due to budgetary and technological reasons, its science case was highly rated. HERO may somehow inherit (part of) its legacy.
Finally, several long-range modified models of gravity [43] imply spherically symmetric modifications of the Newtonian inverse-square law which induce net secular precessions of the pericenter and the mean anomaly at epoch. They would represent further valuable goals for HERO.
The paper is organized as follows. Section 2 is the main body of the paper; it contains a discussion of two possible orbital configurations along with the magnitude of the various pN effects and the size of the corresponding systematic errors due to the current level of mismodeling in the multipolar zonal coefficients of the geopotential. It also deals with the linear combination approach which could be implemented in order to reduce the bias due to the latter ones. Section 3 contains a summary of our findings and offer our conclusions. Appendix A contains the tables and the figures. Appendix B displays the analytically calculated orbital precessions due to the zonal harmonics of the geopotential up to degree = 8 . Appendix C deals in detail with the impact of the atmospheric drag on all the orbital elements of a spherical, passive geodetic satellite in a highly elliptical orbit.

2. Two Different Orbital Configurations for HERO

In Table A2 and Table A6, two different orbital configurations are proposed. They imply highly eccentric orbits, characterized by values of the eccentricity as large as e = 0 . 45 and e = 0 . 82 , respectively, and the critical inclination I crit = arcsin 2 / 5 which allows one to keep the argument of perigee ω essentially fixed over any reasonable time span for an actual data analysis and the longitude of the ascending node Ω circulating at a relatively high pace. Their orbital periods are P b = 4 . 3 h and P b = 21 . 3 h , with perigee heights of h min = 1047 km and h min = 642 km , respectively.
One of the most interesting relativistic features of motion is, perhaps, the relatively large value of the expected semimajor axis increase a ˙ induced by the pN gravitoelectric quadrupolar acceleration of Equation (1); let us recall that it is [1,6,32]
d a d t = 9 a n b 3 R e 2 J 2 e 2 6 + e 2 sin 2 I sin 2 ω 8 c 2 1 e 2 4 ,
where n b = μ / a 3 = 2 / P b is the Keplerian mean motion. The need of a highly eccentric orbit is apparent from Equation (6). The frozen perigee would make the resulting integrated pN shift a secular trend. According to Table A3 and Table A7, its predicted rates for the orbital geometries considered in Table A2 and Table A6 are a ˙ = 3 . 8 cm yr 1 and a ˙ = 11 . 6 cm yr 1 , respectively. Suffice it to say that for the existing passive satellites of the LAGEOS family, whose orbits are essentially circular, secular decay rates have been measured over the last decades with an experimental accuracy of a ˙ 3 cm yr 1 for LAGEOS [44,45,46], and a ˙ = 0 . 7 cm yr 1 for LARES [24]. It is arguable that an active mechanism of compensation of the non-gravitational accelerations would allow the increase of such accuracies, allowing, perhaps, to measure the pN quadrupolar effect at a 1 10 % level, depending on the orbital configuration adopted. As far as possible competing effects of gravitational origin are concerned, neither the static and time-dependent parts of the geopotential nor 3rd-body lunisolar attractions induce nonvanishing averaged perturbations on a. Thus, the reduction of the non-conservative accelerations is of the utmost importance. Among them, a prominent role is played by the atmospheric drag, treated in detail in Appendix C, because it causes a net long-term averaged decay rate of a. By modeling the spacecraft as a LARES-like cannonball geodetic satellite for the sake of simplicity, it can be shown that, in the case of the orbital configuration of Table A2, the average acceleration due to the neutral drag only amounts to A drag = ( 2 . 3 0 . 8 ) × 10 11 m s 2 over one orbital period, so that the resulting effect on a would be as large as a ˙ drag = ( 5 . 1 2 . 0 ) m yr 1 ; see Table A5, and Figure A1 and Figure A2. For the more eccentric orbital configuration of Table A6, we have A drag = ( 1 . 22 7 . 3 ) × 10 11 m s 2 and a ˙ drag = ( 27 . 6 164 . 6 ) m yr 1 , as shown by Table A9, and Figure A3 and Figure A4. An inspection of Figure A1 and Figure A3 reveals that, as expected, most of the disturbing effect occurs around the perigee passage, i.e., for f 0 . This may help in suitably calibrating the counteracting action of the drag-free mechanism.
Table A3 and Table A6 show that all the other Keplerian orbital elements also exhibit non-zero secular rates of change due to Equation (1). This is a potentially important feature since they could be linearly combined, as suggested by2 Shapiro [9] in a different context, in order to decouple the pN effect(s) of interest from the disturbing mismodeled Newtonian secular precessions induced by the Earth’s zonal multipoles. Indeed, contrary to a, the other Keplerian orbital elements do exhibit long-term averaged precessions due to the classical zonal harmonics J , = 2 , 3 , 4 , of the geopotential; they are analytically calculated in Appendix B up to degree = 8 . Depending on the specific orbital geometry, they can be both secular and harmonic, or entirely3 secular. To this aim, we complete the set of the pN orbital effects by analytically calculating the averaged rate of the mean anomaly at epoch due to the pN accelerations considered. It turns out that Equations (1) and (2) induce
η ˙ = μ n b R 2 J 2 80 + 73 e 2 1 + 3 cos 2 I 84 1 + 2 e 2 sin 2 I cos 2 ω 32 c 2 a 3 1 e 2 5 / 2 ,
η ˙ = 9 G S R 2 ε 2 5 cos 3 I + cos I 3 + 10 sin 2 I cos 2 ω 56 a 5 c 2 1 e 2 2 ,
respectively, while the gravitoelectric mass monopole acceleration yields
η ˙ = μ n b 15 + 6 1 e 2 + 9 7 1 e 2 ζ c 2 a 1 e 2 ,
where ζ M m / M + m 2 , and m is the satellite’s mass. Instead, it turns out that there is no net Lense-Thirring effect on η . To the benefit of the reader, we review the linear combination approach, which is a generalization of that proposed explicitly for the first time by Ciufolini et al. [47] to test the pN Lense-Thirring effect in the gravitomagnetic field of the Earth with the artificial satellites of the LAGEOS family. It should be noted that, actually, it is quite general, being not necessarily limited just to the pN spin dipole case. By looking at N orbital elements4 κ ( i ) experiencing classical long-term precessions due to the zonals of the geopotential, the following N linear combinations can be written down
μ pN κ ˙ pN ( i ) + s = 2 N κ ˙ J s ( i ) J s δ J s , i = 1 , 2 , N .
They involve the pN averaged precessions κ ˙ pN ( i ) as predicted by general relativity and scaled by a multiplicative parameter μ pN , and the errors in the computed secular node precessions due to the uncertainties in the first N 1 zonals J s , s = 2 , 3 , N , assumed as mismodeled through δ J s , s = 2 , 3 , N . In the following, we will use the shorthand
κ ˙ . κ ˙ J J
for the partial derivative of the classical averaged precession κ ˙ J with respect to the generic even zonal J of degree ; see Appendix B. Then, the N combinations of Equation (10) are posed equal to the experimental residuals δ κ ˙ ( i ) , i = 1 , 2 , N of each of the N orbital elements considered. In principle, such residuals account for the purposely unmodelled pN effect, the mismodelling of the static and time-varying parts of the geopotential, and the non-gravitational forces. Thus, one gets
δ κ ˙ ( i ) = μ pN κ ˙ pN ( i ) + s = 2 N κ ˙ . s ( i ) δ J s , i = 1 , 2 , N .
If we look at the pN scaling parameter5 μ pN and the mismodeling in the first N 1 zonals δ J s , s = 2 , 3 , N as unknowns, we can interpret Equation (12) as an inhomogenous linear system of N algebraic equations in the N unknowns
μ pN , δ J 2 , δ J 3 δ J N N ,
whose coefficients are
κ ˙ pN ( i ) , κ ˙ . s ( i ) , i = 1 , 2 , N , s = 2 , 3 , N ,
while the constant terms are the N orbital residuals
δ κ ˙ ( i ) , i = 1 , 2 , N .
It turns out that, after some algebraic manipulations, the dimensionless pN scaling parameter, which is 1 in general relativity, can be expressed as
μ pN = C δ C pN .
In Equation (16), the combination of the N orbital residuals
C δ δ κ ˙ ( 1 ) + j = 1 N 1 c j δ κ ˙ ( j + 1 )
is, by construction, independent of the first N 1 zonals, being impacted by the other ones of degree > N along with the non-gravitational perturbations and other possible orbital perturbations which cannot be reduced to the same formal expressions of the first N 1 zonal rates. Instead,
C pN κ ˙ pN ( 1 ) + j = 1 N 1 c j κ ˙ pN ( j + 1 )
combines the N pN orbital precessions as predicted by general relativity. The dimensionless coefficients c j , j = 1 , 2 , N 1 in Equations (17) and (18) depend only on some of the orbital parameters of the satellite(s) involved in such a way that, by construction, C δ = 0 if Equation (17) is calculated by posing
δ κ ˙ ( i ) = κ ˙ . ( i ) δ J , i = 1 , 2 , N
for any of the first N 1 zonals, independently of the value assumed for its uncertainty δ J .
As far as HERO is concerned, the linear combination of the four experimental residuals δ Ω , δ η , δ e , δ ω of the satellite’s node, mean anomaly at epoch, eccentricity and perigee suitably designed to cancel out the secular precessions due to the first three zonal harmonics J 2 , J 3 , J 4 of the geopotential is
C δ = δ Ω + c 1 δ η + c 2 δ e + c 3 δ ω .
The coefficients c 1 , c 2 , c 3 turn out to be
c 1 = Ω ˙ . 2 e ˙ . 3 ω ˙ . 4 e ˙ . 2 Ω ˙ . 3 ω ˙ . 4 Ω ˙ . 2 ω ˙ . 3 e ˙ . 4 + ω ˙ . 2 Ω ˙ . 3 e ˙ . 4 + e ˙ . 2 ω ˙ . 3 Ω ˙ . 4 ω ˙ . 2 e ˙ . 3 Ω ˙ . 4 e ˙ . 2 η ˙ . 3 ω ˙ . 4 η ˙ . 2 e ˙ . 3 ω ˙ . 4 e ˙ . 2 ω ˙ . 3 η ˙ . 4 + ω ˙ . 2 e ˙ . 3 η ˙ . 4 + η ˙ . 2 ω ˙ . 3 e ˙ . 4 ω ˙ . 2 η ˙ . 3 e ˙ . 4 ,
c 2 = Ω ˙ . 2 η ˙ . 3 ω ˙ . 4 + η ˙ . 2 Ω ˙ . 3 ω ˙ . 4 + Ω ˙ . 2 ω ˙ . 3 η ˙ . 4 ω ˙ . 2 Ω ˙ . 3 η ˙ . 4 η ˙ . 2 ω ˙ . 3 Ω ˙ . 4 + ω ˙ . 2 η ˙ . 3 Ω ˙ . 4 e ˙ . 2 η ˙ . 3 ω ˙ . 4 η ˙ . 2 e ˙ . 3 ω ˙ . 4 e ˙ . 2 ω ˙ . 3 η ˙ . 4 + ω ˙ . 2 e ˙ . 3 η ˙ . 4 + η ˙ . 2 ω ˙ . 3 e ˙ . 4 ω ˙ . 2 η ˙ . 3 e ˙ . 4 ,
c 3 = Ω ˙ . 2 e ˙ . 3 η ˙ . 4 + e ˙ . 2 Ω ˙ . 3 η ˙ . 4 + Ω ˙ . 2 η ˙ . 3 e ˙ . 4 η ˙ . 2 Ω ˙ . 3 e ˙ . 4 e ˙ . 2 η ˙ . 3 Ω ˙ . 4 + η ˙ . 2 e ˙ . 3 Ω ˙ . 4 e ˙ . 2 η ˙ . 3 ω ˙ . 4 η ˙ . 2 e ˙ . 3 ω ˙ . 4 e ˙ . 2 ω ˙ . 3 η ˙ . 4 + ω ˙ . 2 e ˙ . 3 η ˙ . 4 + η ˙ . 2 ω ˙ . 3 e ˙ . 4 ω ˙ . 2 η ˙ . 3 e ˙ . 4 .
Their numerical values, computed with the formulas of Appendix B for the orbital configurations of Table A2 and Table A6, are listed in Table A4 and Table A8, respectively. In them, the combined mismodeled classical precessions due to the uncancelled zonals, calculated by assuming the formal, statistical sigmas σ J , = 5 , 6 , of the recent global gravity field solution Tongji-Grace02s [48] as a measure of their uncertainties δ J , = 5 , 6 , , are reported as well. However, caution is in order since the realistic level of mismodeling in the geopotential’s coefficients is usually larger than the mere formal errors released in the models produced by various institutions and publicly available on the Internet at http://icgem.gfz-potsdam.de/tom_longtime. A correct evaluation of the actual uncertainties in the zonal harmonics require great care by suitably comparing different global gravity field models; we will not deal with such a task here. From an inspection of Table A4 and Table A8, it can be noted that the (formal) impact of the uncancelled zonals on the combined pN mass quadrupole effect ( G J 2 / c 2 ) is at the 0 . 6 0 . 07 % level for the proposed orbital configurations of Table A2 and Table A6. If the pN spin dipole Lense-Thirring effect ( G S / c 2 ) is considered, the systematic error due to the mismodeling in J , > 4 is about 0 . 1 0 . 03 % . The pN mass monopole combined precessions ( G M / c 2 ) are affected at the ( 20 5 ) × 10 7 relative level. Instead, it turns out that the the combined mismodelled classical precessions are at the same level of the pN spin octupole trends ( G S ε 2 / c 2 ). It may not be unrealistic to expect that, when the forthcoming global gravity field models based on the analysis of the entire long data records of the dedicated GRACE and GOCE missions will be finally available, the current merely formal level of uncertainties in the geopotential’s zonal harmonics may be considered as realistic. Moreover, in the next years, the mission GRACE-FO (GFO) [49], launched in May 2018, will also contribute to the production of new global Earth’s gravity field models of increased quality.
On the other hand, the size of the coefficients c 1 , c 2 , c 3 amplifies the impact of any non-gravitational perturbations that may affect the spacecraft; thus, they should be effectively counteracted by some active drag-free apparatus. In particular, the coefficient c 2 of the eccentricity is ≃30–40; Table A5 and Table A9 show that the expected secular decrease rate of e due to the atmospheric drag is rather large. Thus, some trade-off may be required among the need of reducing the systematic error of gravitational origin and the actual performance of the drag-free mechanism by looking, e.g., at different linear combinations. It may be interesting to note the case of the mean anomaly at epoch. Indeed, in the case of the high perigee orbital configuration of Table A2, Table A3 and Table A5 tell us that the neutral atmospheric drag would represent just ≃1–2% of the predicted pN G J 2 / c 2 precession on η . On the other hand, the present-day formal mismodeling in the classical J 2 -induced rate is about 19 % of it. If the pN Schwarzschild-like effect is considered, the formal bias due to J 2 is at the 10 5 level, while the impact of the atmospheric drag is as little as 10 6 . The neutral atmospheric drag has a larger impact on the pN precessions of η in the case of the low perigee configuration of Table A6, as shown by Table A7 and Table A9.

3. Summary and Overview

The HERO concept—meant as a hopefully drag-free spacecraft moving in a highly eccentric orbit in a frozen perigee configuration that aimed to perform several tests of relativistic gravity in the Earth’s spacetime—represents, in principle, a promising opportunity to measure a general relativistic effect for the first time, which has never received the same attention as the more well-known Schwarzschild and Lense-Thirring precessions. The post-Newtonian gravitoelectric orbital shifts due to the mass quadrupole moment of the Earth. Indeed, the systematic uncertainty in the combined satellite’s precessions due to the formal, statistical errors in the competing Newtonian mass multipoles of the geopotential—as per one of the most recent global gravity field models—is currently below the per cent level for both the orbital configurations proposed. A unique feature of such a post-Newtonian effect is also that the semimajor axis a undergoes a long-term variation which, for a frozen perigee configuration, resembles a secular trend of the order of ≃4–11 cm yr 1 , depending on the orbital geometry chosen. At present, the secular decay of the semimajor axis of the existing passive geodetic satellite LARES has been measured to an accuracy better than 1 cm yr 1 at 2σ level. As far as the traditional Lense-Thirring and Schwarzschild-like post-Newtonian precessions, the formal systematic bias due to the present-day mismodeling in the classical Earth’s zonal harmonics is currently 0 . 1 % and 0 . 0002 % , respectively, if a suitable linear combination of some of the orbital elements of HERO is adopted. However, it must be stressed that the actual uncertainties in the zonal multipoles of the terrestrial gravity field may usually be (much) worse than the sigmas released in the various global gravity solutions. Nonetheless, it cannot be ruled out that, if and when HERO will fly, our knowledge of the Earth’s gravity field will have reached such levels that today’s only formal uncertainties can finally be considered as truly realistic. In addition to the post-Newtonian accelerations, HERO may perform an accurate test of the gravitational red-shift in view of its high eccentricity. Also several models of modified gravity, which generally affect the perigee and the mean anomaly at epoch with secular precessions, could be fruitfully put to the test. A crucial aspect is represented by the level of compensation of the non-gravitational perturbation which will be practically attainable with some drag-free apparatus; suffice it to say that the nominal size of the competing secular decrease of the semimajor axis due to, for example, the neutral atmospheric drag reaching the ≃ 5–160 m yr 1 level if a passive, cannonball satellite is considered. The variability of atmospheric density is such that, in order to achieve drag-free technology, one would have to measure atmospheric density at the satellite in real-time and transmit the data back to Earth for inclusion and analysis within a drag model. Adding such complexities will, however, add to the actual drag as antennas, batteries, solar panels and so forth, and will increase drag and accompanying errors. Perhaps, as the orbit will have to be determined with Satellite Laser Ranging (SLR), one could design and build a type of corner cube reflector where some aspect of the returned SLR pulse is modulated by the atmospheric pressure signal. SLR systems with kHz pulse rates will be useful here. The investigation of such delicate issues deserve dedicated and detailed analyses.

Funding

This research received no external funding

Acknowledgments

I am grateful to an anonymous referee for insightful observations on the atmospheric drag and the drag-free technology.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Tables and Figures

Table A1. Relevant physical parameters of the Earth [50,51,52]. The zonal harmonics of the geopotential of degree are given by J = 2 + 1 C ¯ , 0 , = 2 , 3 , 4 , , where C ¯ , 0 , = 2 , 3 , 4 , are the fully normalized Stokes coefficients of degree and order m = 0 of the multipolar expansion of the Newtonian part of the Earth’s gravity field. The formal, statistical errors in the first seven Stokes coefficients of the geopotential, along with their nominal values, were retrieved from the global gravity field solution Tongji-Grace02s [48] retrievable on the Internet at http://icgem.gfz-potsdam.de/tom_longtime.
Table A1. Relevant physical parameters of the Earth [50,51,52]. The zonal harmonics of the geopotential of degree are given by J = 2 + 1 C ¯ , 0 , = 2 , 3 , 4 , , where C ¯ , 0 , = 2 , 3 , 4 , are the fully normalized Stokes coefficients of degree and order m = 0 of the multipolar expansion of the Newtonian part of the Earth’s gravity field. The formal, statistical errors in the first seven Stokes coefficients of the geopotential, along with their nominal values, were retrieved from the global gravity field solution Tongji-Grace02s [48] retrievable on the Internet at http://icgem.gfz-potsdam.de/tom_longtime.
Physical ParameterNumerical ValueUnits
Newtonian constant of gravitation G 6 . 67259 × 10 11 kg 1 m 3 s 2
Speed of light in vacuum c 2 . 99792458 × 10 8 m s 1
Gravitational parameter μ 3 . 986004418 × 10 14 m 3 s 2
Angular speed Ψ 7 . 29 × 10 5 s 1
Equatorial radius R e 6378 . 1370 km
Polar radius R p 6356 . 7523 km
Angular momentum S 5 . 86 × 10 33 J s
Normalized Stokes coefficient C ¯ 2 , 0 4 . 84165299806 × 10 4 -
Normalized Stokes coefficient C ¯ 3 , 0 9 . 571989759740 × 10 7 -
Normalized Stokes coefficient C ¯ 4 , 0 5 . 399893295930 × 10 7 -
Normalized Stokes coefficient C ¯ 5 , 0 6 . 86499810446677 × 10 8 -
Normalized Stokes coefficient C ¯ 6 , 0 1 . 49976729587105 × 10 7 -
Normalized Stokes coefficient C ¯ 7 , 0 9 . 05017773295824 × 10 8 -
Normalized Stokes coefficient C ¯ 8 , 0 4 . 94794369681244 × 10 8 -
Formal error σ C ¯ 2 , 0 2 . 98340899705584 × 10 13 -
Formal error σ C ¯ 3 , 0 8 . 39284383652709 × 10 14 -
Formal error σ C ¯ 4 , 0 4 . 07426781903578 × 10 14 -
Formal error σ C ¯ 5 , 0 2 . 57688174349872 × 10 14 -
Formal error σ C ¯ 6 , 0 1 . 89009491873398 × 10 14 -
Formal error σ C ¯ 7 , 0 1 . 50081719867797 × 10 14 -
Formal error σ C ¯ 8 , 0 1 . 27528335995664 × 10 14 -
Table A2. Orbital configuration of the proposed satellite HERO: high perigee case. The orbital motion is rather fast since the orbital period P b is as short as 4 h . The period P ω of the perigee is mainly determined by its secular precession due to J 3 , J 4 because of the critical inclination which makes the secular precession due to J 2 nominally vanishing. Note the relatively short period P Ω of the node, amounting to less than 2 yr .
Table A2. Orbital configuration of the proposed satellite HERO: high perigee case. The orbital motion is rather fast since the orbital period P b is as short as 4 h . The period P ω of the perigee is mainly determined by its secular precession due to J 3 , J 4 because of the critical inclination which makes the secular precession due to J 2 nominally vanishing. Note the relatively short period P Ω of the node, amounting to less than 2 yr .
Orbital and Physical ParameterNumerical ValueUnits
Semimajor axis a13,500km
Orbital period P b 4 . 33 hr
Orbital eccentricity e 0 . 45 -
Perigee height h min 1046.86km
Apogee height h max 13,196.9km
Orbital inclination I 63 . 43 deg
Argument of perigee ω 45deg
Period of the node P Ω 1 . 94 yr
Period of the perigee P ω −1363.4yr
Gravitational redshift Δ U c 2 3 . 7 × 10 10
Table A3. Nominal pN (first four rows from the top) and mismodeled Newtonian (first seven rows from the bottom) rates of change, averaged over one orbital revolution, of the semimajor axis a, the eccentricity e, the inclination I, the longitude of the ascending node Ω , the argument of pericenter ω , and the mean anomaly at epoch η for the ideal (no orbital injection error on I assumed) orbital configuration of Table A2. The units are cm yr 1 for a ˙ , and mas yr 1 for e ˙ , I ˙ , Ω ˙ , ω ˙ , η ˙ . The uncertainties in the classical rates of change due to the geopotential are the formal, statistical errors σ J in J , = 2 , 3 , 8 of the model Tongji-Grace02s [48] quoted in Table A1.
Table A3. Nominal pN (first four rows from the top) and mismodeled Newtonian (first seven rows from the bottom) rates of change, averaged over one orbital revolution, of the semimajor axis a, the eccentricity e, the inclination I, the longitude of the ascending node Ω , the argument of pericenter ω , and the mean anomaly at epoch η for the ideal (no orbital injection error on I assumed) orbital configuration of Table A2. The units are cm yr 1 for a ˙ , and mas yr 1 for e ˙ , I ˙ , Ω ˙ , ω ˙ , η ˙ . The uncertainties in the classical rates of change due to the geopotential are the formal, statistical errors σ J in J , = 2 , 3 , 8 of the model Tongji-Grace02s [48] quoted in Table A1.
Effect a ˙ cm yr 1 e ˙ mas yr 1 I ˙ mas yr 1 Ω ˙ mas yr 1 ω ˙ mas yr 1 η ˙ mas yr 1
J 2 c 2 3 . 8 0 . 42 0 . 02 0 . 82 0 . 14 0 . 87
ε 2 c 2 0 0 . 008 0 . 002 0 0 . 074 0 . 015
S c 2 000 32 . 323 43 . 366 0
M c 2 0000 3237 . 8 9292 . 96
σ J 2 000 0 . 411 0 0 . 164
σ J 3 000 0 . 057 0 . 026 0
σ J 4 0 0 . 002 0 . 0006 0 . 034 0 . 049 0 . 004
σ J 5 0 0 . 005 0 . 001 0 . 010 0 . 036 0 . 004
σ J 6 0 0 . 003 0 . 0009 0 . 002 0 . 025 0 . 002
σ J 7 0 0 . 002 0 . 0007 0 . 002 0 . 015 0 . 002
σ J 8 0 0 . 001 0 . 0004 0 . 004 0 . 006 0 . 001
Table A4. Upper three rows: Numerical values of the coefficients c 1 , c 2 , c 3 of the linear combination of Equation (20) canceling out the classical precessions induced by the first three zonal harmonics J 2 , J 3 , J 4 for the orbital configuration of Table A2. Middle rows: Uncancelled mismodeled precessions due to the zonal harmonics J 5 , J 6 , J 7 , J 8 , in mas yr 1 , linearly combined according to Equation (20); the formal, statistical errors σ J in J , = 5 , 6 , 7 , 8 of the model Tongji-Grace02s [48], quoted in Table A1 were used. Lower four rows: pN precessions, in mas yr 1 , linearly combined according to Equation (20).
Table A4. Upper three rows: Numerical values of the coefficients c 1 , c 2 , c 3 of the linear combination of Equation (20) canceling out the classical precessions induced by the first three zonal harmonics J 2 , J 3 , J 4 for the orbital configuration of Table A2. Middle rows: Uncancelled mismodeled precessions due to the zonal harmonics J 5 , J 6 , J 7 , J 8 , in mas yr 1 , linearly combined according to Equation (20); the formal, statistical errors σ J in J , = 5 , 6 , 7 , 8 of the model Tongji-Grace02s [48], quoted in Table A1 were used. Lower four rows: pN precessions, in mas yr 1 , linearly combined according to Equation (20).
Coefficient of δ η c 1 2 . 51065
Coefficient of δ e c 2 29 . 0889
Coefficient of δ ω c 3 2 . 13813
σ J 5 formal 0 . 06 mas yr 1
σ J 6 formal 0 . 03 mas yr 1
σ J 7 formal 0 . 03 mas yr 1
σ J 8 formal 0 . 02 mas yr 1
J 2 c 2 10 . 75 mas yr 1
ε c 2 0 . 03 mas yr 1
S c 2 60 . 07 mas yr 1
M c 2 30,254.2 mas yr 1
Table A5. Numerically integrated nominal rates of change, averaged over one orbital revolution, of the semimajor axis a, the eccentricity e, the inclination I, the longitude of the ascending node Ω , the argument of pericenter ω , and the mean anomaly at epoch η induced by the neutral atmospheric drag for the orbital configuration of Table A2. The units are m yr 1 for a ˙ , and mas yr 1 for e ˙ , I ˙ , Ω ˙ , ω ˙ , η ˙ . For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s neutral atmospheric density, we adopted r 0 = r min = a 1 e = 1046 . 86 km , ρ 0 = ρ max = 7 . 3 2 . 8 × 10 15 kg m 3 , ρ min = 0 . 001 ρ L , λ = 872 . 87 km 938 . 49 km ; the neutral atmospheric density at the height of LAGEOS is ρ L = 6 . 579 × 10 18 kg m 3 [54]. See Appendix C for details. Neither approximations in e nor in ν Ψ / n b were used. The value of ρ 0 = ρ max was kept fixed over one orbital revolution. Cfr. with the analytical plots in Figure A1 and the numerically produced time series in Figure A2.
Table A5. Numerically integrated nominal rates of change, averaged over one orbital revolution, of the semimajor axis a, the eccentricity e, the inclination I, the longitude of the ascending node Ω , the argument of pericenter ω , and the mean anomaly at epoch η induced by the neutral atmospheric drag for the orbital configuration of Table A2. The units are m yr 1 for a ˙ , and mas yr 1 for e ˙ , I ˙ , Ω ˙ , ω ˙ , η ˙ . For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s neutral atmospheric density, we adopted r 0 = r min = a 1 e = 1046 . 86 km , ρ 0 = ρ max = 7 . 3 2 . 8 × 10 15 kg m 3 , ρ min = 0 . 001 ρ L , λ = 872 . 87 km 938 . 49 km ; the neutral atmospheric density at the height of LAGEOS is ρ L = 6 . 579 × 10 18 kg m 3 [54]. See Appendix C for details. Neither approximations in e nor in ν Ψ / n b were used. The value of ρ 0 = ρ max was kept fixed over one orbital revolution. Cfr. with the analytical plots in Figure A1 and the numerically produced time series in Figure A2.
ρ 0 kg m 3 a ˙ m yr 1 e ˙ mas yr 1 I ˙ mas yr 1 Ω ˙ mas yr 1 ω ˙ mas yr 1 η ˙ mas yr 1
7 . 3 × 10 15 5 . 1 41 0 . 51 0 . 21 0 . 12 0 . 02
2 . 8 × 10 15 2 16 0 . 2 0 . 07 0 . 04 0 . 01
Table A6. Orbital configuration of the proposed satellite HERO: low perigee case. The orbital motion is relatively slow since the orbital period P b is as long as more than 21 h . The period P ω of the perigee is mainly determined by its secular precession due to J 3 , J 4 because of the critical inclination which makes the secular precession due to J 2 nominally vanishing. The period P Ω of the node is rather long, amounting to more than 13 yr .
Table A6. Orbital configuration of the proposed satellite HERO: low perigee case. The orbital motion is relatively slow since the orbital period P b is as long as more than 21 h . The period P ω of the perigee is mainly determined by its secular precession due to J 3 , J 4 because of the critical inclination which makes the secular precession due to J 2 nominally vanishing. The period P Ω of the node is rather long, amounting to more than 13 yr .
Orbital and Physical ParameterNumerical ValueUnits
Semimajor axis a39,000km
Orbital period P b 21 . 29 hr
Orbital eccentricity e 0 . 82 -
Perigee height h min 641 . 86 km
Apogee height h max 64,601.9km
Orbital inclination I 63 . 43 deg
Argument of perigee ω 45deg
Period of the node P Ω 13 . 45 yr
Period of the perigee P ω −8186.71yr
Gravitational redshift Δ U c 2 5 . 7 × 10 10
Table A7. Nominal pN (first four rows from the top) and mismodeled Newtonian (first seven rows from the bottom) rates of change, averaged over one orbital revolution, of the semimajor axis a, the eccentricity e, the inclination I, the longitude of the ascending node Ω , the argument of pericenter ω , and the mean anomaly at epoch η for the ideal (no orbital injection error on I assumed) orbital configuration of Table A6. The units are cm yr 1 for a ˙ , and mas yr 1 for e ˙ , I ˙ , Ω ˙ , ω ˙ , η ˙ . The uncertainties in the classical rates of change due to the geopotential are the formal, statistical errors σ J in J , = 2 , 3 , 8 of the model Tongji-Grace02s [48] quoted in Table A1.
Table A7. Nominal pN (first four rows from the top) and mismodeled Newtonian (first seven rows from the bottom) rates of change, averaged over one orbital revolution, of the semimajor axis a, the eccentricity e, the inclination I, the longitude of the ascending node Ω , the argument of pericenter ω , and the mean anomaly at epoch η for the ideal (no orbital injection error on I assumed) orbital configuration of Table A6. The units are cm yr 1 for a ˙ , and mas yr 1 for e ˙ , I ˙ , Ω ˙ , ω ˙ , η ˙ . The uncertainties in the classical rates of change due to the geopotential are the formal, statistical errors σ J in J , = 2 , 3 , 8 of the model Tongji-Grace02s [48] quoted in Table A1.
Effect a ˙ cm yr 1 e ˙ mas yr 1 I ˙ mas yr 1 Ω ˙ mas yr 1 ω ˙ mas yr 1 η ˙ mas yr 1
J 2 c 2 11 . 6 0 . 115 0 . 010 0 . 100 0 . 022 0 . 092
ε 2 c 2 0 0 . 0006 0 . 0008 0 0 . 0106 0 . 0004
S c 2 000 5 . 09 6 . 83 0
M c 2 0000 555 . 661 1226 . 13
σ J 2 000 0 . 059 0 0 . 015
σ J 3 000 0 . 0128 0 . 006 0
σ J 4 0 0 . 0001 0 . 0002 0 . 005 0 . 007 0 . 0009
σ J 5 0 0 . 0002 0 . 0003 0 . 002 0 . 005 0 . 0006
σ J 6 0 0 . 0002 0 . 0002 0 . 0002 0 . 003 0 . 0003
σ J 7 0 0 . 0001 0 . 0002 0 . 0005 0 . 002 0 . 0002
σ J 8 0 0 . 00008 0 . 0001 0 . 0008 0 . 0007 0 . 00007
Table A8. Upper three rows: Numerical values of the coefficients c 1 , c 2 , c 3 of the linear combination of Equation (20) canceling out the classical precessions induced by the first three zonal harmonics J 2 , J 3 , J 4 for the orbital configuration of Table A6. Middle four rows: Uncancelled mismodeled precessions due to the zonal harmonics J 5 , J 6 , J 7 , J 8 , in mas yr 1 , linearly combined according to Equation (20); the formal, statistical errors σ J in J , = 5 , 6 , 7 , 8 of the model Tongji-Grace02s [48], quoted in Table A1, were used. Lower four rows: pN precessions, in mas yr 1 , linearly combined according to Equation (20).
Table A8. Upper three rows: Numerical values of the coefficients c 1 , c 2 , c 3 of the linear combination of Equation (20) canceling out the classical precessions induced by the first three zonal harmonics J 2 , J 3 , J 4 for the orbital configuration of Table A6. Middle four rows: Uncancelled mismodeled precessions due to the zonal harmonics J 5 , J 6 , J 7 , J 8 , in mas yr 1 , linearly combined according to Equation (20); the formal, statistical errors σ J in J , = 5 , 6 , 7 , 8 of the model Tongji-Grace02s [48], quoted in Table A1, were used. Lower four rows: pN precessions, in mas yr 1 , linearly combined according to Equation (20).
Coefficient of δ η c 1 3 . 91939
Coefficient of δ e c 2 40 . 7154
Coefficient of δ ω c 3 2 . 20981
σ J 5 0 . 003 mas yr 1
σ J 6 0 . 002 mas yr 1
σ J 7 0 . 002 mas yr 1
σ J 8 0 . 001 mas yr 1
J 2 c 2 4 . 373 mas yr 1
ε c 2 0 . 001 mas yr 1
S c 2 9 . 952 mas yr 1
M c 2 6033 . 6 mas yr 1
Table A9. Numerically integrated nominal rates of change, averaged over one orbital revolution, of the semimajor axis a, the eccentricity e, the inclination I, the longitude of the ascending node Ω , the argument of pericenter ω , and the mean anomaly at epoch η induced by the neutral atmospheric drag for the orbital configuration of Table A6. The units are m yr 1 for a ˙ , and mas yr 1 for e ˙ , I ˙ , Ω ˙ , ω ˙ , η ˙ . For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s neutral atmospheric density, we adopted r 0 = r min = a 1 e = 641 . 86 km , ρ 0 = ρ max = 6 . 9 1 . 11 × 10 14 kg m 3 , ρ min = 0 . 0001 ρ L , λ = 3463 . 23 km 3843 . 48 km ; the neutral atmospheric density at the height of LAGEOS is ρ L = 6 . 579 × 10 18 kg m 3 [54]. See Appendix C for details. Neither approximations in e nor in ν Ψ / n b were used. The value of ρ 0 = ρ max was kept fixed over one orbital revolution. Cfr. with the analytical plots in Figure A3 and the numerically produced time series in Figure A4.
Table A9. Numerically integrated nominal rates of change, averaged over one orbital revolution, of the semimajor axis a, the eccentricity e, the inclination I, the longitude of the ascending node Ω , the argument of pericenter ω , and the mean anomaly at epoch η induced by the neutral atmospheric drag for the orbital configuration of Table A6. The units are m yr 1 for a ˙ , and mas yr 1 for e ˙ , I ˙ , Ω ˙ , ω ˙ , η ˙ . For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s neutral atmospheric density, we adopted r 0 = r min = a 1 e = 641 . 86 km , ρ 0 = ρ max = 6 . 9 1 . 11 × 10 14 kg m 3 , ρ min = 0 . 0001 ρ L , λ = 3463 . 23 km 3843 . 48 km ; the neutral atmospheric density at the height of LAGEOS is ρ L = 6 . 579 × 10 18 kg m 3 [54]. See Appendix C for details. Neither approximations in e nor in ν Ψ / n b were used. The value of ρ 0 = ρ max was kept fixed over one orbital revolution. Cfr. with the analytical plots in Figure A3 and the numerically produced time series in Figure A4.
ρ 0 kg m 3 a ˙ m yr 1 e ˙ mas yr 1 I ˙ mas yr 1 Ω ˙ mas yr 1 ω ˙ mas yr 1 η ˙ mas yr 1
6 . 9 × 10 14 164 . 65 152 . 96 2 . 24 0 . 69 0 . 30 0 . 02
1 . 11 × 10 14 27 . 6 25 . 6 0 . 41 0 . 15 0 . 07 0 . 008
Figure A1. Plots of Equations (A49) and (A54) as functions of the true anomaly f from 0 to 2π for the orbital configuration of Table A2. For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s atmospheric density, we adopted r 0 = r min = a 1 e = 1046 . 86 km , ρ 0 = ρ max = 7 . 3 × 10 15 kg m 3 , λ = 872 . 87 km . Neither approximations in e nor in ν Ψ / n b were used. The areas of the regions delimited by the curves and the f axis are the rates of change of the orbital elements averaged over one orbital period P b ; they are numerically calculated and displayed in Table A5. The value of ρ 0 = ρ max was kept fixed over one orbital revolution.
Figure A1. Plots of Equations (A49) and (A54) as functions of the true anomaly f from 0 to 2π for the orbital configuration of Table A2. For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s atmospheric density, we adopted r 0 = r min = a 1 e = 1046 . 86 km , ρ 0 = ρ max = 7 . 3 × 10 15 kg m 3 , λ = 872 . 87 km . Neither approximations in e nor in ν Ψ / n b were used. The areas of the regions delimited by the curves and the f axis are the rates of change of the orbital elements averaged over one orbital period P b ; they are numerically calculated and displayed in Table A5. The value of ρ 0 = ρ max was kept fixed over one orbital revolution.
Universe 05 00165 g0a1
Figure A2. Numerically integrated shifts of the semimajor axis a, the eccentricity e, the longitude of the ascending node Ω , the argument of the perigee ω , and the mean anomaly at epoch η for the orbital configuration of Table A2. For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s atmospheric density, we adopted r 0 = r min = a 1 e = 1046 . 86 km , ρ 0 = ρ max = 7 . 3 × 10 15 kg m 3 , λ = 872 . 87 km . Cfr. with the semi-analytical results of Figure A1 and Table A5.
Figure A2. Numerically integrated shifts of the semimajor axis a, the eccentricity e, the longitude of the ascending node Ω , the argument of the perigee ω , and the mean anomaly at epoch η for the orbital configuration of Table A2. For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s atmospheric density, we adopted r 0 = r min = a 1 e = 1046 . 86 km , ρ 0 = ρ max = 7 . 3 × 10 15 kg m 3 , λ = 872 . 87 km . Cfr. with the semi-analytical results of Figure A1 and Table A5.
Universe 05 00165 g0a2
Figure A3. Plots of Equations (A49) and (A54) as functions of the true anomaly f from 0 to 2π for the orbital configuration of Table A6. For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s atmospheric density, we adopted r 0 = r min = a 1 e = 641 . 86 km , ρ 0 = ρ max = 6 . 9 × 10 14 kg m 3 , λ = 3463 . 23 km . Neither approximations in e nor in ν Ψ / n b were used. The areas of the regions delimited by the curves and the f axis are the rates of change of the orbital elements averaged over one orbital period P b ; they are numerically calculated and displayed in Table A9. The value of ρ 0 = ρ max was kept fixed over one orbital revolution.
Figure A3. Plots of Equations (A49) and (A54) as functions of the true anomaly f from 0 to 2π for the orbital configuration of Table A6. For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s atmospheric density, we adopted r 0 = r min = a 1 e = 641 . 86 km , ρ 0 = ρ max = 6 . 9 × 10 14 kg m 3 , λ = 3463 . 23 km . Neither approximations in e nor in ν Ψ / n b were used. The areas of the regions delimited by the curves and the f axis are the rates of change of the orbital elements averaged over one orbital period P b ; they are numerically calculated and displayed in Table A9. The value of ρ 0 = ρ max was kept fixed over one orbital revolution.
Universe 05 00165 g0a3
Figure A4. Numerically integrated shifts of the semimajor axis a, the eccentricity e, the longitude of the ascending node Ω , the argument of the perigee ω , and the mean anomaly at epoch η for the orbital configuration of Table A6. For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s atmospheric density, we adopted r 0 = r min = a 1 e = 641 . 86 km , ρ 0 = ρ max = 6 . 9 × 10 14 kg m 3 , λ = 3463 . 23 km . Cfr. with the semi-analytical results of Figure A3 and Table A9.
Figure A4. Numerically integrated shifts of the semimajor axis a, the eccentricity e, the longitude of the ascending node Ω , the argument of the perigee ω , and the mean anomaly at epoch η for the orbital configuration of Table A6. For the satellite, assumed spherical in shape and passive, we adopted C D = 3 . 5 , Σ = 2 . 69 × 10 4 as for the existing LARES [53]. In regard to the Earth’s atmospheric density, we adopted r 0 = r min = a 1 e = 641 . 86 km , ρ 0 = ρ max = 6 . 9 × 10 14 kg m 3 , λ = 3463 . 23 km . Cfr. with the semi-analytical results of Figure A3 and Table A9.
Universe 05 00165 g0a4

Appendix B. Classical Long-Term Rates of Change of the Keplerian Orbital Elements up to Degree = 8

Here, we will analytically work out the coefficients κ ˙ . κ ˙ J / J , κ = e , I , Ω , ω , η of the long-term rates of change of the eccentricity e, the inclination I, the longitude of the ascending node Ω , the argument of perigee ω , and the mean anomaly at epoch η induced by the first seven zonal harmonics of the geopotential up to degree = 8 ; the averaged rates of change of the semimajor axis a are all vanishing. We will adopt the Lagrange planetary equations [55,56,57,58] applied to the perturbing potential of degree
Δ U r = μ r R e r J P ξ , = 2 , 3 , 8
where P ξ is the Legendre polynomial of degree , averaged over one orbital period as disturbing function. We will not make any a-priory assumption on the orbital configuration of the satellite. As far as the orientation of the Earth’s spin axis, we will align it to the reference z axis of an equatorial coordinate system. The following formulas include both the genuine secular and the harmonic components with the frequency of the perigee and its integer multiples.

Appendix B.1. The Eccentricity

e ˙ . 2 = 0 ,
e ˙ . 3 = 3 n b R e 3 3 + 5 cos 2 I sin I cos ω 16 a 3 1 e 2 2 ,
e ˙ . 4 = 15 e n b R e 4 5 + 7 cos 2 I sin 2 I sin 2 ω 64 a 4 1 e 2 3 ,
e ˙ . 5 = 15 n b R e 5 2048 a 5 1 e 2 4 × × 28 e 2 7 + 9 cos 2 I sin 3 I cos 3 ω + 4 + 3 e 2 2 sin I + 7 sin 3 I + 3 sin 5 I cos ω ,
e ˙ . 6 = 105 e n b R e 6 sin 2 I 8192 a 6 1 e 2 5 × × 5 2 + e 2 35 + 60 cos 2 I + 33 cos 4 I sin 2 ω + 12 e 2 9 + 11 cos 2 I sin 2 I sin 4 ω ,
e ˙ . 7 = 21 n b R e 7 524288 a 7 1 e 2 6 × × 60 e 2 sin 3 I 3 8 + 3 e 2 189 + 308 cos 2 I + 143 cos 4 I cos 3 ω 44 e 2 11 + 13 cos 2 I sin 2 I cos 5 ω 5 8 + 5 e 2 4 + e 2 25 sin I + 81 sin 3 I + 165 sin 5 I + 429 sin 7 I cos ω ,
e ˙ . 8 = 63 e n b R e 8 sin 2 I 2097152 a 8 1 e 2 7 × × 35 48 + 80 e 2 + 15 e 4 210 + 385 cos 2 I + 286 cos 4 I + 143 cos 6 I sin 2 ω + + 88 e 2 sin 2 I 7 10 + 3 e 2 99 + 156 cos 2 I + 65 cos 4 I sin 4 ω + + 26 e 2 13 + 15 cos 2 I sin 2 I sin 6 ω .

Appendix B.2. The Inclination

I ˙ . 2 = 0 ,
I ˙ . 3 = 3 e n b R e 3 cos I 3 + 5 cos 2 I cos ω 16 a 3 1 e 2 3 ,
I ˙ . 4 = 15 e 2 n b R e 4 cos I 5 + 7 cos 2 I sin I sin 2 ω 64 a 4 1 e 2 4 ,
I ˙ . 5 = 15 e n b R e 5 cot I 2048 a 5 1 e 2 5 × × 28 e 2 7 + 9 cos 2 I sin 3 I cos 3 ω + 4 + 3 e 2 2 sin I + 7 sin 3 I + 3 sin 5 I cos ω ,
I ˙ . 6 = 5 n b R e 6 cot I 65536 a 6 1 e 2 6 × × 840 e 2 2 + e 2 35 + 60 cos 2 I + 33 cos 4 I sin 2 I sin 2 ω 2016 e 4 9 + 11 cos 2 I sin 4 I sin 4 ω ,
I ˙ . 7 = 21 e n b R e 7 cot I 524288 a 7 1 e 2 7 × × 60 e 2 sin 3 I 3 8 + 3 e 2 189 + 308 cos 2 I + 143 cos 4 I cos 3 ω 44 e 2 11 + 13 cos 2 I sin 2 I cos 5 ω 5 8 + 5 e 2 4 + e 2 25 sin I + 81 sin 3 I + 165 sin 5 I + 429 sin 7 I cos ω ,
I ˙ . 8 = 63 e 2 n b R e 8 sin I cos I 2097152 a 8 1 e 2 8 × × 35 48 + 80 e 2 + 15 e 4 210 + 385 cos 2 I + 286 cos 4 I + 143 cos 6 I sin 2 ω + + 88 e 2 sin 2 I 7 10 + 3 e 2 99 + 156 cos 2 I + 65 cos 4 I sin 4 ω + + 26 e 2 13 + 15 cos 2 I sin 2 I sin 6 ω .

Appendix B.3. The Longitude of the Ascending Node

Ω ˙ . 2 = 3 n b R e 2 cos I 2 a 2 1 e 2 2 ,
Ω ˙ . 3 = 3 e n b R e 3 cos I + 15 cos 3 I csc I sin ω 32 a 3 1 e 2 3 ,
Ω ˙ . 4 = 15 n b R e 4 2 + 3 e 2 9 cos I + 7 cos 3 I 2 e 2 5 cos I + 7 cos 3 I cos 2 ω 128 a 4 1 e 2 4 ,
Ω ˙ . 5 = 15 e n b R e 5 2048 a 5 1 e 2 5 × × 4 + 3 e 2 2 cos I + 21 cos 3 I + 5 cos 5 I csc I sin ω + 7 e 2 2 sin 2 I + 15 sin 4 I sin 3 ω ,
Ω ˙ . 6 = 105 n b R e 6 16384 a 6 1 e 2 6 × × 8 + 40 e 2 + 15 e 4 50 cos I + 45 cos 3 I + 33 cos 5 I + + 5 e 2 2 + e 2 70 cos I + 87 cos 3 I + 99 cos 5 I cos 2 ω + + 6 e 4 47 cos I + 33 cos 3 I sin 2 I cos 4 ω ,
Ω ˙ . 7 = 21 e n b R e 7 csc I 524288 a 7 1 e 2 7 × × 5 8 + 5 e 2 4 + e 2 25 cos I + 243 cos 3 I + 825 cos 5 I + 3003 cos 7 I sin ω 30 e 2 8 + 3 e 2 1442 cos I + 1397 cos 3 I + 1001 cos 5 I sin 2 I sin 3 ω 264 e 4 149 cos I + 91 cos 3 I sin 4 I sin 5 ω ,
Ω ˙ . 8 = 63 n b R e 8 2097152 a 8 1 e 2 8 × × 5 16 + 7 e 2 24 + 5 e 2 6 + e 2 1225 cos I + 11 105 cos 3 I + 91 cos 5 I + 65 cos 7 I 70 e 2 48 + 80 e 2 + 15 e 4 105 cos I + 11 11 cos 3 I + 13 cos 5 I + cos 7 I cos 2 ω 616 e 4 10 + 3 e 2 138 cos I + 117 cos 3 I + 65 cos 5 I sin 2 I cos 4 ω 9152 e 6 cos I 2 + 5 cos 2 I sin 4 I cos 6 ω .

Appendix B.4. The Argument of Perigee

ω ˙ . 2 = 3 n b R e 2 3 + 5 cos 2 I 8 a 2 1 e 2 2 ,
ω ˙ . 3 = 3 n b R e 3 64 a 3 e 1 e 2 3 1 3 e 2 4 cos 2 I + 5 1 + 7 e 2 cos 4 I csc I sin ω ,
ω ˙ . 4 = 15 n b R e 4 1024 a 4 1 e 2 4 × × 27 4 + 5 e 2 + 4 cos 2 I 52 63 e 2 + 2 2 + 7 e 2 cos 2 ω + 2 6 + 5 e 2 cos 2 ω + + 7 cos 4 I 28 27 e 2 + 2 2 + 9 e 2 cos 2 ω ,
ω ˙ . 5 = 15 n b R e 5 sin I 4096 a 5 e 1 e 2 5 × × 8 + 74 e 2 + 30 e 4 + 20 + 113 e 2 + 21 e 4 cos 2 I + 14 4 + 5 e 2 9 e 4 cos 4 I 21 4 + 61 e 2 + 33 e 4 cos 6 I csc 2 I sin ω 14 e 2 5 + 7 e 2 + 4 1 + 6 e 2 cos 2 I + 9 + 33 e 2 cos 4 I sin 3 ω ,
ω ˙ . 6 = 105 n b R e 6 65536 a 6 1 e 2 6 × × 5 472 + 1 , 940 e 2 + 675 e 4 cos 2 I + + 3 e 2 2 292 + 99 e 2 cos 4 I + 11 44 + 13 e 2 cos 6 I 5 10 e 2 6 + 7 e 2 + 68 + 254 e 2 + 195 e 4 cos 2 I + 6 4 + 102 e 2 + 55 e 4 cos 4 I + + 33 4 + 34 e 2 + 13 e 4 cos 6 I cos 2 ω + + 2 1128 cos 4 I + 1188 cos 6 I + 25 24 + 100 e 2 + 35 e 4 + 4 cos 2 ω 6 e 2 28 + 45 e 2 + 4 4 + 33 e 2 cos 2 I + 11 4 + 13 e 2 cos 4 I sin 2 I cos 4 ω ,
ω ˙ . 7 = 21 n b R e 7 524288 e a 7 1 e 2 7 × × 5 8 + 156 e 2 + 225 e 4 + 40 e 6 25 sin I + 81 sin 3 I + 165 sin 5 I + 429 sin 7 I sin ω 60 e 2 24 + 95 e 2 + 24 e 4 189 + 308 cos 2 I + 143 cos 4 I sin 3 I sin 3 ω + + 10 e 2 cos I 8 + 5 e 2 4 + e 2 1198 + 2421 cos 2 I 2178 cos 4 I + + 3003 cos 6 I cot I sin ω + + 3 e 2 8 + 3 e 2 523 + 396 cos 2 I + 1001 cos 4 I sin 2 I sin 3 ω + + 528 e 6 cos 2 I 29 + 91 cos 2 I sin 3 I sin 5 ω 528 e 4 5 + 8 e 2 11 + 13 cos 2 I sin 5 I sin 5 ω ,
ω ˙ . 8 = 63 n b R e 8 33554432 a 8 1 e 2 8 × × 5 1225 192 + 35 e 2 48 + 56 e 2 + 9 e 4 280 1664 + 7 e 2 2064 + 2400 e 2 + 385 e 4 cos 2 I 308 1472 + 7 e 2 1776 + 2040 e 2 + 325 e 4 cos 4 I 3432 128 + 7 e 2 144 + 160 e 2 + 25 e 4 cos 6 I 715 704 + 7 e 2 624 + 600 e 2 + 85 e 4 cos 8 I + + 70 35 96 + 208 e 2 + 950 e 4 + 225 e 6 + + 16 384 + 1648 e 2 + 5160 e 4 + 1155 e 6 cos 2 I + + 44 96 + 1360 e 2 + 2870 e 4 + 585 e 6 cos 4 I + + 2288 e 2 48 + 80 e 2 + 15 e 4 cos 6 I + + 143 96 + 1328 e 2 + 1610 e 4 + 255 e 6 cos 8 I cos 2 ω + + 616 e 2 6 280 + 944 e 2 + 363 e 4 + 1960 + 14128 e 2 + 4797 e 4 cos 2 I + + 26 40 + 688 e 2 + 195 e 4 cos 4 I + 65 40 + 208 e 2 + 51 e 4 cos 6 I sin 2 I cos 4 ω + + 4576 e 4 22 + 39 e 2 + 4 2 + 25 e 2 cos 2 I + 5 6 + 17 e 2 cos 4 I sin 4 I cos 6 ω .

Appendix B.5. The Mean Anomaly at Epoch

η ˙ . 2 = 3 n b R e 2 1 + 3 cos 2 I 8 a 2 1 e 2 3 / 2 ,
η ˙ . 3 = 3 n b R e 3 16 e a 3 1 e 2 5 / 2 1 + 4 e 2 3 + 5 cos 2 I sin I sin ω ,
η ˙ . 4 = 15 n b R e 4 1024 a 4 1 e 2 7 / 2 × × 3 e 2 9 + 20 cos 2 I + 35 cos 4 I + 8 2 + 5 e 2 5 + 7 cos 2 I sin 2 I cos 2 ω ,
η ˙ . 5 = 15 n b R e 5 2048 e a 5 1 e 2 9 / 2 × × 4 + 7 e 2 + 18 e 4 2 sin I + 7 sin 3 I + 3 sin 5 I sin ω 28 e 2 1 + 2 e 2 7 + 9 cos 2 I sin 3 I sin 3 ω ,
η ˙ . 6 = 35 n b R e 6 65536 a 6 1 e 2 11 / 2 × × 8 + 20 e 2 + 15 e 4 50 + 105 cos 2 I + 126 cos 4 I + 231 cos 6 I 60 4 + 6 e 2 + 7 e 4 35 + 60 cos 2 I + 33 cos 4 I sin 2 I cos 2 ω 72 e 2 4 + 7 e 2 9 + 11 cos 2 I sin 4 I cos 4 ω ,
η ˙ . 7 = 21 n b R e 7 524288 e a 7 1 e 2 13 / 2 × × 5 8 28 e 2 + 95 e 4 + 40 e 6 25 sin I + 81 sin 3 I + 165 sin 5 I + 429 sin 7 I sin ω + + 180 e 2 8 + 11 e 2 + 8 e 4 189 + 308 cos 2 I + 143 cos 4 I sin 3 I sin 3 ω + + 528 e 4 5 + 8 e 2 11 + 13 cos 2 I sin 5 I sin 5 ω ,
η ˙ . 8 = 63 n b R e 8 33554432 a 8 1 e 2 15 / 2 × × 5 32 + 35 e 4 4 + e 2 1225 + 2520 cos 2 I + 2772 cos 4 I + 3432 cos 6 I + + 6435 cos 8 I + 280 96 80 e 2 + 470 e 4 + 135 e 6 210 + 385 cos 2 I + 286 cos 4 I + + 143 cos 6 I sin 2 I cos 2 ω + + 2464 e 2 40 + 52 e 2 + 27 e 4 99 + 156 cos 2 I + 65 cos 4 I sin 4 I cos 4 ω + + 18304 e 4 2 + 3 e 2 13 + 15 cos 2 I sin 6 I cos 6 ω .

Appendix C. The Atmospheric Drag

The neutral and charged atmospheric drag is potentially a major source of systematic uncertainty since it induces long-term effects on all the satellite’s orbital elements which, for a frozen perigee configuration, may look like secular trends on which the time variability of the atmospheric density is superimposed.
For the sake of simplicity, we will model HERO as a cannonball, passive satellite with the same physical properties of the existing LARES satellite in order to make some quantitative estimates of the disturbing impact of the atmospheric drag on the pN effects of interest. Thus, its perturbing acceleration is customarily modeled as
A D = 1 2 C D Σ ρ V V ,
where C D , Σ , ρ , V are the dimensionless drag coefficient of the spacecraft, its area-to-mass ratio, the atmospheric density at its height, and its velocity with respect to the atmosphere, respectively. By assuming that the atmosphere co-rotates with the Earth, V can be posed equal to
V = v Ψ × r ,
where Ψ is the Earth’s angular velocity. In fact, a decrease of the co-rotation with the height is expected. Membrado & Pacheco [59] modeled it in two scenarios involving a constant and non-constant viscosity. In the first case, the second term in Equation (A38) must be rescaled by R e / r 3 . In general, the atmospheric density may not be considered as constant throughout a highly eccentric orbit covering a wide range of geocentric distances, as in our case; thus, we will model it as
ρ ( r ) = ρ 0 exp r r 0 λ ,
where ρ 0 is the atmospheric density at some reference distance r 0 , while λ is the characteristic scale length. By assuming
r 0 = r min = a 1 e ,
if the atmospheric density is known at the perigee and apogee heights, λ can be determined as
λ = 2 a e ln ρ min ρ max ,
where
ρ min = ρ ( r max ) ,
ρ max = ρ ( r min ) .
In the case of the orbital configuration of Table A2, the perigee height is h min = 1046 . 86 km ; we will determine the corresponding neutral atmospheric density ρ 0 = ρ max as follows. The neutral atmospheric density at the LARES height, which is h LR = 1450 km , amounts to ρ LR = 5 . 644 × 10 16 kg m 3 [53]. According to Brito, Celestino & Moraes [60], the neutral atmospheric density at h 700 = 700 km is 6 . 9 × 10 14 kg m 3 (TD-88 model), or 1 . 11 × 10 14 kg m 3 (NASA model). Thus, it is possible to calculate the characteristic scale length λ valid for the range 700 km < h < 1450 km as
λ 700 / LR = h 700 h LR ln ρ LR ρ 700 = 154 . 422 km 249 . 139 km ,
depending on the value adopted for ρ at 700 km . Since for the orbital configuration of HERO of Table A2 it is h 700 < h min < h LR , one can determine ρ max by using just Equation (A44) in
ρ max = ρ LR exp h LR h min λ 700 / LR = 7 . 3 2 . 8 × 10 15 kg m 3 ,
depending on the value of λ 700 / LR adopted. As far as ρ min is concerned, since h max = 13,169.9 km is much larger than the height of, say, the LAGEOS satellite ( h L = 5891 . 87 km ), for which it is ρ L = 6 . 579 × 10 18 kg m 3 [54], it does not seem unreasonable to assume
ρ min 0 . 001 ρ L ,
or so. Thus, Equation (A41), applied to ρ max , given by Equation (A45), and ρ min , given by Equation (A46), allows to infer the value for λ which must be used in the calculation of the drag effect for HERO. It is
λ = 872 . 87 km 938 . 49 km ,
depending on the value of ρ max adopted. As far as the more eccentric orbital configuration of Table A6 is concerned, by adopting for ρ 0 the two possible values of ρ 700 and, say, ρ min = 0 . 0001 ρ L , it turns out that
λ = 3463 . 23 km 3843 . 48 km .
Actually, even the density at a given height may not be regarded as truly constant because of a variety of geophysical phenomena characterized by quite different time scales. Anyway, in order to have an order-of-magnitude evaluation of the perturbing action of Equation (A37) on the motion of HERO, we calculate the averaged rates of change of its Keplerian orbital elements by keeping ρ 0 = ρ max in Equation (A39) fixed during one orbital period P b ; given its short duration, at least in the case of Table A2, it is not an unreasonable assumption. The large value of the eccentricity makes most of the existing results in the literature unsuitable to the present case; moreover, an exact analytical calculation without recurring to any approximation in both e and ν Ψ / n b is difficult. Thus, we follow two complementary approaches. In one of them, we, first, plot in Figure A1 and Figure A3 the analytical expressions of the ratios of the right-hand-sides of the Gauss perturbing equations for the rates of change of the Keplerian orbital elements, evaluated onto the unperturbed Keplerian ellipse, to P b as functions of the true anomaly f. In the most general case, by assuming the co-rotation of the atmosphere, they are
n b 2 d a d f = C D Σ ρ a n b 1 e 2 V 1 + 2 e cos f + e 2 ν 1 e 2 3 / 2 cos I 2 1 + e cos f 2 ,
n b 2 d e d f = C D Σ ρ n b 1 e 2 3 / 2 V 8 1 + e cos f 4 4 e + cos f 1 + e cos f 2 ν 1 e 2 3 / 2 4 cos f + e 3 + cos 2 f 5 ,
n b 2 d I d f = C D Σ ρ n b ν 1 e 2 3 V sin I cos 2 u 4 1 + e cos f 4 ,
n b 2 d Ω d f = C D Σ ρ n b ν 1 e 2 3 V sin 2 u 8 1 + e cos f 4 ,
n b 2 d ω d f = C D Σ ρ n b 1 e 2 V 8 e 1 + e cos f 4 4 1 e 2 1 + e cos f 2 sin f + + 2 ν 1 e 2 2 cos I 2 sin f + e cos ω sin 2 f + ω ,
n b 2 d η d f = C D Σ ρ n b 1 e 2 2 V sin f 4 e 1 + e cos f 4 2 + 3 e 2 + 2 e 2 + e 2 cos f + e 2 cos 2 f ν 1 e 2 3 / 2 2 + e cos f cos I ,
where
V 2 = 1 ν 2 1 e 2 3 / 2 cos I 1 + e 2 + 2 e cos f + ν 2 1 e 2 3 3 + cos 2 I + 2 sin 2 I cos 2 u 4 1 + e cos f 2 1 + e 2 + 2 e cos f ,
and
u ω + f
is the argument of latitude. It is worthwhile noticing that, in general,
V 2 1 1 ,
being even possible that
V 2 1 1
for some values of f, thus preventing from expanding it in powers of ν . Then, we numerically calculate the areas under the resulting curves, i.e., we numerically integrate the averaged rates of change of the orbital elements for the given orbital configuration: see Table A5 and Table A9. In the second approach, we numerically integrate the equations of motion of the satellite in rectangular Cartesian coordinates, referred to a geocentric equatorial coordinate system, with and without Equation (A37) over 1 year; both the runs share the same initial conditions. Then, we subtract the resulting time series of the orbital elements in order to single out their shifts due to the disturbing acceleration. Finally, we perform a linear fit of them, and look at their slopes; we plot the fitted trends as functions of time t in Figure A2 and Figure A4. Both the methods reciprocally agree well, as shown by Figure A1, Figure A2, Figure A3 and Figure A4 for the neutral drag. A slight reduction turns out to occur if the decrease of the atmospheric co-rotation with height is taken into account as modeled by Membrado & Pacheco [59] by assuming the simpler case of a constant viscosity. We tested our approach by checking that it was able to reproduce the observed features of the semimajor axis decay of LARES recently determined in Pardini et al. [53]. However, it must be stressed that such findings should be deemed just as indicative of the limitations of the scenario considered if a passive spacecraft were to be adopted. Indeed, they were computed preliminarily by assuming the same physical properties of the existing LARES satellite which, in principle, could well be superseded by a new, specifically manufactured spacecraft able to reduce both the drag coefficient C D and the area-to-mass ratio Σ . Moreover, also the actual temporal variability of the atmospheric density over timescales larger than the satellite’s orbital period P b should be taken into account, especially if data will be collected during temporal intervals several years long. To this aim, it is important to note that an inspection of Equations (A49) and (A54) shows that no other sources of long-term modulation are present. Indeed, the circulating node Ω does not enter them, contrary to the perigee ω which, however, is held fixed by the adopted value of the inclination I.

References

  1. Brumberg, V.A. Essential Relativistic Celestial Mechanics; Adam Hilger: Bristol, UK, 1991. [Google Scholar]
  2. Meichsner, J.; Soffel, M.H. Effects on satellite orbits in the gravitational field of an axisymmetric central body with a mass monopole and arbitrary spin multipole moments. Celest. Mech. Dyn. Astr. 2015, 123, 1. [Google Scholar] [CrossRef]
  3. Soffel, M.H. Relativity in Astrometry, Celestial Mechanics and Geodesy; Springer: Heidelberg, Germany, 1989. [Google Scholar]
  4. Will, C.M. Incorporating post-Newtonian effects in -body dynamics. Phys. Rev. D 2014, 89, 044043. [Google Scholar] [CrossRef]
  5. Panhans, M.; Soffel, M.H. Gravito-magnetism of an extended celestial body. Class. Quantum Gravity 2014, 31, 245012. [Google Scholar] [CrossRef]
  6. Soffel, M.; Wirrer, R.; Schastok, J.; Ruder, H.; Schneider, M. Relativistic effects in the motion of artificial satellites: The oblateness of the central body I. Celest. Mech. Dyn. Astr. 1987, 42, 81. [Google Scholar] [CrossRef]
  7. Le Verrier, U.; de Lettre, M.; Le Verrier à, M. Faye sur la théorie de Mercure et sur le mouvement du périhélie de cette planète. Cr. Hebd. Acad. Sci. 1859, 49, 379. [Google Scholar]
  8. Einstein, A. Erklärung der Perihelbewegung des Merkur aus der allgemeinen Relativitätstheorie. Sitzungsberichte der Preußischen Akademie der Wissenschaften 1915, 47, 831–839. [Google Scholar]
  9. Shapiro, I.I. General Relativity and Gravitation; Ashby, N., Bartlett, D.F., Wyss, W., Eds.; Cambridge University Press: Cambridge, UK, 1990; pp. 313–330. [Google Scholar]
  10. Shapiro, I.I.; Pettengill, G.H.; Ash, M.E.; Ingalls, R.P.; Campbell, D.B.; Dyce, R.B. Mercury’s perihelion advance: determination by radar. Phys. Rev. Lett. 1972, 28, 1594. [Google Scholar] [CrossRef]
  11. Anderson, J.D.; Campbell, J.K.; Jurgens, R.F.; Lau, E.L.; Newhall, X.X.; Slade, M.A., III; Standish, E.M., Jr. Recent developments in Solar system tests of general relativity. In Proceedings of the Sixth Marcel Grossmann Meeting on General Relativity; Satō, F., Nakamura, T., Eds.; World Scientific: Singapore, 1992; pp. 353–355. [Google Scholar]
  12. Anderson, J.D.; Keesey, M.S.W.; Lau, E.L.; Standish, E.M., Jr.; Newhall, X.X. Tests of general relativity using astrometric and radio metric observations of the planets. Acta Astronaut. 1978, 5, 43–61. [Google Scholar] [CrossRef]
  13. Shapiro, I.I.; Ash, M.E.; Smith, W.B. Icarus: further confirmation of the relativistic perihelion precession. Phys. Rev. Lett. 1968, 20, 1517. [Google Scholar] [CrossRef]
  14. Shapiro, I.I.; Smith, W.B.; Ash, M.E.; Herrick, S. General relativity and the orbit of Icarus. Astron. J. 1971, 76, 588. [Google Scholar] [CrossRef]
  15. Kramer, M.; Stairs, I.H.; Manchester, R.N.; McLaughlin, M.A.; Lyne, A.G.; Ferdman, R.D.; Burgay, M.; Lorimer, D.R.; Possenti, A.; D’Amico, N.; et al. Tests of general relativity from timing the double pulsar. Science 2006, 314, 97. [Google Scholar] [CrossRef] [PubMed]
  16. Lucchesi, D.M.; Peron, R. Accurate measurement in the field of the earth of the general-relativistic precession of the LAGEOS II pericenter and new constraints on non-newtonian gravity. Phys. Rev. Lett. 2010, 105, 231103. [Google Scholar] [CrossRef] [PubMed]
  17. Lucchesi, D.M.; Peron, R. LAGEOS II pericenter general relativistic precession (1993–2005): Error budget and constraints in gravitational physics. Phys. Rev. D 2014, 89, 082002. [Google Scholar] [CrossRef]
  18. Lense, J.; Thirring, H. Über den Einfluss der Eigenrotation der Zentralkörper auf die Bewegung der Planeten und Monde nach der Einsteinschen Gravitationstheorie. Phys. Z. 1918, 19, 156. [Google Scholar]
  19. Pfister, H. On the history of the so-called Lense-Thirring effect. Gen. Relativ. Gravit. 2007, 39, 1735. [Google Scholar] [CrossRef]
  20. Pfister, H. Eleventh Marcel Grossmann Meeting on Recent Developments in Theoretical and Experimental General Relativity, Gravitation and Relativistic Field Theories; Kleinert, H., Jantzen, R.T., Ruffini, R., Eds.; World Scientific: Singapore, 2008; pp. 2456–2458. [Google Scholar]
  21. Pfister, H. Editorial note to: Hans Thirring, On the formal analogy between the basic electromagnetic equations and Einstein’s gravity equations in first approximation. Gen. Relativ. Gravit. 2012, 44, 3217. [Google Scholar] [CrossRef]
  22. Pfister, H. Relativity and Gravitation; Springer Proceedings in Physics; Bičák, J., Ledvinka, T., Eds.; Springer: Berlin/Heidelberg, Germany, 2014; Volume 157, pp. 191–197. [Google Scholar]
  23. Renzetti, G. History of the attempts to measure orbital frame-dragging with artificial satellites. Open Phys. 2013, 11, 531–544. [Google Scholar] [CrossRef] [Green Version]
  24. Lucchesi, D.M.; Anselmo, L.; Bassan, M.; Magnafico, C.; Pardini, C.; Peron, R.; Pucacco, G.; Visco, M. General Relativity Measurements in the Field of Earth with Laser-Ranged Satellites: State of the Art and Perspectives. Universe 2019, 5, 141. [Google Scholar] [CrossRef]
  25. Pugh, G. Proposal for a Satellite Test of the Coriolis Prediction of General Relativity; Research Memorandum 11, Weapons Systems Evaluation Group; The Pentagon: Washington, DC, USA, 1959. [Google Scholar]
  26. Schiff, L. Gravity probe B: final results of a space experiment to test general relativity. Phys. Rev. Lett. 1960, 4, 215. [Google Scholar] [CrossRef]
  27. Everitt, C.W.F.; DeBra, D.B.; Parkinson, B.W.; Turneaure, J.P.; Conklin, J.W.; Heifetz, M.I.; Keiser, G.M.; Silbergleit, A.S.; Holmes, T.; Kolodziejczak, J.; et al. Gravity probe B: final results of a space experiment to test general relativity. Phys. Rev. Lett. 2011, 106, 221101. [Google Scholar] [CrossRef]
  28. Everitt, C.W.F.; Muhlfelder, B.; DeBra, D.B.; Parkinson, B.W.; Turneaure, J.P.; Silbergleit, A.S.; Acworth, E.B.; Adams, M.; Adler, R.; Bencze, W.J.; et al. The Gravity Probe B test of general relativity. Class. Quantum Gravity 2015, 32, 224001. [Google Scholar] [CrossRef] [Green Version]
  29. Everitt, C.W.F.; Buchman, S.; Debra, D.B.; Keiser, G.M.; Lockhart, J.M.; Muhlfelder, B.; Parkinson, B.W.; Turneaure, J.P. Gyros, Clocks, Interferometers …: Testing Relativistic Gravity in Space. In Lecture Notes in Physics; Lämmerzahl, C., Everitt, C.W.F., Hehl, F.W., Eds.; Springer: Berlin, Germany, 2001; Volume 562, p. 52. [Google Scholar]
  30. Iorio, L. Post-Newtonian direct and mixed orbital effects due to the oblateness of the central body. Int. J. Mod. Phys. D 2015, 24, 1550067. [Google Scholar] [CrossRef]
  31. Iorio, L. The post-Newtonian gravitomagnetic spin-octupole moment of an oblate rotating body and its effects on an orbiting test particle; are they measurable in the Solar system? Mon. Not. R. Astron. Soc. 2019, 484, 4811. [Google Scholar] [CrossRef]
  32. Iorio, L. A possible new test of general relativity with Juno. Class. Quantum Gravity 2013, 30, 195011. [Google Scholar] [CrossRef]
  33. Abbas, U.; Bucciarelli, B.; Lattanzi, M.G. Differential astrometric framework for the Jupiter relativistic experiment with Gaia. Mon. Not. Roy. Astron. Soc. 2019, 485, 1147. [Google Scholar] [CrossRef]
  34. Crosta, M.T.; Mignard, F. Microarcsecond light bending by Jupiter. Class. Quantum Gravity 2006, 23, 4853. [Google Scholar] [CrossRef]
  35. Kopeikin, S.M.; Makarov, V.V. Gravitational bending of light by planetary multipoles and its measurement with microarcsecond astronomical interferometers. Phys. Rev. D 2007, 75, 062002. [Google Scholar] [CrossRef] [Green Version]
  36. Le Poncin-Lafitte, C.; Teyssandier, P. Influence of mass multipole moments on the deflection of a light ray by an isolated axisymmetric body. Phys. Rev. D 2008, 77, 044029. [Google Scholar] [CrossRef] [Green Version]
  37. Angélil, R.; Saha, P.; Bondarescu, R.; Jetzer, P.; Schärer, A.; Lundgren, A. Spacecraft clocks and relativity: Prospects for future satellite missions. Phys. Rev. D 2014, 89, 064067. [Google Scholar] [CrossRef]
  38. Schärer, A.; Angélil, R.; Bondarescu, R.; Jetzer, P.; Lundgren, A. Testing scalar-tensor theories and parametrized post-Newtonian parameters in Earth orbit. Phys. Rev. D 2014, 90, 123005. [Google Scholar] [CrossRef]
  39. Schärer, A.; Bondarescu, R.; Saha, P.; Angélil, R.; Helled, R.; Jetzer, P. Prospects for Measuring Planetary Spin and Frame-Dragging in Spacecraft Timing Signals. Front. Astron. Space Sci. 2017, 4, 11. [Google Scholar] [CrossRef]
  40. Delva, P.; Puchades, N.; Schönemann, E.; Dilssner, F.; Courde, C.; Bertone, S.; Gonzalez, F.; Hees, A.; Le Poncin-Lafitte, C.; Meynadier, F.; et al. Gravitational Redshift Test Using Eccentric Galileo Satellites. Phys. Rev. Lett. 2018, 121, 231101. [Google Scholar] [CrossRef] [PubMed]
  41. Herrmann, S.; Finke, F.; Lülf, M.; Kichakova, O.; Puetzfeld, D.; Knickmann, D.; List, M.; Rievers, B.; Giorgi, G.; Günther, C.; et al. Test of the Gravitational Redshift with Galileo Satellites in an Eccentric Orbit. Phys. Rev. Lett. 2018, 121, 231102. [Google Scholar] [CrossRef] [PubMed]
  42. Altschul, B.; Bailey, Q.G.; Blanchet, L.; Bongs, K.; Bouyer, P.; Cacciapuoti, L.; Capozziello, S.; Gaaloul, N.; Giulini, D.; Hartwi, J.; et al. Quantum tests of the Einstein Equivalence Principle with the STE–QUEST space mission. Adv. Space Res. 2015, 55, 501. [Google Scholar] [CrossRef]
  43. Clifton, T.; Ferreira, P.G.; Padilla, A.; Skordis, C. Modified gravity and cosmology. Phys. Rep. 2012, 513, 1. [Google Scholar] [CrossRef]
  44. Rubincam, D.P. On the secular decrease in the semimajor axis of LAGEOS’s orbit. Celest. Mech. Dyn. Astr. 1982, 26, 361. [Google Scholar] [CrossRef]
  45. Sośnica, K. Determination of Precise Satellite Orbits and Geodetic Parameters Using Satellite Laser Ranging; ETH Zürich, Swiss Geodetic Commission: Zürich, Switzerland, 2014; pp. 94–95. [Google Scholar]
  46. Sośnica, K.; Baumann, C.; Thaller, D.; Jäggi, A.; Dach, R. Combined LARES-LAGEOS Solutions. In Proceedings of the 18th International Workshop on Laser Ranging, Fujiyoshida, Japan, 11–15 November 2013. 13-Po-57. [Google Scholar]
  47. Ciufolini, I.; Lucchesi, D.; Vespe, F.; Mandiello, A. Measurement of dragging of inertial frames and gravitomagnetic field using laser-ranged satellites. Nuovo Cimento A 1996, 109, 575. [Google Scholar] [CrossRef]
  48. Chen, Q.; Shen, Y.; Francis, O.; Chen, W.; Zhang, X.; Hsu, H.J. Tongji-Grace02s and Tongji-Grace02k: High-Precision Static GRACE-Only Global Earth’s Gravity Field Models Derived by Refined Data Processing Strategies. Geophys. Res. 2018, 123, 6111. [Google Scholar] [CrossRef]
  49. Sheard, B.S.; Heinzel, G.; Danzmann, K.; Shaddock, D.A.; Klipstein, W.M.; Folkner, W.M.J. Intersatellite laser ranging instrument for the GRACE follow-on mission. Geodesy 2012, 86, 1083. [Google Scholar] [CrossRef]
  50. Sehnal, L. Effects of the terrestrial infrared radiation pressure on the motion of an artificial satellite. Celest. Mech. Dyn. Astr. 1981, 25, 169. [Google Scholar] [CrossRef]
  51. Durand-Manterola, H.J. Dipolar magnetic moment of the bodies of the solar system and the Hot Jupiters. Planet. Space Sci. 2009, 57, 1405. [Google Scholar] [CrossRef]
  52. Petit, G.; Luzum, B. IERS Conventions (2010). IERS Techn. Note 2010, 36, 1. [Google Scholar]
  53. Pardini, C.; Anselmo, L.; Lucchesi, D.M.; Peron, R. On the secular decay of the LARES semi-major axis. Acta Astronaut. 2017, 140, 469. [Google Scholar] [CrossRef]
  54. Lucchesi, D.M.; Anselmo, L.; Bassan, M.; Pardini, C.; Peron, R.; Pucacco, G.; Visco, M. Testing the gravitational interaction in the field of the Earth via satellite laser ranging and the Laser Ranged Satellites Experiment (LARASE). Class. Quantum Gravity 2015, 32, 155012. [Google Scholar] [CrossRef]
  55. Bertotti, B.; Farinella, P.; Vokrouhlický, D. Physics of the Solar System; Kluwer Academic Press: Dordrecht, The Netherlands, 2003. [Google Scholar]
  56. Capderou, M. Satellites: Orbits and Missions; Springer: Berlin, Germany, 2005. [Google Scholar]
  57. Kopeikin, S.; Efroimsky, M.; Kaplan, G. Relativistic Celestial Mechanics of the Solar System; Wiley-VCH: Weinheim, Germany, 2011. [Google Scholar]
  58. Xu, G. Orbits; Springer: Berlin, Germany, 2008. [Google Scholar]
  59. Membrado, M.; Pacheco, A.F. Decrease of the atmospheric co-rotation with height. Eur. J. Phys. 2010, 31, 1013. [Google Scholar] [CrossRef]
  60. Brito, T.P.; Celestino, C.C.; Moraes, R.V. Study of the decay time of a CubeSat type satellite considering perturbations due to the Earth’s oblateness and atmospheric drag. J. Phys. Conf. Ser. 2015, 641, 012026. [Google Scholar] [CrossRef]
1.
Prior works [37,38,39] on other proposed mission concepts showed that there is a window of perhaps an hour or so over around the perigee, over which these effects would be detectable.
2.
In that case, the aliasing Newtonian effect which must be disentangled from the pN perihelion precessions is due to the quadrupole mass moment J 2 of the Sun.
3.
Indeed, the harmonic perturbations contain also the satellite’s perigee which, for the critical inclination I crit = arcsin 2 / 5 adopted in Table A2 and Table A6, stays essentially constant becoming “frozen”. In this case, also the long-term, harmonic rates become secular trends.
4.
At least one of them must be affected by the pN effect one is looking for. In principle, the N orbital elements κ ( i ) may be different from one another belonging to the same satellite, or some of them may be identical belonging to different spacecraft (e.g., the nodes of two different vehicles).
5.
In general, it is not necessarily one of the parameters of the parameterized post-Newtonian (PPN) formalism, being possibly a combination of some of them.

Share and Cite

MDPI and ACS Style

Iorio, L. A HERO for General Relativity. Universe 2019, 5, 165. https://doi.org/10.3390/universe5070165

AMA Style

Iorio L. A HERO for General Relativity. Universe. 2019; 5(7):165. https://doi.org/10.3390/universe5070165

Chicago/Turabian Style

Iorio, Lorenzo. 2019. "A HERO for General Relativity" Universe 5, no. 7: 165. https://doi.org/10.3390/universe5070165

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop