Wet plume atop of the flattening slab: Insight into intraplate volcanism in East Asia

NASA Astrophysics Data System (ADS)

He, Lijuan

2017-08-01

Geophysical observations imply the intraplate volcanism in East Asia is related to dehydration of slab stagnating in the transition zone. To better understand the dynamics of such process, a thermochemical mantle convection model is constructed to simulate numerically the thermal evolution of slab and the transportation of water in the process of slab downgoing, flattening and stagnation. Equation of water transfer is included, and water effects on density and viscosity are considered. Model results indicate the warming of slab by surrounding mantle is rather slow. Water could be successfully dragged into the transition zone if the reference viscosity of the hydrous layer (with initial water of 2 wt%) is higher than 1017 Pa s and that of mantle is 1021 Pa s. Wet plumes could then originate in the flat-lying part of the slab, relatively far from the trench. Generally, the viscosity of the hydrous layer governs the initiation of wet plume, whereas the viscosity of the overlying mantle wedge controls the activity of the ascending wet plumes - they are more active in the weaker wedge. The complex fluid flow superposed by corner flow and free thermal convection influences greatly the water transport pattern in the upper mantle. Modeling results together with previous modeling infer three stages of water circulation in the big mantle wedge: 1) water is brought into the mantle transition zone by downward subducting slab under some specific thermo-rheological conditions, otherwise water is released at shallow depth near wedge tip; 2) wet plume generates from surface of the flattening slab warmed by surrounding mantle, and 3) water spreads over the big mantle wedge. Wet plume from the flattening Pacific Plate arrives at the lithospheric base and induces melting, which can explain the intraplate Cenozoic volcanoes in East Asia.

Thermal evolution of flattening slab and formation of wet plume: Insight into intraplate volcanism in East Asia

NASA Astrophysics Data System (ADS)

He, L.

2016-12-01

Geophysical observations imply the intraplate volcanism in East Asia is related to dehydration of slab stagnating in the transition zone. To better understand the dynamics of such process, a thermochemical mantle convection model is constructed to simulate numerically the thermal evolution of slab and the transportation of water in the process of subduction. Equation of water transfer is explicitly included, and water effects on density and viscosity are considered. Modeling results indicate that behavior of water transport relates closely to the transient thermal state and viscosities both of the slab and the surrounding mantle. Generally, initiation of wet plume is mainly influenced by the viscosity of the wet layer in the uppermost slab, whereas the horizontal distance of water transport and its ascending rate is affected strongly by the viscosity of the big mantle wedge. Whether water can be carried successfully by slab into the mantle transition zone and trigger wet plume at the surface of flattening slab depends on the viscosity contrast between wet layer and surrounding mantle. The complex fluid flow superposed by corner flow and free thermal convection controls the water transport pattern in the upper mantle. Modeling results together with previous modeling infer three stages of water circulation in the big mantle wedge: 1) water is brought into the mantle transition zone by downward subducting slab when water layer viscosity is much higher than the wedge viscosity, otherwise water is released at shallow depth near wedge tip; 2) wet plume generates from surface of warm flattening slab if containing water, which arrives at the lithospheric base and induces melting; and 3) water spreads all over the big mantle wedge, mantle convection within the big mantle wedge becomes more active, leading to upwelling of asthenosphere and erosion of the overriding continental lithosphere. Wet plume from the flattening Pacific Plate can explain the intraplate Cenozoic volcanoes in East Asia.

Modeling Diverse Pathways to Age Progressive Volcanism in Subduction Zones.

NASA Astrophysics Data System (ADS)

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

2015-12-01

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

Laboratory experiments on subduction-induced circulation in the wedge and the evolution of mantle diapirs

NASA Astrophysics Data System (ADS)

Sylvia, R. T.; Kincaid, C. R.; Behn, M. D.; Zhang, N.

2014-12-01

Circulation in subduction zones involves large-scale, forced-convection by the motion of the down-going slab and small scale, buoyant diapirs of hydrated mantle or subducted sediments. Models of subduction-diapir interaction often neglect large-scale flow patterns induced by rollback, back-arc extension and slab morphology. We present results from laboratory experiments relating these parameters to styles of 4-D wedge circulation and diapir ascent. A glucose fluid is used to represent the mantle. Subducting lithosphere is modeled with continuous rubber belts moving with prescribed velocities, capable of reproducing a large range in downdip relative rollback plate rates. Differential steepening of distinct plate segments simulates the evolution of slab gaps. Back-arc extension is produced using Mylar sheeting in contact with fluid beneath the overriding plate that moves relative to the slab rollback rate. Diapirs are introduced at the slab-wedge interface in two modes: 1) distributions of low density rigid spheres and 2) injection of low viscosity, low density fluid to the base of the wedge. Results from 30 experiments with imposed along-trench (y) distributions of buoyancy, show near-vertical ascent paths only in cases with simple downdip subduction and ratios (W*) of diapir rise velocity to downdip plate rate of W*>1. For W* = 0.2-1, diapir ascent paths are complex, with large (400 km) lateral offsets between source and surfacing locations. Rollback and back-arc extension enhance these offsets, occasionally aligning diapirs from different along-trench locations into trench-normal, age-progressive linear chains beneath the overriding plate. Diapirs from different y-locations may surface beneath the same volcanic center, despite following ascent paths of very different lengths and transit times. In cases with slab gaps, diapirs from the outside edge of the steep plate move 1000 km parallel to the trench before surfacing above the shallow dipping plate. "Dead zones" resulting from lateral and vertical shear in the wedge above the slab gap, produce slow transit times. These 4-D ascent pathways are being incorporated into numerical models on the thermal and melting evolution of diapirs. Models show subduction-induced circulation significantly alters diapir ascent beneath arcs.

Plate and Plume Flux: Constraints for paleomagnetic reference frames and interpretation of deep mantle seismic heterogeneity. (Invited)

NASA Astrophysics Data System (ADS)

Bunge, H.; Schuberth, B. S.; Shephard, G. E.; Müller, D.

2010-12-01

Plate and plume flow are dominant modes of mantle convection, as pointed out by Geoff Davies early on. Driven, respectively, from a cold upper and a hot lower thermal boundary layer these modes are now sufficiently well imaged by seismic tomographers to exploit the thermal boundary layer concept as an effective tool in exploring two long standing geodynamic problems. One relates to the choice of an absolute reference frame in plate tectonic reconstructions. Several absolute reference frames have been proposed over the last decade, including those based on hotspot tracks displaying age progression and assuming either fixity or motion, as well as palaeomagnetically-based reference frames, a subduction reference frame and hybrid versions. Each reference frame implies a particular history of the location of subduction zones through time and thus the evolution of mantle heterogeneity via mixing of subducted slab material in the mantle. Here we compare five alternative absolute plate motion models in terms of their consequences for deep mantle structure. Taking global paleo-plate boundaries and plate velocities back to 140 Ma derived from the new plate tectonic reconstruction software GPlates and assimilating them into vigorous 3-D spherical mantle circulation models, we infer geodynamic mantle heterogeneity and compare it to seismic tomography for each absolute rotation model. We also focus on the challenging problem of interpreting deep mantle seismic heterogeneity in terms of thermal and compositional variations. Using published thermodynamically self-consistent mantle mineralogy models in the pyrolite composition, we find strong plume flux from the CMB, with a high temperature contrast (on the order of 1000 K) across the lower thermal boundary layer is entirely sufficient to explain elastic heterogeneity in the deep mantle for a number of quantitative measures. A high excess temperatures of +1000--1500 K for plumes in the lowermost mantle is particularly important in understanding the strong seismic velocity reduction mapped by tomography in low-velocity bodies of the deep mantle, as this produces significant negative anomalies of shear wave velocity of up to -4%. We note, however, that our results do not account for the curious observation of seismic anti-correlation, which appears difficult to explain in any case. Our results provide important constraints for the integration of plate tectonics and mantle dynamics and their use in forward and inverse geodynamic mantle models.

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.

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.

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.

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.

Numerical modeling the genetic mechanism of Cenozoic intraplate Volcanoes in Northeastern China

NASA Astrophysics Data System (ADS)

Qu, Wulin; Chen, Yongshun John; Zhang, Huai; Jin, Yimin; Shi, Yaolin

2017-04-01

Changbaishan Volcano located about 1400 km west of Japan Trench is an intra continental volcano which having different origin from island arc volcanoes. A number of different mechanisms have been proposed to interpret the origin of intraplate volcanoes, such as deep mantle plumes, back-arc extension and decompressional partial melting, asthenosphere upwelling and decompressional melting, and deep stagnant slab dehydration and partial melting. The recent geophysical research reveals that the slow seismic velocity anomaly extends continuously just below 660 km depth to surface beneath Changbaishan by seismic images and three-dimensional waveform modelling [Tang et al., 2014]. The subduction-induced upwelling occurs within a gap in the stagnant subducted Pacific Plate and produces decompressional melting. Water in deep Earth can reduce viscosity and lower melting temperature and seismic velocity and has effects on many other physical properties of mantle materials. The water-storage capacity of wadsleyite and ringwoodite, which are the main phase in the mantle transition zone, is much greater than that of upper mantle and lower mantle. Geophysical evidences have shown that water content in the mantle transition zone is exactly greater than that of upper mantle and lower mantle [Karato, 2011]. Subducted slab could make mantle transition zone with high water content upward or downward across main phase change surface to release water, and lead to partial melting. We infer that the partial melting mantle and subducted slab materials propagate upwards and form the Cenozoic intraplate Volcanoes in Northeastern China. We use the open source code ASPECT [Kronbichler et al., 2012] to simulate the formation and migration of magma contributing to Changbaishan Volcano. We find that the water entrained by subducted slab from surface has only small proportion comparing to water content of mantle transition zone. Our model provide insights into dehydration melting induced by water transport out of the mantle transition zone associated with dynamic interactions between the subducted slab and surrounding mantle. References Karato, S. (2011), Water distribution across the mantle transition zone and its implications for global material circulation, EARTH PLANET SC LETT, 301(3), 413-423. Kronbichler, M., et al. (2012), High accuracy mantle convection simulation through modern numerical methods, GEOPHYS J INT, 191(1), 12-29. Tang, Y., et al. (2014), Changbaishan volcanism in northeast China linked to subduction-induced mantle upwelling, NAT GEOSCI, 7(6), 470-475.

Bifurcation of the Yellowstone plume driven by subduction-induced mantle flow

NASA Astrophysics Data System (ADS)

Kincaid, C.; Druken, K. A.; Griffiths, R. W.; Stegman, D. R.

2013-05-01

The causes of volcanism in the northwestern United States over the past 20 million years are strongly contested. Three drivers have been proposed: melting associated with plate subduction; tectonic extension and magmatism resulting from rollback of a subducting slab; or the Yellowstone mantle plume. Observations of the opposing age progression of two neighbouring volcanic chains--the Snake River Plain and High Lava Plains--are often used to argue against a plume origin for the volcanism. Plumes are likely to occur near subduction zones, yet the influence of subduction on the surface expression of mantle plumes is poorly understood. Here we use experiments with a laboratory model to show that the patterns of volcanism in the northwestern United States can be explained by a plume upwelling through mantle that circulates in the wedge beneath a subduction zone. We find that the buoyant plume may be stalled, deformed and partially torn apart by mantle flow induced by the subducting plate. Using plausible model parameters, bifurcation of the plume can reproduce the primary volcanic features observed in the northwestern United States, in particular the opposite progression of two volcanic chains. Our results support the presence of the Yellowstone plume in the northwestern United States, and also highlight the power of plume-subduction interactions to modify surface geology at convergent plate margins.

Geophysical signature of hydration-dehydration processes in active subduction zones

NASA Astrophysics Data System (ADS)

Reynard, Bruno

2013-04-01

Seismological and magneto-telluric tomographies are potential tools for imaging fluid circulation when combined with petrophysical models. Recent measurements of the physical properties of serpentine allow refining hydration of the mantle and fluid circulation in the mantle wedge from geophysical data. In the slab lithospheric mantle, serpentinization caused by bending at the trench is limited to a few kilometers below the oceanic crust (<5 km). Double Wadati-Benioff zones, 20-30 km below the crust, are explained by deformation of dry peridotites, not by serpentine dehydration. It reduces the required amount of water stored in solid phases in the slab (Reynard et al., 2010). In the cold (<700°C) fore-arc mantle wedge above the subducting slab, serpentinization is caused by the release of large amounts of hydrous fluids in the cold mantle above the dehydrating subducted plate. Low seismic velocities in the wedge give a time-integrated estimate of hydration and serpentinization. Serpentinization reaches 50-100% in hot subduction, while it is below 10% in cold subduction (Bezacier et al., 2010; Reynard, 2012). Electromagnetic profiles of the mantle wedge reveal high electrical-conductivity bodies. In hot areas of the mantle wedge (> 700°C), water released by dehydration of the slab induces melting of the mantle under volcanic arcs, explaining the observed high conductivities. In the cold melt-free wedge (< 700°C), high conductivities in electromagnetic profiles provide "instantaneous" images of fluid circulation because the measured electrical conductivity of serpentine is below 0.1 mS/m (Reynard et al., 2011). A small fraction (ca. 1% in volume) of connective high-salinity fluids accounts for the highest observed conductivities. Low-salinity fluids (≤ 0.1 m) released by slab dehydration evolve towards high-salinity (≥ 1 m) fluids during progressive serpentinization in the wedge. These fluids can mix with arc magmas at depths and account for high-chlorine melt inclusions in arc lavas. High electrical conductivities up to 1 S/m in the hydrated wedge of the hot subductions (Ryukyu, Kyushu, Cascadia) reflect high fluid concentration, while low to moderate (<0.01 S/m) conductivities in the cold subductions (N-E Japan, Bolivia) reflect low fluid flow. This is consistent with the seismic observations of extensive shallow serpentinization in hot subduction zones, while serpentinization is sluggish in cold subduction zones. Bezacier, L., et al. 2010. Elasticity of antigorite, seismic detection of serpentinites, and anisotropy in subduction zones. Earth and Planetary Science Letters, 289, 198-208. Reynard, B., 2012. Serpentine in active subduction zones. Lithos, http://dx.doi.org/10.1016/j.lithos.2012.10.012. Reynard, B., Mibe, K. & Van de Moortele, B., 2011. Electrical conductivity of the serpentinised mantle and fluid flow in subduction zones. Earth and Planetary Science Letters, 307, 387-394. Reynard, B., Nakajima, J. & Kawakatsu, H., 2010. Earthquakes and plastic deformation of anhydrous slab mantle in double Wadati-Benioff zones. Geophysical Research Letters, 37, L24309.

Machine Learning and Inverse Problem in Geodynamics

NASA Astrophysics Data System (ADS)

Shahnas, M. H.; Yuen, D. A.; Pysklywec, R.

2017-12-01

