Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics

The way in which basal tractions, associated with mantle convection, couples with the lithosphere is a fundamental problem in geodynamics. A successful lithosphere-mantle coupling model for the Earth will satisfy observations of plate motions, intraplate stresses, and the plate boundary zone deformation. We solve the depth integrated three-dimensional force balance equations in a global finite element model that takes into account effects of both topography and shallow lithosphere structure as well as tractions originating from deeper mantle convection. The contribution from topography and lithosphere structure is estimated by calculating gravitational potential energy differences. The basal tractions are derived from a fully dynamic flow model with both radial and lateral viscosity variations. We simultaneously fit stresses and plate motions in order to delineate a best-fit lithosphere-mantle coupling model. We use both the World Stress Map and the Global Strain Rate Model to constrain the models. We find that a strongly coupled model with a stiff lithosphere and 3-4 orders of lateral viscosity variations in the lithosphere are best able to match the observational constraints. Our predicted deviatoric stresses, which are dominated by contribution from mantle tractions, range between 20-70 MPa. The best-fitting coupled models predict strain rates that are consistent with observations. That is, the intraplate areas are nearly rigid whereas plate boundaries and some other continental deformation zones display high strain rates. Comparison of mantle tractions and surface velocities indicate that in most areas tractions are driving, although in a few regions, including western North America, tractions are resistive. Citation: Ghosh, A., W. E. Holt, and L. M. Wen (2013), Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics.

Joint modeling of lithosphere and mantle dynamics: Evaluation of constraints from global tomography models

With the advances in technology, seismological theory, and data acquisition, a number of high-resolution seismic tomography models have been published. However, discrepancies between tomography models often arise from different theoretical treatments of seismic wave propagation, different inversion strategies, and different data sets. Using a fixed velocity-to-density scaling and a fixed radial viscosity profile, we compute global mantle flow models associated with the different tomography models and test the impact of these for explaining surface geophysical observations (geoid, dynamic topography, stress, and strain rates). We use the joint modeling of lithosphere and mantle dynamics approach of Ghosh and Holt (2012) to compute the full lithosphere stresses, except that we use HC for the mantle circulation model, which accounts for the primary flow-coupling features associated with density-driven mantle flow. Our results show that the seismic tomography models of S40RTS and SAW642AN provide a better match with surface observables on a global scale than other models tested. Both of these tomography models have important similarities, including upwellings located in Pacific, Eastern Africa, Iceland, and mid-ocean ridges in the Atlantic and Indian Ocean and downwelling flows mainly located beneath the Andes, the Middle East, and central and Southeast Asia.

Mantle dynamics and seismic tomography

Three-dimensional imaging of the Earth's interior, called seismic tomography, has achieved breakthrough advances in the last two decades, revealing fundamental geodynamical processes throughout the Earth's mantle and core. Convective circulation of the entire mantle is taking place, with subducted oceanic lithosphere sinking into the lower mantle, overcoming the resistance to penetration provided by the phase boundary near 650-km depth that separates the upper and lower mantle. The boundary layer at the base of the mantle has been revealed to have complex structure, involving local stratification, extensive structural anisotropy, and massive regions of partial melt. The Earth's high Rayleigh number convective regime now is recognized to be much more interesting and complex than suggested by textbook cartoons, and continued advances in seismic tomography, geodynamical modeling, and high-pressure–high-temperature mineral physics will be needed to fully quantify the complex dynamics of our planet's interior.

Investigation of Mantle Dynamics from Platinum Group Elements and Rhenium-Osmium Isotope Systematics of Mantle Xenoliths from Oahu, Hawaii

