Observed Gravity Field Research Articles

New 2.75-D gravity modeling reveals the low-density nature of propagator wakes in the Juan de Fuca plate

Propagator wakes within the Juan de Fuca plate represent zones of oceanic crust affected by dynamic readjustments of the spreading centers. They are traditionally mapped from distortions of magnetic anomalies and have been previously interpreted as zones of denser than regular oceanic crust from gravity analysis along several seismic lines in 2-D approximation. In this paper, we utilized a 2.75-D modeling assumption that allowed us to incorporate the effects of the nearby seamounts that were not accounted for in the previous 2-D study. Our analysis reveals that the positive gravity effect of those geological objects was misinterpreted as a need for higher crustal density within propagator wake zones. When seamounts are included in the model, the crust of propagator wakes requires lower densities than the surrounding oceanic crust to explain the observed gravity field. We postulate that the lower densities of propagator wakes relate to faulting during the readjustment of the spreading centers, suggesting that they represent zones of weaker crust.

Gravity modelling of the Tyrrhenian-Calabrian-Ionian subduction system

This study provides a lithospheric-scale model along the Ionian Subduction zone in Southern Italy, contributing to the seismotectonic investigation of a region which is affected by relevant historical seismicity. The study employs gravity forward modelling to build the geo-structural model along a composite, NW-SE trending transect extending from the Ionian to the Tyrrhenian Sea, including the Aeolian arc and the Calabro-Peloritan onshore. Through a multidisciplinary approach, we propose new interpretations of three 2D deep-seismic reflection profiles across the study area. Such interpretative profiles are used as constraints to model the observed Bouguer gravity anomalies providing upper and lower crust geometries. Whilst a tomographic model provides constraints for the lithospheric and asthenospheric modelling. The entire workflow is constrained by literature data about Moho geometry, deep seismicity and tomographic images that are integrated to determine the subducting slab geometry. The proposed model of the entire subducting system reasonably fits the observed gravity field and is coherent with the first-order geological and geophysical constraints. The modelling results in updated Tyrrhenian and Ionian Moho depth, subducting slab geometry and location, and densities of the main units, providing valuable input about the composition and geometry of the Calabrian arc structures.

A Novel Approach for Bathymetry Estimation through Bayesian Gravity Inversion

The bathymetry is the most superficial layer of the Earth’s crust on which it is possible to perform direct measurements. However, it is also well known that water covers more than 70% of the Earth’s surface, so an enormous expenditure of acquisition campaigns should be performed to produce a high-resolution map of this layer. Currently exploiting mainly commercial navigation routes, the sea floor coverage with shipborne sounding is only at 11%, and the remaining part is currently modeled by classical interpolation techniques or satellite-based gravity inversion methods. In the present work, a new method to refine bathymetry modeling based on the exploitation of global gravity field models is presented. In the proposed solution, once modeled and removed from the observed gravity field, the gravitational signals related to the deepest structures, a 3D Bayesian inversion algorithm is used to improve the actual knowledge of bathymetry. The proposed inversion method also enables limiting the solution to shipborne sounding measurements in such a way as to improve the seafloor grid where no local, high-quality information is available. Two test cases are discussed in the Mediterranean Sea region. Promising results are presented, opening the possibility of applying an analogous method to refine the bathymetry modeling at larger scales up to the global one.

3-D density structure of the upper-mantle in the eastern Mediterranean Sea and surrounding region using gravity inversion constrained by seismic velocity model

A 3D density structure of the lithosphere and upper mantle beneath the eastern Mediterranean Sea (EMS) and its adjacent region was constructed based on gravity anomaly inversion constrained by a seismic tomography model. Gravity effects of terrain and crust were removed from the observed gravity field (EIGEN-6C4) to obtain the residual mantle gravity anomaly (RMGA). The density distribution of the lithosphere and upper mantle was investigated. The 3D inversion process was constrained by an initial density model projected from the shear-wave velocity model (SL2013sv). The results show some characteristics of the density distribution in the lithosphere and in the upper mantle that could be related to the tectonic importance of the Mediterranean Sea and the surrounding region. A low-density zone dominates the lithosphere beneath the Sea except for the area around the Arabia Shield and North Anatolian Fault belt. A thinner, high-density layer beneath the southwest of the Sea may be related to the older oceanic lithospheric fragments. The high-density anomalies appear at depths below 280 km beneath the Sea and the Turkish Aegean Sea Plate. However, the low-density anomalies appear on the upper mantle under the trenches of the southwestern part of the Mediterranean Sea, the eastern part of the Aegean Sea, the Red Sea, the Black Sea, and the middle of the Arabia shield. The deep structure under the Eratosthenes seamount in the Mediterranean Sea is the source of the intensity and genesis of tectonic activity. Furthermore, the convergence region of two low-density anomaly zones (Africa-Anatolia) may be interpreted as a significant tectonic unit (Eratosthenes seamount) caused by the arrival of the relatively thick and buoyant Eratosthenes block to its present location south of Cyprus in Holocene time based on the density model interpretation beneath the Mediterranean Sea during the Late Cretaceous and early Tertiary period.

