A. A. Baranov

A new crustal model for Central and Southern Asia

We present a new regional three-layer crustal model for the Central and Southern Asia and surroundings (AsCRUST-08). The model provides Moho boundary, thickness of different layers of consolidated crust and P-velocity distribution in these layers. A large volume of new data on seismic reflections and refractions as well as on surface waves generated by earthquakes or blasts was analyzed. All these data were incorporated into a unified digital 3D integrated model with 1° × 1° resolution. Results are represented as seven numerical maps imaging the distributions of the Moho depth, the thickness of the upper, middle, and lower layers of the consolidated crust, and the P-wave velocities therein.

Exact analytical solutions of the stokes equation for testing the equations of mantle convection with a variable viscosity

Presently, the study of the mantle flow structure is mainly based on numerical modeling. The most important stage of the development of a computer program is its testing. For this purpose, results of various test models of convection flows with a given set of parameters are compared. The solution of the Stokes equation, involving the derivative of viscous stresses, is most difficult. Exact analytical solutions of the Stokes equation are obtained in this work for various cases of special loads. These solutions can be used as benchmarks for testing programs of numerical calculation of viscous flows in both geophysics and engineering. The advantage of this testing technique is the exceptional simplicity of the solution form, the admissibility of any spatial viscosity variations, and the fact that solutions can be compared not for a narrow set of the solution parameters but for any distributions of velocities, viscous stresses, and pressures at all points of the space.

Numerical models of subduction of the oceanic crust with basaltic plateaus

Light continents and islands characterized by a crustal thickness of more than 30 km float over a convective mantle, while the thin basaltic oceanic crust sinks completely in subduction zones. The normal oceanic crust is 7 km thick. However, anomalously thick basaltic plateaus forming as a result of emplacement of mantle plumes into moving oceanic lithospheric plates are also pulled into the mantle. One of the largest basaltic plateaus is the Ontong Java plateau on the Pacific plate, which arose during the intrusion of a giant superplume into the plate ∼100 Myr ago. Notwithstanding its large thickness (averaging ∼30 km), the Ontong Java plateau is still experiencing slow subduction. On the basis of numerical modeling, the paper analyzes the oceanic crust subduction process as a function of the mantle convection vigorousness and the density, thickness, viscosity, and shape of the crust. Even a simplified model of thermocompositional convection in the upper mantle is capable of explaining the observed facts indicating that the oceanic crust and sediments are pulled into the mantle and the continental crust is floating on the mantle.

The effect of spatial variations in viscosity on the structure of mantle flows

According to an opinion widespread in the literature, high viscosity regions (HVRs) in the mantle always affect the structure of mantle flows, changing it in both the HVR itself and the entire mantle. Moreover, a simplified relation is often adopted according to which the flow velocity in the HVR decreases in inverse proportion to viscosity. Therefore, in order to treat a smoother value, some authors introduce a new variable equal to the product of the flow velocity and the viscosity value in a given place. On the basis of numerical modeling, this paper shows that HVRs of two types should be distinguished in the mantle. If an HVR is immobile, mantle flows actually do not penetrate it. If the viscosity increase is more than five orders, the HVR behaves as a solid and flow velocities within it almost vanish. However, if an HVR is free, it moves together with the mantle flow. Then, the general structure of flows changes weakly and flow velocities within the HVR become approximately equal to the average velocity of flows in the absence of the HVR. Horizontal layers and vertical columns differing in viscosity from the mantle behave as regions of the first type, whose flow velocities can differ by a few orders. However, even such large-scale regions as the continental lithosphere, whose viscosity is four to five orders higher than in the surrounding mantle, float together with continents at velocities comparable to mantle flows, i.e., behave as regions of the second type.

Phase Transition Zone Width Implications for Convection Structure

A temperature and pressure increase in the mantle causes phase transitions and related density changes in its material. The transition boundary in the pressure-temperature phase diagram is determined by the curve of phase equilibrium with the slope γ = dp/dT. If the slope is nonzero, a phase transition in hot ascending and cold descending mantle flows occurs at different depths and, therefore, either enhances (γ > 0) or slows down convection (γ < 0). The mantle material has a multicomponent composition. Therefore, phase transitions in the mantle are distributed over an interval of pressures and depths. In this interval, the concentration of one phase smoothly decreases and the concentration of the other increases. The widths of phase transition zones in the Earth’s mantle vary from 3 km for the endothermic transition in olivine at a depth of 660 km to 500 km for the exothermic transition in perovskite, and the high-to-low spin change in the atomic state of iron takes place at a depth of about 1500 km. This work presents results of calculations demonstrating the convection effect of phase transitions as a function of the transition zone width. Transitions of both types with different slopes of the phase curve and different intensities of mantle convection are examined. For the first time, the convection enhancement and an increase in the mass transfer across the phase boundary are quantitatively investigated in the presence of an exothermic phase transition as a function of the slope of the phase curve. The mixing of material under conditions of partially layered convection is examined with the help of markers.

Convection in the Mantle with an Endothermic Phase Transition