Intraplate volcanism that has created the Hawaiian-Emperor seamount chain is generally thought to be formed by a deep-seated mantle plume. While the idea of a Hawaiian plume has not met with substantial opposition, whether or not the Hawaiian plume shows any geochemical signal of receiving materials from the Earth’s Outer Core and how the plume may or may not be reacting with the overriding lithosphere remain debatable issues. In an effort to understand how the Hawaiian plume works I report on the first in-situ sulfides and bulk rock Platinum Group Element (PGE) concentrations, together with Os isotope ratios on well-characterized garnet pyroxenite xenoliths from the island of Oahu in Hawaii. The sulfides are Fe-Ni Monosulfide Solid Solution and show fractionated PGE patterns. Based on the major elements, Platinum Group Elements and experimental data I interpret the Hawaiian sulfides as an immiscible melt that separated from a melt similar to the Honolulu Volcanics (HV) alkali lavas at a pressure-temperature condition of 1530 ± 100OC and 3.1±0.6 GPa., i.e. near the base or slightly below the Pacific lithosphere. The 187Os/188Os ratios of the bulk rock vary from subchondritic to suprachondritic (0.123-0.164); and the 187Os/188Os ratio strongly correlates with major element, High Field Strength Element (HFSE), Rare Earth Element (REE) and PGE abundances. These correlations strongly suggest that PGE concentrations and Os isotope ratios reflect primary mantle processes. I interpret these correlations as the result of melt-mantle reaction at the base of the lithosphere: I suggest that the parental melt that crystallized the pyroxenites selectively picked up radiogenic Os from the grain boundary sulfides, while percolating through the Pacific lithosphere. Thus the sampled pyroxenites essentially represent crystallized melts from different stages of this melt-mantle reaction process at the base of the lithosphere. I further show that the relatively low Pt/Re ratios of the Hawaiian sulfides and the bulk rock pyroxenites suggest that, upon ageing, such pyroxenites plus their sulfides cannot generate the coupled 186Os- 187Os isotope enrichments observed in Hawaiian lavas. Therefore, recycling of mantle sulfides of pyroxenitic parentage is unlikely to explain the enriched Pt-Re-Os isotope systematics of plume-derived lavas.

Three-dimensional models of mantle convection : influence of phase transitions and temperature-dependent viscosity

Two of the most important questions in mantle dynamics are investigated in three separate studies: the influence of phase transitions (studies 1 and 2), and the influence of temperature-dependent viscosity (study 3).

(1) Numerical modeling of mantle convection in a three-dimensional spherical shell incorporating the two major mantle phase transitions reveals an inherently three-dimensional flow pattern characterized by accumulation of cold downwellings above the 670 km discontinuity, and cylindrical 'avalanches' of upper mantle material into the lower mantle. The exothermic phase transition at 400 km depth reduces the degree of layering. A region of strongly-depressed temperature occurs at the base of the mantle. The temperature field is strongly modulated by this partial layering, both locally and in globally-averaged diagnostics. Flow penetration is strongly wavelength-dependent, with easy penetration at long wavelengths but strong inhibition at short wavelengths. The amplitude of the geoid is not significantly affected.

(2) Using a simple criterion for the deflection of an upwelling or downwelling by an endothermic phase transition, the scaling of the critical phase buoyancy parameter with the important lengthscales is obtained. The derived trends match those observed in numerical simulations, i.e., deflection is enhanced by (a) shorter wavelengths, (b) narrower up/downwellings (c) internal heating and (d) narrower phase loops.

(3) A systematic investigation into the effects of temperature-dependent viscosity on mantle convection has been performed in three-dimensional Cartesian geometry, with a factor of 1000-2500 viscosity variation, and Rayleigh numbers of 10^5-10^7. Enormous differences in model behavior are found, depending on the details of rheology, heating mode, compressibility and boundary conditions. Stress-free boundaries, compressibility, and temperature-dependent viscosity all favor long-wavelength flows, even in internally heated cases. However, small cells are obtained with some parameter combinations. Downwelling plumes and upwelling sheets are possible when viscosity is dependent solely on temperature. Viscous dissipation becomes important with temperature-dependent viscosity.

The sensitivity of mantle flow and structure to these various complexities illustrates the importance of performing mantle convection calculations with rheological and thermodynamic properties matching as closely as possible those of the Earth.

Deciphering petrological signatures of reactive melt stagnation and cooling in the oceanic mantle underneath ultraslow-spreading ridges