During the past few decades numerical modeling and traditional HPC have been widely deployed in many diverse fields for problem solutions. However, in recent years the rapid emergence of machine learning (ML), a subfield of the artificial intelligence (AI), in many fields of sciences, engineering, and finance seems to mark a turning point in the replacement of traditional modeling procedures with artificial intelligence-based techniques. The study of the circulation in the interior of Earth relies on the study of high pressure mineral physics, geochemistry, and petrology where the number of the mantle parameters is large and the thermoelastic parameters are highly pressure- and temperature-dependent. More complexity arises from the fact that many of these parameters that are incorporated in the numerical models as input parameters are not yet well established. In such complex systems the application of machine learning algorithms can play a valuable role. Our focus in this study is the application of supervised machine learning (SML) algorithms in predicting mantle properties with the emphasis on SML techniques in solving the inverse problem. As a sample problem we focus on the spin transition in ferropericlase and perovskite that may cause slab and plume stagnation at mid-mantle depths. The degree of the stagnation depends on the degree of negative density anomaly at the spin transition zone. The training and testing samples for the machine learning models are produced by the numerical convection models with known magnitudes of density anomaly (as the class labels of the samples). The volume fractions of the stagnated slabs and plumes which can be considered as measures for the degree of stagnation are assigned as sample features. The machine learning models can determine the magnitude of the spin transition-induced density anomalies that can cause flow stagnation at mid-mantle depths. Employing support vector machine (SVM) algorithms we show that SML techniques can successfully predict the magnitude of the mantle density anomalies and can also be used in characterizing mantle flow patterns. The technique can be extended to more complex problems in mantle dynamics by employing deep learning algorithms for estimation of mantle properties such as viscosity, elastic parameters, and thermal and chemical anomalies.

Effects of compressibility on the temperature jump at the interface of layered, spherical-shell convection

NASA Technical Reports Server (NTRS)

Yen, David A.; Zhang, Shuxia; Langenberger, Sherri E.

1988-01-01

Large temperature jumps at the interface of layered convection are important to the argument used against the likelihood of separate circulations in the upper and lower mantles. This problem was studied within the framework of a compressible, constant viscosity spherical-shell model. Both mechanical and thermal coupling configurations are considered. Although the temperature jumps are reduced by compressibility, their magnitudes remain quite large, in the case of mechanical coupling. For thermal coupling, the temperature jumps become smaller but still are substantial, between 500 to 1000 C. In layered spherical-shell convection, flows in the lower mantle are several times greater than the surface velocities.

Profiling the robustness, efficiency and limits of the forward-adjoint method for 3-D mantle convection modelling

NASA Astrophysics Data System (ADS)

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

2018-02-01

