Evolution of the long-wavelength, subduction-driven topography of South America since 150 Ma

NASA Astrophysics Data System (ADS)

Flament, N. E.; Gurnis, M.; Williams, S.; Bower, D. J.; Seton, M.; Müller, D.

2014-12-01

Subduction to the west of South America spans 6000 km along strike and has been active for over 250 Myr. The influence of the history of subduction on the geodynamics of South America has been profound, driving mountain building and arc volcanism in the Andean Cordillera. Here, we investigate the long-wavelength changes in the topography of South America associated with subduction and plate motion and their interplay with the lithospheric deformation associated with the opening of the South Atlantic. We pay particular attention to the topographic expression of flat-lying subduction zones. We develop time-dependent geodynamic models of mantle flow and lithosphere deformation to investigate the evolution of South American dynamic and total topography since the late Jurassic (150 Ma). Our models are semi-empirical because the computational cost of fully dynamic, evolutionary models is still prohibitive. We impose the kinematics of global plate reconstructions with deforming continents in forward global mantle convection models with compositionally distinct crust and continental lithosphere embedded within the thermal lithosphere. The shallow thermal structure of subducting slabs is imposed, allowing us to investigate the evolution of dynamic topography around flat slab segments in time-dependent models. Multiple cases are used to investigate how the evolution of South American dynamic topography is influenced by mantle viscosity, the kinematics of the opening of the South Atlantic and alternative scenarios for recent and past flat-slab subduction. We predict that the migration of South America over sinking oceanic lithosphere resulted in continental tilt to the west until ~ 45 Ma, inverting to an eastward tilt thereafter. This first-order result is consistent with the reversal of the drainage of the Amazon River system. We investigate which scenarios of flat-slab subduction since the Eocene are compatible with geological constraints on the evolution of the Solimoes Basin, the Chaco Basin, the Sierras Pampeanas and the Central Patagonian Basin. To broadly constrain mantle viscosity, we compare models to the total subsidence inferred from well data offshore Argentina and Brazil, and to mantle tomography, since the initial and boundary conditions are based on independent plate reconstructions.

Insights from 3D numerical simulations on the dynamics of the India-Asia collision zone

NASA Astrophysics Data System (ADS)

Pusok, A. E.; Kaus, B.; Popov, A.

2013-12-01

The dynamics of the India-Asia collision zone remains one of the most remarkable topics of the current research interest: the transition from subduction to collision and uplift, followed by the rise of the abnormally thick Tibetan plateau, and the deformation at its Eastern and Western syntaxes, are processes still not fully understood. Models that have addressed this topic include wholescale underthrusting of Indian lithospheric mantle under Tibet, distributed homogeneous shortening or the thin-sheet model, slip-line field model for lateral extrusion or lower crustal flow models for the exhumation of the Himalayan units and lateral spreading of the Tibetan plateau. Of these, the thin-sheet model has successfully illustrated some of the basic physics of continental collision and has the advantage of a 3D model being reduced to 2D, but one of its major shortcomings is that it cannot simultaneously represent channel flow and gravitational collapse of the mantle lithosphere, since these mechanisms require the lithosphere to interact with the underlying mantle, or to have a vertically non-homogeneous rheology. As a consequence, 3D models are emerging as powerful tools to understand the dynamics of coupled systems. However, because of yet recent developments and various complexities, the current 3D models simulating the dynamics of continent collision zones have relied on certain explicit assumptions, such as replacing part of the asthenosphere with various types of boundary conditions that mimic the effect of mantle flow, in order to focus on the lithospheric/crustal deformation. Here, we employ the parallel 3D code LaMEM (Lithosphere and Mantle Evolution Model), with a finite difference staggered grid solver, which is capable of simulating lithospheric deformation while simultaneously taking mantle flow and a free surface into account. We present qualitative results on lithospheric and upper-mantle scale simulations in which the Indian lithosphere is subducted and/or indented into Asia. We investigate the way deep processes affect continental tectonics at convergent margins, addressing the role the continent subduction and indentation plays on the development of continental tectonics during convergence and we discuss the implications these offer for the Asian tectonics. Acknowledgements: Funding was provided by the European Research Council under the European Community's Seventh Framework Program (FP7/2007-2013) / ERC Grant agreement #258830. Numerical computations have been performed on MOGON (ZDV Mainz computing center) and JUQUEEN (Jülich high-performance computing center).

A Hybrid Approach to Data Assimilation for Reconstructing the Evolution of Mantle Dynamics

NASA Astrophysics Data System (ADS)

Zhou, Quan; Liu, Lijun

2017-11-01

Quantifying past mantle dynamic processes represents a major challenge in understanding the temporal evolution of the solid earth. Mantle convection modeling with data assimilation is one of the most powerful tools to investigate the dynamics of plate subduction and mantle convection. Although various data assimilation methods, both forward and inverse, have been created, these methods all have limitations in their capabilities to represent the real earth. Pure forward models tend to miss important mantle structures due to the incorrect initial condition and thus may lead to incorrect mantle evolution. In contrast, pure tomography-based models cannot effectively resolve the fine slab structure and would fail to predict important subduction-zone dynamic processes. Here we propose a hybrid data assimilation approach that combines the unique power of the sequential and adjoint algorithms, which can properly capture the detailed evolution of the downgoing slab and the tomographically constrained mantle structures, respectively. We apply this new method to reconstructing mantle dynamics below the western U.S. while considering large lateral viscosity variations. By comparing this result with those from several existing data assimilation methods, we demonstrate that the hybrid modeling approach recovers the realistic 4-D mantle dynamics the best.

A Hybrid Forward-Adjoint Data Assimilation Method for Reconstructing the Temporal Evolution of Mantle Dynamics

NASA Astrophysics Data System (ADS)

Zhou, Q.; Liu, L.

2017-12-01

Quantifying past mantle dynamic processes represents a major challenge in understanding the temporal evolution of the solid earth. Mantle convection modeling with data assimilation is one of the most powerful tools to investigate the dynamics of plate subduction and mantle convection. Although various data assimilation methods, both forward and inverse, have been created, these methods all have limitations in their capabilities to represent the real earth. Pure forward models tend to miss important mantle structures due to the incorrect initial condition and thus may lead to incorrect mantle evolution. In contrast, pure tomography-based models cannot effectively resolve the fine slab structure and would fail to predict important subduction-zone dynamic processes. Here we propose a hybrid data assimilation method that combines the unique power of the sequential and adjoint algorithms, which can properly capture the detailed evolution of the downgoing slab and the tomographically constrained mantle structures, respectively. We apply this new method to reconstructing mantle dynamics below the western U.S. while considering large lateral viscosity variations. By comparing this result with those from several existing data assimilation methods, we demonstrate that the hybrid modeling approach recovers the realistic 4-D mantle dynamics to the best.

Effects of grain size evolution on mantle dynamics

NASA Astrophysics Data System (ADS)

Schulz, Falko; Tosi, Nicola; Plesa, Ana-Catalina; Breuer, Doris

2016-04-01

The rheology of planetary mantle materials is strongly dependent on temperature, pressure, strain-rate, and grain size. In particular, the rheology of olivine, the most abundant mineral of the Earth's upper mantle, has been extensively studied in the laboratory (e.g., Karato and Wu, 1993; Hirth and Kohlstedt, 2003). Two main mechanisms control olivine's deformation: dislocation and diffusion creep. While the former implies a power-law dependence of the viscosity on the strain-rate that leads to a non-Newtonian behaviour, the latter is sensitively dependent on the grain size. The dynamics of planetary interiors is locally controlled by the deformation mechanism that delivers the lowest viscosity. Models of the dynamics and evolution of planetary mantles should thus be capable to self-consistently distinguish which of the two mechanisms dominates at given conditions of temperature, pressure, strain-rate and grain size. As the grain size can affect the viscosity associated with diffusion creep by several orders of magnitude, it can strongly influence the dominant deformation mechanism. The vast majority of numerical, global-scale models of mantle convection, however, are based on the use of a linear diffusion-creep rheology with constant grain-size. Nevertheless, in recent studies, a new equation has been proposed to properly model the time-dependent evolution of the grain size (Austin and Evens, 2007; Rozel et al., 2010). We implemented this equation in our mantle convection code Gaia (Hüttig et al., 2013). In the framework of simple models of stagnant lid convection, we compared simulations based on the fully time-dependent equation of grain-size evolution with simulations based on its steady-state version. In addition, we tested a number of different parameters in order to identify those that affects the grain size to the first order and, in turn, control the conditions at which mantle deformation is dominated by diffusion or dislocation creep. References Austin, N. J. and Evans, B. (2007). Geology, 35(4):343. Hirth, G. and Kohlstedt, D. (2003). Geophysical Monograph Series, page 83105. Hüttig, C., Tosi, N., and Moore, W. B. (2013). Physics of the Earth and Planetary Interiors, 220:11-18. Karato, S.-i. and Wu, P. (1993). Science, 260(5109):771778. Rozel, A., Ricard, Y., and Bercovici, D. (2010). Geophysical Journal International, 184(2):719728.

Dynamical effects on the core-mantle boundary from depth-dependent thermodynamical properties of the lower mantle

NASA Technical Reports Server (NTRS)

Zhang, Shuxia; Yuen, David A.

1988-01-01

A common assumption in modeling dynamical processes in the lower mantle is that both the thermal expansivity and thermal conductivity are reasonably constant. Recent work from seismic equation of state leads to substantially higher values for the thermal conductivity and much lower thermal expansivity values in the deep mantle. The dynamical consequences of incorporating depth-dependent thermodynamic properties on the thermal-mechanical state of the lower mantle are examined with the spherical-shell mean-field equations. It is found that the thermal structure of the seismically resolved anomalous zone at the base of the mantle is strongly influenced by these variable properties and, in particular, that the convective distortion of the core-mantle boundary (CMB) is reduced with the decreasing thermal expansivity. Such a reduction of the dynamically induced topography from pure thermal convection would suggest that some other dynamical mechanism must be operating at the CMB.

Hydrogen self-diffusion in single crystal olivine and electrical conductivity of the Earth’s mantle

DOE Office of Scientific and Technical Information (OSTI.GOV)

Novella, Davide; Jacobsen, Benjamin; Weber, Peter K.

Nominally anhydrous minerals formed deep in the mantle and transported to the Earth’s surface contain tens to hundreds of ppm wt H 2O, providing evidence for the presence of dissolved water in the Earth’s interior. Even at these low concentrations, H 2O greatly affects the physico-chemical properties of mantle materials, governing planetary dynamics and evolution. The diffusion of hydrogen (H) controls the transport of H 2O in the Earth’s upper mantle, but is not fully understood for olivine ((Mg, Fe) 2SiO 4) the most abundant mineral in this region. Here we present new hydrogen self-diffusion coefficients in natural olivine singlemore » crystals that were determined at upper mantle conditions (2 GPa and 750–900 °C). Hydrogen self-diffusion is highly anisotropic, with values at 900 °C of 10 -10.9, 10 -12.8 and 10 -11.9 m 2/s along [100], [010] and [001] directions, respectively. Combined with the Nernst-Einstein relation, these diffusion results constrain the contribution of H to the electrical conductivity of olivine to be σH = 10 2.12S/m·C H2O·exp -187kJ/mol/(RT). Comparisons between the model presented in this study and magnetotelluric measurements suggest that plausible H 2O concentrations in the upper mantle (≤250 ppm wt) can account for high electrical conductivity values (10 -2–10 -1 S/m) observed in the asthenosphere.« less

Hydrogen self-diffusion in single crystal olivine and electrical conductivity of the Earth’s mantle

DOE PAGES

Novella, Davide; Jacobsen, Benjamin; Weber, Peter K.; ...

2017-07-13

Nominally anhydrous minerals formed deep in the mantle and transported to the Earth’s surface contain tens to hundreds of ppm wt H 2O, providing evidence for the presence of dissolved water in the Earth’s interior. Even at these low concentrations, H 2O greatly affects the physico-chemical properties of mantle materials, governing planetary dynamics and evolution. The diffusion of hydrogen (H) controls the transport of H 2O in the Earth’s upper mantle, but is not fully understood for olivine ((Mg, Fe) 2SiO 4) the most abundant mineral in this region. Here we present new hydrogen self-diffusion coefficients in natural olivine singlemore » crystals that were determined at upper mantle conditions (2 GPa and 750–900 °C). Hydrogen self-diffusion is highly anisotropic, with values at 900 °C of 10 -10.9, 10 -12.8 and 10 -11.9 m 2/s along [100], [010] and [001] directions, respectively. Combined with the Nernst-Einstein relation, these diffusion results constrain the contribution of H to the electrical conductivity of olivine to be σH = 10 2.12S/m·C H2O·exp -187kJ/mol/(RT). Comparisons between the model presented in this study and magnetotelluric measurements suggest that plausible H 2O concentrations in the upper mantle (≤250 ppm wt) can account for high electrical conductivity values (10 -2–10 -1 S/m) observed in the asthenosphere.« less

Hydrogen self-diffusion in single crystal olivine and electrical conductivity of the Earth's mantle.

PubMed

Novella, Davide; Jacobsen, Benjamin; Weber, Peter K; Tyburczy, James A; Ryerson, Frederick J; Du Frane, Wyatt L

2017-07-13

Nominally anhydrous minerals formed deep in the mantle and transported to the Earth's surface contain tens to hundreds of ppm wt H 2 O, providing evidence for the presence of dissolved water in the Earth's interior. Even at these low concentrations, H 2 O greatly affects the physico-chemical properties of mantle materials, governing planetary dynamics and evolution. The diffusion of hydrogen (H) controls the transport of H 2 O in the Earth's upper mantle, but is not fully understood for olivine ((Mg, Fe) 2 SiO 4 ) the most abundant mineral in this region. Here we present new hydrogen self-diffusion coefficients in natural olivine single crystals that were determined at upper mantle conditions (2 GPa and 750-900 °C). Hydrogen self-diffusion is highly anisotropic, with values at 900 °C of 10 -10.9 , 10 -12.8 and 10 -11.9 m 2 /s along [100], [010] and [001] directions, respectively. Combined with the Nernst-Einstein relation, these diffusion results constrain the contribution of H to the electrical conductivity of olivine to be σ H  = 10 2.12 S/m·C H2O ·exp -187kJ/mol/(RT) . Comparisons between the model presented in this study and magnetotelluric measurements suggest that plausible H 2 O concentrations in the upper mantle (≤250 ppm wt) can account for high electrical conductivity values (10 -2 -10 -1  S/m) observed in the asthenosphere.

