Towards the Moho depth and Moho density contrast along with their uncertainties from seismic and satellite gravity observations

2017 ◽  
Vol 11 (4) ◽  
Author(s):  
M. Abrehdary ◽  
L.E. Sjöberg ◽  
M. Bagherbandi ◽  
D. Sampietro

AbstractWe present a combined method for estimating a new global Moho model named KTH15C, containing Moho depth and Moho density contrast (or shortly Moho parameters), from a combination of global models of gravity (GOCO05S), topography (DTM2006) and seismic information (CRUST1.0 and MDN07) to a resolution of 1° × 1° based on a solution of Vening Meinesz-Moritz’ inverse problem of isostasy. This paper also aims modelling of the observation standard errors propagated from the Vening Meinesz-Moritz and CRUST1.0 models in estimating the uncertainty of the final Moho model. The numerical results yield Moho depths ranging from 6.5 to 70.3 km, and the estimated Moho density contrasts ranging from 21 to 650 kg/m

2021 ◽  
Author(s):  
Lars Erik Sjöberg ◽  
Majid Abrehdary

This chapter describes a theory and application of satellite gravity and altimetry data for determining Moho constituents (i.e. Moho depth and density contrast) with support from a seismic Moho model in a least-squares adjustment. It presents and applies the Vening Meinesz-Moritz gravimetric-isostatic model in recovering the global Moho features. Internal and external uncertainty estimates are also determined. Special emphasis is devoted to presenting methods for eliminating the so-called non-isostatic effects, i.e. the gravimetric signals from the Earth both below the crust and from partly unknown density variations in the crust and effects due to delayed Glacial Isostatic Adjustment as well as for capturing Moho features not related with isostatic balance. The global means of the computed Moho depths and density contrasts are 23.8±0.05 km and 340.5 ± 0.37 kg/m3, respectively. The two Moho features vary between 7.6 and 70.3 km as well as between 21.0 and 650.0 kg/m3. Validation checks were performed for our modeled crustal depths using a recently published seismic model, yielding an RMS difference of 4 km.


2019 ◽  
Vol 13 (4) ◽  
pp. 291-303 ◽  
Author(s):  
M. Abrehdary ◽  
L. E. Sjöberg

Abstract In this research, we present a local Moho model, named MOHV19, including Moho depth and Moho density contrast (or shortly Moho constituents) with corresponding uncertainties, which are mapped from altimetric and gravimetric data (DSNSC08) in addition to seismic tomographic (CRUST1.0) and Earth topographic data (Earth2014) to a resolution of 1° × 1° based on a solution of Vening Meinesz-Moritz’ theory of isostasy. The MOHV19 model covers the area of entire European plate along with the surrounding oceans, bounded by latitudes (30 °N–82 °N) and longitudes (40 °W–70 °E). The article aims to interpret the Moho model resulted via altimetric and gravimetric information from the geological and geophysical perspectives along with investigating the relation between the Moho depth and Moho density contrast. Our numerical results show that estimated Moho depths range from 7.5 to 57.9 km with continental and oceanic averages of 41.3 ± 4.9 km and 21.6 ± 9.2 km, respectively, and an overall average of 30.9 ± 12.3 km. The estimated Moho density contrast ranges from 60.2 to 565.8 kg/m3, with averages of 421.8 ± 57.9 and 284.4 ± 62.9 kg/m3 for continental and oceanic regions, respectively, with a total average of 350.3 ± 91.5 kg/m3. In most areas, estimated uncertainties in the Moho constituents are less than 3 km and 40 kg/m3, respectively, but they reach to much more significant values under Iceland, parts of Gulf of Bothnia and along the Kvitoya Island. Comparing the Moho depths estimated by MOHV19 and those derived by CRUST1.0, MDN07, GRAD09 and MD19 models shows that MOHV19 agree fairly well with CRUST1.0 but rather poor with other models. The RMS difference between the Moho density contrasts estimated by MOHV19 and CRUST1.0 models is 49.45 kg/m3.


2014 ◽  
Vol 58 (2) ◽  
pp. 227-248 ◽  
Author(s):  
Mohammad Bagherbandi ◽  
Robert Tenzer ◽  
Lars E. Sjöberg

2020 ◽  
Author(s):  
Peter Haas ◽  
Joerg Ebbing ◽  
Wolfgang Szwillus ◽  
Philipp Tabelow

