Search

A levelling network was readjusted and a new geoid model compiled within the framework of geodetic vertical datum modernization at the Hong Kong territories. To accomplish all project objectives, the quasigeoid model has to be determinedtoo. A quasigeoid model can be obtained from existing geoid model by applying the geoid-to-quasigeoid separation. Thegeoid-to-quasigeoid separation was traditionally computed as a function of the simple planar Bouguer gravity anomaly, whiledisregarding terrain geometry, topographic density variations, and vertical gravity changes due to mass density heterogeneities below the geoid surface. We applied this approximate method because orthometric heights of levelling benchmarksin Hong Kong were determined only approximately according to Helmert’s theory of orthometric heights. Considering afurther improvement of the accuracy of orthometric heights by applying advanced numerical procedures, we determinedthe geoid-to-quasigeoid separation by applying an accurate method. The comparison of the accurately and approximatelycomputed values of the geoid-to-quasigeoid separation revealed signifcant diferences between them. The approximatevalues are all negative and reach -2.8 cm, whereas values from the accurate method vary between -4.1 and+0.2 cm. In addition, we assessed the efect of anomalous topographic density on the geoid-to-quasigeoid separation by employing a newlydeveloped digital rock density model. According to our estimates the efect of anomalous topographic density reaches amaximum value of 1.6 cm, refecting a predominant presence of light volcanic rocks and sedimentary deposits at the HongKong territories. Our numerical fndings indicate that the conversion between geoid and quasigeoid models should be doneaccurately, even in regions with a moderately elevated topography

The high-resolution Moho depth model is required in various geophysical studies. However, the available models’ resolutions could be improved for this purpose. Large parts of the world still need to be sufficiently covered by seismic data, but existing global Moho models do not fit the present-day requirements for accuracy and resolution. The isostatic models can relatively reproduce a Moho geometry in regions where the crustal structure is in an isostatic equilibrium, but large segments of the tectonic plates are not isostatically compensated, especially along active convergent and divergent tectonic margins. Isostatic models require a relatively good knowledge of the crustal density to correct observed gravity data. To overcome the lack of seismic data and non-uniqueness of gravity inversion, seismic and gravity data should be combined to estimate Moho geometry more accurately. In this study, we investigate the performance of two techniques for combining long- and short-wavelength Moho geometry from seismic and gravity data. Our results demonstrate that both Butterworth and spectral combination techniques can be used to model the Moho geometry. The results show the RMS of Moho depth differences between our model and the reference models are between 1.7 and 4.7 km for the Butterworth filter and between 0.4 and 4.1 km for the spectral combination.

Since marine seismic studies are relatively sparse and unevenly distributed, detailed tomographic images of the Moho geometry under large parts of the world’s oceans and marginal seas are not yet available. Marine gravity data is, therefore, often used to detect the Moho depth in these regions. Alternatively, Airy’s isostatic theory can be applied for this purpose. In this study, we compare different isostatic and gravimetric methods for a Moho recovery under the oceanic crust and continental margins, particularly focusing on a numerical performance of Airy, Vening Meinesz-Moritz (VMM), direct gravity inversion, and generalized (for the Earth’s spherical approximation) Parker-Oldenburg methods. Numerical experiments are conducted to estimate the Moho depth beneath the Indian Ocean. Results reveal that, among these investigated methods, the VMM model is probably the most suitable for a gravimetric Moho recovery beneath the oceanic crust and continental margins, when taking into consideration the lithospheric mantle density information. This method could to some extent model realistically a Moho geometry beneath mid-oceanic spreading ridges, oceanic subductions, most of oceanic volcanic formations, and marine sediment deposits. Nonetheless, this model still cannot fully reproduce a gradual Moho deepening caused by a conductive cooling and a subsequent isostatic rebalance of the oceanic lithosphere, which can functionally be described by a Moho deepening with the increasing ocean-floor age. Results also indicate that the Airy method typically overestimates the Moho depth under oceanic volcanic formations, while the direct gravity inversion and generalized Parker-Oldenburg methods could not reproduce more detailed features in the Moho geometry. Since Pratt’s theory better describes a large-scale isostatic mechanism of the oceanic lithosphere by means of compensation density variations, but does not account for additional changes in compensation depth (i.e., Moho depth)

