Next Article in Journal
Species Abundance Modelling of Arctic-Boreal Zone Ducks Informed by Satellite Remote Sensing
Previous Article in Journal
The Synchrosqueezed Method and Its Theory-Analysis-Based Novel Short-Time Fractional Fourier Transform for Chirp Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamical Model of Rotation and Orbital Coupling for Deimos

1
Yunnan Observatories, Chinese Academy of Sciences, Kunming 650216, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Project Supported by MOE Key Laboratory of TianQin Project, Sun Yat-sen University, Zhuhai 519082, China
4
State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan 430079, China
5
Key Laboratory of Space Object and Debris Observation, PMO, CAS, Nanjing 210023, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(7), 1174; https://doi.org/10.3390/rs16071174
Submission received: 28 February 2024 / Revised: 23 March 2024 / Accepted: 24 March 2024 / Published: 27 March 2024
(This article belongs to the Section Satellite Missions for Earth and Planetary Exploration)

Abstract

:
This paper introduces a novel dynamical model, building upon the existing dynamical model for Deimos in the current numerical ephemerides, which only encompasses the simple libration effects of Deimos. The study comprehensively incorporates the rotational dynamics of Deimos influenced by the torque exerted by the major celestial bodies (Mars, the Sun) in the solar system within the inertial space. Consequently, a full dynamical model is formulated to account for the complete coupling between the rotation and orbit of Deimos. Simultaneously, employing precision orbit determination methods used for artificial satellites, we develop an adjustment model for fitting data to the complete model. The 12-order Adams–Bashforth–Moulton (ABM) integration algorithm is employed to synchronously integrate the 12 state variables of the full model to obtain the orbit of Deimos.The difference in the orbits obtained by integrating the full model over a period of 10 years and those obtained by the simplified model is at the order of 10 km. After precise orbit determination, this difference decreases to below 100 m, so numerical simulation results indicate that the full dynamical model and adjustment model are stable and reliable. Simultaneously, the integration of the Deimos third-order gravity field in the full model over a 10-year period induces only meter-level positional changes. This suggests that when constructing the complete model, the utilization of a second-order gravity field alone is sufficient. Compared to the simple model, the polar axis of Deimos in the inertial space exhibits a more complex oscillation in the full model. Additionally, the full model calculates that the minimum moment of inertia principal axis of Phobos has an amplitude of approximately 0.5 degrees in the longitude direction and does not exceed 2 degrees in the latitude direction. This work further advances the current dynamical model for Deimos and establishes the foundational model for the generation of a new set of precise numerical ephemerides for Deimos.

1. Introduction

Human understanding of the Martian moons began with the discovery by Asaph Hall, a scientist at the United States Naval Observatory, in 1877. He identified the two moons of Mars and named them Phobos and Deimos, respectively. Researchers have conducted comprehensive studies on Phobos [1,2]. Compared to Phobos, there is relatively limited publicly available research on Deimos. Therefore, this paper focuses on Deimos as the primary research target. The physical parameters of Deimos are shown in Table 1 [3]. The study of the Martian moons not only contributes to a deeper understanding of the formation and evolution of terrestrial planet systems but also provides crucial clues in researching the formation and evolution of the solar system. Up to now, humans have launched multiple probes to Mars and conducted extensive exploration activities on its two moons. Approximately 90 years after the initial discovery of the two Martian moons, the first artificial satellite dedicated to Mars, Mariner 4, conducted a flyby of Mars in 1965 [4]. Subsequently, human exploration of Mars has not ceased, and spacecraft have been launched to Mars on multiple occasions. Mariners 6 and 7 also carried out flyby missions of Mars. In 1971, the first Mars orbiter, Mariner 9, was launched. It obtained 214 images of Phobos and Deimos during close flybys, marking the earliest spacecraft observations of the Martian moons. In the late 1970s, the Soviet Union successively launched two orbiters, Viking 1 and Viking 2, to Mars. These spacecraft were equipped with landers. Among them, Mars Express, launched by the European Space Agency (ESA), acquired the most data on the Martian moons. Its impact on understanding the orbits of the Martian moons and the Martian ephemerides has been profound. Notably, data from the Pathfinder Mars lander played a crucial role in the study of Mars’ rotation and gravity field. The gravity field model MRO120F from the Pathfinder mission is utilized in this research [5]. In recent years, both China and Japan have initiated their own Mars exploration missions [6,7,8]. It is anticipated that an abundance of high-precision data will be instrumental in advancing the research on the dynamical models of Martian moons.
In various exploration activities of celestial bodies beyond Earth, investigating their orbital trajectories is a crucial scientific inquiry, holding significant academic importance for a profound understanding of the internal structure and physical parameters of these celestial entities. The data collected from ground-based observations to spacecraft missions provide the foundation for our research on the dynamical models of Martian moons. Multiple international institutions have conducted in-depth research in constructing dynamical models for Martian moons. With the improvement in the volume and precision of human exploration data for Martian moons, the dynamical models for these moons have undergone a developmental process from analytical models and semi-analytical models to numerical models.
In the early stages, dynamical models for Martian moons were typically established using analytical methods. The earliest analytical model was developed by Struve in 1911. He fitted observational data from the United States Naval Observatory collected between 1877 and 1909, calculating numerical values for the Martian moon’s mean longitude, orbital elements, mean motion, and precession rate. It was not until Shor that an analytical theory was introduced that incorporated the influence of tidal acceleration into the model [9]. This modification involved adding an acceleration term to the mean orbital longitude, and fitting was conducted using observational data spanning from 1877 to 1973.
With the continuous improvement in observational precision, analytical methods gradually revealed their limitations, leading to the brief emergence of semi-analytical approaches in dynamical modeling. The theory proposed by Chapront-Touzé [10,11,12] represents a high-precision semi-analytical model. It incorporates perturbations from the Sun and planets, the gravitational field of Mars, and mutual perturbations between Phobos and Deimos. This model is a more comprehensive semi-analytical approach to constructing dynamical models for Martian moons. Through fitting observational data from both ground-based and spacecraft observations spanning from 1877 to 1988, an ephemeris for the Martian moons was established.
With the development of deep space exploration, especially the increase in observational data and the continuous improvement in computing power, numerical methods, characterized by high computational accuracy and simplicity of formulas, have gradually replaced analytical methods as the primary approach in planetary orbit calculations. In 1977, Tolson [13] used numerical methods and Viking spacecraft data to generate the first ephemeris for Phobos. In 2007, Lainey et al. [14] developed the first complete numerical integration ephemeris for the Martian moon. The dynamic model included a 10th-order non-spherical Martian gravity field [15], perturbations from the Sun, Jupiter, Saturn, Earth, and the Moon, mutual perturbations between the moons, tidal effects caused by the two moons on Mars, and the quadrupole moment of Phobos. Subsequently, Jacobson [16,17] further incorporated the librational effects of the moons in the dynamical model, considering perturbations caused by the oscillation of the axis pointing toward Mars in the moon’s body-fixed coordinate system within the orbital plane. Additionally, a new Mars rotation model was adopted, and the model was fitted with the Mars Express dataset, generating the Martian moon ephemeris MARS097.
In Jacobson’s dynamical model, while the libration of the minimum moment of principal axis of Deimos in the orbital plane is considered, it does not incorporate the full rotation of the moon as is the case in the lunar dynamical model. Therefore, the aim of this paper is to establish a rotational model for Deimos in the inertial frame based on the Euler–Liouville rotational equations, following the approach used in the lunar dynamical model. The purpose is to further develop the simple dynamical model for Martian moons and establish a full dynamical model for Deimos in the inertial frame. Simultaneously, based on the rotation model for Deimos, the study incorporates the influence of Deimos’ own gravitational field when it orbits Mars, completing the full dynamical model of the rotational and orbital coupling for Deimos.
The first section of this paper will introduce the establishment process and distinctions between the simple model and the dynamical model considering the coupling of rotation and orbit. In the second section, the least squares method will be employed to introduce the data-fitting generation process and the adjustment model for precise orbit determination using the established dynamical model. Finally, a comparison will be made between the changes in coordinate axes during the orbital integration in inertial space for the two models.

