A finite element analysis of tidal deformation of the entire earth with a discontinuous outer layer

Tidal deformation of the Earth is normally calculated using the analytical solution with some simplified assumptions, such as the Earth is a perfect sphere of continuous media. This paper proposes an alternative way, in which the Earth crust is discontinuous along its boundaries, to calculate the tidal deformation using a finite element method. An in-house finite element code is firstly introduced in brief and then extended here to calculate the tidal deformation. The tidal deformation of the Earth due to the Moon was calculated for an geophysical earth model with the discontinuous outer layer and compared with the continuous case. The preliminary results indicate that the discontinuity could have different effects on the tidal deformation in the local zone around the fault, but almost no effects on both the locations far from the fault and the global deformation amplitude of the Earth. The localized deformation amplitude seems to depend much on the relative orientation between the fault strike direction and the loading direction (i.e. the location of the Moon) and the physical property of the fault.     

Some remarks about the degree-one deformation of the Earth

The degree-one deformation of the Earth (and the induced discrepancy between the figure centre and the mass centre of the Earth) is computed using a theoretical approach (Love numbers formalism) at short timescales (where the Earth has an elastic behaviour) as well as at long timescales (where the Earth has a viscoelastic or quasi-fluid behaviour). For a Maxwell model of rheology, the degree-one relaxation modes associated with the viscoelastic Love numbers have been investigated: the M_o mode does not exist and there is only one transition mode (instead of two) generated by a viscosity discontinuity.
The translations at each interface of the incompressible layers of the earth model [surface, 670 km depth discontinuity, core-mantle boundary (CMB) and inner-core boundary (ICB)] are computed. They are elastic with an order of magnitude of about 1 mm when the excitation source is the atmospheric continental loading or a magnetic pressure acting at the CMB. They are viscoelastic when the earth is submitted to Pleistocene deglaciation, with an order of magnitude of about 1 m. In a quasi-fluid approximation (Newtonian fluid) because of the mantle density heterogeneity their order of magnitude is about 100 m (except for the ICB, which is in quasi-hydrostatic equilibrium at this timescale).     

On azimuthal eigenwavenumbers associated with Laplace's tidal equations

Laplace's tidal equations for the case of an ocean of constant depth bounded by meridians were considered by two authors at a specific frequency as an eigenvalue problem in the azimuthal wavenumber. A finite spectrum of eigenwavenumbers was found. That eigenvalue problem is re-examined by means of asymptotic techniques and numerical integration of the governing equation of the problem. At low frequencies a formula connecting the frequency and the number of eigensolutions is established. It is shown that at a given frequency the spectrum of eigenwavenumbers is wider than that reported, but (for this type of solution) the meridional boundary conditions are satisfied approximately only for the case of very low frequencies.