Knowledge of Earth's past mantle structure is inherently unknown. This lack of knowledge presents problems in many areas of Earth science, including in mantle circulation modelling (MCM). As a mathematical model of mantle convection, MCMs require boundary and initial conditions. While boundary conditions are readily available from sources such as plate reconstructions for the upper surface, and as free slip at the core-mantle boundary, the initial condition is not known. MCMs have historically `created' an initial condition using long `spin up' processes using the oldest available plate reconstruction period available. While these do yield good results when models are run to present day, it is difficult to infer with confidence results from early in a model's history. Techniques to overcome this problem are now being studied in geodynamics, such as by assimilating the known internal structure (e.g. from seismic tomography) of Earth at present day backwards in time. One such method is to use an iterative process known as the forward-adjoint method. While this is an efficient means of solving this inverse problem, it still strains all but the most cutting edge computational systems. In this study we endeavour to profile the effectiveness of this method using synthetic test cases as our known data source. We conclude that savings in terms of computational expense for forward-adjoint models can be achieved by streamlining the time-stepping of the calculation, as well as determining the most efficient method of updating initial conditions in the iterative scheme. Furthermore, we observe that in the models presented, there exists an upper limit on the time interval over which solutions will practically converge, although this limit is likely to be linked to Rayleigh number.

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.

Model of Auctioneer Estimation of Swordtip Squid (Loligo edulis) Quality

NASA Astrophysics Data System (ADS)

Nakamura, Makoto; Matsumoto, Keisuke; Morimoto, Eiji; Ezoe, Satoru; Maeda, Toshimichi; Hirano, Takayuki

The knowledge of experienced auctioneers regarding the circulation of marine products is an essential skill and is necessary for evaluating product quality and managing aspects such as freshness. In the present study, the ability of an auctioneer to quickly evaluate the freshness of swordtip squid (Loligo edulis) at fish markets was analyzed. Evaluation characteristics used by an auctioneer were analyzed and developed using a fuzzy logic model. Forty boxes containing 247 swordtip squid with mantles measuring 220 mm that had been evaluated and assigned to one of five quality categories by an auctioneer were used for the analysis and the modeling. The relationships between the evaluations of appearance, body color, and muscle freshness were statistically analyzed. It was found that a total of four indexes of the epidermis color strongly reflected evaluations of appearance: dispersion ratio of the head, chroma on the head-end mantle and the difference in the chroma and brightness of the mantle. The fuzzy logic model used these indexes for the antecedent-part of the linguistic rules. The results of both simulation and evaluations demonstrate that the model is robust, with the predicted results corresponding with more than 96% of the quality assignments of the auctioneers.

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.

Mantle Degassing and Atmosphere Evolution

NASA Astrophysics Data System (ADS)

Zhang, Y.

2011-12-01

Noble gas isotopes have provided much of our understanding of Earth's early history [1-3]. Various degassing models have been developed, including degassing of the whole mantle, degassing of all gases at similar relative rate [1], solubility-controlled degassing [2], and steady-state degassing models [4]. This report will evaluate various degassing models using recent data. For example, helium outgassing flux has been lowered by more than a factor of two based on sophisticated ocean general circulation models [5], which also impacts on the estimated degassing flux of carbon. Years of measurements and progress have allowed isotopic ratios of various mantle reservoirs being pieced together [6]. For example, 129Xe/130Xe in OIB mantle is found to be lower than that in MORB mantle [7]. Missing Xe has been found to be a non-issue [8]. Nucleogenic 21Ne production rate relative radiogenic 4He has been revised [9-10], which leads to an interesting neon paradox that nucleogenic 21Ne production in the whole silicate Earth is barely enough to supply nucleogenic 21Ne in air. 40Ar/36Ar ratio in BSE seems to be much lower than any OIB samples, another interesting paradox. Although non-nucleogenic mantle neon is solar, nonradiogenic mantle argon is atmospheric [11]. For Kr and Xe, the jury is still out. When mantle degassing models are evaluated using volatile data of the MORB and OIB, solubility-controlled degassing is able to reconcile more data than other degassing models. On the other hand, the vailable data seem to indicate that atmosphere evolution is more than mantle degassing; there may be significant contribution to the atmosphere from impact degassing and other sources. Furthermore, we are now suffering from too many data so that understanding the whole picture is elusive. [1] Allegre et al. (1986/87) EPSL 81, 127-150. [2] Zhang & Zindler (1989) J. Geophys. Res. 94, 13719-13737. [3] Zhang (1998) Geochim. Cosmochim. Acta 62, 3185-3189. [4] Pocelli & Wasserburg (1995) Geochim. Cosmochim. Acta 59, 4921-4937. [5] Bianchi et al. (2010) EPSL 297, 379-386. [6] Jackson et al. (2009) EPSL 287, 519-528. [7] Graham (2002) Rev. Mineral. Geochem. 47, 247-317. [8] Zhang (2002) Earth-Sci. Rev. 59, 235-263. [9] Yatsevich and Honda (1997).

Circulation of carbon dioxide in the mantle: multiscale modeling

NASA Astrophysics Data System (ADS)

Morra, G.; Yuen, D. A.; Lee, S.

2012-12-01

Much speculation has been put forward on the quantity and nature of carbon reservoirs in the deep Earth, because of its involvement in the evolution of life at the surface and inside planetary interiors. Carbon penetrates into the Earth's mantle mostly during subduction of oceanic crust, which contains carbonate deposits [1], however the form that it assumes at lower mantle depths is scarcely understood [2], hampering our ability to estimate the amount of carbon in the entire mantle by orders of magnitude. We present simulations of spontaneous degassing of supercritical CO2 using in-house developed novel implementations of the Fast-Multipole Boundary Element Method suitable for modeling two-phase flow (here mantle mineral and free CO2 fluid) through disordered materials such as porous rocks. Because the mutual interaction of droplets immersed either in a fluid or a solid matrix and their weakening effect to the host rock alters the strength of the mantle rocks, at the large scale the fluid phases in the mantle may control the creeping of mantle rocks [3]. In particular our study focuses on the percolation of supercritical CO2, estimated through the solution of the Laplace equation in a porous system, stochastically generated through a series of random Karhunen-Loeve decomposition. The model outcome is employed to extract the transmissivity of supercritical fluids in the mantle from the lowest scale up to the mantle scale and in combination with the creeping flow of the convecting mantle. The emerging scenarios on the global carbon cycle are finally discussed. [1] Boulard, E., et al., New host for carbon in the deep Earth. Proceedings of the National Academy of Sciences, 2011. 108(13): p. 5184-5187. [2] Walter, M.J., et al., Deep Mantle Cycling of Oceanic Crust: Evidence from Diamonds and Their Mineral Inclusions. Science, 2011. 334(6052): p. 54-57. [3] Morra, G., et al., Ascent of Bubbles in Magma Conduits Using Boundary Elements and Particles. Procedia Computer Science, 2011.; Boundary Element solution of a flow through a porous. Left boxes represent the the matrix associated with the integrals. The flow enters below and emerges at the top, the amount of flow is identical. The flow is spread in the porous and is viscousless (Laplace equation).

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.

Inferring nonlinear mantle rheology from the shape of the Hawaiian swell.

PubMed

Asaadi, N; Ribe, N M; Sobouti, F

2011-05-26

The convective circulation generated within the Earth's mantle by buoyancy forces of thermal and compositional origin is intimately controlled by the rheology of the rocks that compose it. These can deform either by the diffusion of point defects (diffusion creep, with a linear relationship between strain rate and stress) or by the movement of intracrystalline dislocations (nonlinear dislocation creep). However, there is still no reliable map showing where in the mantle each of these mechanisms is dominant, and so it is important to identify regions where the operative mechanism can be inferred directly from surface geophysical observations. Here we identify a new observable quantity--the rate of downstream decay of the anomalous seafloor topography (swell) produced by a mantle plume--which depends only on the value of the exponent in the strain rate versus stress relationship that defines the difference between diffusion and dislocation creep. Comparison of the Hawaiian swell topography with the predictions of a simple fluid mechanical model shows that the swell shape is poorly explained by diffusion creep, and requires a dislocation creep rheology. The rheology predicted by the model is reasonably consistent with laboratory deformation data for both olivine and clinopyroxene, suggesting that the source of Hawaiian lavas could contain either or both of these components.

Comparison of gravimetric and mantle flow solutions for sub-lithopsheric stress modeling and their combination

NASA Astrophysics Data System (ADS)

Eshagh, Mehdi; Steinberger, Bernhard; Tenzer, Robert; Tassara, Andrés

2018-05-01

Based on Hager and O'Connell's solution to mantle flow equations, the stresses induced by mantle convection are determined using the density and viscosity structure in addition to topographic data and a plate velocity model. The solution to mantle flow equations requires the knowledge of mantle properties that are typically retrieved from seismic information. Large parts of the world are, however, not yet covered sufficiently by seismic surveys. An alternative method of modeling the stress field was introduced by Runcorn. He formulated a direct relation between the stress field and gravity data, while adopting several assumptions, particularly disregarding the toroidal mantle flow component and mantle viscosity variations. A possible way to overcome theoretical deficiencies of Runcorn's theory as well as some practical limitations of applying Hager and O'Connell's theory (in the absence of seismic data) is to combine these two methods. In this study, we apply a least-squares analysis to combine these two methods based on the gravity data inversion constraint on mantle flow equations. In particular, we use vertical gravity gradients from the Gravity field and steady state Ocean Circulation Explorer that are corrected for the gravitational contribution of crustal density heterogeneities prior to applying a localized gravity-gradient inversion. This gravitational contribution is estimated based on combining the Vening Meinesz-Moritz and flexural isostatic theories. Moreover, we treat the non-isostatic effect implicitly by applying a band-limited kernel of the integral equation during the inversion. In numerical studies of modeling, the stress field within the South American continental lithosphere we compare the results obtained after applying Runcorn and Hager and O'Connell's methods as well as their combination. The results show that, according to Hager and O'Connell's (mantle flow) solution, the maximum stress intensity is inferred under the northern Andes. Additional large stress anomalies are detected along the central and southern Andes, while stresses under most of old, stable cratonic formations are much less pronounced or absent. A prevailing stress-vector orientation realistically resembles a convergent mantle flow and downward currents under continental basins that separate Andean Orogeny from the Amazonian Shield and adjacent cratons. Runcorn's (gravimetric) solution, on the other hand, reflects a tectonic response of the lithosphere to mantle flow, with the maximum stress intensity detected along the subduction zone between the Nazca and Altiplano plates and along the convergent tectonic margin between the Altiplano and South American plates. The results also reveal a very close agreement between the results obtained from the combined and Hager and O'Connell's solutions.

Water Content of Earth's Continental Mantle Is Controlled by the Circulation of Fluids or Melts

NASA Technical Reports Server (NTRS)

Peslier, Anne; Woodland, Alan B.; Bell, David R.; Lazarov, Marina; Lapen, Thomas J.

2014-01-01

A key mission of the ARES Directorate at JSC is to constrain models of the formation and geological history of terrestrial planets. Water is a crucial parameter to be measured with the aim to determine its amount and distribution in the interior of Earth, Mars, and the Moon. Most of that "water" is not liquid water per se, but rather hydrogen dissolved as a trace element in the minerals of the rocks at depth. Even so, the middle layer of differentiated planets, the mantle, occupies such a large volume and mass of each planet that when it is added at the planetary scale, oceans worth of water could be stored in its interior. The mantle is where magmas originate. Moreover, on Earth, the mantle is where the boundary between tectonic plates and the underlying asthenosphere is located. Even if mantle rocks in Earth typically contain less than 200 ppm H2O, such small quantities have tremendous influence on how easily they melt (i.e., the more water there is, the more magma is produced) and deform (the more water there is, the less viscous they are). These two properties alone emphasize that to understand the distribution of volcanism and the mechanism of plate tectonics, the water content of the mantle must be determined - Earth being a template to which all other terrestrial planets can be compared.

Numerical modelling of volatiles in the deep mantle

NASA Astrophysics Data System (ADS)

Eichheimer, Philipp; Thielmann, Marcel; Golabek, Gregor J.

2017-04-01

The transport and storage of water in the mantle significantly affects several material properties of mantle rocks and thus water plays a key role in a variety of geodynamical processes (tectonics, magmatism etc.). The processes driving transport and circulation of H2O in subduction zones remain a debated topic. Geological and seismological observations suggest different inflow mechanisms of water e.g. slab bending, thermal cracking and serpentinization (Faccenda et al., 2009; Korenaga, 2017), followed by dehydration of the slab. On Earth both shallow and steep subduction can be observed (Li et al., 2011). However most previous models (van Keken et al., 2008; Wilson et al., 2014) did not take different dip angles and subduction velocities of slabs into account. To which extent these parameters and processes influence the inflow of water still remains unclear. We present 2D numerical models simulating the influence of the various water inflow mechanisms on the mantle with changing dip angle and subduction velocity of the slab over time. The results are used to make predictions regarding the rheological behavior of the mantle wedge, dehydration regimes and volcanism at the surface. References: van Keken, P. E., et al. A community benchmark for subduction zone modeling. Phys. Earth Planet. Int. 171, 187-197 (2008). Faccenda, M., T.V. Gerya, and L. Burlini. Deep slab hydration induced by bending-related variations in tectonic pressure. Nat. Geosci. 2, 790-793 (2009). Korenaga, J. On the extent of mantle hydration caused by plate bending. Earth Planet. Sci. Lett. 457, 1-9 (2017). Wilson, C. R., et al. Fluid flow in subduction zones: The role of solid rheology and compaction pressure. Earth Planet. Sci. Lett. 401, 261-274 (2014). Li, Z. H., Z. Q. Xu, and T. V. Gerya. Flat versus steep subduction: Contrasting modes for the formation and exhumation of high- to ultrahigh-pressure rocks in continental collision zones. Earth Planet. Sci. Lett. 301, 65-77 (2011).

Shaping mobile belts by small-scale convection.

PubMed

Faccenna, Claudio; Becker, Thorsten W

2010-06-03

Mobile belts are long-lived deformation zones composed of an ensemble of crustal fragments, distributed over hundreds of kilometres inside continental convergent margins. The Mediterranean represents a remarkable example of this tectonic setting: the region hosts a diffuse boundary between the Nubia and Eurasia plates comprised of a mosaic of microplates that move and deform independently from the overall plate convergence. Surface expressions of Mediterranean tectonics include deep, subsiding backarc basins, intraplate plateaux and uplifting orogenic belts. Although the kinematics of the area are now fairly well defined, the dynamical origins of many of these active features are controversial and usually attributed to crustal and lithospheric interactions. However, the effects of mantle convection, well established for continental interiors, should be particularly relevant in a mobile belt, and modelling may constrain important parameters such as slab coherence and lithospheric strength. Here we compute global mantle flow on the basis of recent, high-resolution seismic tomography to investigate the role of buoyancy-driven and plate-motion-induced mantle circulation for the Mediterranean. We show that mantle flow provides an explanation for much of the observed dynamic topography and microplate motion in the region. More generally, vigorous small-scale convection in the uppermost mantle may also underpin other complex mobile belts such as the North American Cordillera or the Himalayan-Tibetan collision zone.

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 for an accurate prediction of instantaneous flow, but not for a prediction after 10 My of evolution. Therefore, inverse methods (sequential or data assimilation methods) using short-term fully dynamic evolution that predict surface kinematics are promising tools for a better understanding of the state of the Earth's mantle.

Plumes do not Exist: Plate Circulation is Confined to Upper Mantle

NASA Astrophysics Data System (ADS)

Hamilton, W. B.

2002-12-01

Plumes from deep mantle are widely conjectured to define an absolute reference frame, inaugurate rifting, drive plates, and profoundly modify oceans and continents. Mantle properties and composition are assumed to be whatever enables plumes. Nevertheless, purported critical evidence for plume speculation is false, and all data are better interpreted without plumes. Plume fantasies are made ever more complex and ad hoc to evade contradictory data, and have no predictive value because plumes do not exist. All plume conjecture derives from Hawaii and the guess that the Emperor-Hawaii inflection records a 60-degree change in Pacific plate direction at 45 Ma. Paleomagnetic latitudes and smooth Pacific spreading patterns disprove any such change. Rationales for other fixed plumes collapse when tested, and hypotheses of jumping, splitting, and gyrating plumes are specious. Thermal and physical properties of Hawaiian lithosphere falsify plume predictions. Purported tomographic support elsewhere represents artifacts and misleading presentations. Asthenosphere is everywhere near solidus temperature, so melt needs a tensional setting for egress but not local heat. Gradational and inconsistent contrasts between MORB and OIB are as required by depth-varying melt generation and behavior in contrasted settings and do not indicate systematically unlike sources. MORB melts rise, with minimal reaction, through hot asthenosphere, whereas OIB melts react with cool lithosphere, and lose mass, by crystallizing refractories and retaining and assimilating fusibles. The unfractionated lower mantle of plume conjecture is contrary to cosmologic and thermodynamic data, for mantle below 660 km is more refractory than that above. Subduction, due to density inversion by top-down cooling that forms oceanic lithosphere, drives plate tectonics and upper-mantle circulation. It organizes plate motions and lithosphere stress, which controls plate boundaries and volcanic chains. Hinge rollback is the key to kinematics. Arcs advance and collide, fast-spreading Pacific shrinks, etc. A fore-arc basin atop an overriding plate shows that hinge and non-shortening plate front there track together: velocities of rollback and advance are equal. Convergence velocity commonly also equals rollback velocity but often is greater. Slabs sinking broadside push upper mantle back under incoming plates and force rapid Pacific spreading, whereas overriding plates flow forward with retreating hinges. Backarc basins open behind island arcs migrating with hinges. Slabs settle on uncrossable 660-km discontinuity. (Contrary tomographic claims reflect sampling and smearing artifacts, notably due to along-slab raypaths.) Plates advance over sunken slabs and mantle displaced rearward by them, and ridges spread where advancing plates pull away. Ridges migrate over asthenosphere, producing geophysical and bathymetric asymmetry, and tap fresh asthenosphere into which slab material is recycled upward. Sluggish deep-mantle circulation is decoupled from rapid upper-mantle circulation, so plate motions can be referenced to semistable lower mantle. Global plate motions make kinematic sense if Antarctica, almost ringed by departing ridges and varying little in Cenozoic paleomagnetic position, is stationary: hinges roll back, ridges migrate, and directions and velocities of plate rotations accord with subduction, including sliding and crowding of oceanic lithosphere toward free edges, as the dominant drive. (The invalid hotspot and no-net-rotation frames minimize motions of hinges and ridges, and their plate motions lack kinematic sense.) Northern Eurasia also is almost stationary, Africa rotates very slowly counterclockwise toward Aegean and Zagros, Pacific plate races toward surface-exit subduction systems, etc.

Global plate motion frames: Toward a unified model

NASA Astrophysics Data System (ADS)

Torsvik, Trond H.; Müller, R. Dietmar; van der Voo, Rob; Steinberger, Bernhard; Gaina, Carmen

2008-09-01

Plate tectonics constitutes our primary framework for understanding how the Earth works over geological timescales. High-resolution mapping of relative plate motions based on marine geophysical data has followed the discovery of geomagnetic reversals, mid-ocean ridges, transform faults, and seafloor spreading, cementing the plate tectonic paradigm. However, so-called "absolute plate motions," describing how the fragments of the outer shell of the Earth have moved relative to a reference system such as the Earth's mantle, are still poorly understood. Accurate absolute plate motion models are essential surface boundary conditions for mantle convection models as well as for understanding past ocean circulation and climate as continent-ocean distributions change with time. A fundamental problem with deciphering absolute plate motions is that the Earth's rotation axis and the averaged magnetic dipole axis are not necessarily fixed to the mantle reference system. Absolute plate motion models based on volcanic hot spot tracks are largely confined to the last 130 Ma and ideally would require knowledge about the motions within the convecting mantle. In contrast, models based on paleomagnetic data reflect plate motion relative to the magnetic dipole axis for most of Earth's history but cannot provide paleolongitudes because of the axial symmetry of the Earth's magnetic dipole field. We analyze four different reference frames (paleomagnetic, African fixed hot spot, African moving hot spot, and global moving hot spot), discuss their uncertainties, and develop a unifying approach for connecting a hot spot track system and a paleomagnetic absolute plate reference system into a "hybrid" model for the time period from the assembly of Pangea (˜320 Ma) to the present. For the last 100 Ma we use a moving hot spot reference frame that takes mantle convection into account, and we connect this to a pre-100 Ma global paleomagnetic frame adjusted 5° in longitude to smooth the reference frame transition. Using plate driving force arguments and the mapping of reconstructed large igneous provinces to core-mantle boundary topography, we argue that continental paleolongitudes can be constrained with reasonable confidence.

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

NASA Astrophysics Data System (ADS)

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

2015-12-01

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

The proximity of hotspots to convergent and divergent plate boundaries

NASA Technical Reports Server (NTRS)

Weinstein, Stuart A.; Olson, Peter L.

1989-01-01

An analysis of four different hotspot distributions, ranging from Morgan's (1972) original list of 19 to Vogt's (1981) list of 117 reveals that the hotspots are preferentially located near divergent plate boundaries. The probability of this proximity occurring by chance alone is quite remote, less than 0.01 for all four hotspot distributions. The same analysis also reveals that the hotspots are preferentially excluded from regions near convergent plate boundaries. The probability of this exclusion occurring by chance alone is 0.1 or less for three out of the four distributions examined. We interpret this behavior as being a consequence of the effects of large scale convective circulation on ascending mantle plumes. Mantle thermal plumes, the most probable source of hotspots, arise from instabilities in a basal thermal boundary layer. Plumes are suppressed from regions beneath convergent boundaries by descending flow and are entrained into the upwelling flow beneath spreading centers. Plate-scale convective circulation driven by subduction may also advect mantle thermal plumes toward spreading centers.

Continental Subduction: Mass Fluxes and Interactions with the Wider Earth System

NASA Astrophysics Data System (ADS)

Cuthbert, S. J.

2011-12-01

Substantial parts of ultra-high pressure (UHP) terrains probably represent subducted passive continental margins (PCM). This contribution reviews and synthesises research on processes operating in such systems and their implication for the wider Earth system. PCM sediments are large repositories of volatiles including hydrates, nitrogen species, carbonates and hydrocarbons. Sediments and upper/ mid-crustal basement are rich in incompatible elements and are fertile for melting. Lower crust may be more mafic and refractory. Juvenile rift-related mafic rocks also have the potential to generate substantial volumes of granitoid melts, especially if they have been hydrated. Exposed UHP terrains demonstrate the return of continental crust from mantle depths, show evidence for substantial fluxes of aqueous fluid, anatexis and, in entrained orogenic peridotites, metasomatism of mantle rocks by crust- derived C-O-H fluids. However, substantial bodies of continental material may never return to the surface as coherent masses of rock, but remain sequestered in the mantle where they melt or become entrained in the deeper mantle circulation. Hence during subduction, PCM's become partitioned by a range of mechanisms. Mechanical partitioning strips away weaker sediment and middle/upper crust, which circulate back up the subduction channel, while denser, stronger transitional pro-crust and lower crust may "stall" near the base of the lithosphere or be irreversibly subducted to join the global mantle circulation. Under certain conditions sediment and upper crustal basement may reach depths for UHPM. Further partitioning takes place by anatexis, which either aids stripping and exhumation of the more melt-prone rock-masses through mechanical softening, or separates melt from residuum so that melt escapes and is accreted to the upper plate leading to "undercrusting", late-orogenic magmatism and further refinement of the crust. Melt that traverses sections of mantle will interact with it causing metasomatism and refertilisation. Partitioning also takes place by solid-fluid and melt-fluid partitioning. Dehydration may take place both during subduction and exhumation, and fluxes between dehydrating and hydrating rock masses influence the internal fluid budget of the orogen (essential for eclogitisation and densification of mafic lithologies). Ascending granitic melts advect dissolved water to shallow levels, or even the atmosphere. Irreversible subduction of PCM sediment carries water plus nitrogen species to the deeper mantle. Decarbonation of voluminous PCM carbonates depends on thermal regime and may release a pulse of CO2 to the atmosphere, but is limited in colder subduction zones hence transferring large volumes of carbon to the deep mantle. This may ultimately be mobilised by melting or dissolution to form fluid media for diamond formation.

Plate-Tectonic Circulation is Driven by Cooling From the Top and is Closed Within the Upper Mantle

NASA Astrophysics Data System (ADS)

Hamilton, W. B.

2001-12-01

Subduction drives plate tectonics and is due to cooling from the top: circulation is self-organized, and likely is closed above the discontinuity near 660 km. The contrary consensus that plate tectonics is driven by bottom heating and involves the entire mantle combines misunderstood kinematics with flawed concepts of through-the-mantle plumes and subduction. Plume conjecture came from the Emperor-Hawaii progression, the 45 Ma inflection in which was assumed to mark a 60-degree change in direction of that part of the Pacific plate over a fixed plume. Smooth spreading patterns around the east and south margin of the Pacific plate, and paleomagnetic data, disprove such a change. Speculations that plumes move, jump, etc. do not revive falsified conjecture. Geochemical distinctions between enriched island and depleted ridge basalts (which overlap) are expected products of normal upper-mantle processes, not plumes. MORB traverses solidus-T asthenosphere, whereas OIB zone-refines through subsolidus lithosphere and crust, crystallizing refractories to retain T of diminishing melt while assimilating and retaining fusibles. Tomographic inference of deep-mantle subduction is presented misleadingly and may reflect methodological and sampling artifacts (downward smearing, and concentration of recorded body waves in bundles within broad anomalies otherwise poorly sampled). Planetological and other data require hot Earth accretion, and thorough early fractionation, from material much more refractory than primitive meteorites, and are incompatible with the little-fractionated lower mantle postulated to permit whole-mantle circulation. The profound seismic discontinuity near 660 km is a thermodynamic and physical barrier to easy mass transfer in either direction. Refractory lower mantle convects slowly, perhaps in layers, and loses primarily original heat, whereas upper mantle churns rapidly, and the 660 decoupling boundary must have evolved into a compositional barrier also. Plate motions are driven by subduction, the passive falling away of oceanic lithosphere which is negatively buoyant because of top-down cooling. Slabs have top and bottom rolling hinges and sink subvertically (inclinations of slabs mark their positions, not trajectories) into the transition zone, where they are laid down on, and depress, the 660-km discontinuity. Rollback of upper hinges into subducting plates is required by plate behavior at all scales. That fronts of overriding plates advance at rollback velocity is required by common preservation atop their thin leading edges of little-deformed fore-arc basins. Convergence velocity also commonly equals rollback but is faster in some arcs. Steeply-sinking inclined slabs push sublithospheric upper mantle forward into the shrinking ocean from which they came, forcing seafloor spreading therein, and pull overriding plates behind them. Continental plates pass over sunken slabs like tanks above their basal treads, and material from, and displaced rearward by, sunken slabs is cycled into pull-apart oceans opening behind the continents, thus transferring mantle from shrinking to enlarging oceans. Hot mantle displaced above slabs enables backarc spreading. Spreading ridges, in both shrinking and enlarging oceans, are passive byproducts of subduction, and migrate because it is more energy efficient to process new asthenosphere than to get partial melt from increasingly distant sources. A plate-motion framework wherein hinges roll back, ridges migrate, Antarctica is approximately fixed, and intraplate deformation is integrated may approximate an absolute reference to sluggish lower mantle, whereas the hotspot frame is invalid, and the no-net-rotation frame minimizes trench and ridge motions.

Location of eastern Mediterranean hot springs induced by mantle heat flow due to slab roll-back and tearing

NASA Astrophysics Data System (ADS)

Roche, V. M.; Sternai, P.; Guillou-Frottier, L.; Jolivet, L.; Gerya, T.

2016-12-01

The Aegean-Anatolian retreating subduction and collision zones have been investigated through 3D numerical geodynamic models involving slab rollback/tearing/breakoff constrained by, for instance, seismic tomography or anisotropy and geochemical proxies. Here we integrate these investigations by using geothermal anomaly measurements from western Turkey. Such data provides insights into the thermal state of the Aegean-Anatolian region at depth and reflects the development of a widespread active geothermal province that is unlikely to be related only to the Quaternary volcanism because this has a too limited extent in space and time. Firstly, we look for possible connections with larger-scale mantle dynamics and use 3D high-resolution petrological and thermo-mechanical numerical models to quantify the potential contribution of the Aegean-Anatolian subduction dynamics to such measured thermal anomalies. Secondly, the subduction-induced thermal signature at the base of the continental crust is then inserted as the imposed basal thermal condition of 2D models dedicated to the understanding of fluid flow in the shallow crust. These models couple heat transfer and fluid flow equations with appropriate fluid and rock physical properties. Results from the 3D numerical models suggest an efficient control of subduction-related asthenospheric return flow on the regional distribution of thermal anomalies. Results from the 2D numerical models also highlight that low angle normal faults (detachments) in the back-arc region can control the bulk of the heat transport and fluid circulation patterns. Such detachments can drain hot crustal and/or mantellic fluids down to several kilometers depths, thus allowing for or fostering deep fluid circulation.

Pb sbnd Sr sbnd Nd isotopic data of Indian Ocean ridges: new evidence of large-scale mapping of mantle heterogeneities

NASA Astrophysics Data System (ADS)

Hamelin, Bruno; Dupré, Bernard; Allègre, Claude J.

1986-01-01

A Pb sbnd Sr sbnd Nd isotope study of South West and East Indian Ridges confirms that the Indian Ocean belongs to a specific regional isotopic domain, as previously suggested by the results from islands of this ocean. The isotopic domain defined by the Indian MORB is indeed different from that of the North Atlantic and East Pacific Oceans. This demonstrates that the convective circulation of the upper mantle does not allow a rapid homogenization from one region to the other. The isotopic data of the Indian ridges can be interpreted by a contamination model, in which the depleted upper mantle (identical to that under the North Atlantic) is contaminated by two different types of contaminant, one corresponding to the source of the "central Indian Ocean" islands (Amsterdam, St. Paul, Marion, Prince Edward, Réunion, Rodriguez, Mauritius), and the other to a source similar to that of Walvis or Ninety East aseismic ridges. These two contaminants would have contributed to the ridge volcanism in different proportions over time.

Rheological decoupling at the Moho and implication to Venusian tectonics.

PubMed

Azuma, Shintaro; Katayama, Ikuo; Nakakuki, Tomoeki

2014-03-18

Plate tectonics is largely responsible for material and heat circulation in Earth, but for unknown reasons it does not exist on Venus. The strength of planetary materials is a key control on plate tectonics because physical properties, such as temperature, pressure, stress, and chemical composition, result in strong rheological layering and convection in planetary interiors. Our deformation experiments show that crustal plagioclase is much weaker than mantle olivine at conditions corresponding to the Moho in Venus. Consequently, this strength contrast may produce a mechanical decoupling between the Venusian crust and interior mantle convection. One-dimensional numerical modeling using our experimental data confirms that this large strength contrast at the Moho impedes the surface motion of the Venusian crust and, as such, is an important factor in explaining the absence of plate tectonics on Venus.

Rheological decoupling at the Moho and implication to Venusian tectonics

PubMed Central

Azuma, Shintaro; Katayama, Ikuo; Nakakuki, Tomoeki

2014-01-01

Plate tectonics is largely responsible for material and heat circulation in Earth, but for unknown reasons it does not exist on Venus. The strength of planetary materials is a key control on plate tectonics because physical properties, such as temperature, pressure, stress, and chemical composition, result in strong rheological layering and convection in planetary interiors. Our deformation experiments show that crustal plagioclase is much weaker than mantle olivine at conditions corresponding to the Moho in Venus. Consequently, this strength contrast may produce a mechanical decoupling between the Venusian crust and interior mantle convection. One-dimensional numerical modeling using our experimental data confirms that this large strength contrast at the Moho impedes the surface motion of the Venusian crust and, as such, is an important factor in explaining the absence of plate tectonics on Venus. PMID:24638113

REVIEWS OF TOPICAL PROBLEMS: Free convection in geophysical processes

NASA Astrophysics Data System (ADS)

Alekseev, V. V.; Gusev, A. M.

1983-10-01

A highly significant geophysical process, free convection, is examined. Thermal convection often controls the dynamical behavior in several of the earth's envelopes: the atmosphere, ocean, and mantle. Section 2 sets forth the thermohydrodynamic equations that describe convection in a compressible or incompressible fluid, thermochemical convection, and convection in the presence of thermal diffusion. Section 3 reviews the mechanisms for the origin of the global atmospheric and oceanic circulation. Interlatitudinal convection and jet streams are discussed, as well as monsoon circulation and the mean meridional circulation of ocean waters due to the temperature and salinity gradients. Also described are the hypotheses for convective motion in the mantle and the thermal-wave (moving flame) mechanism for inducing global circulation (the atmospheres of Venus and Mars provide illustrations). Eddy formation by convection in a centrifugal force field is considered. Section 4 deals with medium- and small-scale convective processes, including hurricane systems with phase transitions, cellular cloud structure, and convection penetrating into the ocean, with its stepped vertical temperature and salinity microstructure. Self-oscillatory processes involving convection in fresh-water basins are discussed, including effects due to the anomalous (p,T) relation for water.

Geochemical constraints on the origin of serpentinization of oceanic mantle

NASA Astrophysics Data System (ADS)

Li, Z.; Lee, C. A.

2004-12-01

The lower seismic zone of double seismic zones in subducting oceanic lithosphere is suggested to be a result of serpentine or chlorite dehydration in the lithospheric mantle (Hacker et al., 2003). However, the mechanism by which oceanic lithospheric mantle is serpentinized is unclear. One way is through hydrothermal circulation where the lithospheric mantle represents part of the circuit through which seawater passes and then returns to the ocean. Another way is to inject seawater into the lithospheric mantle through fractures in the overlying crust without having a return path of water to the ocean. The two mechanisms differ in that the former is an open system process whereas the latter is a closed system process in which the mantle serves as a ¡°sponge¡± for water. Identifying the dominant process is important. For example, if the mantle is part of a hydrothermal circulation cell, the interaction of seawater with the mantle will influence the composition of seawater. This also has important implications for the heat flow out of seafloor. On the other hand, if serpentinization occurs by a closed system process, there will be no influence on seawater composition. Previous studies have suggested that serpentinization of ophiolite bodies was an isochemical process, hence closed system, but it was not clear in these studies whether serpentinization occurred in situ in the oceanic lithosphere. To better understand serpentinization processes in the oceanic lithosphere, we investigated a continuous transition zone of relatively unaltered harzburgite to completely serpentinized harzburgite in the Feather River Ophiolite in northern California. These samples are highly enriched in Na, K, Rb, Cs, U, and Sr, which strongly suggests that serpentinization occurred while the oceanic lithosphere was beneath the ocean. All samples (n=19) have Al2O3 contents ranging from 0.6 to 2.5 wt.% and have extremely depleted light rare-earth element abundances, indicating that these samples are cpx-free harzburgites, which have experienced roughly 20 to 35% melt extraction. The degree of serpentinization was quantified using the concentration of magnetite, a by-product of serpentinization. The lack of antigorite suggests that serpentinization occurred at temperatures lower than 300 C. By comparing Cr and Cr/Al systematics to that predicted from theoretical partial melting calculations and empirical relationships in unaltered peridotite xenoliths, it is shown that Cr and Al are immobile. Al content was thus used to determine the composition of the protolith, which allows us to estimate the amount of depletion/enrichment of a given element by processes other than melt depletion. Most of the harzburgites show no evidence for mantle metasomatism as evidenced by extreme depletions in LREE elements. Consistent with previous studies, we find no depletions in Mg, Fe, or Ca. As seawater is undersaturated in Mg-bearing minerals, an open system process would yield progressive depletion of Mg as is seen in abyssal peridotites, which have been weathered by seawater at the bottom of the seafloor (e.g., Snow et al. 1995). Collectively, this suggests that, except for the addition of seawater and its constituents, serpentinization of the Feather River Ophiolite, was a closed system process. By combining these observations with the results of our field mapping project, we suggest that serpentinization of the lithospheric mantle occurs by local introduction of seawater through fractures extending from the crust and into the mantle. We find no evidence that serpentinized zones in oceanic lithospheric mantle represents an extremely deep hydrothermal circulation cell.

Relation of major volcanic center concentration on Venus to global tectonic patterns

NASA Technical Reports Server (NTRS)

Crumpler, L. S.; Head, James W.; Aubele, Jayne C.

1993-01-01

Global analysis of Magellan image data indicates that a major concentration of volcanic centers covering about 40 percent of the surface of Venus occurs between the Beta, Atla, and Themis regions. Associated with this enhanced concentration are geological characteristics commonly interpreted as rifting and mantle upwelling. Interconnected low plains in an annulus around this concentration are characterized by crustal shortening and infrequent volcanic centers that may represent sites of mantle return flow and net downwelling. Together, these observations suggest the existence of relatively simple, large-scale patterns of mantle circulation similar to those associated with concentrations of intraplate volcanism on earth.

Constraints on radial anisotropy in the central Pacific upper mantle from the NoMelt OBS array

NASA Astrophysics Data System (ADS)

Russell, J. B.; Gaherty, J. B.; Lin, P. P.; Zebker, M.

2016-12-01

Observations of seismic anisotropy in ocean basins are important for constraining deformation and melting processes in the upper mantle. The NoMelt OBS array was deployed on relatively pristine, 70-Ma seafloor in the central Pacific with the aim of constraining upper-mantle circulation and the evolution of the lithosphere-asthenosphere system. Azimuthal variations in Rayleigh-wave velocity suggest strong anisotropic fabric both in the lithosphere and deep in the asthenosphere, and we aim to evaluate whether radial anisotropy shows a similar pattern. We use a combination of Love waves from earthquakes (20-100 s) as well as high-frequency ambient noise (5-10 s) to estimate VSH in the upper 300 km beneath the NoMelt array. Waveform fitting of the ambient-noise cross spectra provide phase-velocity estimates that are sensitive to the upper 50 km of the mantle. To constrain structure beneath the lid, we employ an array-based approach to measure Love-wave phase velocities across the array using seven shallow-focus events (< 25 km) with high signal-to-noise ratio and diverse azimuthal coverage. The Love wave phase-velocity measurements suggest strong interference of the first overtone for intermediate periods (20-50 s), while longer periods (>60 s) are mostly dominated by fundamental mode energy. Through forward modeling of Love wave Fréchet kernels, we find an extremely strong nonlinearity in individual mode-branch sensitivity that is dependent on the relative velocity difference between the low-velocity zone (LVZ) and the overlying Pacific lid. For the fundamental mode in the presence of a strong LVZ, intermediate periods (20-50 s) have little sensitivity within the lithospheric mantle with peak sensitivity pushed to the base of the low-velocity zone. This peak sensitivity migrates to much shallower depth as the lid/LVZ contrast is reduced. Therefore, we use a Monte Carlo approach to systematically explore the model space and identify the most robust model features required to minimize phase-velocity misfit of the full multimode Love wave arrivals. The resulting VSH model is combined with the NoMelt VSV model to obtain estimates of radial anisotropy for the top 300km of the central Pacific upper-mantle.

A severe combined immunodeficient-hu in vivo mouse model of human primary mantle cell lymphoma.

PubMed

Wang, Michael; Zhang, Liang; Han, Xiaohong; Yang, Jing; Qian, Jianfei; Hong, Sungyoul; Lin, Pei; Shi, Yuankai; Romaguera, Jorge; Kwak, Larry W; Yi, Qing

2008-04-01

To establish a severe combined immunodeficient (SCID)-hu in vivo mouse model of human primary mantle cell lymphoma (MCL) for the study of the biology and novel therapy of human MCL. Primary MCL cells were isolated from spleen, lymph node, bone marrow aspirates, or peripheral blood of six different patients and injected respectively into human bone chips, which had been s.c. implanted in SCID-hu. Circulating human beta(2)-microglobulin in mouse serum was used to monitor the engraftment and growth of patient's MCL cells. H&E staining and immunohistochemical staining with anti-human CD20 and cyclin D1 antibodies were used to confirm the tumor growth and migration. Increasing levels of circulating human beta(2)-microglobulin in mouse serum indicated that the patient's MCL cells were engrafted successfully into human bone chip of SCID-hu mice. The engraftment and growth of patient's MCL cells were dependent on human bone marrow microenvironment. Immunohistochemical staining with anti-human CD20 and cyclin D1 antibodies confirmed that patient's MCL cells were able to not only survive and propagate in the bone marrow microenvironment of the human fetal bone chips, but also similar to the human disease, migrate to lymph nodes, spleen, bone marrow, and gastrointestinal tract of host mice. Treatment of MCL-bearing SCID-hu mice with atiprimod, a novel antitumor compound against the protection of bone marrow stromal cells, induced tumor regression. This is the first human primary MCL animal model that should be useful for the biological and therapeutic research on MCL.

NoMelt Experiment: High-resolution constraints on Pacific upper mantle fabric inferred from radial and azimuthal anisotropy

NASA Astrophysics Data System (ADS)

Russell, J. B.; Gaherty, J. B.; Lin, P. P.; Lizarralde, D.; Collins, J. A.; Hirth, G.; Evans, R. L.

2017-12-01

Observations of seismic anisotropy in the ocean basins are important for constraining deformation and melting processes in the upper mantle. The NoMelt OBS array was deployed on relatively pristine, 70 Ma seafloor in the central Pacific with the aim of constraining upper mantle circulation and the evolution of the lithosphere-asthenosphere system. Surface-waves traversing the array provide a unique opportunity to estimate a comprehensive set of anisotropic parameters. Azimuthal variations in Rayleigh-wave velocity over a period band of 15-180 s suggest strong anisotropic fabric both in the lithosphere and deep in the asthenosphere. High-frequency ambient noise (4-10 s) provides constraints on average VSV and VSH as well as azimuthal variations in both VS and VP in the upper ˜10 km of the mantle. Our best fitting models require radial anisotropy in the uppermost mantle with VSH > VSV by 3 - 7% and as much as 2% radial anisotropy in the crust. Additionally, we find a strong azimuthal dependence for Rayleigh- and Love-wave velocities, with Rayleigh 2θ fast direction parallel to the fossil spreading direction (FSD) and Love 2θ and 4θ fast directions shifted 90º and 45º from the FSD, respectively. These are some of the first direct observations of the Love 2θ and 4θ azimuthal signal, which allows us to directly invert for anisotropic terms G, B, and E in the uppermost Pacific lithosphere, for the first time. Together, these observations of radial and azimuthal anisotropy provide a comprehensive picture of oceanic mantle fabric and are consistent with horizontal alignment of olivine with the a-axis parallel to fossil spreading and having an orthorhombic or hexagonal symmetry.

How supercontinents and superoceans affect seafloor roughness.

PubMed

Whittaker, Joanne M; Müller, R Dietmar; Roest, Walter R; Wessel, Paul; Smith, Walter H F

2008-12-18

Seafloor roughness varies considerably across the world's ocean basins and is fundamental to controlling the circulation and mixing of heat in the ocean and dissipating eddy kinetic energy. Models derived from analyses of active mid-ocean ridges suggest that ocean floor roughness depends on seafloor spreading rates, with rougher basement forming below a half-spreading rate threshold of 30-35 mm yr(-1) (refs 4, 5), as well as on the local interaction of mid-ocean ridges with mantle plumes or cold-spots. Here we present a global analysis of marine gravity-derived roughness, sediment thickness, seafloor isochrons and palaeo-spreading rates of Cretaceous to Cenozoic ridge flanks. Our analysis reveals that, after eliminating effects related to spreading rate and sediment thickness, residual roughness anomalies of 5-20 mGal remain over large swaths of ocean floor. We found that the roughness as a function of palaeo-spreading directions and isochron orientations indicates that most of the observed excess roughness is not related to spreading obliquity, as this effect is restricted to relatively rare occurrences of very high obliquity angles (>45 degrees ). Cretaceous Atlantic ocean floor, formed over mantle previously overlain by the Pangaea supercontinent, displays anomalously low roughness away from mantle plumes and is independent of spreading rates. We attribute this observation to a sub-Pangaean supercontinental mantle temperature anomaly leading to slightly thicker than normal Late Jurassic and Cretaceous Atlantic crust, reduced brittle fracturing and smoother basement relief. In contrast, ocean crust formed above Pacific superswells, probably reflecting metasomatized lithosphere underlain by mantle at only slightly elevated temperatures, is not associated with basement roughness anomalies. These results highlight a fundamental difference in the nature of large-scale mantle upwellings below supercontinents and superoceans, and their impact on oceanic crustal accretion.

Axial high topography and partial melt in the crust and mantle beneath the western Galápagos Spreading Center

USGS Publications Warehouse

Blacic, Tanya M.; Ito, Garrett; Shah, Anjana K.; Canales, Juan Pablo; Lin, Jian

2008-01-01

The hot spot-influenced western Galápagos Spreading Center (GSC) has an axial topographic high that reaches heights of ∼700 m relative to seafloor depth ∼25 km from the axis. We investigate the cause of the unusual size of the axial high using a model that determines the flexural response to loads resulting from the thermal and magmatic structure of the lithosphere. The thermal structure simulated is appropriate for large amounts of cooling by hydrothermal circulation, which tends to minimize the amount of partial melt needed to explain the axial topography. Nonetheless, results reveal that the large axial high near 92°W requires that either the crust below the magma lens contains >35% partial melt or that 20% melt is present in the lower crust and at least 3% in the mantle within a narrow column (<∼10 km wide) extending to depths of 45–65 km. Because melt fractions >35% in the crust are considered unreasonable, it is likely that much of the axial high region of the GSC is underlain by a narrow region of partially molten mantle of widths approaching those imaged seismically beneath the East Pacific Rise. A narrow zone of mantle upwelling and melting, driven largely by melt buoyancy, is a plausible explanation.

Seismic evidence for water transport out of the mantle transition zone beneath the European Alps

NASA Astrophysics Data System (ADS)

Liu, Zhen; Park, Jeffrey; Karato, Shun-ichiro

2018-01-01

The mantle transition zone has been considered a major water reservoir in the deep Earth. Mass transfer across the transition-zone boundaries may transport water-rich minerals from the transition zone into the water-poor upper or lower mantle. Water release in the mantle surrounding the transition zone could cause dehydration melting and produce seismic low-velocity anomalies if some conditions are met. Therefore, seismic observations of low-velocity layers surrounding the transition zone could provide clues of water circulation at mid-mantle depths. Below the Alpine orogen, a depressed 660-km discontinuity has been imaged clearly using seismic tomography and receiver functions, suggesting downwellings of materials from the transition zone. Multitaper-correlation receiver functions show prominent ∼0.5-1.5% velocity reductions at ∼750-800-km depths, possibly caused by partial melting in the upper part of lower mantle. The gap between the depressed 660-km discontinuity and the low-velocity layers is consistent with metallic iron as a minor phase in the topmost lower mantle reported by laboratory studies. Velocity drops atop the 410-km discontinuity are observed surrounding the Alpine orogeny, suggesting upwelling of water-rich rock from the transition zone in response to the downwelled materials below the orogeny. Our results provide evidence that convective penetration of the mantle transition zone pushes hydrated minerals both upward and downward to add hydrogen to the surrounding mantle.

Tomographic and Geodynamic Constraints on Convection-Induced Mixing in Earth's Deep Mantle

NASA Astrophysics Data System (ADS)

Hafter, D. P.; Forte, A. M.; Bremner, P. M.; Glisovic, P.

2017-12-01

Seismological studies reveal two large low-shear-velocity provinces (LLSVPs) in the lowermost mantle (e.g., Su et al. 1994; Wang & Wen 2007; He & Wen 2012), which may represent accumulations of subducted slabs at the CMB (Tan & Gurnis 2005; Christensen & Hoffman 1994) or primordial material generated in the early differentiation of Earth (e.g. Li et al. 2014). The longevity or stability of these large-scale heterogeneities in the deep mantle depends on the vigor and spatial distribution of the convective circulation, which is in turn dependent on the distribution of mantle buoyancy and viscosity (e.g. Glisovic & Forte 2015). Here we explore the state of convective mixing in the mantle using the ASPECT convection code (Kronbichler et al. 2012). A series of experiments are conducted to consider the geochemical and dynamical contributions of LLSVPs to deep-mantle upwellings and corresponding plume-sourced volcanism. The principal feature of these experiments is the use of particle tracers to track geochemical changes in the LLSVPs and mantle plumes in addition to identifying those parts of the mantle that may remain unmixed. We employ 3-D mantle density anomalies derived from joint inversions of seismic, geodynamic and mineral physics constraints and geodynamically-constrained viscosity distributions (Glisovic et al. 2015) to ensure that the predicted flow fields yield a good match to key geophysical constraints (e.g. heat flow, global gravity anomalies and plate velocities).

Migrating Toward Fully 4-D Geodynamical Models of Asthenospheric Circulation and Melt Production at Mid-Ocean Ridges

NASA Astrophysics Data System (ADS)

van Dam, L.; Kincaid, C. R.; Pockalny, R. A.; Sylvia, R. T.; Hall, P. S.

2017-12-01

Lateral migration of mid-ocean ridge spreading centers is a well-documented phenomenon leading to asymmetric melt production and the surficial expressions thereof. This form of plate motion has been difficult to incorporate into both numerical and analogue geodynamical models, and consequently, current estimates of time-dependent flow, material transport, and melting in the mantle beneath ridges are lacking. To address this, we have designed and built an innovative research apparatus that allows for precise and repeatable simulations of mid-ocean ridge spreading and migration. Three pairs of counter-rotating belts with adjustable lateral orientations are scaled to simulate spreading at, and flow beneath, three 600km wide ridge segments with up to 300km transform offsets. This apparatus is attached to a drive system that allows us to test a full range of axis-parallel to axis-normal migration directions, and is suspended above a reservoir of viscous glucose syrup, a scaled analogue for the upper mantle, and neutrally buoyant tracers. We image plate-driven flow in the syrup with high-resolution digital cameras and use particle image velocimetry methods to obtain information about transport pathlines and flow-induced anisotropy. Suites of experiments are run with and without ridge migration to determine the overall significance of migration on spatial and temporal characteristics of shallow mantle flow. Our experiments cover an expansive parameter space by including various spreading rates, migration speeds and directions, degrees of spreading asymmetry, transform-offset lengths, and upper mantle viscosity conditions. Preliminary results highlight the importance of modeling migratory plate forces. Mantle material exhibits a significant degree of lateral transport, particularly between ridge segments and towards the melt triangle. Magma supply to the melting region is highly complex; parcels of material do not necessarily move along fixed streamlines, rather, they can be perturbed upwards and left behind as spreading centers continue to move laterally. These results emphasize that observations of seismic anisotropy should be interpreted in light of intricate flow pathlines, and that melt transport models should consider different paths for melt relative to the solid matrix.

Preliminary Results From the Serpentinite, Extension and Regional Porosity Experiment Across the Nicaraguan Trench (SERPENT)

NASA Astrophysics Data System (ADS)

Key, K. W.; Constable, S.; Evans, R. L.; Naif, S.; Matsuno, T.; Lizarralde, D.

2010-12-01

Water plays an important role in the volcanic processes occurring at convergent margins, as the release of water from the downgoing slab affects the rheology of the mantle, increases melting by lowering the solidus temperature, and alters the chemistry of arc-lavas. Yet, one of the major uncertainties in terms of fluid inputs into the subduction factory concerns the extent of serpentinization of the oceanic upper mantle and the volumes of water that are being carried into the subduction system through this route. In April 2010 we conducted a large-scale marine electromagnetic experiment along a 300 km profile offshore Nicaragua in a region that shows evidence for substantial fault related fluid circulation in the crust and possibly upper mantle, and high Ba/La ratios and water contents in adjacent onshore volcanics that suggest a strong slab fluid input into the arc-melting. Our project is the largest combined controlled-source electromagnetic (CSEM) and magnetotelluric (MT) data set ever collected on an active subduction zone. During the single 28 day research cruise aboard the R/V Melville we collected 54 stations of broadband marine magnetotelluric (MT) data and deep-towed nearly 800 km of controlled-source electromagnetic (CSEM) data. Robust multiple-station array processing of the MT data yields high quality MT responses from 10 to 20,000 s period. The MT responses are fairly 1D over the abyssal plain, showing the effects of a thin veneer of conductive sediments overlying a resistive lithosphere and a deeper conductive mantle. The responses become strongly 2D on the trench outer rise and exhibit large 3D distortions at the bottom of the trench, likely due to a combination of effects from severe topography and seafloor conductivity variations. Two circular CSEM tows of 30 km radius were measured by special long-wire EM (LEM) sensors on the abyssal plain and the outer rise. The LEM data reveals a distinct pattern of electromagnetic polarization that is characteristic of mantle transverse anisotropy. Since the conductive axis is aligned with the fossil ridge-parallel direction and reactivated normal faults in the trench, we interpret this to be caused by conductive serpentinized mantle penetrating faults. Conventional CSEM data recorded at a broad suite of transmission frequencies along the 300 km long profile and a 50 km along strike profile provide constraints on crustal conductivity variations. The analysis of these data is ongoing and will provide a comprehensive picture of the electrical conductivity structure from the seafloor to the upper mantle, representing the entire input into this part of the Central American subduction system. Since conductivity is highly dependent on thermal structure, crack porosity and the presence of serpentinite, our experiment will provide constraints on the depth of active fluid circulation within the oceanic crust and mantle, the variation of fluid circulation with distance from the trench and hence with the degree of plate bending, and the extent of dewatering of the subducting slab in the shallow portion of the mantle wedge.

New Joint Geochemical-Geophysical Record of Time-Dependent Mantle Convection South of Iceland

NASA Astrophysics Data System (ADS)

Jones, S. M.; Murton, B. J.; Fitton, J. G.; White, N. J.; Scientific Team Of Rv Celtic Explorer Cruise Ce0806

2010-12-01

Upper mantle convection is expected to be time-dependent because the mantle Rayleigh number is super-critical by at least 3 orders of magnitude. The North Atlantic V-Shaped Ridges (VSRs) provide a long, spatially extensive and clear record of unsteady mantle convection over time periods of 0.1 to 10 million years. VSRs are diachronous ridges of thick crust formed at the Mid Atlantic Ridge to the north and south of Iceland. Ever since their discovery, it has been generally agreed that this diachronous geometry results from melting anomalies that propagate outward from Iceland within the asthenosphere. The solid earth cycles that generate the VSRs also appear to modulate Atlantic oceanic circulation, since the VSR record correlates well with the % Northern Component Water proxy. Knowledge of geochemical variability associated with the VSRs is lacking in comparison with available geophysical and oceanographic records. Here we report dredged basalt samples obtained during RV Celtic Explorer cruise CE0806 (April-May 2008). Together with an earlier dredged dataset, we now have a record of basalt geochemistry covering intersections of 2 V-Shaped Ridge/Trough pairs with the Mid Atlantic Ridge, as well as a transect along a VSR crest. (1) Along the Mid Atlantic Ridge, spatial variations in trace element concentrations and ratios correlate very well with VSR records from bathymetry and gravity. Minima in incompatible trace element concentrations and in ratios such as La/Sm and Zr/Y correlate with maxima in crustal thickness. There is no associated systematic variation in Nd or Sr isotopes. Observations suggest that VSR melting anomalies are caused by variations in temperature rather than composition of the mantle source. (2) VSR melting anomalies propagate out to at least 1400 km from Iceland, well into the region where fracture zones obscure the topographic and gravitational expressions of the VSRs. A value of >1400 km for the radius of plume head measured at the mid-ocean ridge is >600 km more than the value commonly used in estimating plume volume flux. Both the larger plume head diameter and high propagation speeds implied by VSR planform geometry suggest that some well-known estimates of Icelandic plume flux are too small by up to an order of magnitude. (3) Both geochemical and crustal thickness measurements of VSRs can be matched using a time-dependent mid-ocean ridge melting model atop a simple kinematic model of hotter/cooler pulses travelling radially outward from Iceland within the asthenosphere. (4) A thermal origin for the VSRs strengthens the hypothesis of a direct link between plume pulsing and the fluctuating strength of meridional overturning oceanic circulation. Hotter mantle has a lower density, so locations of higher melt production correspond with locations of increased dynamic support at the mid-ocean ridge. When the same hotter pulses travel beneath the Greenland-Scotland Ridge they cause temporary uplift that affects flow of water between the North Atlantic and Nordic Seas.

Crust-mantle density distribution in the eastern Qinghai-Tibet Plateau revealed by satellite-derived gravity gradients

NASA Astrophysics Data System (ADS)

LI, Honglei; Fang, Jian; Braitenberg, Carla; Wang, Xinsheng

2015-04-01

As the highest, largest and most active plateau on Earth, the Qinghai-Tibet Plateau has a complex crust-mantle structure, especially in its eastern part. In response to the subduction of the lithospheric mantle of the Indian plate, large-scale crustal motion occurs in this area. Despite the many previous studies, geodynamic processes at depth remain unclear. Knowledge of crust and upper mantle density distribution allows a better definition of the deeper geological structure and thus provides critically needed information for understanding of the underlying geodynamic processes. With an unprecedented precision of 1-2 mGal and a spatial resolution better than 100 km, GOCE (Gravity field and steady-state Ocean Circulation Explorer) mission products can be used to constrain the crust-mantle density distribution. Here we used GOCE gravitational gradients at an altitude of 10km after reducing the effects of terrain, sediment thickness variations, and Moho undulations to image the density structures of eastern Tibet up to 200 km depths. We inverted the residual satellite gravitational gradients using a least square approach. The initial density model for the inversion is based on seismic velocities from the tomography. The model is composed of rectangular blocks, having a uniform density, with widths of about 100 km and variable thickness and depths. The thickness of the rectangular cells changes from10 to 60km in accordance with the seismic model. Our results reveal some large-scale, structurally controlled density variations at depths. The lithospheric root defined by higher-density contrast features from southwest to northeast, with shallowing in the central part: base of lithosphere reaches a depth of180 km, less than 100km, and 200 km underneath the Lhasa, Songpan-Ganzi, and Ordos crustal blocks, respectively. However, these depth values only represent a first-order parameterization because they depend on model discretization inherited from the original seismic tomography model. For example, the thickness of the uniform density blocks centered at140 km depth is as large as 60 km. Low-density crustal anomalies beneath the southern Lhasa and Songpan-Ganzi blocks in our model support the idea of weak lower crust and possible crustal flow, as a result of the thermal anomalies caused by the upwelling of hot deep materials. The weak lower crust may cause the decoupling of the upper crust and the mantle. These results are consistent with many other geophysical studies, confirming the effectiveness of the GOCE gravitational gradient data. Using these data in combination with other geodynamic constraints (e.g., gravity and seismic structure and preliminary reference Earth model), an improved dynamic model can be derived.

Seismic Velocity Structure of the Pacific Upper Mantle in the NoMelt Region from Finite-Frequency Traveltime Tomography

NASA Astrophysics Data System (ADS)

Hung, S. H.; Lin, P. Y.; Gaherty, J. B.; Russell, J. B.; Jin, G.; Collins, J. A.; Lizarralde, D.; Evans, R. L.; Hirth, G.

2017-12-01

Surface wave dispersion and magnetotelluric survey from the NoMelt Experiment conducted on 70 Ma central Pacific seafloor revealed an electrically resistive, high shear wave velocity lid of 80 km thick underlain by a non-highly conductive, low-velocity layer [Sarafian et al., 2015; Lin et al., 2016]. The vertical structure of the upper mantle consistent with these observational constraints suggests a plausible convection scenario, where the seismically fast, dehydrated lithosphere preserving very strong fossil spreading fabric moves at a constant plate speed over the hydrated, melt-free athenospheric mantle with the presence of either pressure-driven return flow or thermally-driven small scale circulation. To explore 3-D variations in compressional shear wave velocities related to the lithospheric and asthenospheric mantle dynamics, we employ a multichannel cross correlation method to measure relative traveltime residuals based on the vertical P and traverse S waveforms filtered at 10-33 s from telseismic earthquakes at epicentral distance between 30 and 98 degrees. The obtained P and S residuals show on average peak-to-peak variations of ±0.5 s and ±1 s, respectively, across the NoMelt OBS array. Particularly, the P residuals for most of the events display an asymmetrical pattern with respect to an axis oriented nearly N-S to NE-SW through the array. Preliminary ray-based P tomography results reveal similar asymmetric variations in the uppermost 100 km mantle. To verify the resulting structural features, we will further perform both the P and S traveltime tomography and resolution tests based on a multiscale finite-frequency approach which properly takes into account both the 3D off-path sensitivities of the measured residuals and data-adaptive resolution of the model.

Three-Dimensional Seismic Structure of the Mid-Atlantic Ridge: An Investigation of Tectonic, Magmatic, and Hydrothermal Processes in the Rainbow Area

NASA Astrophysics Data System (ADS)

Dunn, Robert A.; Arai, Ryuta; Eason, Deborah E.; Canales, J. Pablo; Sohn, Robert A.

2017-12-01

To test models of tectonic, magmatic, and hydrothermal processes along slow-spreading mid-ocean ridges, we analyzed seismic refraction data from the Mid-Atlantic Ridge INtegrated Experiments at Rainbow (MARINER) seismic and geophysical mapping experiment. Centered at the Rainbow area of the Mid-Atlantic Ridge (36°14'N), this study examines a section of ridge with volcanically active segments and a relatively amagmatic ridge offset that hosts the ultramafic Rainbow massif and its high-temperature hydrothermal vent field. Tomographic images of the crust and upper mantle show segment-scale variations in crustal structure, thickness, and the crust-mantle transition, which forms a vertical gradient rather than a sharp boundary. There is little definitive evidence for large regions of sustained high temperatures and melt in the lower crust or upper mantle along the ridge axes, suggesting that melts rising from the mantle intrude as small intermittent magma bodies at crustal and subcrustal levels. The images reveal large rotated crustal blocks, which extend to mantle depths in some places, corresponding to off-axis normal fault locations. Low velocities cap the Rainbow massif, suggesting an extensive near-surface alteration zone due to low-temperature fluid-rock reactions. Within the interior of the massif, seismic images suggest a mixture of peridotite and gabbroic intrusions, with little serpentinization. Here diffuse microearthquake activity indicates a brittle deformation regime supporting a broad network of cracks. Beneath the Rainbow hydrothermal vent field, fluid circulation is largely driven by the heat of small cooling melt bodies intruded into the base of the massif and channeled by the crack network and shallow faults.

Quantitative assessments of mantle flow models against seismic observations: Influence of uncertainties in mineralogical parameters

NASA Astrophysics Data System (ADS)

Schuberth, Bernhard S. A.

2017-04-01

One of the major challenges in studies of Earth's deep mantle is to bridge the gap between geophysical hypotheses and observations. The biggest dataset available to investigate the nature of mantle flow are recordings of seismic waveforms. On the other hand, numerical models of mantle convection can be simulated on a routine basis nowadays for earth-like parameters, and modern thermodynamic mineralogical models allow us to translate the predicted temperature field to seismic structures. The great benefit of the mineralogical models is that they provide the full non-linear relation between temperature and seismic velocities and thus ensure a consistent conversion in terms of magnitudes. This opens the possibility for quantitative assessments of the theoretical predictions. The often-adopted comparison between geodynamic and seismic models is unsuitable in this respect owing to the effects of damping, limited resolving power and non-uniqueness inherent to tomographic inversions. The most relevant issue, however, is related to wavefield effects that reduce the magnitude of seismic signals (e.g., traveltimes of waves), a phenomenon called wavefront healing. Over the past couple of years, we have developed an approach that takes the next step towards a quantitative assessment of geodynamic models and that enables us to test the underlying geophysical hypotheses directly against seismic observations. It is based solely on forward modelling and warrants a physically correct treatment of the seismic wave equation without theoretical approximations. Fully synthetic 3-D seismic wavefields are computed using a spectral element method for 3-D seismic structures derived from mantle flow models. This way, synthetic seismograms are generated independent of any seismic observations. Furthermore, through the wavefield simulations, it is possible to relate the magnitude of lateral temperature variations in the dynamic flow simulations directly to body-wave traveltime residuals. The synthetic traveltime data can then be compared - on statistical grounds - to the traveltime variations observed on Earth. Here, we now investigate the influence of uncertainties in the various input parameters that enter our modelling. This is especially important for the material properties at high pressure and high temperature entering the mineralogical models. In particular, this concerns uncertainties that arise from relating measurements in the laboratory to Earth properties on a global scale. As one example, we will address the question on the influence of anelasticity on the variance of global synthetic traveltime residuals. Owing to the differences in seismic frequency content between laboratory measurements (MHz to GHz) and the Earth (mHz to Hz), the seismic velocities given in the mineralogical models need to be adjusted; that is, corrected for dispersion due to anelastic effects. This correction will increase the sensitivity of the seismic velocities to temperature variations. The magnitude of this increase depends on absolute temperature, frequency, the frequency dependence of attenuation and the activation enthalpy of the dissipative process. Especially the latter two are poorly known for mantle minerals and our results indicate that variations in activation enthalpy potentially produce the largest differences in temperature sensitivity with respect to the purely elastic case. We will present new wave propagation simulations and corresponding statistical analyses of traveltime measurements for different synthetic seismic models spanning the possible range of anelastic velocity conversions (while being based on the same mantle circulation model).

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 heterogeneity in the mantle. Our work attempts to reconcile the many proposed current ideas for the deep mantle, giving additional insight from modelling on the latest observations from other Deep Earth disciplines.

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.

Seismic evidence for the depression of the D″ discontinuity beneath the Caribbean: Implication for slab heating from the Earth's core

NASA Astrophysics Data System (ADS)

Ko, Justin Yen-Ting; Hung, Shu-Huei; Kuo, Ban-Yuan; Zhao, Li

2017-06-01

The lowermost 100-300 km of the Earth's mantle commonly regarded as the thermal boundary layer (TBL) of mantle circulation is characterized by its complex physical properties. Beneath the Caribbean this so-called D″ layer features relatively high velocities and abrupt impedance increase at the top (designated as the D″ discontinuity). These seismic characteristics have been attributed to the accumulation of ancient subducted slab material and the phase transition in the major lower mantle mineral of pervoskite. Geodynamic models predict that the blanketing cold slabs may trap enough heat from core to be buoyantly destabilized, and eventually broken apart and entrained into the bottom of the convection cell. Here we explore the D″ structure with unprecedented resolution through modeling traveltimes, amplitudes, and waveform shapes from the USArray. We find an east-to-west asymmetrical undulation of the D″ discontinuity with a V-shaped depression of ∼70-160 km over a lateral distance of 600 km beneath northern South America. The shear velocity perturbations vary in the same trend showing the most pronounced reduction of ∼3-4% below the thinnest D″ layer in close proximity to an intermittently undetected discontinuity. The strong correlation between the D″ topography and velocity variations indicates the phase transition boundary has been perturbed or even disrupted by the large lateral temperature gradient of slab material which has been reheated from the core over extended periods of time.

The survival of geochemical mantle heterogeneities

NASA Astrophysics Data System (ADS)

Albarede, F.

2004-12-01

The last decade witnessed major changes in our perception of the geochemical dynamics of the mantle. Data bases such as PETDB and GEOROC now provide highly constrained estimates of the geochemical properties of dominant rock types and of their statistics, while the new generation of ICP mass spectrometers triggered a quantum leap in the production of high-precision isotopic and elemental data. Such new advances offer a fresh view of mantle heterogeneities and their survival through convective mixing. A vivid example is provided by the new high-density coverage of the Mid-Atlantic ridge by nearly 500 Pb, Nd, and Hf isotopic data. This new data set demonstrates a rich harmonic structure which illustrates the continuing stretching and refolding of subducted plates by mantle convection. Just as for oceanic chemical variability, the survival of mantle geochemical heterogeneities though mantle circulation can be seen as a competition between stirring and renewal. The modern residence (renewal) times of the incompatible lithophile elements in the mantle calculated using data bases vary within a rather narrow range (4-9 Gy). The mantle is therefore not currently at geochemical steady-state and the effect of its primordial layering on modern mantle geochemistry is still strong. Up to 50 percent of incompatible lithophile elements may never have been extracted into the oceanic crust, which generalizes a conclusion reached previously for 40Ar. A balance between the buoyancy flux and viscous dissipation provides frame-independent estimates of the rates of mixing by mantle convection: primordial geochemical anomalies with initial length scales comparable to mantle depths of plate lengths are only marginally visible at the scale of mantle melting underneath mid-ocean ridges (≈~50~km). They may show up, however, in hot spot basalts and even more in melt inclusions. Up to 50 percent primordial material may be present in the mantle, but scattered throughout as small (<~10~km) domains, strongly sheared and refolded, and interlayered with younger recycled material. The exploration of the fine-scale geochemical structure of the mantle and the quest for preserved remnants of very old mantle arise as the strongest priorities of deep Earth geochemistry.

Temporal and Spatial Evolution of Dynamic Support From River Profiles: A Framework for Madagascar and Africa

NASA Astrophysics Data System (ADS)

Paul, J. D.; Roberts, G. G.; White, N.

2012-04-01

It is generally accepted that the surface topography of Africa is a manifestation of convective circulation in the sub-lithospheric mantle. Here, we present an inverse method whereby longitudinal river profiles are interrogated to extract quantitative estimates of spatial and temporal variations in the rate of tectonic uplift. Surface processes can provide an important window into transient convective circulation in the sub-lithospheric mantle. River profiles act as 'tectonic tape recorders': we assume the generation of broad, convex-upward knickzones to represent the effect of tectonic uplift shifting the river system into a state of disequilibrium. Profiles evolve through time primarily via the headward retreat of these knickzones. We use a conjugate gradient inverse algorithm to minimise the misfit between observed river profiles - derived from a regional Digital Elevation Model (DEM) - and calculated profiles obtained by varying the uplift rate history. We jointly invert a total of 98 Malagasy and 570 African river profiles to obtain a history of the cumulative tectonic uplift through geological time. We show that Africa has undergone two phases of rapid uplift: first in Eocene times; secondly, since 10 Ma. While the first gave rise to broad, long wavelength topography, the second led to more localised domal swells of high relief. We propose the existence of two wavelengths of dynamic support, reflecting a change in the style of convection in the upper mantle since 50 Ma. Our results correlate strongly with independent geological estimates of uplift across Africa and Madagascar, while our calculated landscape surface following 50 Myr of uplift corresponds closely to a surface fit across present-day drainage divides. Finally we calculate the solid sediment flux delivered to major African deltas as a function of time. This onshore record provides an important indirect constraint on the history of vertical motions at the surface, and agrees well with the offshore flux record, obtained from mapping the thickness of chronostratigraphic sediment packages at the deltas.

Seawater Circulation and Thermal Sink at OCEAN Ridges - FIELD Evidence in Oman Ophiolite

NASA Astrophysics Data System (ADS)

Nicolas, A. A.; Boudier, F. I.; Cathles, L. M.; Buck, W. R.; Celerier, B. P.

2014-12-01

Exceptionally, the lowermost gabbros in the Oman ophiolite are black and totally fresh, except for minute traces of impregnation by seawater fluids at very high temperature (~1000°C). These black gabbros sharply contrast with normal, whitish gabbros altered down to Low-T~500-350°C. These hydrous alterations are ascribed to an unconventional model of seawater circulation and cooling of the permanent magma chambers of fast spreading ocean ridges. In this model, gabbros issued from the magma chamber cross a ~100 m thick thermal boundary layer (TBL) before reaching a narrow, Low-T high permeability channel where the heated return seawater is flowing towards black smokers and the local gabbros are altered. Uprising mantle diapirs in Oman diverge at ~5 km on each side of the palaeo-ridge axis and feed an overlying magma chamber that closes at this distance from axis. Preservation of black gabbros along the Moho implies that the loop of seawater alteration locally does not reach Moho beyond this ~5km distance (otherwise black gabbros would be altered in whitish gabbros). This defines an internal "thermal sink" within ~5 km to the ridge axis. There, the sink is efficiently cooled by the active hydrothermal convection that is ridge transverse. This has been documented near the Galapagos ridge by marine geophysical data, within the same distance. Beyond this critical distance, the cooling system becomes dominantly conductive and ridge-parallel. The TBL and attached return flow channels must be rising into the overcooled, accreted crust. Beyond the thermal sink, the 500°C isotherm rebounds into the crust. It is only after ~ 1My of crustal drift that this isotherm penetrates into the uppermost mantle in a sustained fashion, developing serpentinites at the expense of peridotites.

Formation and modification of chromitites in the mantle

NASA Astrophysics Data System (ADS)

Arai, Shoji; Miura, Makoto

2016-11-01

Podiform chromitites have long supplied us with unrivaled information on various mantle processes, including the peridotite-magma reaction, deep-seated magmatic evolution, and mantle dynamics. The recent discovery of ultrahigh-pressure (UHP) chromitites not only sheds light on a different aspect of podiform chromitites, but also changes our understanding of the whole picture of podiform chromitite genesis. In addition, new evidence was recently presented for hydrothermal modification/formation chromite/chromitite in the mantle, which is a classical but innovative issue. In this context, we present here an urgently needed comprehensive review of podiform chromitites in the upper mantle. Wall-rock control on podiform chromitite genesis demonstrates that the peridotite-magma reaction at the upper mantle condition is an indispensable process. We may need a large system in the mantle, far larger than the size of outcrops or mining areas, to fulfill the Cr budget requirement for podiform chromitite genesis. The peridotite-magma reaction over a large area may form a melt enriched with Na and other incompatible elements, which mixes with a less evolved magma supplied from the depth to create chromite-oversaturated magma. The incompatible-element-rich magma trapped by the chromite mainly precipitates pargasite and aspidolite (Na analogue of phlogopite), which are stable under upper mantle conditions. Moderately depleted harzburgites, which contain chromite with a moderate Cr# (0.4-0.6) and a small amount of clinopyroxene, are the best reactants for the chromitite-forming reaction, and are the best hosts for podiform chromitites. Arc-type chromitites are dominant in ophiolites, but some are of the mid-ocean ridge type; chromitites may be common beneath the ocean floor, although it has not yet been explored for chromitite. The low-pressure (upper mantle) igneous chromitites were conveyed through mantle convection or subduction down to the mantle transition zone to form ultrahigh-pressure chromitites. Some of these reappear at the shallower mantle, and can coexist with newly formed low-pressure igneous chromitites. High-temperature hydrothermal fluids can dissolve and precipitate chromite, and hydrothermal chromitites (chromitites precipitated from aqueous fluids) are possibly formed within the mantle where the circulation of hydrous fluid is available, e.g., at the mantle wedge.

It's the little things that matter most: The role of volatiles in volcanoes and their magmatic roots

NASA Astrophysics Data System (ADS)

Keller, T.; Suckale, J.

2017-12-01

Many volcanic eruptions are driven by volatiles - mostly H2O and CO2 - that degas from magmas rising up beneath the volcano. Gas expands during ascent, thus frequently creating lavas with upward of 50% vesicularity. That is a particularly compelling observation considering that volatiles are only present at concentrations of order 100 ppm in the mantle source. Yet, even at these small concentrations, volatiles significantly lower the peridotite solidus. That leads to the production of reactive volatile-rich melts at depth, which has important consequences for melt transport in the asthenosphere. Thus, volatiles have a pivotal role both at the beginning and the end of the magmatic storyline. A growing amount of observational evidence provides various perspectives on these systems. Volcanic products are characterised increasingly well by geochemical and petrological data. And, volcano monitoring now often provides continuous records of degassing flux and composition. What is missing to better interpret these data are coupled fluid mechanic and thermodynamic models that link melt production and reactive transport in the mantle and crust with degassing-driven volcanic activity at the surface. Such models need to describe the deformation and segregation of multiple material phases (liquids, solids, gases) and track the reactive transport of diverse chemical components (major elements, trace elements, volatiles). I will present progress towards a generalization of existing two-phase model for melt transport in the mantle, extending them to three-phase flows appropriate for magma circulation and degassing in volcanoes. What sets the two environments apart is the presence of a compressible vapor in volcanoes. Also, volcanic degassing may occur by convecting suspensions as well as porous segregation. The model framework we are developing for these processes is based on mixture theory. Uncovering the underlying physics that connects these diverse expressions of magma transport will provide an opportunity to gain deeper insights into magmatic and volcanic phenomena as related rather than separate processes. In time we may thus come to more fully understand how it is that the little things that are mantle volatiles do matter most in volcanoes and their magmatic roots.

Noble gas models of mantle structure and reservoir mass transfer

NASA Astrophysics Data System (ADS)

Harrison, Darrell; Ballentine, Chris J.

Noble gas observations from different mantle samples have provided some of the key observational data used to develop and support the geochemical "layered" mantle model. This model has dominated our conceptual understanding of mantle structure and evolution for the last quarter of a century. Refinement in seismic tomography and numerical models of mantle convection have clearly shown that geochemical layering, at least at the 670 km phase change in the mantle, is no longer tenable. Recent adaptations of the mantle-layering model that more successfully reconcile whole-mantle convection with the simplest data have two common features: (i) the requirement for the noble gases in the convecting mantle to be sourced, or "fluxed", by a deep long-lived volatile-rich mantle reservoir; and (ii) the requirement for the deep mantle reservoirs to be seismically invisible. The fluxing requirement is derived from the low mid-ocean ridge basalt (MORB)-source mantle 3He concentration, in turn calculated from the present day 3He flux from mid-ocean ridges into the oceans (T½ ˜ 1,000 yr) and the ocean crust generation rate (T½ ˜ 108 yr). Because of these very different residence times we consider the 3He concentration constraint to be weak. Furthermore, data show 3He/22Ne ratios derived from different mantle reservoirs to be distinct and require additional complexities to be added to any model advocating fluxing of the convecting mantle from a volatile-rich mantle reservoir. Recent work also shows that the convecting mantle 20Ne/22Ne isotopic composition is derived from an implanted meteoritic source and is distinct from at least one plume source system. If Ne isotope heterogeneity between convecting mantle and plume source mantle is confirmed, this result then excludes all mantle fluxing models. While isotopic heterogeneity requires further quantification, it has been shown that higher 3He concentrations in the convecting mantle, by a factor of 3.5, remove the need for the noble gases in the convecting mantle to be sourced from such a deep hidden reservoir. This "zero paradox" concentration [Ballentine et al., 2002] is then consistent with the different mantle source 3He/22Ne and 20Ne/22Ne heterogeneities. Higher convecting mantle noble gas concentrations also eliminate the requirement for a hidden mantle 40Ar rich-reservoir and enables the heat/4He imbalance to be explained by temporal variance in the different mechanisms of heat vs. He removal from the mantle system—two other key arguments for mantle layering. Confirmation of higher average convecting mantle noble gas concentrations remains the key test of such a concept.

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.

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 dynamics to assess the influence of realistic mineralogical density on mantle convection.

Dynamics of Compressible Convection and Thermochemical Mantle Convection

NASA Astrophysics Data System (ADS)

Liu, Xi

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

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.

Water pumping in mantle shear zones

PubMed Central

Précigout, Jacques; Prigent, Cécile; Palasse, Laurie; Pochon, Anthony

2017-01-01

Water plays an important role in geological processes. Providing constraints on what may influence the distribution of aqueous fluids is thus crucial to understanding how water impacts Earth's geodynamics. Here we demonstrate that ductile flow exerts a dynamic control on water-rich fluid circulation in mantle shear zones. Based on amphibole distribution and using dislocation slip-systems as a proxy for syn-tectonic water content in olivine, we highlight fluid accumulation around fine-grained layers dominated by grain-size-sensitive creep. This fluid aggregation correlates with dislocation creep-accommodated strain that localizes in water-rich layers. We also give evidence of cracking induced by fluid pressure where the highest amount of water is expected. These results emphasize long-term fluid pumping attributed to creep cavitation and associated phase nucleation during grain size reduction. Considering the ubiquitous process of grain size reduction during strain localization, our findings shed light on multiple fluid reservoirs in the crust and mantle. PMID:28593947

Crustal accretion at fast spreading ridges and implications for hydrothermal circulation

NASA Astrophysics Data System (ADS)

Theissen-Krah, S.; Rupke, L.; Hasenclever, J.

2015-12-01

Oceanic crust is continuously created at mid-ocean ridges, but the location of lower crust crystallization continues to be debated since the proposal of the gabbro glacier and many sills end-member models. Geophysical and geochemical studies find evidence for either of the models. The crust is cooled by a combination of heat diffusion and advection, and hydrothermal circulation is thought to play a key role in distinguishing between both models. We use our numerical model for joint modeling of crustal accretion and hydrothermal circulation1 to test different accretion and hydrothermal cooling scenarios. The results match the seismic and structural observations from the East Pacific Rise2 and the Oman Ophiolite3, with a shallow melt lens at the correct location overlaying a narrow volume of partially molten rocks. Our results show that no more than 25-50% of the lower crust crystallizes in situ and that deep circulation is likely to occur at fast and intermediate spreading ridges. The occurrence of deep hydrothermal cooling however does not rule out that a major portion of the lower crust is formed in the shallow melt lens; our simulations rather suggest that it is necessary independent of where in the lower crust crystallization takes place. 1 Theissen-Krah, S., Iyer, K., Rupke, L. H. & Morgan, J. P. Coupled mechanical and hydrothermal modeling of crustal accretion at intermediate to fast spreading ridges. Earth and Planetary Science Letters 311, 275-286, doi:10.1016/j.epsl.2011.09.018 (2011). 2 Dunn, R. A., Toomey, D. R. & Solomon, S. C. Three-dimensional seismic structure and physical properties of the crust and shallow mantle beneath the East Pacific Rise at 9 degrees 30'N. Journal of Geophysical Research-Solid Earth 105, 23537-23555 (2000). 3 Nicolas, A. & Boudier, F. Structural contribution from the Oman ophiolite to processes of crustal accretion at the East Pacific Rise. Terra Nova 27, 77-96, doi:10.1111/ter.12137 (2015).

The Upper Mantle Flow Field around South-Africa as Reflected by Isotopic Provinciality

NASA Astrophysics Data System (ADS)

Meyzen, C.; Blichert-Toft, J.; Ludden, J.; Humler, E.; Mevel, C.; Albarede, F.

2006-12-01

Isotopic studies of MORB have established the existence of broad isotopic provinces within the underlying asthenosphere, such as in the Indian Ocean (DUPAL). How these features relate to mantle circulation is, however, still unknown. The steepness of the transition between such isotopic provinces will define the geometry of the velocity field in the upper mantle. In this respect, the transition between the Indian and South Atlantic provinces, two domains that are isotopically contrasted, should be readily identifiable over this long ridge segment. Here, we present Hf isotope data for 60 samples dredged along the SWIR between 35° and 69°E. The new Hf isotope data show that the Indian asthenosphere does not spill directly into the South Atlantic upper mantle: the general decreasing southward gradient observed for ^{176}Hf/^{177}Hf down the mid- Atlantic Ridge, and also for Sr isotopes and model Th/U ratios (derived from Pb isotopes), is overprinted by material with radiogenic Sr, unradiogenic Hf and high Th/U. The Indian domain grades into the South Atlantic around Bouvet, while the South Atlantic collides with the Atlantic province around Tristan. We interpret these features to represent fronts between three adjacent isotopic provinces similar to what has been suggested for the Australian-Antarctic Discordance. The common DUPAL signature of MORB and OIB from the Indian province and the geochemistry of Gulf of Aden MORB and the Afar plume suggest that the source of this distinctive mantle component is deep and lies to the north of the province. This is also what the three-dimensional flow field computed by Behn et al. (2004) from shear-wave splitting shows with a major lower mantle upwelling radiating at the base of the asthenosphere under the Afar plume. Lower mantle gushing out from this source flows southward unimpeded along the Indian ridges, whereas it only reaches the South Atlantic ridge after first having been deflected under the deep roots of the South African Archean cratons. Erosion of these roots by the asthenospheric drift confers a distinct continental signature on the source of South Atlantic MORB. This pattern is also consistent with the observation that the lowest He isotope values occur, on average, along the South Atlantic ridge. To some extent, the dynamics of the North Atlantic upper mantle mirrors the Indian situation: the flow field of Behn et al. (2004) shows that the North Atlantic asthenosphere also fills up through deep mantle upwellings, which is consistent with the Dupal-like isotopic signature of the Arctic ridges. M.D. Behn, C.P. Conrad and P.G. Silver (2004), Detection of upper mantle flow associated with the African Superplume, Earth. Planet. Sci. Lett., 224, 259-274.

Modes of mantle convection and the removal of heat from the earth's interior

NASA Technical Reports Server (NTRS)

Spohn, T.; Schubert, G.

1982-01-01

Thermal histories for two-layer and whole-mantle convection models are calculated and presented, based on a parameterization of convective heat transport. The model is composed of two concentric spherical shells surrounding a spherical core. The models were constrained to yield the observed present-day surface heat flow and mantle viscosity, in order to determine parameters. These parameters were varied to determine their effects on the results. Studies show that whole-mantle convection removes three times more primordial heat from the earth interior and six times more from the core than does two-layer convection (in 4.5 billion years). Mantle volumetric heat generation rates for both models are comparable to that of a potassium-depleted chondrite, and thus surface heat-flux balance does not require potassium in the core. Whole and two-layer mantle convection differences are primarily due to lower mantle thermal insulation and the lower heat removal efficiency of the upper mantle as compared with that of the whole mantle.

A mantle plume model for the Equatorial Highlands of Venus

NASA Technical Reports Server (NTRS)

Kiefer, Walter S.; Hager, Bradford H.

1991-01-01

The possibility that the Equatorial Highlands are the surface expressions of hot upwelling mantle plumes is considered via a series of mantle plume models developed using a cylindrical axisymmetric finite element code and depth-dependent Newtonian rheology. The results are scaled by assuming whole mantle convection and that Venus and the earth have similar mantle heat flows. The best model fits are for Beta and Atla. The common feature of the allowed viscosity models is that they lack a pronounced low-viscosity zone in the upper mantle. The shape of Venus's long-wavelength admittance spectrum and the slope of its geoid spectrum are also consistent with the lack of a low-viscosity zone. It is argued that the lack of an asthenosphere on Venus is due to the mantle of Venus being drier than the earth's mantle. Mantle plumes may also have contributed to the formation of some smaller highland swells, such as the Bell and Eistla regions and the Hathor/Innini/Ushas region.

Experimental constraints on the serpentinization rate of fore-arc peridotites: Implications for the upwelling condition of the slab-derived fluid

NASA Astrophysics Data System (ADS)

Nakatani, T.; Nakamura, M.

2016-08-01

To constrain the water circulation in subduction zones, the hydration rates of peridotites were investigated experimentally in fore-arc mantle conditions. Experiments were conducted at 400-580°C and 1.3 and 1.8 GPa, where antigorite is expected to form as a stable serpentine phase. Crushed powders of olivine ± orthopyroxene and orthopyroxene + clinopyroxene were reacted with 15 wt % distilled water for 4-19 days. The synthesized serpentine varieties were lizardite and aluminous lizardite (Al-lizardite) in all experimental conditions except those of 1.8 GPa and 580°C in the olivine + orthopyroxene system, in which antigorite was formed. In the olivine + orthopyroxene system, the reactions were interface-controlled except for the reaction at 400°C, which was transport-controlled. The corresponding reaction rates were 7.0 × 10-12 to 1.5 × 10-11 m s-1 at 500-580°C and 7.5 × 10-16 m2 s-1 at 400°C for the interface and transport-controlled reactions, respectively. Based on a simple reaction-transport model including these hydration rates, we infer that penetration of the slab-derived fluid all the way through a water-unsaturated fore-arc mantle is allowed only when focused flow occurs with a spacing larger than 77-229 km in hot subduction zones such as Nankai and Cascadia. However, the necessary spacing is only 2.3-4.6 m in intermediate-temperature subduction zones such as Kyushu and Costa Rica. These calculations imply that fluid leakage in hot subduction zones may occur after the fore-arc mantle is totally hydrated, whereas in intermediate-temperature subduction zones, leakage through a water-unsaturated fore-arc mantle may be facilitated.

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.

Tomography and Dynamics of Western-Pacific Subduction Zones

NASA Astrophysics Data System (ADS)

Zhao, D.

2012-01-01

We review the significant recent results of multiscale seismic tomography of the Western-Pacific subduction zones and discuss their implications for seismotectonics, magmatism, and subduction dynamics, with an emphasis on the Japan Islands. Many important new findings are obtained due to technical advances in tomography, such as the handling of complex-shaped velocity discontinuities, the use of various later phases, the joint inversion of local and teleseismic data, tomographic imaging outside a seismic network, and P-wave anisotropy tomography. Prominent low-velocity (low-V) and high-attenuation (low-Q) zones are revealed in the crust and uppermost mantle beneath active arc and back-arc volcanoes and they extend to the deeper portion of the mantle wedge, indicating that the low-V/low-Q zones form the sources of arc magmatism and volcanism, and the arc magmatic system is related to deep processes such as convective circulation in the mantle wedge and dehydration reactions in the subducting slab. Seismic anisotropy seems to exist in all portions of the Northeast Japan subduction zone, including the upper and lower crust, the mantle wedge and the subducting Pacific slab. Multilayer anisotropies with different orientations may have caused the apparently weak shear-wave splitting observed so far, whereas recent results show a greater effect of crustal anisotropy than previously thought. Deep subduction of the Philippine Sea slab and deep dehydration of the Pacific slab are revealed beneath Southwest Japan. Significant structural heterogeneities are imaged in the source areas of large earthquakes in the crust, subducting slab and interplate megathrust zone, which may reflect fluids and/or magma originating from slab dehydration that affected the rupture nucleation of large earthquakes. These results suggest that large earthquakes do not strike anywhere, but in only anomalous areas that may be detected with geophysical methods. The occurrence of deep earthquakes under the Japan Sea and the East Asia margin may be related to a metastable olivine wedge in the subducting Pacific slab. The Pacific slab becomes stagnant in the mantle transition zone under East Asia, and a big mantle wedge (BMW) has formed above the stagnant slab. Convective circulations and fluid and magmatic processes in the BMW may have caused intraplate volcanism (e.g., Changbai and Wudalianchi), reactivation of the North China craton, large earthquakes, and other active tectonics in East Asia. Deep subduction and dehydration of continental plates (such as the Eurasian plate, Indian plate and Burma microplate) are also found, which have caused intraplate magmatism (e.g., Tengchong) and geothermal anomalies above the subducted continental plates. Under Kamchatka, the subducting Pacific slab shortens toward the north and terminates near the Aleutian-Kamchatka junction. The slab loss was induced by friction with the surrounding asthenosphere, as the Pacific plate rotated clockwise 30 Ma ago, and then it was enlarged by the slab-edge pinch-off by the asthenospheric flow. The stagnant slab finally collapses down to the bottom of the mantle, which may trigger upwelling of hot mantle materials from the lower mantle to the shallow mantle. Suggestions are also made for future directions of the seismological research of subduction zones.

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-3 and degrees 4-9 geoid, either a gradually increase of viscosity in the upper part of lower mantle or larger thermal expansivity in the lower mantle is needed. We also consider thermo-chemical models to examine the effects of the African and Pacific thermochemical piles (i.e., LLSVPSs) on the geoid and the inferred mantle viscosity and Ra.

A kinematic model for the late Cenozoic development of southern California crust and upper mantle

NASA Technical Reports Server (NTRS)

Humphreys, Eugene D.; Hager, Bradford H.

1990-01-01

A model is developed for the young and ongoing kinematic deformation of the southern California crust and upper mantle. The kinematic model qualitatively explains both the overall seismic structure of the upper mantle and much of the known geological history of the late Cenozoic as consequences of ongoing convection beneath southern California. In this model, the high-velocity upper-mantle anomaly of the Transverse ranges is created through the convergence and sinking of the entire thickness of subcrustal lihtosphere, and the low-velocity upper-mantle anomaly beneath the Salton Trough region is attributed to high temperatures and 1-4 percent partial melt related to adiabatic decompression during mantle upwelling.

Compositional mantle layering revealed by slab stagnation at ~1000-km depth

PubMed Central

Ballmer, Maxim D.; Schmerr, Nicholas C.; Nakagawa, Takashi; Ritsema, Jeroen

2015-01-01

Improved constraints on lower-mantle composition are fundamental to understand the accretion, differentiation, and thermochemical evolution of our planet. Cosmochemical arguments indicate that lower-mantle rocks may be enriched in Si relative to upper-mantle pyrolite, whereas seismic tomography images suggest whole-mantle convection and hence appear to imply efficient mantle mixing. This study reconciles cosmochemical and geophysical constraints using the stagnation of some slab segments at ~1000-km depth as the key observation. Through numerical modeling of subduction, we show that lower-mantle enrichment in intrinsically dense basaltic lithologies can render slabs neutrally buoyant in the uppermost lower mantle. Slab stagnation (at depths of ~660 and ~1000 km) and unimpeded slab sinking to great depths can coexist if the basalt fraction is ~8% higher in the lower mantle than in the upper mantle, equivalent to a lower-mantle Mg/Si of ~1.18. Global-scale geodynamic models demonstrate that such a moderate compositional gradient across the mantle can persist can in the presence of whole-mantle convection. PMID:26824060

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

PubMed Central

Lee, Changyeol; Wada, Ikuko

2017-01-01

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

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

PubMed

Lee, Changyeol; Wada, Ikuko

2017-06-29

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

Electrical conductivity imaging in the western Pacific subduction zone

NASA Astrophysics Data System (ADS)

Utada, Hisashi; Baba, Kiyoshi; Shimizu, Hisayoshi

2010-05-01

Oceanic plate subduction is an important process for the dynamics and evolution of the Earth's interior, as it is regarded as a typical downward flow of the mantle convection that transports materials from the near surface to the deep mantle. Recent seismological study showed evidence suggesting the transportation of a certain amount of water by subduction of old oceanic plate such as the Pacific plate down to 150-200 km depth into the back arc mantle. However it is not well clarified how deep into the mantle the water can be transported. The electromagnetic induction method to image electrical conductivity distribution is a possible tool to answer this question as it is known to be sensitive to the presence of water. Here we show recent result of observational study from the western Pacific subduction zone to examine the electrical conductivity distribution in the upper mantle and in the mantle transition zone (MTZ), which will provide implications how water distributes in the mantle. We take two kinds of approach for imaging the mantle conductivity, (a) semi-global and (b) regional induction approaches. Result may be summarized as follows: (a) Long (5-30 years) time series records from 8 submarine cables and 13 geomagnetic observatories in the north Pacific region were analyzed and long period magnetotelluric (MT) and geomagnetic deep sounding (GDS) responses were estimated in the period range from 1.7 to 35 days. These frequency dependent response functions were inverted to 3-dimensional conductivity distribution in the depth range between 350 and 850 km. Three major features are suggested in the MTZ depth such as, (1) a high conductivity anomaly beneath the Philippine Sea, (2) a high conductivity anomaly beneath the Hawaiian Islands, and (3) a low conductivity anomaly beneath and in the vicinity of northern Japan. (b) A three-year long deployment of ocean bottom electro-magnetometers (OBEM's) was conducted in the Philippine Sea and west Pacific Ocean from 2005 to 2008. As a preliminary investigation, MT response functions from 20 sites in the Philippine Sea and 4 sites in the west Pacific basin in the period range between 300 and 80000 sec were respectively inverted to one-dimensional (1-D) profile of electrical conductivity by quantitatively considering the effect of the heterogeneous conductivity distribution (ocean and lands) at the surface. The resultant 1-D models show three main features: (1) Strong contrast in the conductivity for the shallower 200 km of the upper mantle depths is recognized between the two regions, which is qualitatively consistent with the difference in lithospheric age. (2) The conductivity at 200-300 km depth is more or less similar to each other at about 0.3 S /m. (3) The conductivity around the MTZ depth is higher for the Philippine Sea mantle than for the Pacific mantle, which is consistent with the implication obtained from a semi-global approach (a). As already suggested in our previous work, the high conductivity in the MTZ below the Philippine Sea can be explained by the excess conduction due to the presence of hydrogen (water) in wadesleyite or in ringwoodite. Therefore, it implies a large scale circulation of water in the back arc mantle not only in the upper mantle but also down to the MTZ depth. However, our interpretation indicates that the high conductivity of the Philippine Sea uppermost upper mantle cannot be explained only by the effect of hydrogen conduction in olivine, but that additional conduction enhancement such as the presence of partial melt is required.

Origin of geochemical mantle components: Role of spreading ridges and thermal evolution of mantle

NASA Astrophysics Data System (ADS)

Kimura, Jun-Ichi; Gill, James B.; van Keken, Peter E.; Kawabata, Hiroshi; Skora, Susanne

2017-02-01

We explore the element redistribution at mid-ocean ridges (MOR) using a numerical model to evaluate the role of decompression melting of the mantle in Earth's geochemical cycle, with focus on the formation of the depleted mantle component. Our model uses a trace element mass balance based on an internally consistent thermodynamic-petrologic computation to explain the composition of MOR basalt (MORB) and residual peridotite. Model results for MORB-like basalts from 3.5 to 0 Ga indicate a high mantle potential temperature (Tp) of 1650-1500°C during 3.5-1.5 Ga before decreasing gradually to ˜1300°C today. The source mantle composition changed from primitive (PM) to depleted as Tp decreased, but this source mantle is variable with an early depleted reservoir (EDR) mantle periodically present. We examine a two-stage Sr-Nd-Hf-Pb isotopic evolution of mantle residues from melting of PM or EDR at MORs. At high-Tp (3.5-1.5 Ga), the MOR process formed extremely depleted DMM. This coincided with formation of the majority of the continental crust, the subcontinental lithospheric mantle, and the enriched mantle components formed at subduction zones and now found in OIB. During cooler mantle conditions (1.5-0 Ga), the MOR process formed most of the modern ocean basin DMM. Changes in the mode of mantle convection from vigorous deep mantle recharge before ˜1.5 Ga to less vigorous afterward is suggested to explain the thermochemical mantle evolution.

Seismic-geodynamic constraints on three-dimensional structure, vertical flow, and heat transfer in the mantle

USGS Publications Warehouse

Forte, A.M.; Woodward, R.L.

1997-01-01

Joint inversions of seismic and geodynamic data are carried out in which we simultaneously constrain global-scale seismic heterogeneity in the mantle as well as the amplitude of vertical mantle flow across the 670 km seismic discontinuity. These inversions reveal the existence of a family of three-dimensional (3-D) mantle models that satisfy the data while at the same time yielding predictions of layered mantle flow. The new 3-D mantle models we obtain demonstrate that the buoyancy forces due to the undulations of the 670 km phase-change boundary strongly inhibit the vertical flow between the upper and lower mantle. The strong stabilizing effect of the 670 km topography also has an important impact on the predicted dynamic topography of the Earth's solid surface and on the surface gravity anomalies. The new 3-D models that predict strongly or partially layered mantle flow provide essentially identical fits to the global seismic data as previous models that have, until now, predicted only whole-mantle flow. The convective vertical transport of heat across the mantle predicted on the basis of the new 3-D models shows that the heat flow is a minimum at 1000 km depth. This suggests the presence at this depth of a globally defined horizon across which the pattern of lateral heterogeneity changes rapidly. Copyright 1997 by the American Geophysical Union.

Predicting SKS-splitting from 35 Myr of subduction and mantle flow evolution in the western Mediterranean

NASA Astrophysics Data System (ADS)

Chertova, Maria; Spakman, Wim; Faccenda, Manuele

2017-04-01

We investigate the development of mantle anisotropy associated with the evolution of the Rif-Gibraltar-Betic (RGB) slab of the western Mediterranean and predict SKS-splitting directions for comparison with the recent observations compiled in Diaz and Gallart (2014). Our numerical model of slab evolution starts at 35 Ma and builds on our on recent work (Chertova et al., 2014) with the extension of imposing mantle flow velocities on the side boundaries of the model (Chertova et al., 2017). For the calculation of the evolution of finite strain deformation from the mantle flow field and for prediction of SKS-splitting directions we use the modified D-Rex program of Faccenda (2014). We test the predicted splitting observations against present-day shear wave splitting observations for subduction models with open boundary conditions (Chertova, 2014) and for models with various prescribed mantle flow conditions on the model side boundaries. The latter are predicted time-dependent (1 Myr time steps) velocity boundary conditions computed from back-advection of a temperature and density model of the present-day mantle scaled from a global seismic tomography model (Steinberger et al., 2015). These boundary conditions where used recently to demonstrate the relative insensitivity of RGB slab position and overall slab morphology for external mantle flow (Chertova et al., 2017). Using open boundaries only we obtain a poor to moderate fit between predicted and observed splitting directions after 35 Myr of slab and mantle flow evolution. In contrast, a good fit is obtained when imposing the computed mantle flow velocities on the western, southern, and northern boundaries during 35 Myr of model evolution. This successful model combines local slab-driven mantle flow with remotely forced mantle flow. We are in the process to repeat these calculations for shorter periods of mantle flow evolution to determine how much of past mantle flow is implicitly recorded in present-day observation of SKS splitting. In combination with our recent work on the influence of external mantle flow on RGB slab evolution (Chertova et al., 2017) we have also demonstrated that (1) the preferred slab evolution model of Chertova et al. (2014; their "Scenario 1" in which RGB subduction starts at the Baleares margin some 35 Myr ago and then rolls back southward to Africa and next to the W and finally to NW to create the future Rif-Gibraltar-Betics cordillera), is robust with respect to the impact of global mantle flow for the past 35 Myr and that (2) only the combination of global flow with local slab-induced flow leads to mantle anisotropy prediction that consistent with present-day observations of present-day SKS splitting. Steinberger, B., W.Spakman, P.Japsen and T.H.Torsvik (2015), The key role of global solid Earth processes in the late Cenozoic intensification of Greenland glaciation. Terra Nova, 27 Chertova, M.V., W.Spakman, T. Geenen, A.P. van den Berg, D.J.J. van Hinsbergen (2014), Underpinning tectonic reconstructions of the western Mediterranean region with dynamic slab evolution from 3-D numerical modeling. J. Geophys. Res. Solid Earth Chertova, M., W.Spakman and B.Steinberger (2017), Mantle flow influence on subduction evolution, submitted to J. Geophys. Res. Solid Earth Faccenda, M. (2014), Mid mantle seismic anisotropy around subduction zones, Physics of the Earth and Planetary Interiors Diaz, J., and J. Gallart (2014) Seismic anisotropy from the Variscan core of Iberia to the western African Craton: New constraints on upper mantle flow at regional scale. Earth and Planetary Science Letters

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

NASA Astrophysics Data System (ADS)

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

2017-07-01

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