2. The Establishment of the Simple Model and Full Model

2.1. The Rotation of the Mars

When describing the motion of Martian moons, the coordinate system chosen is the International Celestial Reference Frame (ICRF), while the gravitational field is represented in the Martian body-fixed coordinate system. Therefore, calculating the acceleration induced by the Martian gravity field on the moons requires an appropriate rotation model. Konopliv et al. conducted an in-depth study of the Mars rotation based on data from the Pathfinder lander, establishing the most accurate current rotation model for Mars [5]. According to this model, the method for converting positions from the Mars-fixed frame to the ICRF is as follows:
R I N = R Z ( N ) R X ( J ) R Z ( ψ ) R X ( I ) R Z ( ϕ ) R Y ( x p ) R X ( y p ) R B F
where R I N and R B F represent the coordinates of Deimos in the inertial frame and the body-fixed frame, respectively. N is the angle in the plane of the Earth-mean-equator of J2000 (EME2000) from the vernal equinox (intersection of the mean ecliptic and EME2000 planes) to the intersection of the Mars-mean-orbit of J2000 and EME2000 planes; J is the inclination of the Mars-mean-orbit plane relative to the EME2000 plane; ψ is the angle in the Mars-mean-orbit plane from the intersection of the EME2000 and Mars-mean-orbit planes to the intersection of the Mars-mean-orbit and Mars-true-equator of date planes; I is the inclination of the Mars true equator relative to the Mars-mean-orbit plane; ϕ is the angle in the plane of the Mars true equator from the intersection of the Mars-mean-orbit plane and Mars true equator to the Mars prime meridian; and x p , y p represents the polar motion.

2.2. Mars Gravity Field

Due to the irregular shape and uneven mass distribution of Mars, when calculating the gravity on Deimos, it cannot be treated as an ideal homogeneous sphere. Non-spherical perturbations should be considered. The gravitational potential function is expanded into spherical harmonics as follows:
V = G M r n = 2 m = 0 n R mars r n P ¯ n m ( sin ϕ ) C ¯ n m cos m λ + S ¯ n m sin m λ
In this expression, G M is the gravitational constant of Mars, R mars is the radius of Mars, P ¯ n m represents the normalized associated Legendre function, and C ¯ n m and S ¯ n m are the normalized spherical harmonic coefficients of the Mars gravity field model. λ , ϕ , and r denote the longitude, latitude, and distance to the center of Mars in the fixed coordinate system of Deimos on Mars. When using the gravity field formula, it is necessary to follow the coordinate transformation method described in the previous section. First, the coordinates of Deimos in the inertial frame should be converted to the Mars-fixed frame. Subsequently, the integration calculations are performed by transforming from the Mars-fixed frame back to the inertial frame.
When computing gravitational field acceleration, a recursive method is commonly employed to avoid complex differentiation processes. For the Deimos orbit calculation, taking the example of a 10th-order Martian gravity field, the Cunningham [18] algorithm is used to calculate the gravitational field and its partial derivatives.
When calculating the Martian gravity potential, Cunningham provided a single recursion formula that can be used to calculate the gravity potential and the corresponding acceleration in Cartesian coordinates. The formula is defined as follows:
V n m = R r n + 1 · P n m ( sin ϕ ) · cos m λ W n m = R r n + 1 · P n m ( sin ϕ ) · sin m λ
where the gravitational field can be represented as (non-normalized coefficients are used here):
U = G M R n = 0 m = 0 n C n m V n m + S n m W n m
Define the relationship:
V m m = ( 2 m 1 ) x R r 2 V m 1 , m 1 y R r 2 W m 1 , m 1 W m m = ( 2 m 1 ) x R r 2 W m 1 , m 1 + y R r 2 V m 1 , m 1 V n m = 2 n 1 n m · z R r 2 · V n 1 , m n + m 1 n m · R 2 r 2 · V n 2 , m W n m = 2 n 1 n m · z R r 2 · W n 1 , m n + m 1 n m · R 2 r 2 · V n 2 , m
Then, the acceleration of the gravitational field can be represented as:
x ¨ n m = G M R 2 C n 0 V n + 1 , 1 , m = 0 G M R 2 · 1 2 · C n m V n + 1 , m + 1 S n m W n + 1 , m + 1 + ( n m + 2 ) ! ( n m ) ! · C n m V n + 1 , m 1 + S n m W n + 1 , m 1 , m > 0
y ¨ n m = G M R 2 C n 0 W n + 1 , 1 , m = 0 G M R 2 · 1 2 · C n m W n + 1 , m + 1 + S n m V n + 1 , m + 1 + ( n m + 2 ) ! ( n m ) ! · C n m V n + 1 , m 1 + S n m W n + 1 , m 1 , m > 0
z ¨ n m = G M R 2 · ( n m + 1 ) · C n m V n + 1 , m S n m V n + 1 , m