An endothermic phase transition at a depth of 660 km in the mantle partially slows down mantle flows. Many models considering the possibility of temporary layering of flows with separation of convection in the upper and lower mantle have been constructed over the past two decades. The slowing-down effect of the endothermic phase transition is very sensitive to the slope of the phase-equilibrium curve. However, laboratory measurements contain considerable uncertainties admitting both a partial convection layering and only an insignificant slowing down of a part of downgoing mantle flows. In this work, we present results of calculations of mantle flows within a wide range of phase-transition parameter values, determine ranges of one-and two-layer convection, and derive dependences of the amplitude and period of oscillations on phase-transition parameters.

Influence of an endothermic phase transition on mass transfer between the upper and the lower mantle

Geochemical data indicate that two major reservoirs 1–2 Ga in age are present in the mantle. The upper mantle, feeding mid-ocean ridges, is depleted in chemical elements carried away into the continental crust. The lower mantle, feeding hotspot plumes, is close in composition to primordial matter. The 660-km depth of an endothermic phase transition in olivine has been considered over the last two decades as a possible boundary between the reservoirs. In this period, many models of mantle convection were constructed that used values of the phase transition parameters, which led to temporal (up to 1 Gyr long) convection layerings and periodic avalanche-induced mantle intermixing events. However, laboratory measurements with new improved instrumentation give other values of the phase transition parameters that require a revision of the majority of the existence of large-scale avalanches in the Earth’s history becomes disputable. The paper is devoted to comprehensive study of the phase transition effect on the structure of mantle flows with different values of phase transition parameters and Rayleigh numbers; in particular, the mass transfer through the phase boundary is calculated for different regimes of steady-state convection.

The structure of 2D mantle convection and stress fields: Effects of viscosity distribution

Spatial fields of temperature, velocity, overlithostatic pressure, and horizontal stresses in the Earth’s mantle are studied in two-dimensional (2D) numerical Cartesian models of mantle convection with variable viscosity. The calculations are carried out for three different patterns of the viscosity distribution in the mantle: (a) an isoviscous model, (b) a four-layer viscosity model, and (c) a temperature- and pressure-dependent viscosity model. The pattern of flows, the stresses, and the surface heat flow are strongly controlled by the viscosity distribution. This is connected with the formation of a cold highly viscous layer on the surface, which is analogous to the oceanic lithosphere and impedes the heat transfer. For the Rayleigh number Ra = 107, the Nusselt number, which characterizes the heat transfer, is Nu = 34, 28, and 15 in models with constant, four-layered, and p, T-dependent viscosity, respectively. In all three models, the values of overlithostatic pressure and horizontal stresses σxx in a vast central region of the mantle, which occupies the bulk of the entire volume of the computation domain, are approximately similar, varying within ±5 MPa (±50 bar). This follows from the fact that the dimensionless mantle viscosity averaged over volume is almost similar in all these models. In the case of temperature- and pressure-dependent viscosity, the overlithostatic pressure and stress σxx fields exhibit much stronger concentration towards the horizontal boundaries of the computation domain compared to the isoviscous model. This effect occurs because the upwellings and downwellings in a highly viscous region experience strong variations in both amplitude and direction of flow velocity near the horizontal boundaries. In the models considered with the parameters used, the stresses in the upper and lower mantle are approximately identical, that is, there is no denser concentration of stresses in the upper or lower mantle. In contrast to the overlithostatic pressure field, the fields of horizontal stresses σxx in all models do not exhibit deep roots of highly viscous downwelling flows.

Horizontal stresses in the mantle and in the moving continent for the model of two-dimensional convection with varying viscosity

Numerical experiments on studying the spatial fields and evolution of viscous overlithostatic horizontal stresses and pressure in the mantle and in the moving continent are carried out. The continent moves consistently with time-dependent forces, which act from the viscous mantle. By introducing the varying viscosity, we gain the possibility for taking into account the oceanic lithosphere and the difference between the viscosity of the upper and the lower mantle in the context of a purely viscous model. The typical overlithostatic horizontal stresses in the main part of the mantle are ±(7–9) MPa (70–90 bar); in the highly viscous regions and, particularly, in the subduction zones they are at least three times larger. The descending mantle flows in the depth interval from approximately 50 km to about 300 km are more sharply pronounced in the pressure field than in the field of horizontal stresses. At the considered stages of motion and in different parts, the continent is characterized by the following typical values of stresses: the overlithostatic pressure ranges from −5 to +15 MPa; the horizontal overlithostatic tensile stress amounts up to −4MPa (−40 bar); and the compressive stress in case of the overriding of the subduction zone attains +35 MPa (350 bar).

Moho Depth in Antarctica from Seismic Data

A refined digital model of the Moho depth is constructed for the Antarctica on a uniform grid with resolution of 1° × 1°. The model is based on seismic data. Results are presented as a digital table that defines the Moho depth (the Moho, or M) at each point. A large volume of new data on reflection, refraction, converted and surface waves, as well as receiver functions and data on subglacial relief, were analyzed. The new model provides far more precise and detailed information about the Moho than the previous model. The difference in the crustal thickness between these two models may amount up to −10–±24 km.