<p>We present a new inverse approach to invert satellite gravity gradients for the Moho depth under consideration of a laterally varying density contrast between crust and mantle. The inverse problem is linearized and solved with the classical Gauss-Newton algorithm in a spherical geometry. To ensure stable solutions, the Jacobian is smoothed with second-order Tikhonov regularization. During the inversion, the Moho depth is discretized into tesseroids by reference Moho depth and density contrast, from which the gravitational effect can be calculated. As a computational benefit, the Jacobian is calculated only once and afterwards weighted with the laterally varying density contrast. We look for a Moho depth model that simultaneously explains the gravity gradient field and a least misfit to existing seismic Moho depth determinations. We perform the inversion both on regional and global scale.</p><p>The laterally varying density contrast is based on different tectonic units, which are defined by independent global geological and geophysical data, such as regionalization of dispersion curves. This is beneficial in remote areas, where seismic investigations are very sparse and the crustal structure is to a large extent unknown. Applying the inversion to the Amazonian Craton and its surroundings shows a lower density contrast at the Moho depth for the continental interior compared to oceanic domains. This is in accordance with the tectono-thermal architecture of the lithosphere. The inverted values of the density vary between 300-450 kg/m<sup>3</sup>. The inverted Moho depth shows a clear separation between the Sao Francisco Craton and shallower Amazonian Craton.</p><p>Gravity inversion with a laterally varying density contrast requires a uniform reference Moho depth. On a global scale, we utilize our inversion to estimate a reference Moho depth that is in accordance with crustal buoyancy. The inverted density contrasts show a similar trend like the regional study area. The inverted Moho depth shows expected tectonic features. Our method of computing the Jacobian once and weighting with lateral variable density contrasts is a valuable optimization of standard gravity inversion.</p>


2021 ◽  
Vol 225 (3) ◽  
pp. 1952-1962
Author(s):  
M Abrehdary ◽  
L E Sjöberg

SUMMARY As recovering the crust–mantle/Moho density contrast (MDC) significantly depends on the properties of the Earth's crust and upper mantle, varying from place to place, it is an oversimplification to define a constant standard value for it. It is especially challenging in Antarctica, where almost all the bedrock is covered with a thick layer of ice, and seismic data cannot provide a sufficient spatial resolution for geological and geophysical applications. As an alternative, we determine the MDC in Antarctica and its surrounding seas with a resolution of 1° × 1° by the Vening Meinesz-Moritz gravimetric-isostatic technique using the XGM2019e Earth Gravitational Model and Earth2014 topographic/bathymetric information along with CRUST1.0 and CRUST19 seismic crustal models. The numerical results show that our model, named HVMDC20, varies from 81 kg m−3 in the Pacific Antarctic mid-oceanic ridge to 579 kg m−3 in the Gamburtsev Mountain Range in the central continent with a general average of 403 kg m−3. To assess our computations, we compare our estimates with those of some other gravimetric as well as seismic models (KTH11, GEMMA12C, KTH15C and CRUST1.0), illustrating that our estimates agree fairly well with KTH15C and CRUST1.0 but rather poor with the other models. In addition, we compare the geological signatures with HVMDC20, showing how the main geological structures contribute to the MDC. Finally, we study the remaining glacial isostatic adjustment effect on gravity to figure out how much it affects the MDC recovery, yielding a correlation of the optimum spectral window (7≤ n ≤12) between XGM2019e and W12a GIA models of the order of ∼0.6 contributing within a negligible $ \pm 14$ kg m−3 to the MDC.


2013 ◽  
Vol 43 (2) ◽  
pp. 83-98
Author(s):  
Robert Tenzer

Abstract The gravimetric inverse problem for finding the Moho density contrast is formulated in this study. The solution requires that the crust density structure and the Moho depths are a priori known, for instance, from results of seismic studies. The relation between the isostatic gravity data (i.e., the complete-crust stripped isostatic gravity disturbances) and the Moho density contrast is defined by means of the Fredholm integral equation of the first kind. The closed analytical solution of the integral equation is given. Alternative expressions for solving the inverse problem of isostasy are defined in frequency domain. The isostatic gravity data are computed utilizing methods for a spherical harmonic analysis and synthesis of the gravity field. For this purpose, we define various spherical functions, which define the crust density structures and the Moho interface globally.