2.3. General Relativity Perturbation Model

During the motion of Deimos, it is mainly influenced by the General Relativity effect generated by the mass of Mars. Since the eccentricity of Deimos’ orbit is very small, it can be considered a nearly circular orbit. In this case, the acceleration induced by the General Relativity effect of Mars on Deimos is given by:
a r e l = G M | r | 3 r 3 v 2 c 2
where G M is the gravitational constant of Mars, v is the velocity of Deimos in the Mars inertial frame, r is the position vector from the Mars to Deimos, and c is the speed of light.

2.4. Solid Tide Perturbation

Due to the influence of the gravitational forces from other celestial bodies in the solar system, Mars experiences changes in its gravitational field, causing solid tidal effects. For a celestial body with a mass M j and a radius a j , the solid tidal effect (according to elastic model pointed in [19] and also [20]) on Deimos’ gravitational on Deimos’ gravitational potential function is given by [21]:
V tide = k 2 2 G M j r j 3 a j 5 | r | 3 3 r j · r 2 1
where r j = r r m   r j is the position vector from the perturbing body to Deimos, r is the position vector from Mars to Deimos, r m is the position vector from Mars to the perturbing body, and the Love number of Mars has a value of 0.169 [20]. The solid tide model used in this study does not incorporate the tidal energy dissipation effects of Mars. We referred to previous studies on tidal energy dissipation [22,23,24,25,26], and integrating the established dynamical model over ten years in numerical experiments reveals that tidal energy dissipation effects only cause differences of less than a few meters. Therefore, within the time scale of this study (10 years), we have not considered the influence of tidal dissipation.

2.5. Simple Model

The current Deimos dynamical model is referred to as a simplified model. In the simple Deimos dynamical model, the librations of Deimos are considered. Jacobson [16] introduced the specific expression for this effect:
F = 3 2 μ R 2 r 4 J 2 + 6 C 22 cos 2 θ r ^ + 4 C 22 sin 2 θ t ^
θ = f M + 2 e 1 C 3 ( B A ) sin M
where μ represents the gravitational constant of Mars; R is the radius of Deimos; r represents the distance from Mars to Deimos; J 2 is the second-degree sectoral harmonic coefficient of Deimos; C 22 is the second-degree tesseral harmonic coefficient of Deimos; f and M are the true anomaly and mean anomaly of Deimos’ orbit, respectively; e is the eccentricity of the orbit; A, B, and C are the moments of inertia with respect to the principal axes of Deimos; r ^ is the unit vector from Mars to Deimos; and t ^ is the unit vector in the Deimos orbital plane that is perpendicular to the velocity vector of Deimos. This model is the dynamical model used in the current ephemeris.

2.6. Full Model