Lithospheric Structure of the Red Sea Based on 3D Density Modeling: A Contrasting Rift Architecture

AbstractThe Red Sea is an ideal location for studying rifting processes, offering a young and active intra‐continental rift at the final stages of breakup. We explore the 3D lithospheric structure of the Red Sea by analyzing the gravity response of four end‐member models of rift architecture, including two end‐member types margin architecture Type I—narrow margins and exhumed lithospheric mantle, and Type II—(ultra)wide margins and removal of lithospheric mantle. Additionally, we test two options for the distribution of the oceanic crust (a) limited, that is, confined only to regions of magnetic stripes, and (b) extended, that is, available in vast areas within the basin. South of latitude 23°N, our results suggest the presence of ultrawide margins and limited oceanic crust. North of this latitude, the model of exhumed mantle and limited oceanic crust has minimized residuals compared to the observed gravity field, and agrees with a previously published regional tomographic model. Additionally, we find evidence for the presence of a high‐density body along the southern Arabian coast, probably associated with magmatic underplating. We discuss the lithospheric structure of the Red Sea with respect to the nearby Afar plume, agreeing that the close proximity of the central‐southern regions to the plume promoted a reduction in the strength of the lithosphere, and led to the development of ultrawide margins in these areas.

Choosing a Regularization Parameter in the Problem of Analytical Continuation of Gravitational Fields (Separation of Anomalies Generated by Shallow and Deep Sources)

Abstract—We present a filtration tomography technique for isolating components of the gravity field anomalies generated by inhomogeneities of the horizontally layered density model. The filtration algorithm of field separation relies on the solution of the forward and inverse problems of analytical continuation of harmonic functions through the horizontal boundary plane. We applied the regularizing algorithms to analytical continuation of the gravity field “down” to its generating sources. The fields successively recalculated upward and downward relative to preset depths allowed us to partition the initial (total) field as the sum of the fields generated in the layers based on the properly selected adaptive regularization parameter α. For the sake of stability of the inverse problem solution in the analytical continuation of the observed gravity field to a certain depth, we used the Lavrentiev regularization scheme involving the L-curve method (for selecting the adaptive regularization parameter). The smoothing regularization parameter values obtained from the preset successive depth intervals and grid step for the observed field are shown to be optimal for dividing the observed field into components corresponding to different depths. The developed algorithms for massively parallel computing systems and their application to a group of different heights were numerically implemented on the Uran supercomputer.

Marine Geophysical Data and Its Application to Assessment of Crustal Structure along the Northern Egyptian Passive Continental Margin

The Egyptian passive continental margin is considered a remarkable geologic setting, in addition to being an occupation and manufacturing locality in north Egypt. This work used accessible potential field data, such as marine gravity data, to provide a wider vision of the potential field of the area under investigation. The results of the two-dimensional (2D) gravity modelling revealed a good agreement between the modelled gravity and observed gravity fields, including known regional structures found in the investigated area. The findings revealed that crustal modelling was affected by the tectonic structure and the huge thickness of sedimentary layers, which act as barriers to the crystal crust. The results revealed that the crustal thickness and density are spread among the deposited layer and the inferior mantle in the Moho range. Furthermore, it was found that the basement extent lies nearly 6–9 km lower in the northern Egyptian coastline to approximately 13 km under the Herodotus abyssal plain. Moreover, it was shown that the thickness of the sedimentary layers deposit increases near the East Mediterranean Ridge.

ТРЕХМЕРНАЯ ПЛОТНОСТНАЯ МОДЕЛЬ ВЕРХНЕЙ КОРЫ В ОБЛАСТИ СОЧЛЕНЕНИЯ ЛОСЕВСКОГО И ВОРОНЦОВСКОГО ТЕРРЕЙНОВ (ВОРОНЕЖСКИЙ КРИСТАЛЛИЧЕСКИЙ МАССИВ)

