Search

In global studies investigating the Earth's lithospheric structure, the spectral expressions for the gravimetric forward and inverse modeling of the global gravitational and crustal structure models are preferably used, because of their numerical efficiency. In regional studies, the applied numerical schemes typically utilize the expressions in spatial form. Since the gravity-gradient observations have a more localized support than the gravity measurements, the gravity-gradient data (such as products from the Gravity field and steady-state Ocean Circulation Explorer - GOCE - gravity-gradiometry satellite mission) could preferably be used in regional studies, because of reducing significantly the spatial data-coverage required for a regional inversion or interpretation. In this study, we investigate this aspect in context of a regional Moho recovery. In particular, we compare the numerical performance of solving the Vening Meinesz-Moritz's (VMM) inverse problem of isostasy in spectral and spatial domains from the gravity and (vertical) gravity-gradient data. We demonstrate that the VMM spectral solutions from the gravity and gravity-gradient data are (almost) the same, while the VMM spatial solutions differ from the corresponding spectral solutions, especially when using the gravity-gradient data. The validation of the VMM solutions, however, reveals that the VMM spatial solution from the gravity-gradient data has a slightly better agreement with seismic models. A more detailed numerical analysis shows that the VMM spatial solution formulated for the gravity gradient is very sensitive to horizontal spatial variations of the vertical gravity gradient, especially in vicinity of the computation point. Consequently, this solution provides better results in regions with a relatively well-known crustal structure, while suppressing errors caused by crustal model uncertainties from distant zones. Based on these findings we argue that the gravity-gradient data are more suitable than the gravity data for a regional Moho recovery. [ABSTRACT FROM AUTHOR]

In this study we investigate the lithospheric stresses computed from the gravity and lithospheric structure models. The functional relation between the lithospheric stress tensor and the gravity field parameters is formulated based on solving the boundary-value problem of elasticity in order to determine the propagation of stresses inside the lithosphere, while assuming the horizontal shear stress components (computed at the base of the lithosphere) as lower boundary values for solving this problem. We further suppress the signature of global mantle flow in the stress spectrum by subtracting the long-wavelength harmonics (below the degree of 13). This numerical scheme is applied to compute the normal and shear stress tensor components globally at the Moho interface. The results reveal that most of the lithospheric stresses are accumulated along active convergent tectonic margins of oceanic subductions and along continent-to-continent tectonic plate collisions. These results indicate that, aside from a frictional drag caused by mantle convection, the largest stresses within the lithosphere are induced by subduction slab pull forces on the side of subducted lithosphere, which are coupled by slightly less pronounced stresses (on the side of overriding lithospheric plate) possibly attributed to trench suction. Our results also show the presence of (intra-plate) lithospheric loading stresses along Hawaii islands. The signature of ridge push (along divergent tectonic margins) and basal shear traction resistive forces is not clearly manifested at the investigated stress spectrum (between the degrees from 13 to 180). [ABSTRACT FROM AUTHOR]

A regional recovery of the Earth's gravity field from satellite observables has become particularly important in various geoscience studies in order to better localize stochastic properties of observed data, while allowing the inversion of a large amount of data, collected with a high spatial resolution only over the area of interest. One way of doing this is to use observables, which have a more localized support. As acquired in recent studies related to a regional inversion of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) data, the satellite gravity-gradient observables have a more localized support than the gravity observations. Following this principle, we compare here the performance of the second- and third-order derivatives of the gravitational potential in context of a regional gravity modeling, namely estimating the gravity anomalies. A functional relation between these two types of observables and the gravity anomalies is formulated by means of the extended Stokes' integral formula (or more explicitly its second- and third-order derivatives) while the inverse solution is carried out by applying a least-squares technique and the ill-posed inverse problem is stabilized by applying Tikhonov's regularization. Our results reveal that the third-order radial derivatives of the gravitational potential are the most suitable among investigated input data types for a regional gravity recovery, because these observables preserve more information on a higher-frequency part of the gravitational spectrum compared to the vertical gravitational gradients. We also demonstrate that the higher-order horizontal derivatives of the gravitational potential do not necessary improve the results. We explain this by the fact that most of the gravity signal is comprised in its radial component, while the horizontal components are considerably less sensitive to spatial variations of the gravity field. [ABSTRACT FROM AUTHOR]

We develop and apply the algorithm for a regional gravimetric Moho recovery in spatial domain. The functional relation between the (known) refined gravity field and the (unknown and sought) Moho depth is defined by means of solving the Fredholm integral equation of the first kind. Linearization is applied to define the Moho depth corrections with respect to the mean Moho depth. The system of linearized observation equations is solved to find the Moho depth corrections, and Tikhonov’s regularization is applied to stabilize this ill-posed inverse problem. Developed algorithm is applied to model the Moho depth regionally under the Tibetan Plateau and Himalayas, while adopting the uniform and variable Moho density contrast models, and gravimetric results are validated using the CRUST1.0 seismic model. Our results show a relatively good agreement between gravimetric and seismic models, without presence of a significant systematic bias. Our result, however, indicates that for a regional study the variable Moho density contrast might not improve the results especially when density structure of the lower crust and uppermost mantle is not know accurately. In that case, the use of the uniform Moho density contrast for a regional Moho recovery is more appropriate. [ABSTRACT FROM AUTHOR]

