Can mantle convection be self-regulated?

PubMed Central

Korenaga, Jun

2016-01-01

The notion of self-regulating mantle convection, in which heat loss from the surface is constantly adjusted to follow internal radiogenic heat production, has been popular for the past six decades since Urey first advocated the idea. Thanks to its intuitive appeal, this notion has pervaded the solid earth sciences in various forms, but approach to a self-regulating state critically depends on the relation between the thermal adjustment rate and mantle temperature. I show that, if the effect of mantle melting on viscosity is taken into account, the adjustment rate cannot be sufficiently high to achieve self-regulation, regardless of the style of mantle convection. The evolution of terrestrial planets is thus likely to be far from thermal equilibrium and be sensitive to the peculiarities of their formation histories. Chance factors in planetary formation are suggested to become more important for the evolution of planets that are more massive than Earth. PMID:27551689

Numerical Mantle Convection Models of Crustal Formation in an Oceanic Environment in the Early Earth

NASA Astrophysics Data System (ADS)

van Thienen, P.; van den Berg, A. P.; Vlaar, N. J.

2001-12-01

The generation of basaltic crust in the early Earth by partial melting of mantle rocks, subject to investigation in this study, is thought to be a first step in the creation of proto-continents (consisting largely of felsic material), since partial melting of basaltic material was probably an important source for these more evolved rocks. In the early Archean the earth's upper mantle may have been hotter than today by as much as several hundred degrees centigrade. As a consequence, partial melting in shallow convective upwellings would have produced a layering of basaltic crust and underlying depleted (lherzolitic-harzburgitic) mantle peridotite which is much thicker than found under modern day oceanic ridges. When a basaltic crustal layer becomes sufficiently thick, a phase transition to eclogite may occur in the lower parts, which would cause delamination of this dense crustal layer and recycling of dense eclogite into the upper mantle. This recycling mechanism may have contributed significantly to the early cooling of the earth during the Archean (Vlaar et al., 1994). The delamination mechanism which limits the build-up of a thick basaltic crustal layer is switched off after sufficient cooling of the upper mantle has taken place. We present results of numerical modelling experiments of mantle convection including pressure release partial melting. The model includes a simple approximate melt segregation mechanism and basalt to eclogite phase transition, to account for the dynamic accumulation and recycling of the crust in an upper mantle subject to secular cooling. Finite element methods are used to solve for the viscous flow field and the temperature field, and lagrangian particle tracers are used to represent the evolving composition due to partial melting and accumulation of the basaltic crust. We find that this mechanism creates a basaltic crust of several tens of kilometers thickness in several hundreds of million years. This is accompanied by a cooling of

Facilitating atmosphere oxidation through mantle convection

NASA Astrophysics Data System (ADS)

Lee, K. K. M.; Gu, T.; Creasy, N.; Li, M.; McCammon, C. A.; Girard, J.

2017-12-01

Earth's mantle connects the surface with the deep interior through convection, and the evolution of its redox state will affect the distribution of siderophile elements, recycling of refractory isotopes, and the oxidation state of the atmosphere through volcanic outgassing. While the rise of oxygen in the atmosphere, i.e., the Great Oxidation Event (GOE) occurred 2.4 billion years ago (Ga), multiple lines of evidence point to oxygen production in the atmosphere well before 2.4 Ga. In contrast to the fluctuations of atmospheric oxygen, vanadium in Archean mantle lithosphere suggests that the mantle redox state has been constant for 3.5 Ga. Indeed, the connection between the redox state of the deep Earth and the atmosphere is enigmatic as is the effect of redox state on mantle dynamics. Here we show a redox-induced density contrast affects mantle convection and may potentially cause the oxidation of the upper mantle. We compressed two synthetic enstatite chondritic samples with identical bulk compositions but formed under different oxygen fugacities (fO2) to lower mantle pressures and temperatures and find Al2O3 forms its own phase separate from the dominant bridgmanite phase in the more reduced composition, in contrast to a more Al-rich, bridgmanite-dominated assemblage for a more oxidized starting composition. As a result, the reduced material is 1-1.5% denser than the oxidized material. Subsequent experiments on other plausible mantle compositions, which differ only in redox state of the starting glass materials, show similar results: distinct mineral assemblages and density contrasts up to 4%. Our geodynamic simulations suggest that such a density contrast causes a rapid ascent and accumulation of oxidized material in the upper mantle, with descent of the denser reduced material to the core-mantle boundary. The resulting heterogeneous redox conditions in Earth's interior may have contributed to the large low-shear velocity provinces in the lower mantle and the

Mantle convection and the distribution of geochemical reservoirs in the silicate shell of the Earth

NASA Astrophysics Data System (ADS)

Walzer, Uwe; Hendel, Roland

2010-05-01

We present a dynamic 3-D spherical-shell model of mantle convection and the evolution of the chemical reservoirs of the Earth`s silicate shell. Chemical differentiation, convection, stirring and thermal evolution constitute an inseparable dynamic system. Our model is based on the solution of the balance equations of mass, momentum, energy, angular momentum, and four sums of the number of atoms of the pairs 238U-206Pb, 235U-207Pb, 232Th-208Pb, and 40K-40Ar. Similar to the present model, the continental crust of the real Earth was not produced entirely at the start of the evolution but developed episodically in batches [1-7]. The details of the continental distribution of the model are largely stochastic, but the spectral properties are quite similar to the present real Earth. The calculated Figures reveal that the modeled present-day mantle has no chemical stratification but we find a marble-cake structure. If we compare the observational results of the present-day proportion of depleted MORB mantle with the model then we find a similar order of magnitude. The MORB source dominates under the lithosphere. In our model, there are nowhere pure unblended reservoirs in the mantle. It is, however, remarkable that, in spite of 4500 Ma of solid-state mantle convection, certain strong concentrations of distributed chemical reservoirs continue to persist in certain volumes, although without sharp abundance boundaries. We deal with the question of predictable and stochastic portions of the phenomena. Although the convective flow patterns and the chemical differentiation of oceanic plateaus are coupled, the evolution of time-dependent Rayleigh number, Rat , is relatively well predictable and the stochastic parts of the Rat(t)-curves are small. Regarding the juvenile growth rates of the total mass of the continents, predictions are possible only in the first epoch of the evolution. Later on, the distribution of the continental-growth episodes is increasingly stochastic

Water in geodynamical models of mantle convection and plate tectonics

NASA Astrophysics Data System (ADS)

Rodríguez-González, J.; Van Hunen, J.; Chotalia, K.; Lithgow-Bertelloni, C. R.; Rozel, A.; Tackley, P. J.; Nakagawa, T.

2017-12-01

The presence of water in the the mantle has a significant effect in the dynamical and thermal evolution of Earth, which partially explains the differences with other planets and is a key factor for the presence of life on Earth. First, a small amount of water can decrease the mantle viscosity by a several orders of magnitude, thereby changing the convection regime and affecting the thermal evolution. Second, the presence of water significantly changes the solidus curve, with crucial implications for melting. Third, water in the mantle can change the Clapeyron slope of mantle materials, which changes the depth at which phase transitions take place. The thermal and dynamical evolution of Earth under the presence of water in the mantle has been the focus of recent studies, but many questions remain unanswered. In this project we intend to investigate how the maximum water capacity of different mantle regions affects water transport and Earth's convective regime. We will study the effect phase transitions under the presence of water, which can change the buoyancy of slabs in the transition zone. We present preliminary results numerical models of global mantle convection for the whole history of earth using the numerical geodynamics software tool StagYY. We will use a new parametrisation of dehydration processes, obtained from high-resolution numerical simulations, to implement a more accurate description of the water released from the slab as it travels through the mantle. We have integrated recent experimental results of the water capacity of deep mantle minerals to study the water circulation and the total water budget. We use data from the most recent experiments and ab-inito calculations to implement a realistic rheology.

Osmium isotopes and mantle convection.

PubMed

Hauri, Erik H

2002-11-15

The decay of (187)Re to (187)Os (with a half-life of 42 billion years) provides a unique isotopic fingerprint for tracing the evolution of crustal materials and mantle residues in the convecting mantle. Ancient subcontinental mantle lithosphere has uniquely low Re/Os and (187)Os/(188)Os ratios due to large-degree melt extraction, recording ancient melt-depletion events as old as 3.2 billion years. Partial melts have Re/Os ratios that are orders of magnitude higher than their sources, and the subduction of oceanic or continental crust introduces into the mantle materials that rapidly accumulate radiogenic (187)Os. Eclogites from the subcontinental lithosphere have extremely high (187)Os/(188)Os ratios, and record ages as old as the oldest peridotites. The data show a near-perfect partitioning of Re/Os and (187)Os/(188)Os ratios between peridotites (low) and eclogites (high). The convecting mantle retains a degree of Os-isotopic heterogeneity similar to the lithospheric mantle, although its amplitude is modulated by convective mixing. Abyssal peridotites from the ocean ridges have low Os isotope ratios, indicating that the upper mantle had undergone episodes of melt depletion prior to the most recent melting events to produce mid-ocean-ridge basalt. The amount of rhenium estimated to be depleted from the upper mantle is 10 times greater than the rhenium budget of the continental crust, requiring a separate reservoir to close the mass balance. A reservoir consisting of 5-10% of the mantle with a rhenium concentration similar to mid-ocean-ridge basalt would balance the rhenium depletion of the upper mantle. This reservoir most likely consists of mafic oceanic crust recycled into the mantle over Earth's history and provides the material that melts at oceanic hotspots to produce ocean-island basalts (OIBs). The ubiquity of high Os isotope ratios in OIB, coupled with other geochemical tracers, indicates that the mantle sources of hotspots contain significant quantities

Long-term preservation of early formed mantle heterogeneity by mobile lid convection: Importance of grainsize evolution

NASA Astrophysics Data System (ADS)

Foley, Bradford J.; Rizo, Hanika

2017-10-01

The style of tectonics on the Hadean and Archean Earth, particularly whether plate tectonics was in operation or not, is debated. One important, albeit indirect, constraint on early Earth tectonics comes from observations of early-formed geochemical heterogeneities: 142Nd and 182W anomalies recorded in Hadean to Phanerozoic rocks from different localities indicate that chemically heterogeneous reservoirs, formed during the first ∼500 Myrs of Earth's history, survived their remixing into the mantle for over 1 Gyrs. Such a long mixing time is difficult to explain because hotter mantle temperatures, expected for the early Earth, act to lower mantle viscosity and increase convective vigor. Previous studies found that mobile lid convection typically erases heterogeneity within ∼100 Myrs under such conditions, leading to the hypothesis that stagnant lid convection on the early Earth was responsible for the observed long mixing times. However, using two-dimensional Cartesian convection models that include grainsize evolution, we find that mobile lid convection can preserve heterogeneity at high mantle temperature conditions for much longer than previously thought, because higher mantle temperatures lead to larger grainsizes in the lithosphere. These larger grainsizes result in stronger plate boundaries that act to slow down surface and interior convective motions, in competition with the direct effect temperature has on mantle viscosity. Our models indicate that mobile lid convection can preserve heterogeneity for ≈0.4-1 Gyrs at early Earth mantle temperatures when the initial heterogeneity has the same viscosity as the background mantle, and ≈1-4 Gyrs when the heterogeneity is ten times more viscous than the background mantle. Thus, stagnant lid convection is not required to explain long-term survival of early formed geochemical heterogeneities, though these heterogeneities having an elevated viscosity compared to the surrounding mantle may be essential for their

Whole-mantle convection with tectonic plates preserves long-term global patterns of upper mantle geochemistry.

PubMed

Barry, T L; Davies, J H; Wolstencroft, M; Millar, I L; Zhao, Z; Jian, P; Safonova, I; Price, M

2017-05-12

The evolution of the planetary interior during plate tectonics is controlled by slow convection within the mantle. Global-scale geochemical differences across the upper mantle are known, but how they are preserved during convection has not been adequately explained. We demonstrate that the geographic patterns of chemical variations around the Earth's mantle endure as a direct result of whole-mantle convection within largely isolated cells defined by subducting plates. New 3D spherical numerical models embedded with the latest geological paleo-tectonic reconstructions and ground-truthed with new Hf-Nd isotope data, suggest that uppermost mantle at one location (e.g. under Indian Ocean) circulates down to the core-mantle boundary (CMB), but returns within ≥100 Myrs via large-scale convection to its approximate starting location. Modelled tracers pool at the CMB but do not disperse ubiquitously around it. Similarly, mantle beneath the Pacific does not spread to surrounding regions of the planet. The models fit global patterns of isotope data and may explain features such as the DUPAL anomaly and long-standing differences between Indian and Pacific Ocean crust. Indeed, the geochemical data suggests this mode of convection could have influenced the evolution of mantle composition since 550 Ma and potentially since the onset of plate tectonics.

Towards driving mantle convection by mineral physics

NASA Astrophysics Data System (ADS)

Piazzoni, A. S.; Bunge, H.; Steinle-Neumann, G.

2005-12-01

Models of mantle convection have become increasingly sophisticated over the past decade, accounting, for example, for 3 D spherical geometry, and changes in mantle rheology due to variations in temperature and stress. In light of such advances it is surprising that growing constraints on mantle structure derived from mineral physics have not yet been fully brought to bear on mantle convection models. In fact, despite much progress in our understanding of mantle mineralogy a partial description of the equation of state is often used to relate density changes to pressure and temperature alone, without taking into account compositional and mineralogical models of the mantle. Similarly, for phase transitions an incomplete description of thermodynamic constraints is often used, resulting in significant uncertainties in model behavior. While a number of thermodynamic models (some with limited scope) have been constructed recently, some lack the rigor in thermodynamics - for example with respect to the treatment of solid solution - that is needed to make predictions about mantle structure. Here we have constructed a new thermodynamic database for the mantle and have coupled the resulting density dynamically with mantle convection models. The database is build on a self-consistent Gibb's free energy minimization of the system MgO-FeO-SiO2-CaO-Al2O3 that is appropriate for standard (dry) chemical models of the Earth's mantle for relevant high pressure and temperature phases. We have interfaced the database with a high-resolution 2-D convection code (2DTERRA), dynamically coupling the thermodynamic model (density) with the conservation equations of mantle flow. The coupled model is run for different parameterizations of viscosity, initial temperature conditions, and varying the internal vs. external heating. We compare the resulting flow and temperature fields to cases with the Boussinesq approximation and other classical descriptions of the equation of state in mantle

Thermal Evolution of Earth's Mantle During the Accretion

NASA Astrophysics Data System (ADS)

Arkani-Hamed, J.; Roberts, J. H.

2017-12-01

Earth is likely formed by accreting Moon to Mars size embryos. The impact heating by an embryo melts the embryo and the upper mantle of the Earth beneath the impact site. The iron core of the embryo sinks and merges with the core of the Earth, while the mantle of the embryo mixes with the upper mantle of the Earth, producing a buoyant molten/partially molten magma pond. Strong but localized mantle dynamics results in fast lithostatic adjustment that pours out a huge amount of molten and partially molten magma which spread on the Earth, and together with impact ejecta creates a globe encircling magma ocean. The lithostatic adjustment diminishes as the magma ocean becomes globe encircling within 104 to 105 yr. The major part of the thermal evolution of Earth's mantle after an impact takes place in the presence of a thick and hot magma ocean, which hampers heat loss from the mantle and suppresses global mantle dynamics. Because the impact velocity of an embryo increases as the Earth grows, a given magma ocean is hotter than the previous ones. We investigated this scenario using 25 Moon to Mars size embryos. Due to random geographic impact sites we considered vertical impacts since no information is available about the impact angles. This may over estimate the impact heating by a factor of 1.4 with respect to the most probable impact angle of 45o. The thermal structure of the Earth at the end of accretion is layered, aside from the localized magma ponds that are distributed randomly due to the random geographic impact sites. We also take into account the impact heating of the solid lower mantle, the heating of the lower mantle by the gravitational energy released through sinking of an embryo's core. We then follow the thermal evolution of the mantle of a growing Earth using a 3D convection model. The Earth grows due to merging of the impactor iron core with the Earth's core, and the accumulating magma ocean on the surface. The growth enhances the lithostatic pressure

Global-scale water circulation in the Earth's mantle: Implications for the mantle water budget in the early Earth

NASA Astrophysics Data System (ADS)

Nakagawa, Takashi; Spiegelman, Marc W.

2017-04-01

We investigate the influence of the mantle water content in the early Earth on that in the present mantle using numerical convection simulations that include three processes for redistribution of water: dehydration, partitioning of water into partially molten mantle, and regassing assuming an infinite water reservoir at the surface. These models suggest that the water content of the present mantle is insensitive to that of the early Earth. The initial water stored during planetary formation is regulated up to 1.2 OMs (OM = Ocean Mass; 1.4 ×1021 kg), which is reasonable for early Earth. However, the mantle water content is sensitive to the rheological dependence on the water content and can range from 1.2 to 3 OMs at the present day. To explain the evolution of mantle water content, we computed water fluxes due to subducting plates (regassing), degassing and dehydration. For weakly water dependent viscosity, the net water flux is almost balanced with those three fluxes but, for strongly water dependent viscosity, the regassing dominates the water cycle system because the surface plate activity is more vigorous. The increased convection is due to enhanced lubrication of the plates caused by a weak hydrous crust for strongly water dependent viscosity. The degassing history is insensitive to the initial water content of the early Earth as well as rheological strength. The degassing flux from Earth's surface is calculated to be approximately O (1013) kg /yr, consistent with a coupled model of climate evolution and mantle thermal evolution.

Probing Seismically Melting Induced Mantle Heterogeneities in Thermal-chemical Convection Models

NASA Astrophysics Data System (ADS)

Heck, H. V.; Davies, H.; Nowacki, A.; Wookey, J. M.

2015-12-01

Two regions at the base of the Earth's mantle (the Large Low-Shear Velocity Provinces) pose a fundamental problem in understanding large-scale mantle dynamics and history. Are they dense piles of (possibly primordial) material separated from mantle circulation, or large-scale thermal features which are part of global mantle convection? Or some combination of the two? We use our numerical 3D spherical mantle convection code to perform simulations of the Earths mantle dynamical evolution. We drive the surface velocity of the model according to 200 Ma plate motion reconstructions, to arrive at Earth-like structures in the mantle at present day. Variations in bulk chemistry will be tracked in two ways: 1) by starting the calculations with a (primordial) dense layer at the base of the mantle, and 2) by tracking basalt fraction which is fractionated upon melting close to the surface. The resulting distribution of chemical heterogeneity and temperature will be converted to seismic velocities. This will be done with a thermodynamical database (Stixrude & Lithgow-Bertelloni, GJI, 2005, 2011), allowing us to compare the model with previous observations of triplications and waveform complexity near the margins of the LLSVPs. These observations have been taken as proof that strong chemical variations are present; our simulations can be used to show whether this is true, or if purely thermal convection can also cause these features. We simulate finite-frequency, 3D seismograms at ~5 s period and compare these with previous studies.

Mantle convection and plate tectonics: toward an integrated physical and chemical theory

PubMed

Tackley

2000-06-16

Plate tectonics and convection of the solid, rocky mantle are responsible for transporting heat out of Earth. However, the physics of plate tectonics is poorly understood; other planets do not exhibit it. Recent seismic evidence for convection and mixing throughout the mantle seems at odds with the chemical composition of erupted magmas requiring the presence of several chemically distinct reservoirs within the mantle. There has been rapid progress on these two problems, with the emergence of the first self-consistent models of plate tectonics and mantle convection, along with new geochemical models that may be consistent with seismic and dynamical constraints on mantle structure.

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.

Two-Body Convection in the Mantle of the Earth: E/W Asymmetry, Under Astronomically Determined Tilt in g

NASA Astrophysics Data System (ADS)

Bostrom, R. C.

2002-12-01

Under purely geocentric gravity, over time displacement under mantle convection is globally symmetrical, resulting in zero net lithosphere rotation. The effect is here explored of substituting the asymmetric Earth-Moon field, gconv, prevalent in actuality. The gravity responsible for mantle convection is defined as the vector sum of a vertical component and the day-averaged attraction of masses lagging tidal equilibrium. The increasingly accurately measured lunar recession may then be used to delimit the internal field in terms of the secular luni-tidal interval of the Earth as a whole, some 600 seconds [1], without having to identify tidal components i.e. separate marine from body tides. In context the astronomic phase-lag may be viewed as a global isostatic anomaly, in which the longitude circles marking Earth's gravimetric figure are located east of those describing its perpetually unattained equilibrium figure by some 89 km at the Equator. Reference the hydrostatic ellipsoid gconv is tilted by the astronomically delimited amount, albeit that the phase lag is attributable in part to the convection itself. As with the convection, the tectonic significance of its asymmetry is determinable geodetically. Using present art-state a strategically located GPS grid [2] would provide continuously more precise separation of the asymmetric component of surface displacement. In developing plate-motion models including members of the Nuvel series, it would be logical to follow up rather than discard the set permitting minor asymmetrical convection sans net torque, such as an element of net-lithosphere-rotation relative to plumes. To conserve system angular-momentum, this may be the only valid set. Characteristics of the convection to be expected accord with 'paradoxical' features of plate tectonics under purely radial gravity, including: difficulty in closing plate-motion circuits; net-lithosphere-rotation refce. hot-spots, sans net torque; geotectonic maps ranging from

Tectonic predictions with mantle convection models

NASA Astrophysics Data System (ADS)

Coltice, Nicolas; Shephard, Grace E.

2018-04-01

Over the past 15 yr, numerical models of convection in Earth's mantle have made a leap forward: they can now produce self-consistent plate-like behaviour at the surface together with deep mantle circulation. These digital tools provide a new window into the intimate connections between plate tectonics and mantle dynamics, and can therefore be used for tectonic predictions, in principle. This contribution explores this assumption. First, initial conditions at 30, 20, 10 and 0 Ma are generated by driving a convective flow with imposed plate velocities at the surface. We then compute instantaneous mantle flows in response to the guessed temperature fields without imposing any boundary conditions. Plate boundaries self-consistently emerge at correct locations with respect to reconstructions, except for small plates close to subduction zones. As already observed for other types of instantaneous flow calculations, the structure of the top boundary layer and upper-mantle slab is the dominant character that leads to accurate predictions of surface velocities. Perturbations of the rheological parameters have little impact on the resulting surface velocities. We then compute fully dynamic model evolution from 30 and 10 to 0 Ma, without imposing plate boundaries or plate velocities. Contrary to instantaneous calculations, errors in kinematic predictions are substantial, although the plate layout and kinematics in several areas remain consistent with the expectations for the Earth. For these calculations, varying the rheological parameters makes a difference for plate boundary evolution. Also, identified errors in initial conditions contribute to first-order kinematic errors. This experiment shows that the tectonic predictions of dynamic models over 10 My are highly sensitive to uncertainties of rheological parameters and initial temperature field in comparison to instantaneous flow calculations. Indeed, the initial conditions and the rheological parameters can be good enough

Plate Tectonic Cycling and Whole Mantle Convection Modulate Earth's 3He/22Ne Ratio

NASA Astrophysics Data System (ADS)

Dygert, N. J.; Jackson, C.; Hesse, M. A.; Tremblay, M. M.; Shuster, D. L.; Gu, J.

2016-12-01

3He and 22Ne are not produced in the mantle or fractionated by partial melting, and neither isotope is recycled back into the mantle by subduction of oceanic basalt or sediment. Thus, it is a surprise that large 3He/22Ne variations exist within the mantle and that the mantle has a net elevated 3He/22Ne ratio compared to volatile-rich planetary precursor materials. Depleted subcontinental lithospheric mantle and mid-ocean ridge basalt (MORB) mantle have distinctly higher 3He/22Ne compared to ocean island basalt (OIB) sources ( 4-12.5 vs. 2.5-4.5, respectively) [1,2]. The low 3He/22Ne of OIBs approaches chondritic ( 1) and solar nebula values ( 1.5). The high 3He/22Ne of the MORB mantle is not similar to solar sources or any known family of meteorites, requiring a mechanism for fractionating He from Ne in the mantle and suggesting isolation of distinct mantle reservoirs throughout geologic time. We model the formation of a MORB source with elevated and variable 3He/22Ne though diffusive exchange between dunite channel-hosted basaltic liquids and harzburgite wallrock beneath mid-ocean ridges. Over timescales relevant to mantle upwelling beneath spreading centers, He may diffuse tens to hundreds of meters into wallrock while Ne is relatively immobile, producing a regassed, depleted mantle lithosphere with elevated 3He/22Ne. Subduction of high 3He/22Ne mantle would generate a MORB source with high 3He/22Ne. Regassed, high 3He/22Ne mantle lithosphere has He concentrations 2-3 orders of magnitude lower than undegassed mantle. To preserve the large volumes of high 3He/22Ne mantle required by the MORB source, mixing between subducted and undegassed mantle reservoirs must have been limited throughout geologic time. Using the new 3He/22Ne constraints, we ran a model similar to [3] to quantify mantle mixing timescales, finding they are on the order of Gyr assuming physically reasonable seafloor spreading rates, and that Earth's convecting mantle has lost >99% of its primordial

Large-scale compositional heterogeneity in the Earth's mantle

NASA Astrophysics Data System (ADS)

Ballmer, M.

2017-12-01

Seismic imaging of subducted Farallon and Tethys lithosphere in the lower mantle has been taken as evidence for whole-mantle convection, and efficient mantle mixing. However, cosmochemical constraints point to a lower-mantle composition that has a lower Mg/Si compared to upper-mantle pyrolite. Moreover, geochemical signatures of magmatic rocks indicate the long-term persistence of primordial reservoirs somewhere in the mantle. In this presentation, I establish geodynamic mechanisms for sustaining large-scale (primordial) heterogeneity in the Earth's mantle using numerical models. Mantle flow is controlled by rock density and viscosity. Variations in intrinsic rock density, such as due to heterogeneity in basalt or iron content, can induce layering or partial layering in the mantle. Layering can be sustained in the presence of persistent whole mantle convection due to active "unmixing" of heterogeneity in low-viscosity domains, e.g. in the transition zone or near the core-mantle boundary [1]. On the other hand, lateral variations in intrinsic rock viscosity, such as due to heterogeneity in Mg/Si, can strongly affect the mixing timescales of the mantle. In the extreme case, intrinsically strong rocks may remain unmixed through the age of the Earth, and persist as large-scale domains in the mid-mantle due to focusing of deformation along weak conveyor belts [2]. That large-scale lateral heterogeneity and/or layering can persist in the presence of whole-mantle convection can explain the stagnation of some slabs, as well as the deflection of some plumes, in the mid-mantle. These findings indeed motivate new seismic studies for rigorous testing of model predictions. [1] Ballmer, M. D., N. C. Schmerr, T. Nakagawa, and J. Ritsema (2015), Science Advances, doi:10.1126/sciadv.1500815. [2] Ballmer, M. D., C. Houser, J. W. Hernlund, R. Wentzcovitch, and K. Hirose (2017), Nature Geoscience, doi:10.1038/ngeo2898.

Archean greenstone-tonalite duality: Thermochemical mantle convection models or plate tectonics in the early Earth global dynamics?

NASA Astrophysics Data System (ADS)

Kerrich, Robert; Polat, Ali

2006-03-01

Mantle convection and plate tectonics are one system, because oceanic plates are cold upper thermal boundary layers of the convection cells. As a corollary, Phanerozoic-style of plate tectonics or more likely a different version of it (i.e. a larger number of slowly moving plates, or similar number of faster plates) is expected to have operated in the hotter, vigorously convecting early Earth. Despite the recent advances in understanding the origin of Archean greenstone-granitoid terranes, the question regarding the operation of plate tectonics in the early Earth remains still controversial. Numerical model outputs for the Archean Earth range from predominantly shallow to flat subduction between 4.0 and 2.5 Ga and well-established steep subduction since 2.5 Ga [Abbott, D., Drury, R., Smith, W.H.F., 1994. Flat to steep transition in subduction style. Geology 22, 937-940], to no plate tectonics but rather foundering of 1000 km sectors of basaltic crust, then "resurfaced" by upper asthenospheric mantle basaltic melts that generate the observed duality of basalts and tonalities [van Thienen, P., van den Berg, A.P., Vlaar, N.J., 2004a. Production and recycling of oceanic crust in the early earth. Tectonophysics 386, 41-65; van Thienen, P., Van den Berg, A.P., Vlaar, N.J., 2004b. On the formation of continental silicic melts in thermochemical mantle convection models: implications for early Earth. Tectonophysics 394, 111-124]. These model outputs can be tested against the geological record. Greenstone belt volcanics are composites of komatiite-basalt plateau sequences erupted from deep mantle plumes and bimodal basalt-dacite sequences having the geochemical signatures of convergent margins; i.e. horizontally imbricated plateau and island arc crust. Greenstone belts from 3.8 to 2.5 Ga include volcanic types reported from Cenozoic convergent margins including: boninites; arc picrites; and the association of adakites-Mg andesites- and Nb-enriched basalts. Archean cratons

Evolution of the earliest mantle caused by the magmatism-mantle upwelling feedback: Implications for the Moon and the Earth

NASA Astrophysics Data System (ADS)

Ogawa, M.

2017-12-01

The two most important agents that cause mantle evolution are magmatism and mantle convection. My earlier 2D numerical models of a coupled magmatism-mantle convection system show that these two agents strongly couple each other, when the Rayleigh number Ra is sufficiently high: magmatism induced by a mantle upwelling flow boosts the upwelling flow itself. The mantle convection enhanced by this positive feedback (the magmatism-mantle upwelling, or MMU, feedback) causes vigorous magmatism and, at the same time, strongly stirs the mantle. I explored how the MMU feedback influences the evolution of the earliest mantle that contains the magma ocean, based on a numerical model where the mantle is hot and its topmost 1/3 is partially molten at the beginning of the calculation: The evolution drastically changes its style, as Ra exceeds the threshold for onset of the MMU feedback, around 107. At Ra < 107, basaltic materials generated by the initial widespread magmatism accumulate in the deep mantle to form a layer; the basaltic layer is colder than the overlying shallow mantle. At Ra > 107, however, the mantle remains compositionally more homogeneous in spite of the widespread magmatism, and the deep mantle remains hotter than the shallow mantle, because of the strong convective stirring caused by the feedback. The threshold value suggests that the mantle of a planet larger than Mars evolves in a way substantially different from that in the Moon does. Indeed, in my earlier models, magmatism makes the early mantle compositionally stratified in the Moon, but the effects of strong convective stirring overwhelms that of magmatism to keep the mantle compositionally rather homogeneous in Venus and the Earth. The MMU feedback is likely to be a key to understanding why vestiges of the magma ocean are so scarce in the Earth.

Constraints from Earth's heat budget on mantle dynamics

NASA Astrophysics Data System (ADS)

Kellogg, L. H.; Ferrachat, S.

2006-12-01

Recent years have seen an increase in the number of proposed models to explain Earth's mantle dynamics: while two end-members, pure layered convection with the upper and lower mantle convecting separately from each other, and pure, whole mantle convection, appear not to satisfy all the observations, several addition models have been proposed. These models include and attempt to characterize least one reservoir that is enriched in radiogenic elements relative to the mid-ocean ridge basalt (MORB) source, as is required to account for most current estimates of the Earth's heat budget. This reservoir would also be responsible for the geochemical signature in some ocean island basalts (OIBs) like Hawaii, but must be rarely sampled at the surface. Our current knowledge of the mass- and heat-budget for the bulk silicate Earth from geochemical, cosmochemical and geodynamical observations and constraints enables us to quantify the radiogenic heat enrichment required to balance the heat budget. Without assuming any particular model for the structure of the reservoir, we first determine the inherent trade-off between heat production rate and mass of the reservoir. Using these constraints, we then investigate the dynamical inferences of the heat budget, assuming that the additional heat is produced within a deep layer above the core-mantle boundary. We carry out dynamical models of layered convection using four different fixed reservoir volumes, corresponding to deep layers of thicknesses 150, 500 1000 and 1600 km, respectively, and including both temperature-dependent viscosity and an instrinsic viscosity jump between upper and lower mantle. We then assess the viability of these cases against 5 criteria: stability of the deep layer through time, topography of the interface, effective density profile, intrinsic chemical density and the heat flux at the CMB.

Investigating melting induced mantle heterogeneities in plate driven mantle convection models

NASA Astrophysics Data System (ADS)

Price, M.; Davies, H.; Panton, J.

2017-12-01

Observations from geochemistry and seismology continue to suggest a range of complex heterogeneity in Earth's mantle. In the deep mantle, two large low velocity provinces (LLVPs) have been regularly observed in seismic studies, with their longevity, composition and density compared to the surrounding mantle debated. The cause of these observed LLVPs is equally uncertain, with previous studies advocating either thermal or thermo-chemical causes. There is also evidence that these structures could provide chemically distinct reservoirs within the mantle, with recent studies also suggesting there may be additional reservoirs in the mantle, such as bridgmanite-enriched ancient mantle structures (BEAMS). One way to test these hypotheses is using computational models of the mantle, with models that capture the full 3D system being both complex and computationally expensive. Here we present results from our global mantle model TERRA. Using our model, we can track compositional variations in the convecting mantle that are generated by self-consistent, evolving melting zones. Alongside the melting, we track trace elements and other volatiles which can be partitioned during melting events, and expelled and recycled at the surface. Utilising plate reconstruction models as a boundary condition, the models generate the tectonic features observed at Earth's surface, while also organising the lower mantle into recognisable degree-two structures. This results in our models generating basaltic `oceanic' crusts which are then brought into the mantle at tectonic boundaries, providing additional chemical heterogeneity in the mantle volume. Finally, by utilising thermodynamic lookup tables to convert the final outputs from the model to seismic structures, together with resolution filters for global tomography models, we are able to make direct comparisons between our results and observations. By varying the parameters of the model, we investigate a range of current hypotheses for