The results of detailed three-dimensional density modeling of the upper part of the Earth's crust for the area located in the northeastern part of the Voronezh crystalline massif in the area of the junction of the Losevsky and Vorontsovsky terranes are considered. The model obtained by solving the inverse gravimetry problem characterizes the density distribution of the rocks of the crystalline basement of the territory to a depth of 16 km. An important role in such constructions is given to the starting model of the studied environment, which was formed taking into account the available a priori geological and geophysical information: the regional density model of the East European Platform lithosphere and the gravity field, the data on the density of rocks of the sedimentary cover and the crystalline basement, estimates of the “gravity” layer thickness, the geological map of the crystalline basement, and the digital elevation models of the region. The three-dimensional density model allowed specifying the geological structure of the upper crust of the study area and tracing the deep continuation of the surface structures of the basement. The reliability of the obtained results is determined by the consistency of the accepted a priori data on the environment and the consistency of the model and observed gravity fields of the study area.

Numerical Investigation of Lunar Basin Formation Constrained by Gravity Signature

AbstractImpact basins on the Moon can serve as a benchmark for timing and intensity of the impact flux in the inner solar system. The basin morphology and morphometry depend on impactor size, mass, and velocity as well as the thermal state of the lunar lithosphere which is a function of the cooling history. Erosion by superimposed impact bombardment has altered the surface expression of basin structures over time, making it difficult to determine the size unequivocally solely based on topographic expression. The gravity signature of basins is thought to be a less altered measure of the size of impact structures. By a systematic study of basin formation using the iSALE‐2D shock‐physics code, we investigate the influence of the lunar thermal state and different impactor properties on the transient crater and final basin size and on the resulting gravity anomaly. As constraints we use gravity data of 16 farside basins and their assumed formation ages to estimate the subsurface temperature related to the cooling history of the Moon. Our modeling results confirm that the thermal state affects the basin formation process and the basin sizes significantly. We provide quantitative relationships between the observed gravity signal, the different basin sizes, and the impactor diameter considering the thermal state of the Moon upon impact, which correlates with the formation ages or periods in the literature. Our study allows for estimating the impactor size from the observed gravity field if the formation age and, thus, the thermal state of the lithosphere is approximately known.

Forward Modelling of Bouguer Anomalies along a transect of the Southern Apennines and the Southern Tyrrhenian Sea (Italy)

In the present study, we perform a gravity modelling at crustal scale along the trace of the CROP-04 (on-shore) and M-6B (off-shore) deep seismic reflection profiles, crossing the Southern Apennines and the Southern Tyrrhenian Sea (Italy). Along the 321 km-long modelled profile, we investigate the crustal-scale sources for the observed gravity anomalies through a simplified model of the crust and upper mantle across both onshore and offshore areas.After a compelling review of the published Moho geometries in the area, that were retrieved from either active or passive seismic methods, we test them in the observed gravity field through forward modelling of the Bouguer gravity anomalies. The comparison between the different Moho interpretations highlights the major contributors to the observed Bouguer gravity at the crustal scale, defining a set of starting values of these parameters for our final model.The proposed model locates the westward flexure of the Adriatic Moho, mimicking the subduction of the Adriatic lithosphere beneath the Peri-Tyrrhenian block and locates the step between the western (Tyrrhenian) and the eastern (Adriatic) Moho beneath the Apennines range providing a valuable geometrical and compositional model at the crustal scale. The model depicts a typical oceanic-to-continental crust transition in the Tyrrhenian domain and represents a solid starting base for further detailed modelling across the area.

Integrated geological and gravity modelling to improve 3-D model harmonization—Methods and benefits for the Saxony-Anhalt/Brandenburg cross-border region (North German Basin)