To improve the numerical efficiency, the Fast Fourier Transform (FFT) technique was facilitated in Parker-Oldenburg's method for a regional gravimetric Moho recovery, which assumes the Earth's planar approximation. In this study, we extend this definition for global applications while assuming a spherical approximation of the Earth. In particular, we utilize the FFT technique for a global Moho recovery, which is practically realized in two numerical steps. The gravimetric forward modeling is first applied, based on methods for a spherical harmonic analysis and synthesis of the global gravity and lithospheric structure models, to compute the refined gravity field, which comprises mainly the gravitational signature of the Moho geometry. The gravimetric inverse problem is then solved iteratively in order to determine the Moho depth. The application of FFT technique to both numerical steps reduces the computation time to a fraction of that required without applying this fast algorithm. The developed numerical producers are used to estimate the Moho depth globally, and the gravimetric result is validated using the global (CRUST1.0) and regional (ESC) seismic Moho models. The comparison reveals a relatively good agreement between the gravimetric and seismic models, with the RMS of differences (of 4-5 km) at the level of expected uncertainties of used input datasets, while without the presence of significant systematic bias. [ABSTRACT FROM AUTHOR]

We investigate the horizontal stress field beneath crustal structures of central Eurasia. The numerical procedure applied for a simultaneous determination of the sub-crustal stress and the crustal thickness from the global gravity, terrain and crustal structure models is based on solving Navier-Stokes' problem which incorporates the inverse solution to the Vening Meinesz- Moritz's problem of isostasy. The numerical results reveal that a spatial distribution of the sub-crustal stress in this study area closely resembles the regional tectonic configuration comprising parts of the Eurasian, Indian and Arabian lithospheric plates. The maximum shear stress intensity is generated by a subduction of the Indian plate beneath the Tibetan block. The intra-plate tectonic configuration is marked by the stress anomalies distributed along major active strike-slip fault systems and sections of subduction which separate the Tibetan and Iranian blocks from the rest of the Eurasian plate. The most pronounced intra-plate stress anomalies are related with a subduction of the Eurasian plate beneath the Tibetan block. We also demonstrate that a prevailing convergent orientation of stress vectors agree with the compressional tectonism of orogenic formations (Himalaya and Tibet Plateau, Than Shan, Zargos and Iranian Plateau), while the extensional tectonism of continental basins (Tarim, Ganges-Brahmaputra, Sichuan) is manifested by a divergence of stress vectors. [ABSTRACT FROM AUTHOR]

The occurrence of subduction earthquakes usually leads to considerable localized mass migration changes. To improve the detection of earthquake subductions as well as the constraint of fault parameters, we derive expressions that describe changes of the gravitational curvatures (GC), i.e., the third-order derivatives of the Earth's gravitational potential, caused by a point dislocation while adopting a spherical symmetric Earth model. As a 3-D tensor matrix, the GC have twenty-seven components of those seven are independent. First, we investigate the dislocation Love numbers of the Earth's gravitational potential and derive the Green's functions of GC caused by four independent point sources in a spherical inhomogeneous Earth model. We then present the GC changes in a half-space Earth model. Furthermore, we conduct a sensitivity study by using three physical quantities that involve gravitation, gravitational gradients, and GC to compare their abilities in a seismic source depth detection. Our numerical results reveal that changes in the GC are more sensitive to a medium information about the field source compared to gravitation and gravitational gradients. This finding indicates that GC measurements could provide a more detailed information about slip fault parameters when considering a heterogeneous slip. Despite a widespread application of gravity gradients in Earth science, especially after launching the Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) satellite mission, measurements of the third-order derivatives of the gravitational potential have an enormous potential in the study of the solid Earth, although further work is needed in terms of instrument design and development. [Display omitted] • Expressions describing changes of the third-order derivatives of the Earth's gravitational potential caused by point dislocation are derived. • Changes in third-order derivatives of the Earth's gravitational potential provide a more detailed information on slip fault parameters. [ABSTRACT FROM AUTHOR]

The results of global and regional studies often show significant disagreement between the Moho depths determined using seismic and isostatic models. In this study, we estimate the differences between these two models in central Eurasia. The Vening Meinesz-Moritz (VMM) inverse problem of isostasy is utilized to determine the isostatic Moho depths. The estimated VMM Moho depths are then corrected for the sediment density contrast. The application of this correction improves the agreement between the isostatic and seismic Moho models. The existing discrepancies between the isostatic and seismic models are finally modeled by applying the non-isostatic correction, which accounts for the unmodelled mantle density heterogeneities and other geodynamic processes, which are not taken into account in classical isostatic models. Our results reveal that the non-isostatic correction still cannot fully describe mechanisms affecting the Moho geometry along the convergent continent-tocontinent tectonic plate boundaries occurring beneath Himalayas despite an overall good performance of the applied method. [ABSTRACT FROM AUTHOR]