Supercontinent cycles, true polar wander, and very long-wavelength mantle convection

NASA Astrophysics Data System (ADS)

Zhong, Shijie; Zhang, Nan; Li, Zheng-Xiang; Roberts, James H.

2007-09-01

We show in this paper that mobile-lid mantle convection in a three-dimensional spherical shell with observationally constrained mantle viscosity structure, and realistic convective vigor and internal heating rate is characterized by either a spherical harmonic degree-1 planform with a major upwelling in one hemisphere and a major downwelling in the other hemisphere when continents are absent, or a degree-2 planform with two antipodal major upwellings when a supercontinent is present. We propose that due to modulation of continents, these two modes of mantle convection alternate within the Earth's mantle, causing the cyclic processes of assembly and breakup of supercontinents including Rodinia and Pangea in the last 1 Ga. Our model suggests that the largely degree-2 structure for the present-day mantle with the Africa and Pacific antipodal superplumes, is a natural consequence of this dynamic process of very long-wavelength mantle convection interacting with supercontinent Pangea. Our model explains the basic features of true polar wander (TPW) events for Rodinia and Pangea including their equatorial locations and large variability of TPW inferred from paleomagnetic studies. Our model also suggests that TPW is expected to be more variable and large during supercontinent assembly, but small after a supercontinent acquires its equatorial location and during its subsequent dispersal.

Mantle dynamics in super-Earths: Post-perovskite rheology and self-regulation of viscosity

NASA Astrophysics Data System (ADS)

Tackley, P. J.; Ammann, M.; Brodholt, J. P.; Dobson, D. P.; Valencia, D.

2013-07-01

The discovery of extra-solar "super-Earth" planets with sizes up to twice that of Earth has prompted interest in their possible lithosphere and mantle dynamics and evolution. Simple scalings suggest that super-Earths are more likely than an equivalent Earth-sized planet to be undergoing plate tectonics. Generally, viscosity and thermal conductivity increase with pressure while thermal expansivity decreases, resulting in lower convective vigour in the deep mantle, which, if extralopated to the largest super-Earths might, according to conventional thinking, result in no convection in their deep mantles due to the very low effective Rayleigh number. Here we evaluate this. First, as the mantle of a super-Earth is made mostly of post-perovskite we here extend the density functional theory (DFT) calculations of post-perovskite activation enthalpy of to a pressure of 1 TPa, for both slowest diffusion (upper-bound rheology) and fastest diffusion (lower-bound rheology) directions. Along a 1600 K adiabat the upper-bound rheology would lead to a post-perovskite layer of a very high (˜1030 Pa s) but relatively uniform viscosity, whereas the lower-bound rheology leads to a post-perovskite viscosity increase of ˜7 orders of magnitude with depth; in both cases the deep mantle viscosity would be too high for convection. Second, we use these DFT-calculated values in statistically steady-state numerical simulations of mantle convection and lithosphere dynamics of planets with up to ten Earth masses. The models assume a compressible mantle including depth-dependence of material properties and plastic yielding induced plate-like lithospheric behaviour. Results confirm the likelihood of plate tectonics for planets with Earth-like surface conditions (temperature and water) and show a self-regulation of deep mantle temperature. The deep mantle is not adiabatic; instead feedback between internal heating, temperature and viscosity regulates the temperature such that the viscosity has the

Linking the geological record for large igneous provinces and hotspots with tomography-based numerical models of thermal convection in the Earth's mantle

NASA Astrophysics Data System (ADS)

Glisovic, P.; Forte, A. M.; Rowley, D. B.; Simmons, N. A.; Grand, S. P.

2013-12-01

Current tomographic imaging of the 3-D structure in Earth's interior reveals several large-scale anomalies of strongly reduced seismic velocity in the deep lower mantle, in particular beneath the Perm region in Western Siberia, the East Pacific Rise, the West Pacific (Caroline Islands), the Southwest Indian Ocean, as well as under Western and Southern Africa. We have carried out mantle dynamic simulations (Glisovic et al., GJI 2012) of the evolution of these large-scale structures that directly incorporate robust constraints provided by joint seismic-geodynamic inversions of mantle density structure with further constraints provided by mineral physics data (Simmons et al., GJI 2009, JGR 2010). These tomography-based convection simulations also incorporate constraints on mantle viscosity inferred by inversion of a suite of convection-related and glacial isostatic adjustment data sets (Mitrovica & Forte, EPSL 2004) and are characterized by Earth-like Rayleigh numbers. The convection simulations provide a detailed insight into the very-long-time evolution of the buoyancy of these lower-mantle anomalies. We find, in particular, that the buoyancy associated with the 'Perm Anomaly' generates a very long-lived hot upwelling or 'superplume' that is connected to the paleomagnetic location of the Siberian Traps (Smirnov & Tarduno, EPSL 2010) and also to location of North Atlantic Igneous Provinces (i.e., the opening of North Atlantic Ocean). These convection simulations (both backwards and forwards in time) also reveal stable and long-lived plume-like upwellings under the East Pacific Rise, as previously identified by Rowley et al. (AGU 2011, Nature - in review), in particular beneath the Easter & Pitcairn hotspots. Finally we also provide detailed reconstructions of the 65 Myr evolution of the 'Reunion plume' that gave rise to the Deccan Traps.

Chondritic xenon in the Earth's mantle.

PubMed

Caracausi, Antonio; Avice, Guillaume; Burnard, Peter G; Füri, Evelyn; Marty, Bernard

2016-05-05

Noble gas isotopes are powerful tracers of the origins of planetary volatiles, and the accretion and evolution of the Earth. The compositions of magmatic gases provide insights into the evolution of the Earth's mantle and atmosphere. Despite recent analytical progress in the study of planetary materials and mantle-derived gases, the possible dual origin of the planetary gases in the mantle and the atmosphere remains unconstrained. Evidence relating to the relationship between the volatiles within our planet and the potential cosmochemical end-members is scarce. Here we show, using high-precision analysis of magmatic gas from the Eifel volcanic area (in Germany), that the light xenon isotopes identify a chondritic primordial component that differs from the precursor of atmospheric xenon. This is consistent with an asteroidal origin for the volatiles in the Earth's mantle, and indicates that the volatiles in the atmosphere and mantle originated from distinct cosmochemical sources. Furthermore, our data are consistent with the origin of Eifel magmatism being a deep mantle plume. The corresponding mantle source has been isolated from the convective mantle since about 4.45 billion years ago, in agreement with models that predict the early isolation of mantle domains. Xenon isotope systematics support a clear distinction between mid-ocean-ridge and continental or oceanic plume sources, with chemical heterogeneities dating back to the Earth's accretion. The deep reservoir now sampled by the Eifel gas had a lower volatile/refractory (iodine/plutonium) composition than the shallower mantle sampled by mid-ocean-ridge volcanism, highlighting the increasing contribution of volatile-rich material during the first tens of millions of years of terrestrial accretion.

The influence of water on mantle convection and plate tectonics

NASA Astrophysics Data System (ADS)

Brändli, S.; Tackley, P. J.

2017-12-01

Water has a significant influence to mantle rheology and therefore also to the convection of the mantle and the plate tectonics. The viscosity of the mantle can be decreased by up to two orders of magnitude when water is present in the mantle. Another effect of the water is the change in the solidus of the mantle and therefore the melting regime. This two effects of water in the mantle have a significant influence to mantle convection and plate tectonics. The influx of water to the mantle is driven by plate tectonics as wet oceanic lithosphere is subducted into the mantle and then brought back to the lithosphere and the surface by MOR-, arc- and hotspot volcanism. Studies show that the amount of water in the mantle is about three times bigger than the amount of water in the oceans. To model this water cycle multiple additions to StagYY are necessary. With the enhanced code we calculated multiple steady state models with a wide range of parameters to study the effect of water on the mantle rheology and the behavior of the lithosphere. The results will help us to understand the earths interior and its reaction and behavior under partially hydrated conditions.

On the formation of continental silicic melts in thermochemical mantle convection models: implications for early Earth

NASA Astrophysics Data System (ADS)

van Thienen, P.; van den Berg, A. P.; Vlaar, N. J.

2004-12-01

Important constituents of Archean cratons, formed in the early and hot history of the Earth, are Tonalite-Trondhjemite-Granodiorite (TTG) plutons and greenstone belts. The formation of these granite-greenstone terrains is often ascribed to plate-tectonic processes. Buoyancy considerations, however, do not allow plate tectonics to take place in a significantly hotter Earth. We therefore propose an alternative mechanism for the coeval and proximate production of TTG plutons and greenstone-like crustal successions. That is, when a locally anomalously thick basaltic crust has been produced by continued addition of extrusive or intrusive basalts due to partial melting of the underlying convecting mantle, the transition of a sufficient amount of basalt in the lower crust to eclogite may trigger a resurfacing event, in which a complete crustal section of over 1000 km long sinks into the mantle in less than 2 million years. Pressure release partial melting in the complementary upwelling mantle produces large volumes of basaltic material replacing the original crust. Partial melting at the base of this newly produced crust may generate felsic melts which are added as intrusives and/or extrusives to the generally mafic crustal succession, adding to what resembles a greenstone belt. Partial melting of metabasalt in the sinking crustal section produces a significant volume of TTG melt which is added to the crust directly above the location of 'subduction', presumably in the form of a pluton. This scenario is self-consistently produced by numerical thermochemical mantle convection models, presented in this paper, including partial melting of mantle peridotite and crustal (meta)basalt. The metamorphic p, T conditions under which partial melting of metabasalt takes place in this scenario are consistent with geochemical trace element data for TTGs, which indicate melting under amphibolite rather than eclogite facies. Other geodynamical settings which we have also investigated

Thermal and chemical convection in planetary mantles

NASA Technical Reports Server (NTRS)

Dupeyrat, L.; Sotin, C.; Parmentier, E. M.

1995-01-01

Melting of the upper mantle and extraction of melt result in the formation of a less dense depleted mantle. This paper describes series of two-dimensional models that investigate the effects of chemical buoyancy induced by these density variations. A tracer particles method has been set up to follow as closely as possible the chemical state of the mantle and to model the chemical buoyant force at each grid point. Each series of models provides the evolution with time of magma production, crustal thickness, surface heat flux, and thermal and chemical state of the mantle. First, models that do not take into account the displacement of plates at the surface of Earth demonstrate that chemical buoyancy has an important effect on the geometry of convection. Then models include horizontal motion of plates 5000 km wide. Recycling of crust is taken into account. For a sufficiently high plate velocity which depends on the thermal Rayleigh number, the cell's size is strongly coupled with the plate's size. Plate motion forces chemically buoyant material to sink into the mantle. Then the positive chemical buoyancy yields upwelling as depleted mantle reaches the interface between the upper and the lower mantle. This process is very efficient in mixing the depleted and undepleted mantle at the scale of the grid spacing since these zones of upwelling disrupt the large convective flow. At low spreading rates, zones of upwelling develop quickly, melting occurs, and the model predicts intraplate volcanism by melting of subducted crust. At fast spreading rates, depleted mantle also favors the formation of these zones of upwelling, but they are not strong enough to yield partial melting. Their rapid displacement toward the ridge contributes to faster large-scale homogenization.

Time-dependent convection models of mantle thermal structure constrained by seismic tomography and geodynamics: implications for mantle plume dynamics and CMB heat flux

NASA Astrophysics Data System (ADS)

Glišović, P.; Forte, A. M.; Moucha, R.

2012-08-01

One of the outstanding problems in modern geodynamics is the development of thermal convection models that are consistent with the present-day flow dynamics in the Earth's mantle, in accord with seismic tomographic images of 3-D Earth structure, and that are also capable of providing a time-dependent evolution of the mantle thermal structure that is as 'realistic' (Earth-like) as possible. A successful realization of this objective would provide a realistic model of 3-D mantle convection that has optimal consistency with a wide suite of seismic, geodynamic and mineral physical constraints on mantle structure and thermodynamic properties. To address this challenge, we have constructed a time-dependent, compressible convection model in 3-D spherical geometry that is consistent with tomography-based instantaneous flow dynamics, using an updated and revised pseudo-spectral numerical method. The novel feature of our numerical solutions is that the equations of conservation of mass and momentum are solved only once in terms of spectral Green's functions. We initially focus on the theory and numerical methods employed to solve the equation of thermal energy conservation using the Green's function solutions for the equation of motion, with special attention placed on the numerical accuracy and stability of the convection solutions. A particular concern is the verification of the global energy balance in the dissipative, compressible-mantle formulation we adopt. Such validation is essential because we then present geodynamically constrained convection solutions over billion-year timescales, starting from present-day seismically constrained thermal images of the mantle. The use of geodynamically constrained spectral Green's functions facilitates the modelling of the dynamic impact on the mantle evolution of: (1) depth-dependent thermal conductivity profiles, (2) extreme variations of viscosity over depth and (3) different surface boundary conditions, in this case mobile

The role of upper mantle mineral phase transitions on the current structure of large-scale Earth's mantle convection.

NASA Astrophysics Data System (ADS)

Thoraval, C.

2017-12-01

Describing the large-scale structures of mantle convection and quantifying the mass transfer between upper and lower mantle request to account for the role played by mineral phase transitions in the transition zone. We build a density distribution within the Earth mantle from velocity anomalies described by global seismic tomographic models. The density distribution includes thermal anomalies and topographies of the phase transitions at depths of 410 and 660 km. We compute the flow driven by this density distribution using a 3D spherical circulation model, which account for depth-dependent viscosity. The dynamic topographies at the surface and at the CMB and the geoid are calculated as well. Within the range of viscosity profiles allowing for a satisfying restitution of the long wavelength geoid, we perform a parametric study to decipher the role of the characteristics of phase diagrams - mainly the Clapeyron's slopes - and of the kinetics of phase transitions, which may modify phase transition topographies. Indeed, when a phase transition is delayed, the boundary between two mineral phases is both dragged by the flow and interfere with it. The results are compared to recent estimations of surface dynamic topography and to the phase transition topographies as revealed by seismic studies. The consequences are then discussed in terms of structure of mantle flow. Comparisons between various tomographic models allow us to enlighten the most robust features. At last, the role played by the phase transitions on the lateral variations of mass transfer between upper and lower mantle are quantified by comparison to cases with no phase transitions and confronted to regional tomographic models, which reflect the variability of the behaviors of the descending slabs in the transition zone.

Using the heterogeneity distribution in Earth's mantle to study structure and flow

NASA Astrophysics Data System (ADS)

Rost, S.; Frost, D. A.; Bentham, H. L.

2016-12-01

The Earth's interior contains heterogeneities on many scale-lengths ranging from continent sized structures such as Large-Low Shear Velocity Provinces (LLSVPs) to grain-sized anomalies resolved using geochemistry. Sources of heterogeneity in Earth's mantle are for example the recycling of crustal material through the subduction process as well as partial melting and compositional variations. The subduction and recycling of oceanic crust throughout Earth's history leads to strong heterogeneities in the mantle that can be detected using seismology and geochemistry. Current models of mantle convection show that the subducted crustal material can be long-lived and is transported passively throughout the mantle by convective flows. Settling and entrainment is dependent on the density structure of the heterogeneity. Imaging heterogeneities throughout the mantle therefore allows imaging mantle flow especially in areas of inhibited flow due to e.g. viscosity changes or changes in composition or dynamics. The short-period seismic wavefield is dominated by scattered seismic energy partly originating from scattering at small-scale heterogeneities in Earth's mantle. Using specific raypath configurations we are able to sample different depths throughout Earth's mantle for the existence and properties of heterogeneities. These scattering probes show distinct variations in energy content with frequency indicating dominant heterogeneity length-scales in the mantle. We detect changes in heterogeneity structure both in lateral and radial directions. The radial heterogeneity structure requires changes in mantle structure at depths of 1000 km and 1800 to 2000 km that could indicate a change in viscosity structure in the mid mantle partly changing the flow of subducted crustal material into the deep mantle. Lateral changes in heterogeneity structure close to the core mantle boundary indicate lateral transport inhibited by the compositional anomalies of the LLSVPs.

Episodic tectonic plate reorganizations driven by mantle convection

NASA Astrophysics Data System (ADS)

King, Scott D.; Lowman, Julian P.; Gable, Carl W.

2002-10-01

Periods of relatively uniform plate motion were interrupted several times throughout the Cenozoic and Mesozoic by rapid plate reorganization events [R. Hey, Geol. Soc. Am. Bull. 88 (1977) 1404-1420; P.A. Rona, E.S. Richardson, Earth Planet. Sci. Lett. 40 (1978) 1-11; D.C. Engebretson, A. Cox, R.G. Gordon, Geol. Soc. Am. Spec. Pap. 206 (1985); R.G. Gordon, D.M. Jurdy, J. Geophys. Res. 91 (1986) 12389-12406; D.A. Clague, G.B. Dalrymple, US Geol. Surv. Prof. Pap. 1350 (1987) 5-54; J.M. Stock, P. Molnar, Nature 325 (1987) 495-499; C. Lithgow-Bertelloni, M.A. Richards, Geophys. Res. Lett. 22 (1995) 1317-1320; M.A. Richards, C. Lithgow-Bertelloni, Earth Planet. Sci. Lett. 137 (1996) 19-27; C. Lithgow-Bertelloni, M.A. Richards, Rev. Geophys. 36 (1998) 27-78]. It has been proposed that changes in plate boundary forces are responsible for these events [M.A. Richards, C. Lithgow-Bertelloni, Earth Planet. Sci. Lett. 137 (1996) 19-27; C. Lithgow-Bertelloni, M.A. Richards, Rev. Geophys. 36 (1998) 27-78]. We present an alternative hypothesis: convection-driven plate motions are intrinsically unstable due to a buoyant instability that develops as a result of the influence of plates on an internally heated mantle. This instability, which has not been described before, is responsible for episodic reorganizations of plate motion. Numerical mantle convection experiments demonstrate that high-Rayleigh number convection with internal heating and surface plates is sufficient to induce plate reorganization events, changes in plate boundary forces, or plate geometry, are not required.

