Mikhail K. Kaban

EuCRUST-07: A new reference model for the European crust

[1] We present a new digital model (EuCRUST-07) for the crust of Western and Central Europe and surroundings (35N–71N, 25W–35E). Available results of seismic reflection, refraction and receiver functions studies are assembled in an integrated model at a uniform grid (15 0 � 15 0 ). The model consists of three layers: sediments and two layers of the crystalline crust. Besides depth to the boundaries, we provide average P-wave velocities in the upper and lower parts of the crystalline crust. The new model demonstrates large differences in the Moho depth compared to previous compilations, over ±10 km in some specific areas (e.g. the Baltic Shield). Furthermore, the velocity structure of the crust is much more heterogeneous than in previous maps. EuCRUST-07 offers a starting point for numerical modeling of deeper structures by allowing correction for crustal effects beforehand and to resolve trade-off with mantle heterogeneities. Citation: Tesauro, M., M. K. Kaban, and S. A. P. L. Cloetingh (2008), EuCRUST-07: A new reference model for the European crust, Geophys. Res. Lett., 35, L05313, doi:10.1029/2007GL032244.

Density of the continental roots: compositional and thermal contributions

The origin and evolution of cratonic roots has been debated for many years. Precambrian cratons are underlain by cold lithospheric roots that are chemically depleted. Thermal and petrologic data indicate that Archean roots are colder and more chemically depleted than Proterozoic roots. This observation has led to the hypothesis that the degree of depletion in a lithospheric root depends mostly on its age. Here we test this hypothesis using gravity, thermal, petrologic, and seismic data to quantify differences in the density of cratonic roots globally. In the first step in our analysis we use a global crustal model to remove the crustal contribution to the observed gravity. The result is the mantle gravity anomaly field, which varies over cratonic areas from 3100 to +100 mGal. Positive mantle gravity anomalies are observed for cratons in the northern hemisphere: the Baltic shield, East European Platform, and the Siberian Platform. Negative anomalies are observed over cratons in the southern hemisphere: Western Australia, South America, the Indian shield, and Southern Africa. This indicates that there are significant differences in the density of cratonic roots, even for those of similar age. Root density depends on temperature and chemical depletion. In order to separate these effects we apply a lithospheric temperature correction using thermal estimates from a combination of geothermal modeling and global seismic tomography models. Gravity anomalies induced by temperature variations in the uppermost mantle range from 3200 to +300 mGal, with the strongest negative anomalies associated with mid-ocean ridges and the strongest positive anomalies associated with cratons. After correcting for thermal effects, we obtain a map of density variations due to lithospheric compositional variations. These maps indicate that the average density decrease due to the chemical depletion within cratonic roots varies from 1.1% to 1.5%, assuming the chemical boundary layer has the same thickness as the thermal boundary layer. The maximal values of the density drop are in the range 1.7^2.5%, and correspond to the Archean portion of each craton. Temperatures within cratonic roots vary strongly, and our analysis indicates that density variations in the roots due to temperature are larger than the variations due to chemical differences. ; 2003 Elsevier Science B.V. All rights reserved.

Deep Europe today: Geophysical synthesis of the upper mantle structure and lithospheric processes over 3.5 Ga

abstract We present a summary of geophysical models of the subcrustal lithosphere of Europe. This includes the results from seismic (reflection and refraction profiles, P- and S-wave tomography, mantle anisotropy), gravity, thermal, electromagnetic, elastic and petro-logical studies of the lithospheric mantle. We discuss major tectonic processes as reflected in the lithospheric structure of Europe, from Precambrian terrane accretion and subduction to Phanerozoic rifting, volcanism, subduction and continent-continent collision. The differences in the lithospheric structure of Precambrian and Phanerozoic Europe, as illustrated by a comparative analysis of different geophysical data, are shown to have both a compositional and a thermal origin. We propose an integrated model of physical properties of the European subcrustal lithosphere, with emphasis on the depth intervals around 150 and 250 km. At these depths, seismic velocity models, constrained by body- and surface-wave continent-scale tomography, are compared with mantle temperatures and mantle gravity anomalies. This comparison provides a framework for discussion of the physical or chemical origin of the major lithospheric anomalies and their relation to large-scale tectonic processes, which have formed the present lithosphere of Europe.

