Andrei A. Dmitrovskii

Ocean tide loading (OTL) brings about recurring deformation of the Earth’s surface. Some of the OTL harmonics, e.g. M2, O1, Mf, cause sufficiently large surface displacement to be registered by the Global Navigation Satellite Systems (GNSS). These displacements are sensitive to the interior structure of the planet in a broad range of temporal and spatial scales making them a potentially unique source of information about the planet’s response at low frequencies. Comparison between observations and predictions for 1D elastic Earth models result in discrepancies of up to 3 mm (Bos et al., 2015, Martens et al., 2016). Spatial coherency of these discrepancies hints to 3D interior structure as one of the main sources of such residuals.In this context, we present a framework to invert OTL observations for 3D crustal and mantle structure based on a trust-region Newton-type iterative algorithm. Furthermore, we resort to the adjoint approach as an efficient means of computing the gradient for the high-dimensional model space. Focusing on the design of the inverse algorithm, we constrain ourselves to deformations of an isotropic elastic planet, which are governed by a self-adjoint forward operator. In order to assess the robustness of the method, we perform a suite of 3D synthetic inversions that mimic the distribution of the GNSS stations in South America. Preliminary results indicate enhanced sensitivities to the crustal and upper mantle density and elastic properties in the vicinity of the coastlines.

Deformations of the solid Earth as a response to ocean tidal loading (OTL) are sensitive to the material properties of Earth’s interior across a broad range of spatial and temporal scales. Studying tidal response can provide constraints on the interior structure, which are complementary to seismic tomography and particularly important to explore the interior response to low frequency loads. Although seismic tomography is widely used to constrain the Earth’s interior, it is prone to be only slightly sensitive to the density distribution in the interior with an increase of the sensitivity towards the long period signal. Whereas previous research (e.g. Ito & Simons, 2011, Martens et al., 2016) has shown that the tidal surface displacements may be sensitive to elastic properties of the interior to the same extent as to the mass distribution in the lithosphere and the mantle. The latter are of massive interest to all fields of geophysics and especially geodynamics. We present a numerical approach to simulate the elastic and gravitational responses of the solid Earth that relies on the spectral-element method. Modeling the governing equations in a 3-D Earth using a coupled system of the elastostatic and Poisson’s equations enables us to include effects like topography or lateral variations in Earth structure. The adjoint method is a powerful technique to simultaneously compute sensitivity with respect to all material parameters, e.g., density and elastic moduli, by solving an auxiliary linear system. We introduce a recipe for computing adjoint-based sensitivities of the complex-valued amplitude of surface displacement by two simulations for the real and imaginary part of the surface load. Those two simulations are independent under assumption of negligible attenuation.