Thermal structure and thickness of the lithospheric mantle underlying the Siberian Craton from the kraton and kimberlit superlong seismic profiles

Abstract

A self-consistent approach is proposed for the investigation of the thermal conditions, chemical composition, and internal structure of the upper mantle of the Earth. Using this approach, the thermal state of the lithospheric mantle beneath the Siberian Craton (SC) is reconstructed from P velocities, taking into account the phase transitions, anharmonicity, and the effects of anelasticity. The velocities of seismic waves are more sensitive to temperature than to the composition of the mantle rocks, which allows the velocity models to be effectively used for reconstruction of the thermal regime of the mantle. The temperature at depths 100–300 km is reconstructed by inversion of the Kraton and Kimberlit superlong seismic profiles for compositions of the garnet harzburgite, lherzolite, and intermediate composition of garnet peridotite. The averaged temperature in the normal continental mantle is reconstructed by inversion of the IASP91 reference model for depleted and fertile substance. One-dimensional models and two-dimensional thermal fields undergo a substantial fall in temperature (∼300–600°C) beneath the Siberian Craton as compared to the temperatures of the continental mantle and paleotemperatures inferred from the thermobarometry of xenoliths. Temperature profiles of the Siberian Craton deduced from seismic data lie between the conductive geotherms of 32.5–40.0 mW/m2 and below the P(H)-T values obtained for low- and high-temperature xenoliths from the Mir, Udachnaya, and Obnazhennaya kimberlite pipes. The thickness of the thermal lithosphere estimated from the intersection with the potential adiabat is 300–320 km, which is consistent with the data on heat flows and seismotomographic observations. This provides grounds for the assumption that the low-temperature anomalies (thermal roots of continents) penetrate down to a depth of 300 km. The analysis of the sensitivity of seismic velocity and density to the variations in temperature, pressure, and chemical and phase composition of petrological models shows that recognition of fine differences in chemical composition of the lithospheric rocks by seismic methods is impossible.