Subcrustal density inhomogeneities of Northern Eurasia as derived from the gravity data and isostatic models of the lithosphere

For the territory of Northern Eurasia (6°E–165°W; 30–75°N) the distribution of anomalous masses in the lithosphere has been estimated in accordance with the lithosphere isostatic model. The method of model construction is based on the admittance technique. The experimental admittance presents a relation between the part of the outer load uncompensated by the Moho undulations and the residual gravity field and is used to select the best model. The 1 × 1° averaged values of topography elevations, basement and Moho depths, sedimentary cover density and gravity anomalies have been used as initial data. According to the correlation equation relating the outer load and Moho depths, the mean density contrast between the lower crust and the subcrustal lithosphere is 0.43 g/cm3, but the Moho undulation can not provide complete isostatic equilibrium. In some areas, the part of the outer load uncompensated by Moho undulations may be as large as 107 kg/m2 and the residual gravity field is as intensive as + 260 mGal. Assuming that for loads of wavelength > 200 km, local isostatic compensation is valid, in accordance with the admittance analysis, the anomalous masses compensating for the part of the outer load, which is not compensated by Moho undulations, have to be located partly in the lower crust and in the subcrustal layer. The regional trend of anomalous compensating masses is negative under Western Europe, the Mediterranean, Eastern Asia and adjacent marginal seas, and positive under the East European Platform and Western and Central Asia. The local compensating masses correspond to particular tectonic units. The isostatic gravity anomalies of Northern Eurasia have been determined and the long-wave component of the field reflecting anomalous masses under the isostatic compensation level has been evaluated.

Crust and mantle of the Tien Shan from data of the receiver function tomography

A 3-D velocity model of the Tien Shan crust and upper mantle is constructed through the inversion of the receiver functions of P and S waves together with teleseismic traveltime anomalies at nearly 40 local seismic stations. It is found that in the vast central region, where no strong earthquakes have been known over the past century, the S wave velocity at depths of 10–35 km is lower than in adjacent regions by up to 10%. These data are evidence for mechanical weakness of the crust preventing the accumulation of elastic energy. Apparently, the lower velocity and the weakness of the crust are due to the presence of water. The weakness of the crust is one of the possible reasons for the strain localization responsible for the formation of the present Tien Shan but can also be due in part to the young orogenesis. The crustal thickness is largest (about 60 km) in the Tarim-Tien Shan junction zone. The crust-mantle boundary in this region descends by a jump as a result of an increase in the lower crust thickness. This is probably due to the underthrusting of the Tien Shan by the Tarim lithosphere. This causes the mechanically weak lower crust of the Tarim to delaminate and accumulate in nearly the same way as an accretionary prism during the subduction of oceanic lithosphere. In the upper mantle, the analysis has revealed a low velocity anomaly, apparently related to basaltic outflows of the Upper Cretaceous-Early Paleogene. The Cenozoic Bachu uplift in the northern Tarim depression is also associated with the low velocity anomaly. The Naryn depression is characterized by a high velocity in the upper mantle and can be interpreted as a fragment of an ancient platform.

Density inhomogeneities, isostasy and flexural rigidity of the lithosphere in the Transcaspian region

