Moho Depth Model Research Articles

Moho depth model of South America from a machine learning approach

Crustal structure models play an important role in the characterization of seismogenic zones, in the regional delimitation of geological provinces and particularly, in understanding the genesis and evolution of sedimentary basins. Despite the increasing number of new seismic surveys and seismographic networks in South America, crustal thickness measurements are still scarce and irregularly sampled, reducing the resolution of crustal model maps. To overcome these challenges, a novel approach based on machine learning techniques is proposed, in order to explore higher resolution gravity datasets in the interpolation of crustal thickness measurement points obtained from previous seismic/seismological compilations. The algorithm used in this study is based on Gaussian processes prediction methods, which allowed inferring the depth of Moho to South America. The prediction error of the model obtained from the training and testing database was 3.48 km, which is compatible with the uncertainties derived from the H-k stacking analysis. The depth range varied from 69.8 km beneath the Andes to 4.3 km in oceanic regions. The average Moho depth for the South American Platform is 39.1 km, allowing a spatial correlation of deeper and shallower Moho regions with different types of continental basins. Compared to other models, the model resulting from this study presents fine-scale features highlighting the limits of the main tectonic domains and a good agreement with the suture zones. Overall, this study demonstrates the potential of applying machine learning tools in crustal-scale imaging using sparse datasets, providing new advances in Moho modeling of the South America, as well as new perspectives on its the history and tectonic evolution.

A three-dimensional Moho depth model beneath the Yemeni highlands and rifted volcanic margins of the Red Sea and Gulf of Aden, Southwest Arabia

Knowing Moho discontinuity undulation is fundamental to understanding mechanisms of lithosphere-asthenosphere interaction, extensional tectonism and crustal deformation in volcanic passive margins such as the study area, which is located in the southwestern corner of the Arabian Peninsula bounded by the Red Sea and the Gulf of Aden. In this work, a 3D Moho depth model of the study area is constructed for the first time by inverting gravity data from the Earth Gravitational Model (EGM2008) using the Parker-Oldenburg algorithm. This model indicates the shallow zone is situated at depths of 20 km to 24 km beneath coastal plains, whereas the deep zone is located below the plateau at depths of 30 km to 35 km and its deepest part coincides mainly with the Dhamar-Rada’a Quaternary volcanic field. The results also indicate two channels of hot magmatic materials joining both the Sana’a-Amran Quaternary volcanic field and the Late Miocene Jabal An Nar volcanic area with the Dhamar-Rada’a volcanic field. This conclusion is supported by the widespread geothermal activity (of mantle origin) distributed along these channels, isotopic data, and the upper mantle low velocity zones indicated by earlier studies.©2023 China Geology Editorial Office.

Moho Inversion by Gravity Anomalies in the South China Sea: Updates and Improved Iteration of the Parker‐Oldenburg Algorithm

AbstractThe Moho is the interface between crust and mantle, and accurate location of the Moho is important for both resource exploration and deep earth condition and structural change investigations. The theory of the traditional Parker‐Oldenburg (P‐O) method is quite simple and it is widely applied in the frequency domain of Moho depth inversion. However, Moho fluctuation simulations using the P‐O method are not reliable because of the lack of field geographic data constraints during the inversion process and the excessively smoothing of data details caused by using filters to correct the source data signals. To solve those problems, we propose an improved iteration P‐O method with a variable density contrast model, the iterative process is constrained by geological data in the inversion parameters, and the variable depth of the gravity interface is iterated using an equivalent form of upward continuation in the Fourier domain, which is more stable and convergent than downward continuation term in original P‐O method. Synthetic experiments indicate that improved method has the better consistency among the simulations than original method, and our improved method has the smallest root mean square (RMS) of 0.59 km. In a real case, we employed the improved method to invert for the Moho depth of the South China Sea (SCS), and the RMS between our Moho depth model and the seismological data is the smallest value of 2.27 km. The synthetic experiments and application of the model to the SCS further prove that our method is practical and efficient.

Papua New Guinea Moho inversion based on XGM 2019e gravity field model

