GEODYNAMICS

Purpose. To investigate the features of the algorithm implementation for finding the derivatives of the spatial distribution function of the planet's masses with the use of high-order Stokes constants and, on the basis of this, to find its analytical expression. According to the given methodology, to carry out calculations with the help of which to carry on the study of dynamic phenomena occurring inside an ellipsoidal planet. The proposed method involves the determination of the derivatives of the mass distribution function by the sum, the coefficients of which are obtained from the system of equations, which is incorrect. In order to solve it, an error-resistant method for calculating unknowns was used. The implementation of the construction is carried out in an iterative way, while for the initial approximation we take the three-dimensional function of the density of the Earth's masses, built according to Stokes constants up to the second order inclusive, by dynamic compression by the one-dimensional density distribution, and we determine the expansion coefficients of the derivatives of the function in the variables to the third order inclusive. They are followed by the corresponding density function, which is then taken as the initial one. The process is repeated until the specified order of approximation is reached. To obtain a stable result, we use the Cesaro summation method (method of means).. The calculations performed with the help of programs that implement the given algorithm, while the achieved high (ninth) order of obtaining the terms of the sum of calculations. The studies of the convergence of the sum of the series have been carried out, and on this basis, a conclusion has been made about the advisability of using the generalized finding of the sums based on the Cesaro method. The optimal number of contents of the sum terms has been chosen, provides convergence both for the mass distribution function and for its derivatives. Calculations of the deviations of mass distribution from the mean value ("inhomogeneities") for extreme points of the earth's geoid, which basically show the total compensation along the radius of the Earth, have been performed. For such three-dimensional distributions, calculations were performed and schematic maps were constructed according to the taken into account values of deviations of three-dimensional distributions of the mean ("inhomogeneities") at different depths reflecting the general structure of the Earth's internal structure. The presented vector diagrams of the horizontal components of the density gradient at characteristic depths (2891 km - core-mantle, 700 km - middle of the mantle, also the upper mantle - 200, 100 km) allow us to draw preliminary conclusions about the global movement of masses. At the same time, a closed loop is observed on the “core-mantle” edge, which is an analogy of a closed electric circuit. For shallower depths, differentiation of vector motions is already taking place, which gives hope for attracting these vector-grams to the study of dynamic motions inside the Earth. In fact, the vertical component (derivative with respect to the z variable) is directed towards the center of mass and confirms the main property of mass distributions - growth when approaching the center of mass. The method of stable solution of incorrect linear systems is applied, by means of which the vector-gram of the gradient of the mass distribution function is constructed. The nature of such schemes provides a tool for possible causes of mass redistribution in the middle of the planet and to identify possible factors of tectonic processes in the middle of the Earth, i.e indirectly confirms the gravitational convection of masses. The proposed technique can be used to create detailed models of density functions and its characteristics (derivatives) of the planet's interior, and the results of numerical experiments - to solve tectonics problems.

Population Properties of Neutron Stars in the Coalescing Compact Binaries

We perform a hierarchical Bayesian inference to investigate the population properties of the coalescing compact binaries involving at least one neutron star (NS). With the current gravitational-wave (GW) observation data, we can rule out none of the double Gaussian, single Gaussian, and uniform NS mass distribution models, though a specific double Gaussian model inferred from the Galactic NSs is found to be slightly more preferred. The mass distribution of black holes (BHs) in the neutron star–black hole (NSBH) population is found to be similar to that in the Galactic X-ray binaries. Additionally, the ratio of the merger rate densities between NSBHs and BNSs is estimated to be ∼3:7. The spin properties of the binaries, though constrained relatively poorly, play a nontrivial role in reconstructing the mass distribution of NSs and BHs. We find that a perfectly aligned spin distribution can be ruled out, while a purely isotropic distribution of spin orientation is still allowed. To evaluate the feasibility of reliably determining the population properties of NSs in the coalescing compact binaries with upcoming GW observations, we perform simulations with a mock population. We find that with 100 detections (including BNSs and NSBHs) the mass distribution of NSs can be well determined, and the fraction of BNSs can also be accurately estimated.

MMX geodesy investigations: science requirements and observation strategy

In order to investigate the origin of Phobos and Deimos, the Japanese Martian Moons eXploration (MMX) mission is scheduled for launch in 2024. MMX will make comprehensive remote-sensing measurements of both moons and return regolith samples from Phobos to Earth. Geodetic measurements of gravity, shape, and rotation parameter of a body provides constraints on its internal structure reflecting its origin and evolution. Moments of inertia are important parameters to constrain the internal mass distribution, but they have not been well determined for the Martian moons yet. We discuss the mission requirements related to the moments of inertia to detect a potential heterogeneity of the mass distribution inside Phobos. We introduce mission instruments and operational strategies to meet the mission requirements. We present a preliminary imaging strategy from a quasi-satellite orbit for a base shape model that is expected to be created at the early stage of the mission. Geodetic products including ephemeris, gravity field, rotation parameter of Phobos, and spacecraft orbit are of importance not only for the geodetic study, but also for interpreting data from various mission instruments and selecting possible landing sites. Graphical Abstract