The 3-D lithospheric-density model for the southeastern part of the Caspian Sea and the Transcaspian area, practically coinciding with the territory of the Turkmen Republic, has been constructed based on geophysical data and in accordance with the principle of isostasy. From the model selected the anomalous density of the subcrustal layer between the Moho discontinuity and the 100-km depth level is found to be — 100 kg/m3 under the Tien-Shan, − 50 kg/m3 under the Kopet-Dag mountain area, + 80 kg/m3 under the central region of the South Caspian basin, −50 kg/m3 under the eastern part of the basin, known as the West Turkmenian depression, and + 45 kg/m3 under the Murgab depression. Significant disturbances of the local isostasy are determined both in the northern and central areas of the South Caspian basin and also in the area of the Kara-Bogaz swell of the Turan platform and for the Kopet-Dag foredeep. indicating a high level of stresses in the lithosphere. The shape of the Turan plate determined by the seismic profiling is accounted for by elastic deformation resulting from the forces acting on the southern edge of the plate in the area of the Turan plate-Kopet-Dag collision. The elastic thickness of the Turan plate is estimated as 25 ± 5 km. The results obtained seem to confirm the idea that the decomposition of the Turan plate has taken place in the zone of the plates interaction and the decomposed material is situated under the Kopet-Dag ridge. We propose that the Kara-Bogaz swell is supported by the mantle material upwelling whereas the subsidence of the adjacent part of the South Caspian basin may be due to the downgoing mantle flow i.e., a small convection cell is suggested in that area.

Density structure of the lithosphere in the southwestern United States and its tectonic significance

We calculate a density model of the lithosphere of the southwestern United States through an integrated analysis of gravity, seismic refraction, drill hole, and geological data. Deviations from the average upper mantle density are as much as ±3%. A comparison with tomographic images of seismic velocities indicates that a substantial part (>50%) of these density variations is due to changes in composition rather than temperature. Pronounced mass deficits are found in the upper mantle under the Basin and Range Province and the northern part of the California Coast Ranges and adjacent ocean. The density structure of the northern and central/southern Sierra Nevada is remarkably different. The central/southern part is anomalous and is characterized by a relatively light crust underlain by a higher-density upper mantle that may be associated with a cold, stalled subducted plate. High densities are also determined within the uppermost mantle beneath the central Transverse Ranges and adjoining continental slope. The average density of the crystalline crust under the Great Valley and western Sierra Nevada is estimated to be up to 200 kg m−3 higher than the regional average, consistent with tectonic models for the obduction of oceanic crust and uppermost mantle in this region.

Density, temperature and composition of the North American lithosphere: new insights from a joint analysis of seismic, gravity and mineral physics data: 1. Density structure of the crust and upper mantle.

We introduce a new method to construct integrated 3-D models of density, temperature, and compositional variations of the crust and upper mantle based on a combined analysis of gravity, seismic, and tomography data with mineral physics constraints. The new technique is applied to North America. In the first stage, we remove the effect of the crust from the observed gravity field and topography, using a new crustal model (NACr2014). In the second step, the residual mantle gravity field and residual topography are inverted to obtain a 3-D density model of the upper mantle. The inversion technique accounts for the notion that these fields are controlled by the same factors but in a different way, e.g., depending on depth and horizontal dimension. This enables us to locate the position of principal density anomalies in the upper mantle. Afterward, we estimate the thermal contribution to the density structure by inverting two tomography models for temperature (NA07 and SL2013sv), assuming a laterally and vertically uniform “fertile” mantle composition. Both models show the cold internal part and the hot western margin of the continent, while in some Proterozoic regions (e.g., Grenville province) NA07 at a depth of 100 km is >200°C colder than SL2013sv. After removing this effect from the total mantle anomalies, the residual “compositional” fields are obtained. Some features of the composition density distribution, which are invisible in the seismic tomography data, are detected for the first time in the upper mantle. These results serve as a basis for the second part of the study, in which we improve the thermal and compositional models by applying an iterative approach to account for the effect of composition on the thermal model.

Density, temperature, and composition of the North American lithosphere—New insights from a joint analysis of seismic, gravity, and mineral physics data: 2. Thermal and compositional model of the upper mantle