The main difference between the full model and the simple model lies in the consideration of the complete rotation of Deimos in inertial space in the full model. There have been some research efforts on the rotation of Phobos [27,28]. In this section, we will establish the rotation model of Deimos based on the theory of rigid body rotation. Based on the established rotation model, it is possible to calculate the gravitational force of Deimos on Mars and the acceleration of Deimos caused by the reaction force it experiences due to the gravitational interaction with Mars. Thus, a dynamical model for the coupled rotation and orbit of Deimos has been established. The gravitational field coefficients for Deimos are listed in Table 2 [3]. This study further advances the current dynamical model that only considers the simple librations of Deimos. It comprehensively calculates the rotation of Deimos in inertial space due to torques from celestial bodies such as Mars and the Sun.
The differential equations governing the rotational motion of a rigid body are classically represented by the Euler–Liouville equations, and their form in an inertial frame is as follows:
d ( I ω ) d t + ω × I ω = N
where N is the torque acting on Deimos, ω represents the rotation of Deimos, and I is inertia tensor matrix.
When describing the rotation of Deimos, we take the center of mass of Deimos as the origin, the three principal axes of Deimos as the coordinate axes, and establish the Deimos Principal Axis (PA) reference frame rotating with Deimos. Simultaneously, we can obtain the celestial reference frame of the Deimos center of mass with coordinates aligned with ICRS. Assuming Deimos as an ideal rigid body, the transformation between the Deimos body-fixed coordinate system and the inertial frame is determined by three Euler angles ( ϕ , θ , ψ ) describing the rotation of the rigid body. This transformation process is expressed as:
R I N = R Z ( ϕ ) R X ( θ ) R Z ( ψ ) R B F
ω represents the rotation of the rigid-body moon under the action of torque in the inertial coordinate system, and the relationship between ω and the three Euler angles and the rates of change of the three Euler angles can be obtained:
ω x = ϕ ˙ sin θ sin ψ + θ ˙ cos ψ ω y = ϕ ˙ sin θ cos ψ θ ˙ sin ψ ω z = ϕ ˙ cos θ + ψ ˙
Solving for the first derivative of ω , we can obtain the accelerations of the three Euler angles:
ϕ ¨ = ω ˙ x   sin   ψ   +   ω ˙ x   cos   ψ   +   θ ˙   ( ψ ˙     ϕ ˙   cos   θ ) sin   θ θ ¨ = ω ˙ x cos ψ ω ˙ y sin ψ ϕ ˙ ψ ˙ sin θ ψ ¨ = ω ˙ z ϕ ¨ cos θ + ϕ ˙ ψ sin θ
However, when employing Equation (16), one may encounter singular situations. In such cases, using Wisdom angles instead of Euler angles might be more beneficial [29]. According to the Euler–Liouville equation, the value of ω ˙ depends on the torque acting on the Deimos and the moment of inertia. Taking the first derivative of Equation (13) yields the expression:
ω ˙ = I 1 ( N ω × I ω )
I is the inertia tensor matrix, which can be expressed as:
I = A 0 0 0 B 0 0 0 C
The torque acting on Deimos during its rotation can be divided into three components: (1) the torque produced by Mars as a point mass in the Deimos gravitational field; (2) the torque exerted by the Sun, treating it as a point mass due to the presence of the Deimos gravitational field; and (3) the torque generated in the Deimos due to the irregular shape of Mars, under the influence of the Deimos gravitational field [30]. The formulas can be expressed as:
N = N + N + N F
where N and N represent the torques generated by the point masses Mars and the Sun, and N F represents the torque produced by the shape effect of Mars [31]. N F is generated by the second-order term of the Martian gravitational field and is mainly related to the value of Mars. In numerical experiments, the torque generated by Phobos is a magnitude smaller than that of Mars and the Sun. Therefore, the model does not consider the impact of the torque generated by Phobos. The force exerted by Mars in the moon’s gravitational field is expressed as:
F = M U d
The torque produced by this force in the Deimos is expressed as:
N = M r d × U d
M represents the mass of Mars, r d represents the coordinates of Deimos in the Mars-centered inertial frame, and U d represents the gravity field of Deimos. And N F represents the torque induced by the oblateness of Mars and the shape effect of Deimos. The expression is given by [30]:
N F m D = 15 μ M R M 2 J 2 M 2 r M D 5 1 7 r ^ M D · p ^ 2 r ^ M D × I m D r ^ M D + 2 r ^ M D · p ^ r ^ M D × I m D p ^ + p ^ × I m D r ^ M D 2 5 p ^ × I m D p ^
m D is the mass of Deimos, J 2 M represents the value of the Mars flattening, r M D is the distance between Deimos and Mars, r ^ is the unit vector pointing from Deimos to Mars, p ^ is the direction of Mars’s polar axis, and all vectors are described in the Deimos body-fixed coordinate system.
According to the established model, the initial Euler angles of Deimos at the J2000 epoch are provided by Brent et al. [32]. The integration of the established rotational differential equations using the ABM integration algorithm yields the results for the three Euler angles as shown in Figure 1 and Figure 2.

2.7. Result of the Model

At this point, we have considered factors influencing the motion of Deimos, utilizing a second-order gravity field for Deimos. We have established a dynamical model for the coupled rotation and orbit of Deimos. The initial value of Deimos’ position at the J2000 epoch is provided by MARS097. The integration orbit for both the simple dynamical model and the full model with rotation–orbit coupling established in this study are compared in Figure 3.
The experimental results indicate that, considering the J2000 epoch as the integration starting point with a time span of 10 years, the difference in the computed results between the two dynamical models is around 10 km. The difference increases linearly over time, primarily due to the inclusion of rotation and the gravitational field of Deimos itself in the dynamical model. This discrepancy suggests that considering the complete rotation effect of Deimos in the dynamical model is essential for improving its accuracy.
When considering the third-order gravity field of Deimos, the integrated orbit of Deimos shows differences compared to the integrated orbit of the dynamical model considering only the second-order gravity field of Deimos as shown in Figure 4.
The model results indicate that, after integrating for 10 years, the difference between the second- and third-order gravity fields of Deimos is only on the order of meters, which is an extremely small value. Therefore, when constructing a comprehensive model for Deimos, considering gravity field orders up to the second order is sufficient.

3. Adjustment Model

In this section, we conducted a precise orbit determination for the established full dynamical model of Deimos, based on the methodology employed in the precise orbit determination of artificial satellites [33]. When selecting data for the fitting process, we initially utilized observations from the ephemeris MARS097 as the basis for refining the simple dynamical model. The resulting fitted orbit serves as the input data for the full model. This approach ensures consistency in the parameters between the two dynamical models, with the integrated results highlighting differences arising solely from distinct rotational modes between the simple and full models.

3.1. Variation Equations