Temporal variations in the Earth’s gravity field can be used for monitoring of lithospheric deformations. The network of continuously operating gravity stations is required for this purpose but a global coverage by such network is currently extremely sparse. Temporal variations in long-wavelength part of the Earth’s gravity field have been, however, observed by two satellite missions, namely the Gravity Recovery And Climate Experiment (GRACE) and the GRACE Follow-On (GRACE-FO). These satellite gravity observations can be used to study long-wavelength deformations of the lithosphere. Consequently, considering the lithosphere as a spherical elastic shell and solving the partial differential equation of elasticity for it, the stress, strain and displacement inside the lithosphere can be estimated. The lower boundary of this shell is assumed to be stressed by mantle convection, which has a direct relation to the Earth’s gravity field according to Runcorn’s theory. Changes in gravity field lead to changes in the sublithospheric stress and the stress propagated throughout the lithosphere. In this study, we develop mathematical models in spherical coordinates for describing the stress propagation from the sublithosphere through the lithosphere. We then organize a system of observation equations for finding a special solution to the boundary-value problem of elasticity in the way that provides a stable solution. In contrast, models presented in previously published studies are ill-posed. Furthermore, we use constants of the solution determined from the boundary stresses to determine the strain and displacements leading to these stresses, while in previous studies only the stress has been considered according to rheological properties of the lithosphere. We demonstrate a practical applicability of this theoretical model to estimate the stress–strain redistribution caused by the Sar-e-Pol Zahab 2018 earthquake in Iran by using the GRACE-FO monthly solutions.

In geophysical and geodetic studies, gravity inversion is typically performed such that observed gravity values are first continued downward onto a regular (planar, spherical or spheroidal) surface by solving an inverse integral transform, which originates from a classical solution to the first boundary-value problem in potential theory. A typical example is continuing gravity observed at the topographic surface down to the mean sea level (geoid). Nowadays, gravity-dedicated satellite missions and aerial gravimetry provide gravity data above the topographic surface in addition to classical terrestrial gravity observations. For specific purposes (such as gravity data combination and validation, or quasigeoid determination), satellite and aerial gravity observations have to be continued to the irregular topographic surface. In this study, we address this issue by formulating a functional model for a spectral downward continuation of selected gravitational field quantities to an irregular topographic surface. Moreover, we generalize this functional model to allow for transformation between different types of gravitational field quantities. In particular, we derive spectral weights for estimation of the disturbing potential or disturbing/anomalous gravity at the Earth’s surface by combining the first-, second- and third-order radial gradients of the disturbing potential (disturbing gradients). The correctness of the developed combined spectral estimator is verified in a closed-loop test based on synthetic satellite disturbing gradients. The combined spectral estimator is applied to simulated satellite disturbing gradients polluted by a realistic Gaussian noise. Results of the numerical experiments show that the combined spectral estimator puts the highest importance on the least polluted disturbing gradient, while the contribution of the least accurate disturbing gradient is negligible. An important advantage of this spectral combination method is that no matrix inversi

We estimate the (effective) elastic thickness of the Iranian lithosphere (and adjoining tectonic plates) by using the approach that combines the Vening Meinesz-Moritz's (VMM) regional isostatic principle with the isostatic flexural model formulated based on solving a flexural differential equation for a thin elastic shell. To model the response on a load more realistically, we also consider the lithospheric density structure. The resulting expression describes a functional relation that links gravity field and mechanical properties of the lithosphere. The Young modulus and the Poisson ratio are computed from seismic velocity data in prior of estimating the lithospheric elastic thickness. The presented results reveal that the estimated elastic thickness closely resembles a regional tectonic configuration associated with the extensional tectonism along the Red Sea-Gulf Rift System, the continental collision of the Arabian and Eurasian plates, and the subduction along the Makran Subduction Zone. Seismically and volcanically active convergent tectonic margins of the Zagros and Kopeh Dagh Fold and Thrust Belts further extending along the Makran Accretionary Complex are characterised by a low lithospheric strength, with the elastic thickness typically less than ∼30 km. These small values of the elastic thickness are in a striking contrast to much larger values within most of the Central Iranian Blocks. According to our estimate, local maxima there reach ∼70 km in the Tabas micro-block. The elastic thickness of the Turan and Arabian Platforms reaches maxima of ∼100 km. These results generally support the hypothesis that tectonically active zones and orogens have a relatively low strength, resulting in a significant response of the lithosphere on various tectonic loads, compared to a significant strength of old cratonic formations. Interestingly, however, we observe a striking contrast between a low strength of the Arabian Shield compared to a high strength of the Arabian P