Theoretical calculation of tidal Love numbers of the Moon with a new spectral element method

The tidal Love numbers of the Moon are a set of nondimensional parameters that describe the deformation responses of the Moon to the tidal forces of external celestial bodies.They play an important role in the theoretical calculation of the Moon’s tidal deformation and the inversion of its internal structure.In this study,we introduce the basic theory for the theoretical calculation of the tidal Love numbers and propose a new method of solving the tidal Love numbers:the spectral element method.Moreover,we explain the mathematical theory and advantages of this method.On the basis of this new method,using 10 published lunar internal structure reference models,the lunar surface and lunar internal tidal Love numbers were calculated,and the influence of different lunar models on the calculated Love numbers was analyzed.Results of the calculation showed that the difference in the second-degree lunar surface Love numbers among different lunar models was within 8.5%,the influence on the maximum vertical displacement on the lunar surface could reach±8.5 mm,and the influence on the maximum gravity change could reach±6μGal.Regarding the influence on the Love numbers inside the Moon,different lunar models had a greater impact on the Love numbers h_(2) and l_(2) than on k_(2) in the lower lunar mantle and core.

Three-Dimensional Multiferroic Structures under Time-Harmonic Loading

Magneto-electro-elastic(MEE)materials are a specific class of advanced smart materials that simultaneouslymanifest the coupling behavior under electric,magnetic,and mechanical loads.This unique combination ofproperties allows MEE materials to respond to mechanical,electric,and magnetic stimuli,making them versatile forvarious applications.This paper investigates the static and time-harmonic field solutions induced by the surface loadin a three-dimensional(3D)multilayered transversally isotropic(TI)linear MEE layered solid.Green’s functionscorresponding to the applied uniform load(in both horizontal and vertical directions)are derived using the FourierBessel series(FBS)system of vector functions.By virtue of this FBS method,two sets of first-order ordinarydifferential equations(i.e.,N-type and LM-type)are obtained,with the expansion coefficients being Love numbers.It is noted that the LM-type system corresponds to the MEE-coupled P-,SV-,and Rayleigh waves,while the N-typecorresponds to the purely elastic SH-and Love waves.By applying the continuity conditions across interfaces,the solutions for each layer of the structure(from the bottom to the top)are derived using the dual-variable andposition(DVP)method.This method(i.e.,DVP)is unconditionally stable when propagating solutions throughdifferent layers.Numerical examples illustrate the impact of load types,layering,and frequency on the response ofthe structure,as well as the accuracy and convergence of the proposed approach.The numerical results are usefulin designing smart devices made of MEE solids,which are applicable to engineering fields like renewable energy.

Time-space characteristics of viscoelastic post-seismic deformations corresponding to different rheology models

On a long time(>1 a)scale,the viscoelastic properties of mantle media significantly affect post-seismic deformation.The stress field disturbance in viscoelastic medium caused by fault slip gradually relax,and the relaxation process and its temporal-spatial characteristics are determined by the viscoelastic model.In this paper,we assume that the mantle media are types of common linear rheological models,i.e.,the Burgers body,the standard linear solid,and the Maxell body,and we calculate the dislocation Love number and Green function for a spherically symmetric,non-rotating,viscoelastic,and isotropic(SNRVEI)Earth model.The characteristics of the post-seismic relaxation deformations corresponding to the different models are compared.Our results show that for a short time period,the Burgers body and standard linear solid are similar;while for the long time period,the Burgers body and Maxwell body are similar.This suggests that the observations of post-seismic deformation on the surface have a great potential for the inversion of underground viscoelastic structures.However,the potential of using surface displacement to distinguish different rheological models is limited when the observation period is not long enough.

Characteristics of gravity signal and loading effect in China

The complex geographical environment in China makes its gravity signals miscellaneous.This work gives a comprehensive representation and explanation in secular trend of gravity change in different regions,the key features of which include positive trend in inner Tibet Plateau and South China and negative trend in North China plain and high mountain Asia（HMA）.We also present the patterns of amplitudes and phases of annual and semiannual change.The mechanism underlying the semiannual period is explicitly discussed.The displacement in three directions expressed in terms of geo-potential spherical coefficients and load Love numbers are given.A case study applied with these equations is presented.The results show that Global Positioning System（GPS） observations can be used to compare with Gravity Recovery and Climate Experiment（GRACE） derived displacement and the vertical direction has a signal-noise-ratio of about one order of magnitude larger than the horizontal directions.

Shape and gravitational field of the ellipsoidal satellites

The shape and gravitational field of ellipsoidal satellites are studied by using the tidal theory. For ellipsoidal satellites, the following conclusions were obtained: Firstly, in the early stage of the satellite formation, strong tidal friction allowed the satellites move in a synchronous orbit and evolve into a triaxial ellipsoidal shape. Because the tidal potential from the associated primary and the centrifugal potential from the satellite spin are nearly fixed at the surface, the early satellites are the viscoelastic celestial body, and their surfaces are nearly in the hydrostatic equilibrium state. The deformation is fixed in the surface of the satellite. By using the related parameters of primary and satellite, the tidal height and the theoretical lengths of three primary radii of the ellipsoidal satellite are calculated. Secondly, the current ellipsoidal satellites nearly maintain their ellipsoidal shape from solidification, which happened a few billion years ago. According to the satellite shape, we estimated the orbital period and spinning angular velocity, and then determined the evolution of the orbit. Lastly, assuming an ellipsoidal satellite originated in the hydrostatic equilibrium state, the surface shape could be determined by tidal, rotation, and additional potentials. However, the shape of the satellite's geoid differs from its surface shape. The relationship between these shapes is discussed and a formula for the gravitational harmonic coefficients is presented.