Global forward modelling of the Earth's gravitational potential, a classical problem in geophysics and geodesy, is relevant for a range of applications such as gravity interpretation, isostatic hypothesis testing or combined gravity field modelling with high and ultra-high resolution. This study presents spectral forward modelling with volumetric mass layers to degree 2190 for the first time based on two different levels of approximation. In spherical approximation, the mass layers are referred to a sphere, yielding the spherical topographic potential. In ellipsoidal approximation where an ellipsoid of revolution provides the reference, the ellipsoidal topographic potential (ETP) is obtained. For both types of approximation, we derive a mass layer concept and study it with layered data from the Earth2014 topography model at 5-arc-min resolution. We show that the layer concept can be applied with either actual layer density or density contrasts w.r.t. a reference density, without discernible differences in the computed gravity functionals. To avoid aliasing and truncation errors, we carefully account for increased sampling requirements due to the exponentiation of the boundary functions and consider all numerically relevant terms of the involved binominal series expansions. The main outcome of our work is a set of new spectral models of the Earth's topographic potential relying on mass layer modelling in spherical and in ellipsoidal approximation. We compare both levels of approximations geometrically, spectrally and numerically and quantify the benefits over the frequently used rock-equivalent topography (RET) method. We show that by using the ETP it is possible to avoid any displacement of masses and quantify also the benefit of mapping-free modelling. The layer-based forward modelling is corroborated by GOCE satellite gradiometry, by in-situ gravity observations from recently released Antarctic gravity anomaly grids and degree correlations with spectral models of the Earth's observed geopotential. As the main conclusion of this work, the mass layer approach allows more accurate modelling of the topographic potential because it avoids 10-20-mGal approximation errors associated with RET techniques. The spherical approximation is suited for a range of geophysical applications, while the ellipsoidal approximation is preferable for applications requiring high accuracy or high resolution. [ABSTRACT FROM AUTHOR]

Seismic data are primarily used in studies of the Earth's lithospheric structure including the Moho geometry. In regions, where seismic data are sparse or completely absent, gravimetric or combined gravimetric-seismic methods could be applied to determine the Moho depth. In this study, we derive and present generalized expressions for solving the Vening Meinesz- Moritz's (VMM) inverse problem of isostasy for a Moho depth determination from gravity and vertical gravity-gradient data. By solving the (non-linear) Fredholm's integral equation of the first kind, the linearized observation equations, which functionally relate the (given) gravity/gravity-gradient data to the (unknown) Moho depth, are derived in the spectral domain. The VMM gravimetric results are validated by using available seismic and gravimetric Moho models. Our results show that the VMM Moho solutions obtained by solving the VMM problem for gravity and gravity-gradient data are almost the same. This finding indicates that in global applications, using the global gravity/gravity-gradient data coverage, the spherical harmonic expressions for the gravimetric forward and inverse modelling yield (theoretically) the same results. Globally, these gravimetric solutions have also a relatively good agreement with the CRUST1.0 and GEMMA GOCE models in terms of their rms Moho differences (4.7 km and 4.1 km, respectively). [ABSTRACT FROM AUTHOR]

The geoid-to-quasigeoid separation is often computed only approximately as a function of the simple planar Bouguer gravity anomaly and the height of the computation point while disregarding the contributions of terrain geometry and anomalous topographic density as well as the sub-geoid masses. In this study we demonstrate that these contributions are significant and, therefore, should be taken into consideration when investigating the relation between the normal and orthometric heights particularly in the mountainous, polar and geologically complex regions. These contributions are evaluated by applying the spectral expressions for gravimetric forward modelling and using the EIGEN-6C4 gravity model, the Earth2014 datasets of terrain, ice thickness and inland bathymetry and the CRUST1.0 sediment and (consolidated) crustal density data. Since the global crustal density models currently available (e.g. CRUST1.0) have a limited accuracy and resolution, the comparison of individual density contributions is-for consistency-realized with a limited spectral resolution up to a spherical harmonic degree 360 (or 180). The results reveal that the topographic contribution globally varies between $$-$$ 0.33 and 0.57 m, with maxima in Himalaya and Tibet. The contribution of ice considerably modifies the geoid-to-quasigeoid separation over large parts of Antarctica and Greenland, where it reaches $$\sim $$ 0.2 m. The contributions of sediments and bedrock are less pronounced, with the values typically varying only within a few centimetres. These results, however, have still possibly large uncertainties due to the lack of information on the actual sediment and bedrock density. The contribution of lakes is mostly negligible; its maxima over the Laurentian Great Lakes and the Baikal Lake reach only several millimetres. The contribution of the sub-geoid masses is significant. It is everywhere negative and reaches extreme values of $$-$$ 4.43 m. According to our estimates, the geoid-to-quasigeoid separation globally varies within $$-$$ 4.19 and 0.26 m while the corresponding values computed according to a classical definition are only negative and reach extreme values of $$-$$ 3.5 m. A comparison of these results reveals that inaccuracies caused by disregarding the terrain geometry and mass density heterogeneities distributed within the topography and below the geoid surface can reach $$\pm $$ 2 m or more in the mountainous regions. [ABSTRACT FROM AUTHOR]