For a specified epoch, changes in the initial values of the state vector result in variations in the position and velocity at subsequent epochs, represented by the state transition matrix. The variation equations can be established from the state transition matrix:
d ( Φ , S ) d x = 0 3 × 3 1 3 × 3 a r a v 6 × 6 · ( Φ , S ) + 0 3 × 6 0 3 × n p 0 3 × 6 a p 6 × 6 + n p
For the full model, involving both the rotation and orbital parameters (12 state variables in total), the structure of the state transition matrix is outlined in Equation (24):
x y z v x v y v z ϕ θ ψ ϕ ˙ θ ˙ ψ ˙ v x 0 0 0 1 0 0 0 0 0 0 0 0 v y 0 0 0 0 1 0 0 0 0 0 0 0 v z 0 0 0 0 0 1 0 0 0 0 0 0 a x a x x a x y a x z a x v x a x v y a x v z a x ϕ a x θ a x ψ 0 0 0 a y a y x a y y a y z a y v x a y v y a y v z a y ϕ ˙ a y θ a y ψ 0 0 0 a z a z x a z y a z z a z v x a z v y a z v z a z ϕ a z θ a z ψ 0 0 0 ϕ ˙ 0 0 0 0 0 0 0 0 0 1 0 0 θ ˙ 0 0 0 0 0 0 0 0 0 0 1 0 ψ ˙ 0 0 0 0 0 0 0 0 0 0 0 1 ϕ ¨ ϕ ¨ x ϕ ¨ y ϕ ¨ z 0 0 0 ϕ ¨ ϕ ˙ ϕ ¨ θ ϕ ¨ ψ ϕ ¨ ϕ ˙ ϕ ¨ θ ˙ ϕ ¨ ψ ˙ θ ¨ θ ¨ x θ ¨ y θ ¨ z 0 0 0 θ ¨ ϕ θ ¨ θ θ ¨ ψ θ ¨ ϕ ˙ θ ¨ θ ˙ θ ¨ ψ ˙ ψ ¨ ψ ¨ x ψ ¨ y ψ ¨ z 0 0 0 ψ ¨ ϕ ψ ¨ θ ψ ¨ ψ ψ ¨ ϕ ˙ ψ ¨ θ ˙ ψ ¨ ψ ˙
Let φ ¨ F 1 , θ ¨ F 2 , ψ ¨ F 3 , and calculate the first derivative of F 1 , F 2 , F 3 :
F 1 ϕ = csc θ sin ψ ω ˙ x ϕ + cos ψ ω ˙ y ϕ F 2 ϕ = cos ψ ω ˙ x ϕ + sin ψ ω ˙ y ϕ F 3 ϕ = ω ˙ z ϕ cos θ F 1 ϕ
F 1 θ = csc θ sin ψ ω ˙ x θ + cos ψ ω ˙ y θ cot θ csc θ ω ˙ x sin ψ + ω ˙ y cos ψ θ ˙ ψ ˙ cot θ csc θ + θ ˙ ϕ ˙ cot θ csc 2 θ F 2 θ = cos ψ ω ˙ x θ sin ψ ω ˙ y θ θ ˙ ϕ ˙ cos θ F 3 θ = ω ˙ z θ cos ψ F 1 θ + sin θ F 1 + θ ˙ ϕ ˙ cos θ
F 1 ψ = csc θ ω ˙ x cos ψ + sin ψ ω ˙ x ψ ω ˙ y sin ψ + cos ψ ω ˙ y ψ F 2 ψ = cos ψ ω ˙ x ψ sin ψ ω ˙ x cos ψ ω ˙ y sin ψ ω ˙ y ψ F 3 ψ = ω ˙ z ψ cos θ F 1 ψ
According to this method, the expression for the first derivative with respect to Euler angular rates F A ˙ can also be derived. According to the Euler–Liouville equation, the expression for the partial derivative of ω ˙ with respect to Euler angles ( ω ˙ A ) is given by:
ω ˙ A = I 1 N A ω A × I ω ω × I ω A
N A = M b r b A × U b + r b × U b A
Coordinates in the body-fixed coordinate system of Deimos are denoted by ( ξ 1 , ξ 2 , ξ 3 ) . According to the transformation matrix from inertial frame to body-fixed frame, the expression for r A can be solved as:
x ϕ = ( cos ϕ sin ψ sin ψ cos θ cos ϕ ) ξ 1 + ( cos ϕ cos ψ sin ψ cos θ sin ϕ ) ξ 2 x θ = sin ψ sin θ sin ϕ ξ 1 sin ψ sin θ cos ϕ ξ 2 + sin ψ sin θ ξ 3 x ψ = ( cos ϕ sin ψ sin ψ cos θ cos ψ ) ξ 1 + ( sin ϕ sin ψ + cos ψ cos θ sin ϕ ) ξ 2 + cos ψ cos θ ξ 3
y ϕ = ( cos ψ sin ϕ cos ψ cos θ cos ϕ ) ξ 1 + ( cos ϕ sin ψ sin ϕ cos θ cos ψ ) ξ 3 y θ = cos ψ sin θ sin ϕ ξ 1 cos ψ sin θ cos ϕ ξ 2 + cos ψ cos θ ξ 3 y ψ = ( cos ϕ cos ψ + sin ψ cos θ sin ϕ ) ξ 1 + ( sin ϕ cos ψ + sin ψ cos θ cos ϕ ) ξ 2 sin ψ sin θ ξ 3
z ϕ = sin ψ cos ϕ ξ 1 + sin θ sin ϕ ξ 3 z θ = cos θ sin ϕ ξ 1 cos θ cos ϕ ξ 2 sin θ ξ 3 z ψ = 0
Finally, employing the precision orbit determination method, the computed state transition matrix is synchronously integrated with the orbit. Subsequently, utilizing a least-squares fit to observational data, the calculated results are compared with the MARS097 ephemeris as depicted in Figure 5 and Figure 6.
With the J2000 epoch as the integration starting point, the maximum discrepancy between the final results of the data fitting using the full dynamical model for Deimos and the MARS097 ephemeris is approximately 50 m. This demonstrates the stability and reliability of the established dynamical model and adjustment approach.

3.2. Variation of the Deimos Principal Axes in the Inertial Frame