The global mid-ocean ridge system creates oceanic crust and lithosphere that covers more than two-thirds of the Earth. Basalts are volumetrically the most important rock type sampled at mid-ocean ridges. For this reason, our present understanding of upper mantle dynamics and the chemical evolution of the earth is strongly influenced by the study of mid-ocean ridge basalts (MORB). However, MORB are aggregates of polybarically generated small melt increments that can undergo a variety of physical and chemical processes during their ascent and consequently affect their derivative geochemical composition. Therefore, MORB do not represent “direct” windows to the underlying upper mantle. Abyssal peridotites, upper mantle rocks recovered from the ocean floor, are the residual complement to MORB melting and provide essential information on melt extraction from the upper mantle. In this study, abyssal peridotites are examined to address these overarching questions posed by previous studies of MORB: How are basaltic melts formed in the mantle, how are they extracted from the mantle and what physical and chemical processes control mantle melting? The number of studies on abyssal peridotites is small compared to those on basalts, in part because seafloor exposures of abyssal peridotites are relatively rare. For this reason, abyssal peridotite characteristics need to be considered in the context of subaerially exposed peridotites associated with ophiolites, orogenic peridotite bodies and basalt-hosted xenoliths. However, orogenic peridotite bodies are mainly associated with passive continental margins, most ophiolites are formed in supra-subduction zone settings, and peridotite xenoliths are often contaminated by their host magma. Therefore, studies of abyssal peridotites are essential to understanding the primary characteristics of the oceanic upper mantle free from the influence of continental rifting, subduction and tectonic emplacement processes. Nevertheless, numerous processes such as melt stagnation and cooling-induced, inter-mineral exchange can affect residual abyssal peridotite compositions after the cessation of melting. The aim of this study is to address these post-melting modifications of abyssal peridotites from a petrological-geochemical perspective. The samples in this study were dredged along the axis of the ultraslow-spreading Gakkel Ridge in the Arctic Ocean within the “Sparsely Magmatic Zone”, a 100 km ridge section where only mantle rocks are exposed. During two expeditions (ARK XVII-2 in 2001 and ARK XX-2 in 2004), exceptionally fresh peridotites were recovered. The boulders and cobbles collected cover a range of mantle rock compositions, with most characterized as plagioclase-free spinel peridotites or plagioclase- spinel peridotites. This thesis investigates melt stagnation and cooling processes in the upper mantle and is divided into two parts. The first part focuses on processes in the stability field of spinel peridotites (>10 kb) such as melt refertilization and cooling related trace element exchange, while the second part investigates processes in the stability field of plagioclase peridotites (< 10 kb) such as reactive melt migration and melt stagnation. The dissertation chapters are organized to follow the theoretical ascent of a mantle parcel upwelling beneath the location where the samples were collected.

Anisotropy model for mantle convection

The paper presents a new theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including folding and kinking in multi-layered visco-elastic rock. The orientation of slip planes in the context of crystallographic slip is determined by the normal vector, the so-called director of these surfaces. The model is applied to simulate anisotropic natural mantle convection. We compare the evolution of the director and approximately steady states of isotropic and anisotropic convection. The isotropic case has a simple steady state solution, whereas the orthotropic convection model produces a continuously evolving patterning in tile core of the convection cell which makes only a near-steady condition possible, in which the thermal boundary layer appears to be well aligned with the flow and hence as observed in seismic tomomgraphy strong anistropic.

利用宽频带台阵资料研究华北地区地壳及上地幔间断面结构

Intense tectonic renovation has occurred in the eastern continent of china since Mesozoic, as evidenced by the high heat flow, widespread magma extrusion and volcanic activities, and development of large sedimentary basins. To explain the cause and mechanism for the tectonic process in this period, some researchers have put forward various models, such as mantle plume, subduction of the Pacific slab, Yangtze Block-North China Block collision, etc. Their seismological evidence, however, is still scarce..During the period from 2000 to 2003, large temporary seismic arrays were established in North China by the Institute of the Geology and Geophysics, Chinese Academy of Sciences. Total 129 portable seismic stations were linearly emplaced across the western and eastern boundaries of the Bohai Bay Basin, and accumulated a large amount of high-quality data. Moreover, abundant data were also collected at the capital digital seismic network established in the ninth five-year period of national economic and social development. These provide an unprecedented opportunity for us to study the deep structure and associated geodynamic mechanism of lithospheric processes in North China using seismological techniques.Seismology is a kind of observation-based science. The development of seismic observations greatly promotes the improvement of seismologic theory and methodology. At the beginning of this thesis, I review the history of seismic observation progress, and present some routine processing techniques used in the array seismology. I also introduce two popular seismic imaging methods (receiver function method and seismic tomography).Receiver function method has been widely used to study the crustal and upper mantle structures, and many relevant research results have been published. In this thesis I elaborate the theory of this method, including the basic concept of receiver functions and the methodology for data pre-processing, stacking and migration. I also address some problems often encountered in practical applications of receiver function imaging.By using the teleseismic data collected at the temporary seismic arrays in North China, in particular, the traveltime information of P-to-S conversion and multiple reverberations of the Moho discontinuity, I obtain the distributions of the crustal thickness and the poisson ratio at the northwest boundary area of the Bohai Bay Basin and discuss the geological implications of the results.Through detailed intestigations on the crustal structural feature around the middle part of the Tanlu fault, considerable disparity in poisson ratios is found in the western and eastern sides of the Tanlu fault. Moreover, an obvious Moho offset is coincidently observed at the same surface location. A reasonable density model for the Tanlu fault area is also derived by simulating the observed gravity variations. Both receiver function study and gravity anomaly modeling suggest that the crustal difference between the western and eastern sides of the Tanlu fault is mainly resulted from their different compositions.With common conversion point imaging of receiver functions, I estimate the depths of the upper and lower boundaries of the mantle transition zone, i.e., the 410 and 660 km discontinuities, beneath most part of the North China continent The thickness of the transition zone (TTZ) in the study area is calculated by subtracting the depth of .410 km discontinuity from that of the 660km discontinuity. The resultant TTZ is 10-15 km larger in the east than in the west of the study area. Phase transitions at the 410 km and the 660 km discontinuities are known to have different Clapeyron slopes. Therefore, the TTZ is sensitive to the temperature changes in the transition zone. Previous studies have shown that the TTZ would be smaller in the mantle plume areas and become larger when the remnants of subducted slabs are present The hypothesis of mantle plume cannot give a reasonable interpretation to the observed TTZ beneath North China, Instead, the receiver function imaging results favor a dynamic model that correlates the thermal structure of the mantle transition zone and associated upper mantle dynamics of North China to the Pacific plate subduction process.