The construction of the high-resolution Moho depth model is significant for studying the characteristics of the complex tectonic movement (seafloor spreading, plate subduction phenomena) in Papua New Guinea. We calculate the region’s Moho relief and lithosphere thinning factor using the XGM 2019e gravity field model and nonlinear fast gravity inversion method under the GEMMA Moho depth model’s constraint considering the influence of lithosphere thermal gravity anomaly. The calculation result shows that the Moho depth is between 6—34 km, forming two large depressions in Woodlark Basin (WB) and Solomon Sea Plate (SSP) with deep scattered islands. In addition, the findings suggest that Significant differences exist in the shape and tectonic movement intensity of the North and South oceanic crust at the WB. Nevertheless, the lithosphere extends evenly in Manus Basin (MB). WB collided with the Solomon Islands at a higher angle than the SSP subducted under Bismarck Sea Plate (BSP); strong earthquakes may frequently occur on both sides and in deeper positions at West New Britain Trench in the future.

Three-dimensional Moho depth model of the eastern Indian shield and its isostatic implications

Three-dimensional Moho depth model of the eastern Indian shield and its isostatic implications

A High-Resolution Global Moho Model from Combining Gravimetric and Seismic Data by Using Spectral Combination Methods

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.

A refined Moho depth model from a joint analysis of gravity and seismic data of the South China Sea basin and its tectonic implications

As a typical marginal sea in the western Pacific, the South China Sea (SCS) has experienced a complex evolutionary process, and mapping its Moho depth can help us understand the accretion process of the oceanic crust. In this study, we present a refined Moho depth model in the SCS based on gravity inversion constrained by seismic data. For a more reliable gravity-derived Moho depth, an improved sediment thickness grid is constructed from a large volume of seismic reflection/refraction data. The gravity-derived Moho depths are in good agreement with the seismically determined ones, with the majority of differences falling within ±2.5 km. Our results reveal a significant crustal thickness asymmetry across the mid-ocean ridge of the East Sub-basin of the SCS basin since the ridge jump at 23.6 Ma. We attribute this asymmetry to higher residual mantle temperature and more magma supply north of the spreading center. The higher mantle temperature to the north is evidenced by the shallower Curie depth estimated from magnetic anomalies than to the south of the spreading ridge.

MOHV21: a least squares combination of five global Moho depth models

The purpose of this study is to determine MOHV21, a Moho depth model based on an optimal combination of five global seismic and gravimetric-isostatic models of Moho depth by a weighted least squares approach at a resolution of 1° × 1°. For proper weighting among the data, the study starts with determining (mostly missing) standard errors and correlations among the models. The standard errors among the input models range from 1.0 (in Brazil) to 6.8 km (in Peru) and from 0.1 (in Huna Bay) to 6.0 km (in East Pacific Ridge) for Moho depth on land and ocean, respectively. The correlations among the five models range between − 0.99 and + 0.90. The Moho depths for MOHV21 at land regions vary between 14.5 (at the Horn of Africa) and 75 km (in the Himalayas) and between 6.6 (in the Greenland Sea) and 51.8 (in the Gulf of Bothnia) for land and ocean regions, respectively (However, note that, the Gulf of Bothnia belongs to continental crust, while the oceanic crust is generally within 20 km). The standard errors are generally within a few km but reaches 6.8 km (9%) in the highest mountains. The shallow Moho depths along mid-ocean ridges are well exposed in the model. Notable regional Moho highs are visualized in the Tarim basin in NW China of 59 ± 6.5 km and in Central Finland of 57 ± 4.7 km. A comparison of MOHV21 with a mosaic of regional models shows large differences reaching ± 25 km in Africa, Antarctic, and parts of S. America, while the differences are relatively modest in those parts of oceans that are available in the regional models.

Lithospheric Structure Near Jiuyishan, South China: Implications for Asthenospheric Upwelling and Lithospheric Modification