Continuous reorientation of synchronous terrestrial planets due to mantle convection

NASA Astrophysics Data System (ADS)

Leconte, Jérémy

2018-03-01

Many known rocky exoplanets are thought to have been spun down by tidal interactions to a state of synchronous rotation, in which a planet's period of rotation is equal to that of its orbit around its host star. Investigations into atmospheric and surface processes occurring on such exoplanets thus commonly assume that day and night sides are fixed with respect to the surface over geological timescales. Here we use an analytical model to show that true polar wander—where a planetary body's spin axis shifts relative to its surface because of changes in mass distribution—can continuously reorient a synchronous rocky exoplanet. As occurs on Earth, we find that even weak mantle convection in a rocky exoplanet can produce density heterogeneities within the mantle sufficient to reorient the planet. Moreover, we show that this reorientation is made very efficient by the slower rotation rate of a synchronous planet when compared with Earth, which limits the stabilizing effect of rotational and tidal deformations. Furthermore, a relatively weak lithosphere limits its ability to support remnant loads and stabilize against reorientation. Although uncertainties exist regarding the mantle and lithospheric evolution of these worlds, we suggest that the axes of smallest and largest moment of inertia of synchronous exoplanets with active mantle convection change continuously over time, but remain closely aligned with the star-planet and orbital axes, respectively.

The longevity of Archean mantle residues in the convecting upper mantle and their role in young continent formation

NASA Astrophysics Data System (ADS)

Liu, Jingao; Scott, James M.; Martin, Candace E.; Pearson, D. Graham

2015-08-01

The role played by ancient melt-depleted lithospheric mantle in preserving continental crust through time is critical in understanding how continents are built, disrupted and recycled. While it has become clear that much of the extant Archean crust is underpinned by Archean mantle roots, reports of Proterozoic melt depletion ages for peridotites erupted through Phanerozoic terranes raise the possibility that ancient buoyant lithospheric mantle acts as a "life-raft" for much of the Earth's continental crust. Here we report the largest crust-lithospheric mantle age decoupling (∼2.4 Ga) so far observed on Earth and examine the potential cause for such extreme age decoupling. The Phanerozoic (<300 Ma) continental crust of West Otago, New Zealand, is intruded by Cenozoic diatremes that have erupted cratonic mantle-like highly depleted harzburgites and dunites. These peridotites have rhenium depletion Os model ages that vary from 0.5 to 2.7 Ga, firmly establishing the record of an Archean depletion event. However, the vast range in depletion ages does not correlate with melt depletion or metasomatic tracer indices, providing little support for the presence of a significant volume of ancient mantle root beneath this region. Instead, the chemical and isotopic data are best explained by mixing of relict components of Archean depleted peridotitic mantle residues that have cycled through the asthenosphere over Ga timescales along with more fertile convecting mantle. Extensive melt depletion associated with the "docking" of these melt residues beneath the young continental crust of the Zealandia continent explains the decoupled age relationship that we observe today. Hence, the newly formed lithospheric root incorporates a mixture of ancient and modern mantle derived from the convecting mantle, cooled and accreted in recent times. We argue that in this case, the ancient components played no earlier role in continent stabilization, but their highly depleted nature along with

Multipathing within LLSVPs in models of thermal and thermochemical mantle convection

NASA Astrophysics Data System (ADS)

Nowacki, A.; Walpole, J.; Davies, R.; Heck, H. V.; Wookey, J. M.; Davies, H.

2016-12-01

Two regions at the base of Earth's mantle (the Large Low Shear Velocity Provinces, or LLSVPs) are seismically slow and comprise a large proportion of the lowermost few hundred km of the mantle. It is debated whether these regions might be the remnants of a basal magma ocean or other Earth-forming processes, in which case the regions may provide information about Earth's history. However, it is still uncertain what the current physical properties of the LLSVPs are. Is the cause of the LLSVPs' seismic signature primarily thermal or chemical? One argument for a largely chemical origin is that seismically `sharp sides' can be inferred from waves which exhibit `multipathing' (the arrival of more than one wave due to refraction) when traversing these regions. This implies strong gradients in velocity which are seemingly unlikely to be sustained in a purely thermal situation, where diffusive processes and convection would act to equilibrate temperatures over short length scales (of a few tens of km). We address this by simulating mantle convection with Earth-like parameters in a 3D spherical geometry for two end-member cases: an isochemical (T) mantle, and a thermochemical (TC) case where a global, dense layer exists initially at the base of the mantle. We impose 200 Ma of plate motion history, track the location of the dense material, and convert both models to seismic velocity using a thermodynamical database (Stixrude & Lithgow-Bertelloni, GJI, 2005, 2011). Previous work has shown the cases are not easily distinguishable tomographically, so we seek to reproduce observations of `sharp sides' by creating finite-frequency synthetics at relatively high frequencies ( 0.2 Hz), using the spectral element method. We find that in a number of regions in both T and TC models, we observe multipathing in Sdiff waves which traverse the LLSVPs. Events beneath Tonga recorded in southern Africa yield strongly azimuth-dependent arrival times, as expected by features in the convection

Linking TERRA and DRex to relate mantle convection and seismic anisotropy

NASA Astrophysics Data System (ADS)

Walker, Andrew; Davies, Huw; Davies, Rhodri; Wookey, James

2015-04-01

Seismic anisotropy caused by flow induced alignment of the olivine crystals in Earth's upper mantle provides a powerful way to test our ideas of mantle convection. We have been working to directly combine computer simulations of mantle dynamics, using fluid mechanics at the continuum scale, with models of rock deformation to capture fabric evolution at the grain scale. By combining models of deformation at these two scales we hope to be able to rigorously test hypothesis linking mantle flow to seismic anisotropy in regions as diverse as subduction zones, the lithosphere-asthenosphere boundary, and the transition zone. We also intend to permit feedback, for example via geometrical softening, from the model of fabric development into the material properties used in the convection simulation. We are building a flexible framework for this approach which we call Theia. Our initial implementation uses the TERRA convection code (Baumgardner, J. Stat. Phys. 39:501-511, 1985; Davies et al. Geosci. Model Dev. 6:1095-1107, 2013) to drive DRex (Kaminski et al. Geophys. J. Int. 158:744-752, 2004), which is used to predict the evolution of crystallographic preferred orientation in the upper mantle. Here we describe our current implementation which makes use of the ability of TERRA to track markers, or particles, through the evolving flow field. These tracers have previously been used to track attributes such as the bulk chemical composition or trace element ratios. Our modification is to use this technology to track a description of the current state of the texture and microstructure (encompassing an orientation distribution function, grain size parameters and dislocation density) such that we can advance models of polycrystalline deformation for many simultaneous DRex instances alongside and in sync with models of mantle convection. We will also describe initial results from our first use of the Theia framework where we are investigating the effect of asthenospheric viscosity

Mantle convection: concensus and queries (Augustus Love Medal Lecture)

NASA Astrophysics Data System (ADS)

Ricard, Y.

2012-04-01