Predicting and testing continental vertical motion histories since the Paleozoic

NASA Astrophysics Data System (ADS)

Zhang, Nan; Zhong, Shijie; Flowers, Rebecca M.

2012-02-01

Dynamic topography at the Earth's surface caused by mantle convection can affect a range of geophysical and geological observations including bathymetry, sea-level change, continental flooding, sedimentation and erosion. These observations provide important constraints on and test of mantle dynamic models. Based on global mantle convection models coupled with the surface plate motion history, we compute dynamic topography and its history for the last 400 Ma associated with Pangea assembly and breakup, with particular focus on cratonic regions. We propose that burial-unroofing histories of cratons inferred from thermochronology data can be used as a new diagnostic to test dynamic topography and mantle dynamic models. Our models show that there are currently two broad dynamic topography highs in the Pacific and Africa for the present-day Earth that are associated with the broad, warm structures (i.e., superplumes) in the deep mantle, consistent with previous proposals of dynamical support for the Pacific and African superswells. Our models reveal that Pangea assembly and breakup, by affecting subduction and mantle upwelling processes, have significant effects on continental vertical motions. Our models predict that the Slave craton in North America subsides before Pangea assembly at 330 Ma but uplifts significantly from 330 Ma to 240 Ma in response to pre-Pangea subduction and post-assembly mantle warming. The Kaapvaal craton of Africa is predicted to undergo uplift from ~180 Ma to 90 Ma after Pangea breakup, but its dynamic topography remains stable for the last 90 Ma. The predicted histories of elevation change for the Slave and Kaapvaal cratons compare well with the burial-unroofing histories inferred from thermochronology studies, thus supporting our dynamic models including the development of the African superplume mantle structure. The vertical motion histories for other cratons can provide further tests of and constraints on our mantle dynamic models.

Predicting and testing continental vertical motion histories since the Paleozoic

NASA Astrophysics Data System (ADS)

Zhang, N.; Zhong, S.; Flowers, R. M.

2011-12-01

Dynamic topography at the Earth's surface caused by mantle convection can affect a range of geophysical and geological observations including bathymetry, sea-level change, continental flooding, sedimentation and erosion. These observations provide important constraints on and test of mantle dynamic models. Based on global mantle convection models coupled with the surface plate motion history, we compute dynamic topography and its history for the last 400 Ma associated with Pangea assembly and breakup, with particular focus on continental cratonic regions. We propose that burial-unroofing histories of continental cratons inferred from thermochronology data can be used as a new diagnostic to test dynamic topography and mantle dynamic models. Our models show that there are currently two broad dynamic topography highs in the Pacific and Africa for the present-day Earth that are associated with the broad, warm structures (i.e., superplumes) in the deep mantle, consistent with previous proposals of dynamical support for the Pacific and African superswells. Our models reveal that Pangea assembly and breakup, by affecting subduction and mantle upwelling processes, have significant effects on continental vertical motions. Our models predict that the Slave craton in North America subsides before Pangea assembly at 330 Ma but uplifts significantly from 330 Ma to 240 Ma in response to pre-Pangea subduction and post-assembly mantle warming. The Kaapvaal craton of Africa is predicted to undergo uplift from ~180 Ma to 90 Ma after Pangea breakup, but its dynamic topography remains stable for the last 90 Ma. The predicted histories of elevation change for the Slave and Kaapvaal cratons compare well with the burial-unroofing histories inferred from thermochronology studies, thus supporting our dynamic models including the development of the African superplume mantle structure. The vertical motion histories for other cratons can provide further tests and constraints on our mantle dynamic models.

Orogenic delamination - dynamics, effects, and geological expression

NASA Astrophysics Data System (ADS)

Ueda, Kosuke; Gerya, Taras

2010-05-01

Unbundling of continental lithosphere and removal of its mantle portion have been described by two mutually rather exclusive models, convective thinning and integral delamination. Either disburdens the remaining lithosphere, weakens the remainder, and causes uplift and extension. Increased heat flux is likely to promote high-degree crustal melting, and has been viewed as a source for voluminous granitic intrusions in late or collapsing orogenic settings. Collapse may be driven by any of gravitational potential differences from orogen to foreland, by stress inversion in the unburdened domain, or by suction of a retreating trench. In this study, we investigate prerequisites, mechanism, and development paths for orogeny-related mantle lithosphere removal. Our experiments numerically reproduce delamination which self-consistently results from the dynamics of a decoupling collision zone. In particular, it succeeds without a seed facilitating initial separation of layers. External shortening of a continent - ocean - continent assembly, such as to initiate oceanic subduction, is lifted before the whole oceanic part is consumed, leaving slab pull to govern further convergence. Once buoyant continental crust enters, the collision zone locks, and convergence diminishes. Under favourable conditions, delamination then initiates close to the edge of the mantle wedge and at deep crustal levels. While it initially separates upper crust from lower crust according to the weakness minimum in the lithospheric strength profile, the lower crust is eventually also delaminated from the subducting lithospheric mantle, owing to buoyancy differences. The level of delamination within the lithosphere seems thus first rheology-controlled, then density-controlled. Subduction-coupled delamination is contingent on retreat and decoupling of the subducting slab, which in turn is dependent on effective rheological weakening of the plate contact. Weakening is a function of shear-heating and hereby of collision rate, melting and hydration, the latter two incorporating the effects of sediment subduction and phase changes. The drag available for slab retreat scales with the age of the descending oceanic lithosphere; integrated strength of the lithosphere and activation volume for mantle creep additionally control angle and depth of the descent. Fully developed delamination is observed from between 10 to 15 Ma after collision ceases, with following trenchward migration of the delamination front. Consequently, the main maximum extension migrates, while local, partly intermittent compression can be observed on smaller scale. Across the orogen, extension thus has a strongly diachronous main component. We track common surface observables such as heat flow, partially melted rocks (domal migmatites), and predicted geo-/thermochronological ages over the evolving plate boundary. Geochemical projections of our observations confirm potential contamination of reservoirs - although the net delamination level follows the Moho, some crustal remnants along the old slab still sink through the 660-discontinuity. On the other hand, the base of the delaminated domain is not as plain a contact as in concept. Where the contact of asthenosphere with delaminated crust is the location of high-degree melting, also traces of original lithospheric mantle can be entangled. Our results do not fully support the conceptual distinction between convective thinning and blockwise delamination. While the foundering portion initially retains a fairly coherent, slab-like perimeter, the actual separation of layers in a limited process-zone occurs in smaller -scale eddies. Also, convection of the whole uprising asthenosphere wedge is dynamically not discernible from the latter and crucial for the removal of lithospheric mantle. The removed lithosphere does initially not convect, but subsequently shows an increasing tendency to drip down. In the presented case, extension in the axial zone of the orogen is not (only) caused by unsupported gravitational potential of the core domain itself, but actively driven by slab retreat with a shallow mantle dynamic contribution.

An Interconnected Network of Core-Forming Melts Produced by Shear Deformation

NASA Technical Reports Server (NTRS)

Bruhn, D.; Groebner, N.; Kohlstedt, D. L.

2000-01-01

The formation mechanism of terrestrial planetary is still poorly understood, and has been the subject of numerous experimental studies. Several mechanisms have been proposed by which metal-mainly iron with some nickel-could have been extracted from a silicate mantle to form the core. Most recent models involve gravitational sinking of molten metal or metal sulphide through a partially or fully molten mantle that is often referred to as a'magma ocean. Alternative models invoke percolation of molten metal along an interconnected network (that is, porous flow) through a solid silicate matrix. But experimental studies performed at high pressures have shown that, under hydrostatic conditions, these melts do not form an interconnected network, leading to the widespread assumption that formation of metallic cores requires a magma ocean. In contrast, here we present experiments which demonstrate that shear deformation to large strains can interconnect a significant fraction of initially isolated pockets of metal and metal sulphide melts in a solid matrix of polycrystalline olivine. Therefore, in a dynamic (nonhydrostatic) environment, percolation remains a viable mechanism for the segregation and migration of core-forming melts in a solid silicate mantle.

Geochemical Constraints on the Size of the Moon-Forming Giant Impact

NASA Astrophysics Data System (ADS)

Piet, Hélène; Badro, James; Gillet, Philippe

2017-12-01

Recent models involving the Moon-forming giant impact hypothesis have managed to reproduce the striking isotopic similarity between the two bodies, albeit using two extreme models: one involves a high-energy small impactor that makes the Moon out of Earth's proto-mantle; the other supposes a gigantic collision between two half-Earths creating the Earth-Moon system from both bodies. Here we modeled the geochemical influence of the giant impact on Earth's mantle and found that impactors larger than 15% of Earth mass result in mantles always violating the present-day concentrations of four refractory moderately siderophile trace elements (Ni, Co, Cr, and V). In the aftermath of the impact, our models cannot further discriminate between a fully and a partially molten bulk silicate Earth. Then, the preservation of primordial geochemical reservoirs predating the Moon remains the sole argument against a fully molten mantle after the Moon-forming impact.

Variation of the subsidence parameters, effective thermal conductivity, and mantle dynamics

NASA Astrophysics Data System (ADS)

Adam, C.; King, S. D.; Vidal, V.; Rabinowicz, M.; Jalobeanu, A.; Yoshida, M.

2015-09-01

The subsidence of young seafloor is generally considered to be a passive phenomenon related to the conductive cooling of the lithosphere after its creation at mid-oceanic ridges. Recent alternative theories suggest that the mantle dynamics plays an important role in the structure and depth of the oceanic lithosphere. However, the link between mantle dynamics and seafloor subsidence has still to be quantitatively assessed. Here we provide a statistical study of the subsidence parameters (subsidence rate and ridge depth) for all the oceans. These parameters are retrieved through two independent methods, the positive outliers method, a classical method used in signal processing, and through the MiFil method. From the subsidence rate, we compute the effective thermal conductivity, keff, which ranges between 1 and 7 W m-1 K-1. We also model the mantle flow pattern from the S40RTS tomography model. The density anomalies derived from S40RTS are used to compute the instantaneous flow in a global 3D spherical geometry. We show that departures from the keff = 3 Wm-1K-1 standard value are systematically related to mantle processes and not to lithospheric structure. Regions characterized by keff > 3 Wm-1K-1 are associated with mantle uplifts (mantle plumes or other local anomalies). Regions characterized by keff < 3 Wm-1K-1 are related to large-scale mantle downwellings such as the Australia-Antarctic Discordance (AAD) or the return flow from the South Pacific Superswell to the East Pacific Rise. This demonstrates that mantle dynamics plays a major role in the shaping of the oceanic seafloor. In particular, the parameters generally considered to quantify the lithosphere structure, such as the thermal conductivity, are not only representative of this structure but also incorporate signals from the mantle convection occurring beneath the lithosphere. The dynamic topography computed from the S40RTS tomography model reproduces the subsidence pattern observed in the bathymetry. Overall we find a good correlation between the subsidence parameters derived from the bathymetry and the dynamic topography. This demonstrates that these parameters are strongly dependent on mantle dynamics.

An interconnected network of core-forming melts produced by shear deformation

PubMed

Bruhn; Groebner; Kohlstedt

2000-02-24

The formation mechanism of terrestrial planetary cores is still poorly understood, and has been the subject of numerous experimental studies. Several mechanisms have been proposed by which metal--mainly iron with some nickel--could have been extracted from a silicate mantle to form the core. Most recent models involve gravitational sinking of molten metal or metal sulphide through a partially or fully molten mantle that is often referred to as a 'magma ocean'. Alternative models invoke percolation of molten metal along an interconnected network (that is, porous flow) through a solid silicate matrix. But experimental studies performed at high pressures have shown that, under hydrostatic conditions, these melts do not form an interconnected network, leading to the widespread assumption that formation of metallic cores requires a magma ocean. In contrast, here we present experiments which demonstrate that shear deformation to large strains can interconnect a significant fraction of initially isolated pockets of metal and metal sulphide melts in a solid matrix of polycrystalline olivine. Therefore, in a dynamic (non-hydrostatic) environment, percolation remains a viable mechanism for the segregation and migration of core-forming melts in a solid silicate mantle.

Modeling Earth's surface topography: decomposition of the static and dynamic components

NASA Astrophysics Data System (ADS)

Guerri, M.; Cammarano, F.; Tackley, P. J.

2017-12-01

Isolating the portion of topography supported by mantle convection, the so-called dynamic topography, would give us precious information on vigor and style of the convection itself. Contrasting results on the estimate of dynamic topography motivate us to analyse the sources of uncertainties affecting its modeling. We obtain models of mantle and crust density, leveraging on seismic and mineral physics constraints. We use the models to compute isostatic topography and residual topography maps. Estimates of dynamic topography and associated synthetic geoid are obtained by instantaneous mantle flow modeling. We test various viscosity profiles and 3D viscosity distributions accounting for inferred lateral variations in temperature. We find that the patterns of residual and dynamic topography are robust, with an average correlation coefficient of 0.74 and 0.71, respectively. The amplitudes are however poorly constrained. For the static component, the considered lithospheric mantle density models result in topographies that differ, on average, 720 m, with peaks reaching 1.7 km. The crustal density models produce variations in isostatic topography averaging 350 m, with peaks of 1 km. For the dynamic component, we obtain peak-to-peak topography amplitude exceeding 3 km for all the tested mantle density and viscosity models. Such values of dynamic topography produce geoid undulations that are not in agreement with observations. Assuming chemical heterogeneities in the lower mantle, in correspondence with the LLSVPs (Large Low Shear wave Velocity Provinces), helps to decrease the amplitudes of dynamic topography and geoid, but reduces the correlation between synthetic and observed geoid. The correlation coefficients between the residual and dynamic topography maps is always less than 0.55. In general, our results indicate that, i) current knowledge of crust density, mantle density and mantle viscosity is still limited, ii) it is important to account for all the various sources of uncertainties when computing static and dynamic topography. In conclusion, a multidisciplinary approach, which involves multiple geophysics observations and constraints from mineral physics, is necessary for obtaining robust density models and, consequently, for properly estimating the dynamic topography.

Coldspots and hotspots - Global tectonics and mantle dynamics of Venus

NASA Technical Reports Server (NTRS)

Bindschadler, Duane L.; Schubert, Gerald; Kaula, William M.

1992-01-01

Based on geologic observations provided by Magellan's first cycle of data collection and recent models of mantle convection in spherical shells and crustal deformation, the major topographic and geologic features of Venus are incorporated into a model of global mantle dynamics. Consideration is given to volcanic rises, such as Beta Regio and Atla Regio, plateau-shaped highlands dominated by complex ridged terrain (e.g., Ovda Regio and Alpha Regio), and circular lowland regions, such as Atalanta Planitia. Each of these features is related to either mantle plumes (hotspots) or mantle downwellings (coldspots).