AbstractIn this study, a high‐resolution 3D lithospheric S wave velocity model of the Jiuyishan region is constructed by joint inversion of receiver function and ambient noise. Our Moho depth model reveals relatively thick crust (∼35 km) beneath the Yangtze block and thin crust (∼27 km) beneath the Cathaysia block. The model depicts a complex boundary between these blocks, roughly following the Chenzhou‐Linwu fault in the northern part but gradually curving westward in the southern part between the two blocks. The 3D S wave velocity model shows low‐velocity anomalies in the middle crust, as well as in the upper mantle, interpreted as evidence of asthenospheric upwellings in the upper mantle. These upwellings modify lithosphere structures by changing their temperatures and compositions, and may provide a heat source through the Chenzhou‐Linwu fault for the low‐velocity anomaly in the shallow part beneath Jiuyishan region.

Two-step Gravity Inversion Reveals Variable Architecture of African Cratons

The lithospheric build-up of the African continent is still to a large extent unexplored. In this contribution, we present a new Moho depth model to discuss the architecture of the three main African cratonic units, which are: West African Craton, Congo Craton, and Kalahari Craton. Our model is based on a two-step gravity inversion approach that allows variable density contrasts across the Moho depth. In the first step, the density contrasts are varied for all non-cratonic units, in the second step for the three cratons individually. The lateral extension of the tectonic units is defined by a regionalization map, which is calculated from a recent continental seismic tomography model. Our Moho depth is independently constrained by pointwise active seismics and receiver functions. Treating the constraints separately reveals a variable range of density contrasts and different trends in the estimated Moho depth for the three cratons. Some of the estimated density contrasts vary substantially, caused by sparse data coverage of the seismic constraints. With a density contrast of Δρ= 200 kg/m3the Congo Craton features a cool and undisturbed lithosphere with smooth density contrasts across the Moho. The estimated Moho depth shows a bimodal pattern with average Moho depth of 39–40 km for the Kalahari and Congo Cratons and 33–34 km for the West African Craton. We link our estimated Moho depth with the cratonic extensions, imaged by seismic tomography, and with topographic patterns. The results indicate that cratonic lithosphere is not necessarily accompanied by thick crust. For the West African Craton, the estimated thin crust, i.e. shallow Moho, contrasts to thick lithosphere. This discrepancy remains enigmatic and requires further studies.

A local lithospheric structure model for Vietnam derived from a high-resolution gravimetric geoid

High-resolution Moho and lithosphere–asthenosphere boundary depth models for Vietnam and its surrounding areas are determined based on a recently released geoid model constructed from surface and satellite gravity data (GEOID_LSC_C model) and on 3ʹʹ resolution topography data (mixed SRTM model). A linear density gradient for the crust and a temperature-dependent density for the lithospheric mantle were used to determine the lithospheric structure under the assumption of local isostasy. In a first step, the impact of correcting elevation data from sedimentary basins to estimate Moho depth has been evaluated using CRUST1.0 model. Results obtained from a test area where seismic data are available, which demonstrated that the sedimentary effect should be considered before the inversion process. The geoid height and elevation-corrected sedimentary layer were filtered to remove signals originating below the lithosphere. The resulting Moho and lithosphere–asthenosphere boundary depth models computed at 1ʹ resolution were evaluated against seismic data as well as global and local lithospheric models available in the study region. These comparisons indicate a consistency of our Moho depth estimation with the seismic data within 1.5 km in standard deviation for the whole Vietnam. This new Moho depth model for the study region represents a significant improvement over the global models CRUST1.0 and GEMMA, which have standard deviations of 3.2 and 3.3 km, respectively, when compared to the seismic data. Even if a detailed geological interpretation of the results is out of scope of this paper, a joint analysis of the obtained models with the high-resolution Bouguer gravity anomaly is finally discussed in terms of the main geological patterns of the study region. The high resolution of our Moho and lithosphere–asthenosphere boundary depth models contribute to better constrain the lithospheric structure as well as tectonic and geodynamic processes of this region. The differences in Moho depth visible in the northeast and southwest sides of the Red River Fault Zone confirmed that the Red River Fault Zone may be considered the boundary between two continental blocks: South China and Indochina blocks. However, no remarkable differences in lithosphere–asthenosphere boundary depth were obtained from our results. This suggests that the Red River Fault Zone developed within the crust and remained a crustal fault.