Thermal convection driven by surface cooling and internal heat production is the cause of endogenic activity of all planets, expressed as tectonic activity and volcanism for solid planets. The sluggish convection of the silicated mantle also controls the activity of the metallic core and the possibility of an active dynamo. A glimpse of the internal structure of Earth's mantle is provided by seismic tomography. However, both the limited resolution of seismic methods and the complexity of the relations between seismic velocities and the thermo-mechanical parameters (mostly temperature and density), leave to the geodynamicist a large degree of interpretation. At first order, a very simple model of mantle heterogeneities, only built from the paleogeographic positions of Cenozoic and Mesozoic slabs, explains the pattern and amplitude of Earth's plate motions and gravity field, while being in agreement with long wavelength tomography. This indicates that the mantle dynamics is mostly controlled by thermal anomalies and by the dynamics of the top boundary layer, the lithosphere. However, the presence of various complexities due to variations in elemental composition and to phase transitions is required by seismology, mineralogy and geochemistry. I will review how these complexities affect the dynamics of the transition zone and of the deep mantle and discuss the hypothesis on their origins, either primordial or as a consequence of plate tectonics. The rheologies that are used in global geodynamic models for the mantle and the lithosphere remain very simplistic. Some aspects of plate tectonics (e.g., the very existence of plates, their evolution, the dynamics of one-sided subductions...) are now reproduced by numerical simulations. However the rheologies implemented and their complexities remain only remotely related to that of solid minerals as observed in laboratories. The connections between the quantities measured at microscopic scale (e.g., mineralogy, grainsize

Thermal-chemical Mantle Convection Models With Adaptive Mesh Refinement

NASA Astrophysics Data System (ADS)

Leng, W.; Zhong, S.

2008-12-01

In numerical modeling of mantle convection, resolution is often crucial for resolving small-scale features. New techniques, adaptive mesh refinement (AMR), allow local mesh refinement wherever high resolution is needed, while leaving other regions with relatively low resolution. Both computational efficiency for large- scale simulation and accuracy for small-scale features can thus be achieved with AMR. Based on the octree data structure [Tu et al. 2005], we implement the AMR techniques into the 2-D mantle convection models. For pure thermal convection models, benchmark tests show that our code can achieve high accuracy with relatively small number of elements both for isoviscous cases (i.e. 7492 AMR elements v.s. 65536 uniform elements) and for temperature-dependent viscosity cases (i.e. 14620 AMR elements v.s. 65536 uniform elements). We further implement tracer-method into the models for simulating thermal-chemical convection. By appropriately adding and removing tracers according to the refinement of the meshes, our code successfully reproduces the benchmark results in van Keken et al. [1997] with much fewer elements and tracers compared with uniform-mesh models (i.e. 7552 AMR elements v.s. 16384 uniform elements, and ~83000 tracers v.s. ~410000 tracers). The boundaries of the chemical piles in our AMR code can be easily refined to the scales of a few kilometers for the Earth's mantle and the tracers are concentrated near the chemical boundaries to precisely trace the evolvement of the boundaries. It is thus very suitable for our AMR code to study the thermal-chemical convection problems which need high resolution to resolve the evolvement of chemical boundaries, such as the entrainment problems [Sleep, 1988].

Predicting lower mantle heterogeneity from 4-D Earth models

NASA Astrophysics Data System (ADS)

Flament, Nicolas; Williams, Simon; Müller, Dietmar; Gurnis, Michael; Bower, Dan J.

2016-04-01

The Earth's lower mantle is characterized by two large-low-shear velocity provinces (LLSVPs), approximately ˜15000 km in diameter and 500-1000 km high, located under Africa and the Pacific Ocean. The spatial stability and chemical nature of these LLSVPs are debated. Here, we compare the lower mantle structure predicted by forward global mantle flow models constrained by tectonic reconstructions (Bower et al., 2015) to an analysis of five global tomography models. In the dynamic models, spanning 230 million years, slabs subducting deep into the mantle deform an initially uniform basal layer containing 2% of the volume of the mantle. Basal density, convective vigour (Rayleigh number Ra), mantle viscosity, absolute plate motions, and relative plate motions are varied in a series of model cases. We use cluster analysis to classify a set of equally-spaced points (average separation ˜0.45°) on the Earth's surface into two groups of points with similar variations in present-day temperature between 1000-2800 km depth, for each model case. Below ˜2400 km depth, this procedure reveals a high-temperature cluster in which mantle temperature is significantly larger than ambient and a low-temperature cluster in which mantle temperature is lower than ambient. The spatial extent of the high-temperature cluster is in first-order agreement with the outlines of the African and Pacific LLSVPs revealed by a similar cluster analysis of five tomography models (Lekic et al., 2012). Model success is quantified by computing the accuracy and sensitivity of the predicted temperature clusters in predicting the low-velocity cluster obtained from tomography (Lekic et al., 2012). In these cases, the accuracy varies between 0.61-0.80, where a value of 0.5 represents the random case, and the sensitivity ranges between 0.18-0.83. The largest accuracies and sensitivities are obtained for models with Ra ≈ 5 x 107, no asthenosphere (or an asthenosphere restricted to the oceanic domain), and a

Development of a Mantle Convection Physical Model to Assist with Teaching about Earth's Interior Processes

NASA Astrophysics Data System (ADS)

Glesener, G. B.; Aurnou, J. M.

2010-12-01

The Modeling and Educational Demonstrations Laboratory (MEDL) at UCLA is developing a mantle convection physical model to assist educators with the pedagogy of Earth’s interior processes. Our design goal consists of two components to help the learner gain conceptual understanding by means of visual interactions without the burden of distracters, which may promote alternative conceptions. Distracters may be any feature of the conceptual model that causes the learner to use inadequate mental artifact to help him or her understand what the conceptual model is intended to convey. The first component, and most important, is a psychological component that links properties of “everyday things” (Norman, 1988) to the natural phenomenon, mantle convection. Some examples of everyday things may be heat rising out from a freshly popped bag of popcorn, or cold humid air falling from an open freezer. The second component is the scientific accuracy of the conceptual model. We would like to simplify the concepts for the learner without sacrificing key information that is linked to other natural phenomena the learner will come across in future science lessons. By taking into account the learner’s mental artifacts in combination with a simplified, but accurate, representation of what scientists know of the Earth’s interior, we expect the learner to have the ability to create an adequate qualitative mental simulation of mantle convection. We will be presenting some of our prototypes of this mantle convection physical model at this year’s poster session and invite constructive input from our colleagues.

Dynamical Generation of the Transition Zone in the Earth's Mantle

NASA Astrophysics Data System (ADS)

Hansen, U.; Stemmer, K.

2005-12-01

The internal structure of the Earth is made up by a series of layers, though it is unclear how many layers exist and if there are layers invisible to remote sensing techniques. The transition zone is likely to be a boundary layer separating the convective systems in the lower and upper mantle. It seems likely that currently there is some mass exchange across this boundary, rather than the two systems beeing strictly separated.a Double-diffusive convection(d.d.c) is a vital mechanism which can generate layered structure and may thus be an important mmical machinery behind the formation of the transition zone. Double-diffusive convection determines the dynamics of systems whose density is influenced by at least two components with different molecular diffusivities.In the mantle, composition and temperature play the role of those two components. By means of numerical experiments we demonstrate that under mantle relevant conditions d.d.c typically leads to the formation of a transition zone. The calculations encompass two- and three dimensional Cartesian geometries as well as fully 3D spherical domains. We have further included strongly temperature dependent viscosity and find that this leads to even more pronounced layering. In most cases a layered flow pattern emerges, where two layers with a transition zone in between resembles a quasistationary state. Thus, the transition zone can be the result of a self organization process of the convective flow in the mantle. The presence of a phase transition further helps to stabilize the boundary against overturning, even on a time scale on the order of the age of the Earth.

The role of thermodynamics in mantle convection: is mantle-layering intermittent?

NASA Astrophysics Data System (ADS)

Stixrude, L. P.; Cagney, N.; Lithgow-Bertelloni, C. R.

2016-12-01

We examine the thermal evolution of the Earth using a 1D model in which mixing length theory is used to characterise the role of thermal convection. Unlike previous work, our model accounts for the complex role of thermodynamics and phase changes through the use of HeFESTo (Stixrude & Lithgow-Bertelloni, Geophys. J. Int. 184, 2011), a comprehensive thermodynamic model that enables self-consistent computation of phase equilibria, physical properties (e.g. density, thermal expansivity etc.) and mantle isentropes. Our model also accounts for the freezing of the inner core, radiogenic heating and Arrhenius rheology, and is validated by comparing our results to observations, including the present-day size of the inner core and the heat flux at the surface.If phase changes and the various thermodynamic effects on mantle properties are neglected, the results are weakly dependent on the initial conditions, as has been observed in several previous studies. However, when these effects are accounted for, the initial temperature profile has a strong influence on the thermal evolution of the mantle, because small changes in the temperature and phase-assemblage can lead to large changes in the local physical properties and the adiabatic gradient.The inclusion of thermodynamic effects leads to some new and interesting insights. We demonstrate that the Clapeyron slope and the thermal gradient at the transition zone both vary significantly with time; this causes the mantle to switch between a layered state, in which convection across the transition zone is weak or negligible, and an un-layered state, in which there is no resistance to mass transfer between the upper and lower mantles.Various plume models describe plumes either rising directly from the CMB to the lithosphere, or stalling at the transition zone before spawning new plumes in the upper mantle. The observance of switching behaviour indicates that both models may be applicable depending on the state of the mantle: plumes

Supercontinent Formation in 3-D Spherical Mantle Convection Models With Multiple Continental Blocks

NASA Astrophysics Data System (ADS)

Zhang, N.; Zhong, S.; McNamara, A.

2007-12-01

Much of the large-scale tectonics on the Earth in the last Ga is predominated by the assembly and breakup of supercontinents Rodinia and Pangea. However, the mechanism that is responsible for supercontinent formation remains poorly understood. Zhong et al [2007] recently showed that mantle convection with moderately strong lithosphere and lower mantle is characterized by a largely degree-1 planform in which one hemisphere is predominated by upwellings while the other by downwellings. They further suggested that the downwellings should attract all the continental blocks to merge in the downwelling hemisphere, thus leading to supercontinent formation there. However, Zhong et al. [2007] did not consider drifting and collision processes of continents. In this study, we explore the supercontinent formation mechanisms by including drifting and collision processes of multiple continental blocks in 3-D spherical mantle convection models. We use thermochemical CitcomS code to model 3-D spherical mantle convection with continental blocks. In our models, particles are used to represent continents and to track their motions. We found that for models with mantle viscosity (i.e., moderately strong lithosphere and lower mantle) that leads to degree-1 convection as reported in Zhong et al. [2007], initially evenly- distributed continental blocks always merge to form a supercontinent on a time-scale of about 6 transit times (i.e., corresponding to about 300 Ma). The hemisphere where a supercontinent is formed is predominated by downwellings as continents merge towards there, while the other hemisphere by upwellings. However, after the supercontinent formation, upwellings are generated beneath the supercontinent. This scenario is qualitatively consistent with what Zhong et al. [2007] proposed. We also found that while some convection models with intrinsically small-scale planforms may also lead to formation of a supercontinent, some other models may fail to produce a supercontinent

Long-Term Stability of Plate-Like Behavior Caused by Hydrous Mantle Convection and Water Absorption in the Deep Mantle

NASA Astrophysics Data System (ADS)

Nakagawa, Takashi; Iwamori, Hikaru

2017-10-01

We investigate the cycling of water (regassing, dehydration, and degassing) in mantle convection simulations as a function of the strength of the oceanic lithosphere and its influence on the evolution of mantle water content. We also consider pseudo-plastic yielding with a friction coefficient for simulating brittle behavior of the plates and the water-weakening effect of mantle materials. This model can generate long-term plate-like behavior as a consequence of the water-weakening effect of mantle minerals. This finding indicates that water cycling plays an essential role in generating tectonic plates. In vigorous plate motion, the mantle water content rapidly increases by up to approximately 4-5 ocean masses, which we define as the "burst" effect. A burst is related to the mantle temperature and water solubility in the mantle transition zone. When the mantle is efficiently cooled down, the mantle transition zone can store water transported by the subducted slabs that can pass through the "choke point" of water solubility. The onset of the burst effect is strongly dependent on the friction coefficient. The burst effect of the mantle water content could have significantly influenced the evolution of the surface water if the burst started early, in which case the Earth's surface cannot preserve the surface water over the age of the Earth.

The source location of mantle plumes from 3D spherical models of mantle convection

NASA Astrophysics Data System (ADS)

Li, Mingming; Zhong, Shijie

2017-11-01

Mantle plumes are thought to originate from thermal boundary layers such as Earth's core-mantle boundary (CMB), and may cause intraplate volcanism such as large igneous provinces (LIPs) on the Earth's surface. Previous studies showed that the original eruption sites of deep-sourced LIPs for the last 200 Myrs occur mostly above the margins of the seismically-observed large low shear velocity provinces (LLSVPs) in the lowermost mantle. However, the mechanism that leads to the distribution of the LIPs is not clear. The location of the LIPs is largely determined by the source location of mantle plumes, but the question is under what conditions mantle plumes form outside, at the edges, or above the middle of LLSVPs. Here, we perform 3D geodynamic calculations and theoretical analyses to study the plume source location in the lowermost mantle. We find that a factor of five decrease of thermal expansivity and a factor of two increase of thermal diffusivity from the surface to the CMB, which are consistent with mineral physics studies, significantly reduce the number of mantle plumes forming far outside of thermochemical piles (i.e., LLSVPs). An increase of mantle viscosity in the lowermost mantle also reduces number of plumes far outside of piles. In addition, we find that strong plumes preferentially form at/near the edges of piles and are generally hotter than that forming on top of piles, which may explain the observations that most LIPs occur above LLSVP margins. However, some plumes originated at pile edges can later appear above the middle of piles due to lateral movement of the plumes and piles and morphologic changes of the piles. ∼65-70% strong plumes are found within 10 degrees from pile edges in our models. Although plate motion exerts significant controls over the large-scale mantle convection in the lower mantle, mantle plume formation at the CMB remains largely controlled by thermal boundary layer instability which makes it difficult to predict geographic

Convection experiments in a centrifuge and the generation of plumes in a very viscous fluid. [for earth mantle models

NASA Technical Reports Server (NTRS)

Nataf, H.-C.; Hager, B. H.; Scott, R. F.

1984-01-01

In this paper, experiments are described for which inertial effects are negligible. A small aspect-ratio tank filled with a very viscous fluid (Pr = 10 to the 6th) is used to observe the behavior of convection for Rayleigh numbers up to 6.3 x 10 to the 5th. These high values are reached by conducting the experiment in a centrifuge which provides a 130-fold increase in apparent gravity. Rotational effects are small, but cannot be totally dismissed. In this geometry, thermal boundary layer instabilities are indeed observed, and are found to be very similar to their lower Prandtl number counterparts. It is tentatively concluded that once given a certain degree of 'vulnerability' convection can develop 'plume' like instabilities, even when the Prandtl number is infinite. The concept is applied to the earth's mantle and it is speculated that 'plumes' could well be the dominant mode of small-scale convection under the lithospheric plates.

Numerical Mantle Convection Models With a Flexible Thermodynamic Interface

NASA Astrophysics Data System (ADS)

van den Berg, A. P.; Jacobs, M. H.; de Jong, B. H.

2001-12-01

Accurate material properties are needed for deep mantle (P,T) conditions in order to predict the longterm behavior of convection planetary mantles. Also the interpretation of seismological observations concerning the deep mantle in terms of mantle flow models calls for a consistent thermodynamical description of the basic physical parameters. We have interfaced a compressible convection code using the anelastic liquid approach based on finite element methods, to a database containing a full thermodynamic description of mantle silicates (Ita and King, J. Geophys. Res., 99, 15,939-15,940, 1994). The model is based on high resolution (P,T) tables of the relevant thermodynamic properties containing typically 50 million (P,T) table gridpoints to obtain resolution in (P,T) space of 1 K and an equivalent of 1 km. The resulting model is completely flexible such that numerical mantle convection experiments can be performed for any mantle composition for which the thermodynamic database is available. We present results of experiments for 2D cartesian models using a data base for magnesium-iron silicate in a pyrolitic composition (Stixrude and Bukowinski, Geoph.Monogr.Ser., 74, 131-142, 1993) and a recent thermodynamical model for magnesium silicate for the complete mantle (P,T) range, (Jacobs and Oonk, Phys. Chem. Mineral, 269, inpress 2001). Preliminary results of bulksound velocity distribution derived in a consistent way from the convection results and the thermodynamic database show a `realistic' mantle profile with bulkvelocity variations decreasing from several percent in the upper mantle to less than a percent in the deep lower mantle.

Constraints on mantle viscosity from convection models with plate motion history

NASA Astrophysics Data System (ADS)

Mao, W.; Zhong, S.

2017-12-01

The Earth's long-wavelength geoid and dynamic topography are mainly controlled by the mantle buoyancy and viscosity structure. Previous dynamical models for the geoid provide constraints on the 1-D mantle viscosity, using mantle buoyancy derived from seismic topography models. However, it is a challenge in these studies on how to convert seismic velocity to density anomalies and mantle buoyancy. Furthermore, these studies provide constraints only on relative viscosity variations but not on absolute magnitude of viscosity. In this study, we formulate time-dependent 3-D spherical mantle convection models with imposed plate motion history and seek constraints on mantle viscosity structure for both its radial relative variations and its absolute magnitude (i.e., Rayleigh number), using the geoid from the convection models. We found that the geoid at intermediate wavelengths of degrees 4-9 is mainly controlled by the subducted slabs in the upper mantle and the upper part of lower mantle that result from subduction from the last 50 Myr or the Cenozoic. To fit the degrees 4-9 geoid, we need viscosity contrast β defined as the ratio of the lower mantle viscosity and the asthenospheric viscosity to be larger than 2000 and Ra to be 1e8 (defined by the Earth's radius). The best fit model leads to 57% variance reduction and 76% correlation between the model and the observations. However, the long-wavelength geoid at degrees 2-3 is controlled by the lower mantle structure which requires much longer time scale to develop, as seen from our modeling. The preferred viscosity structure and Rayleigh number as constrained by the Cenozoic plate motion and the degrees 4-9 geoid no longer provide adequate fit to the geoid in models with the plate motion history for the last 450 Myr. The degrees 4-9 geoid amplitude is smaller for the models with longer plate motion history and a smaller Ra is required to fit the observation. In order to satisfy the relative amplitude between degrees 2

Gravitational Core-Mantle Coupling and the Acceleration of the Earth

NASA Technical Reports Server (NTRS)

Rubincam, David Parry; Smith, David E. (Technical Monitor)

2001-01-01

Gravitational core-mantle coupling may be the cause of the observed variable acceleration of the Earth's rotation on the 1000 year timescale. The idea is that density inhomogeneities which randomly come and go in the liquid outer core gravitationally attract density inhomogeneities in the mantle and crust, torquing the mantle and changing its rotation state. The corresponding torque by the mantle on the core may also explain the westward drift of the magnetic field of 0.2 deg per year. Gravitational core-mantle coupling would stochastically affect the rate of change of the Earth's obliquity by just a few per cent. Its contribution to polar wander would only be about 0.5% the presently observed rate. Tidal friction is slowing down the rotation of the Earth, overwhelming a smaller positive acceleration from postglacial rebound. Coupling between the liquid outer core of the Earth and the mantle has long been a suspected reason for changes in the length-of-day. The present investigation focuses on the gravitational coupling between the density anomalies in the convecting liquid outer core and those in the mantle and crust as a possible cause for the observed nonsecular acceleration on the millenial timescale. The basic idea is as follows. There are density inhomogeneities caused by blobs circulating in the outer core like the blobs in a lava lamp; thus the outer core's gravitational field is not featureless. Moreover, these blobs will form and dissipate somewhat randomly. Thus there will be a time variability to the fields. These density inhomogeneities will gravitationally attract the density anomalies in the mantle.

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.

Internally heated mantle convection and the thermal and degassing history of the earth

NASA Technical Reports Server (NTRS)

Williams, David R.; Pan, Vivian

1992-01-01

An internally heated model of parameterized whole mantle convection with viscosity dependent on temperature and volatile content is examined. The model is run for 4l6 Gyr, and temperature, heat flow, degassing and regassing rates, stress, and viscosity are calculated. A nominal case is established which shows good agreement with accepted mantle values. The effects of changing various parameters are also tested. All cases show rapid cooling early in the planet's history and strong self-regulation of viscosity due to the temperature and volatile-content dependence. The effects of weakly stress-dependent viscosity are examined within the bounds of this model and are found to be small. Mantle water is typically outgassed rapidly to reach an equilibrium concentration on a time scale of less than 200 Myr for almost all models, the main exception being for models which start out with temperatures well below the melting temperature.

Constraints on Thermochemical Convection of the Mantle from Plume-related Observations

NASA Astrophysics Data System (ADS)

Zhong, S.

2005-05-01

Although geochemical observations have long suggested a layered mantle with more enriched mantle material in the bottom layer to provide a significant amount of heat to the top layer, the nature of such a layering remains unclear. An important observation that has been used to argue against the conventional layered mantle model (i.e., the layering at the 670 km depth) was the plume heat flux [Davies, 1999]. Plume heat flux is estimated as ~ 3.5 TW, or 10% of the surface heat flux [Davies, 1988; Sleep, 1990]. In this study, we demonstrate with 3-D spherical models of mantle convection with depth- and temperature-dependent viscosity that observed plume heat flux, plume excess temperature (<350°C), and upper mantle temperature (~ 1300°C) can pose important constraints on the layered mantle convection. We show that for a purely thermal convection model (i.e., a whole mantle convection), the observations of plume heat flux, plume excess temperature, and upper mantle temperature can be simultaneously explained only when internal heating rate is about 65%. For smaller internal heating rate, plume heat flux and plume excess temperature would be too large, and upper mantle temperature would be too small, compared with the observed. This suggests that for a whole mantle convection the CMB heat flux needs to be > 10 TW. For a core with no significant heat producing elements, such large CMB heat flux may lead to too rapid cooling of the core or a too young inner core. A layered mantle convection may help reduce the CMB heat flux. For layered convection models, we found that the top layer needs to be ~70% internally heated to explain the upper mantle temperature and plume-related observations, and this required internal heating ratio is insensitive to the layer thickness for the bottom layer (we used ~600 km and 1100 km thicknesses). This result suggests that heat generation rate for the bottom layer cannot be significantly larger (< a factor of 2) than that for the top layer

A new conceptual model for whole mantle convection and the origin of hotspot plumes

NASA Astrophysics Data System (ADS)

Yoshida, Masaki

2014-08-01

A new conceptual model of mantle convection is constructed for consideration of the origin of hotspot plumes, using recent evidence from seismology, high-pressure experiments, geodynamic modeling, geoid inversion studies, and post-glacial rebound analyses. This conceptual model delivers several key points. Firstly, some of the small-scale mantle upwellings observed as hotspots on the Earth's surface originate at the base of the mantle transition zone (MTZ), in which the Archean granitic continental material crust (TTG; tonalite-trondhjemite-granodiorite) with abundant radiogenic elements is accumulated. Secondly, the TTG crust and the subducted oceanic crust that have accumulated at the base of MTZ could act as thermal or mechanical insulators, leading to the formation of a hot and less viscous layer just beneath the MTZ; which may enhance the instability of plume generation at the base of the MTZ. Thirdly, the origin of some hotspot plumes is isolated from the large low shear-wave velocity provinces (LLSVPs) under Africa and the South Pacific. I consider that the conceptual model explains why almost all the hotspots around Africa are located above the margins of the African LLSVP. Because a planetary-scale trench system surrounding a “Pangean cell” has been spatially stable throughout the Phanerozoic, a large amount of the oceanic crustal layer is likely to be trapped in the MTZ under the Pangean cell. Therefore, under Africa, almost all of the hotspot plumes originate from the base of the MTZ, where a large amount of TTG and/or oceanic crusts has accumulated. This conceptual model may explain the fact that almost all the hotspots around Africa are located on margins above the African LLSVP. It is also considered that some of the hotspot plumes under the South Pacific thread through the TTG/oceanic crusts accumulated around the bottom of the MTZ, and some have their roots in the South Pacific LLSVP while others originate from the MTZ. The numerical simulations

Early formation and long-term stability of continents resulting from decompression melting in a convecting mantle

NASA Astrophysics Data System (ADS)

De Smet, J.; Van den Berg, A. P.; Vlaar, N. J.

2000-07-01

The origin of stable old continental cratonic roots is still debated. We present numerical modelling results which show rapid initial formation during the Archaean of continental roots of ca. 200 km thick. These results have been obtained from an upper mantle thermal convection model including differentiation by pressure release partial melting of mantle peridotite. The upper mantle model includes time-dependent radiogenic heat production and thermal coupling with a heat reservoir representing the Earth's lower mantle and core. This allows for model experiments including secular cooling on a time-scale comparable to the age of the Earth. The model results show an initial phase of rapid continental root growth of ca. 0.1 billion year, followed by a more gradual increase of continental volume by addition of depleted material produced through hot diapiric, convective upwellings which penetrate the continental root from below. Within ca. 0.6 Ga after the start of the experiment, secular cooling of the mantle brings the average geotherm below the peridotite solidus thereby switching off further continental growth. At this time the thickness of the continental root has grown to ca. 200 km. After 1 Ga of secular cooling small scale thermal instabilities develop at the bottom of the continental root causing continental delamination without breaking up the large scale layering. This delaminated material remixes with the deeper layers. Two more periods, each with a duration of ca. 0.5 Ga and separated by quiescent periods were observed when melting and continental growth was reactivated. Melting ends at 3 Ga. Thereafter secular cooling proceeds and the compositionally buoyant continental root is stabilized further through the increase in mechanical strength induced by the increase of the temperature dependent mantle viscosity. Fluctuating convective velocity amplitudes decrease to below 10 mma -1 and the volume average temperature of the sub-continental convecting mantle has

Early Earth slab stagnation

NASA Astrophysics Data System (ADS)

Agrusta, R.; Van Hunen, J.

2016-12-01

At present day, the Earth's mantle exhibits a combination of stagnant and penetrating slabs within the transition zone, indicating a intermittent convection mode between layered and whole-mantle convection. Isoviscous thermal convection calculations show that in a hotter Earth, the natural mode of convection was dominated by double-layered convection, which may imply that slabs were more prone to stagnate in the transition zone. Today, slab penetration is to a large extent controlled by trench mobility for a plausible range of lower mantle viscosity and Clapeyron slope of the mantle phase transitions. Trench mobility is, in turn, governed by slab strength and density and upper plate forcing. In this study, we systematically investigate the slab-transition zone internation in the Early Earth, using 2D self-consistent numerical subduction models. Early Earth's higher mantle temperature facilitates decoupling between the plates and the underlying asthenosphere, and may result in slab sinking almost without trench retreat. Such behaviour together with a low resistance of a weak lower mantle may allow slabs to penetrate. The ability of slab to sink into the lower mantle throughout Earth's history may have important implications for Earth's evolution: it would provide efficient mass and heat flux through the transition zone therefore provide an efficient way to cool and mix the Earth's mantle.

The Elephants' Graveyard: Constraints from Mantle Plumes on the Fate of Subducted Slabs and Implications for the Style of Mantle Convection

NASA Astrophysics Data System (ADS)

Lassiter, J. C.

2007-12-01

The style of mantle convection (e.g., layered- vs. whole-mantle convection) is one of the most hotly contested questions in the Geological Sciences. Geochemical arguments for and against mantle layering have largely focused on mass-balance evidence for the existence of "hidden" geochemical reservoirs. However, the size and location of such reservoirs are largely unconstrained, and most geochemical arguments for mantle layering are consistent with a depleted mantle comprising most of the mantle mass and a comparatively small volume of enriched, hidden material either within D" or within seismically anomalous "piles" beneath southern Africa and the South Pacific. The mass flux associated with subduction of oceanic lithosphere is large and plate subduction is an efficient driver of convective mixing in the mantle. Therefore, the depth to which oceanic lithosphere descends into the mantle is effectively the depth of the upper mantle in any layered mantle model. Numerous geochemical studies provide convincing evidence that many mantle plumes contain material which at one point resided close to the Earth's surface (e.g., recycled oceanic crust ± sediments, possibly subduction-modified mantle wedge material). Fluid dynamic models further reveal that only the central cores of mantle plumes are involved in melt generation. The presence of recycled material in the sources of many ocean island basalts therefore cannot be explained by entrainment of this material during plume ascent, but requires that recycled material resides within or immediately above the thermo-chemical boundary layer(s) that generates mantle plumes. More recent Os- isotope studies of mantle xenoliths from OIB settings reveal the presence not only of recycled crust in mantle plumes, but also ancient melt-depleted harzburgite interpreted to represent ancient recycled oceanic lithosphere [1]. Thus, there is increasing evidence that subducted slabs accumulate in the boundary layer(s) that provide the source

The role of solid-solid phase transitions in mantle convection

NASA Astrophysics Data System (ADS)

Faccenda, Manuele; Dal Zilio, Luca

2017-01-01

With changing pressure and temperature conditions, downwelling and upwelling crustal and mantle rocks experience several solid-solid phase transitions that affect the mineral physical properties owing to structural changes in the crystal lattice and to the absorption or release of latent heat. Variations in density, together with phase boundary deflections related to the non-null reaction slope, generate important buoyancy forces that add to those induced by thermal perturbations. These buoyancy forces are proportional to the density contrast between reactant and product phases, their volume fraction, the slope and the sharpness of the reaction, and affect the style of mantle convection depending on the system composition. In a homogeneous pyrolitic mantle there is little tendency for layered convection, with slabs that may stagnate in the transition zone because of the positive buoyancy caused by post-spinel and post-ilmenite reactions, and hot plumes that are accelerated by phase transformations in the 600-800 km depth range. By adding chemical and mineralogical heterogeneities as on Earth, phase transitions introduce bulk rock and volatiles filtering effects that generate a compositional gradient throughout the entire mantle, with levels that are enriched or depleted in one or more of these components. Phase transitions often lead to mechanical softening or hardening that can be related to a different intrinsic mechanical behaviour and volatile solubility of the product phases, the heating or cooling associated with latent heat, and the transient grain size reduction in downwelling cold material. Strong variations in viscosity would enhance layered mantle convection, causing slab stagnation and plume ponding. At low temperatures and relatively dry conditions, reactions are delayed due to the sluggish kinetics, so that non-equilibrium phase aggregates can persist metastably beyond the equilibrium phase boundary. Survival of low-density metastable olivine

Revealing the Earth's mantle from the tallest mountains using the Jinping Neutrino Experiment.

PubMed

Šrámek, Ondřej; Roskovec, Bedřich; Wipperfurth, Scott A; Xi, Yufei; McDonough, William F

2016-09-09

The Earth's engine is driven by unknown proportions of primordial energy and heat produced in radioactive decay. Unfortunately, competing models of Earth's composition reveal an order of magnitude uncertainty in the amount of radiogenic power driving mantle dynamics. Recent measurements of the Earth's flux of geoneutrinos, electron antineutrinos from terrestrial natural radioactivity, reveal the amount of uranium and thorium in the Earth and set limits on the residual proportion of primordial energy. Comparison of the flux measured at large underground neutrino experiments with geologically informed predictions of geoneutrino emission from the crust provide the critical test needed to define the mantle's radiogenic power. Measurement at an oceanic location, distant from nuclear reactors and continental crust, would best reveal the mantle flux, however, no such experiment is anticipated. We predict the geoneutrino flux at the site of the Jinping Neutrino Experiment (Sichuan, China). Within 8 years, the combination of existing data and measurements from soon to come experiments, including Jinping, will exclude end-member models at the 1σ level, define the mantle's radiogenic contribution to the surface heat loss, set limits on the composition of the silicate Earth, and provide significant parameter bounds for models defining the mode of mantle convection.

Modeling mantle convection in the spherical annulus

NASA Astrophysics Data System (ADS)

Hernlund, John W.; Tackley, Paul J.

2008-12-01

Most methods for modeling mantle convection in a two-dimensional (2D) circular annular domain suffer from innate shortcomings in their ability to capture several characteristics of the spherical shell geometry of planetary mantles. While methods such as rescaling the inner and outer radius to reduce anomalous effects in a 2D polar cylindrical coordinate system have been introduced and widely implemented, such fixes may have other drawbacks that adversely affect the outcome of some kinds of mantle convection studies. Here we propose a new approach that we term the "spherical annulus," which is a 2D slice that bisects the spherical shell and is quantitatively formulated at the equator of a spherical polar coordinate system after neglecting terms in the governing equations related to variations in latitude. Spherical scaling is retained in this approximation since the Jacobian function remains proportional to the square of the radius. We present example calculations to show that the behavior of convection in the spherical annulus compares favorably against calculations performed in other 2D annular domains when measured relative to those in a fully three-dimensional (3D) spherical shell.

Cold cratonic roots and thermal blankets: How continents affect mantle convection

USGS Publications Warehouse

Trubitsyn, V.P.; Mooney, W.D.; Abbott, D.H.

2003-01-01

Two-dimensional convection models with moving continents show that continents profoundly affect the pattern of mantle convection. If the continents are wider than the wavelength of the convection cells (???3000 km, the thickness of the mantle), they cause neighboring deep mantle thermal upwellings to coalesce into a single focused upwelling. This focused upwelling zone will have a potential temperature anomaly of about 200??C, much higher than the 100??C temperature anomaly of upwelling zones generated beneath typical oceanic lithosphere. Extensive high-temperature melts (including flood basalts and late potassic granites) will be produced, and the excess temperature anomaly will induce continental uplift (as revealed in sea level changes) and the eventual breakup of the supercontinent. The mantle thermal anomaly will persist for several hundred million years after such a breakup. In contrast, small continental blocks (<1000 km diameter) do not induce focused mantle upwelling zones. Instead, small continental blocks are dragged to mantle downwelling zones, where they spend most of their time, and will migrate laterally with the downwelling. As a result of sitting over relatively cold mantle (downwellings), small continental blocks are favored to keep their cratonic roots. This may explain the long-term survival of small cratonic blocks (e.g., the Yilgarn and Pilbara cratons of western Australia, and the West African craton). The optimum size for long-term stability of a continental block is <3000 km. These results show that continents profoundly affect the pattern of mantle convection. These effects are illustrated in terms of the timing and history of supercontinent breakup, the production of high-temperature melts, and sea level changes. Such two-dimensional calculations can be further refined and tested by three-dimensional numerical simulations of mantle convection with moving continental and oceanic plates.

Subduction History and the Evolution of Earth's Lower Mantle

NASA Astrophysics Data System (ADS)

Bull, Abigail; Shephard, Grace; Torsvik, Trond

2016-04-01

Understanding the complex structure, dynamics and evolution of the deep mantle is a fundamental goal in solid Earth geophysics. Close to the core-mantle boundary, seismic images reveal a mantle characterised by (1) higher than average shear wave speeds beneath Asia and encircling the Pacific, consistent with sub ducting lithosphere beneath regions of ancient subduction, and (2) large regions of anomalously low seismic wavespeeds beneath Africa and the Central Pacific. The anomalously slow areas are often referred to as Large Low Shear Velocity Provinces (LLSVPs) due to the reduced velocity of seismic waves passing through them. The origin, composition and long-term evolution of the LLSVPs remain enigmatic. Geochemical inferences of multiple chemical reservoirs at depth, strong seismic contrasts, increased density, and an anticorrelation of shear wave velocity to bulk sound velocity in the anomalous regions imply that heterogeneities in both temperature and composition may be required to explain the seismic observations. Consequently, heterogeneous mantle models place the anomalies into the context of thermochemical piles, characterised by an anomalous component whose intrinsic density is a few percent higher relative to that of the surrounding mantle. Several hypotheses have arisen to explain the LLSVPs in the context of large-scale mantle convection. One end member scenario suggests that the LLSVPs are relatively mobile features over short timescales and thus are strongly affected by supercontinent cycles and Earth's plate motion history. In this scenario, the African LLSVP formed as a result of return flow in the mantle due to circum-Pangean subduction (~240 Ma), contrasting a much older Pacific LLSVP, which may be linked to the Rodinia supercontinent and is implied to have remained largely unchanged since Rodinian breakup (~750-700 Ma). This propounds that Earth's plate motion history plays a controlling role in LLSVP development, suggesting that the location

Numerical study of the origin and stability of chemically distinct reservoirs deep in Earth's mantle

NASA Astrophysics Data System (ADS)

van Thienen, P.; van Summeren, J.; van der Hilst, R. D.; van den Berg, A. P.; Vlaar, N. J.

Seismic tomography is providing mounting evidence for large scale compositional heterogeneity deep in Earth's mantle; also, the diverse geochemical and isotopic signatures observed in oceanic basalts suggest that the mantle is not chemically homogeneous. Isotopic studies on Archean rocks indicate that mantle inhomogeneity may have existed for most of the Earth's history. One important component may be recycled oceanic crust, residing at the base of the mantle. We investigate, by numerical modeling, if such reservoirs may have been formed in the early Earth, before plate tectonics (and subduction) were possible, and how they have survived—and evolved—since then. During Earth's early evolution, thick basaltic crust may have sunk episodically into the mantle in short but vigorous diapiric resurfacing events. These sections of crust may have resided at the base of the mantle for very long times. Entrainment of material from the enriched reservoirs thus produced may account for enriched mantle and high-μ signatures in oceanic basalts, whereas deep subduction events may have shaped and replenished deep mantle reservoirs. Our modeling shows that (1) convective instabilities and resurfacing may have produced deep enriched mantle reservoirs before the era of plate tectonics; (2) such formation is qualitatively consistent with the geochemical record, which shows multiple distinct ocean island basalt sources; and (3) reservoirs thus produced may be stable for billions of years.

Interaction between Edge-Driven Convection and Mantle Plumes

NASA Astrophysics Data System (ADS)

Manjón-Cabeza Córdoba, A.; Ballmer, M.

2017-12-01

Intraplate volcanism can occur in a variety of geodynamic settings. Its characteristics can inform about the underlying mantle dynamics. A non-negligible number of intraplate oceanic volcanoes are located close to continental shelves (e.g. Bermuda, Canary Islands, Cape Verde…). In these regions, any putative plumes would interact with Edge-Driven Convection (EDC), a mode of Small-Scale Convection that is triggered along steps of lithospheric thickness. We have systematically explored 2-D geodynamic models of EDC, varying e.g. the viscosity of the mantle, geometry of the edge, potential temperature, etc. In addition, we study the influence of a mantle plume with variable excess temperature and buoyancy flux at a given distance to the edge. The mantle-convection code is coupled with a new melting parameterization that considers the depletion effect on productivity. We apply this parameterization not only to predict the extent of melting for a given lithology, but also the major-element composition of extracted melts for comparison with geochemical data. Results show that the first EDC upwellings are always localized in the oceanic domain at a distance from the continental margin that depends on mantle viscosity. The initial geometry of the edge does not have a significant influence on the "steady-state" shape of EDC. Depending on the distance of the plume from the edge and plume vigor, the plume is either deflected or enhanced by EDC. The mix of materials that melts in the mantle, as well as the amount of melting, is controlled by the interaction of the plume with EDC (e.g., with melting restricted to fertile heterogeneities in the end-member EDC case). Because several model parameters affect this interaction and related melting, a joint analysis of major-element and trace-element composition of hotspot lavas is required to constrain mantle processes.

Searching for Hysteresis in Models of Mantle Convection with Grain-Damage

NASA Astrophysics Data System (ADS)

Lamichhane, R.; Foley, B. J.

2017-12-01

The mode of surface tectonics on terrestrial planets is determined by whether mantle convective forces are capable of forming weak zones of localized deformation in the lithosphere, which act as plate boundaries. If plate boundaries can form then a plate tectonic mode develops, and if not convection will be in the stagnant lid regime. Episodic subduction or sluggish lid convection are also possible in between the nominal plate tectonic and stagnant lid regimes. Plate boundary formation is largely a function of the state of the mantle, e.g. mantle temperature or surface temperature, and how these conditions influence both mantle convection and the mantle rheology's propensity for forming weak, localized plate boundaries. However, a planet's tectonic mode also influences whether plate boundaries can form, as the driving forces for plate boundary formation (e.g. stress and viscous dissipation) are different in a plate tectonic versus stagnant lid regime. As a result, tectonic mode can display hysteresis, where convection under otherwise identical conditions can reach different final states as a result of the initial regime of convection. Previous work has explored this effect in pseudoplastic models, finding that it is more difficult to initiate plate tectonics starting from a stagnant lid state than it is to sustain plate tectonics when already in a mobile lid regime, because convective stresses in the lithosphere are lower in a stagnant lid regime than in a plate tectonic regime. However, whether and to what extent such hysteresis is displayed when alternative rheological models for lithospheric shear localization are used is unknown. In particular, grainsize reduction is commonly hypothesized to be a primary cause of shear localization and plate boundary formation. We use new models of mantle convection with grain-size evolution to determine how the initial mode of surface tectonics influences the final convective regime reached when convection reaches statistical

On retrodictions of global mantle flow with assimilated surface velocities

NASA Astrophysics Data System (ADS)

Colli, Lorenzo; Bunge, Hans-Peter; Schuberth, Bernhard S. A.

2016-04-01

Modeling past states of Earth's mantle and relating them to geologic observations such as continental-scale uplift and subsidence is an effective method for testing mantle convection models. However, mantle convection is chaotic and two identical mantle models initialized with slightly different temperature fields diverge exponentially in time until they become uncorrelated, thus limiting retrodictions (i.e., reconstructions of past states of Earth's mantle obtained using present information) to the recent past. We show with 3-D spherical mantle convection models that retrodictions of mantle flow can be extended significantly if knowledge of the surface velocity field is available. Assimilating surface velocities produces in some cases negative Lyapunov times (i.e., e-folding times), implying that even a severely perturbed initial condition may evolve toward the reference state. A history of the surface velocity field for Earth can be obtained from past plate motion reconstructions for time periods of a mantle overturn, suggesting that mantle flow can be reconstructed over comparable times.

On retrodictions of global mantle flow with assimilated surface velocities

NASA Astrophysics Data System (ADS)

Colli, Lorenzo; Bunge, Hans-Peter; Schuberth, Bernhard S. A.

2015-10-01

Modeling past states of Earth's mantle and relating them to geologic observations such as continental-scale uplift and subsidence is an effective method for testing mantle convection models. However, mantle convection is chaotic and two identical mantle models initialized with slightly different temperature fields diverge exponentially in time until they become uncorrelated, thus limiting retrodictions (i.e., reconstructions of past states of Earth's mantle obtained using present information) to the recent past. We show with 3-D spherical mantle convection models that retrodictions of mantle flow can be extended significantly if knowledge of the surface velocity field is available. Assimilating surface velocities produces in some cases negative Lyapunov times (i.e., e-folding times), implying that even a severely perturbed initial condition may evolve toward the reference state. A history of the surface velocity field for Earth can be obtained from past plate motion reconstructions for time periods of a mantle overturn, suggesting that mantle flow can be reconstructed over comparable times.

Self-Consistent Generation of Continental Crust in Global Mantle Convection Models

NASA Astrophysics Data System (ADS)

Jain, C.; Rozel, A. B.; Tackley, P.

2016-12-01

Numerical modelling commonly shows that mantle convection and continents have strong feedbacks on each other (Philips and Coltice, JGR 2010; Heron and Lowman, JGR 2014), but the continents are always inserted a priori while basaltic (oceanic) crust is generated self-consistently in such models (Rolf et al., EPSL 2012). We aim to implement self-consistent generation of continental crust in global models of mantle convection using StagYY (Tackley, PEPI 2008). The silica-rich continental crust appears to have been formed by fractional melting and crystallization in episodes of relatively rapid growth from late Archaean to late Proterozoic eras (3-1 Ga) (Hawkesworth & Kemp, Nature 2006). It takes several stages of differentiation to generate continental crust. First, the basaltic magma is extracted from the pyrolitic mantle. Second, it goes through eclogitic transformation and then partially melts to form Na-rich Tonalite-Trondhjemite-Granodiorite (TTG) which rise to form proto-continents (Rudnick, Nature 1995; Herzberg & Rudnick, Lithos 2012). TTGs dominate the grey gneiss complexes which make up most of the continental crust. Based on the melting conditions proposed by Moyen (Lithos, 2011), we parameterize TTG formation and henceforth, the continental crust. Continental crust can also be destroyed by subduction or delamination. We will investigate continental growth and destruction history in the models spanning the age of the Earth.

Self-consistent generation of continental crust in global mantle convection models

NASA Astrophysics Data System (ADS)

Jain, Charitra; Rozel, Antoine; Tackley, Paul

2017-04-01

Numerical modeling commonly shows that mantle convection and continents have strong feedbacks on each other (Philips and Coltice, JGR 2010; Heron and Lowman, JGR 2014), but the continents are always inserted a priori while basaltic (oceanic) crust is generated self-consistently in such models (Rolf et al., EPSL 2012). We aim to implement self-consistent generation of continental crust in global models of mantle convection using StagYY (Tackley, PEPI 2008). The silica-rich continental crust appears to have been formed by fractional melting and crystallization in episodes of relatively rapid growth from late Archean to late Proterozoic eras (3-1 Ga) (Hawkesworth & Kemp, Nature 2006). It takes several stages of differentiation to generate continental crust. First, the basaltic magma is extracted from the pyrolitic mantle. Second, it goes through eclogitic transformation and then partially melts to form Na-rich Tonalite-Trondhjemite-Granodiorite (TTG) which rise to form proto-continents (Rudnick, Nature 1995; Herzberg & Rudnick, Lithos 2012). TTGs dominate the grey gneiss complexes which make up most of the continental crust. Based on the melting conditions proposed by Moyen (Lithos, 2011), we parameterize TTG formation and henceforth, the continental crust. Continental crust can also be destroyed by subduction or delamination. We will investigate continental growth and destruction history in the models spanning the age of the Earth.

ISS COLUMBUS laboratory experiment `GeoFlow I and II' -fluid physics research in microgravity environment to study convection phenomena inside deep Earth and mantle

NASA Astrophysics Data System (ADS)

Futterer, Birgit; Egbers, Christoph; Chossat, Pascal; Hollerbach, Rainer; Breuer, Doris; Feudel, Fred; Mutabazi, Innocent; Tuckerman, Laurette

Overall driving mechanism of flow in inner Earth is convection in its gravitational buoyancy field. A lot of effort has been involved in theoretical prediction and numerical simulation of both the geodynamo, which is maintained by convection, and mantle convection, which is the main cause for plate tectonics. Especially resolution of convective patterns and heat transfer mechanisms has been in focus to reach the real, highly turbulent conditions inside Earth. To study specific phenomena experimentally different approaches has been observed, against the background of magneto-hydrodynamic but also on the pure hydrodynamic physics of fluids. With the experiment `GeoFlow' (Geophysical Flow Simulation) instability and transition of convection in spherical shells under the influence of central-symmetry buoyancy force field are traced for a wide range of rotation regimes within the limits between non-rotating and rapid rotating spheres. The special set-up of high voltage potential between inner and outer sphere and use of a dielectric fluid as working fluid induce an electro-hydrodynamic force, which is comparable to gravitational buoyancy force inside Earth. To reduce overall gravity in a laboratory this technique requires microgravity conditions. The `GeoFlow I' experiment was accomplished on International Space Station's module COLUM-BUS inside Fluid Science Laboratory FSL und supported by EADS Astrium, Friedrichshafen, User Support und Operations Centre E-USOC in Madrid, Microgravity Advanced Research and Support Centre MARS in Naples, as well as COLUMBUS Control Center COL-CC Munich. Running from August 2008 until January 2009 it delivered 100.000 images from FSL's optical diagnostics module; here more precisely the Wollaston shearing interferometry was used. Here we present the experimental alignment with numerical prediction for the non-rotating and rapid rotation case. The non-rotating case is characterized by a co-existence of several stationary supercritical

Mantle dynamics and seismic tomography

PubMed Central

Tanimoto, Toshiro; Lay, Thorne

2000-01-01

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

Boundary-modulated Thermal Convection Model in the Mantle

NASA Astrophysics Data System (ADS)

Kurita, K.; Kumagai, I.

2008-12-01

Analog experiments have played an important role in the constructing ideas of mantle dynamics. The series of experiments by H. Ramberg is one of the successful examples. Recently, however the realm of the analog experiments seems to be overwhelmed by steady progress of computer simulations. Is there still room for the analog experiments? This might be a main and hidden subject of this session. Here we propose a working hypothesis how the convecting mantle behaves based on the analog experiments in the system of viscous fluid and particles. The essential part is the interaction of convecting flow with heterogeneities existing in the boundaries. It is proposed the preexisting topographical heterogeneity in the boundary could control the flow pattern of convecting fluid. If this kind of heterogeneity can be formed as a consequence of convective motion and mobilized by the flow, the convection also can control the heterogeneity. We can expect interactions in two ways, by which the system behaves in a self-organize fashion. To explore the mutual interactions between convection flow and heterogeneity the system of viscous fluid and particles with slightly higher density is selected as 2D Rayleigh-Benard type convection. The basic structure consists of a basal particulate layer where permeable convection transports heat and an upper viscous fluid layer. By reducing the magnitude of the density difference the convective flow can mobilize the particles and can erode the basal layer. The condition of this erosion can be identified in the phase diagram of the particle Shields"f and the Rayleigh numbers. At Ra greater than 107 the convection style drastically changed before and after the erosion. Before the erosion where the flat interface of the boundary is maintained small scaled turbulent convection pattern is dominant. After the erosion where the interface becomes bumpy the large scale convective motion is observed. The structure is coherent to that of the boundary. This

Assessing the role of slab rheology in coupled plate-mantle convection models

NASA Astrophysics Data System (ADS)

Bello, Léa; Coltice, Nicolas; Tackley, Paul J.; Dietmar Müller, R.; Cannon, John

2015-11-01

Reconstructing the 3D structure of the Earth's mantle has been a challenge for geodynamicists for about 40 yr. Although numerical models and computational capabilities have substantially progressed, parameterizations used for modeling convection forced by plate motions are far from being Earth-like. Among the set of parameters, rheology is fundamental because it defines in a non-linear way the dynamics of slabs and plumes, and the organization of lithosphere deformation. In this study, we evaluate the role of the temperature dependence of viscosity (variations up to 6 orders of magnitude) and the importance of pseudo-plasticity on reconstructing slab evolution in 3D spherical models of convection driven by plate history models. Pseudo-plasticity, which produces plate-like behavior in convection models, allows a consistent coupling between imposed plate motions and global convection, which is not possible with temperature-dependent viscosity alone. Using test case models, we show that increasing temperature dependence of viscosity enhances vertical and lateral coherence of slabs, but leads to unrealistic slab morphologies for large viscosity contrasts. Introducing pseudo-plasticity partially solves this issue, producing thin laterally and vertically more continuous slabs, and flat subduction where trench retreat is fast. We evaluate the differences between convection reconstructions employing different viscosity laws to be very large, and similar to the differences between two models with the same rheology but using two different plate histories or initial conditions.

Towards a New Framework for Interpreting Relations Between Mantle Dynamics and Processes at the Earth's Surface: A Case Study Involving the Deccan Traps

NASA Astrophysics Data System (ADS)

Glisovic, P.; Forte, A. M.

2017-12-01

An outstanding challenge in modern geodynamics is the utilization of mantle convection models and geophysical data to successfully explain geological events and processes that alter Earth's biosphere, climate, and surface. A key challenge in this modelling is the determination of the initial (and unknown) configuration of mantle heterogeneity in the geological past. The first step in addressing this challenge is recognizing that seismic tomography is our most powerful tool for mapping the present-day, internal structure of the mantle. We, therefore, implemented a new back-and-forth iterative method for time-reversed, tomography-based convection modelling to reconstruct Earth's internal 3-D structure and dynamics over the Cenozoic [Glisovic & Forte 2016 (JGR)]. This backward convection modelling also includes another key input - the depth variation of mantle viscosity inferred from joint inversions of the global convection-related observables and a suite of glacial isostatic adjustments (GIA) data [Mitrovica & Forte 2004 (EPSL), Forte et al. 2010 (EPSL)]. This state-of-the-art, time-reversed convection model is able to show that massive outpourings of basalt in west-central India, known as the Deccan Traps, about 65 million years ago can be directly linked to the presence of two different deep-mantle hotspots: Réunion and Comores [Glisovic & Forte 2017 (Science)]. This work constitutes case study showing how time-reversed convection modelling provides a new framework for interpreting the relations between mantle dynamics and changing paleogeography and it provides a roadmap for a new series of studies that will elucidate these linkages.

Satellite Detection of the Convection Generated Stresses in Earth

NASA Technical Reports Server (NTRS)

Liu, Han-Shou; Kolenkiewicz, Ronald; Li, Jin-Ling; Chen, Jiz-Hong

2003-01-01

We review research developments on satellite detection of the convection generated stresses in the Earth for seismic hazard assessment and Earth resource survey. Particular emphasis is laid upon recent progress and results of stress calculations from which the origin and evolution of the tectonic features on Earth's surface can be scientifically addressed. An important aspect of the recent research development in tectonic stresses relative to earthquakes is the implications for earthquake forecasting and prediction. We have demonstrated that earthquakes occur on the ring of fire around the Pacific in response to the tectonic stresses induced by mantle convection. We propose a systematic global assessment of the seismic hazard based on variations of tectonic stresses in the Earth as observed by satellites. This space geodynamic approach for assessing the seismic hazard is unique in that it can pinpoint the triggering stresses for large earthquakes without ambiguities of geological structures, fault geometries, and other tectonic properties. Also, it is distinct from the probabilistic seismic hazard assessment models in the literature, which are based only on extrapolations of available earthquake data.

The Earth's mantle in a microwave oven: thermal convection driven by a heterogeneous distribution of heat sources

NASA Astrophysics Data System (ADS)

Fourel, Loïc; Limare, Angela; Jaupart, Claude; Surducan, Emanoil; Farnetani, Cinzia G.; Kaminski, Edouard C.; Neamtu, Camelia; Surducan, Vasile

2017-08-01

Convective motions in silicate planets are largely driven by internal heat sources and secular cooling. The exact amount and distribution of heat sources in the Earth are poorly constrained and the latter is likely to change with time due to mixing and to the deformation of boundaries that separate different reservoirs. To improve our understanding of planetary-scale convection in these conditions, we have designed a new laboratory setup allowing a large range of heat source distributions. We illustrate the potential of our new technique with a study of an initially stratified fluid involving two layers with different physical properties and internal heat production rates. A modified microwave oven is used to generate a uniform radiation propagating through the fluids. Experimental fluids are solutions of hydroxyethyl cellulose and salt in water, such that salt increases both the density and the volumetric heating rate. We determine temperature and composition fields in 3D with non-invasive techniques. Two fluorescent dyes are used to determine temperature. A Nd:YAG planar laser beam excites fluorescence, and an optical system, involving a beam splitter and a set of colour filters, captures the fluorescence intensity distribution on two separate spectral bands. The ratio between the two intensities provides an instantaneous determination of temperature with an uncertainty of 5% (typically 1K). We quantify mixing processes by precisely tracking the interfaces separating the two fluids. These novel techniques allow new insights on the generation, morphology and evolution of large-scale heterogeneities in the Earth's lower mantle.

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)

On the temporal evolution of long-wavelength mantle structure of the Earth since the early Paleozoic

NASA Astrophysics Data System (ADS)

Zhong, Shijie; Rudolph, Maxwell L.

2015-05-01

The seismic structure of the Earth's lower mantle is characterized by a dominantly degree-2 pattern with the African and Pacific large low shear velocity provinces (i.e., LLSVP) that are separated by circum-Pacific seismically fast anomalies. It is important to understand the origin of such a degree-2 mantle structure and its temporal evolution. In this study, we investigated the effects of plate motion history and mantle viscosity on the temporal evolution of the lower mantle structure since the early Paleozoic by formulating 3-D spherical shell models of thermochemical convection. For convection models with realistic mantle viscosity and no initial structure, it takes about ˜50 Myr to develop dominantly degree-2 lower mantle structure using the published plate motion models for the last either 120 Ma or 250 Ma. However, it takes longer time to develop the mantle structure for more viscous mantle. While the circum-Pangea subduction in plate motion history models promotes the formation of degree-2 mantle structure, the published pre-Pangea plate motions before 330 Ma produce relatively cold lower mantle in the African hemisphere and significant degree-1 structure in the early Pangea (˜300 Ma) or later times, even if the lower mantle has an initially degree-2 structure and a viscosity as high as 1023 Pas. This suggests that the African LLSVP may not be stationary since the early Paleozoic. With the published plate motion models and lower mantle viscosity of 1022 Pas, our mantle convection models suggest that the present-day degree-2 mantle structure may have largely been formed by ˜200 Ma.

An 1-2-1 Cyclic Model for the Evolution of Mantle Structure

NASA Astrophysics Data System (ADS)

Zhong, S.; Zhang, N.

2006-12-01

The present-day Earth`s mantle is predominated by long-wavelength structures including circum-Pacific subducted slabs and Africa and Pacific super-plumes. These long-wavelength structures are largely controlled by the history of plate tectonic motion. Although it dictates the evolution of mantle structure, global plate tectonic history prior to 120 Ma is poorly constrained except for continental motions that can be reliably traced back to >1 Ga. An important observation of continental motions in the last 1 Ga is the two episodes of formation and breakup of super-continents Pangea and Rodinia. We formulated 3D global models of mantle convection with temperature- and depth-dependent viscosity to study the formation of mantle structure. We found that for the upper mantle with 30 times smaller viscosity than the lower mantle, in the absence of continents, mantle convection is characterized by a hemispherically asymmetric structure in which one hemisphere is largely upwellings, while the other hemisphere contains downwellings (i.e., degree-1 convection). This is the first study in which degree-1 mantle convection is observed in mobile-lid/plate-tectonic convection regime at high Rayleigh number. This result suggests that degree-1 convection is a dynamically preferred state for the Earth`s mantle. We suggest that the evolution of mantle structure is controlled by a cyclic process of formation and breakdown of degree-1 convection modulated strongly by continents. The formation and breakup of supercontinents are surface manifestation of this cyclic process. During the degree-1 convection state, the upwellings in one hemisphere push all continents into the other hemisphere with the downwellings to form a supercontinent. The non-subducting nature of continents dictates that subduction in the downwelling hemisphere occurs along the edge of the supercontinent upon its formation. The insulating effect of a supercontinent and return flow from the circum-supercontinent subduction

Spontaneous development of arcuate single-sided subduction in global 3-D mantle convection models with a free surface

NASA Astrophysics Data System (ADS)

Crameri, Fabio; Tackley, Paul

2014-05-01

plate surface and weak oceanic crust produce single-sided subduction on Earth, Geophys. Res. Lett., 39(3), L03,306. Crameri, F., H. Schmeling, G. J. Golabek, T. Duretz, R. Orendt, S. J. H. Buiter, D. A. May, B. J. P. Kaus, T. V. Gerya, and P. J. Tackley (2012b), A comparison of numerical surface topography calculations in geodynamic modelling: an evaluation of the 'sticky air' method, Geophys. J. Int., 189(1), 38-54. Matsumoto, T., and Y. Tomoda (1983), Numerical simulation of the initiation of subduction at the fracture zone, J. Phys. Earth, 31(3), 183-194. Schmeling, H., A. Babeyko, A. Enns, C. Faccenna, F. Funiciello, T. Gerya, G. Golabek, S. Grigull, B. Kaus, G. Morra, S. Schmalholz, and J. van Hunen (2008), A benchmark comparison of spontaneous subduction models-Towards a free surface, Phys. Earth Planet. Int., 171(1-4), 198-223. Tackley, P. J. (2008), Modelling compressible mantle convection with large viscosity contrasts in a three-dimensional spherical shell using the yin-yang grid, Phys. Earth Planet. Int., 171(1-4), 7-18.

Numerical Modeling of Deep Mantle Flow: Thermochemical Convection and Entrainment

NASA Astrophysics Data System (ADS)

Mulyukova, Elvira; Steinberger, Bernhard; Dabrowski, Marcin; Sobolev, Stephan

2013-04-01

One of the most robust results from tomographic studies is the existence of two antipodally located Large Low Shear Velocity Provinces (LLSVPs) at the base of the mantle, which appear to be chemically denser than the ambient mantle. Results from reconstruction studies (Torsvik et al., 2006) infer that the LLSVPs are stable, long-lived, and are sampled by deep mantle plumes that rise predominantly from their margins. The origin of the dense material is debated, but generally falls within three categories: (i) a primitive layer that formed during magma ocean crystallization, (ii) accumulation of a dense eclogitic component from the recycled oceanic crust, and (iii) outer core material leaking into the lower mantle. A dense layer underlying a less dense ambient mantle is gravitationally stable. However, the flow due to thermal density variations, i.e. hot rising plumes and cold downwelling slabs, may deform the layer into piles with higher topography. Further deformation may lead to entrainment of the dense layer, its mixing with the ambient material, and even complete homogenisation with the rest of the mantle. The amount of the anomalous LLSVP-material that gets entrained into the rising plumes poses a constraint on the survival time of the LLSVPs, as well as on the plume buoyancy, on the lithospheric uplift associated with plume interaction and geochemical signature of the erupted lavas observed at the Earth's surface. Recent estimates for the plume responsible for the formation of the Siberian Flood Basalts give about 15% of entrained dense recycled oceanic crust, which made the hot mantle plume almost neutrally buoyant (Sobolev et al., 2011). In this numerical study we investigate the mechanics of entrainment of a dense basal layer by convective mantle flow. We observe that the types of flow that promote entrainment of the dense layer are (i) upwelling of the dense layer when it gets heated enough to overcome its stabilizing chemical density anomaly, (ii

Passive margins getting squeezed in the mantle convection vice

NASA Astrophysics Data System (ADS)

Yamato, Philippe; Husson, Laurent; Becker, Thorsten W.; Pedoja, Kevin

2013-12-01

margins often exhibit uplift, exhumation, and tectonic inversion. We speculate that the compression in the lithosphere gradually increased during the Cenozoic, as seen in the number of mountain belts found at active margins during that period. Less clear is how that compression increase affects passive margins. In order to address this issue, we design a 2-D viscous numerical model wherein a lithospheric plate rests above a weaker mantle. It is driven by a mantle conveyor belt, alternatively excited by a lateral downwelling on one side, an upwelling on the other side, or both simultaneously. The lateral edges of the plate are either free or fixed, representing the cases of free convergence, and collision (or slab anchoring), respectively. This distinction changes the upper mechanical boundary condition for mantle circulation and thus, the stress field. Between these two regimes, the flow pattern transiently evolves from a free-slip convection mode toward a no-slip boundary condition above the upper mantle. In the second case, the lithosphere is highly stressed horizontally and deforms. For a constant total driving force, compression increases drastically at passive margins if upwellings are active. Conversely, if downwellings alone are activated, compression occurs at short distances from the trench and extension prevails elsewhere. These results are supported by Earth-like models that reveal the same pattern, where active upwellings are required to excite passive margins compression. Our results substantiate the idea that compression at passive margins is in response to the underlying mantle flow that is increasingly resisted by the Cenozoic collisions.

Numerical studies on convective stability and flow pattern in three-dimensional spherical mantle of terrestrial planets

NASA Astrophysics Data System (ADS)

Yanagisawa, Takatoshi; Kameyama, Masanori; Ogawa, Masaki

2016-09-01

We explore thermal convection of a fluid with a temperature-dependent viscosity in a basally heated 3-D spherical shell using linear stability analyses and numerical experiments, while considering the application of our results to terrestrial planets. The inner to outer radius ratio of the shell f assumed in the linear stability analyses is in the range of 0.11-0.88. The critical Rayleigh number Rc for the onset of thermal convection decreases by two orders of magnitude as f increases from 0.11 to 0.88, when the viscosity depends sensitively on the temperature, as is the case for real mantle materials. Numerical simulations carried out in the range of f = 0.11-0.55 show that a thermal boundary layer (TBL) develops both along the surface and bottom boundaries to induce cold and hot plumes, respectively, when f is 0.33 or larger. However, for smaller f values, a TBL develops only on the bottom boundary. Convection occurs in the stagnant-lid regime where the root mean square velocity on the surface boundary is less than 1 per cent of its maximum at depth, when the ratio of the viscosity at the surface boundary to that at the bottom boundary exceeds a threshold that depends on f. The threshold decreases from 106.5 at f = 0.11 to 104 at f = 0.55. If the viscosity at the base of the convecting mantle is 1020-1021 Pa s, the Rayleigh number exceeds Rc for Mars, Venus and the Earth, but does not for the Moon and Mercury; convection is unlikely to occur in the latter planets unless the mantle viscosity is much lower than 1020 Pa s and/or the mantle contains a strong internal heat source.

On global gravity anomalies and two-scale mantle convection

NASA Technical Reports Server (NTRS)

Marsh, B. D.; Marsh, J. G.

1976-01-01

The two-scale model of mantle convection developed by Richter and Parsons (1975) predicts that if the depth of the convective layer is about 600 km, then for a plate moving at 10 cm/yr, longitudinal convective rolls will be produced in about 50 million years, and the strike of these rolls indicates the direction of motion of the plate relative to the upper mantle. The paper tests these predictions by examining a new global free air gravity model complete to the 30th degree and order. The free air gravity map developed shows a series of linear positive and negative anomalies (with transverse wavelengths of about 2000 km) spanning the Pacific Ocean, crossing the Pacific rise and striking parallel to the Hawaiian seamounts. It is suggested that the pattern of these anomalies may indicate the presence of longitudinal convective rolls beneath the Pacific plates, a result which tends to support the predictions of Richter and Parsons.

Self-Consistent Generation of Primordial Continental Crust in Global Mantle Convection Models

NASA Astrophysics Data System (ADS)

Jain, C.; Rozel, A.; Tackley, P. J.

2017-12-01

We present the generation of primordial continental crust (TTG rocks) using self-consistent and evolutionary thermochemical mantle convection models (Tackley, PEPI 2008). Numerical modelling commonly shows that mantle convection and continents have strong feedbacks on each other. However in most studies, continents are inserted a priori while basaltic (oceanic) crust is generated self-consistently in some models (Lourenco et al., EPSL 2016). Formation of primordial continental crust happened by fractional melting and crystallisation in episodes of relatively rapid growth from late Archean to late Proterozoic eras (3-1 Ga) (Hawkesworth & Kemp, Nature 2006) and it has also been linked to the onset of plate tectonics around 3 Ga. It takes several stages of differentiation to generate Tonalite-Trondhjemite-Granodiorite (TTG) rocks or proto-continents. First, the basaltic magma is extracted from the pyrolitic mantle which is both erupted at the surface and intruded at the base of the crust. Second, it goes through eclogitic transformation and then partially melts to form TTGs (Rudnick, Nature 1995; Herzberg & Rudnick, Lithos 2012). TTGs account for the majority of the Archean continental crust. Based on the melting conditions proposed by Moyen (Lithos 2011), the feasibility of generating TTG rocks in numerical simulations has already been demonstrated by Rozel et al. (Nature, 2017). Here, we have developed the code further by parameterising TTG formation. We vary the ratio of intrusive (plutonic) and extrusive (volcanic) magmatism (Crisp, Volcanol. Geotherm. 1984) to study the relative volumes of three petrological TTG compositions as reported from field data (Moyen, Lithos 2011). Furthermore, we systematically vary parameters such as friction coefficient, initial core temperature and composition-dependent viscosity to investigate the global tectonic regime of early Earth. Continental crust can also be destroyed by subduction or delamination. We will investigate

Mantle convection pattern and subcrustal stress field under South America

NASA Technical Reports Server (NTRS)

Liu, H.-S.

1980-01-01

The tectonic, igneous and metallogenic features of South America are discussed in terms of the crustal deformation associated with stresses due to mantle convection as inferred from the high degree harmonics in the geopotential field. The application of Runcorn's model for the laminar viscous flows in the upper mantle to satellite and gravity data results in a convection pattern which reveals the ascending flows between the descending Nazca plate and the overlying South American plate as well as segments of the descending Nazca plate beneath South America. The arc volcanism in South America is shown apparently to be related to the upwelling of high-temperature material induced by the subduction of the Nazca plate, with the South American basin systems associated with downwelling mantle flows. The resulting tensional stress fields are shown to be regions of structural kinship characterized by major concentrations of ore deposits and related to the cordillera, shield and igneous systems and the upward Andean movements. It is suggested that the upwelling convection flows in the upper mantle, coupled with crustal tension, have provided an uplift mechanism which has forced the hydrothermal systems in the basement rocks to the surface.

Continental Growth and Mantle Hydration as Earth System Feedback Cycles and possible Effects of the Biosphere

NASA Astrophysics Data System (ADS)

Höning, D.; Spohn, T.

2016-12-01

The evolution of Earth is charcterized by intertwined feedback cycles. We focus on two feedback cycles that include the mantle water budget and the continental crust and study possible effects of the Earth's biosphere. The first feedback loop includes cycling of water into the mantle at subduction zones and outgassing at volcanic chains and mid-ocean ridges. Water will reduce the viscosity of mantle rock, and therefore the speed of mantle convection and plate subduction will increase with the mantle water concentration, eventually enhancing the rates of mantle water regassing and outgassing. A second feedback loop includes the production and erosion of continental crust. Continents grow by volcanism above subduction zones, whose total length is determined by the total area of the continents. Furthermore, the erosion rate of the continents is proportional to the total surface area of continental crust. The rate of sediment subduction affects the rate of transport of water to the mantle and the production rate of new continental crust. We present a model that includes both cycles and show how the system develops stable and unstable fixed points in a plane defined by mantle water concentration and surface are of continents. The stable points represent either an Earth mostly covered by continents and a wet mantle or an Earth mostly covered by oceans with a dry mantle. The presently observed Earth is inbetween these extreme states but the state is intrinsically unstable. We couple the feedback model to a parameterized thermal evolution model. We show how Earth evolved towards its present unstable state. We argue that other feedback cycles such as the carbonate silicate cycle may act to stabilize the present state, however. By enhancing continental weathering and erosion, and eventually the sediment transport into subduction zones, the biosphere impacts both feedback cycles and might play a crucial role in regulating Earth's system and keep continental crust coverage and

Deep mantle roots and continental hypsometry: implications for whole-Earth elemental cycling, long-term climate, and the Cambrian explosion

NASA Astrophysics Data System (ADS)

Lee, C. T.

2016-12-01

Most of Earth's continents today are above sea level, but the presence of thick packages of ancient sediments on top of the stable cores of continents indicates that continents must have been submerged at least once in their past. Elevations of continents are controlled by the interplay between crustal thickness, mantle root thickness and the temperature of the ambient convecting mantle. The history of a continent begins with mountain building through magmatic or tectonic crustal thickening, during which exhumation of deep-seated igneous and metamorphic rocks are highest. Mountain building is followed by a long interval of subsidence as a result of continued, but decreasing erosion and thermal relaxation, the latter in the form of a growing thermal boundary layer. Subsidence is manifest first as a boring interval in which no sedimentary record is preserved, followed by continent-scale submergence wherein sediments are deposited directly on deep-seated igneous/metamorphic basement, generating a major disconformity. The terminal resting elevation of a mature continent, however, is defined by the temperature of the ambient convecting mantle: below sea level when the mantle is hot and above sea level when the mantle is cold. Using thermobarometric constraints on secular cooling of Earth's mantle, our results suggest that Earth, for most of its history, must have been a water world, with regions of land confined to narrow orogenic belts and oceans characterized by deep basins and shallow continental seas, the latter serving as repositories of sediments and key redox-sensitive biological nutrients, such as phosphorous. Cooling of the Earth led to the gradual and irreversible rise of the continents, culminating in rapid emergence, through fits and starts and possible instabilities in climate, between 500-1000 Ma. Such emergence fundamentally altered marine biogeochemical cycling, continental weathering and the global hydrologic cycle, defining the backdrop for the

Passive margins getting squeezed in the mantle convection vice

NASA Astrophysics Data System (ADS)

Yamato, Philippe; Husson, Laurent; Becker, Thorsten W.; Pedoja, Kevin

2014-05-01

Passive margins often exhibit uplift, exhumation and tectonic inversion. We speculate that the compression in the lithosphere gradually increased during the Cenozoic. In the same time, the many mountain belts at active margins that accompany this event seem readily witness this increase. However, how that compression increase affects passive margins remains unclear. In order to address this issue, we design a 2D viscous numerical model wherein a lithospheric plate rests above a weaker mantle. It is driven by a mantle conveyor belt, alternatively excited by a lateral downwelling on one side, an upwelling on the other side, or both simultaneously. The lateral edges of the plate are either free or fixed, representing the cases of free convergence, and collision or slab anchoring, respectively. This distinction changes the upper boundary condition for mantle circulation and, as a consequence, the stress field. Our results show that between these two regimes, the flow pattern transiently evolves from a free-slip convection mode towards a no-slip boundary condition above the upper mantle. In the second case, the lithosphere is highly stressed horizontally and deforms. For an equivalent bulk driving force, compression increases drastically at passive margins provided that upwellings are active. Conversely, if downwellings alone are activated, compression occurs at short distances from the trench and extension prevails elsewhere. These results are supported by Earth-like 3D spherical models that reveal the same pattern, where active upwellings are required to excite passive margins compression. These results support the idea that compression at passive margins, is the response to the underlying mantle flow, that is increasingly resisted by the Cenozoic collisions.

Continental margin subsidence from shallow mantle convection: Example from West Africa

NASA Astrophysics Data System (ADS)

Lodhia, Bhavik Harish; Roberts, Gareth G.; Fraser, Alastair J.; Fishwick, Stewart; Goes, Saskia; Jarvis, Jerry

2018-01-01

Spatial and temporal evolution of the uppermost convecting mantle plays an important role in determining histories of magmatism, uplift, subsidence, erosion and deposition of sedimentary rock. Tomographic studies and mantle flow models suggest that changes in lithospheric thickness can focus convection and destabilize plates. Geologic observations that constrain the processes responsible for onset and evolution of shallow mantle convection are sparse. We integrate seismic, well, gravity, magmatic and tomographic information to determine the history of Neogene-Recent (<23 Ma) upper mantle convection from the Cape Verde swell to West Africa. Residual ocean-age depths of +2 km and oceanic heat flow anomalies of +16 ± 4 mW m-2 are centered on Cape Verde. Residual depths decrease eastward to zero at the fringe of the Mauritania basin. Backstripped wells and mapped seismic data show that 0.4-0.8 km of water-loaded subsidence occurred in a ∼500 × 500 km region centered on the Mauritania basin during the last 23 Ma. Conversion of shear wave velocities into temperature and simple isostatic calculations indicate that asthenospheric temperatures determine bathymetry from Cape Verde to West Africa. Calculated average excess temperatures beneath Cape Verde are > + 100 °C providing ∼103 m of support. Beneath the Mauritania basin average excess temperatures are < - 100 °C drawing down the lithosphere by ∼102 to 103 m. Up- and downwelling mantle has generated a bathymetric gradient of ∼1/300 at a wavelength of ∼103 km during the last ∼23 Ma. Our results suggest that asthenospheric flow away from upwelling mantle can generate downwelling beneath continental margins.

Seismic Heterogeneity Caused by Oceanic Crust Differentiation and Segregation in the Convecting Mantle

NASA Astrophysics Data System (ADS)

Xie, S.; Tackley, P. J.

2003-12-01

This presentation focuses on the seismic signature of mantle heterogeneity associated with crustal differentiation and segregation in the lower mantle. Segregation of subducted oceanic crust above the CMB has often been invoked as a way of explaining the isotopic signature of OIB geochemical endmembers such as HIMU. Here a mantle convection model that includes melting-induced differentiation and plate tectonics is run for billions of years and the resulting thermo-chemical heterogeneity is studied. Statistical diagnostics such as radial correlation functions (Jordan et al., 1993) and spectral heterogeneity maps (Tackley et al., 1994) are used to characterize the observational signature of the thermo-chemical structures and compare them to global seismic tomographic models. In the reference case, crust is denser than the background mantle at the CMB. Due to this density contrast, the crustal material forms a thick and dense layer at the bottom of the mantle, although the layer interface is not sharp as is commonly obtained in models where a layer is inserted a priori. An enormous amount of long-wavelength volumetric heterogeneity is found in the lower mantle. The presence of oceanic crust near the surface also contributes to heterogeneity at the top of the mantle. In cases where the subducted crust is neutrally buoyant or buoyant in the deepest mantle, a large amount of heterogeneity instead exists in the the mid-mantle region, which is not observed in tomographic models of the real Earth. Unlike the reference case with a thick layer at the bottom of the mantle, these cases have heterogeneity right at the CMB, and this strong heterogeneity exists at both short and long wavelength. When comparing to Earth, it appears that models in which dense subducted crust settles into a layer above the CMB are consistent with constraints from seismic tomography; such a layer is therefore a viable location for the storage of geochemical endmembers.

African hot spot volcanism: small-scale convection in the upper mantle beneath cratons.

PubMed

King, S D; Ritsema, J

2000-11-10

Numerical models demonstrate that small-scale convection develops in the upper mantle beneath the transition of thick cratonic lithosphere and thin oceanic lithosphere. These models explain the location and geochemical characteristics of intraplate volcanos on the African and South American plates. They also explain the presence of relatively high seismic shear wave velocities (cold downwellings) in the mantle transition zone beneath the western margin of African cratons and the eastern margin of South American cratons. Small-scale, edge-driven convection is an alternative to plumes for explaining intraplate African and South American hot spot volcanism, and small-scale convection is consistent with mantle downwellings beneath the African and South American lithosphere.

Mantle convection with plates and mobile, faulted plate margins.

PubMed

Zhong, S; Gurnis, M

1995-02-10

A finite-element formulation of faults has been incorporated into time-dependent models of mantle convection with realistic rheology, continents, and phase changes. Realistic tectonic plates naturally form with self-consistent coupling between plate and mantle dynamics. After the initiation of subduction, trenches rapidly roll back with subducted slabs temporarily laid out along the base of the transition zone. After the slabs have penetrated into the lower mantle, the velocity of trench migration decreases markedly. The inhibition of slab penetration into the lower mantle by the 670-kilometer phase change is greatly reduced in these models as compared to models without tectonic plates.

Mantle circulation models with variational data assimilation: Inferring past mantle flow and structure from plate motion histories and seismic tomography

NASA Astrophysics Data System (ADS)

Bunge, Hans-Peter

2002-08-01

Earth's mantle overturns itself about once every 200 Million years (myrs). Prima facie evidence for this overturn is the motion of tectonic plates at the surface of the Earth driving the geologic activity of our planet. Supporting evidence also comes from seismic tomograms of the Earth's interior that reveal the convective currents in remarkable clarity. Much has been learned about the physics of solid state mantle convection over the past two decades aided primarily by sophisticated computer simulations. Such simulations are reaching the threshold of fully resolving the convective system globally. In this talk we will review recent progress in mantle dynamics studies. We will then turn our attention to the fundamental question of whether it is possible to explicitly reconstruct mantle flow back in time. This is a classic problem of history matching, amenable to control theory and data assimilation. The technical advances that make such approach feasible are dramatically increasing compute resources, represented for example through Beowulf clusters, and new observational initiatives, represented for example through the US-Array effort that should lead to an order-of-magnitude improvement in our ability to resolve Earth structure seismically below North America. In fact, new observational constraints on deep Earth structure illustrate the growing importance of of improving our data assimilation skills in deep Earth models. We will explore data assimilation through high resolution global adjoint models of mantle circulation and conclude that it is feasible to reconstruct mantle flow back in time for at least the past 100 myrs.

Mantle Circulation Models with variational data assimilation: Inferring past mantle flow and structure from plate motion histories and seismic tomography

NASA Astrophysics Data System (ADS)

Bunge, H.; Hagelberg, C.; Travis, B.

2002-12-01

EarthScope will deliver data on structure and dynamics of continental North America and the underlying mantle on an unprecedented scale. Indeed, the scope of EarthScope makes its mission comparable to the large remote sensing efforts that are transforming the oceanographic and atmospheric sciences today. Arguably the main impact of new solid Earth observing systems is to transform our use of geodynamic models increasingly from conditions that are data poor to an environment that is data rich. Oceanographers and meteorologists already have made substantial progress in adapting to this environment, by developing new approaches of interpreting oceanographic and atmospheric data objectively through data assimilation methods in their models. However, a similarly rigorous theoretical framework for merging EarthScope derived solid Earth data with geodynamic models has yet to be devised. Here we explore the feasibility of data assimilation in mantle convection studies in an attempt to fit global geodynamic model calculations explicitly to tomographic and tectonic constraints. This is an inverse problem not quite unlike the inverse problem of finding optimal seismic velocity structures faced by seismologists. We derive the generalized inverse of mantle convection from a variational approach and present the adjoint equations of mantle flow. The substantial computational burden associated with solutions to the generalized inverse problem of mantle convection is made feasible using a highly efficient finite element approach based on the 3-D spherical fully parallelized mantle dynamics code TERRA, implemented on a cost-effective topical PC-cluster (geowulf) dedicated specifically to large-scale geophysical simulations. This dedicated geophysical modeling computer allows us to investigate global inverse convection problems having a spatial discretization of less than 50 km throughout the mantle. We present a synthetic high-resolution modeling experiment to demonstrate that mid

Deep mantle structure as a reference frame for movements in and on the Earth.

PubMed

Torsvik, Trond H; van der Voo, Rob; Doubrovine, Pavel V; Burke, Kevin; Steinberger, Bernhard; Ashwal, Lewis D; Trønnes, Reidar G; Webb, Susan J; Bull, Abigail L

2014-06-17

Earth's residual geoid is dominated by a degree-2 mode, with elevated regions above large low shear-wave velocity provinces on the core-mantle boundary beneath Africa and the Pacific. The edges of these deep mantle bodies, when projected radially to the Earth's surface, correlate with the reconstructed positions of large igneous provinces and kimberlites since Pangea formed about 320 million years ago. Using this surface-to-core-mantle boundary correlation to locate continents in longitude and a novel iterative approach for defining a paleomagnetic reference frame corrected for true polar wander, we have developed a model for absolute plate motion back to earliest Paleozoic time (540 Ma). For the Paleozoic, we have identified six phases of slow, oscillatory true polar wander during which the Earth's axis of minimum moment of inertia was similar to that of Mesozoic times. The rates of Paleozoic true polar wander (<1°/My) are compatible with those in the Mesozoic, but absolute plate velocities are, on average, twice as high. Our reconstructions generate geologically plausible scenarios, with large igneous provinces and kimberlites sourced from the margins of the large low shear-wave velocity provinces, as in Mesozoic and Cenozoic times. This absolute kinematic model suggests that a degree-2 convection mode within the Earth's mantle may have operated throughout the entire Phanerozoic.

Anisotropic rheology of a polycrystalline aggregate and convection in planetary mantles

NASA Astrophysics Data System (ADS)

Pouilloux, L. S.; Labrosse, S.; Kaminski, E.

2011-12-01

Observations of seismic anisotropy in the Earth mantle is often related to the crystal preferred orientation of polycrystalline aggregates. In this case, the physical properties depends on the direction and require the use of tensors to be fully described. In particular, the viscosity must be defined as a fourth order tensor whereas the thermal conductivity is a 2nd order tensor. However, the dynamical implications of such physical properties have received little attention until now. In this work, we present the mathematical formulation for an anisotropic medium and the relationship with dislocation creep deformation. We explore extensively the problem of the onset of Rayleigh-Bénard convection with such anisotropic properties. We finally presents some numerical results on the time-dependent problem using an orthotropic law for an ice polycrystal. Geophysical implications of this work related to the dynamics of planetary mantles are discussed, especially the potential of anisotropic rheology to localize deformation.

Self-Organized Mantle Layering After the Magma-Ocean Period

NASA Astrophysics Data System (ADS)

Hansen, U.; Dude, S.

2017-12-01

The thermal history of the Earth, it's chemical differentiation and also the reaction of the interior with the atmosphere is largely determined by convective processes within the Earth's mantle. A simple physical model, resembling the situation, shortly after core formation, consists of a compositionally stable stratified mantle, as resulting from fractional crystallization of the magma ocean. The early mantle is subject to heating from below by the Earth's core and cooling from the top through the atmosphere. Additionally internal heat sources will serve to power the mantle dynamics. Under such circumstances double diffusive convection will eventually lead to self -organized layer formation, even without the preexisting jumps is material properties. We have conducted 2D and 3D numerical experiments in Cartesian and spherical geometry, taking into account mantle realistic values, especially a strong temperature dependent viscosity and a pressure dependent thermal expansivity . The experiments show that in a wide parameter range. distinct convective layers evolve in this scenario. The layering strongly controls the heat loss from the core and decouples the dynamics in the lower mantle from the upper part. With time, individual layers grow on the expense of others and merging of layers does occur. We observe several events of intermittent breakdown of individual layers. Altogether an evolution emerges, characterized by continuous but also spontaneous changes in the mantle structure, ranging from multiple to single layer flow. Such an evolutionary path of mantle convection allows to interpret phenomena ranging from stagnation of slabs at various depth to variations in the chemical signature of mantle upwellings in a new framework.

Broad plumes rooted at the base of the Earth's mantle beneath major hotspots.

PubMed

French, Scott W; Romanowicz, Barbara

2015-09-03

Plumes of hot upwelling rock rooted in the deep mantle have been proposed as a possible origin of hotspot volcanoes, but this idea is the subject of vigorous debate. On the basis of geodynamic computations, plumes of purely thermal origin should comprise thin tails, only several hundred kilometres wide, and be difficult to detect using standard seismic tomography techniques. Here we describe the use of a whole-mantle seismic imaging technique--combining accurate wavefield computations with information contained in whole seismic waveforms--that reveals the presence of broad (not thin), quasi-vertical conduits beneath many prominent hotspots. These conduits extend from the core-mantle boundary to about 1,000 kilometres below Earth's surface, where some are deflected horizontally, as though entrained into more vigorous upper-mantle circulation. At the base of the mantle, these conduits are rooted in patches of greatly reduced shear velocity that, in the case of Hawaii, Iceland and Samoa, correspond to the locations of known large ultralow-velocity zones. This correspondence clearly establishes a continuous connection between such zones and mantle plumes. We also show that the imaged conduits are robustly broader than classical thermal plume tails, suggesting that they are long-lived, and may have a thermochemical origin. Their vertical orientation suggests very sluggish background circulation below depths of 1,000 kilometres. Our results should provide constraints on studies of viscosity layering of Earth's mantle and guide further research into thermochemical convection.

Structure and dynamics of Earth's lower mantle.

PubMed

Garnero, Edward J; McNamara, Allen K

2008-05-02

Processes within the lowest several hundred kilometers of Earth's rocky mantle play a critical role in the evolution of the planet. Understanding Earth's lower mantle requires putting recent seismic and mineral physics discoveries into a self-consistent, geodynamically feasible context. Two nearly antipodal large low-shear-velocity provinces in the deep mantle likely represent chemically distinct and denser material. High-resolution seismological studies have revealed laterally varying seismic velocity discontinuities in the deepest few hundred kilometers, consistent with a phase transition from perovskite to post-perovskite. In the deepest tens of kilometers of the mantle, isolated pockets of ultralow seismic velocities may denote Earth's deepest magma chamber.

Onset of solid state mantle convection and mixing during magma ocean solidification

NASA Astrophysics Data System (ADS)

Maurice, Maxime; Tosi, Nicola; Samuel, Henri; Plesa, Ana-Catalina; Hüttig, Christian; Breuer, Doris

2017-04-01

The fractional crystallization of a magma ocean can cause the formation of a compositional layering that can play a fundamental role for the subsequent long-term dynamics of the interior, for the evolution of geochemical reservoirs, and for surface tectonics. In order to assess to what extent primordial compositional heterogeneities generated by magma ocean solidification can be preserved, we investigate the solidification of a whole-mantle Martian magma ocean, and in particular the conditions that allow solid state convection to start mixing the mantle before solidification is completed. To this end, we performed 2-D numerical simulations in a cylindrical geometry. We treat the liquid magma ocean in a parametrized way while we self-consistently solve the conservation equations of thermochemical convection in the growing solid cumulates accounting for pressure-, temperature- and, where it applies, melt-dependent viscosity as well as parametrized yield stress to account for plastic yielding. By testing the effects of different cooling rates and convective vigor, we show that for a lifetime of the liquid magma ocean of 1 Myr or longer, the onset of solid state convection prior to complete mantle crystallization is likely and that a significant part of the compositional heterogeneities generated by fractionation can be erased by efficient mantle mixing.

Discontinuous Galerkin (DG) Method for solving time dependent convection-diffusion type temperature equation : Demonstration and Comparison with Other Methods in the Mantle Convection Code ASPECT

NASA Astrophysics Data System (ADS)

He, Y.; Puckett, E. G.; Billen, M. I.; Kellogg, L. H.

2016-12-01

For a convection-dominated system, like convection in the Earth's mantle, accurate modeling of the temperature field in terms of the interaction between convective and diffusive processes is one of the most common numerical challenges. In the geodynamics community using Finite Element Method (FEM) with artificial entropy viscosity is a popular approach to resolve this difficulty, but introduce numerical diffusion. The extra artificial viscosity added into the temperature system will not only oversmooth the temperature field where the convective process dominates, but also change the physical properties by increasing the local material conductivity, which will eventually change the local conservation of energy. Accurate modeling of temperature is especially important in the mantle, where material properties are strongly dependent on temperature. In subduction zones, for example, the rheology of the cold sinking slab depends nonlinearly on the temperature, and physical processes such as slab detachment, rollback, and melting all are sensitively dependent on temperature and rheology. Therefore methods that overly smooth the temperature may inaccurately represent the physical processes governing subduction, lithospheric instabilities, plume generation and other aspects of mantle convection. Here we present a method for modeling the temperature field in mantle dynamics simulations using a new solver implemented in the ASPECT software. The new solver for the temperature equation uses a Discontinuous Galerkin (DG) approach, which combines features of both finite element and finite volume methods, and is particularly suitable for problems satisfying the conservation law, and the solution has a large variation locally. Furthermore, we have applied a post-processing technique to insure that the solution satisfies a local discrete maximum principle in order to eliminate the overshoots and undershoots in the temperature locally. To demonstrate the capabilities of this new

Mantle Plumes and Geologically Recent Volcanism on Mars

NASA Astrophysics Data System (ADS)

Kiefer, W. S.

2013-12-01

Despite its small size, Mars has remained volcanically active until the geologically recent past. Crater retention ages on the volcanos Arsia Mon, Olympus Mons, and Pavonis Mons indicate significant volcanic activity in the last 100-200 million years. The radiometric ages of many shergottites, a type of igneous martian meteorite, indicate igneous activity at about 180 million years ago. These ages correspond to the most recent 2-4% of the age of the Solar System. The most likely explanation for this young martian volcanism is adiabatic decompression melting in upwelling mantle plumes. Multiple plumes may be active at any time, with each of the major volcanos in the Tharsis region being formed by a separate plume. Like at least some terrestrial mantle plumes, mantle plumes on Mars likely form via an instability of the thermal boundary layer at the base of the mantle. Because Mars operates in the stagnant lid convection regime, the temperature difference between mantle and core is lower than on Earth. This reduces the temperature contrast between mantle and core, resulting in mantle plumes on Mars that are about 100 K hotter than the average mantle. The chemical composition of the martian meteorites indicates that the martian mantle is enriched in both iron and sodium relative to Earth's mantle. This lowers the dry solidus on early Mars by 30-40 K relative to Earth. Migration of sodium to the crust over time decreases this difference in solidus temperature to about 15 K at present, but that is sufficient to increase the current plume magma production rate by a factor of about 2. Hydrous phases in the martian meteorites indicate the presence of a few hundred ppm water in the mantle source region, roughly the same as Earth. Finite element simulations of martian plumes using temperature-dependent viscosity and realistic Rayleigh numbers can reproduce the geologically recent magma production rate that is inferred from geologic mapping and the melt fraction inferred from

Osmium isotope variations in the Pacific mantle: implications for the distribution of heterogeneity in the convecting mantle

NASA Astrophysics Data System (ADS)

Ishikawa, A.; Senda, R.; Suzuki, K.; Tani, K.; Ishii, T.

2015-12-01

Recent accumulation of Os isotope data obtained either from abyssal peridotites or from ocean island peridotite xenoliths has clearly demonstrated that the modern convecting mantle is substantially heterogeneous in Os-isotope composition. Unlike other radiogenic isotope heterogeneities observed in oceanic basalts, largely controlled by incorporation of recycled crustal materials, it seems likely that the observed range of Os-isotope compositions in oceanic peridotites directly reflect varying degrees of ancient melt extraction from peridotitic mantle. Hence, global variations of Os-isotope compositions in oceanic peridotites may provide an important piece of information in unraveling the geochemical and geodynamic evolution of the convecting mantle. Here we present the Os-isotope variations in peridotite-serpentinite recovered from the Pacific area because the number of data available is yet scarce when compared with data from other oceans (Atlantic, Arctic and Indian Ocean). Our primary purpose is to test whether mantle domains underlying four major oceans are distinct in terms of Os isotope variations, reflecting the pattern of mantle convection or mixing efficiency. We examined 187Os/188Os ratios and highly siderophile element concentrations in serpentinized harzburgite recovered from Hess Deep in the East Pacific Rise, a mantle section in the Taitao ophiolite, Chile (Schulte et al., 2009), serpentinized harzburgite bodies in the Izu-Ogasawara and Tonga forearc (Parkinson et al., 1998), peridotite xenoliths from the Pali-Kaau vent in O'ahu island, Hawaii (Bizimis et al., 2007), and low-temperature type peridotite xenoliths from Malaita, Solomon Islands (Ishikawa et al., 2011). The results demonstrate that samples from each area display very similar Os-isotope variations with a pronounced peak in 187Os/188Os = 0.125-0.128. Moreover, the relatively larger datasets obtained from Hess Deep, Taitao and Malaita clearly exhibit the presence of secondary peak in 187Os

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.

Mantle convection patterns reveal the enigma of the Red Sea rifting

NASA Astrophysics Data System (ADS)

Petrunin, Alexey; Kaban, Mikhail; El Khrepy, Sami; Al-Arifi, Nassir

2017-04-01

Initiation and further development of the Red Sea rift (RSR) is usually associated with the Afar plume at the Oligocene-Miocene separating the Arabian plate from the rest of the continent. Usually, the RSR is divided into three parts with different geological, tectonic and geophysical characteristics, but the nature of this partitioning is still debatable. To understand origin and driving forces responsible for the tectonic partitioning of the RSR, we have developed a global mantle convection model based on the refined density model and viscosity distribution derived from tectonic, rheological and seismic data. The global density model of the upper mantle is refined for the Middle East based on the high-resolution 3D model (Kaban et al., 2016). This model based on a joint inversion of the residual gravity and residual topography provides much better constraints on the 3D density structure compared to the global model based on seismic tomography. The refined density model and the viscosity distribution based on a homologous temperature approach provide an initial setup for further numerical calculations. The present-day snapshot of the mantle convection is calculated by using the code ProSpher 3D that allows for strong lateral variations of viscosity (Petrunin et al., 2013). The setup includes weak plate boundaries, while the measured GPS velocities are used to constrain the solution. The resulting mantle flow patterns show clear distinctions among the mantle flow patterns below the three parts of the RSR. According to the modeling results, tectonics of the southern part of the Red Sea is mainly determined by the Afar plume and the Ethiopian rift opening. It is characterized by a divergent mantle flow, which is connected to the East African Rift activity. The rising mantle flow is traced down to the transition zone and continues in the lower mantle for a few thousand kilometers south-west of Afar. The hot mantle anomaly below the central part of the RSR can be

Reduced Lattice Thermal Conductivity of Fe-bearing Bridgmanite in Earth's Deep Mantle

NASA Astrophysics Data System (ADS)

Hsieh, W. P.; Deschamps, F.; Okuchi, T.; Lin, J. F.

2017-12-01

Complex seismic and thermo-chemical features have been revealed in Earth's lowermost mantle. Particularly, possible iron enrichments in the large low shear-wave velocity provinces (LLSVPs) could influence thermal transport properties of the constituting minerals in this region, which, in turn, may alter the lower mantle dynamics and heat flux across core-mantle boundary (CMB). Thermal conductivity of bridgmanite is expected to partially control the thermal evolution and dynamics of Earth's lower mantle. Importantly, the pressure-induced lattice distortion in bridgmanite could affect its lattice thermal conductivity, but this effect remains largely unknown. Here we report our measurements of the lattice thermal conductivity of Fe-bearing and (Fe,Al)-bearing bridgmanites to 120 GPa using optical pump-probe spectroscopy. The thermal conductivity of Fe-bearing bridgmanite increases monotonically with pressure, but drops significantly around 45 GPa presumably due to pressure-induced lattice distortion on iron sites. Our findings indicate that lattice thermal conductivity at lowermost mantle conditions is twice smaller than previously thought. The decrease in the thermal conductivity of bridgmanite in mid-lower mantle and below would promote mantle flow against a potential viscosity barrier, facilitating slabs crossing over the 1000-km depth. Modeling of our results applied to the LLSVPs shows that variations in iron and bridgmanite fractions induce a significant thermal conductivity decrease, which would enhance internal convective flow. Our CMB heat flux modeling indicates that, while heat flux variations are dominated by thermal effects, variations in thermal conductivity also play a significant role. The CMB heat flux map we obtained is substantially different from those assumed so far, which may influence our understanding of the geodynamo.

Subduction and volatile recycling in Earth's mantle

NASA Technical Reports Server (NTRS)

King, S. D.; Ita, J. J.; Staudigel, H.

1994-01-01

The subduction of water and other volatiles into the mantle from oceanic sediments and altered oceanic crust is the major source of volatile recycling in the mantle. Until now, the geotherms that have been used to estimate the amount of volatiles that are recycled at subduction zones have been produced using the hypothesis that the slab is rigid and undergoes no internal deformation. On the other hand, most fluid dynamical mantle flow calculations assume that the slab has no greater strength than the surrounding mantle. Both of these views are inconsistent with laboratory work on the deformation of mantle minerals at high pressures. We consider the effects of the strength of the slab using two-dimensional calculations of a slab-like thermal downwelling with an endothermic phase change. Because the rheology and composition of subducting slabs are uncertain, we consider a range of Clapeyron slopes which bound current laboratory estimates of the spinel to perovskite plus magnesiowustite phase transition and simple temperature-dependent rheologies based on an Arrhenius law diffusion mechanism. In uniform viscosity convection models, subducted material piles up above the phase change until the pile becomes gravitationally unstable and sinks into the lower mantle (the avalanche). Strong slabs moderate the 'catastrophic' effects of the instabilities seen in many constant-viscosity convection calculations; however, even in the strongest slabs we consider, there is some retardation of the slab descent due to the presence of the phase change.

Benchmarking FEniCS for mantle convection simulations

NASA Astrophysics Data System (ADS)

Vynnytska, L.; Rognes, M. E.; Clark, S. R.

2013-01-01

This paper evaluates the usability of the FEniCS Project for mantle convection simulations by numerical comparison to three established benchmarks. The benchmark problems all concern convection processes in an incompressible fluid induced by temperature or composition variations, and cover three cases: (i) steady-state convection with depth- and temperature-dependent viscosity, (ii) time-dependent convection with constant viscosity and internal heating, and (iii) a Rayleigh-Taylor instability. These problems are modeled by the Stokes equations for the fluid and advection-diffusion equations for the temperature and composition. The FEniCS Project provides a novel platform for the automated solution of differential equations by finite element methods. In particular, it offers a significant flexibility with regard to modeling and numerical discretization choices; we have here used a discontinuous Galerkin method for the numerical solution of the advection-diffusion equations. Our numerical results are in agreement with the benchmarks, and demonstrate the applicability of both the discontinuous Galerkin method and FEniCS for such applications.

Stability and growth of continental shields in mantle convection models including recurrent melt production

NASA Astrophysics Data System (ADS)

de Smet, J. H.; van den Berg, A. P.; Vlaar, N. J.

1998-10-01

The long-term growth and stability of compositionally layered continental upper mantle has been investigated by numerical modelling. We present the first numerical model of a convecting mantle including differentiation through partial melting resulting in a stable compositionally layered continental upper mantle structure. This structure includes a continental root extending to a depth of about 200 km. The model covers the upper mantle including the crust and incorporates physical features important for the study of the continental upper mantle during secular cooling of the Earth since the Archaean. Among these features are: a partial melt generation mechanism allowing consistent recurrent melting, time-dependent non-uniform radiogenic heat production, and a temperature- and pressure-dependent rheology. The numerical results reveal a long-term growth mechanism of the continental compositional root. This mechanism operates through episodical injection of small diapiric upwellings from the deep layer of undepleted mantle into the continental root which consists of compositionally distinct depleted mantle material. Our modelling results show the layered continental structure to remain stable during at least 1.5 Ga. After this period mantle differentiation through partial melting ceases due to the prolonged secular cooling and small-scale instabilities set in through continental delamination. This stable period of 1.5 Ga is related to a number of limitations in our model. By improving on these limitations in the future this stable period will be extended to more realistic values.

Fluctuations in seafloor spreading predicted by tectonic reconstructions and mantle convection models

NASA Astrophysics Data System (ADS)

Coltice, Nicolas; Seton, Maria; Rolf, Tobias; Müller, R. Dietmar; Tackley, Paul J.

2013-04-01

The theory of plate tectonics theory has enabled possible the reconstruction of the ancient seafloor and paleogeography. Over 50 years of data collection and kinematic reconstruction efforts, plate models have improved significantly (Seton et al., 2012) although reconstructions of ancient seafloor are naturally limited by the limited preservation of of very old seafloor. It is challenging to reconstruct ancient ocean basins and associated plate boundaries for times earlier than 200 Ma, since seafloor of this age is not preserved. This means we can merely reconstruct only 5% of the history of the planet in this fashion. However, geodynamic models can now help evaluate how seafloor spreading may evolve over longer time periods, since recent developments of numerical models of mantle convection with pseudo-plasticity can generate long-term solutions that simulate a form of seafloor spreading (Moresi and Solomatov, 1998; Tackley, 2000a; Tackley, 2000b). The introduction of models of continental lithosphere further improves the quality of the predictions: the computed distribution of seafloor ages reproduces the consumption of young seafloor as observed on the present-day Earth (Coltice et al., 2012). The time-dependence of the production of new seafloor has long been debated and there is no consensus on how much it has varied in the past 150My, and how it could have fluctuated over longer time-scales. Using plate reconstructions, Parsons (1982) and Rowley (2002) proposed the area vs. age distribution of the seafloor could have experienced limited fluctuations in the past 150My while others suggest stronger variations would fit the observations equally well (Seton et al., 2009. Here we propose to investigate the global dynamics of seafloor spreading using state-of-the-art plate reconstructions and geodynamic models. We focus on the evolution of the distribution of seafloor ages because fundamental geophysical observations like mantle heat flow or sea level provide

Tidal tomography constrains Earth's deep-mantle buoyancy.

PubMed

Lau, Harriet C P; Mitrovica, Jerry X; Davis, James L; Tromp, Jeroen; Yang, Hsin-Ying; Al-Attar, David

2017-11-15

Earth's body tide-also known as the solid Earth tide, the displacement of the solid Earth's surface caused by gravitational forces from the Moon and the Sun-is sensitive to the density of the two Large Low Shear Velocity Provinces (LLSVPs) beneath Africa and the Pacific. These massive regions extend approximately 1,000 kilometres upward from the base of the mantle and their buoyancy remains actively debated within the geophysical community. Here we use tidal tomography to constrain Earth's deep-mantle buoyancy derived from Global Positioning System (GPS)-based measurements of semi-diurnal body tide deformation. Using a probabilistic approach, we show that across the bottom two-thirds of the two LLSVPs the mean density is about 0.5 per cent higher than the average mantle density across this depth range (that is, its mean buoyancy is minus 0.5 per cent), although this anomaly may be concentrated towards the very base of the mantle. We conclude that the buoyancy of these structures is dominated by the enrichment of high-density chemical components, probably related to subducted oceanic plates or primordial material associated with Earth's formation. Because the dynamics of the mantle is driven by density variations, our result has important dynamical implications for the stability of the LLSVPs and the long-term evolution of the Earth system.

Source of Volatiles in Earth's Deep Mantle from Neon Isotope Systematics in the South Atlantic

NASA Astrophysics Data System (ADS)

Williams, C. D.; Mukhopadhyay, S.

2016-12-01

The noble gases play an important role in understanding Earth's accretion and subsequent evolution. Neon isotopes in particular have the potential to distinguish between distinct sources of Earth's volatiles e.g., acquisition of nebular gas, solar wind implanted materials or chondritic meteorites and their components. The neon isotopic composition of the deep mantle remains subject to debate with the majority of mantle-derived basalts displaying maximum 20Ne/22Ne ratios less than 12.5, similar to values determined for the convective mantle (20Ne/22Ne = 12.49 +/- 0.04; [1]). These values are also much lower than those of solar wind (20Ne/22Ne = 13.8; [2,3]) and estimates of the nebular gas (20Ne/22Ne = 13.4; [4]) but comparable to solar wind implanted meteoritic materials (20Ne/22Ne = 12.5-12.7; [5]). Here we determine the neon isotopic composition of mantle-derived materials from the south Atlantic. These samples display strong linear correlations in 20Ne/22Ne-21Ne/22Ne space with maximum 20Ne/22Ne ratios that are resolvable from and higher than materials derived from the convecting mantle as well as models of solar wind implantation. These results supplement a growing database of mantle materials characterized by 20Ne/22Ne ratios greater than 12.5, challenging the notion that the entire mantle acquired volatiles from solar wind implanted meteoritic materials. In this presentation we will explore alternative origins for these volatiles and provide testable predictions for each scenario. [1] G. Holland, C.J. Ballentine.. Nature 441 (2006), 186-191. [2] A. Gimberg et al. GCA 72 (2008), 626-645. [3] V.S. Heber et al. GCA 73 (2009), 7414-7432. [4] V. S. Heber et al. ApJ 759 (2012), 121. [5] R. Wieler in: D. Porcelli, C.J. Ballentine, R. Wieler (Eds.), Reviews in Mineralogy and Geochemistry 47 (2002), 21-70.

Three-dimensional instabilities of mantle convection with multiple phase transitions

NASA Technical Reports Server (NTRS)

Honda, S.; Yuen, D. A.; Balachandar, S.; Reuteler, D.

1993-01-01

The effects of multiple phase transitions on mantle convection are investigated by numerical simulations that are based on three-dimensional models. These simulations show that cold sheets of mantle material collide at junctions, merge, and form a strong downflow that is stopped temporarily by the transition zone. The accumulated cold material gives rise to a strong gravitational instability that causes the cold mass to sink rapidly into the lower mantle. This process promotes a massive exchange between the lower and upper mantles and triggers a global instability in the adjacent plume system. This mechanism may be cyclic in nature and may be linked to the generation of superplumes.

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.

Water and hydrogen are immiscible in Earth's mantle.

PubMed

Bali, Enikő; Audétat, Andreas; Keppler, Hans

2013-03-14

In the deep, chemically reducing parts of Earth's mantle, hydrous fluids contain significant amounts of molecular hydrogen (H2). Thermodynamic models of fluids in Earth's mantle so far have always assumed that molecular hydrogen and water are completely miscible. Here we show experimental evidence that water and hydrogen can coexist as two separate, immiscible phases. Immiscibility between water and hydrogen may be the cause of the formation of enigmatic, ultra-reducing domains in the mantle that contain moissanite (SiC) and other phases indicative of extremely reducing conditions. Moreover, the immiscibility between water and hydrogen may provide a mechanism for the rapid oxidation of Earth's upper mantle immediately following core formation.

Decreasing µ142Nd Variation in the Archean Convecting Mantle from 4.0 to 2.5 Ga: Heterogeneous Domain Mixing or Crustal Recycling?

NASA Astrophysics Data System (ADS)

Brandon, A. D.; Debaille, V.

2014-12-01

The 146Sm-142Nd (t1/2=68 Ma) chronometer can be used to examine silicate differentiation in the first 400 Ma of Earth history. Early fractionation between Sm and Nd is recorded in cratonic Archean rocks in their 142Nd/144Nd ratios that that deviate up to ±20 ppm, or μ142Nd - ppm deviation relative to the present-day convecting mantle at 0. These values likely record early extraction of incompatible trace element (ITE) enriched material with -μ142Nd, either as crust or late stage residual melt from a magma ocean, and resulting in a complimentary ITE depleted residual mantle with +μ142Nd. If this early-formed ITE-enriched material was re-incorporated rapidly back into the convecting mantle, both ITE-enriched and ITE-depleted mantle domains would have been established in the Hadean. Alternatively, if it was early-formed crust that remained stable it could have slowly eroded and progressively remixed into the convecting mantle as subducted sediment during the Archean. Each of these scenarios could potentially explain the decrease in the maximum variation in µ142Nd from ±20 at 4.0 Ga to 0 at 2.5 Ga [1,2,3]. In the scenario where these variations reflect mixing of mantle domains, this implies long mantle mixing times of greater than 1 Ga in the Archean in order to preserve the early-formed heterogeneities. This can be achieved in a stagnant lid tectonic regime in the Archean with sporadic and short subduction cycles [2]. This scenario would also indicate that mixing times in the convecting mantle were much slower than the previously proposed 100 Ma in the Hadean and Archean. In the alternative scenario, sediment with -µ142Nd was progressively mixed into the mantle via subduction in the Archean [3]. This scenario doesn't require slow mantle mixing times or a stagnant-lid regime. It requires crustal resident times of up to 750 Ma to maintain a steady supply of ancient sediment recycling over the Archean. Each of these scenarios evoke very contrasting conditions for

Effects from equation of state and rheology in dissipative heating in compressible mantle convection

NASA Technical Reports Server (NTRS)

Yuen, David A.; Quareni, Francesca; Hong, H.-J.

1987-01-01

The effects of compressibility on mantle convection are considered, incorporating the effects of equations of state and rheology in the dissipative heating term of the energy equation. The ways in which compression may raise the interior mantle temperature are explicitly demonstrated, and it is shown how this effect can be used to constrain some of the intrinsic parameters associated with the equation of state in the mantle. It is concluded that the coupling between variable viscosity and equation of state in dissipative heating is potentially an important mechanism in mantle convection. These findings emphasize that rheology, equation of state, and radiogenic heating are all linked to each other by nonlinear thermomechanical couplings.

Earth's Fiercely Cooling Core - 24 TW

NASA Astrophysics Data System (ADS)

Morgan, Jason P.; Vannucchi, Paola

2014-05-01

Earth's mantle and core are convecting planetary heat engines. The mantle convects to lose heat from slow cooling, internal radioactivity, and core heatflow across its base. Its convection generates plate tectonics, volcanism, and the loss of ~35 TW of mantle heat through Earth's surface. The core convects to lose heat from slow cooling, small amounts of internal radioactivity, and the freezing-induced growth of a compositionally denser inner core. Core convection produces the geodynamo generating Earth's geomagnetic field. The geodynamo was thought to be powered by ~4 TW of heatloss across the core-mantle boundary, a rate sustainable (cf. Gubbins et al., 2003; Nimmo, 2007) by freezing a compositionally denser inner core over the ~3 Ga that Earth is known to have had a strong geomagnetic field (cf. Tarduno, 2007). However, recent determinations of the outer core's thermal conductivity(Pozzo et al., 2012; Gomi et al., 2013) indicate that >15 TW of power should conduct down its adiabat. Conducted power is unavailable to drive thermal convection, implying that the geodynamo needs a long-lived >17 TW power source. Core cooling was thought too weak for this, based on estimates for the Clapeyron Slope for high-pressure freezing of an idealized pure-iron core. Here we show that the ~500-1000 kg/m3 seismically-inferred jump in density between the liquid outer core and solid inner core allows us to directly infer the core-freezing Clapeyron Slope for the outer core's actual composition which contains ~8±2% lighter elements (S,Si,O,Al, H,…) mixed into a Fe-Ni alloy. A PREM-like 600 kg/m3 - based Clapeyron Slope implies there has been ~774K of core cooling during the freezing and growth of the inner core, releasing ~24 TW of power during the past ~3 Ga. If so, core cooling can easily power Earth's long-lived geodynamo. Another major implication of ~24 TW heatflow across the core-mantle boundary is that the present-day mantle is strongly 'bottom-heated', and diapiric mantle

A sequential data assimilation approach for the joint reconstruction of mantle convection and surface tectonics

NASA Astrophysics Data System (ADS)

Bocher, M.; Coltice, N.; Fournier, A.; Tackley, P. J.

2016-01-01

With the progress of mantle convection modelling over the last decade, it now becomes possible to solve for the dynamics of the interior flow and the surface tectonics to first order. We show here that tectonic data (like surface kinematics and seafloor age distribution) and mantle convection models with plate-like behaviour can in principle be combined to reconstruct mantle convection. We present a sequential data assimilation method, based on suboptimal schemes derived from the Kalman filter, where surface velocities and seafloor age maps are not used as boundary conditions for the flow, but as data to assimilate. Two stages (a forecast followed by an analysis) are repeated sequentially to take into account data observed at different times. Whenever observations are available, an analysis infers the most probable state of the mantle at this time, considering a prior guess (supplied by the forecast) and the new observations at hand, using the classical best linear unbiased estimate. Between two observation times, the evolution of the mantle is governed by the forward model of mantle convection. This method is applied to synthetic 2-D spherical annulus mantle cases to evaluate its efficiency. We compare the reference evolutions to the estimations obtained by data assimilation. Two parameters control the behaviour of the scheme: the time between two analyses, and the amplitude of noise in the synthetic observations. Our technique proves to be efficient in retrieving temperature field evolutions provided the time between two analyses is ≲10 Myr. If the amplitude of the a priori error on the observations is large (30 per cent), our method provides a better estimate of surface tectonics than the observations, taking advantage of the information within the physics of convection.

Convective instability within the Tibetan Lithospheric Mantle (Invited)

NASA Astrophysics Data System (ADS)

Houseman, G. A.; Molnar, P. H.; Evans, L.; England, P. C.

2013-12-01

Studies of seismic surface waves in Asia show that shear-wave speeds at depths of ~120-250km beneath the Tibetan Plateau are higher than is generally observed for continents, other than beneath Archaean cratons. The high-speed layer has been interpreted as continental lithosphere that was thickened during the convergence between India and Asia. This interpretation contradicts conceptual models in which gravitational instabilities remove a significant fraction of the mantle lithosphere beneath Tibet during that convergence. In contrast, the suggestion of relatively recent (post-early-Miocene) surface uplift of the Plateau, inferred from the onset of normal faulting across the plateau, synchronous increased rates of compressional deformation in the surroundings of the the plateau, and widespread volcanism in the northern part of the plateau, implies action of a mechanism that increased the gravitational potential energy of, and temperatures within, the Tibetan lithosphere in a way that would not occur if the mantle lithosphere had simply thickened continually throughout the India-Asia convergence. A resolution to this paradox is suggested by the observation that, while shear-wave speeds are indeed high at depths of 120-250 km beneath the Tibetan plateau, they are anomalously low at shallower depths, implying a temperature inversion that is hard to reconcile with uninterrupted lithospheric thickening. We suggest that the ensemble of observations may be explained by the convective overturn of a lithospheric root that is depleted in iron such that it remains buoyant with respect to normal upper mantle. The increased rate of strain within the Tibetan lithosphere once convergence began reduced its effective viscosity, and continuing convergence thickened the lithospheric root. These conditions led to convective overturn, similar to the original conceptual models, with the difference that the overturn was confined within the root, which remains buoyant with respect to the

A Comparison of Methods for Modeling Geochemical Variability in the Earth's Mantle

NASA Astrophysics Data System (ADS)

Kellogg, J. B.; Tackley, P. J.

2004-12-01

Numerial models of isotopic and chemical heterogeneity of the Earth's mantle fall into three categories, in decreasing order of computational demand. First, several authors have used chemical tracers within a full thermo-chemical convection calculation (e.g., Christensen and Hofmann, 1994, van Keken and Ballentine, 1999; Xie and Tackley, 2004). Second, Kellogg et al. (2002) proposed an extension of the traditional geochemical box model calculations in which numerous subreservoirs were tracked within the bulk depleted mantle reservoir. Third, Allègre and Lewin (1995) described a framework in which the variance in chemical and isotopic ratios were treated as quantities intrinsic to the bulk reservoirs, complete with sources and sinks. Results from these three methods vary, particularly with respect to conclusions drawn about the meaning of the Pb-Pb pseudo-isochron. We revisit these methods in an attempt to arrive at a common understanding. By considering all three we better identify the strengths and weaknesses of each approach and allow each to inform the other. Finally, we present results from a new hybrid model that combines the complexity and regional-scale variability of the thermochemical convection models with the short length-scale sensitivity of the Kellogg et al. approach.

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.

The generation of plate tectonics from mantle convection

NASA Astrophysics Data System (ADS)

Bercovici, David

2003-01-01

In the last decade, significant progress has been made toward understanding how plate tectonics is generated from mantle dynamics. A primary goal of plate-generation studies has been the development of models that allow the top cold thermal boundary layer of mantle convection, i.e. the lithosphere, to develop broad and strong plate-like segments separated by narrow, weak and rapidly deforming boundaries; ideally, such models also permit significant strike-slip (toroidal) motion, passive ridges (i.e. pulled rather than pried apart), and self-consistent initiation of subduction. A major outcome of work so far is that nearly all aspects of plate generation require lithospheric rheologies and shear-localizing feedback mechanisms that are considerably more exotic than rheologies typically used in simple fluid-dynamical models of mantle flow. The search for plate-generating behavior has taken us through investigations of the effects of shear weakening ('stick-slip') and viscoplastic rheologies, of melting at ridges and low-viscosity asthenospheres, and of grain-size dependent rheologies and damage mechanics. Many such mechanisms, either by themselves or in combination, have led to self-consistent fluid-mechanical models of mantle flow that are remarkably plate-like, which is in itself a major accomplishment. However, many other important problems remain unsolved, such as subduction intiation and asymmetry, temporal evolution of plate geometry, rapid changes in plate motion, and the Archaean initiation of the plate-tectonic mode of convection. This paper presents a brief review of progress made in the plate-generation problem over the last decade, and discusses unresolved issues and future directions of research in this important area.

Simulation of active tectonic processes for a convecting mantle with moving continents

USGS Publications Warehouse

Trubitsyn, V.; Kaban, M.; Mooney, W.; Reigber, C.; Schwintzer, P.

2006-01-01

Numerical models are presented that simulate several active tectonic processes. These models include a continent that is thermally and mechanically coupled with viscous mantle flow. The assumption of rigid continents allows use of solid body equations to describe the continents' motion and to calculate their velocities. The starting point is a quasi-steady state model of mantle convection with temperature/ pressure-dependent viscosity. After placing a continent on top of the mantle, the convection pattern changes. The mantle flow subsequently passes through several stages, eventually resembling the mantle structure under present-day continents: (a) Extension tectonics and marginal basins form on boundary of a continent approaching to subduction zone, roll back of subduction takes place in front of moving continent; (b) The continent reaches the subduction zone, the extension regime at the continental edge is replaced by strong compression. The roll back of the subduction zone still continues after closure of the marginal basin and the continent moves towards the upwelling. As a result the ocean becomes non-symmetric and (c) The continent overrides the upwelling and subduction in its classical form stops. The third stage appears only in the upper mantle model with localized upwellings. ?? 2006 The Authors Journal compilation ?? 2006 RAS.

Sulfur in Earth's Mantle and Its Behavior During Core Formation

NASA Technical Reports Server (NTRS)

Chabot, Nancy L.; Righter,Kevin

2006-01-01

The density of Earth's outer core requires that about 5-10% of the outer core be composed of elements lighter than Fe-Ni; proposed choices for the "light element" component of Earth's core include H, C, O, Si, S, and combinations of these elements [e.g. 1]. Though samples of Earth's core are not available, mantle samples contain elemental signatures left behind from the formation of Earth's core. The abundances of siderophile (metal-loving) elements in Earth's mantle have been used to gain insight into the early accretion and differentiation history of Earth, the process by which the core and mantle formed, and the composition of the core [e.g. 2-4]. Similarly, the abundance of potential light elements in Earth's mantle could also provide constraints on Earth's evolution and core composition. The S abundance in Earth's mantle is 250 ( 50) ppm [5]. It has been suggested that 250 ppm S is too high to be due to equilibrium core formation in a high pressure, high temperature magma ocean on early Earth and that the addition of S to the mantle from the subsequent accretion of a late veneer is consequently required [6]. However, this earlier work of Li and Agee [6] did not parameterize the metalsilicate partitioning behavior of S as a function of thermodynamic variables, limiting the different pressure and temperature conditions during core formation that could be explored. Here, the question of explaining the mantle abundance of S is revisited, through parameterizing existing metal-silicate partitioning data for S and applying the parameterization to core formation in Earth.

Deformation of "stable" continental interiors by mantle convection: Implications for intraplate stress in the New Madrid Seismic Zone

NASA Astrophysics Data System (ADS)

Forte, A. M.; Moucha, R.; Simmons, N. A.; Grand, S. P.; Mitrovica, J. X.

2011-12-01

The enigmatic origin of large-magnitude earthquakes far from active plate boundaries, especially those occurring in so-called "stable" continental interiors, is a source of continuing controversy that has eluded a satisfactory explanation using past geophysical models of intraplate deformation and faulting. One outstanding case of such major intraplate earthquakes is the 1811-1812 series of events in the New Madrid Seismic Zone (NMSZ). We contend that the origin of some of these enigmatic intraplate events is due to regional variations in the pattern of tectonic stress generated by mantle convective flow acting on the overlying lithosphere and crust. Mantle convection affects the entire surface of the planet, irrespective of the current configuration of surface plate boundaries. In addition, it must be appreciated that plate tectonics is not a 2-D process, because the convective flow that drives the observed horizontal motions of the tectonic plates also drives vertical displacements of the crust across distances as great as 2 to 3 km. This dynamic topography is directly correlated with convection-driven stress field variations in the crust and lithosphere and these stresses can be locally focussed if the mantle rheology below the lithosphere is characterised by sufficiently low viscosities. We have developed global models of convection-driven mantle flow [Forte et al. 2009,2010] that are based on recent high-resolution 3-D tomography models derived from joint inversions of seismic, geodynamic and mineral physics data [Simmons et al. 2007,2008,2010]. These tomography-based mantle convection models also include a full suite of surface geodynamic (postglacial rebound and convection) constraints on the depth-dependent average viscosity of the mantle [Mitrovica & Forte 2004]. Our latest tomography-based and geodynamically-constrained convection calculations reveal that mantle flow under the central US are driven by density anomalies within the lower mantle associated

Iron-carbonate interaction at Earth's core-mantle boundary

NASA Astrophysics Data System (ADS)

Dorfman, S. M.; Badro, J.; Nabiei, F.; Prakapenka, V.; Gillet, P.

2015-12-01

Carbon storage and flux in the deep Earth are moderated by oxygen fugacity and interactions with iron-bearing phases. The amount of carbon stored in Earth's mantle versus the core depends on carbon-iron chemistry at the core-mantle boundary. Oxidized carbonates subducted from Earth's surface to the lowermost mantle may encounter reduced Fe0 metal from disproportionation of Fe2+ in lower mantle silicates or mixing with the core. To understand the fate of carbonates in the lowermost mantle, we have performed experiments on sandwiches of single-crystal (Ca0.6Mg0.4)CO3 dolomite and Fe foil in the laser-heated diamond anvil cell at lower mantle conditions of 49-110 GPa and 1800-2500 K. Syntheses were conducted with in situ synchrotron X-ray diffraction to identify phase assemblages. After quench to ambient conditions, samples were sectioned with a focused Ga+ ion beam for composition analysis with transmission electron microscopy. At the centers of the heated spots, iron melted and reacted completely with the carbonate to form magnesiowüstite, iron carbide, diamond, magnesium-rich carbonate and calcium carbonate. In samples heated at 49 and 64 GPa, the two carbonates exhibit a eutectoid texture. In the sample heated at 110 GPa, the carbonates form rounded ~150-nm-diameter grains with a higher modal proportion of interspersed diamonds. The presence of reduced iron in the deep lower mantle and core-mantle boundary region will promote the formation of diamonds in carbonate-bearing subducted slabs. The complete reaction of metallic iron to oxides and carbides in the presence of mantle carbonate supports the formation of these phases at the Earth's core-mantle boundary and in ultra-low velocity zones.

Modelling the possible interaction between edge-driven convection and the Canary Islands mantle plume

NASA Astrophysics Data System (ADS)

Negredo, A. M.; Rodríguez-González, J.; Fullea, J.; Van Hunen, J.

2017-12-01

The close location between many hotspots and the edges of cratonic lithosphere has led to the hypothesis that these hotspots could be explained by small-scale mantle convection at the edge of cratons (Edge Driven Convection, EDC). The Canary Volcanic Province hotspot represents a paradigmatic example of this situation due to its close location to the NW edge of the African Craton. Geochemical evidence, prominent low seismic velocity anomalies in the upper and lower mantle, and the rough NE-SW age-progression of volcanic centers consistently point out to a deep-seated mantle plume as the origin of the Canary Volcanic Province. It has been hypothesized that the plume material could be affected by upper mantle convection caused by the thermal contrast between thin oceanic lithosphere and thick (cold) African craton. Deflection of upwelling blobs due to convection currents would be responsible for the broader and more irregular pattern of volcanism in the Canary Province compared to the Madeira Province. In this study we design a model setup inspired on this scenario to investigate the consequences of possible interaction between ascending mantle plumes and EDC. The Finite Element code ASPECT is used to solve convection in a 2D box. The compositional field and melt fraction distribution are also computed. Free slip along all boundaries and constant temperature at top and bottom boundaries are assumed. The initial temperature distribution assumes a small long-wavelength perturbation. The viscosity structure is based on a thick cratonic lithosphere progressively varying to a thin, or initially inexistent, oceanic lithosphere. The effects of assuming different rheologies, as well as steep or gradual changes in lithospheric thickness are tested. Modelling results show that a very thin oceanic lithosphere (< 30 km) is needed to generate partial melting by EDC. In this case partial melting can occur as far as 700 km away from the edge of the craton. The size of EDC cells is

The tungsten isotopic composition of the Earth's mantle before the terminal bombardment.

PubMed

Willbold, Matthias; Elliott, Tim; Moorbath, Stephen

2011-09-07

Many precious, 'iron-loving' metals, such as gold, are surprisingly abundant in the accessible parts of the Earth, given the efficiency with which core formation should have removed them to the planet's deep interior. One explanation of their over-abundance is a 'late veneer'--a flux of meteorites added to the Earth after core formation as a 'terminal' bombardment that culminated in the cratering of the Moon. Some 3.8 billion-year-old rocks from Isua, Greenland, are derived from sources that retain an isotopic memory of events pre-dating this cataclysmic meteorite shower. These Isua samples thus provide a window on the composition of the Earth before such a late veneer and allow a direct test of its importance in modifying the composition of the planet. Using high-precision (less than 6 parts per million, 2 standard deviations) tungsten isotope analyses of these rocks, here we show that they have a isotopic tungsten ratio (182)W/(184)W that is significantly higher (about 13 parts per million) than modern terrestrial samples. This finding is in good agreement with the expected influence of a late veneer. We also show that alternative interpretations, such as partial remixing of a deep-mantle reservoir formed in the Hadean eon (more than four billion years ago) or core-mantle interaction, do not explain the W isotope data well. The decrease in mantle (182)W/(184)W occurs during the Archean eon (about four to three billion years ago), potentially on the same timescale as a notable decrease in (142)Nd/(144)Nd (refs 3 and 6). We speculate that both observations can be explained if late meteorite bombardment triggered the onset of the current style of mantle convection.

Thermal evolution of the earth

NASA Technical Reports Server (NTRS)

Spohn, T.

1984-01-01

The earth's heat budget and models of the earth's thermal evolution are discussed. Sources of the planetary heat are considered and modes of heat transport are addressed, including conduction, convection, and chemical convection. Thermal and convectional models of the earth are covered, and models of thermal evolution are discussed in detail, including changes in the core, the influence of layered mantle convection on the thermal evolution, and the effect of chemical differentiation on the continents.

Formation of plate boundaries: The role of mantle volatilization

NASA Astrophysics Data System (ADS)

Seno, Tetsuzo; Kirby, Stephen H.

2014-02-01

In the early Earth, convection occurred with the accumulation of thick crust over a weak boundary layer downwelling into the mantle (Davies, G.F., 1992. On the emergence of plate tectonics. Geology 20, 963-966.). This would have transitioned to stagnant-lid convection as the mantle cooled (Solomatov, V.S., Moresi, L.-N., 1997. Three regimes of mantle convection with non-Newtonian viscosity and stagnant lid convection on the terrestrial planets. Geophys. Res. Lett. 24, 1907-1910.) or back to a magma ocean as the mantle heated (Sleep, N., 2000. Evolution of the mode of convection within terrestrial planets. J. Geophys. Res. 105(E7): 17563-17578). Because plate tectonics began operating on the Earth, subduction must have been initiated, thus avoiding these shifts. Based on an analogy with the continental crust subducted beneath Hindu Kush and Burma, we propose that the lithosphere was hydrated and/or carbonated by H2O-CO2 vapors released from magmas generated in upwelling plumes and subsequently volatilized during underthrusting, resulting in lubrication of the thrust above, and subduction of the lithosphere along with the overlying thick crust. Once subduction had been initiated, serpentinized forearc mantle may have formed in a wedge-shaped body above a dehydrating slab. In relict arcs, suture zones, or rifted margins, any agent that warms and dehydrates the wedge would weaken the region surrounding it, and form various types of plate boundaries depending on the operating tectonic stress. Thus, once subduction is initiated, formation of plate boundaries might be facilitated by a major fundamental process: weakening due to the release of pressurized water from the warming serpentinized forearc mantle.

Linking the Earth's surface with the deep-mantle plume beneath a region from Iceland to the city of Perm

NASA Astrophysics Data System (ADS)

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

2014-05-01

Current tomography models consistently reveal three large-scale regions of strongly reduced seismic velocity in the lowermost mantle under the Pacific, Africa and a region that extends from below Iceland to the city of Perm (the Perm Anomaly). We have carried out mantle dynamic simulations (Glišović et al., GJI 2012; Glišović & Forte, EPSL 2014) of the evolution of these large-scale structures that directly incorporate: 1) robust constraints provided by 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. The convection simulations provide a detailed insight into the very-long-time evolution of the buoyancy of these lower-mantle anomalies. We find, in particular, that the buoyancy associated with the Perm Anomaly generates a very long-lived superplume that is connected to the paleomagnetic location of the Siberian Traps at the time of their eruption (Smirnov & Tarduno, EPSL 2010) and also to location of North Atlantic Igneous Provinces (i.e., the opening of North Atlantic Ocean).

Magnesium stable isotope composition of Earth's upper mantle

NASA Astrophysics Data System (ADS)

Handler, Monica R.; Baker, Joel A.; Schiller, Martin; Bennett, Vickie C.; Yaxley, Gregory M.

2009-05-01

The mantle is Earth's largest reservoir of Mg containing > 99% of Earth's Mg inventory. However, no consensus exists on the stable Mg isotope composition of the Earth's mantle or how variable it is and, in particular, whether the mantle has the same stable Mg isotope composition as chondrite meteorites. We have determined the Mg isotope composition of olivine from 22 mantle peridotites from eastern Australia, west Antarctica, Jordan, Yemen and southwest Greenland by pseudo-high-resolution MC-ICP-MS on Mg purified to > 99%. The samples include fertile lherzolites, depleted harzburgites and dunites, cryptically metasomatised ('dry') peridotites and modally metasomatised apatite ± amphibole-bearing harzburgites and wehrlites. Olivine from these samples of early Archaean through to Permian lithospheric mantle have δ25Mg DSM-3 = - 0.22 to - 0.08‰. These data indicate the bulk upper mantle as represented by peridotite olivine is homogeneous within current analytical uncertainties (external reproducibility ≤ ± 0.07‰ [2 sd]). We find no systematic δ25Mg variations with location, lithospheric age, peridotite fertility, or degree or nature of mantle metasomatism. Although pyroxene may have slightly heavier δ25Mg than coexisting olivine, any fractionation between mantle pyroxene and olivine is also within current analytical uncertainties with a mean Δ25Mg pyr-ol = +0.06 ± 0.10‰ (2 sd; n = 5). Our average mantle olivine δ25Mg DSM-3 = - 0.14 ± 0.07‰ and δ26Mg DSM-3 = - 0.27 ± 0.14‰ (2 sd) are indistinguishable from the average of data previously reported for terrestrial basalts, confirming that basalts have stable Mg isotope compositions representative of the mantle. Olivine from five pallasite meteorites have δ25Mg DSM-3 = - 0.16 to - 0.11‰ that are identical to terrestrial olivine and indistinguishable from the average δ25Mg previously reported for chondrites. These data provide no evidence for measurable heterogeneity in the stable Mg isotope

Episodic large-scale overturn of two-layer mantles in terrestrial planets

NASA Astrophysics Data System (ADS)

Herrick, D. L.; Parmentier, E. M.

1994-01-01

It is usually assumed that the upper and lower mantles of a chemically stratified planet are arranged so that the upper mantle is chemically less dense and that these layers convect separately. Possible buoyant overturn of the two mantle layers has not previously been considered. Such overturn would initially occur when thermal expansion of a chemically denser lower mantle more than offsets the compositional density difference between the layers, reversing the relative sense of buoyancy. Once overturn has occurred, the chemically denser, but thermally less dense upper mantle cools more efficiently than the lower mantle and loses its relative thermal buoyancy. If mixing is slow, this leads to repeated overturns that result in thermal histories that differ radically from those obtained without this large-scale overturning. Thermal evolution calculations, for a two-layer mantle over a wide range of parameter space, show that large-scale overturn occurs cyclically with a well-defined period. This period depends most strongly on the viscosity of the lower mantle, to which it is approximately proportional. Geologically interesting overturn periods on the order of 107 to 109 years result for lower mantle viscosities of 1022 to 1024 Pa s for the Earth and Venus, and 1021 to 1023 Pa s for Mars. The mantles of Mercury and the Moon are too thin to permit two-layer convection, and therefore the model is not appropriate for them. Overturn cannot occur on Earth or Venus if the compositional density difference between the layers exceeds about 4%, or on Mars if it exceeds about 2%. Large-scale mantle overturn could have significant tectonic consequences such as the initiation of a new plate tectonic cycle on the Earth or a major resurfacing event on Mars or Venus. Such episodic events in the evolution of a planet are not easily explained by whole mantle thermal convection.

Chondritic Xenon in the Earth's mantle: new constrains on a mantle plume below central Europe

NASA Astrophysics Data System (ADS)

Caracausi, Antonio; Avice, Guillaume; Bernard, Peter; Furi, Evelin; Marty, Bernard

2016-04-01

Due to their inertness, their low abundances, and the presence of several different radiochronometers in their isotope systematics, the noble gases are excellent tracers of mantle dynamics, heterogeneity and differentiation with respect to the atmosphere. Xenon deserves particular attention because its isotope systematic can be related to specific processes during terrestrial accretion (e.g., Marty, 1989; Mukhopadhyay, 2012). The origin of heavy noble gases in the Earth's mantle is still debated, and might not be solar (Holland et al., 2009). Mantle-derived CO2-rich gases are particularly powerful resources for investigating mantle-derived noble gases as large quantities of these elements are available and permit high precision isotope analysis. Here, we report high precision xenon isotopic measurements in gases from a CO2 well in the Eifel volcanic region (Germany), where volcanic activity occurred between 700 ka and 11 ka years ago. Our Xe isotope data (normalized to 130Xe) show deviations at all masses compared to the Xe isotope composition of the modern atmosphere. The improved analytical precision of the present study, and the nature of the sample, constrains the primordial Xe end-member as being "chondritic", and not solar, in the Eifel mantle source. This is consistent with an asteroidal origin for the volatile elements in Earth's mantle and it implies that volatiles in the atmosphere and in the mantle originated from distinct cosmochemical sources. Despite a significant fraction of recycled atmospheric xenon in the mantle, primordial Xe signatures still survive in the mantle. This is also a demonstration of a primordial component in a plume reservoir. Our data also show that the reservoir below the Eifel region contains heavy-radiogenic/fissiogenic xenon isotopes, whose ratios are typical of plume-derived reservoirs. The fissiogenic Pu-Xe contribution is 2.26±0.28 %, the UXe contribution is negligible, the remainder being atmospheric plus primordial. Our

Understanding the Earth's Mantle Through Advanced Elasticity Measurements

NASA Astrophysics Data System (ADS)

Marquardt, Hauke; Schulze, Kirsten; Kurnosov, Alexander; Buchen, Johannes; Frost, Daniel; Boffa Ballaran, Tiziana; Marquardt, Katharina; Kawazoe, Takaaki

2017-04-01

Constraints on the inner structure, chemical and mineralogical composition as well as dynamics of Earth's mantle can be derived through comparison of laboratory elasticity data to seismological observables. A quantitative knowledge of the elastic properties of mantle minerals, and their variations with chemical composition, at pressure and temperature conditions of Earth's mantle is key to construct reliable synthetic mineral physics-based seismic velocity models to be compared to seismic observables. We will discuss results of single-crystal elasticity measurements on Earth mantle minerals that have been conducted using the combined Brillouin scattering and x-ray diffraction (XRD) system at BGI Bayreuth in combination with advanced sample preparation using the focused ion beam (FIB) technique [1] that allows for tailoring sizes and shapes of tiny single-crystals. In our experiments, multiple FIB-prepared single-crystals were loaded in a single sample chamber of a resistively-heated diamond-anvil cell (DAC). The possiblity to measure simultaneously acoustic wave velocities and density (unit-cell parameters) in the DAC in combination with the multi-sample approach facilitates direct quantification of the effects of chemical substitution on the elasticity and seismic wave velocities at non-ambient conditions. Our experimental approach eliminates uncertainties arising from the combination of data collected under (potentially) different conditions in several DAC runs, in different laboratories and/or from using different pressure-temperature sensors. We will present our recent experiments on the elasticity of single-crystal Fe-Al-bearing bridgmanite in the lower mantle and discuss implications for the composition and oxidation state of Earth's lower mantle. We will further discuss our laboratory data on the effects of 'water' and iron on the seismic wave velocities of ringwoodite in Earth's transition zone and outline implications for mapping 'water' in the transition

A benchmark initiative on mantle convection with melting and melt segregation

NASA Astrophysics Data System (ADS)

Schmeling, Harro; Dohmen, Janik; Wallner, Herbert; Noack, Lena; Tosi, Nicola; Plesa, Ana-Catalina; Maurice, Maxime

2015-04-01

In recent years a number of mantle convection models have been developed which include partial melting within the asthenosphere, estimation of melt volumes, as well as melt extraction with and without redistribution at the surface or within the lithosphere. All these approaches use various simplifying modelling assumptions whose effects on the dynamics of convection including the feedback on melting have not been explored in sufficient detail. To better assess the significance of such assumptions and to provide test cases for the modelling community we initiate a benchmark comparison. In the initial phase of this endeavor we focus on the usefulness of the definitions of the test cases keeping the physics as sound as possible. The reference model is taken from the mantle convection benchmark, case 1b (Blanckenbach et al., 1989), assuming a square box with free slip boundary conditions, the Boussinesq approximation, constant viscosity and a Rayleigh number of 1e5. Melting is modelled assuming a simplified binary solid solution with linearly depth dependent solidus and liquidus temperatures, as well as a solidus temperature depending linearly on depletion. Starting from a plume free initial temperature condition (to avoid melting at the onset time) three cases are investigated: Case 1 includes melting, but without thermal or dynamic feedback on the convection flow. This case provides a total melt generation rate (qm) in a steady state. Case 2 includes batch melting, melt buoyancy (melt Rayleigh number Rm), depletion buoyancy and latent heat, but no melt percolation. Output quantities are the Nusselt number (Nu), root mean square velocity (vrms) and qm approaching a statistical steady state. Case 3 includes two-phase flow, i.e. melt percolation, assuming a constant shear and bulk viscosity of the matrix and various melt retention numbers (Rt). These cases should be carried out using the Compaction Boussinseq Approximation (Schmeling, 2000) or the full compaction

Water in Earth's mantle: Hydrogen analysis of mantle olivine, pyroxenes and garnet using the SIMS

NASA Technical Reports Server (NTRS)

Kurosawa, Masanori; Yurimoto, Hisayoshi; Sueno, Shigeho

1993-01-01

Hydrogen (or water) in the Earth's interior plays a key role in the evolution and dynamics of the planet. However, the abundance and the existence form of the hydrogen have scarcely been clear in practice. Hydrogen in the mantle was incorporated in the interior during the formation of the Earth. The incorporated hydrogen was hardly possible to concentrate locally inside the Earth considering its high mobility and high reactivity. The hydrogen, preferably, could be distributed homogeneously over the mantle and the core by the subsequent physical and chemical processes. Therefore, hydrogen in the mantle could be present in the form of trace hydrogen in nominally anhydrous mantle minerals. The hydrogen and the other trace elements in mantle olivines, orthopyroxenes, clinopyroxenes, and garnets were determined using secondary ion mass spectrometry (SIMS) for elucidating (1) the exact hydrogen contents, (2) the correlation between the hydrogen and the other trace elements, (3) the dependence of the hydrogen contents on the depth, and (4) the dependence of the whole rock water contents on the depth.

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

Eutectic melting temperature of the lowermost Earth's mantle

NASA Astrophysics Data System (ADS)

Andrault, D.; Lo Nigro, G.; Bolfan-Casanova, N.; Bouhifd, M.; Garbarino, G.; Mezouar, M.

2009-12-01

Partial melting of the Earth's deep mantle probably occurred at different stages of its formation as a consequence of meteoritic impacts and seismology suggests that it even continues today at the core-mantle boundary. Melts are important because they dominate the chemical evolution of the different Earth's reservoirs and more generally the dynamics of the whole planet. Unfortunately, the most critical parameter, that is the temperature profile inside the deep Earth, remains poorly constrained accross the planet history. Experimental investigations of the melting properties of materials representative of the deep Earth at relevant P-T conditions can provide anchor points to refine past and present temperature profiles and consequently determine the degree of melting at the different geological periods. Previous works report melting relations in the uppermost lower mantle region, using the multi-anvil press [1,2]. On the other hand, the pyrolite solidus was determined up to 65 GPa using optical observations in the laser-heated diamond anvil cell (LH-DAC) [3]. Finally, the melting temperature of (Mg,Fe)2SiO4 olivine is documented at core-mantle boundary (CMB) conditions by shock wave experiments [4]. Solely based on these reports, experimental data remain too sparse to draw a definite melting curve for the lower mantle in the relevant 25-135 GPa pressure range. We reinvestigated melting properties of lower mantle materials by means of in-situ angle dispersive X-ray diffraction measurements in the LH-DAC at the ESRF [5]. Experiments were performed in an extended P-T range for two starting materials: forsterite and a glass with chondrite composition. In both cases, the aim was to determine the onset of melting, and thus the eutectic melting temperatures as a function of pressure. Melting was evidenced from drastic changes of diffraction peak shape on the image plate, major changes in diffraction intensities in the integrated pattern, disappearance of diffraction rings

The origin of volatiles in the Earth's mantle

NASA Astrophysics Data System (ADS)

Hier-Majumder, Saswata; Hirschmann, Marc M.

2017-08-01

The Earth's deep interior contains significant reservoirs of volatiles such as H, C, and N. Due to the incompatible nature of these volatile species, it has been difficult to reconcile their storage in the residual mantle immediately following crystallization of the terrestrial magma ocean (MO). As the magma ocean freezes, it is commonly assumed that very small amounts of melt are retained in the residual mantle, limiting the trapped volatile concentration in the primordial mantle. In this article, we show that inefficient melt drainage out of the freezing front can retain large amounts of volatiles hosted in the trapped melt in the residual mantle while creating a thick early atmosphere. Using a two-phase flow model, we demonstrate that compaction within the moving freezing front is inefficient over time scales characteristic of magma ocean solidification. We employ a scaling relation between the trapped melt fraction, the rate of compaction, and the rate of freezing in our magma ocean evolution model. For cosmochemically plausible fractions of volatiles delivered during the later stages of accretion, our calculations suggest that up to 77% of total H2O and 12% of CO2 could have been trapped in the mantle during magma ocean crystallization. The assumption of a constant trapped melt fraction underestimates the mass of volatiles in the residual mantle by more than an order of magnitude.