Hello Guest User, You are visiting this website from a computer with an IP address of 172.70.42.41 with the name of '?' since Thu Mar 28, 2024 at 2:44:06 PM PT for approx. 0 minutes now.
The Sheahan Diamond Literature Reference Compilation - Scientific and Media Articles based on Major Keyword - Convection
The Sheahan Diamond Literature Reference Compilation is compiled by Patricia Sheahan who publishes on a monthly basis a list of new scientific articles related to diamonds as well as media coverage and corporate announcements called the Sheahan Diamond Literature Service that is distributed as a free pdf to a list of followers. Pat has kindly agreed to allow her work to be made available as an online digital resource at Kaiser Research Online so that a broader community interested in diamonds and related geology can benefit. The references are for personal use information purposes only; when available a link is provided to an online location where the full article can be accessed or purchased directly. Reproduction of this compilation in part or in whole without permission from the Sheahan Diamond Literature Service is strictly prohibited. Return to Diamond Keyword Index
Sheahan Diamond Literature Reference Compilation - Scientific Articles by Author for all years
Each article reference in the SDLRC is tagged with one or more key words assigned by Pat Sheahan to highlight the main topics of the article. In an effort to make it easier for users to track down articles related to a specific topic, KRO has extracted these key words and developed a list of major key words presented in this Key Word Index to which individual key words used in the article reference have been assigned. In most of the individual Key Word Reports the references are in crhonological order, though in some such as Deposits the order is first by key word and then chronological. Only articles classified as "technical" (mainly scientific journal articles) and "media" (independent media articles) are included in the Key Word Index. References that were added in the most recent monthly update are highlighted in yellow.
Convection is the manner by which heat generated inside the earth's core through radioactive decay causes hot rocks within the mantle (astheosphere) to rise while cooler rock sinks, creating a heat circulation that drives the plates which rest on the mantle (lithosphere). Convection is the engine that drivers plate tectonics and generates magmas which work their way to the surface. Convection is relevant to diamonds because it affects craton construction and eruption of kimberlites.
Three dimensional mantle convection simulations with a low viscosity asthenosphere and the relationship between heat flow and the horizontal length scale
Geophysical Research Letters, Vol. 35, 10, May 28, L10304
Global trench migration velocities and slab migration induced upper mantle volume fluxes: constraints to find an Earth reference frame based on minimizing viscous dissipation.
Earth Science Reviews, Vol. 88, 1-2, May pp. 118-144.
Osmium isotopes in Baffin Island and West Greenland picrites: implications for the 187 Os and 188 Os composition of the convection mantle and nature 3He/4he
Earth and Planetary Interiors, Vol. 278, 3-4, pp. 267-277.
Dynamic role of the rheological contrast between cratonic and oceanic lithospheres in the longevity of cratonic lithosphere: a three dimensional numerical study.
On the stability of thermal stratification of highly compressible fluids with depth dependent physical properties: implications for the mantle convection.
Geophysical Journal International, Vol. 195, 3, pp. 1443-1454.
Geochemistry, Geophysics, Geosystems: G3, Vol. 16, 2, pp. 366-379.
United States, Colorado Plateau
Convection
Abstract: Although volcanism in the southwestern United States has been studied extensively, its origin remains controversial. Various mechanisms such as mantle plumes, upwelling in response to slab sinking, and small-scale convective processes have been proposed, but have not been evaluated within the context of rapidly shearing asthenosphere that is thought to underlie this region. Using geodynamic models that include this shear, we here explore spatiotemporal patterns of mantle melting and volcanism near the Colorado Plateau. We show that the presence of viscosity heterogeneity within an environment of asthenospheric shearing can give rise to decompression melting along the margins of the Colorado Plateau. Our models indicate that eastward shear flow can advect pockets of anomalously low viscosity toward the edges of thickened lithosphere beneath the plateau, where they can induce decompression melting in two ways. First, the arrival of the pockets critically changes the effective viscosity near the plateau to trigger small-scale edge-driven convection. Second, they can excite shear-driven upwelling (SDU), in which horizontal shear flow becomes redirected upward as it is focused within the low-viscosity pocket. We find that a combination of “triggered” edge-driven convection and SDU can explain volcanism along the margins of the Colorado Plateau, its encroachment toward the plateau's southwestern edge, and the association of volcanism with slow seismic anomalies in the asthenosphere. Geographic patterns of intraplate volcanism in regions of vigorous asthenospheric shearing may thus directly mirror viscosity heterogeneity of the sublithospheric mantle.
Geological Society of America Special Paper, No. 514, pp. SPE514-05.
Mantle
Convection
Abstract: The role of decoupling in the low-velocity zone is crucial for understanding plate tectonics and mantle convection. Mantle convection models fail to integrate plate kinematics and thermodynamics of the mantle. In a first gross estimate, we computed at >300 km3/yr the volume of the plates lost along subduction zones. Mass balance predicts that slabs are compensated by broad passive upwellings beneath oceans and continents, passively emerging at oceanic ridges and backarc basins. These may correspond to the broad low-wavespeed regions found in the upper mantle by tomography. However, west-directed slabs enter the mantle more than three times faster (?232 km3/yr) than in the opposite east- or northeast-directed subduction zones (?74 km3/yr). This difference is consistent with the westward drift of the outer shell relative to the underlying mantle, which accounts for the steep dip of west-directed slabs, the asymmetry between flanks of oceanic ridges, and the directions of ridge migration. The larger recycling volumes along west-directed subduction zones imply asymmetric cooling of the underlying mantle and that there is an "easterly" directed component of the upwelling replacement mantle. In this model, mantle convection is tuned by polarized decoupling of the advecting and shearing upper boundary layer. Return mantle flow can result from passive volume balance rather than only by thermal buoyancy-driven upwelling.
Geochemistry, Geophysics, Geosystems: G3, Vol. 16, 11, Nov. pp. 3924-3945.
Mantle
Convection
Abstract: Numerical models show that small-scale convection (SSC) occurring atop a mantle plume is a plausible mechanism to rejuvenate the lithosphere. The triggering of SSC depends on the density contrast and on the rheology of the unstable layer underlying the stagnant upper part of the thermal boundary layer (TBL). Partial melting may change both properties. We analyze, using 2-D numerical simulations, how partial melting influences the dynamics of time-dependent SSC instabilities and the resulting thermo-mechanical rejuvenation of an oceanic plate moving atop of a plume. Our simulations show a complex behavior, with acceleration, no change, or delay of the SSC onset, due to competing effects of the latent heat of partial melting, which cools the plume material, and of the buoyancy increase associated with both melt retention and depletion of residue following melt extraction. The melt-induced viscosity reduction is too localized to affect significantly SSC dynamics. Faster SSC triggering is promoted for low melting degrees (low plume temperature anomalies, thick lithosphere, or fast moving plates), which limit both the temperature reduction due to latent heat of melting and the accumulation of depleted buoyant residue in the upper part of the unstable layer. In contrast, high partial melting degrees lead to a strong temperate decrease due to latent heat of melting and development of a thick depleted layer within the sublithospheric convecting layer, which delay the development of gravitational instabilities. Despite differences in SSC dynamics, the thinning of the lithosphere is not significantly enhanced relatively to simulations that neglect partial melting.
Abstract: The chemical composition of Earth's mantle can tell us how our planet formed and how subsequent mantle dynamics have since homogenized the mantle through convective processes. Most terrestrial rocks have a similar tungsten (W) isotope composition (1), but some rocks that have been dated at 2.8 Ga (billion years old) (2), 3.8 Ga (3), and 3.96 Ga (4) have elevated 182W/184W ratios. This is reported as µ182W, in parts per million (ppm) deviation from the bulk silicate Earth. Until now, the outliers have included only these ancient rock samples with a small µ182W excess (?15 ppm) that can be attributed to the final ?0.5% of Earth's mass that accreted late in its accretion history. On page 809 of this issue, Rizo et al. (5) report W isotope data from young mantle-derived rocks with µ182W excesses of 10 to 48 ppm. This result is spectacular because the range of µ182W values in mantle-derived rocks is larger than can be accommodated by late accretion; the implication is that remnants of Earth's earliest mantle have been preserved over the entirety of Earth's history.
Progress in Earth and Planetary Science, Open access
Mantle
Convection, geodynamics
Abstract: Subduction initiation is a key in understanding the dynamic evolution of the Earth and its fundamental difference to all other rocky planetary bodies in our solar system. Despite recent progress, the question about how a stiff, mostly stagnant planetary lid can break and become part in the global overturn of the mantle is still unresolved. Many mechanisms, externally or internally driven, are proposed in previous studies. Here, we present the results on subduction initiation obtained by dynamically self-consistent, time-dependent numerical modelling of mantle convection. We show that the stress distribution and resulting deformation of the lithosphere are strongly controlled by the top boundary formulation: A free surface enables surface topography and plate bending, increases gravitational sliding of the plates and leads to more realistic, lithosphere-scale shear zones. As a consequence, subduction initiation induced by regional mantle flow is demonstrably favoured by a free surface compared to the commonly applied, vertically fixed (i.e. free-slip) surface. In addition, we present global, three-dimensional mantle convection experiments that employ basal heating that leads to narrow mantle plumes. Narrow mantle plumes impinging on the base of the plate cause locally weak plate segments and a large topography at the lithosphere-asthenosphere boundary. Both are shown to be key to induce subduction initiation. Finally, our model self-consistently reproduces an episodic lid with a fast global overturn due to the hotter mantle developed below a former stagnant lid. We conclude that once in a stagnant-lid mode, a planet (like Venus) might preferentially evolve by temporally discrete, global overturn events rather than by a continuous recycling of lid and that this is something worth testing more rigorously in future studies.
Abstract: Series of high-resolution numerical simulations of three-dimensional mantle convection were performed to examine the interaction between the drifting continental lithospheres and the underlying mantle structure for 250 m.y. from the present, and to predict the configuration of the future supercontinent. The density anomaly of the mantle interior was determined by the seismic velocity anomaly from global seismic tomography data sets, which contain well-resolved subducting slabs. The present-day plate motion was imposed for the first stage of the simulation as a velocity boundary condition at the top surface boundary, instead of a shear stress-free condition. The switching time from the plate motion boundary to shear stress-free conditions was taken as a free parameter. The results revealed that Australia, Eurasia, North America, and Africa will merge together in the Northern Hemisphere to form a new supercontinent within ?250 m.y. from the present. The continental drift was assumed to be realized by plate-scale mantle flow, rather than large-scale upwelling plumes. That is, continuously moving plates at the surface for the first stage of the simulation are mechanically coupled with the subducting slabs in the mantle; this enhances the underlying mantle downwelling flow. As a result, persistent continental drift can be reproduced for long future time periods even though top surface boundary conditions may switch in response to shear stress-free conditions. The configuration of the numerically reproduced future supercontinent in this study is broadly consistent with the hypothetical model of Amasia as indicated by previous findings from geological correlations and a paleogeographic reconstruction.
Abstract: Lateral variations in lithosphere thickness are observed in many continental regions, especially at the boundary between the ancient cratonic core and the adjacent more juvenile lithosphere. In some places, such as the North America craton margin in western Canada and the Sorgenfrei-Tornquist Zone in northern Europe, the transition in lithosphere thickness has a steep gradient (>45°) and it appears to be a long-lived feature (at least 50 Ma). We use thermal-mechanical numerical models to address the dynamics of lithospheric thickness changes on timescales of 100 Ma. Models start with the juxtaposition of 60 km thick lithosphere ("mobile belt") and 160 km thick lithosphere ("craton"). In the reference model, all mantle materials have a damp olivine rheology and a density comparable to primitive mantle. With this configuration, edge-driven mantle convection occurs at the craton boundary, resulting in a lateral smoothing of the thickness transition. The density and rheology of the craton mantle lithosphere are then varied to approximate changes in composition and water content. For all densities, a steep transition is maintained only if the craton strength is 5-50 times stronger than the reference damp olivine. If dry olivine is an upper limit on strength, only cratonic mantle with moderate compositional buoyancy (20-40 kg/m3 less dense than primitive mantle) remains stable. At higher densities, the thick lithosphere is eroded through downwellings, and the craton margin migrates inboard. Conversely, a compositionally buoyant craton destabilises through lateral spreading below the mobile belt.
Abstract: High-Mg lavas are characteristic of the mid-Miocene volcanism in Inner Asia. In the Vitim Plateau, small volume high-Mg volcanics erupted at 16-14 Ma, and were followed with voluminous moderate-Mg lavas at 14-13 Ma. In the former unit, we have recorded a sequence of (1) initial basaltic melts, contaminated by crustal material, (2) uncontaminated high-Mg basanites and basalts of transitional (K-Na-K) compositions, and (3) picrobasalts and basalts of K series; in the latter unit a sequence of (1) initial basalts and basaltic andesites of transitional (Na-K-Na) compositions and (2) basalts and trachybasalts of K-Na series. From pressure estimation, we infer that the high-Mg melts were derived from the sub-lithospheric mantle as deep as 150 km, unlike the moderate-Mg melts that were produced at the shallow mantle. The 14-13 Ma rock sequence shows that initial melts equilibrated in a garnet-free mantle source with subsequently reduced degree of melting garnet-bearing material. No melting of relatively depleted lithospheric material, evidenced by mantle xenoliths, was involved in melting, however. We suggest that the studied transition from high- to moderate-Mg magmatism was due to the mid-Miocene thermal impact on the lithosphere by hot sub-lithospheric mantle material from the Transbaikalian low-velocity (melting) domain that had a potential temperature as high as 1510 °?. This thermal impact triggered rifting in the lithosphere of the Baikal Rift Zone.
Abstract: On the basis of quantum-chemical calculations of the linear to isomeric bent transition of the SiO2 molecule, it is suggested that the bent to linear transition of SiO2 forms can occur in melted mantle minerals of the lower mantle. This may be important for the formation of the peculiarities of mantle convection and origination of plumes.
Geochemistry, Geophysics, Geosystems: G3, Vol. 18, 7, pp. 2727-2747.
Canada, Somerset Island, Saskatchewan, United States, Kansas
magmatism, convection, diamond genesis
Abstract: Thirty new high-precision U-Pb perovskite and zircon ages from kimberlites in central North America delineate a corridor of mid-Cretaceous (115–92 Ma) magmatism that extends ?4000 km from Somerset Island in Arctic Canada through central Saskatchewan to Kansas, USA. The least contaminated whole rock Sr, Nd, and Hf isotopic data, coupled with Sr isotopic data from groundmass perovskite indicates an exceptionally limited range in Sr-Nd-Hf isotopic compositions, clustering at the low ?Nd end of the OIB array. These isotopic compositions are distinct from other studied North American kimberlites and point to a sublithospheric source region. This mid-Cretaceous kimberlite magmatism cannot be related to mantle plumes associated with the African or Pacific large low-shear wave velocity province (LLSVP). All three kimberlite fields are adjacent to strongly attenuated lithosphere at the edge of the North American craton. This facilitated edge-driven convection, a top-down driven processes that caused decompression melting of the transition zone or overlying asthenosphere. The inversion of ringwoodite and/or wadsleyite and release of H2O, with subsequent metasomatism and synchronous wet partial melting generates a hot CO2 and H2O-rich protokimberlite melt. Emplacement in the crust is controlled by local lithospheric factors; all three kimberlite fields have mid-Cretaceous age, reactivated major deep-seated structures that facilitated kimberlite melt transit through the lithosphere.
Journal of Asian Earth Sciences, Vol. 145, pt. B, pp. 334-348.
Mantle
convection, tectonics
Abstract: During Solar System condensation, the early Earth formed through planetesimal accretion, including collision of a Mars-sized asteroid. These processes rapidly increased the overall thermal budget and partial fusion of the planet. Aided by heat supplied by radioactivity and infall of the Fe-Ni core, devolatilization and chemical-density stratification attended planetary growth. After the thermal maximum at ?4.4 Ga, terrestrial temperatures gradually declined as an early Hadean magma ocean solidified. By ?4.3-4.2 Ga, H2O oceans + a dense CO2-rich atmosphere blanketed the terrestrial surface. Near-surface temperatures had fallen well below the low-P solidi of dry peridotite, basalt, and granite, ?1300, ?1120, and ?950 °C, respectively. At less than half their melting T, rocky materials existed as thin lithospheric platelets in the surficial Hadean Earth. Upper mantle stagnant-lid convection may have operated locally, but was rapidly overwhelmed by heat build-up-induced asthenospheric circulation, rifting and subduction, because massive heat transfer required vigorous mantle overturn in the early, hot planet. Bottom-up mantle overturn, involving abundant plume ascent, brought deep-seated heat to the surface. It decreased over time as cooling, plate enlargement, and top-down plate descent increased. Thickening, lateral extension, and contraction typified the post-Hadean lithosphere. Geologic evolutionary stages included: (a) ?4.5-4.4 Ga, the magma ocean solidified, generating ephemeral, ductile platelets; (b) ?4.4-2.7 Ga, small oceanic and continental plates were produced, then were destroyed by mantle return flow before ?4.0 Ga; eventually, continental material began to accumulate as largely subsea, sialic crust-capped lithospheric collages; (c) ?2.7-1.0 Ga, progressive suturing of old shields and younger orogenic belts led to cratonal plates typified by emerging continental freeboard, intense sedimentary differentiation, and episodic glaciation during transpolar plate drift; temporally limited stagnant-lid mantle convection occurred beneath growing supercontinents; (d) ?1.0 Ga-present, laminar-flowing mantle cells are capped by giant, stately moving plates. Near-restriction of komatiitic lavas to the Archean, and formation of multicycle sediments, ophiolite complexes ± alkaline igneous rocks, and high-pressure/ultrahigh-pressure (HP/UHP) metamorphic belts in youngest Proterozoic and Phanerozoic orogens reflect increasing density of cool oceanic plates, but decreasing subductability of enlarging, more buoyant continental plates. Attending assembly of supercontinents, negative buoyancy of thickening oceanic lithosphere began to control the overturn of suboceanic mantle as cold, top-down convection. The scales and dynamics of hot asthenospheric upwelling versus plate foundering and mantle return flow (bottom-up plume ascent versus top-down plate subduction) evolved gradually, due to planetary cooling. After accretion of the Earth, heat transfer through mantle convection has resulted in the existence of surficial rocky plates or platelets, and vigorous, lithosphere-coupled mantle overturn since ?4.4 Ga. Thus plate-tectonic processes have typified the Earth’s thermal history since Hadean time.
Abstract: Mantle convection is a fundamental planetary process. Its plate mode is established and expressed by plate tectonics. Its plume mode also is established and expressed by interregional geological patterns. We developed both an event-based stratigraphic framework to illustrate the surface effects predicted by the plume model of Griffiths et al. (1989) and Griffiths and Campbell (1990) and a methodology to analyze continent-scale geological maps based on unconformities and hiatuses. The surface expression of ascending plumes lasts for tens-of-millions-of-years and rates vary over a few million years. As the plume ascends, its surface expression narrows, but increases in amplitude, leaving distinct geological and stratigraphic patterns in the geologic record, not only above the plume-head center, but also above its margins and in distal regions a few thousands-of-kilometers from the center. To visualize these patterns, we constructed sequential geological maps, chronostratigraphic sections, and hiatus diagrams. Dome-uplift with erosion (?engör, 2001) and the flood basalts (Duncan and Richards, 1991; Ernst and Buchan, 2001a) are diagnostic starting points for plume-stratigraphic analyses. Mechanical collapse of the dome results in narrow rifting (Burke and Dewey, 1973), drainage-network reorganization (Cox, 1989), and flood-basalt eruption. In the marginal region, patterns of vertical movement, deformation and surface response are transient and complex. At first, the plume margin is uplifted together with the central region, but then it subsides as the plume ascents farther; With plume-head flattening, the plume margin experiences renewed outward-migrating surface uplift, erosion, broad crustal faulting, and drainage reorganization. Knickpoint migration occurs first inward-directed at ˝ the rate of plume ascent and later outward-directed at the rate of asthenospheric flow. Interregional-scale unconformity-bounded stratigraphic successions document the two inversions. The distal regions, which did not experience any plume-related uplift, yield complete sedimentary records of the event; Event-related time gaps (hiatuses) in the sedimentary record increase towards the center, but the event horizon is best preserved in the distal region; it may be recognized by tracing its contacts from the center outwards. We extracted system- and series-hiatuses from interregional geological maps and built hiatus maps as proxies for paleo-dynamic topography and as a basis for comparison with results from numerical models. Interregional-scale geological maps are well suited to visualize plume-related geological records of dynamic topography in continental regions. However, geological records and hiatus information at the resolution of stages will be needed at interregional scales. The plume-stratigraphic framework is event-based, interregional, but not global, with time-dependent amplitudes that are significantly larger than those of global eustatic sea-level fluctuations. Global stratigraphic syntheses require integration of plate- and plume-stratigraphic frameworks before eustatic contributions may be assessed.
Nonlinear Processes Geophysics, Vol. 25, pp. 99-123. pdf
Mantle
convection
Abstract: Recent advances in mantle convection modeling led to the release of a new generation of convection codes, able to self-consistently generate plate-like tectonics at their surface. Those models physically link mantle dynamics to surface tectonics. Combined with plate tectonic reconstructions, they have the potential to produce a new generation of mantle circulation models that use data assimilation methods and where uncertainties in plate tectonic reconstructions are taken into account. We provided a proof of this concept by applying a suboptimal Kalman filter to the reconstruction of mantle circulation (Bocher et al., 2016). Here, we propose to go one step further and apply the ensemble Kalman filter (EnKF) to this problem. The EnKF is a sequential Monte Carlo method particularly adapted to solve high-dimensional data assimilation problems with nonlinear dynamics. We tested the EnKF using synthetic observations consisting of surface velocity and heat flow measurements on a 2-D-spherical annulus model and compared it with the method developed previously. The EnKF performs on average better and is more stable than the former method. Less than 300 ensemble members are sufficient to reconstruct an evolution. We use covariance adaptive inflation and localization to correct for sampling errors. We show that the EnKF results are robust over a wide range of covariance localization parameters. The reconstruction is associated with an estimation of the error, and provides valuable information on where the reconstruction is to be trusted or not.
Earth and Planetary Science Letters, Vol. 496, 1, pp. 142-158.
Asia, Borneo
convection
Abstract: Most, but not all, geodynamic models predict 1-2 km of mantle convective draw-down of the Earth's surface in a region centered on Borneo within southeast Asia. Nevertheless, there is geomorphic, geologic and geophysical evidence which suggests that convective uplift might have played some role in sculpting Bornean physiography. For example, a long wavelength free-air gravity anomaly of +60 mGal centered on Borneo coincides with the distribution of Neogene basaltic magmatism and with the locus of sub-plate slow shear wave velocity anomalies. Global positioning system measurements, an estimate of elastic thickness, and crustal isostatic considerations suggest that regional shortening does not entirely account for kilometer-scale regional elevation. Here, we explore the possible evolution of the Bornean landscape by extracting and modeling an inventory of 90 longitudinal river profiles. Misfit between observed and calculated river profiles is minimized by smoothly varying uplift rate as a function of space and time. Erosional parameters are chosen by assuming that regional uplift post-dates Eocene deposition of marine carbonate rocks. The robustness of this calibration is tested against independent geologic observations such as thermochronometric measurements, offshore sedimentary flux calculations, and the history of volcanism. A calculated cumulative uplift history suggests that kilometer-scale Bornean topography grew rapidly during Neogene times. This suggestion is corroborated by an offshore Miocene transition from carbonate to clastic deposition. Co-location of regional uplift and slow shear wave velocity anomalies immediately beneath the lithospheric plate implies that regional uplift could have been at least partly generated and maintained by temperature anomalies within an asthenospheric channel.
Abstract: The linear structures of seismically fast anomalies, often interpreted as subducted slabs, in the southern Asia and circum-Pacific lower mantle provided strong evidence for the whole mantle convection model. However, recent seismic studies have consistently shown that subducted slabs are deflected horizontally for large distances in mantle transition zone in the western Pacific and other subduction zones, suggesting that the slabs meet significant resistance to their descending motion and become stagnant in the transition zone. This poses challenges to the whole mantle convection model and also brings the origin of stagnant slabs into question. Here, using a global mantle convection model with realistic spine-post-spinel phase change (?2 MPa K?ą Clapeyron slope) and plate motion history, we demonstrate that the observed stagnant slabs in the transition zone and other slab structures in the lower mantle can be explained by the presence of a thin, weak layer at the phase change boundary that was suggested by mineral physics and geoid modelling studies. Our study also shows that the stagnant slabs mostly result from subduction in the past 20-30 million years, confirming the transient nature of slab stagnation and phase change dynamics on timescales of tens of millions of years from previous studies.
Abstract: The existence of undulations of the geoid, gravity and bathymetry in ocean basins, as well as anomalies in heat flow, point to the existence of small scale convection beneath tectonic plates. The instabilities that could develop at the base of the lithosphere are sufficiently small scale (< 500 km) that they remain mostly elusive from seismic detection. We take advantage of 3D spherical numerical geodynamic models displaying plate-like behavior to study the interaction between large-scale flow and small-scale convection. We find that finger-shaped instabilities develop at seafloor ages > 60 Ma. They form networks that are shaped by the plate evolution, slabs, plumes and the geometry of continental boundaries. Plumes impacting the boundary layer from below have a particular influence through rejuvenating the thermal lithosphere. They create a wake in which new instabilities form downstream. These wakes form channels that are about 1000 km wide, and thus are possibly detectable by seismic tomography. Beneath fast plates, cold sinking instabilities are tilted in the direction opposite to plate motion, while they sink vertically for slow plates. These instabilities are too small to be detected by usual seismic methods, since they are about 200 km in lateral scale. However, this preferred orientation of instabilities below fast plates could produce a pattern of large-scale azimuthal anisotropy consistent with both plate motions and the large scale organisation of azimuthal anisotropy obtained from recent surface wave models.
Abstract: Three?dimensional spherical mantle convection was simulated to predict future continental motion and investigate the driving force of continental motion. Results show that both the time required (?300 Ma from the present) and the process for the next supercontinent formation are sensitive to the choice of critical rheological parameters for mantle dynamics, such as a viscosity contrast between the upper and lower mantles and a yield strength of the lithosphere. From all the numerical models studied herein, mantle drag force by horizontal mantle flow beneath the continents may mostly act as a resistance force for the continental motion in the process of forming a new supercontinent. The maximum absolute magnitude of the tensional and compressional stress acting at the base of the moving continents is in the order of 100 MPa, which is comparable to a typical value of the slab pull force.
Abstract: Subducting slabs are the primary drivers of plate tectonics and mantle circulation, but can also undergo various instabilities that cause dramatic adjustments in tectonic evolution and motion. Slab rollback or trench retreat is possibly a dominant form of time dependence in the plate-mantle system, causing plates to shrink and the mantle to undergo complex flow patterns. Likewise, slab detachment can induce abrupt adjustments in both plate motions and vertical displacement of continents. The arrival or accumulation of continental crust over a subduction zone induces high stresses on the plate and slab that can trigger either rollback or detachment or both. However, these processes necessarily interact because of how stress is relieved and plate motions altered. Here we present a simple boundary-layer like model of coupled trench retreat and slab detachment, induced by continent accumulation, and with slab necking augmented by grain-damage self-weakening (to allow for abrupt necking). With this model we find that, with continental accumulation, initial rollback is at first modest. However, as the stress from continental accumulation peaks, it triggers abrupt slab detachment. The subsequent slab loss causes the plate to lose its primary motive force and to thus undergo a more dramatic and rapid rollback event. After the larger rollback episode, the contracted continental mass re-expands partially. Plausible grain-damage parameters and 40?km thick crust cause abrupt detachment and major rollback to occur after a few hundred million years, which means the plates remain stable for that long, in agreement with the typical age for most large plates. While the complexity of some field areas with a well documented history of detachment and rollback, such as the Mediterranean, taxes the sophistication of our toy model, other simpler geological examples, such as on the western North American plate, show that episodes of rollback can follow detachment.
Abstract: The structure of the lithosphere is key to reconciling the dynamic topography predicted by mantle convection models with residual topography derived from observations, suggest analyses of both models and data.
Physics of the Earth and Planetary Interiors, in press 10.1016/j.pepi.2019.106299 18p. Pdf
Mantle
convection
Abstract: The solid Earth system is characterized by plate tectonics, a low viscosity zone beneath plates (the asthenosphere), and long wavelength flow in the convecting mantle. We use suites of numerical experiments to show: 1) How long wavelength flow and the operation of plate tectonics can generate and maintain an asthenosphere, and 2) How an asthenosphere can maintain long wavelength flow and plate tectonics. Plate subduction generates a sub-adiabatic temperature gradient in the mantle which, together with temperature-dependent viscosity, leads to a viscosity increase from the upper to the lower mantle. This allows mantle flow to channelize in a low viscosity region beneath plates (an asthenosphere forms dynamically). Flow channelization, in turn, stabilizes long wavelength convection. The degree of dynamic viscosity variations from the upper to the lower mantle increases with the wavelength of convection and drops toward zero if the system transitions from plate tectonics to a single plate planet. The plate margin strength needed to initiate that transition increases for long wavelength cells (long wavelength flow allows plate tectonics to exist over a wider range of plate margin strength). The coupled feedbacks allow for a linked causality between plates, the asthenosphere, and the wavelength of mantle flow, with none being more fundamental than the others and the existence of each depending on the others. Under this hypothesis, the asthenosphere is defined by an active process, plate tectonics, which maintains it and is maintained by it and plate tectonics is part of an emergent, self-sustaining flow system that bootstraps itself into existence.
Abstract: The exchange of volatile species—water, carbon dioxide, nitrogen and halogens—between the mantle and the surface of the Earth has been a key driver of environmental changes throughout Earth’s history. Degassing of the mantle requires partial melting and is therefore linked to mantle convection, whose regime and vigour in the Earth’s distant past remain poorly constrained1,2. Here we present direct geochemical constraints on the flux of volatiles from the mantle. Atmospheric xenon has a monoisotopic excess of 129Xe, produced by the decay of extinct 129I. This excess was mainly acquired during Earth’s formation and early evolution3, but mantle degassing has also contributed 129Xe to the atmosphere through geological time. Atmospheric xenon trapped in samples from the Archaean eon shows a slight depletion of 129Xe relative to the modern composition4,5, which tends to disappear in more recent samples5,6. To reconcile this deficit in the Archaean atmosphere by mantle degassing would require the degassing rate of Earth at the end of the Archaean to be at least one order of magnitude higher than today. We demonstrate that such an intense activity could not have occurred within a plate tectonics regime. The most likely scenario is a relatively short (about 300 million years) burst of mantle activity at the end of the Archaean (around 2.5 billion years ago). This lends credence to models advocating a magmatic origin for drastic environmental changes during the Neoarchaean era, such as the Great Oxidation Event.
Geochemical Perspectives Letters, Vol. 9, pp. 21-25.
Mantle
convection
Abstract: Noble gases serve as unique tracers of the origin and evolution of Earth’s volatile reservoirs owing to their inert nature and contribution from extinct and extant radioactivities. However, noble gases are low in abundance relative to many other elements, particularly in the Earth’s mantle. Additionally, mantle-derived samples show large post-eruptive atmospheric contamination, rendering the determination of the primary mantle composition challenging. The sources of mantle krypton and xenon remain debated due to their partially resolvable excess, if any, relative to the atmosphere. Atmospheric noble gases also appear to be recycled into the mantle via subduction, progressively overprinting the initial mantle signature. Here we develop a new protocol to accumulate non-contaminated mantle-derived xenon, in particular the low abundant 124-126-128Xe. The results show the highest excesses in 124-126-128Xe ever measured in the mantle relative to the atmosphere and point toward a chondritic origin for mantle xenon. The fissiogenic isotopes 131-132-134-136Xe allow the onset of efficient xenon recycling in the mantle to be constrained at around 3 Gyr ago, implying that volatile recycling before 3 Ga would have been negligible.
Geochimica et Cosmochimica Acta, Vol. 244, pp. 56-85.
Mantle
convection
Abstract: Atmospheric xenon is strongly mass fractionated, the result of a process that apparently continued through the Archean and perhaps beyond. Previous models that explain Xe fractionation by hydrodynamic hydrogen escape cannot gracefully explain how Xe escaped when Ar and Kr did not, nor allow Xe to escape in the Archean. Here we show that Xe is the only noble gas that can escape as an ion in a photo-ionized hydrogen wind, possible in the absence of a geomagnetic field or along polar magnetic field lines that open into interplanetary space. To quantify the hypothesis we construct new 1-D models of hydrodynamic diffusion-limited hydrogen escape from highly-irradiated CO2-H2-H atmospheres. The models reveal three minimum requirements for Xe escape: solar EUV irradiation needs to exceed that of the modern Sun; the total hydrogen mixing ratio in the atmosphere needs to exceed 1% (equiv. to CH4); and transport amongst the ions in the lower ionosphere needs to lift the Xe ions to the base of the outflowing hydrogen corona. The long duration of Xe escape implies that, if a constant process, Earth lost the hydrogen from at least one ocean of water, roughly evenly split between the Hadean and the Archean. However, to account for both Xe’s fractionation and also its depletion with respect to Kr and primordial 244Pu, Xe escape must have been limited to small apertures or short episodes, which suggests that Xe escape was restricted to polar windows by a geomagnetic field, or dominated by outbursts of high solar activity, or limited to transient episodes of abundant hydrogen, or a combination of these. Xenon escape stopped when the hydrogen (or methane) mixing ratio became too small, or EUV radiation from the aging Sun became too weak, or charge exchange between Xe+ and O2 rendered Xe neutral. In our model, Xe fractionation attests to an extended history of hydrogen escape and Earth oxidation preceding and ending with the Great Oxidation Event (GOE).
Abstract: Two end-member scenarios have been proposed for the tectonic situation along the eastern margins of Gondwanaland before Zealandia was formed ca. 100 million years ago (Ma), namely: (1) A subduction zone located far from the eastern margin of Zealandia, wherein Zealandia may have separated from Gondwanaland by plume push of an active hotspot plume.; (2) A subduction zone located along the eastern margin of Gondwanaland, wherein Zealandia possibly separated from Gondwanaland via trench/subduction retreat. Assuming that the thermal structure of the deep mantle and source of hotspot plumes remained relatively stationary over the last hundred million years, major hotspot plumes with a large buoyancy flux did not exist under Zealandia; the eastern margins of Gondwanaland were far from two large low-shear-velocity provinces under the Africa-Atlantic and South Pacific regions. Herein, through numerical studies of three-dimensional global mantle convection, we examined the mantle convection and surface tectonic patterns at ~100 Ma. The present model considered the real configuration of Gondwanaland at the model surface to observe long-term variations of mantle convection and the resulting surface tectonic conditions. The results demonstrate that the extensive subduction zones developed preferentially along the eastern margin of Gondwanaland when the temperature anomaly of the lower mantle was primarily dominated by high-temperature regions under present-day Africa-Atlantic and South Pacific regions. The results of this study support one of the proposed hypotheses, where the breakup at the eastern margins of Gondwanaland at ~100 Ma occurred via trench/subduction retreat.
Abstract: Paleo-temperature data indicates that the Earth's mantle did not cool at a constant rate over geologic time. Post magma ocean cooling was slow with an onset of more rapid mantle cooling between 2.5 and 3.0 Gyr. We explore the hypothesis that this multi-stage cooling is a result of deep water cycling coupled to thermal mantle convection. As warm mantle ascends, producing melt, the mantle is dehydrated. This tends to stiffens the mantle, which slows convective vigor causing mantle heating. At the same time, an increase in temperature tends to lower mantle viscosity which acts to increase convective vigor. If these two tendencies are in balance, then mantle cooling can be weak. If the balance is broken, by a switch to a net rehydration of the mantle, then the mantle can cool more rapidly. We use coupled water cycling and mantle convection models to test the viability of this hypothesis. We test models with different parameterizations to allow for variable degrees of plate margin strength. We also perform a layered uncertainty analysis on all the models to account for input, parameter, and structural model uncertainties. Within model and data uncertainty, the hypothesis that deep water cycling, together with a combination of plate strength and mantle viscosity resisting mantle overturn, can account for paleo data constraints on mantle cooling.
Abstract: Conventional wisdom holds that the motion of tectonic plates drives motion in the Earth’s rocky interior (i.e., in the Earth’s asthenosphere). Recent seismological observations have brought this view into question as they indicate that the velocity of the asthenosphere can exceed tectonic plate velocity. This suggests that interior motions can drive plate motions. We explore models of coupled plate tectonics and interior motions to address this discrepancy. The models reveal that the coupling between plates and the asthenosphere is not an issue of plates drive asthenosphere motion or asthenosphere motion drives plates. Both factors work in tandem with the balance being a function of plate margins strength and asthenosphere rheology. In particular, a power-law viscosity allows pressure gradients to generate interior flow that can locally drive plate motion. The models also reveal a hysteresis effect that allows different tectonic states (plate tectonics versus a single plate planet) to exist at the same parameter conditions. This indicates that history and initial conditions can play a role in determining if a planet will or will not have plate tectonics.
Geophysical Research Letters, 10.1029/2020/GL089556 11p. Pdf
Mantle
convection
Abstract: It is generally thought that tectonic plates drive motion in the Earth's rocky interior. Recent observations have challenged this view as they indicate that interior motion can drive tectonic plates. Models of coupled tectonics and interior flow are used to address this discrepancy. The models reveal that the question of “does plate tectonics drive interior flow or does interior flow drive plate tectonics” may be ill founded as both possibilities may be active at the same time. The balance between the two drivers is found to depend on plate margin strength. The models also reveal that different tectonic modes can exist under the same physical conditions. This indicates a planet's initial state can determine if it will or will not have plate tectonics.