A New Moho Depth Model for Fennoscandia with Special Correction for the Glacial Isostatic Effect

In this study, we present a new Moho depth model in Fennoscandia and its surroundings. The model is tailored from data sets of XGM2019e gravitationl field, Earth2014 topography and seismic crustal model CRUST1.0 using the Vening Meinesz-Moritz model based on isostatic theory to a resolution of 1° × 1°. To that end, the refined Bouguer gravity disturbance is determined by reducing the observed field for gravity effect of topography, density heterogeneities related to bathymetry, ice, sediments, and other crustal components. Moreover, stripping of non-isostatic effects of gravity signals from mass anomalies below the crust due to crustal thickening/thinning, thermal expansion of the mantle, Delayed Glacial Isostatic Adjustment (DGIA), i.e., the effect of future GIA, and plate flexure has also been performed. As Fennoscandia is a key area for GIA research, we particularly investigate the DGIA effect on the gravity disturbance and gravimetric Moho depth determination in this area. One may ask whether the DGIA effect is sufficiently well removed in the application of the general non-isostatic effects in such an area, and to answer this question, the Moho depth is determined both with and without specific removal of the DGIA effect prior to non-isostatic effect and Moho depth determinations. The numerical results yield that the RMS difference of the Moho depth from our model HVMD19 vs. the seismic CRUST19 and GRAD09 models are 3.8/4.2 km and 3.7/4.0 km when the above strategy for removing the DGIA effect is/is not applied, respectively, and the mean value differences are 1.2/1.4 km and 0.98/1.4 km, respectively. Hence, our study shows that the specific correction for the DGIA effect on gravity disturbance is slightly significant, resulting in individual changes in the gravimetric Moho depth up to − 1.3 km towards the seismic results. On the other hand, our study shows large discrepancies between gravimetric and seismic Moho models along the Norwegian coastline, which might be due to uncompensated non-isostatic effects caused by tectonic motions.

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Victoria Land, Antarctica: An Improved Geodynamic Interpretation Based on the Strain Rate Field of the Current Crustal Motion and Moho Depth Model

In Antarctica, the severe climatic conditions and the thick ice sheet that covers the largest and most internal part of the continent make it particularly difficult to systematically carry out geophysical and geodetic observations on a continental scale. It prevents the comprehensive understanding of both the onshore and offshore geology as well as the relationship between the inner part of East Antarctica (EA) and the coastal sector of Victoria Land (VL). With the aim to reduce this gap, in this paper multiple geophysical dataset collected since the 1980s in Antarctica by Programma Nazionale di Ricerche in Antartide (PNRA) were integrated with geodetic observations. In particular, the analyzed data includes: (i) Geodetic time series from Trans Antarctic Mountains DEFormation (TAMDEF), and Victoria Land Network for DEFormation control (VLNDEF) GNSS stations installed in Victoria Land; (ii) the integration of on-shore (ground points data and airborne) gravity measurements in Victoria Land and marine gravity surveys performed in the Ross Sea and the narrow strip of Southern Ocean facing the coasts of northern Victoria Land. Gravity data modelling has improved the knowledge of the Moho depth of VL and surrounding the offshore areas. By the integration of geodetic and gravitational (or gravity) potential results it was possible to better constrain/identify four geodynamic blocks characterized by homogeneous geophysical signature: the Southern Ocean to the N, the Ross Sea to the E, the Wilkes Basin to the W, and VL in between. The last block is characterized by a small but significant clockwise rotation relative to East Antarctica. The presence of a N-S to NNW-SSE 1-km step in the Moho in correspondence of the Rennick Geodynamic Belt confirms the existence of this crustal scale discontinuity, possibly representing the tectonic boundary between East Antarctica and the northern part of VL block, as previously proposed by some geological studies.

Sensitivity analysis of gravity gradient inversion of the Moho depth—a case example for the Amazonian Craton