Geochemical cycling, mass balance, and Earth’s dynamic structure (Invited)

NASA Astrophysics Data System (ADS)

Allegre, C. J.

2009-12-01

The use of radiogenic tracers is now established as one of the standard methods (together with seismic tomography) to constrain Earth models. One of the advantages of radiogenic tracers is that they constrain not only present day structures, but also their historical evolution since the earliest time. A mandatory condition is to use not only one but all of the available isotopic systems, with their diverse cycling properties and different radioactive half-lives. Mass balances calculations are the basic tool which allows one to use quantitatively the various tracers. However, the inverse method approach, when applied to mass balance, is strongly non-linear, particularly for tracers with intermediate to short half lives. 142Nd, 82W, and 129Xe allows one to define a consistent model for the Present and past evolutions. Results point to a 3-layer model for the present mantle: an upper mantle with two layers separated by the 450 km discontinuity, and a lower mantle below the 670 km discontinuity. These layers all convect independently though none of them is fully isolated. They exchange matter and energy, not necessarily through the same process. For instance, most hotspots (though not all of them!) are boundary layer instabilities generated at the 670 km discontinuity, but there is important heat transfer from the lower mantle through megablobs. Recent results on the Hadean period are in agreement with core-mantle- atmosphere differentiation at 4444 My. The early crust was mostly anorthositic but has been recycled into the upper mantle within the first billion years. Formation of granitic continents started at 4300 My. 80% of their material was already present on the surface as acidic rocks by 2500 My. The reworking process involved in continental development increased continuously with time and today is the dominant process. Continents are now in a stage of steady state, with general addition of mantle and subducted continental material.

The importance of grain size to mantle dynamics and seismological observations

NASA Astrophysics Data System (ADS)

Gassmoeller, R.; Dannberg, J.; Eilon, Z.; Faul, U.; Moulik, P.; Myhill, R.

2017-12-01

Grain size plays a key role in controlling the mechanical properties of the Earth's mantle, affecting both long-timescale flow patterns and anelasticity on the timescales of seismic wave propagation. However, dynamic models of Earth's convecting mantle usually implement flow laws with constant grain size, stress-independent viscosity, and a limited treatment of changes in mineral assemblage. We study grain size evolution, its interplay with stress and strain rate in the convecting mantle, and its influence on seismic velocities and attenuation. Our geodynamic models include the simultaneous and competing effects of dynamic recrystallization resulting from dislocation creep, grain growth in multiphase assemblages, and recrystallization at phase transitions. They show that grain size evolution drastically affects the dynamics of mantle convection and the rheology of the mantle, leading to lateral viscosity variations of six orders of magnitude due to grain size alone, and controlling the shape of upwellings and downwellings. Using laboratory-derived scaling relationships, we convert model output to seismologically-observable parameters (velocity, attenuation) facilitating comparison to Earth structure. Reproducing the fundamental features of the Earth's attenuation profile requires reduced activation volume and relaxed shear moduli in the lower mantle compared to the upper mantle, in agreement with geodynamic constraints. Faster lower mantle grain growth yields best fit to seismic observations, consistent with our re-examination of high pressure grain growth parameters. We also show that ignoring grain size in interpretations of seismic anomalies may underestimate the Earth's true temperature variations.

Dynamics of Compressible Convection and Thermochemical Mantle Convection

NASA Astrophysics Data System (ADS)

Liu, Xi

The Earth's long-wavelength geoid anomalies have long been used to constrain the dynamics and viscosity structure of the mantle in an isochemical, whole-mantle convection model. However, there is strong evidence that the seismically observed large low shear velocity provinces (LLSVPs) in the lowermost mantle are chemically distinct and denser than the ambient mantle. In this thesis, I investigated how chemically distinct and dense piles influence the geoid. I formulated dynamically self-consistent 3D spherical convection models with realistic mantle viscosity structure which reproduce Earth's dominantly spherical harmonic degree-2 convection. The models revealed a compensation effect of the chemically dense LLSVPs. Next, I formulated instantaneous flow models based on seismic tomography to compute the geoid and constrain mantle viscosity assuming thermochemical convection with the compensation effect. Thermochemical models reconcile the geoid observations. The viscosity structure inverted for thermochemical models is nearly identical to that of whole-mantle models, and both prefer weak transition zone. Our results have implications for mineral physics, seismic tomographic studies, and mantle convection modelling. Another part of this thesis describes analyses of the influence of mantle compressibility on thermal convection in an isoviscous and compressible fluid with infinite Prandtl number. A new formulation of the propagator matrix method is implemented to compute the critical Rayleigh number and the corresponding eigenfunctions for compressible convection. Heat flux and thermal boundary layer properties are quantified in numerical models and scaling laws are developed.

Convective Differentiation of the Earth's Mantle

NASA Astrophysics Data System (ADS)

Hansen, U.; Schmalzl, J.; Stemmer, K.

2007-05-01

The differentiation of the Earth is likely to be influenced by convective motions within the early mantle. Double- diffusive convection (d.d.c), driven by thermally and compositionally induced density differences is considered as a vital mechanism behind the dynamic differentiation of the early mantle.. We demonstrate that d.d.c can lead to layer formation on a planetary scale in the diffusive regime where composition stabilizes the system whil heat provides the destabilizing force. Choosing initial conditions in which a stable compositional gradient overlies a hot reservoir we mimic the situation of a planet in a phase after core formation. Differently from earlier studies we fixed the temperature rather than the heat flux at the lower boundary, resembling a more realistic condition for the core-mantle boundary. We have carried out extended series of numerical experiments, ranging from 2D calculations in constant viscosity fluids to fully 3D experiments in spherical geometry with strongly temperature dependent viscosity. The buoyancy ratio R and the Lewis number Le are the important dynamical parameters. In all scenarios we could identify a parameter regime where the non-layered initial structure developed into a state consisting of several, mostly two layers. Initially plumes from the bottom boundary homogenize a first layer which subsequently thickens. The bottom layer heats up and then convection is initiated in the top layer. This creates dynamically (i.e. without jump in the material behavior) a stack of separately convecting layers. The bottom layer is significantly thicker than the top layer. Strongly temperature dependent viscosity leads to a more complex evolution The formation of the bottom layer is followed by the generation of several layers on top. Finally the uppermost layer starts to convect. In general, the multilayer structure collapses into a two layer system. We employed a numerical technique, allowing for a diffusion free treatment of the compositional field. In each case a similar evolution has been observed. This indicates that a temporary formation of layered structures in planetary interiors is a typical phenomenon. Moreover, in this scenario, plate tectonics appears only in later stages of the evolution.

Quantitative Restoration of the Evolution of Mantle Structures Using Data Assimilation

NASA Astrophysics Data System (ADS)

Ismail-Zadeh, A.; Schubert, G.; Tsepelev, I.

2008-12-01

Rapid progress in imaging deep Earth structures and in studies of physical and chemical properties of mantle rocks facilitates research in assimilation of data related to mantle dynamics. We present a quantitative approach to assimilation of geophysical and geodetic data, which allows for incorporating observations and unknown initial conditions for mantle temperature and flow into a three-dimensional dynamic model in order to determine the initial conditions in the geological past. Once the conditions are determined the evolution of mantle structures can be restore backward in time. We apply data assimilation techniques to model the evolution of mantle plumes and lithospheric slabs. We show that the geometry of the mantle structures changes with time diminishing the degree of surface curvature of the structures, because the heat conduction smoothes the complex thermal surfaces of mantle bodies with time. Present seismic tomography images of mantle structures do not allow definition of the sharp shapes of these structures. Assimilation of mantle temperature and flow to the geological past instead provides a quantitative tool to restore thermal shapes of prominent structures in the past from their diffusive shapes at present.

Quantifying mantle structure and dynamics using plume tracing in seismic tomography

NASA Astrophysics Data System (ADS)

O'Farrell, K. A.; Eakin, C. M.; Jackson, M. G.; Jones, T. D.; Lekic, V.; Lithgow-Bertelloni, C. R.

2017-12-01

Directly linking deep mantle processes with surface features and dynamics is a complex problem. Hotspot volcanism gives us surface observables of mantle signatures, but the depth and source of the mantle plumes feeding these hotspots are highly debated. To address these issues, it is necessary to consider the entire journey of a plume through the mantle. By analyzing the behavior of mantle plumes we can constrain the vigor of mantle convection, the net rotation of the mantle and the role of thermal versus chemical anomalies as well as the bulk physical properties such as the viscosity profile. To do this, we developed a new algorithm to trace plume-like features in shear-wave (Vs) seismic tomography models based on picking local minima in the velocity and searching for continuous features with depth. We applied this method to recent tomographic models and find 60+ continuous plume conduits that are > 750 km long. Approximately a third of these can be associated with known hotspots at the surface. We analyze the morphology of these continuous conduits and infer large scale mantle flow patterns and properties. We find the largest lateral deflections in the conduits occur near the base of the lower mantle and in the upper mantle (near the thermal boundary layers). The preferred orientation of the plume deflections show large variability at all depths and indicate no net mantle rotation. Plate by plate analysis shows little agreement in deflection below particular plates, indicating these deflected features might be long lived and not caused by plate shearing. Changes in the gradient of plume deflection are inferred to correspond with viscosity contrasts in the mantle and found below the transition zone as well as at 1000 km depth. From this inferred viscosity structure, we explore the dynamics of a plume through these viscosity jumps. We also retrieve the Vs profiles for the conduits and compare with the velocity profiles predicted for different mantle adiabat temperatures. We are able to constrain the average temperature anomaly of the conduits to be around 150 K. We use these thermal anomalies in conjunction with our measured plume tilts/deflections to further explore the dynamics of plume conduits in the lower mantle and transition zone.

Magnetohydrodynamic Convection in the Outer Core and its Geodynamic Consequences

NASA Technical Reports Server (NTRS)

Kuang, Weijia; Chao, Benjamin F.; Fang, Ming

2004-01-01

The Earth's fluid outer core is in vigorous convection through much of the Earth's history. In addition to generating and maintaining Earth s time-varying magnetic field (geodynamo), the core convection also generates mass redistribution in the core and a dynamical pressure field on the core-mantle boundary (CMB). All these shall result in various core-mantle interactions, and contribute to surface geodynamic observables. For example, electromagnetic core-mantle coupling arises from finite electrically conducting lower mantle; gravitational interaction occurs between the cores and the heterogeneous mantle; mechanical coupling may also occur when the CMB topography is aspherical. Besides changing the mantle rotation via the coupling torques, the mass-redistribution in the core shall produce a spatial-temporal gravity anomaly. Numerical modeling of the core dynamical processes contributes in several geophysical disciplines. It helps explain the physical causes of surface geodynamic observables via space geodetic techniques and other means, e.g. Earth's rotation variation on decadal time scales, and secular time-variable gravity. Conversely, identification of the sources of the observables can provide additional insights on the dynamics of the fluid core, leading to better constraints on the physics in the numerical modeling. In the past few years, our core dynamics modeling efforts, with respect to our MoSST model, have made significant progress in understanding individual geophysical consequences. However, integrated studies are desirable, not only because of more mature numerical core dynamics models, but also because of inter-correlation among the geophysical phenomena, e.g. mass redistribution in the outer core produces not only time-variable gravity, but also gravitational core-mantle coupling and thus the Earth's rotation variation. They are expected to further facilitate multidisciplinary studies of core dynamics and interactions of the core with other components of the Earth.