The fundamental distinction between the full model and the simple dynamical model lies in the incorporation of the Deimos’s full rotation in the inertial frame. In the simple model, the moon’s polar axis direction is solely perpendicular to the orbital plane. Nevertheless, when accounting for the complete rotation, the moon’s polar axis undergoes oscillation in the inertial frame as the moon rotates. Figure 7 provides a comparison of the directional deviation of the polar axis under the two motion modes.
Starting from the J2000 epoch, Figure 7 illustrates the evolution of the angle between the polar axis and the vector (0, 0, 1) in the inertial frame over a 50-year period for both motion models. In the simplified model, the polar axis exhibits relative stability in space, primarily influenced by the precession effect. However, in the full model, accounting for the rotational forces from Mars and the Sun, the polar axis undergoes intricate motion in the inertial frame. Figure 8 and Figure 9 represent the angles between the minimum moment of inertia axis of the moon in the moon-fixed coordinate system and the Mars position vector in the longitude and latitude directions, respectively. The results reveal that the moon’s amplitude in the longitude direction is approximately 0.5 degrees, and in the latitude direction, it does not exceed 2 degrees.

4. Conclusions

This study presents a novel model describing the motion of Deimos in an inertial frame. Building upon existing models, it thoroughly calculates the rotational effects induced by the irregular shape of Deimos itself and the torques exerted by large celestial bodies (Mars, Sun) in the solar system in inertial space. Furthermore, it considers the coupling effect resulting from rotation and orbit, establishing a full dynamical model for the rotation of Deimos coupled with its orbit. The model is integrated using the Adams–Bashforth–Moulton (ABM) integration algorithm to concurrently solve the 12 elements governing rotation and orbital dynamics in the differential equations.
To validate the model, precise orbit determination techniques for artificial satellites are employed to fit data to the full model. The results demonstrate the stability and reliability of the established model. Additionally, the study derives an analytical expression for the solution of the coupling effect between Deimos’ rotation and orbit. The variational equation, starting from the J2000 epoch, is integrated over ten years. The final outcomes reveal differences of less than 50 m when compared to the current Deimos ephemeris MARS097, showcasing its practical applicability in engineering.
As various exploration missions targeting the Martian moons continue in the future, it is anticipated that a considerable amount of high-precision data on these moons will be generated. Leveraging these data in conjunction with the full dynamical model established in this study, there is potential to create a new generation of more accurate ephemerides for the Martian moons. Additionally, by utilizing the full model established in this paper along with future observational data, there is the prospect of quantitatively studying orbital resonances among celestial bodies in the solar system [34,35] and the Lense–Thirring effect [36,37,38] in orbital mechanics. This, in turn, holds the potential to further enhance the accuracy of the Deimos’ dynamical model.
The model established in this paper lays the foundation for further constructing an accurate ephemeris for Deimos. It demonstrates versatility in solving the dynamical equations for celestial bodies similar to Deimos and can serve as a reference for developing numerical ephemerides for major celestial bodies in the solar system.

Author Contributions

This paper is a collaborative work by all of the authors. K.H.: Conceptualization, methodology, software, and writing—original draft preparation; L.Z.: software; Y.Y.: Conceptualization, methodology; M.Y.: methodology, software; Y.L.: Conceptualization. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the National Natural Science Foundation of China (12033009, 12103087), the National Key Research and Development Program of China (2021YFA0715101), the International Partnership Program of Chinese Academy of Sciences (020GJHZ2022034FN), the Yunnan Province Foundation (202201AU070225, 202301AT070328, 202401AT070141), the Yunnan Revitalization Talent Support Program, grant from Key Laboratory of Lunar and Deep Space Exploration, Chinese Academy of Sciences (LDSE202004), open research fund of state key laboratory of information engineering in surveying, mapping and remote sensing, Wuhan University (21P02), and supported by Key Laboratory of TianQin Project (Sun Yat-sen University), Ministry of Education. M. Ye is supported by Natural Science Foundation of Hubei Province, China, No. 2022CFB123, Space Optoelectronic Measurement and Perception Laboratory, Beijing Institute of Control Engineering (No. LabSOMP-2023-09).

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