SUMMARYWe present a novel approach for linearized gravity inversion to estimate the Moho depth, which allows the use of any gravitational component instead of the vertical gravity component only. The inverse problem is solved with the Gauss–Newton algorithm and the gravitational field of the undulating Moho depth is calculated with tesseroids. Hereby, the density contrast can be laterally variable by using information from seismological regionalization. Our approach is illustrated with a synthetic example, which we use to explore different regularization parameters. The vertical gravity gradient gzz provides the most reasonable results with appropriate parameters. As a case example, we invert for the Moho depth of the Amazonian Craton and its surroundings. The results are constrained by estimates from active seismic measurements. Our new Moho depth model correlates to tectonic domains and is in agreement with previous models. The estimated density contrasts of the tectonic domains agree well with the lithospheric architecture and show with 300–450 kg m–3 lower density contrasts for continental domains, whereas the oceans reveal a density contrast of 450–500 kg m–3. The wider range of estimated density contrast for the continent reflects uncertainties in Precambrian Fold Belts that arise from its small gravity signal. Our results demonstrate that a variable density contrast at the Moho depth is a valuable enhancement for gravity inversion.

Moho depth of the British Isles: a probabilistic perspective

SUMMARYWe present a new Moho depth model of the British Isles and surrounding areas from the most up-to-date compilation of Moho depth estimates obtained from refraction, reflection and receiver function data. We use a probabilistic, trans-dimensional and hierarchical approach for the surface reconstruction of Moho topography. This fully data-driven approach allows for adaptive parametrization, assessment of relative importance between different data-types and uncertainties quantification on the reconstructed surface. Our results confirm the first order features of the Moho topography obtained in previous work such as deeper Moho (29–36 km) in continental areas (e.g. Ireland and Great Britain) and shallower Moho (12–22 km) offshore (e.g. in the Atlantic Ocean, west of Ireland). Resolution is improved by including recent available data, especially around the Porcupine Basin, onshore Ireland and Great Britain. NE trending features in Moho topography are highlighted above the Rockall High (about 28 km) and the Rockall Trough (with a NE directed deepening from 12 to about 20 km). A perpendicular SE oriented feature (Moho depth 26–28 km) is located between the Orkney and the Shetland, extending further SW in the North Sea. Onshore, our results highlight the crustal thinning towards the N in Ireland and an E–W oriented transition between deep (34 km) and shallow (about 28 km) Moho in Scotland. Our probabilistic results are compared with previous models showing overall differences around ±2 km, within the posterior uncertainties calculated with our approach. Bigger differences are located where different data are used between models or in less constrained areas where posterior uncertainties are high.

Moho Depth Variations From Receiver Function Imaging in the Northeastern North China Craton and Its Tectonic Implications

AbstractA detailed knowledge of the crustal thickness in the northeastern North China Craton (NCC) is important for understanding the unusual Phanerozoic destruction of the craton. We achieve this goal by employing a 2‐D wave equation‐based migration method to P receiver functions from 198 broadband seismic stations, using Ps conversions and surface‐reflected multiples. By combining receiver function images along 19 profiles, we constructed a high‐resolution Moho depth model for the northeastern NCC. The results present dominant E‐W Moho depth variations similar to previous observations and new regional N‐S variations beneath both sides of the North‐South Gravity Lineament. To the west, while a deeper Moho (∼42 km) appears in the interior of the Trans‐North China Orogen, a relatively shallow Moho (∼38 km) exists in the northern margin of the Trans‐North China Orogen to western NCC. To the east, the crust beneath the Yan Mountains in the marginal area is thicker (∼32–40 km) than that (∼26–32 km) beneath the Bohai Bay Basin in the craton interior, and the Moho further shallows from NE (∼32 km) to SW (∼26 km) within the basin. Along with other observations, we conclude that the dominant E‐W difference may have been associated with the Paleo‐Pacific plate subduction under eastern Asia since the Mesozoic. The newly observed complex N‐S variations may have reflected the structural heterogeneity of the cratonic lithosphere inherited since the formation of the NCC in the Paleoproterozoic, or spatially uneven effects on the cratonic lithosphere of subsequent thermotectonic events during the long‐term evolution of the craton, or both.