Mantle plumes - A boundary layer approach for Newtonian and non-Newtonian temperature-dependent rheologies. [modeling for island chains and oceanic aseismic ridges

NASA Technical Reports Server (NTRS)

Yuen, D. A.; Schubert, G.

1976-01-01

Stress is placed on the temperature dependence of both a linear Newtonian rheology and a nonlinear olivine rheology in accounting for narrow mantle flow structures. The boundary-layer theory developed incorporates an arbitrary temperature-dependent power-law rheology for the medium, in order to facilitate the study of mantle plume dynamics under real conditions. Thermal, kinematic, and dynamic structures of mantle plumes are modelled by a two-dimensional natural-convection boundary layer rising in a fluid with a temperature-dependent power-law relationship between shear stress and strain rate. An analytic similarity solution is arrived at for upwelling adjacent to a vertical isothermal stress-free plane. Newtonian creep as a deformation mechanism, thermal anomalies resulting from chemical heterogeneity, the behavior of plumes in non-Newtonian (olivine) mantles, and differences in the dynamics of wet and dry olivine are discussed.

Origin and thermal evolution of Mars

NASA Technical Reports Server (NTRS)

Schubert, Gerald; Soloman, S. C.; Turcotte, D. L.; Drake, M. J.; Sleep, N. H.

1990-01-01

The thermal evolution of Mars is governed by subsolidus mantle convection beneath a thick lithosphere. Models of the interior evolution are developed by parameterizing mantle convective heat transport in terms of mantle viscosity, the superadiabatic temperature rise across the mantle, and mantle heat production. Geological, geophysical, and geochemical observations of the compositon and structure of the interior and of the timing of major events in Martian evolution are used to constrain the model computations. Such evolutionary events include global differentiation, atmospheric outgassing, and the formation of the hemispherical dichotomy and Tharsis. Numerical calculations of fully three-dimensional, spherical convection in a shell the size of the Martian mantle are performed to explore plausible patterns of Martian mantel convection and to relate convective features, such as plumes, to surface features, such as Tharsis. The results from the model calculations are presented.

Preface: Deep Slab and Mantle Dynamics

NASA Astrophysics Data System (ADS)

Suetsugu, Daisuke; Bina, Craig R.; Inoue, Toru; Wiens, Douglas A.

2010-11-01

We are pleased to publish this special issue of the journal Physics of the Earth and Planetary Interiors entitled "Deep Slab and Mantle Dynamics". This issue is an outgrowth of the international symposium "Deep Slab and Mantle Dynamics", which was held on February 25-27, 2009, in Kyoto, Japan. This symposium was organized by the "Stagnant Slab Project" (SSP) research group to present the results of the 5-year project and to facilitate intensive discussion with well-known international researchers in related fields. The SSP and the symposium were supported by a Grant-in-Aid for Scientific Research (16075101) from the Ministry of Education, Culture, Sports, Science and Technology of the Japanese Government. In the symposium, key issues discussed by participants included: transportation of water into the deep mantle and its role in slab-related dynamics; observational and experimental constraints on deep slab properties and the slab environment; modeling of slab stagnation to constrain its mechanisms in comparison with observational and experimental data; observational, experimental and modeling constraints on the fate of stagnant slabs; eventual accumulation of stagnant slabs on the core-mantle boundary and its geodynamic implications. This special issue is a collection of papers presented in the symposium and other papers related to the subject of the symposium. The collected papers provide an overview of the wide range of multidisciplinary studies of mantle dynamics, particularly in the context of subduction, stagnation, and the fate of deep slabs.

Effects of iron on the lattice thermal conductivity of Earth's deep mantle and implications for mantle dynamics

NASA Astrophysics Data System (ADS)

Hsieh, Wen-Pin; Deschamps, Frédéric; Okuchi, Takuo; Lin, Jung-Fu

2018-04-01

Iron may critically influence the physical properties and thermochemical structures of Earth's lower mantle. Its effects on thermal conductivity, with possible consequences on heat transfer and mantle dynamics, however, remain largely unknown. We measured the lattice thermal conductivity of lower-mantle ferropericlase to 120 GPa using the ultrafast optical pump-probe technique in a diamond anvil cell. The thermal conductivity of ferropericlase with 56% iron significantly drops by a factor of 1.8 across the spin transition around 53 GPa, while that with 8–10% iron increases monotonically with pressure, causing an enhanced iron substitution effect in the low-spin state. Combined with bridgmanite data, modeling of our results provides a self-consistent radial profile of lower-mantle thermal conductivity, which is dominated by pressure, temperature, and iron effects, and shows a twofold increase from top to bottom of the lower mantle. Such increase in thermal conductivity may delay the cooling of the core, while its decrease with iron content may enhance the dynamics of large low shear-wave velocity provinces. Our findings further show that, if hot and strongly enriched in iron, the seismic ultralow velocity zones have exceptionally low conductivity, thus delaying their cooling.

Effects of iron on the lattice thermal conductivity of Earth's deep mantle and implications for mantle dynamics.

PubMed

Hsieh, Wen-Pin; Deschamps, Frédéric; Okuchi, Takuo; Lin, Jung-Fu

2018-04-17

Iron may critically influence the physical properties and thermochemical structures of Earth's lower mantle. Its effects on thermal conductivity, with possible consequences on heat transfer and mantle dynamics, however, remain largely unknown. We measured the lattice thermal conductivity of lower-mantle ferropericlase to 120 GPa using the ultrafast optical pump-probe technique in a diamond anvil cell. The thermal conductivity of ferropericlase with 56% iron significantly drops by a factor of 1.8 across the spin transition around 53 GPa, while that with 8-10% iron increases monotonically with pressure, causing an enhanced iron substitution effect in the low-spin state. Combined with bridgmanite data, modeling of our results provides a self-consistent radial profile of lower-mantle thermal conductivity, which is dominated by pressure, temperature, and iron effects, and shows a twofold increase from top to bottom of the lower mantle. Such increase in thermal conductivity may delay the cooling of the core, while its decrease with iron content may enhance the dynamics of large low shear-wave velocity provinces. Our findings further show that, if hot and strongly enriched in iron, the seismic ultralow velocity zones have exceptionally low conductivity, thus delaying their cooling.

Global geodynamic models constrained by tectonic reconstructions including plate deformation

NASA Astrophysics Data System (ADS)

Gurnis, M.; Flament, N.; Spasojevic, S.; Williams, S.; Seton, M.; Müller, R. D.

2011-12-01

In order to investigate the effect of mantle flow on the Earth's surface, imposing the kinematics predicted by plate reconstructions in global convection models has become common practice. Such models are valuable to investigate the effect of the mantle flow beneath the lithosphere on surface topography. Changes in surface topography due to lithospheric deformation are so far not part of top-down tectonic models in which plates are treated as rigid in traditional tectonic reconstructions. We introduce a new generation of geodynamic models that are based on tectonic reconstructions with deforming plates at both passive and convergent margins. These models allow us to investigate the relationships between lithospheric deformation and mantle flow, and their combined effects on surface topography. In traditional tectonic reconstructions, continents are represented as rigid blocks that either overlap or are separated by gaps in full-fit reconstructions. Reconstructions that include a global network of topological plate polygons avoid continental overlaps and gaps, but velocities are still derived on the basis of the Euler poles for rigid blocks. To resolve these issues, we developed a series of deforming plate models using the open source plate modeling software GPlates. For a given area, our methodology requires the relative motions between major rigid continental blocks, and a definition of the regions in which continental lithosphere deformed between these blocks. We use geophysical and geological data to define the limit between rigid and deforming areas, and the deformation history of non-rigid blocks. The velocity field predicted by these reconstructions is then used as a time-dependent surface boundary condition in global 3-D geodynamic models. To incorporate the continental lithosphere in our global models, we embed compositionally distinct crust and continental lithosphere within the thermal lithosphere. We define three isostatic columns of different thickness and buoyancy based on the tectonothermal age of the continents: Archean, Proterozoic and Phanerozoic. In the fourth isostatic column, the oceans, the thickness of the thermal lithosphere is assimilated using the half-space cooling model. We also use this capacity to define the thickness of the thermal lithosphere for different continental types, with the exception of the deforming areas that are fully dynamic. Finally, we introduce a new slab assimilation method in which the thermal structure of the slab, derived analytically, is progressively assimilated in the upper mantle into the dynamic models. This method not only improves the continuity of slabs in our models, but it also allows us to model flat slab segments that are particularly relevant for dynamic topography. This new generation of models allows us to analyse the contributions of continental deformation and of mantle flow to surface topography. We compare our results to geological and geophysical data, including stratigraphy, paleo-altimetry, paleo-environment and mantle tomography. This allows us to place constraints on key model parameters and to refine our knowledge of plate-mantle interactions during continental deformation.

Complex interactions between diapirs and 4-D subduction driven mantle wedge circulation.

NASA Astrophysics Data System (ADS)

Sylvia, R. T.; Kincaid, C. R.

2015-12-01

Analogue laboratory experiments generate 4-D flow of mantle wedge fluid and capture the evolution of buoyant mesoscale diapirs. The mantle is modeled with viscous glucose syrup with an Arrhenius type temperature dependent viscosity. To characterize diapir evolution we experiment with a variety of fluids injected from multiple point sources. Diapirs interact with kinematically induced flow fields forced by subducting plate motions replicating a range of styles observed in dynamic subduction models (e.g., rollback, steepening, gaps). Data is collected using high definition timelapse photography and quantified using image velocimetry techniques. While many studies assume direct vertical connections between the volcanic arc and the deeper mantle source region, our experiments demonstrate the difficulty of creating near vertical conduits. Results highlight extreme curvature of diapir rise paths. Trench-normal deflection occurs as diapirs are advected downward away from the trench before ascending into wedge apex directed return flow. Trench parallel deflections up to 75% of trench length are seen in all cases, exacerbated by complex geometry and rollback motion. Interdiapir interaction is also important; upwellings with similar trajectory coalesce and rapidly accelerate. Moreover, we observe a new mode of interaction whereby recycled diapir material is drawn down along the slab surface and then initiates rapid fluid migration updip along the slab-wedge interface. Variability in trajectory and residence time leads to complex petrologic inferences. Material from disparate source regions can surface at the same location, mix in the wedge, or become fully entrained in creeping flow adding heterogeneity to the mantle. Active diapirism or any other vertical fluid flux mechanism employing rheological weakening lowers viscosity in the recycling mantle wedge affecting both solid and fluid flow characteristics. Many interesting and insightful results have been presented based upon 2-D, steady-state thermal and flow regimes. We reiterate the importance of 4-D time evolution in subduction models. Analogue experiments allow added feedbacks and complexity improving intuition and providing insight for further investigation.

Supercontinent Pangea, Mantle Dynamics, and Reference Frame of Global Plate Motions

NASA Astrophysics Data System (ADS)

Zhong, S.; Rudolph, M. L.; Liu, X.

2014-12-01

Arguably the most important and challenging goal in geodynamics is to understand the two-way dynamics between tectonic plates and mantle convection. While it has long been recognized that the present-day degree-2 mantle structure as imaged seismically is closely related to the plate motions (Hager and O'Connell, 1981) and their history (<119 Ma) (Ricard et al., 1993; McNamara and Zhong, 2005), recent studies have expanded this concept, from two different perspectives, by seeking connections between Pangea assembly and breakup and mantle structure and dynamics. First, it has been proposed that the large igneous provinces (LIPs) and kimberlite volcanism erupted mainly along the edges of the two major seismically slow anomalies above the core-mantle boundary (often referred to as the Africa and Pacific LLSVPs) (Torsvik et al, 2010). This has led to the proposal that the present-day degree-2 mantle structure has existed for >500 Ma (Torsvik et al., 2014), although its statistical significance has been challenged (Austermann et al., 2013). The proposals of the spatially stable Africa and Pacific LLSVPs and of the LIP eruptions along their edges have also been exploited in attempts to build global plate motion models since the Pangea assembly by providing a plate motion reference frame or inferring true polar wander (TPW) corrections to the plate motions (Torsvik et al., 2014). Second, mantle dynamics studies indicate that degree-1 mantle convection, which is expected with realistic lithospheric and mantle viscosity, may be needed for assembly of a supercontinent (e.g., Pangea) (Zhong et al., 2007). This suggests that the present degree-2 mantle structure may have been formed only after the Pangea assembly from an initially degree-1 structure - a scenario that is consistent with convection calculations with a proxy plate motion model that considers Pangea process (Zhang et al., 2010). In this presentation, in addition to critically reviewing these arguments, we will discuss calculations of long-wavelength geoid for the mantle with thermochemical piles and LLSVPs and their potential effects on TPW determinations and hence reconstruction of plate motion (i.e., net lithospheric rotation). We will also present additional calculations of mantle structure evolution using different plate motion history models.

Teaching machines to find mantle composition

NASA Astrophysics Data System (ADS)

Atkins, Suzanne; Tackley, Paul; Trampert, Jeannot; Valentine, Andrew

2017-04-01

The composition of the mantle affects many geodynamical processes by altering factors such as the density, the location of phase changes, and melting temperature. The inferences we make about mantle composition also determine how we interpret the changes in velocity, reflections, attenuation and scattering seen by seismologists. However, the bulk composition of the mantle is very poorly constrained. Inferences are made from meteorite samples, rock samples from the Earth and inferences made from geophysical data. All of these approaches require significant assumptions and the inferences made are subject to large uncertainties. Here we present a new method for inferring mantle composition, based on pattern recognition machine learning, which uses large scale in situ observations of the mantle to make fully probabilistic inferences of composition for convection simulations. Our method has an advantage over other petrological approaches because we use large scale geophysical observations. This means that we average over much greater length scales and do not need to rely on extrapolating from localised samples of the mantle or planetary disk. Another major advantage of our method is that it is fully probabilistic. This allows us to include all of the uncertainties inherent in the inference process, giving us far more information about the reliability of the result than other methods. Finally our method includes the impact of composition on mantle convection. This allows us to make much more precise inferences from geophysical data than other geophysical approaches, which attempt to invert one observation with no consideration of the relationship between convection and composition. We use a sampling based inversion method, using hundreds of convection simulations run using StagYY with self consistent mineral physics properties calculated using the PerpleX package. The observations from these simulations are used to train a neural network to make a probabilistic inference for major element oxide composition of the mantle. We find we can constrain bulk mantle FeO molar percent, FeO/MgO and FeO/SiO2 using observations of the temperature and density structure of the mantle in convection simulations.

Report of the panel on earth structure and dynamics, section 6

NASA Technical Reports Server (NTRS)

Dziewonski, Adam M.; Mcadoo, David C.; Oconnell, Richard J.; Smylie, Douglas E.; Yoder, Charles F.

1991-01-01

The panel identified problems related to the dynamics of the core and mantle that should be addressed by NASA programs. They include investigating the geodynamo based on observations of the Earth's magnetic field, determining the rheology of the mantle from geodetic observations of post-glacial vertical motions and changes in the gravity field, and determining the coupling between plate motions and mantle flow from geodetic observations of plate deformation. Also emphasized is the importance of support for interdisciplinary research to combine various data sets with models which couple rheology, structure and dynamics.

Blending geological observations and convection models to reconstruct mantle dynamics

NASA Astrophysics Data System (ADS)

Coltice, Nicolas; Bocher, Marie; Fournier, Alexandre; Tackley, Paul

2015-04-01

Knowledge of the state of the Earth mantle and its temporal evolution is fundamental to a variety of disciplines in Earth Sciences, from the internal dynamics to its many expressions in the geological record (postglacial rebound, sea level change, ore deposit, tectonics or geomagnetic reversals). Mantle convection theory is the centerpiece to unravel the present and past state of the mantle. For the past 40 years considerable efforts have been made to improve the quality of numerical models of mantle convection. However, they are still sparsely used to estimate the convective history of the solid Earth, in comparison to ocean or atmospheric models for weather and climate prediction. The main shortcoming is their inability to successfully produce Earth-like seafloor spreading and continental drift self-consistently. Recent convection models have begun to successfully predict these processes. Such breakthrough opens the opportunity to retrieve the recent dynamics of the Earth's mantle by blending convection models together with advanced geological datasets. A proof of concept will be presented, consisting in a synthetic test based on a sequential data assimilation methodology.

Reconstructing mantle flow and long-wavelength dynamic topography since the Jurassic Period (GD Division Outstanding ECS Award Lecture)

NASA Astrophysics Data System (ADS)

Flament, Nicolas

2017-04-01

Global tectonic reconstructions can be used as boundary conditions of forward mantle convection models to simulate past mantle flow and long-wavelength dynamic topography. The predictions of such models can be compared to seismic tomography, to estimates of residual topography and to geological indicators of past vertical motions. Here we present models that reproduce the present-day structure of the lower mantle, including two large structures that resemble the Pacific and African Large Low Shear Velocity Provinces (LLSVPs, ˜15,000 km in diameter) and a smaller structure that resembles the recently discovered Perm Anomaly (˜1,000 km in diameter). The match between predicted and seismically inferred lower mantle structure is quantified across a series of mantle flow and tomography models. In the models, the Perm-like anomaly forms in isolation within a closed and long-lived subduction network (East Asia, Northern Tethys and Mongol-Okhotsk) ˜22,000 km in circumference before migrating ˜1,500 km westward at an average rate of 1 cm yr-1 since 150 million years ago. These results indicate a greater mobility of deep mantle structures than previously recognized, and illustrate that the predictive power of mantle flow models has significantly increased over the last thirty years. We suggest that the mobile Perm Anomaly could be linked to the ˜258 Ma Emeishan volcanics, in contrast to the previously proposed ˜251 Ma Siberian Traps. We also compare the present-day dynamic topography predicted by forward mantle flow models to residual topography models, and show that radial and lateral viscosity variations significantly influence the distribution of power of predicted dynamic topography as a function of spherical harmonic degree. We finally show how past vertical motions preserved in the geological record and the present-day position of slabs in the mantle inferred from seismic tomography may be used to constrain tectonic reconstructions and mantle rheology, including examples focusing on the large-scale topographic asymmetry of the South Atlantic domain and on the uplift history of the eastern highlands of Australia.

Density Anomalies in the Mantle and the Gravitational Core-Mantle Interaction

NASA Technical Reports Server (NTRS)

Kuang, Weijia; Liu, Lanbo

2003-01-01

Seismic studies suggest that the bulk of the mantle is heterogeneous, with density variations in depth as well as in horizontal directions (latitude and longitude). This density variation produces a three- dimensional gravity field throughout the Earth. On the other hand, the core density also varies in both time and space, due to convective core flow. Consequently, the fluid outer core and the solid mantle interact gravitationally due to the mass anomalies in both regions. This gravitational core-mantle interaction could play a significant role in exchange of angular momentum between the core and the mantle, and thus the change in Earth's rotation on time scales of decades and longer. Aiming at estimating the significance of the gravitational core-mantle interaction on Earth's rotation variation, we introduce in our MoSST core dynamics model a heterogeneous mantle, with a density distribution derived from seismic results. In this model, the core convection is driven by the buoyancy forces. And the density variation is determined dynamically with the convection. Numerical simulation is carried out with different parameter values, intending to extrapolate numerical results for geophysical implications.

Is Earth coming out of the recent ice house age in the long-term? - constraints from probable mantle CO2-degassing reconstructions

NASA Astrophysics Data System (ADS)

Hartmann, Jens; Li, Gaojun; West, A. Joshua

2017-04-01

Enhanced partial melting of mantle material probably started when the subduction motor started around 3.2 Ga ago as evidenced by the formation history of the continental crust. Carbon is degassing due partial melting as it is an incompatible element. Therefore, mantle carbon degassing rates would change with time proportionally to the reservoir mantle concentration evolution and the ocean crust production rate, causing a distinct CO2-degassing rate change with time. The evolution of the mantle degassing rate has some implications for the reconstruction of the carbon cycle and therefore climate and Earth surface processes rates, as CO2-degassing rates are used to constrain or to balance the atmosphere-ocean-crust carbon cycle system. It will be shown that compilations of CO2-degassing from relevant geological sources are probably exceeding the established CO2-sink terrestrial weathering, which is often used to constrain long-term mantle degassing rates to close the carbon cycle on geological time scales. In addition, the scenarios for the degassing dynamics from the mantle sources suggest that the mantle is depleting its carbon content since 3 Ga. This has further implications for the long-term CO2-sink weathering. Results will be compared with geochemical proxies for weathering and weathering intensity dynamics, and will be set in context with snow ball Earth events and long-term emplacement dynamics of mafic areas as Large Igneous Provinces. Decreasing mantle degassing rates since about 2 Ga suggest a constraint for the evolution of the carbon cycle and recycling potential of the amount of subducted carbon. If the given scenarios hold further investigation, the contribution of mantle degassing to climate forcing (directly and via recycling) will decrease further.

Earth Structure, Ice Mass Changes, and the Local Dynamic Geoid

NASA Astrophysics Data System (ADS)

Harig, C.; Simons, F. J.

2014-12-01

Spherical Slepian localization functions are a useful method for studying regional mass changes observed by satellite gravimetry. By projecting data onto a sparse basis set, the local field can be estimated more easily than with the full spherical harmonic basis. We have used this method previously to estimate the ice mass change in Greenland from GRACE data, and it can also be applied to other planetary problems such as global magnetic fields. Earth's static geoid, in contrast to the time-variable field, is in large part related to the internal density and rheological structure of the Earth. Past studies have used dynamic geoid kernels to relate this density structure and the internal deformation it induces to the surface geopotential at large scales. These now classical studies of the eighties and nineties were able to estimate the mantle's radial rheological profile, placing constraints on the ratio between upper and lower mantle viscosity. By combining these two methods, spherical Slepian localization and dynamic geoid kernels, we have created local dynamic geoid kernels which are sensitive only to density variations within an area of interest. With these kernels we can estimate the approximate local radial rheological structure that best explains the locally observed geoid on a regional basis. First-order differences of the regional mantle viscosity structure are accessible to this technique. In this contribution we present our latest, as yet unpublished results on the geographical and temporal pattern of ice mass changes in Antarctica over the past decade, and we introduce a new approach to extract regional information about the internal structure of the Earth from the static global gravity field. Both sets of results are linked in terms of the relevant physics, but also in being developed from the marriage of Slepian functions and geoid kernels. We make predictions on the utility of our approach to derive fully three-dimensional rheological Earth models, to be used for corrections for glacio-isostatic adjustment, as necessary for the interpretation of time-variable gravity observations in terms of ice sheet mass-balance studies.

Retrodicting the Cenozoic evolution of the mantle: Implications for dynamic surface topography

NASA Astrophysics Data System (ADS)

Glišović, Petar; Forte, Alessandro; Rowley, David; Simmons, Nathan; Grand, Stephen

2014-05-01

Seismic tomography is the essential starting ingredient for constructing realistic models of the mantle convective flow and for successfully predicting a wide range of convection-related surface observables. However, the lack of knowledge of the initial thermal state of the mantle in the geological past is still an outstanding problem in mantle convection. The resolution of this problem requires models of 3-D mantle evolution that yield maximum consistency with a wide suite of geophysical constraints. Quantifying the robustness of the reconstructed thermal evolution is another major concern. We have carried out mantle dynamic simulations (Glišović & Forte, EPSL 2014) using a pseudo-spectral solution for compressible-flow thermal convection in 3-D spectral geometry that directly incorporate: 1) joint seismic-geodynamic inversions of mantle density structure with constraints provided by mineral physics data (Simmons et al., GJI 2009); and 2) constraints on mantle viscosity inferred by inversion of a suite of convection-related and glacial isostatic adjustment data sets (Mitrovica & Forte, EPSL 2004) characterised by Earth-like Rayleigh numbers. These time-reversed convection simulations reveal how the buoyancy associated with hot, active upwellings is a major driver of the mantle-wide convective circulation and the changes in dynamic topography at the Earth's surface. These simulations reveal, for example, a stable and long-lived superplume under the East Pacific Rise (centred under the Easter and Pitcairn hotspots) that was previously identified by Rowley et al. (AGU 2011, Nature in review) on the basis of plate kinematic data. We also present 65 Myr reconstructions of the Reunion plume that gave rise to the Deccan Traps.