2020 ◽  
Author(s):  
Majid Abrehdary ◽  
Lars Sjöberg

<p>Seismic data are the preliminary information for investigating Earth’s interior structure. Since large parts of the world are not yet sufficiently covered by such data, products from Earth satellite gravity and altimetry missions can be used as complimentary for this purpose. This is particularly the case in most of the ocean areas, where seismic data are sparse. One important information of Earth’s interior is the crustal/Moho depth, which is widely mapped from seismic surveys. In this study, we aim at presenting a new Moho depth model from satellite altimetry derived gravity and seismic data in Fennoscandia and its surroundings using the Vening Meinesz-Moritz (VMM) model based on isostatic theory. To that end, the refined Bouguer gravity disturbance (reduced for gravity of topography, density heterogeneities related to bathymetry, ice, sediments, and other crustal components by applying so-called stripping gravity corrections) is corrected for so-called non-isostatic effects (NIEs) of nuisance gravity signals from mass anomalies below the crust due to crustal thickening/thinning, thermal expansion of the mantle, Delayed Glacial Isostatic Adjustment (DGIA) and plate flexure. As Fennoscandia is a key area for GIA research, we particularly investigate the DGIA effect on the gravity disturbance and Moho depth determination from gravity in this area. To do so, the DGIA effect is removed and restored from the NIEs prior to the application of the VMM formula. The numerical results show that the RMS difference of the Moho depth from the (mostly) seismic CRUST1.0 model is 3.6/4.3 km when the above strategy for removing the DGIA effect is/is not applied, respectively. Also, the mean value differences are 0.9 and 1.5 km, respectively. Hence, our study shows that our method of correcting for the DGIA effect on gravity disturbance is significant, resulting in individual changes in Moho depth up to several kilometres.</p>


2019 ◽  
Vol 13 (1) ◽  
pp. 33-40 ◽  
Author(s):  
M. Abrehdary ◽  
L. E. Sjöberg ◽  
D. Sampietro

Abstract The determination of the oceanic Moho (or crust-mantle) density contrast derived from seismic acquisitions suffers from severe lack of data in large parts of the oceans, where have not yet been sufficiently covered by such data. In order to overcome this limitation, gravitational field models obtained by means of satellite altimetry missions can be proficiently exploited, as they provide global uniform information with a sufficient accuracy and resolution for such a task. In this article, we estimate a new Moho density contrast model named MDC2018, using the marine gravity field from satellite altimetry in combination with a seismic-based crustal model and Earth’s topographic/bathymetric data. The solution is based on the theory leading to Vening Meinesz-Moritz’s isostatic model. The study results in a high-accuracy Moho density contrast model with a resolution of 1° × 1° in oceanic areas. The numerical investigations show that the estimated density contrast ranges from 14.2 to 599.7 kg/m3 with a global average of 293 kg/m3. In order to evaluate the accuracy of the MDC2018 model, the result was compared with some published global models, revealing that our altimetric model is able to image rather reliable information in most of the oceanic areas. However, the differences between this model and the published results are most notable along the coastal and polar zones, which are most likely due to that the quality and coverage of the satellite altimetry data are worsened in these regions.


2021 ◽  
Author(s):  
Willi Freeden

<p>The lecture highlights arguments that, coming from multiscale mathematics, have fostered the advancement of gravimetry, as well as those that, generated by gravimetric problems, have contributed to the enhancement in constructive approximation and numerics. Inverse problems in gravimetry are delt with multiscale mollifier decorrelation strategies. Two examples are studied in more detail: (i) Vening Meinesz multiscale surface mollifier regularization to determine locally the Earth's disturbing potential from deflections of vertical, (ii) Newton multiscale volume mollifier regularization of the inverse gravimetry problem to derive locally the density contrast distribution from functionals of the Newton integral and to detect fine particulars of geological relevance. All in all, the Vening Meinesz medal  lecture is meant as an  \lq \lq appetizer'' served to enjoy the tasty meal "Mathematical Geoscience Today'' to be shared by geoscientists and mathematicians in the field of gravimetry. It provides innovative concepts and locally relevant applications presented in a monograph to be published by Birkhäuser in the book series “Geosystems Mathematics” (2021).</p>


Sign in / Sign up

Export Citation Format

Share Document