Deimos’ ephemerides files can be downloaded from https://naif.jpl.nasa.gov/pub/naif/INSIGHT/kernels/spk/mar097s.bsp/ (accessed on 1 May 2023) (Mars097), and https://ftp.imcce.fr/pub/ephem/satel/NOE/JUPITER/ (accessed on 1 May 2023) (NOE-4-2020). And the EPM ephemeris is obtained from the Online Ephemeris Service http://iaaras.ru/en/dept/ephemeris/online/ (accessed on 1 May 2023).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Kamada, A.; Kuroda, T.; Terada, N.; Kobayashi, M.; Nakagawa, H.; Miyamoto, H. Modeling 4.3 billion years of water history on Phobos. Icarus 2024, 410, 115916. [Google Scholar] [CrossRef]
  2. Kuramoto, K.; Kawakatsu, Y.; Fujimoto, M.; Araya, A.; Barucci, M.A.; Genda, H.; Hirata, N.; Ikeda, H.; Imamura, T.; Helbert, J.; et al. Martian moons exploration MMX: Sample return mission to Phobos elucidating formation processes of habitable planets. Earth Planets Space 2021, 74, 12. [Google Scholar] [CrossRef]
  3. Rubincam, D.P.; Chao, B.F.; Thomas, P.C. The Gravitational Field of Deimos. Icarus 1995, 114, 63–67. [Google Scholar] [CrossRef]
  4. Anderson, H.R. Mariner iv measurements near mars—Initial results. spacecraft description and encounter sequence. Science 1965, 149, 1226–1228. [Google Scholar] [CrossRef]
  5. Konopliv, A.S.; Park, R.S.; Folkner, W.M. An improved JPL Mars gravity field and orientation from Mars orbiter and lander tracking data. Icarus 2016, 274, 253–260. [Google Scholar] [CrossRef]
  6. Zezhou, S.; Wei, R.; Yang, J.; Chuang, W.; Jie, D.; Baichao, C. Key Mission System Design of Tianwen-1 Mars Probe. Aerosp. Control. Appl. 2021, 47, 9–16. [Google Scholar]
  7. Kuramoto, K.; Kawakatsu, Y.; Fujimoto, M.; Bibring, J.P.; Lawrence, D.J.; Genda, H.; Hirata, N.; Imamura, T.; Kameda, S.; Kobayashi, M.; et al. Martian Moons eXploration MMX: Current Status Report 2020. In Proceedings of the JpGU—AGU Joint Meeting 2020, Chiba, Japan, 24–28 May 2020. [Google Scholar]
  8. Croswell, K. An ambitious new mission seeks the origin of the moons of Mars. Proc. Natl. Acad. Sci. USA 2023, 120, e2302307120. [Google Scholar] [CrossRef] [PubMed]
  9. Shor, V.A. The motion of the Martian satellites. Celest. Mech. 1975, 12, 61–75. [Google Scholar] [CrossRef]
  10. Chapront-Touzé, M. ESAPHO: A semi-analytical theory for the orbital motion of Phobos. Astron. Astrophys. 1988, 200, 255–268. [Google Scholar]
  11. Chapront-Touzé, M. Orbits of the Martian satellites from ESAPHO and ESADE theories. Astron. Astrophys. 1990, 240, 159–172. [Google Scholar]
  12. Chapront-Touzé, M. Phobos’ physical libration and complements to the ESAPHO theory for the orbital motion of Phobos. Astron. Astrophys. 1990, 235, 447–458. [Google Scholar]
  13. Tolson, R.H.; Blackshear, W.T.; Mason, M.L.; Kelly, G.M. The mass of Phobos. Geophys. Res. Lett. 1977, 4, 551–554. [Google Scholar] [CrossRef]
  14. Lainey, V.; Dehant, V.; Pätzold, M. First numerical ephemerides of the Martian moons. Astron. Astrophys. 2007, 465, 1075–1084. [Google Scholar] [CrossRef]
  15. Lainey, V.; Duriez, L.; Vienne, A. New accurate ephemerides for the Galilean satellites of Jupiter. I. Numerical integration of elaborated equations of motion. Astron. Astrophys. 2004, 420, 1171–1183. [Google Scholar] [CrossRef]
  16. Jacobson, R.A. The Orbits and Masses of the Martian Satellites and the Libration of Phobos. Astron. J. 2010, 139, 668–679. [Google Scholar] [CrossRef]
  17. Lainey, V.; Pasewaldt, A.; Robert, V.; Rosenblatt, P.; Jaumann, R.; Oberst, J.; Roatsch, T.; Willner, K.; Ziese, R.; Thuillot, W. Mars moon ephemerides after 14 years of Mars Express data. Astron. Astrophys. 2020, 650, A64. [Google Scholar] [CrossRef]
  18. Cunningham, L.E. On the computation of the spherical harmonic terms needed during the numerical integration of the orbital motion of an artificial satellite. Celest. Mech. 1970, 2, 207–216. [Google Scholar] [CrossRef]
  19. Yoder, C.F.; Konopliv, A.S.; Yuan, D.; Standish, E.M.; Folkner, W.M. Fluid Core Size of Mars from Detection of the Solar Tide. Science 2003, 300, 299–303. [Google Scholar] [CrossRef]
  20. Konopliv, A.S.; Park, R.S.; Rivoldini, A.; Baland, R.; Maistre, S.L.; Hoolst, T.V.; Yseboodt, M.; Dehant, V. Detection of the Chandler Wobble of Mars From Orbiting Spacecraft. Geophys. Res. Lett. 2020, 47, e2020GL090568. [Google Scholar] [CrossRef]
  21. Goossens, S.; Matsumoto, K. Lunar degree 2 potential Love number determination from satellite tracking data. Geophys. Res. Lett. 2008, 35. [Google Scholar] [CrossRef]
  22. Mignard, F. The evolution of the lunar orbit revisited. I. Moon Planets 1979, 20, 301–315. [Google Scholar] [CrossRef]
  23. Efroimsky, M.; Lainey, V. Physics of Bodily Tides in Terrestrial Planets and the Appropriate Scales of Dynamical Evolution. J. Geophys. Res. 2007, 112. [Google Scholar] [CrossRef]
  24. Efroimsky, M.; Makarov, V.V. Tidal Friction and Tidal Lagging. Applicability Limitations of a Popular Formula for the Tidal Torque. Astrophys. J. 2012, 764, 10. [Google Scholar] [CrossRef]
  25. Emelyanov, N.V. Influence of tides in viscoelastic bodies of planet and satellite on the satellite’s orbital motion. Mon. Not. R. Astron. Soc. 2018, 479, 1278–1286. [Google Scholar] [CrossRef]
  26. Pou, L.; Nimmo, F.; Rivoldini, A.; Khan, A.; Bagheri, A.; Gray, T.R.; Samuel, H.; Lognonné, P.; Plesa, A.; Gudkova, T.; et al. Tidal Constraints on the Martian Interior. J. Geophys. Res. Planets 2022, 127, e2022JE007291. [Google Scholar] [CrossRef]
  27. Rambaux, N.; Castillo-Rogez, J.C.; Maistre, S.L.; Rosenblatt, P. Rotational motion of Phobos. Astron. Astrophys. 2012, 548, A14. [Google Scholar] [CrossRef]
  28. Yang, Y.; Yan, J.; Guo, X.; He, Q.; Barriot, J. An elastic model of Phobos’ libration. Astron. Astrophys. 2020, 636, A27. [Google Scholar] [CrossRef]
  29. Wisdom, J.; Peale, S.J.; Mignard, F. The chaotic rotation of Hyperion. Icarus 1984, 58, 137–152. [Google Scholar] [CrossRef]
  30. Pavlov, D.; Williams, J.G.; Suvorkin, V.V. Determining parameters of Moon’s orbital and rotational motion from LLR observations using GRAIL and IERS-recommended models. Celest. Mech. Dyn. Astron. 2016, 126, 61–88. [Google Scholar] [CrossRef]
  31. Breedlove, W.J. A Numerical Study of the Effects of Fourth Degree Terms in the Earth-Moon Mutual Potential on Lunar Physical Librations. In Proceedings of the a Symposium, Austin, TX, USA, 8–10 June 1976. [Google Scholar]
  32. Archinal, B.A.; Acton, C.H.; A’Hearn, M.F.; Conrad, A.; Consolmagno, G.J.; Duxbury, T.C.; Hestroffer, D.; Hilton, J.L.; Kirk, R.L.; Klioner, S.A.; et al. Report of the IAU Working Group on Cartographic Coordinates and Rotational Elements: 2015. Celest. Mech. Dyn. Astron. 2018, 130, 1–46. [Google Scholar] [CrossRef]
  33. Beletskiy, V.V. The Libration of a Satellite; NASA Technical Reports Server (NTRS): Washington, DC, USA, 1959; pp. 13–31. [Google Scholar]
  34. Peale, S.J. Orbital Resonances in the Solar System. Annu. Rev. Astron. Astrophys. 1976, 14, 215–246. [Google Scholar] [CrossRef]
  35. Goldreich, P.; Peale, S.J. Spin-orbit coupling in the solar system. Astron. J. 1966, 71, 425. [Google Scholar] [CrossRef]
  36. Ciufolini, I.; Pavlis, E.C. A confirmation of the general relativistic prediction of the Lense–Thirring effect. Nature 2004, 431, 958–960. [Google Scholar] [CrossRef] [PubMed]
  37. Iorio, L. On the Lense-Thirring test with the Mars Global Surveyor in the gravitational field of Mars. Cent. Eur. J. Phys. 2007, 8, 509–513. [Google Scholar] [CrossRef]
  38. Iorio, L. Juno, the angular momentum of Jupiter and the Lense–Thirring effect. J. End-to-End Test. 2008, 15, 554–560. [Google Scholar] [CrossRef]