SUMMARYAs 3-D geological models become more numerous and widely available, the opportunity arises to combine them into large regional compilations. One of the biggest challenges facing these compilations is the connection and alignment of individual models, especially in less explored areas or across political borders. In this regard, gravity modelling is suitable for revealing additional subsurface information that can support a harmonization of structural models. Here, we present an integrated geological and gravity modelling approach to support the harmonization process of two geological 3-D models of the North German Basin in the cross-border region between the federal states of Saxony-Anhalt and Brandenburg. Gravity gradient calculation, filtering and Euler deconvolution are utilized to reveal new insights into the local fault system and gravity anomaly sources. The independent models are merged and harmonized during 3-D forward and inverse gravity modelling. Herein, density gradients for individual layers are incorporated in the framework of model parametrization. The resulting geological 3-D model consists of harmonized interfaces and is consistent with the observed gravity field. To demonstrate the plausibility of the derived model, we discuss the new geophysical findings on the sedimentary and crustal structures of the cross-border region in the context of the regional geological setting. The cross-border region is dominated by an NW–SE oriented fault system that coincides with the Elbe Fault System. We interpret a low-density zone within the basement of the Mid-German Crystalline Rise as a northward continuation of the Pretzsch–Prettin Crystalline Complex into the basement of the North German Basin. Additionally, we observe two types of anticlines within the basin, which we link to provinces of contrasting basement rigidity. Our gravity modelling implies that the Zechstein salt has mostly migrated into the deeper parts of the basin west of the Seyda Fault. Finally, we identify a pronounced syncline that accommodates a narrow and up to 800 m deep Cenozoic basin.

Moho Beneath Tibet Based on a Joint Analysis of Gravity and Seismic Data

AbstractWe use the improved Parker‐Oldenburg's formulas that include a reference depth into the exponential term and employ the Gauss‐fast Fourier transform method to determine Moho depth beneath the Tibetan Plateau. The synthetic models demonstrate that the improved Parker's formula has high accuracy with the maximum absolute error less than 0.25 mGal compared to the analytical solution. Two inversion parameters, that is, the reference depth and the density contrast, are essential for the Moho estimation based on the gravity field, and they need to be determined in advance to obtain correct results. Therefore, the Moho estimates derived from existing seismic studies are used to reduce the nonuniqueness of the gravity inversion and to determine these parameters by searching for the maximum correlation between the gravity‐inverted and seismic‐derived Moho depths. Another critical issue is to remove beforehand the gravity effects of other factors, which affect the observed gravity field besides Moho variations. In addition to the topography, the gravity effects of the sedimentary layer and crystalline crust are removed based on existing crustal models, while the upper mantle impact is determined based on the seismic tomography model. The inversion results show that the Moho structure under the Tibetan plateau is very complex with the depths varying from about 30–40 km in the surrounding basins (e.g., Ganges basin, Sichuan basin, and Tarim basin) to 60–80 km within the plateau. This considerable difference up to 40 km on the Moho depth reveals the substantial uplift and thickening of the crust in the Tibetan Plateau. Furthermore, two visible “Moho depression belts” are observed within the plateau with the maximum Moho deepening along the Indus‐Tsangpo Suture and along the northern margin of Tibet bounding the Tarim basin with the relatively shallow Moho in central Tibet between them. The southern “belt” is likely formed in compressional environment, where the Indian plate underthrusts northward beneath the Tibetan Plateau, while the northern one could be formed by the southward underthrust of the Asian lithosphere beneath Tibet.

Density distribution across the Alpine lithosphere constrained by 3-D gravity modelling and relation to seismicity and deformation

Abstract. The Alpine orogen formed as a result of the collision between the Adriatic and European plates. Significant crustal heterogeneity exists within the region due to the long history of interplay between these plates, other continental and oceanic blocks in the region, and inherited crustal features from earlier orogenies. Deformation relating to the collision continues to the present day. Here, a seismically constrained, 3-D structural and density model of the lithosphere of the Alps and their respective forelands, derived from integrating numerous geoscientific datasets, was adjusted to match the observed gravity field. It is shown that the distribution of seismicity and deformation within the region correlates well to thickness and density changes within the crust, and that the present-day Adriatic crust is both thinner and denser (22.5 km, 2800 kg m−3) than the European crust (27.5 km, 2750 kg m−3). Alpine crust derived from each respective plate is found to show the same trend, with zones of Adriatic provenance (Austro-Alpine unit and Southern Alps) found to be denser and those of European provenance (Helvetic zone and Tauern Window) to be less dense. This suggests that the respective plates and related terranes had similar crustal properties to the present-day ones prior to orogenesis. The model generated here is available for open-access use to further discussions about the crust in the region.

A numerical model for the gravimetric recovery of sub-lithospheric mantle structures