Temperature and compositional variations of the North American (NA) lithospheric mantle are estimated using a new inversion technique introduced in Part 1, which allows us to jointly interpret seismic tomography and gravity data, taking into account depletion of the lithospheric mantle beneath the cratonic regions. The technique is tested using two tomography models (NA07 and SL2013sv) and different lithospheric density models. The first density model (Model I) reproduces the typical compositionally stratified lithospheric mantle, which is consistent with xenolith samples from the central Slave craton, while the second one (Model II) is based on the direct inversion of the residual gravity and residual topography. The results obtained, both in terms of temperature and composition, are more strongly influenced by the input models derived from seismic tomography, rather than by the choice of lithospheric density Model I versus Model II. The final temperatures estimated in the Archean lithospheric root are up to 150°C higher than in the initial thermal models obtained using a laterally and vertically uniform “fertile” compositional model and are in agreement with temperatures derived from xenolith data. Therefore, the effect of the compositional variations cannot be neglected when temperatures of the cratonic lithospheric mantle are estimated. Strong negative compositional density anomalies ( 92, characterize the lithospheric mantle of the northwestern part of the Superior craton and the central part of the Slave and Churchill craton, according to both tomographic models. The largest discrepancies between the results based on different tomography models are observed in the Proterozoic regions, such as the Trans Hudson Orogen (THO), Rocky Mountains, and Colorado Plateau, which appear weakly depleted (>−0.025 g/cm3 corresponding to Mg # ∼91) when model NA07 is used, or locally characterized by high-density bodies when model SL2013sv is used. The former results are in agreement with those based on the interpretation of xenolith data. The high-density bodies might be interpreted as fragments of subducted slabs or of the advection of the lithospheric mantle induced from the eastward-directed flat slab subduction. The selection of a seismic tomography model plays a significant role when estimating lithospheric density, temperature, and compositional heterogeneity. The consideration of the results of more than one model gives a more complete picture of the possible compositional variations within the NA lithospheric mantle.

Variations of the lithospheric strength and elastic thickness in North America

We evaluate the effect of temperature variations on strength and effective elastic thickness (Te) of the lithosphere of the North American (NA) continent. To this purpose, we use two thermal models that are corrected for compositional variations and anelasticity effects in the upper mantle. These thermal models are obtained from a joint inversion of gravity data and two recent seismic tomography models (NA07 and SL2013sv). The crustal rheology was defined using NACr14, the most recent NA crustal model. This model specifies seismic velocities and thickness for a three-layer model of the crystalline crust. Strength in the lithosphere and in the crust has similar distributions, indicating that local geotherms play a dominant role in determining strength rather than crustal composition. A pronounced contrast is present in strength between cratonic and off-cratonic regions. Lithospheric strength in the off-cratonic regions is prevalently localized within the crust and Te shows low values ( 150 km). In contrast to previous results, our models indicate that Phanerozoic regions located close to the edge of the cratons, as the Appalachians, are characterized by low strength. We also find that locally weak zones exist within the cratons (e.g., beneath the intracratonic Illinois Basin and Midcontinent rift). Seismic tomography models NA07 and SL2013sv differ mainly in some peripheral parts of the cratons, as the Proterozoic Canadian Platform, the Grenville, and the western part of the Yavapai-Mazatzal province, where the integrated strength for the model NA07 is 10 times larger than in model SL2013sv due to a temperature difference (>200°C) in the uppermost mantle. The differences in Te between the two models are less pronounced. In both models, Proterozoic regions reactivated by Meso-Cenozoic tectonics (e.g., Rocky Mountains and the Mississippi Embayment) are characterized by a weak lithosphere due to the absence of the mechanically strong part of the mantle lithospheric layer. Intraplate earthquakes are distributed along the edges of the cratons, demonstrating that tectonic stress accumulates there, while the cores of the cratons remain undeformed. In both models, intraplate earthquakes occur in weak lithosphere (∼0.5 × 1013 Pa s, Te ∼ 15 km) or near the edges of strong cratonic blocks, characterized by pronounced contrasts of strength and Te.