Pb-, Sr-, and Nd-isotopic composition of basaltic rocks from the Lau Basin

New Pb, Sr, and Nd isotope data are presented for 64 samples from the six backarc sites drilled during Leg 135. Systematic changes in Pb and Sr compositions illustrate significant isotopic variations between and within sites as well as provide two key pieces of information. First, a recent influx of asthenosphere with Indian Ocean mantle affinities has occurred and has successfully displaced older "Pacific" asthenosphere from the mantle underlying the backarc region. Second, clear evidence exists for mixing between these two asthenospheric end-members and at least one "arc-like" component. The latter was not the same as most material currently erupting in the Tofua Arc, but it must have had a more radiogenic Pb-isotope signature, perhaps similar to rocks analyzed from the islands of Tafahi, and Niuatoputapu. A comparison between the isotopic variations and the tectonic setting of the drill sites reveals consistent and important information regarding the mantle dynamics beneath the evolving backarc basin. We propose a model in which the source of upwelling magmas changes from Pacific to Indian Ocean asthenosphere with the propagation of seafloor spreading, a model with important implications for the rate of mantle influx into this region. Although the chemistries of backarc magmas have been profoundly influenced by this process, an additional consequence is the advection of Indian Ocean asthenosphere into the sub-arc mantle source. The isotopic compositions of arc rocks from the vicinity have been reevaluated on the basis of the proposed mantle advection model. We suggest that the slab-derived flux of trace elements into the arc wedge has remained relatively uniform with time (i.e., ~40 Ma), so that the change in arc chemistry results from mantle source substitution, rather than from differences in the composition of the downgoing plate.

Geochemistry of highly depleted peridotites of the Mid-Atlantic Ridge

During ODP Leg 209, a magma-starved area of the Mid-Atlantic Ridge (MAR) was drilled in the vicinity of the Fifteen-Twenty Fracture Zone (FZ) that offsets one of the slowest portions of the spreading ridge. We present here the results of a bulk rock multi-elemental study of 27 peridotites drilled at Sites 1272 and 1274 (to the south and the north of the FZ, respectively). The peridotites comprise mainly of harzburgites with minor dunites. Clinopyroxene (Cpx), which is interstitial and interpreted as secondary, is observed in Site 1274 peridotites. Sites 1272 and 1274 peridotites have low Al2O3 contents (<1 anhydrous wt.%), high Mg# (>91.5), and bulk rock trace elements compositions mostly below 0.1X primitive mantle (PM). These peridotites, and in particular Site 1272 peridotites, represent the most depleted peridotites yet sampled at a slow spreading ridge. Their compositions indicate high degrees of partial melting and melt extraction. A single open-system melting event (melting plus percolation of melts produced within upwelling mantle) can explain their highly depleted yet linear chondrite-normalized REE patterns, characterized by a steady depletion from HREE to LREE. Late melt-rock reactions and precipitation of Cpx explains the slightly less depleted compositions of Site 1274 peridotites. Hence, the differences in composition between Sites 1272 and 1274 peridotites do not provide evidence for regional variations in the degrees of partial melting from the south to the north of the FZ. The occurrence of highly refractory peridotites in the Fifteen-Twenty area suggests we sampled a more actively convecting mantle than generally supposed below slow spreading centers.