It is a well-known fact that the long-wavelength terrestrial geoid undulations are mainly attributed to deep mantle density heterogeneities, while more detailed features in the geoid geometry are associated with the topography and the lithospheric density structure. To enhance a gravitational signature of mantle density heterogeneities below the lithosphere, the gravitational contributions of topography and lithospheric density heterogeneities should be modelled and subsequently removed from Earth's gravity field. The refined gravity field obtained after this numerical procedure is more suitable for a recovery of a mantle density structure (below the lithosphere). Following this idea, methods for a spherical harmonic analysis and synthesis of gravity field and lithospheric density structures are presented, and a theoretical relation between gravity field and mass density structure is formulated. Since a gravimetric recovery of inner density structure has a non-unique solution, we propose an alternative method based on a conversion of seismic velocities to mass densities. A forward modelling approach is then employed to find the mantle density configuration that generates the gravitational field that best approximates the corresponding refined gravitational field obtained from observed gravity field after subtracting the gravitational signal of the lithosphere.

3-D crustal density model of the Sea of Marmara

Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible “end-member” solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kg m−3) and an intermediate to mafic lower crystalline crust (average density of 2890 kg m−3) appear to be cross-cut by two large, dome-shaped mafic high-density bodies (density of 2890 to 3150 kg m−3) of considerable thickness above a rather uniform lithospheric mantle (3300 kg m−3). The spatial correlation between two major bends of the main Marmara fault and the location of the high-density bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.

Imitation modelling technology for gravity inversion cases

Formulation of the problem. A gravity method is aimed at prospecting and exploration of mineral resources which are based on the study of the geological section structure. The task of quantitative interpretation of the gravimetric materials, which uses methods for solving direct and inverse gravity problems, is the modelling of a gravity field (direct problem) and geological media’s density structure (inverse problem). The important features of methods for density structure modelling of complex geological media are geological content, consistency with a priori data and its subordination to geological hypotheses. It is proposed to analyze these properties by a imitation technique. The purpose of the article is to describe the imitation gravimetric modelling method, based on the construction of an informal sequence of equivalent solutions. The purpose of imitation modelling is to study the properties of gravity inversion in general formulation as well as to assess the degree of detail and reliability of the methodology and technologies of gravity modelling, which is claimed to be an effective solution to geological problems. Methods. Imitation modelling technology and methods of solving gravity direct and inverse problems for geodensity model of complex geological environment. Results. The examples of density and structural simulation testing of the informal sequence of equivalent solutions and its computer technologies show that complex interpretation of wells, seismic and gravity data enables to create detailed density models of geological medium. Studies have also been conducted of ways to increase the reliability of gravitational modelling. Scientific novelty and practical significance. It is revealed that the best approximation of the regional background is an inclined plane, which approximates the observed gravity field along characteristic pickets over the research areas that are better studied. Also, an increase in the reliability of modelling can be achieved by rebuilding near side zones in structural type models in an interactive process of solving structural inverse gravity problems. Substantive modelling depends primarily on the experience of the interpreter, since computer technologies for solving direct and inverse gravity problems are only an interpretation tool.

A map of absolute values for gravity fields of Ukraine and some aspects of its possible interpretation

An electronic map of absolute values of gravity field ( G н ) has been plotted for the territory of Ukraine for the first time with application of combined electronic matrices of Buge anomalies in case of the density of intermediate layer 2.3 g/cm 3 and the values of altitude marks (the Baltic system). The size of unit cell of the matrix is 1´1 km, the error obtained for the map interpolation G н for flat country is ±2,1 mGal. It has been shown that the maps of G н give more complete information on tectonic processes, taking part in the Earth’s crust formation because they reflect combined effects of density and relief-producing factors, which are not taken into account while maps of Buge anomalies are being analyzed. Examples have been given of large lineaments of tectonic origin, clearly visible on G н map and practically unnoticeable on the map of Buge anomalies. Comparison of anomalies of G н and Buge (2,3 g/cm 3 ) in mountain areas (Crimea, Carpathians) gives an idea of hypogenous nature of gravity objects, and large regional differences of anomalies G н , which are as big as 950 mGal within the limits of Ukraine must be taken into account while problems of geophysical ecology and compiling meteorology forecasts are being solved. A necessity is being indicated to introduce maps of G н field to a set of gravimetric materials used for geological-geophysical interpretation. In this connection we propose to introduce in the running regulations for gravity works a demand on presenting a map of observed gravity field as a report.

High-resolution regional gravity field model of Pakistan based on best residual terrain model (RTM) of topography and interpolation techniques

