Physics – Geophysics
Scientific paper
May 2010
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=2010eguga..12.6660z&link_type=abstract
EGU General Assembly 2010, held 2-7 May, 2010 in Vienna, Austria, p.6660
Physics
Geophysics
Scientific paper
We deal with the equations and boundary conditions describing deformation and gravitational potential of prestressed spherically symmetric elastic bodies by decomposing governing equations into a series of boundary value problems (BVP) for ordinary differential equations (ODE) of the second order. In contrast to traditional Runge-Kutta integration techniques, highly accurate pseudospectral schemes are employed to directly discretize the BVP on Chebyshev grids and a set of linear algebraic equations with an almost block diagonal matrix is derived. As a consequence of keeping the governing ODEs of the second order instead of the usual first-order equations, the resulting algebraic system is half-sized but derivatives of the model parameters are required. Moreover, they can be easily evaluated for models, where structural parametres are piecewise polynomially dependent. Both accuracy and efficiency of the method are tested by evaluating the tidal Love numbers for the Earth's model PREM. Finally, we also derive complex Love numbers for models with the Maxwell viscoelastic rheology, where viscosity is a depth-dependent function. The method is applied to evaluation of the tidal Love numbers for models of Mars and Venus. The Love numbers of the two Martian models - the former optimized to cosmochemical data and the latter to the moment of inertia (Sohl and Spohn, 1997) - are h2=0.172 (0.212) and k2=0.093 (0.113). For Venus, the value of k2=0.295 (Konopliv and Yoder, 1996), obtained from the gravity-field analysis, is consistent with the results for our model with the liquid-core radius of 3110 km (Zábranová et al., 2009). Together with rapid evaluation of free oscillation periods by an analogous method, this combined matrix approach could by employed as an efficient numerical tool in structural studies of planetary bodies. REFERENCES Konopliv, A. S. and Yoder, C. F., 1996. Venusian k2 tidal Love number from Magellan and PVO tracking data, Geophys. Res. Lett., 23, 1857-1860. Sohl, F., and Spohn, T., 1997. The interior structure of Mars: Implications from SNC meteorites, J. Geophys. Res., 102, 1613-1635. Zabranova, E., Hanyk L. and Matyska, C.: Matrix Pseudospectral Method for Elastic Tides Modeling. In: Holota P. (Ed.): Mission and Passion: Science. A volume dedicated to Milan Bursa on the occasion of his 80th birthday. Published by the Czech National Committee of Geodesy and Geophysics. Prague, 2009, pp. 243-260.
Hanyk Ladidslav
Matyska Ctirad
Zabranova Eliska
No associations
LandOfFree
Matrix Pseudospectral Method for (Visco)Elastic Tides Modeling of Planetary Bodies does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Matrix Pseudospectral Method for (Visco)Elastic Tides Modeling of Planetary Bodies, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Matrix Pseudospectral Method for (Visco)Elastic Tides Modeling of Planetary Bodies will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-917584