Figure 1. The change of Euler angles with respect to the integration time.
Figure 1. The change of Euler angles with respect to the integration time.
Remotesensing 16 01174 g001
Figure 2. The change of Euler angle rates with respect to the integration time.
Figure 2. The change of Euler angle rates with respect to the integration time.
Remotesensing 16 01174 g002
Figure 3. The difference in position between the Deimos full dynamics model and simple model.
Figure 3. The difference in position between the Deimos full dynamics model and simple model.
Remotesensing 16 01174 g003
Figure 4. The difference in position in using the second-order gravity field and the third-order gravity field in the full model for Deimos.
Figure 4. The difference in position in using the second-order gravity field and the third-order gravity field in the full model for Deimos.
Remotesensing 16 01174 g004
Figure 5. The discrepancy between the precise orbit determination of the full dynamics model and the position of MARS097.
Figure 5. The discrepancy between the precise orbit determination of the full dynamics model and the position of MARS097.
Remotesensing 16 01174 g005
Figure 6. The discrepancy in XYZ coordinates between the precise orbit determination of the full dynamics model and the position of MARS097.
Figure 6. The discrepancy in XYZ coordinates between the precise orbit determination of the full dynamics model and the position of MARS097.
Remotesensing 16 01174 g006
Figure 7. The orientation of the pole for Deimos in the inertial frame in both the simple model and the full model.
Figure 7. The orientation of the pole for Deimos in the inertial frame in both the simple model and the full model.
Remotesensing 16 01174 g007
Figure 8. The angle between the minimum moment of inertia axis and the Mars position vector in the longitude direction in the complete model.
Figure 8. The angle between the minimum moment of inertia axis and the Mars position vector in the longitude direction in the complete model.
Remotesensing 16 01174 g008
Figure 9. The angle between the minimum moment of inertia axis and the Mars position vector in the latitude direction in the full model.
Figure 9. The angle between the minimum moment of inertia axis and the Mars position vector in the latitude direction in the full model.
Remotesensing 16 01174 g009
Table 1. Major physics parameters of Deimos.
Table 1. Major physics parameters of Deimos.
ParametersSymbolValue
Mass parameters GM m 3 / s 2 9.68 × 10 4
Radius R   ( m ) 6250
Second-degree and order tesseral harmonic C 22 4.773 × 10 2
Second-degree zonal harmonic J 2 0.108
The moments of inertia A 0.338
B 0.461
C 0.508
Table 2. Normalized gravitational coefficients of Deimos in units of 10 2 .
Table 2. Normalized gravitational coefficients of Deimos in units of 10 2 .
( l , m ) C n m ( × 10 2 ) S n m ( × 10 2 )
(1, 0)00
(1, 1)00
(2, 0)−4.8270
(2, 1)−0.07270.362
(2, 2)4.7730.122
(3, 0)0.9700
(3, 1)1.4300.0521
(3, 2)−1.108−0.499
(3, 3)−0.1210.728
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Huang, K.; Zhang, L.; Yang, Y.; Ye, M.; Li, Y. Dynamical Model of Rotation and Orbital Coupling for Deimos. Remote Sens. 2024, 16, 1174. https://doi.org/10.3390/rs16071174

AMA Style

Huang K, Zhang L, Yang Y, Ye M, Li Y. Dynamical Model of Rotation and Orbital Coupling for Deimos. Remote Sensing. 2024; 16(7):1174. https://doi.org/10.3390/rs16071174

Chicago/Turabian Style

Huang, Kai, Lijun Zhang, Yongzhang Yang, Mao Ye, and Yuqiang Li. 2024. "Dynamical Model of Rotation and Orbital Coupling for Deimos" Remote Sensing 16, no. 7: 1174. https://doi.org/10.3390/rs16071174

APA Style

Huang, K., Zhang, L., Yang, Y., Ye, M., & Li, Y. (2024). Dynamical Model of Rotation and Orbital Coupling for Deimos. Remote Sensing, 16(7), 1174. https://doi.org/10.3390/rs16071174

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