This study focuses on the development of absolute gravity model for Pakistan based on best possible residual terrain model of gravity using residual terrain modeling technique. The datasets used for the development of model are observed gravity, global gravity models, and Shuttle Radar Topographic Mission (SRTM30) elevation data. The residual terrain modeling technique has been used in the remove-restore procedure for smoothing the observed gravity field. Different topographic elevation models were tested in the model selection and one best possible model with minimum mean and standard deviation was selected for residual terrain effects. Least square collocation technique has been used for quality control and error estimates. The best possible covariance model was established from residual gravity for onward prediction of gravity anomalies at the earth surface for error and prediction analysis. The residual terrain effect of gravity, value of free air anomaly from EGM96, and observed free air anomaly are added to normal gravity to compute the absolute gravity at earth surface. The prediction of these parameters is made by employing Lagrange interpolation with least square adjustment. The results are compared with ∼5% randomly selected data points not utilized for the development of covariance function and/or model development. Spline interpolation technique has also been used for the prediction of gravity field-related parameters. Lagrange interpolation exhibits relatively superior results over spline-based interpolation. This is as per expectation due to the reason that additional gridding for spline interpolation filters the signal part as well. This fact is evident from the results of spline interpolation of Grid-I and Grid-II with relatively better prediction results in Grid-I. This version of the model is capable of prediction having limiting error of 30 mGal. The predicted results show that 96.16% of prediction data falls within above-mentioned limit with Lagrange interpolation technique with least square adjustment for whole Pakistan area. The adverse effect of gridding is absent in case of Grid-I due to relatively flat areas and predicted data matches totally with control values for both spline as well as Lagrange interpolations. However, in case of Grid-II which includes high mountains of Himalaya, gridding effect is present and the accuracy of the predicted results falls to ∼92%. The computed results have been compared with absolute values predicted using EGM96 and EGM2008 models as well. The gravity field recovered with PAKGM model is much better, i.e., ∼ 96.16%, than both with EGM96 and EGM2008 which is about 85% only.

Integrated crustal–geological cross-section of Ellesmere Island

Abstract The crustal seismic velocity model (based on receiver functions) of Ellesmere Island and the structural geological cross-section of Ellesmere Island, both documented and discussed elsewhere in this volume, are here integrated into a crustal-scale transect crossing all the main tectonic domains. The velocity model satisfies much of the observed gravity field, but implies minor modifications with potentially important implications for characterizing the lower crust over the transect. The crust of the Pearya Terrane includes a high-velocity and high-density lower crustal body, suggested to represent a mafic underplate linked to the emplacement of the High Arctic Large Igneous Province. A similar body also lies directly beneath the Hazen Plateau, but this is more likely to be inherited from earlier tectonic stages than to be linked to the High Arctic Large Igneous Province. A large-scale basement-involving thrust, possibly linked to a deep detachment of Ellesmerian age, lies immediately south of the Pearya Terrane and forms the northern backdrop to a crustal-scale pop-up structure that accommodates Eurekan-aged shortening in northern Ellesmere Island. The thickest crust and deepest Moho along the transect are below the Central Ellesmerian fold belt, where the Moho is flexured downwards to the north to a depth of about 48 km beneath the load of the structurally thickened supracrustal strata of the fold belt.

The principal moments of inertia calculated with the hydrostatic equilibrium figure of the Earth

As an indication of the Earth's mass distribution, the principal moments of inertia (PMOI, i.e., A, B, C) of the Earth are the basic parameters in studies of the global dynamics of the earth, like earth nutation, and the geophysics. From the aspect of observation, the PMOI can be calculated from the spherical coefficients of observed gravity field. In this paper, the PMOI are calculated directly according to its definition with the figures of the Earth's interior derived by a generalized theory of the hydrostatic equilibrium figure of the Earth. We obtain that the angle between the principal axis of the maximum moment of PMOI and the rotational axis is 0.184°, which means that the other two principal axes are very closely in the equatorial plane. Meanwhile, B-A is 1.60 × 10−5MR2, and the global dynamical flattening (H) is calculated to be 3.29587 × 10−3, which is 0.67% different from the latest observation derived value Hobs(3.273795 × 10−3) (Petit and Luzum, 2010), and this is a significant improvement from the 1.1% difference between the value of H derived from traditional theories of the figure of the Earth and the value of Hobs. It shows that we can calculate the PMOI and H with an appropriate accuracy by a generalized theory of the hydrostatic equilibrium figure of the Earth.