Lunar surface temperature is one of the fundamental thermophysical parameters of the lunar regolith, which is of great significance to the interpretation of remote-sensing thermal data. In this study, a daytime surface temperature model is established focusing on the lunar superficial layer with high spatial-temporal resolution. The physical parameters at the time of interest are adopted, including effective solar irradiance, lunar libration, large-scale topographic shading, and surrounding diffuse reflection. Thereafter, the 1/64° temperature distributions at five local times are quantitatively generated and analyzed in Sinus Iridum. Also, combined with Chang’E-2 microwave radiometer (CELMS) data and Diviner thermal infrared (TIR) data, the spectral emissivity distributions are estimated as a potential geological application of the simulated surface temperature. The results are as follows: (1) daytime surface temperature in Sinus Iridum is significantly affected by the local topography and observation time, and the influence of diffuse reflection energy is obvious; (2) the emissivity distributions provide a new way to understand the thermophysical properties difference of lunar regolith at different depths; (3) the influence of lunar orbiting revolution and precession on surface temperature should be analyzed carefully, which shows the importance of using the parameters at the time of interest.

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.

Context. Study the rotation of a celestial body is an efficient way to infer its interior structure, and then may give information of its origin and evolution. In this study, based on the latest shape model of Phobos from Mars Express (MEX) mission, the polyhedron approximation approach was used to simulate the gravity field of Phobos. Then, the gravity information was combined with the newest geophysical parameters such as GM and k₂ to construct the numerical model of Phobos’ rotation. And with an appropriate angles transformation, we got the librational series respect to Martian mean equator of date. Aims. The purpose of this paper is to develop a numerical model of Phobos’ rotational motion that includes the elastic properties of Phobos. The frequencies analysis of the librational angles calculated from the numerical integration results emphasize the relationship between geophysical properties and dynamics of Phobos. This work will also be useful for a future space mission dedicated to Phobos. Methods. Based on the latest shape model of Phobos from MEX mission, we firstly modeled the gravity field of Phobos, then the gravity coefficients were combined with some of the newest geophysical parameters to simulate the rotational motion of Phobos. To investigate how the elastic properties of Phobos affect its librational motion, we adopted various k₂ into our numerical integration. Then the analysis was performed by iterating a frequency analysis and linear least-squares fit of Phobos’ physical librations. From this analysis, we identified the influence of k₂ on the largest librational amplitude and its phase. Results. We showed the first ten periods of the librational angles and found that they agree well with the previous numerical results which Phobos was treated as a perfectly rigid body. We also found that the maximum amplitudes of the three parameters of libration are also close to the results from a rigid model, which is mainly due to the inclination of Phobos and moments of inertia. The other amplitudes are slightly different, since the physics contained in our model is different to that of a previous study, specifically, the different low-degree gravity coefficients and ephemeris. The libration in longitude τ has the same quadratic term with previous numerical study, which is consistent with the secular acceleration of Phobos falling onto Mars. We investigated the influence of the tidal Love number k₂ on Phobos’ rotation and found a detectable amplitude changes (0.0005°) expected in the future space mission on τ, which provided a potential possibility to constrain the k₂ of Phobos by observing its rotation. We also studied the influence of Phobos’ orbit accuracy on its libration and suggested a simultaneous integration of orbit and rotation in future work.

This paper primarily investigates the effect of the tilt of corner cube reflector (CCR) arrays on lunar laser ranging (LLR). A mathematical model was established to study the random errors caused by the tilt of the CCR arrays. The study found that, ideally, when the laser ranging pulse width is 10 picoseconds or less, it is possible to distinguish from which specific corner cubes within the CCR array each peak in the echo signal originates. Consequently, partial data from the echo can be extracted for signal processing, significantly reducing random errors and improving the single-shot precision of LLR. The distance obtained by extracting part of the echo can be reduced to the center position of the array, thereby providing multiple higher-precision ranging results from each measurement. This not only improves the precision of LLR but also increases the data volume. A simulation experiment based on the 1.2 m laser ranging system at Yunnan Observatories was conducted. By extracting one peak for signal processing, the single-shot precision improved from 32.24 mm to 2.52 mm, validating the theoretical analysis results. Finally, an experimental laser ranging system based on a 53 cm binocular telescope system was established for ground experiments. The experimental results indicated that the echo signal could identify the tilt state of the CCR array. By extracting the peak returned by the central CCR for signal processing, the ranging precision was greatly improved. Through theoretical analyses, simulation experiments, and ground experiments, a solution to reduce the random errors caused by the tilt of the CCR array was provided. This offers an approach to enhance the single-shot precision of future LLR and provides a reference for upgrading ground-based equipment at future laser ranging stations.

Context. High-precision ephemerides are not only useful in supporting space missions, but also in investigating the physical nature of celestial bodies. This paper reports an update to the orbit and rotation model of the Martian moon Phobos. In contrast to earlier numerical models, this paper details a dynamical model that fully considers the rotation of Phobos. Here, Phobos’ rotation is first described by Euler’s rotational equations and integrated simultaneously with the orbital motion equations. We discuss this dynamical model, along with the differences with respect to the model now in use. Aims. This work is aimed at updating the physical model embedded in the ephemerides of Martian moons, considering improvements offered by exploiting high-precision observations expected from future missions (e.g., Japanese Martian Moons exploration, MMX), which fully supports future studies of the Martian moons. Methods. The rotational motion of Phobos can be expressed by Euler’s rotational equations and integrated in parallel with the equations of the orbital motion of Phobos around Mars. In order to investigate the differences between the two models, we first reproduced and simulated the dynamical model that is now used in the ephemerides, but based on our own parameters. We then fit the model to the newest Phobos ephemeris published by Institut de Mécanique Céleste et de Calcul des Éphémérides (IMCCE). Based on our derived variational equations, the influence of the gravity field, the Love number, k₂, and the rotation behavior were studied by fitting the full model to the simulated simple model. Our revised dynamic model for Phobos was constructed as a general method that can be extended with appropriate corrections (mainly rotation) to systems other than Phobos, such as the Saturn and Jupiter systems. Results. We present the variational equation for Phobos’ rotation employing the symbolic Maple computation software. The adjustment test simulations confirm the latitude libration of Phobos, suggesting gravity field coefficients obtained using a shape model and homogeneous density hypothesis should be re-examined in the future in the context of dynamics. Furthermore, the simulations with different k₂ values indicate that it is difficult to determine k₂ efficiently using the current data.