Dynamical consequences of mantle heterogeneity in two-phase models of mid-ocean ridges

NASA Astrophysics Data System (ADS)

Katz, R. F.

2010-12-01

The mid-ocean ridge system, over 50,000 km in length, samples the magmatic products of a large swath of the asthenosphere. It provides our best means to assess the heterogeneity structure of the upper mantle. Interpretation of the diverse array of observations of MOR petrology, geochemistry, tomography, etc requires models that can map heterogeneity structure onto predictions testable by comparison with these observations. I report on progress to this end; in particular, I describe numerical models of coupled magma/mantle dynamics at mid-ocean ridges [1,2]. These models incorporate heterogeneity in terms of a simple, two-component thermochemical system with specified amplitude and spatial distribution. They indicate that mantle heterogeneity has significant fluid-dynamical consequences for both mantle and magmatic flow. Models show that the distribution of enrichment can lead to asymmetry in the strength of upwelling across the ridge-axis and channelised magmatic transport to the axis. Furthermore, heterogeneity can cause off-axis upwelling of partially molten diapirs, trapping of enriched melts off-axis, and re-fertilization of the mantle by pooled and refrozen melts. Predicted consequences of geochemical heterogeneity may also be considered. References: [1] Katz, RF, (2008); Magma dynamics with the Enthalpy Method: Benchmark Solutions and Magmatic Focusing at Mid-ocean Ridges. Journal of Petrology, doi: 10.1093/petrology/egn058. [2] Katz RF, (2010); Porosity-driven convection and asymmetry beneath mid-ocean ridges. Submitted to G3.

Early evolution and dynamics of Earth from a molten initial stage

NASA Astrophysics Data System (ADS)

Lourenço, Diogo L.; Tackley, Paul J.

2014-05-01

It is now well established that most of the terrestrial planets underwent a magma ocean stage during their accretion. On Earth, it is probable that at the end of accretion, giant impacts like the hypothesised Moon-forming impact, together with other sources of heat such as conversion of gravitational energy of formation into heat, heat losses from the core at the core-mantle boundary, radioactive decay, electromagnetic induction heating and tidal heating, melted a substantial part of the mantle. The thermal and chemical evolution of the resulting magma ocean most certainly had dramatic consequences on the history of the planet, influencing the chemical composition of the mantle after differentiation, the style of tectonic regime prevailing in the solid-state mantle and its habitability. Considerable research has been done on magma oceans using 1-D models (e.g.: Abe, PEPI 1997; Solomatov, Treat. Geophys. 2007; Elkins-Tanton, EPSL 2008). However, its dynamics, evolution from a molten state to the present day solid state, and crystallisation are still not fully understood and are more complex than a 1-D formulation. Recent advances in computational methods and resources allow us to address numerically more complex problems, with higher resolution and multiple physics incorporated. Moreover, new developments in mineral physics that indicate that melt can be denser than solid at high pressures (e.g.: de Koker et al., EPSL 2013) can have very important impacts on the classical views of the solidification of magma oceans (Labrosse et al., Nature 2007). The goal of our study is to understand and characterise the influence of melting on the long-term thermo-chemical evolution of rocky planet interiors, starting from an initial molten state (magma ocean). Our approach is to test existing published 1-D parameterisations of magma ocean dynamics and extend them into 2-D models. We will address this problem using the numerical code StagYY (Tackley, PEPI 2008), which uses a finite-volume scheme for advection of temperature, a multigrid solver to obtain a velocity-pressure solution at each timestep, tracers to track composition, and a treatment of partial melting and crustal formation. Additional enhancements are needed in the code and are related to the physics and parameterisation of melting.

Topography: dusting for the fingerprints of mantle dynamics

NASA Astrophysics Data System (ADS)

Faccenna, C.; Becker, T. W.

2016-12-01

The surface of the Earth is an ever-changing expression of the dynamic processes occurring deep in the mantle and at and above its surface, but our ability to "read" landscapes in terms of their underlying tectonic or climatic forcing is rudimentary. During the last decade, particular attention has been drawn to the deep, convection-related component of topography, induced by the stress produced at the base of the lithosphere by mantle flow, and its relevance compared to the (iso)static component. Despite much progress, several issues, including the magnitude and rate of this dynamic component, remain open. Here, we use key sites from convergent margins (e.g., the Apennines) and from intraplate settings (e.g., Ethiopia) to estimate the amplitude and rate of topography change and to disentangle the dynamic from the static component. On the base of those and other examples, we introduce the concept of a Topographic Fingerprint: any combination of mantle, crustal and surface processes that will result in a distinctive, thus predictable, topographic expression.

Impact of lithosphere rheology on the dynamic topography

NASA Astrophysics Data System (ADS)

Burov, Evgueni; Gerya, Taras; Koptev, Alexander

2014-05-01

Dynamic topography is a key observable signature of the Earth's and planetary (e.g. Venus) mantle dynamics. In general view, it reflects complex mantle flow patterns, and hence is supposed to correlate at different extent with seismic tomography, SKS fast orientations, geodetic velocity fields and geoid anomalies. However, identification of dynamic topography had no systematic success, specifically in the Earth's continents. Here we argue that lithosphere rheology, in particular, rheological stratification of continents, results in modulation of dynamic topography, converting commonly expected long-wavelength/small amplitude undulations into short-wavelength surface undulations with wide amplitude spectrum, superimposed onto "tectonic" topography. These ideas are explored in 3D using unprecedentedly high resolution numerical experiments (grid step size 2-3 km for 1500x1500x600 km computational area) incorporating realistic rheologically stratified lithosphere. Such high resolution is actually needed to resolve small-scale crustal faulting and inter-layer coupling/uncoupling that shape surface topography. The results reveal strikingly discordant, counterintuitive features of 3D dynamic topography, going far beyond the inferences from previous models. In particular, even weak anisotropic tectonic stress field results both in large-scale small-amplitude dynamic topography and in strongly anisotropic short-wavelength (at least in one direction) dynamic topography with wide amplitude range (from 100 to 2000-3000 m), including basins and ranges and large-scale linear normal and strike-slip faults. Even very slightly pre-stressed strong lithosphere yields and localizes deformation much easier , than un-prestressed one, in response to plume impact and mantle flow. The results shed new light on the importance of lithosphere rheology and active role of lithosphere in mantle-lithosphere interactions as well as on the role of mantle flow and far-field stresses in tectonic-scale deformation. We show, for example, that crustal fault patterns initiated by plume impact are rapidly re-organized in sub-linear rifts and spreading centers, which orientation is largely dictated (e.g., perpendicular to) by the direction of the tectonic far-field stress field, as well as the plume-head material soon starts to flow along the sub-linear rifted shear zones in crustal and mantle lithosphere further amplifying their development. The final surface deformation and mantle flow patterns rapidly loose the initial axisymmetric character and take elongated sub-linear shapes whereas brittle deformation at surface is amplified and stabilized by coherent flow of mantle/plume-head material from below. These "tectonically" looking dynamic topography patterns are quite different from those expected from conventional models as well as from those directly observed, for example, on Venus where plume-lithosphere interactions produce only axisymmetric coronae domal-shaped features with radiating extensional rifts, suggesting that the Venusian lithosphere is rheologically too weak , and its crust is too thin, to produce any significant impact on the dynamic topography.

Mantle dynamics following supercontinent formation

NASA Astrophysics Data System (ADS)

Heron, Philip J.