The Crustal Structure of the North–South Earthquake Belt in China Revealed from Deep Seismic Soundings and Gravity Data

The North–South earthquake belt (NSEB) is one of the major earthquake regions in China. The studies of crustal structure play a great role in understanding tectonic evolution and in evaluating earthquake hazards in this region. However, some fundamental crustal parameters, especially crustal interface structure, are not clear in this region. In this paper, we reconstructed the crustal interface structure around the NSEB based on both the deep seismic sounding (DSS) data and the gravity data. We firstly reconstructed the crustal structure of crystalline basement (interface G), interface between upper and lower crusts (interface C) and Moho in the study area by compiling the results of 38 DSS profiles published previously. Then, we forwardly calculated the gravity anomalies caused by the interfaces G and C, and then subtracted them from the complete Bouguer gravity anomalies, yielding the regional gravity anomalies mainly due to the Moho interface. We then utilized a lateral-variable density interface inversion technique with constraints of the DSS data to invert the regional anomalies for the Moho depth model in the study area. The reliability of our Moho depth model was evaluated by comparing with other Moho depth models derived from other gravity inversion technique and receiver function analysis. Based on our Moho depth model, we mapped the crustal apparent density distribution in the study area for better understanding the geodynamics around the NSEB.

Moho depth model for the Central Asian Orogenic Belt from satellite gravity gradients

AbstractThe main purpose of this study is to construct a new 3‐D model of the Central Asian Orogenic Belt (CAOB) crust, which can be used as a starting point for future lithospheric studies. The CAOB is a Paleozoic accretionary orogen surrounded by the Siberian Craton to the north and the North China and Tarim Cratons to the south. This area is of great interest due to its enigmatic and still not completely understood geodynamic evolution. First, we estimate an initial crustal thickness by inversion of the vertical gravity component of the Gravity Field and Steady‐State Ocean Circulation Explorer (GOCE) and DTU10 models. Second, 3‐D forward modeling of the GOCE gravity gradients is performed, which determines the topography of the Moho, the geometry, and the density distribution of the deeper parts of the CAOB and its surroundings, taking into account the lateral and vertical density variations of the crust. The model is constrained by seismic refraction, reflection, and receiver function studies and geological studies. In addition, we discuss the isostatic implications of the differences between the seismic Moho and the resulting 3‐D gravity Moho, complemented by the analysis of the lithostatic load distribution at the upper mantle level. Finally, the correlation between the contrasting tectonic domains and the thickness of the crust reveals the inheritance of Paleozoic and Mesozoic geodynamics, particularly the magmatic provinces and the orocline which preserve their crustal features.

High-resolution Moho model for Greenland from EIGEN-6C4 gravity data

The crust–mantle boundary (the Moho) is a first order interface in the Earth and the depth to the Moho is therefore well studied in most regions. However, below regions which are covered by large ice sheets, such as Greenland and Antarctica, the Moho is only partly known and seismic data are difficult to obtain. Here, we take advantage of the global gravity model EIGEN-6C4, together with the Parker-Oldenburg algorithm, to estimate the depth to the Moho beneath Greenland and surroundings. The available free-air gravity data are corrected for the topographic effect and the effect of sedimentary basins. We also correct for the effect on gravity due to the weight of the ice sheet and the accompanying deflection of the Earth's surface, which has not previously been taken into account in gravity studies of currently glaciated regions. Our final Moho depth model for Greenland has an associated uncertainty of ±4.5km for areas with sedimentary basins and ±4km for areas without sedimentary basins. The model shows maximum Moho depths below east Greenland of up to 55km and values less than 20km offshore east Greenland. There is a marked increase in Moho depth of 10–15km from northern to central Greenland, indicating a significant change in geology. A deep Moho at the northern coast of Greenland towards Ellesmere Island might be related to the location of the hot-spot track. Our Moho model is consistent with previous models, but has a higher lateral resolution of 0.1° and covers the entire area of on- and offshore Greenland.