Geothermal heat flux under the Antarctic ice is one of the least known parameters. Different methods (based on e.g., magnetic or seismic data) have been applied in recent years to quantify the ...thermal structure and the geothermal heat flux, resulting in vastly different estimates. In this study, we use a Bayesian Monte-Carlo-Markov-Chain approach to explore the consistency of such models and to which degree lateral variations of the thermal parameters are required. Hereby, we evaluate the input from different lithospheric models and how they influence surface heat flux. We demonstrate that both Curie isotherm and heat production are dominating parameters for the thermal calculation and that use of incorrect models or sparsely available data lead to unreliable results. As an alternative approach, geological information should be coupled with geophysical data analysis, as we demonstrate for the Antarctic Peninsula.
SUMMARY
We present a new global model for the Earth’s lithosphere and upper mantle (LithoRef18) obtained through a formal joint inversion of 3-D gravity anomalies, geoid height, satellite-derived ...gravity gradients and absolute elevation complemented with seismic, thermal and petrological prior information. The model includes crustal thickness, average crustal density, lithospheric thickness, depth-dependent density of the lithospheric mantle, lithospheric geotherms, and average density of the sublithospheric mantle down to 410 km depth with a surface discretization of 2° × 2°. Our results for lithospheric thickness and sublithospheric density structure are in excellent agreement with estimates from recent seismic tomography models. A comparison with higher resolution regional studies in a number of regions around the world indicates that our values of crustal thickness and density are an improvement over a number of previous global crustal models. Given the strong similarity with recent tomography models down to 410 km depth, LithoRef18 can be readily merged with these seismic models to include seismic velocities as part of the reference model. We include several analyses of robustness and reliability of input data, method and results. We also provide easy-to-use codes to interrogate the model and use its predictions for the development of higher-resolution models.
Considering the model‘s features and data fitting statistics, LithoRef18 will be useful in a wide range of geophysical and geochemical applications by serving as a reference or initial lithospheric model for (i) higher-resolution gravity, seismological and/or integrated geophysical studies of the lithosphere and upper mantle, (ii) including far-field effects in gravity-based regional studies, (iii) global circulation/convection models that link the lithosphere with the deep Earth, (iv) estimating residual, static and dynamic topography, (v) thermal modelling of sedimentary basins and (vi) studying the links between the lithosphere and the deep Earth, among others. Several avenues for improving the reliability of LithoRef18’s predictions are also discussed. Finally, the inversion methodology presented in this work can be applied in other planets for which potential field data sets are either the only or major constraints to their internal structures (e.g. Moon, Venus, etc.).
Display omitted
•New method of global gravity gradient inversion for the Moho depth has been developed.•Estimated Moho depth is highlighted in a Gondwana reconstruction.•Secular trend shows ...increasing crustal thickness towards end of Archean.
We present a global study of the crustal structure with emphasis on the twelve main cratons of the earth. In an inverse scheme, satellite gravity gradient data from the GOCE mission are inverted for the Moho depth, exploiting laterally variable density contrasts based on seismic tomography. Our results are constrained by an active source seismic data base, as well as a tectonic regionalization map, derived from seismic tomography. For the global analysis, we implement a moving window approach to perform the gravity inversion, followed by interpolating the estimated density contrasts of common tectonic units with a flood-fill algorithm. The Moho depth model reveals variable patterns, with average values ranging between 33–40 km for the individual cratons. We observe low density contrasts for the cratons in the northern hemisphere, indicating old cratonic lithosphere with only smooth density gradients, and variable density contrasts for the other cratons. The estimated Moho depth is implemented in a Gondwana reconstruction, highlighting the characteristics of once connected cratons with regard to their tectonic evolution. Further, we investigate the estimated Moho depth together with the stabilization age of the individual cratons. We identify a secular trend of increasing crustal thickness proceeding from the Archean that reaches its turning point at the Archean-Proterozoic boundary. The Paleoproterozoic cratons are characterized by shallow average Moho depth, reflecting exhumation of the lower crust during orogenic events. We propose that magmatic underplating and isostatic adjustment of elevated topography after craton stabilization cause deviations from the observed secular trend of Moho depth variation. This is partly reflected by an anomalous deep Moho depth. We emphasize that our estimated Moho depth and density contrasts are suitable for a wide range of applications, not only for solid earth community, but also for interdisciplinary global scale studies, as well as local studies in the oceanic domain.
SUMMARY
We 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.
SUMMARY
Estimating the depth to magnetic bottom (DTB) from magnetic data is one of the most important and difficult potential field inversion problems. Since DTB can often be linked to the Curie ...isotherm depth of magnetite (∼580 °C), it could provide crucial constraints on heat flow, even in remote or inaccessible areas. Spectral methods are the most popular approach to estimate DTB, but their reliability has been challenged on many grounds. In contrast, space-domain methods have received relatively little attention, even though they might avoid some of the limitations of spectral methods. Furthermore, many DTB estimation methods are to some extent ad hoc, which makes uncertainty estimation and effective communication of the results difficult. In this work, we develop a Bayesian approach to estimate susceptibility and DTB from magnetic data. We describe the subsurface in terms of tesseroids and use a two-step inversion procedure that consists of a Monte Carlo Markov Chain hyperparameter optimization and a linearized inversion. This way, the uncertainties due to unknown hyperparameter are rigorously propagated to the final maps of susceptibility and DTB. Additionally, pointwise constraints based on heat flow measurements can be easily included into the inversion. Synthetic tests are used to determine the accuracy and reliability of the new algorithm. We find that heat flow constraints are necessary to achieve reliable results, although already a small number of points is sufficient. Finally, we apply the algorithm to the Australian continent and demonstrate applicability to real data.
Curvature components derived from satellite gravity gradients provide new global views of Earth's structure. The satellite gravity gradients are based on the GOCE satellite mission and we illustrate ...by curvature images how the Earth is seen differently compared to seismic imaging. Tectonic domains with similar seismic characteristic can exhibit distinct differences in satellite gravity gradients maps, which points to differences in the lithospheric build-up. This is particularly apparent for the cratonic regions of the Earth. The comparisons demonstrate that the combination of seismological, and satellite gravity gradient imaging has significant potential to enhance our knowledge of Earth's structure. In remote frontiers like the Antarctic continent, where even basic knowledge of lithospheric scale features remains incomplete, the curvature images help unveil the heterogeneity in lithospheric structure, e.g. between the composite East Antarctic Craton and the West Antarctic Rift System.
Active source seismology provides a critical constraint on the global crustal structure. However, the heterogeneous data coverage means that interpolation is necessary to fill the gap between seismic ...profiles. This has the potential to cause large uncertainties especially if the data are interpolated over a large distance. In previous models, geological intuition was often employed to ensure reasonable results. To investigate crustal model uncertainty, we apply geostatistical analysis to a database of active seismic investigations. Unlike previous models, our workflow in the construction of the crustal model is completely transparent. Apart from the points from the database, we only use an a priori separation in oceanic and continental domains. We calculate global maps of Moho depth and average P wave velocity in the crystalline crust. Additionally, we obtain the interpolation error and error covariance. Overall, our results agree with previous global crustal models such as Crust1.0. Our uncertainty estimates show that the Moho depth uncertainty in the most well studied areas such as North America and Europe is less than 4 km but can reach 10 km or more in frontier regions such as most of Africa. P wave velocity shows the same pattern, but is less accurate overall, due to more small‐scale variation. We demonstrate the benefit of having a numerical estimate of uncertainty by propagating the uncertainty to the residual topography. We see two main uses for our crustal model in the geophysical research community: (1) as a starting model for inversions focusing on the crust and upper mantle and (2) as a starting point for including other pointwise information about crustal structure, for example, from passive seismology.
Key Points
A data‐driven approach is presented to determine crustal structure and its uncertainty on a global scale
Uncertainty of crustal thickness is less than 4 km in well‐studied areas and reaches 12 km in poorly studied regions
The Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite aimed at determining the Earth's mean gravity field. GOCE delivered gravity gradients containing directional ...information, which are complicated to use because of their error characteristics and because they are given in a rotating instrument frame indirectly related to the Earth. We compute gravity gradients in grids at 225 km and 255 km altitude above the reference ellipsoid corresponding to the GOCE nominal and lower orbit phases respectively, and find that the grids may contain additional high-frequency content compared with GOCE-based global models. We discuss the gradient sensitivity for crustal depth slices using a 3D lithospheric model of the North-East Atlantic region, which shows that the depth sensitivity differs from gradient to gradient. In addition, the relative signal power for the individual gradient component changes comparing the 225 km and 255 km grids, implying that using all components at different heights reduces parameter uncertainties in geophysical modelling. Furthermore, since gravity gradients contain complementary information to gravity, we foresee the use of the grids in a wide range of applications from lithospheric modelling to studies on dynamic topography, and glacial isostatic adjustment, to bedrock geometry determination under ice sheets.
Geothermal heat flow is an important boundary condition for ice sheets, affecting, for example, basal melt rates, but for ice-covered regions, we only have sparse heat flow observations with partly ...high uncertainty of up to 30 m W m.sup.-2 . In this study, we first investigate the agreement between such pointwise heat flow observations and solid Earth models, applying a 1D steady-state approach to perform a statistical analysis for the entire Arctic region. We find that most of the continental heat flow observations have a high reliability and agreement to solid Earth models, except a few data points, such as, for example, the NGRIP (North Greenland Ice Core Project) point in central Greenland.