This thesis presents mantle convection numerical simulations of supercontinent formation. Approximately 300 million years ago, through the large-scale subduction of oceanic sea floor, continental material amalgamated to form the supercontinent Pangea. For 100 million years after its formation, Pangea remained relatively stationary, and subduction of oceanic material featured on its margins. The present-day location of the continents is due to the rifting apart of Pangea, with supercontinent dispersal being characterized by increased volcanic activity linked to the generation of deep mantle plumes. The work presented here investigates the thermal evolution of mantle dynamics (e.g., mantle temperatures and sub-continental plumes) following the formation of a supercontinent. Specifically, continental insulation and continental margin subduction are analyzed. Continental material, as compared to oceanic material, inhibits heat flow from the mantle. Previous numerical simulations have shown that the formation of a stationary supercontinent would elevate sub-continental mantle temperatures due to the effect of continental insulation, leading to the break-up of the continent. By modelling a vigorously convecting mantle that features thermally and mechanically distinct continental and oceanic plates, this study shows the effect of continental insulation on the mantle to be minimal. However, the formation of a supercontinent results in sub-continental plume formation due to the re-positioning of subduction zones to the margins of the continent. Accordingly, it is demonstrated that continental insulation is not a significant factor in producing sub-supercontinent plumes but that subduction patterns control the location and timing of upwelling formation. A theme throughout the thesis is an inquiry into why geodynamic studies would produce different results. Mantle viscosity, Rayleigh number, continental size, continental insulation, and oceanic plate boundary evolution are explored in over 600 2D and over 20 3D numerical simulations to better understand how modelling method affects conclusions on mantle convection studies. The results from this thesis show that the failure to model tectonic plates, a high vigour of convection, and a (pseudo) temperature-dependent viscosity would distort the role of mantle plumes, continent insulation, and subduction in the thermal evolution of mantle dynamics.

Petrological Geodynamics of Mantle Melting II. AlphaMELTS + Multiphase Flow: Dynamic Fractional Melting

NASA Astrophysics Data System (ADS)

Tirone, Massimiliano

2018-03-01

In this second installment of a series that aims to investigate the dynamic interaction between the composition and abundance of the solid mantle and its melt products, the classic interpretation of fractional melting is extended to account for the dynamic nature of the process. A multiphase numerical flow model is coupled with the program AlphaMELTS, which provides at the moment possibly the most accurate petrological description of melting based on thermodynamic principles. The conceptual idea of this study is based on a description of the melting process taking place along a 1-D vertical ideal column where chemical equilibrium is assumed to apply in two local sub-systems separately on some spatial and temporal scale. The solid mantle belongs to a local sub-system (ss1) that does not interact chemically with the melt reservoir which forms a second sub-system (ss2). The local melt products are transferred in the melt sub-system ss2 where the melt phase eventually can also crystallize into a different solid assemblage and will evolve dynamically. The main difference with the usual interpretation of fractional melting is that melt is not arbitrarily and instantaneously extracted from the mantle, but instead remains a dynamic component of the model, hence the process is named dynamic fractional melting (DFM). Some of the conditions that may affect the DFM model are investigated in this study, in particular the effect of temperature, mantle velocity at the boundary of the mantle column. A comparison is made with the dynamic equilibrium melting (DEM) model discussed in the first installment. The implications of assuming passive flow or active flow are also considered to some extent. Complete data files of most of the DFM simulations, four animations and two new DEM simulations (passive/active flow) are available following the instructions in the supplementary material.

Isotopic evidence for a large-scale plume-derived mantle domain between the Indian and Pacific mantles beneath the Southern Ocean.

NASA Astrophysics Data System (ADS)

Park, S. H.; Langmuir, C. H.; Scott, S. R.; Sims, K. W. W.; Lin, J.; Kim, S. S.; Blichert-Toft, J.; Choi, H.; Yang, Y. S.; Michael, P. J.

2017-12-01

Earth's upper mantle is characterized by Indian- and Pacific-type domains with distinctive isotope characteristics. The boundary between these two mantle regions has been hypothesized to be located at the Australian-Antarctic-Discordance (AAD), where regions west and east of the AAD are Indian- and Pacific-type, respectively. It was further posited that the Pacific mantle feeds into the Indian mantle as the boundary is moving westward. These scenarios have important implications for the dynamics of mantle convection in the area. In the present model, regions east of the AAD are assumed to be entirely Pacific-type mantle, but our recent recovery of basalts from a 2,000-km sampling gap along the Australian-Antarctic Ridge (AAR), located east of the AAD on the Pacific side, challenges this picture. Here we show that the Hf, Nd, Pb, and Sr isotopic compositions of AAR MORB are distinct from those of Pacific and Indian MORB. Rather, the AAR lavas show mixing relationships with volcanoes from the Hikurangi seamounts, the Balleney and Scott Islands, the West Antarctic Rift System, New Zealand, and east Australia. According to tectonic reconstruction models, these volcanoes are related to super-plume activity that caused Gondwana to break up at 90 Ma. These results imply that a large-scale plume-derived mantle domain exists between the Indian and Pacific mantle domains, and that mantle dynamics along the AAD should be reinterpreted in light of interaction with a super-plume.

The Earth's Mantle.

ERIC Educational Resources Information Center

McKenzie, D. P.

1983-01-01

The nature and dynamics of the earth's mantle is discussed. Research indicates that the silicate mantle is heated by the decay of radioactive isotopes and that the heat energizes massive convention currents in the upper 700 kilometers of the ductile rock. These currents and their consequences are considered. (JN)

Thermally-Driven Mantle Plumes Reconcile Hot-spot Observations

NASA Astrophysics Data System (ADS)

Davies, D.; Davies, J.

2008-12-01

Hot-spots are anomalous regions of magmatism that cannot be directly associated with plate tectonic processes (e.g. Morgan, 1972). They are widely regarded as the surface expression of upwelling mantle plumes. Hot-spots exhibit variable life-spans, magmatic productivity and fixity (e.g. Ito and van Keken, 2007). This suggests that a wide-range of upwelling structures coexist within Earth's mantle, a view supported by geochemical and seismic evidence, but, thus far, not reproduced by numerical models. Here, results from a new, global, 3-D spherical, mantle convection model are presented, which better reconcile hot-spot observations, the key modification from previous models being increased convective vigor. Model upwellings show broad-ranging dynamics; some drift slowly, while others are more mobile, displaying variable life-spans, intensities and migration velocities. Such behavior is consistent with hot-spot observations, indicating that the mantle must be simulated at the correct vigor and in the appropriate geometry to reproduce Earth-like dynamics. Thermally-driven mantle plumes can explain the principal features of hot-spot volcanism on Earth.

A dynamic model of Venus's gravity field

NASA Technical Reports Server (NTRS)

Kiefer, W. S.; Richards, M. A.; Hager, B. H.; Bills, B. G.

1984-01-01

Unlike Earth, long wavelength gravity anomalies and topography correlate well on Venus. Venus's admittance curve from spherical harmonic degree 2 to 18 is inconsistent with either Airy or Pratt isostasy, but is consistent with dynamic support from mantle convection. A model using whole mantle flow and a high viscosity near surface layer overlying a constant viscosity mantle reproduces this admittance curve. On Earth, the effective viscosity deduced from geoid modeling increases by a factor of 300 from the asthenosphere to the lower mantle. These viscosity estimates may be biased by the neglect of lateral variations in mantle viscosity associated with hot plumes and cold subducted slabs. The different effective viscosity profiles for Earth and Venus may reflect their convective styles, with tectonism and mantle heat transport dominated by hot plumes on Venus and by subducted slabs on Earth. Convection at degree 2 appears much stronger on Earth than on Venus. A degree 2 convective structure may be unstable on Venus, but may have been stabilized on Earth by the insulating effects of the Pangean supercontinental assemblage.

Effects of upper mantle heterogeneities on the lithospheric stress field and dynamic topography

NASA Astrophysics Data System (ADS)

Osei Tutu, Anthony; Steinberger, Bernhard; Sobolev, Stephan V.; Rogozhina, Irina; Popov, Anton A.

2018-05-01

The orientation and tectonic regime of the observed crustal/lithospheric stress field contribute to our knowledge of different deformation processes occurring within the Earth's crust and lithosphere. In this study, we analyze the influence of the thermal and density structure of the upper mantle on the lithospheric stress field and topography. We use a 3-D lithosphere-asthenosphere numerical model with power-law rheology, coupled to a spectral mantle flow code at 300 km depth. Our results are validated against the World Stress Map 2016 (WSM2016) and the observation-based residual topography. We derive the upper mantle thermal structure from either a heat flow model combined with a seafloor age model (TM1) or a global S-wave velocity model (TM2). We show that lateral density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting dynamic topography that appears more sensitive to such heterogeneities. The modeled stress field directions, using only the mantle heterogeneities below 300 km, are not perturbed much when the effects of lithosphere and crust above 300 km are added. In contrast, modeled stress magnitudes and dynamic topography are to a greater extent controlled by the upper mantle density structure. After correction for the chemical depletion of continents, the TM2 model leads to a much better fit with the observed residual topography giving a good correlation of 0.51 in continents, but this correction leads to no significant improvement of the fit between the WSM2016 and the resulting lithosphere stresses. In continental regions with abundant heat flow data, TM1 results in relatively small angular misfits. For example, in western Europe the misfit between the modeled and observation-based stress is 18.3°. Our findings emphasize that the relative contributions coming from shallow and deep mantle dynamic forces are quite different for the lithospheric stress field and dynamic topography.

Modeling Diverse Pathways to Age Progressive Volcanism in Subduction Zones.

NASA Astrophysics Data System (ADS)

Kincaid, C. R.; Szwaja, S.; Sylvia, R. T.; Druken, K. A.

2015-12-01

One of the best, and most challenging clues to unraveling mantle circulation patterns in subduction zones comes in the form of age progressive volcanic and geochemical trends. Hard fought geological data from many subduction zones, like Tonga-Lau, the Cascades and Costa-Rica/Nicaragua, reveal striking temporal patterns used in defining mantle flow directions and rates. We summarize results from laboratory subduction models showing a range in circulation and thermal-chemical transport processes. These interaction styles are capable of producing such trends, often reflecting apparent instead of actual mantle velocities. Lab experiments use a glucose working fluid to represent Earth's upper mantle and kinematically driven plates to produce a range in slab sinking and related wedge transport patterns. Kinematic forcing assumes most of the super-adiabatic temperature gradient available to drive major downwellings is in the tabular slabs. Moreover, sinking styles for fully dynamic subduction depend on many complicating factors that are only poorly understood and which can vary widely even for repeated parameter combinations. Kinematic models have the benefit of precise, repeatable control of slab motions and wedge flow responses. Results generated with these techniques show the evolution of near-surface thermal-chemical-rheological heterogeneities leads to age progressive surface expressions in a variety of ways. One set of experiments shows that rollback and back-arc extension combine to produce distinct modes of linear, age progressive melt delivery to the surface through a) erosion of the rheological boundary layer beneath the overriding plate, and deformation and redistribution of both b) mantle residuum produced from decompression melting and c) formerly active, buoyant plumes. Additional experiments consider buoyant diapirs rising in a wedge under the influence of rollback, back-arc spreading and slab-gaps. Strongly deflected diapirs, experiencing variable rise rates, also commonly surface as linear, age progressive tracks. Applying these results to systems like the Cascades and Tonga-Lau suggest there are multiple ways to produce timing trends due both to linear flows and waves of heterogeneity obliquely impacting surface plates.

Melting and reactive flow of a volatilized mantle beneath mid-ocean ridges: theory and numerical models

NASA Astrophysics Data System (ADS)

Keller, Tobias; Katz, Richard F.

2015-04-01

Laboratory experiments indicate that even small concentrations volatiles (H2O or CO2) in the upper mantle significantly affect the silicate melting behavior [HK96,DH06]. The presence of volatiles stabilizes volatile-rich melt at high pressure, thus vastly increasing the volume of the upper mantle expected to be partially molten [H10,DH10]. These small-degree melts have important consequences for chemical differentiation and could affect the dynamics of mantle flow. We have developed theory and numerical implementation to simulate thermo-chemically coupled magma/mantle dynamics in terms of a two-phase (rock+melt), three component (dunite+MORB+volatilized MORB) physical model. The fluid dynamics is based on McKenzie's equations [McK84], while the thermo-chemical formulation of the system is represented by a novel disequilibrium multi-component melting model based on thermo-dynamic theory [RBS11]. This physical model is implemented as a parallel, two-dimensional, finite-volume code that leverages tools from the PETSc toolkit. Application of this simulation code to a mid-ocean ridge system suggests that the methodology captures the leading-order features of both hydrated and carbonated mantle melting, including deep, low-degree, volatile-rich melt formation. Melt segregation leads to continuous dynamic thermo-chemical dis-equilibration, while phenomenological reaction rates are applied to continually move the system towards re-equilibration. The simulations will be used first to characterize volatile extraction from the MOR system assuming a chemically homogeneous mantle. Subsequently, simulations will be extended to investigate the consequences of heterogeneity in lithology [KW12] and volatile content. These studies will advance our understanding of the role of volatiles in the dynamic and chemical evolution of the upper mantle. Moreover, they will help to gauge the significance of the coupling between the deep carbon cycle and the ocean/atmosphere system. REFERENCES HK96 Hirth & Kohlstedt (1996), Earth Planet Sci Lett DH06 Dasgupta & Hirschmann (2006), doi:10.1038/nature04612. H10 Hirschmann (2010), doi:10.1016/j.pepi.2009.12.003. DH10 Dasgupta & Hirschmann (2010), doi:10.1016/j.epsl.2010.06.039. McK84 McKenzie (1984), J Pet KW12 Katz & Weatherley (2012), doi: 10.1016/j.epsl.2012.04.042. RBS11 Rudge, Bercovici & Spiegelman (2011), doi: 10.1111/j.1365-246X.2010.04870.x

Reconciling Long-Wavelength Dynamic Topography, Geoid Anomalies and Mass Distribution on Earth

NASA Astrophysics Data System (ADS)

Hoggard, M.; Richards, F. D.; Ghelichkhan, S.; Austermann, J.; White, N.

2017-12-01

