Numerical simulation of geomechanical processes in the lithosphere in terms of heat flow.

In this paper, the states of stress and strain in the deep structures in the Yakutsk-Vilyui large igneous province (LIP) are numerically analyzed. The "Craton-1980" profile with 60 km depth was chosen because it is characteristic for the Yakutsk-Vilyui LIP and contains deformed crystalline foundation and the upper mantle layer. This geological profile was obtained by deep seismic sounding surveys. Using the "Craton-1980" profile data a 2D computer structural model was constructed. To solve the posed problem, a model of the continental yield stress envelopes (YSE) was used in which the strength properties of Earth's crust layers change in accordance with depth and due to the influence of heat flow from the mantle layer. To describe the inelastic properties, the use of Nikolaevskiy's elastoplastic model was made. Numerical simulation allowed for revealing the regions in the lithosphere responsible for the formation of stress concentrations and inelastic strain localization. The location of these regions was shown to be determined not only by the heterogeneous layers in Earth's crust with curved interfaces but also by pronounced non-monotonic dependence of strength with depth. [ABSTRACT FROM AUTHOR]