Since the first satellite observations in the late 1950s, we have known that that the Earth's non-hydrostatic geoid is dominated by spherical harmonic degree 2 (wavelengths of 16,000 km). Peak amplitudes are approximately ± 100 m, with highs centred on the Pacific Ocean and Africa, encircled by lows in the vicinity of the Pacific Ring of Fire and at the poles. Initial seismic tomography models revealed that the shear-wave velocity, and therefore presumably the density structure, of the lower mantle is also dominated by degree 2. Anti-correlation of slow, probably low density regions beneath geoid highs indicates that the mantle is affected by large-scale flow. Thus, buoyant features are rising and exert viscous normal stresses that act to deflect the surface and core-mantle boundary (CMB). Pioneering studies in the 1980s showed that a viscosity jump between the upper and lower mantle is required to reconcile these geoid and tomographically inferred density anomalies. These studies also predict 1-2 km of dynamic topography at the surface, dominated by degree 2. In contrast to this prediction, a global observational database of oceanic residual depth measurements indicates that degree 2 dynamic topography has peak amplitudes of only 500 m. Here, we attempt to reconcile observations of dynamic topography, geoid, gravity anomalies and CMB topography using instantaneous flow kernels. We exploit a density structure constructed from blended seismic tomography models, combining deep mantle imaging with higher resolution upper mantle features. Radial viscosity structure is discretised, and we invert for the best-fitting viscosity profile using a conjugate gradient search algorithm, subject to damping. Our results suggest that, due to strong sensitivity to radial viscosity structure, the Earth's geoid seems to be compatible with only ± 500 m of degree 2 dynamic topography.

Upper mantle anisotropy beneath Peru from SKS splitting: Constraints on flat slab dynamics and interaction with the Nazca Ridge

NASA Astrophysics Data System (ADS)

Eakin, Caroline M.; Long, Maureen D.; Wagner, Lara S.; Beck, Susan L.; Tavera, Hernando

2015-02-01

The Peruvian flat slab is by far the largest region of flat subduction in the world today, but aspects of its structure and dynamics remain poorly understood. In particular, questions remain over whether the relatively narrow Nazca Ridge subducting beneath southern Peru provides dynamic support for the flat slab or it is just a passive feature. We investigate the dynamics and interaction of the Nazca Ridge and the flat slab system by studying upper mantle seismic anisotropy across southern Peru. We analyze shear wave splitting of SKS, sSKS, and PKS phases at 49 stations distributed across the area, primarily from the PerU Lithosphere and Slab Experiment (PULSE). We observe distinct spatial variations in anisotropic structure along strike, most notably a sharp transition from coherent splitting in the north to pervasive null (non-split) arrivals in the south, with the transition coinciding with the northern limit of the Nazca Ridge. For both anisotropic domains there is evidence for complex and multi-layered anisotropy. To the north of the ridge our *KS splitting measurements likely reflect trench-normal mantle flow beneath the flat slab. This signal is then modified by shallower anisotropic layers, most likely in the supra-slab mantle, but also potentially from within the slab. To the south the sub-slab mantle is similarly anisotropic, with a trench-oblique fast direction, but widespread nulls appear to reflect dramatic heterogeneity in anisotropic structure above the flat slab. Overall the regional anisotropic structure, and thus the pattern of deformation, appears to be closely tied to the location of the Nazca Ridge, which further suggests that the ridge plays a key role in the mantle dynamics of the Peruvian flat slab system.

Constraints on dynamic topography from asymmetric subsidence of the mid-ocean ridges

NASA Astrophysics Data System (ADS)

Watkins, C. Evan; Conrad, Clinton P.

2018-02-01

Stresses from mantle convection deflect Earth's surface vertically, producing dynamic topography that is important for continental dynamics and sea-level change but difficult to observe due to overprinting by isostatic topography. For long wavelengths (∼104 km), the amplitude of dynamic topography is particularly uncertain, with mantle flow models typically suggesting larger amplitudes (>1000 m) than direct observations. Here we develop a new constraint on the amplitude of long-wavelength dynamic topography by examining asymmetries in seafloor bathymetry across mid-ocean ridges. We compare bathymetric profiles across the Mid-Atlantic Ridge (MAR) and the East Pacific Rise (EPR) and we find that the South American flank of both ridges subsides faster than its opposing flank. This pattern is consistent with dynamic subsidence across South America, supported by downwelling in the lower mantle. To constrain the amplitude of dynamic topography, we compare bathymetric profiles across both ridges after correcting bathymetry for several different models of dynamic topography with varying amplitudes and spatial patterns. We find that long-wavelength dynamic topography with an amplitude of only ∼500 m explains the observed asymmetry of the MAR. A similar model can explain EPR asymmetry but is complicated by additional asymmetrical topography associated with tectonic, crustal thickness, and/or asthenospheric temperature asymmetries across the EPR. After removing 500 m of dynamic topography, both the MAR and EPR exhibit a slower seafloor subsidence rate (∼280-290 m/Myr1/2) than previously reported. Our finding of only ∼500 m of long-wavelength dynamic topography may indicate the importance of thermochemical convection and/or large viscosity variations for lower mantle dynamics.

Modelling of Equilibrium Between Mantle and Core: Refractory, Volatile, and Highly Siderophile Elements

NASA Technical Reports Server (NTRS)

Righter, K.; Danielson, L.; Pando, K.; Shofner, G.; Lee, C. -T.

2013-01-01

Siderophile elements have been used to constrain conditions of core formation and differentiation for the Earth, Mars and other differentiated bodies [1]. Recent models for the Earth have concluded that the mantle and core did not fully equilibrate and the siderophile element contents of the mantle can only be explained under conditions where the oxygen fugacity changes from low to high during accretion and the mantle and core do not fully equilibrate [2,3]. However these conclusions go against several physical and chemical constraints. First, calculations suggest that even with the composition of accreting material changing from reduced to oxidized over time, the fO2 defined by metal-silicate equilibrium does not change substantially, only by approximately 1 logfO2 unit [4]. An increase of more than 2 logfO2 units in mantle oxidation are required in models of [2,3]. Secondly, calculations also show that metallic impacting material will become deformed and sheared during accretion to a large body, such that it becomes emulsified to a fine scale that allows equilibrium at nearly all conditions except for possibly the length scale for giant impacts [5] (contrary to conclusions of [6]). Using new data for D(Mo) metal/silicate at high pressures, together with updated partitioning expressions for many other elements, we will show that metal-silicate equilibrium across a long span of Earth s accretion history may explain the concentrations of many siderophile elements in Earth's mantle. The modeling includes refractory elements Ni, Co, Mo, and W, as well as highly siderophile elements Au, Pd and Pt, and volatile elements Cd, In, Bi, Sb, Ge and As.

Iron Spin Crossover in the New Hexagonal Aluminous (NAL) Phase

NASA Astrophysics Data System (ADS)

Hsu, H.

2017-12-01

The new hexagonal aluminous (NAL) phase, chemical formula AB2C6O12 (A = Na+, K+, Ca2+; B = Mg2+, Fe2+, Fe3+; C = Al3+, Si4+, Fe3+), is considered a major component ( 20 vol%) of mid-ocean ridge basalt (MORB) at lower-mantle conditions. Given that MORB can be transported back into the Earth's lower mantle via subduction, a thorough knowledge of the NAL phase is essential to fully understand the fate of subducted MORB and its role in mantle dynamics and heterogeneity. In this presentation, the complicated spin crossover of the Fe-bearing NAL phase will be discussed based on a series of first-principles calculations [1], in which the local density approximation + self-consistent Hubbard U (LDA+Usc) method was adopted. As revealed by these calculations, only the ferric iron (Fe3+) substituting Al/Si in the octahedral (C) site undergoes a crossover from the high-spin (HS) to the low-spin (LS) state at 40 GPa, while iron substituting Mg in the trigonal-prismatic (B) site remains in the HS state, regardless of its oxidation state (Fe2+ or Fe3+). The volume/elastic anomalies, iron nuclear quadrupole splittings, and crystal field spltting determined by calculations are in great agreement with experiments [2,3]. The calculations further predict that the HS-LS transition pressure of the NAL phase barely increases with temperature due to the three nearly degenerate LS states of Fe3+, suggesting that the elastic anomalies of this mineral can occur at the top lower mantle. [1] H. Hsu, Phys. Rev. B 95, 020406(R) (2017). [2] Y. Wu et al. Earth Planet. Sci. Lett. 434, 91-100 (2016). [3] S. S. Lobanov et al., J. Geophys. Res. Solid Earth 122, 3565 (2017).

Topographic asymmetry of the South Atlantic from global models of mantle flow and lithospheric stretching

NASA Astrophysics Data System (ADS)

Flament, Nicolas; Gurnis, Michael; Williams, Simon; Seton, Maria; Skogseid, Jakob; Heine, Christian; Dietmar Müller, R.

2014-02-01

The relief of the South Atlantic is characterized by elevated passive continental margins along southern Africa and eastern Brazil, and by the bathymetric asymmetry of the southern oceanic basin where the western flank is much deeper than the eastern flank. We investigate the origin of these topographic features in the present and over time since the Jurassic with a model of global mantle flow and lithospheric deformation. The model progressively assimilates plate kinematics, plate boundaries and lithospheric age derived from global tectonic reconstructions with deforming plates, and predicts the evolution of mantle temperature, continental crustal thickness, long-wavelength dynamic topography, and isostatic topography. Mantle viscosity and the kinematics of the opening of the South Atlantic are adjustable parameters in thirteen model cases. Model predictions are compared to observables both for the present-day and in the past. Present-day predictions are compared to topography, mantle tomography, and an estimate of residual topography. Predictions for the past are compared to tectonic subsidence from backstripped borehole data along the South American passive margin, and to dynamic uplift as constrained by thermochronology in southern Africa. Comparison between model predictions and observations suggests that the first-order features of the topography of the South Atlantic are due to long-wavelength dynamic topography, rather than to asthenospheric processes. The uplift of southern Africa is best reproduced with a lower mantle that is at least 40 times more viscous than the upper mantle.

Topographic asymmetry of the South Atlantic from global models of mantle flow and lithospheric stretching

NASA Astrophysics Data System (ADS)

Flament, Nicolas; Gurnis, Michael; Williams, Simon; Seton, Maria; Skogseid, Jakob; Heine, Christian; Müller, Dietmar

2014-05-01

The relief of the South Atlantic is characterized by elevated passive continental margins along southern Africa and eastern Brazil, and by the bathymetric asymmetry of the southern oceanic basin where the western flank is much deeper than the eastern flank. We investigate the origin of these topographic features in the present and over time since the Jurassic with a model of global mantle flow and lithospheric deformation. The model progressively assimilates plate kinematics, plate boundaries and lithospheric age derived from global tectonic reconstructions with deforming plates, and predicts the evolution of mantle temperature, continental crustal thickness, long-wavelength dynamic topography, and isostatic topography. Mantle viscosity and the kinematics of the opening of the South Atlantic are adjustable parameters in multiple model cases. Model predictions are compared to observables both for the present-day and in the past. Present-day predictions are compared to topography, mantle tomography, and an estimate of residual topography. Predictions for the past are compared to tectonic subsidence from backstripped borehole data along the South American passive margin, and to dynamic uplift as constrained by thermochronology in southern Africa. Comparison between model predictions and observations suggests that the first-order features of the topography of the South Atlantic are due to long-wavelength dynamic topography, rather than to asthenospheric processes. We find the uplift of southern Africa to be best reproduced with a lower mantle that is at least 40 times more viscous than the upper mantle.

Clustering of arc volcanoes caused by temperature perturbations in the back-arc mantle

PubMed Central

Lee, Changyeol; Wada, Ikuko

2017-01-01

Clustering of arc volcanoes in subduction zones indicates along-arc variation in the physical condition of the underlying mantle where majority of arc magmas are generated. The sub-arc mantle is brought in from the back-arc largely by slab-driven mantle wedge flow. Dynamic processes in the back-arc, such as small-scale mantle convection, are likely to cause lateral variations in the back-arc mantle temperature. Here we use a simple three-dimensional numerical model to quantify the effects of back-arc temperature perturbations on the mantle wedge flow pattern and sub-arc mantle temperature. Our model calculations show that relatively small temperature perturbations in the back-arc result in vigorous inflow of hotter mantle and subdued inflow of colder mantle beneath the arc due to the temperature dependence of the mantle viscosity. This causes a three-dimensional mantle flow pattern that amplifies the along-arc variations in the sub-arc mantle temperature, providing a simple mechanism for volcano clustering. PMID:28660880

Clustering of arc volcanoes caused by temperature perturbations in the back-arc mantle.

PubMed

Lee, Changyeol; Wada, Ikuko

2017-06-29

Clustering of arc volcanoes in subduction zones indicates along-arc variation in the physical condition of the underlying mantle where majority of arc magmas are generated. The sub-arc mantle is brought in from the back-arc largely by slab-driven mantle wedge flow. Dynamic processes in the back-arc, such as small-scale mantle convection, are likely to cause lateral variations in the back-arc mantle temperature. Here we use a simple three-dimensional numerical model to quantify the effects of back-arc temperature perturbations on the mantle wedge flow pattern and sub-arc mantle temperature. Our model calculations show that relatively small temperature perturbations in the back-arc result in vigorous inflow of hotter mantle and subdued inflow of colder mantle beneath the arc due to the temperature dependence of the mantle viscosity. This causes a three-dimensional mantle flow pattern that amplifies the along-arc variations in the sub-arc mantle temperature, providing a simple mechanism for volcano clustering.

Dynamics of cratons in an evolving mantle

NASA Astrophysics Data System (ADS)

O'Neill, C. J.; Lenardic, A.; Griffin, W. L.; O'Reilly, Suzanne Y.

2008-04-01

The tectonic quiescence of cratons on a tectonically active planet has been attributed to their physical properties such as buoyancy, viscosity, and yield strength. Previous modelling has shown the conditions under which cratons may be stable for the present, but cast doubt on how they survived in a more energetic mantle of the past. Here we incorporate an endothermic phase change at 670 km, and a depth-dependent viscosity structure consistent with post-glacial rebound and geoid modelling, to simulate the dynamics of cratons in an "Earth-like" convecting system. We find that cratons are unconditionally stable in such systems for plausible ranges of viscosity ratios between the root and asthenosphere (50-150) and the root/oceanic lithosphere yield strength ratio (5-30). Realistic mantle viscosity structures have limited effect on the average background cratonic stress state, but do buffer cratons from extreme stress excursions. An endothermic phase change at 670 km introduces an additional time-dependence into the system, with slab breakthrough into the lower mantle associated with 2-3 fold stress increases at the surface. Under Precambrian mantle conditions, however, the dominant effect is not more violent mantle avalanches, or faster mantle/plate velocities, but rather the drastic viscosity drop which results from hotter mantle conditions in the past. This results in a large decrease in the cratonic stress field, and promotes craton survival under the evolving mantle conditions of the early Earth.

On evolutionary climate tracks in deep mantle volatile cycle computed from numerical mantle convection simulations and its impact on the habitability of the Earth-like planets

NASA Astrophysics Data System (ADS)

Nakagawa, T.; Tajika, E.; Kadoya, S.

2017-12-01

Discussing an impact of evolution and dynamics in the Earth's deep interior on the surface climate change for the last few decades (see review by Ehlmann et al., 2016), the mantle volatile (particularly carbon) degassing in the mid-oceanic ridges seems to play a key role in understanding the evolutionary climate track for Earth-like planets (e.g. Kadoya and Tajika, 2015). However, since the mantle degassing occurs not only in the mid-oceanic ridges but also in the wedge mantle (island arc volcanism) and hotspots, to incorporate more accurate estimate of mantle degassing flux into the climate evolution framework, we developed a coupled model of surface climate-deep Earth evolution in numerical mantle convection simulations, including more accurate deep water and carbon cycle (e.g. Nakagawa and Spiegelman, 2017) with an energy balance theory of climate change. Modeling results suggest that the evolution of planetary climate computed from a developed model is basically consistent with an evolutionary climate track in simplified mantle degassing model (Kadoya and Tajika, 2015), but an occurrence timing of global (snowball) glaciation is strongly dependent on mantle degassing rate occurred with activities of surface plate motions. With this implication, the surface plate motion driven by deep mantle dynamics would play an important role in the planetary habitability of such as the Earth and Earth-like planets over geologic time-scale.

Estimating long-wavelength dynamic topographic change of passive continental margins since the Early Cretaceous

NASA Astrophysics Data System (ADS)

Müller, Dietmar; Hassan, Rakib; Gurnis, Michael; Flament, Nicolas; Williams, Simon

2017-04-01

The influence of mantle convection on dynamic topographic change along continental margins is difficult to unravel, because their stratigraphic record is dominated by tectonic subsidence caused by rifting. Yet, dynamic topography can potentially introduce significant depth anomalies along passive margins, influencing their water depth, sedimentary environments and geohistory. Here we follow a three-fold approach to estimate changes in dynamic topography along both continental interiors and passive margins based on a set of seven global mantle convection models. These models include different methodologies (forward and hybrid backward-forward methods), different plate reconstructions and alternative mantle rheologies. We demonstrate that a geodynamic forward model that includes adiabatic heating in addition to internal heating from radiogenic sources, and a mantle viscosity profile with a gradual increase in viscosity below the mantle transition zone, provides a greatly improved match to the spectral range of residual topography end-members as compared with previous models at very long wavelengths (spherical degrees 2-3). We combine global sea level estimates with predicted surface dynamic topography to evaluate the match between predicted continental flooding patterns and published paleo-coastlines by comparing predicted versus geologically reconstructed land fractions and spatial overlaps of flooded regions for individual continents since 140 Ma. Modelled versus geologically reconstructed land fractions match within 10% for most models, and the spatial overlaps of inundated regions are mostly between 85% and 100% for the Cenozoic, dropping to about 75-100% in the Cretaceous. We categorise the evolution of modelled dynamic topography in both continental interiors and along passive margins using cluster analysis to investigate how clusters of similar dynamic topography time series are distributed spatially. A subdivision of four clusters is found to best reveal end-members of dynamic topography evolution along passive margins and their hinterlands, differentiating topographic stability, long-term pronounced subsidence, initial stability over a dynamic high followed by moderate subsidence and regions that are relatively proximal to subduction zones with varied dynamic topography histories. Along passive continental margins the most commonly observed process is a gradual move from dynamic highs towards lows during the fragmentation of Pangea, reflecting that many passive margins now overly slabs sinking in the lower mantle. Our best-fit model results in up to 500 ±150 m of total dynamic subsidence of continental interiors while along passive margins the maximum predicted dynamic topographic change over 140 million years is about 350 ±150 m of subsidence. Models with plumes exhibit clusters of transient passive margin uplift of about 200 ±200m. The good overall match between predicted dynamic topography and geologically mapped paleo-coastlines makes a convincing case that mantle-driven topographic change is a critical component of relative sea level change, and one of the main driving forces generating the observed geometries and timings of large-scale shifts in paleo-coastlines.

Trench dynamics: Effects of dynamically migrating trench on subducting slab morphology and characteristics of subduction zones systems

NASA Astrophysics Data System (ADS)

Yoshida, Masaki

2017-07-01

Understanding the mechanisms of trench migration (retreat or advance) is crucial to characterizing the driving forces of Earth's tectonics plates, the origins of subducting slab morphologies in the deep mantle, and identifying the characteristics of subduction zones systems, which are among the fundamental issues of solid Earth science. A series of numerical simulations of mantle convection, focusing on plate subduction in a three-dimensional (3-D) regional spherical shell coordinate system, was performed to examine subduction zone characteristics, including geodynamic relationships among trench migration, back-arc stress, and slab morphology. The results show that a subducting slab tends to deflect around the base of the mantle transition zone and form a sub-horizontal slab because its front edge (its 'toe') is subject to resistance from the highly viscous lower mantle. As the sub-horizontal slab starts to penetrate into the lower mantle from its 'heel,' the toe of the slab is drawn into the lower mantle. The results for models with dynamically migrating trenches suggest that trench retreat is the dynamically self-consistent phenomenon in trench migration. The reason for this is that the strong lateral mantle flow that is generated as a sequence of events leading from corner flow at the subduction initiation to return flow of the formation of a sub-horizontal slab in the shallower part of mantle wedge produces the retreat of the subducting slab. In fact, a 'mantle suction force,' which is generated in the mantle wedge to fill space left by the retreating subducting plate, is enhanced by the subsequent trench retreat. Even when upwelling flow with significant positive buoyancy originates just above a mantle phase boundary at a depth of 410 km (as inferred from independent seismic tomographic, geodynamic, geochemical, and mineral physics), reaches the base of the overriding plate, and the overriding plate is slightly thinned, lithospheric stress tends to be compressed above the upwelling flow. The reason for this is that the strong lateral mantle flow originating from the upwelling flow generates resistance drag force at the base of the overriding plates. This situation may apply to a case of East Asia, under which the typical morphology of sub-horizontal slabs can be seen by seismic tomography. The strong lateral velocity observed in the shallower mantle wedge in the present numerical simulation may account for both the compressional subduction tectonics and back arc compression in the Japan-Kuril-Kamchatka, Aleutian, and South Chile trenches, as well as for weak plate-slab coupling, strong seismic coupling, and the possibility of great earthquakes along these trenches.

Mantle viscosity structure constrained by joint inversions of seismic velocities and density

NASA Astrophysics Data System (ADS)

Rudolph, M. L.; Moulik, P.; Lekic, V.

2017-12-01

The viscosity structure of Earth's deep mantle affects the thermal evolution of Earth, the ascent of mantle upwellings, sinking of subducted oceanic lithosphere, and the mixing of compositional heterogeneities in the mantle. Modeling the long-wavelength dynamic geoid allows us to constrain the radial viscosity profile of the mantle. Typically, in inversions for the mantle viscosity structure, wavespeed variations are mapped into density variations using a constant- or depth-dependent scaling factor. Here, we use a newly developed joint model of anisotropic Vs, Vp, density and transition zone topographies to generate a suite of solutions for the mantle viscosity structure directly from the seismologically constrained density structure. The density structure used to drive our forward models includes contributions from both thermal and compositional variations, including important contributions from compositionally dense material in the Large Low Velocity Provinces at the base of the mantle. These compositional variations have been neglected in the forward models used in most previous inversions and have the potential to significantly affect large-scale flow and thus the inferred viscosity structure. We use a transdimensional, hierarchical, Bayesian approach to solve the inverse problem, and our solutions for viscosity structure include an increase in viscosity below the base of the transition zone, in the shallow lower mantle. Using geoid dynamic response functions and an analysis of the correlation between the observed geoid and mantle structure, we demonstrate the underlying reason for this inference. Finally, we present a new family of solutions in which the data uncertainty is accounted for using covariance matrices associated with the mantle structure models.

Large-scale retreat and advance of shallow seas in Southeast Asia driven by mantle flow

NASA Astrophysics Data System (ADS)

Zahirovic, Sabin; Flament, Nicolas; Dietmar Müller, R.; Seton, Maria; Gurnis, Michael

2016-04-01

The Indonesian islands and surrounding region represent one of the most submerged, low-lying continental areas on Earth. Almost half of this region, known as Sundaland, is presently inundated by a shallow sea. The role of mantle convection in driving long-wavelength topography and vertical motion of the lithosphere in this region has largely been ignored when interpreting regional stratigraphic sections, despite a consensus that Southeast Asia presently situated on a "dynamic topography low" resulting from long-term post-Pangea subduction. However, dynamic topography is typically described as a temporally and spatially transient process, implying that Sundaland may have experienced significant vertical motions in the geological past, and thus must be considered when interpreting relative sea level changes and the paleogeographic indicators of advancing and retreating shallow seas. Although the present-day low regional elevation has been attributed to the massive volume of oceanic slabs sinking in the mantle beneath Southeast Asia, a Late Cretaceous to Eocene regional unconformity indicates that shallow seas retreated following regional flooding during the mid-Cretaceous sea level highstand. During the Eocene, less than one fifth of Sundaland was submerged, despite global sea level being ~200 m higher than at present. The regional nature of the switch from marine to terrestrial environments, that is out-of-sync with eustatic sea levels, suggests that broad mantle-driven dynamic uplift may have led to the emergence of Sundaland in the Late Cretaceous and Paleocene. We use numerical forward modelling of plate tectonics and mantle convection, and compare the predicted trends of dynamic topography with evidence from regional paleogeography and eustasy to determine the extent to which mantle-driven vertical motions of the lithosphere have influenced regional basin histories in Southeast Asia. A Late Cretaceous collision of Gondwana-derived terranes with Sundaland choked the active margin, leading to slab breakoff and a weakened mantle down-welling acting on the overriding plate, which resulted in regional dynamic uplift and emergence from a ~10-15 Myr-long subduction hiatus along the Sunda active margin. This explains the absence of sediment deposition across Sundaland and the emergence of Sundaland between ~80-60 Ma. Renewed subduction from ~60 Ma reinitiated dynamic subsidence of Sundaland, leading to submergence from ~40 Ma despite falling long-term global sea levels. Our results highlight a complete 'down-up-down' dynamic topography cycle experienced by Sundaland over 100 million years, with the transience of topography revealed in sedimentary basin stratigraphy punctuated with regional unconformities. Subduction-driven mantle convection models are now able to transform the geological record of basins into a dynamic surface history, enabling a deeper understanding of mechanisms that control landscape evolution across spatial and temporal scales.

Anomalous Subsidence at Rifted Continental Margins: Distinguishing Mantle Dynamic Topography from Anomalous Oceanic Crustal Thickness

NASA Astrophysics Data System (ADS)

Cowie, L.; Kusznir, N. J.

2012-12-01

It has been proposed that some continental rifted margins have anomalous subsidence histories and that at breakup they were elevated at shallower bathymetries than the isostatic response of classical rift models (McKenzie 1978) would predict. The existence of anomalous syn or post breakup subsidence of this form would have important implications for our understanding of the geodynamics of continental breakup and rifted continental margin formation, margin subsidence history and the evolution of syn and post breakup depositional systems. We have investigated three rifted continental margins; the Gulf of Aden, Galicia Bank and the Gulf of Lions, to determine whether the oceanic crust in the ocean-continent transition of these margins has present day anomalous subsidence and if so, whether it is caused by mantle dynamic topography or anomalous oceanic crustal thickness. Residual depth anomalies (RDA) corrected for sediment loading, using flexural backstripping and decompaction, have been calculated by comparing observed and age predicted oceanic bathymetries in order to identify anomalous oceanic bathymetry and subsidence at these margins. Age predicted bathymetric anomalies have been calculated using the thermal plate model predictions from Crosby & McKenzie (2009). Non-zero sediment corrected RDAs may result from anomalous oceanic crustal thickness with respect to the global average, or from mantle dynamic uplift. Positive RDAs may result from thicker than average oceanic crust or mantle dynamic uplift; negative RDAs may result from thinner than average oceanic crust or mantle dynamic subsidence. Gravity inversion incorporating a lithosphere thermal gravity anomaly correction and sediment thickness from 2D seismic data has been used to determine Moho depth and oceanic crustal basement thickness. The reference Moho depths used in the gravity inversion have been calibrated against seismic refraction Moho depths. The gravity inversion crustal basement thicknesses together with Airy isostasy have been used to predict a "synthetic" gravity derived RDA. Sediment corrected RDA for oceanic crust in the Gulf of Aden are positive (+750m) indicating anomalous uplift with respect to normal subsidence. Gravity inversion predicts normal thickness oceanic crust and a zero "synthetic" gravity derived RDA in the oceanic domain. The difference between the positive sediment corrected RDA and the zero "synthetic" gravity derived RDA, implies that the anomalous subsidence reported in the Gulf of Aden is the result of mantle dynamic uplift. For the oceanic crust outboard of Galicia Bank both the sediment corrected RDA and the "synthetic" gravity derived RDA are negative (-800m) and of similar magnitude, indicating anomalous subsidence, which is the result of anomalously thin oceanic crust, not mantle dynamic topography. We conclude that there is negligible mantle dynamic topography influencing the Galicia Bank region. In the Gulf of Lions, gravity inversion predicts thinner than average oceanic crust. Both sediment corrected RDA (-1km) and "synthetic" gravity derived RDA (-500m) are negative. The more negative sediment corrected RDA compared with the "synthetic" gravity derived RDA implies that the anomalous subsidence in the Gulf of Lions is the result of mantle dynamic subsidence as well as thinner than average oceanic crust.

The effect of ilmenite viscosity on the dynamics and evolution of an overturned lunar cumulate mantle

NASA Astrophysics Data System (ADS)

Zhang, Nan; Dygert, Nick; Liang, Yan; Parmentier, E. M.

2017-07-01

Lunar cumulate mantle overturn and the subsequent upwelling of overturned mantle cumulates provide a potential framework for understanding the first-order thermochemical evolution of the Moon. Upwelling of ilmenite-bearing cumulates (IBCs) after the overturn has a dominant influence on the dynamics and long-term thermal evolution of the lunar mantle. An important parameter determining the stability and convective behavior of the IBC is its viscosity, which was recently constrained through rock deformation experiments. To examine the effect of IBC viscosity on the upwelling of overturned lunar cumulate mantle, here we conduct three-dimensional mantle convection models with an evolving core superposed by an IBC-rich layer, which resulted from mantle overturn after magma ocean solidification. Our modeling shows that a reduction of mantle viscosity by 1 order of magnitude, due to the presence of ilmenite, can dramatically change convective planform and long-term lunar mantle evolution. Our model results suggest a relatively stable partially molten IBC layer that has surrounded the lunar core to the present day.