Mantle Flow Models Research Articles

Dynamics of the North American Plate: Large‐Scale Driving Mechanism From Far‐Field Slabs and the Interpretation of Shallow Negative Seismic Anomalies

AbstractWith a small fraction of marginal subduction zones, the driving mechanism for the North American plate motion is in debate. We construct global mantle flow models simultaneously constrained by geoid and plate motions to investigate the driving forces for the North American plate motion. By comparing the model with only near‐field subducting slabs and that with global subducting slabs, we find that the contribution to the motion of the North American plate from the near‐field Aleutian, central American, and Caribbean slabs is small. In contrast, other far‐field slabs, primarily the major segments around western Pacific subduction margins, provide the dominant large‐scale driving forces for the North American plate motion. The coupling between far‐field slabs and the North American plate suggests a new form of active plate interactions within the global self‐organizing plate tectonic system. We further evaluate the extremely slow seismic velocity anomalies associated with the shallow partial melt around the southwestern North America. Interpreting these negative seismic shear‐velocity anomalies as purely thermal origin generates considerably excessive resistance to the North American plate motion. A significantly reduced velocity‐to‐density scaling for these negative seismic shear‐velocity anomalies must be incorporated into the construction of the buoyancy field to predict the North American plate motion. We also examine the importance of lower mantle buoyancy including the ancient descending Kula‐Farallon plates and the active upwelling below the Pacific margin of the North American plate. Lower mantle buoyancy primarily affects the amplitudes, as opposed to the patterns of both North American and global plate motions.

Strong ULVZ and Slab Interaction at the Northeastern Edge of the Pacific LLSVP Favors Plume Generation

AbstractStrong waveform complexities, including multipathing of the S diffracted phase and rapid changes in differential ScS‐S times, are observed for multiple deep Fiji earthquakes recorded at the USArray. The complexities occur at the northeastern edge of the Pacific Large Low Shear Velocity Province (LLSVP), about 12 degrees southeast of present‐day Hawaiʻi. Waveform modeling of the multipathing provides good constraints on an ultra‐low velocity zone (ULVZ) with a width of 5 degree located near the inner edge of the LLSVP. Based on the mineralogical‐modeling of the ULVZ as a solid iron‐rich magnesiowüstite‐bearing assemblage with compatible morphology predicted from geodynamical simulations, a ULVZ model with a thickness of 30 km and a shear wave velocity reduction of 18% is preferred. The rapid change in differential ScS‐S travel time is best explained by having both the aforementioned ULVZ and an adjacent high velocity structure near the LLSVP. Furthermore, a low‐velocity plume‐like structure could potentially explain the observed S travel time delay independent of ScS. These seismic features are proposed to be a ULVZ driven toward the edge of the LLSVP while potentially pushed by a subducted slab. This configuration may trigger plume generation due to strong thermal instabilities and is in the same vicinity where mantle flow models place the present‐day Hawaiian plume source. Multiple ScS can potentially be used to verify vertical plume structure in tomographic models but the accuracy of upper mantle structure, which is a key reflection point, needs to be considered.

The evolution of basal mantle structure in response to supercontinent aggregation and dispersal

Seismic studies have revealed two Large Low-Shear Velocity Provinces (LLSVPs) in the lowermost mantle. Whether these structures remain stable over time or evolve through supercontinent cycles is debated. Here we analyze a recently published mantle flow model constrained by a synthetic plate motion model extending back to one billion years ago, to investigate how the mantle evolves in response to changing plate configurations. Our model predicts that sinking slabs segment the basal thermochemical structure below an assembling supercontinent, and that this structure eventually becomes unified due to slab push from circum-supercontinental subduction. In contrast, the basal thermochemical structure below the superocean is generally coherent due to the persistence of a superocean in our imposed plate reconstruction. The two antipodal basal thermochemical structures exchange material several times when part of one of the structures is carved out and merged with the other one, similarly to “exotic” tectonic terranes. Plumes mostly rise from thick basal thermochemical structures and in some instances migrate from the edges towards the interior of basal thermochemical structures due to slab push. Our results suggest that the topography of basal structures and distribution of plumes change over time due to the changing subduction network over supercontinent cycles.

High 3He/4He in central Panama reveals a distal connection to the Galápagos plume

It is well established that mantle plumes are the main conduits for upwelling geochemically enriched material from Earth's deep interior. The fashion and extent to which lateral flow processes at shallow depths may disperse enriched mantle material far (>1,000 km) from vertical plume conduits, however, remain poorly constrained. Here, we report He and C isotope data from 65 hydrothermal fluids from the southern Central America Margin (CAM) which reveal strikingly high 3He/4He (up to 8.9RA) in low-temperature (≤50 °C) geothermal springs of central Panama that are not associated with active volcanism. Following radiogenic correction, these data imply a mantle source 3He/4He >10.3RA (and potentially up to 26RA, similar to Galápagos hotspot lavas) markedly greater than the upper mantle range (8 ± 1RA). Lava geochemistry (Pb isotopes, Nb/U, and Ce/Pb) and geophysical constraints show that high 3He/4He values in central Panama are likely derived from the infiltration of a Galápagos plume-like mantle through a slab window that opened ∼8 Mya. Two potential transport mechanisms can explain the connection between the Galápagos plume and the slab window: 1) sublithospheric transport of Galápagos plume material channeled by lithosphere thinning along the Panama Fracture Zone or 2) active upwelling of Galápagos plume material blown by a "mantle wind" toward the CAM. We present a model of global mantle flow that supports the second mechanism, whereby most of the eastward transport of Galápagos plume material occurs in the shallow asthenosphere. These findings underscore the potential for lateral mantle flow to transport mantle geochemical heterogeneities thousands of kilometers away from plume conduits.

Modelling the role of dynamic topography and eustasy in the evolution of the Great Artesian Basin

AbstractWidespread flooding of the Australian continent during the Early Cretaceous, referred to as the Eromanga Sea, deposited extensive shallow marine sediments throughout the Great Artesian Basin (GAB). This event had been considered ‘out of sync’ with eustatic sea level and was instead solely attributed to dynamic subsidence associated with Australia's passage over eastern Gondwanan subducted material. However, mantle convection models previously used to explain this event have since been shown to overestimate dynamic topography amplitude by a factor of two compared with residual topography estimates. Previous models were also based on a Cretaceous eustatic sea level peak at ca. 90 Ma in conventional eustatic sea level curves; however, more recent estimates of global sea level from ocean basin volume (OBV) suggest this peak may have occurred earlier at ca. 120 Ma. Our work links time‐dependent erosion and deposition with dynamic topography and eustasy to test their contribution to basin development using the landscape evolution code pyBadlands. Our results show that the lower amplitude estimates of dynamic topography derived from pseudo‐compressible mantle flow models better reflect the Cretaceous vertical motions of the Australian continent (ca. 100 m) compared with their incompressible counterparts (ca. 200–400 m). Additionally, our models include the Neogene north‐eastward tilting of Australia, elusive in most previously published geodynamic models. In conjunction with an OBV‐derived sea level curve, our preferred landscape evolution model broadly matches the Cretaceous inundation patterns and first‐order sedimentary sequences in the GAB. The results highlight that the Early Cretaceous inundation of the Australian continent is likely a combination of high global sea levels and the regional effects of dynamic subsidence. Our work provides a framework for a new generation of evolving paleogeographic models at continental scales, while also providing key insights into the viability of existing sea level curves and dynamic topography estimates for reproducing topographic and basin evolution.

Spatial Characteristics of Recycled and Primordial Reservoirs in the Deep Mantle

AbstractThe spatial distribution of the geochemical domains hosting recycled crust and primordial (high‐3He/4He) reservoirs, and how they are linked to mantle convection, are poorly understood. Two continent‐sized seismic anomalies located near the core‐mantle boundary—called the Large Low Shear Wave Velocity Provinces (LLSVPs)—are potential geochemical reservoir hosts. It has been suggested that high‐3He/4He hotspots are spatially confined to the LLSVPs, hotspots sampling recycled continental crust are associated with only one of the LLSVPs, and recycled continental crust shows no relationship with latitude. We reevaluate the links between LLSVPs and isotopic signatures of hotspot lavas using improved mantle flow models including plume conduit advection. While most hotspots with the highest‐3He/4He can indeed be traced to the LLSVP interiors, at least one high‐3He/4He hotspot, Yellowstone, is located outside of the LLSVPs. This suggests high‐3He/4He is not geographically confined to the LLSVPs. Instead, a positive correlation between hotspot buoyancy flux and maximum hotspot 3He/4He suggests that it is plume dynamics (i.e., buoyancy), not geography, which determines whether a dense, deep, and possibly widespread high‐3He/4He reservoir is entrained. We also show that plume‐fed EM hotspots (enriched mantle, with low‐143Nd/144Nd), signaling recycled continental crust, are spatially linked to both LLSVPs, and located primarily in the southern hemisphere. Lastly, we confirm that hotspots sampling HIMU (“high‐μ,” or high 238U/204Pb) domains are not spatially limited to the LLSVPs. These findings clarify and advance our understanding of deep mantle reservoir distributions, and we discuss how continental and oceanic crust subduction is consistent with the spatial decoupling of EM and HIMU.

Controls on the present-day dynamic topography predicted from mantle flow models since 410 Ma

SUMMARYMantle convection induces dynamic topography, the lithosphere's surface deflections driven by the vertical stresses from sublithospheric mantle convection. Dynamic topography has important influences on a range of geophysical and geological observations. Here, we studied controls on the Earth's dynamic topography through 3-D spherical models of mantle convection, which use reconstructed past 410 Myr global plate motion history as time-dependent surface mechanical boundary condition. The numerical model assumes the extended-Boussinesq approximation and includes strongly depth- and temperature-dependent viscosity and phase changes in the mantle. Our results show that removing the chemical layer above the core–mantle boundary (CMB) and including depth-dependent thermal expansivity have both a limited influence on the predicted present-day dynamic topography. Considering phase transitions in our models increases the predicted amplitude of dynamic topography, which is mainly influenced by the 410 km exothermic phase transition. The predicted dynamic topography is very sensitive to shallow temperature-induced lateral viscosity variations (LVVs) and Rayleigh number. The preservation of LVVs significantly increases the negative dynamic topography at subduction zones. A decrease (or increase) of Rayleigh number increases (or decreases) the predicted present-day dynamic topography considerably. The dynamic topography predicted from the model considering LVVs and with a Rayleigh number of 6 × 108 is most compatible with residual topography models. This Rayleigh number is consistent with the convective vigor of the Earth as supported by generating more realistic lower mantle structure, slab sinking rate and surface and CMB heat fluxes. The evolution of the surface heat flux pattern is similar to the long-term eustatic sea level change. Before the formation of Pangea, large negative dynamic topography formed between the plate convergence region of Gondwana and Laurussia. The predicted dynamic topography similar to that of present-day has already emerged by about 262 Ma. Powers for degrees 1–3 dynamic topography at 337 and 104 Ma which correspond to times of higher plate velocities and higher surface heat fluxes are larger.

Effects of basal drag on subduction dynamics from 2D numerical models

Abstract. Subducting slabs are an important driver of plate motions, yet the relative importance of different forces in governing subduction motions and styles remains incompletely understood. Basal drag has been proposed to be a minor contributor to subduction forcing because of the lack of correlation between plate size and velocity in observed and reconstructed plate motions. Furthermore, in single subduction system models, low basal drag leads to subduction behaviour most consistent with the observation that trench migration velocities are generally low compared to convergence velocities. By contrast, analytical calculations and global mantle flow models indicate basal drag can be substantial. In this study, we revisit this problem by examining the drag at the base of the lithosphere, for a single subduction system, in 2D models with a free trench and composite non-linear rheology. We compare the behaviour of short and long plates for a range of asthenospheric and lithospheric rheologies. We reproduce results from previous modelling studies, including low ratios of trench over plate motions. However, we also find that any combination of asthenosphere and lithosphere viscosity that produces Earth-like subduction behaviour leads to a correlation of velocities with plate size, due to the role of basal drag. By examining Cenozoic plate motion reconstructions, we find that slab age and plate size are positively correlated: higher slab pull for older plates tends to be offset by higher basal drag below these larger plates. This, in part, explains the lack of plate velocity–size correlation in observations, despite the important role of basal drag in the subduction force balance.

The influence of mantle flow on intracontinental basins: Three examples from Australia

AbstractDuring the Paleozoic, sedimentary basins developed within Gondwana without evolving to diverging plate boundaries. Such intracontinental basins present long subsidence histories with multiple phases of accelerated subsidence that are not always easily explained by far‐field tectonic forces, and may be driven by processes other than rifting and thermal subsidence. Here we investigate the subsidence of Paleozoic Australian intracontinental basins by comparing one‐dimensional backstripped tectonic subsidence histories from the western Australian Canning and Southern Carnarvon Basins and the central Australian Cooper Basin to forward subsidence models for pure shear lithospheric thinning. We make the hypothesis that differences between observed and model subsidence may be explained by mantle‐flow driven topography, in addition to tectonic forces. To test this hypothesis, we compute dynamic topography from the first geodynamic models of mantle flow spanning the entire Phanerozoic Eon, and we analyse the relationship between dynamic topography and anomalous basin subsidence to dynamic topography and mantle flow. Although reconstructions of mantle flow in deep geological times are uncertain, our results suggest that long‐wavelength dynamic topography could explain aspects of the complex tectonic histories intracontinental basins. In the presented reconstruction of mantle flow, topographic rebound following the sinking of a Cambrian aged slab resulted in a minor phase of dynamic uplift in the Cooper Basin in middle Permian times. Throughout Carboniferous‐Triassic times Australia was positioned above a mantle upwelling driven by a hot structure at the base of the mantle. Structural uplift in the Canning and Southern Carnarvon basins during the Triassic‐Jurassic interval was augmented by dynamic uplift produced by that large‐scale upwelling, and possibly augmented by a focused active mantle plume during the Permo‐Triassic. In Late Jurassic‐Cretaceous times, Australia drifted east away from the mantle upwelling, resulting in a period of subsidence in the Canning and Southern Carnarvon basins. During the Cretaceous the Cooper Basin moved over a downwelling produced by long‐lived subduction along the east Australian margin, resulting in a period of accelerated subsidence.

Lithospheric mantle anisotropy from local events beneath the Sunda–Banda arc transition and its geodynamic implications

Shear wave splitting analysis to characterise lithospheric mantle anisotropy has been performed to provide better knowledge about lithospheric deformation and mantle flow beneath the Sunda–Banda arc transition, Indonesia. The tectonic setting of the study area is very complex characterised by the transition from subduction along Sunda arc to collision in Banda arc. The splitting measurements show lateral and vertical variation in the fast directions of the S-waves in this region. When the splitting results are analysed through 2D delay-time tomography and spatial averaging, systematic patterns in delay times and fast polarisation become more visible. In the subduction domain, the spatial averages of fast directions are dominated by two distinct fast polarisations: perpendicular and parallel to the plate motion for shallow and deep events, respectively. The results suggest that anisotropy in this area is not only controlled by anisotropic source related to the simple mantle flow model, but also by anisotropic fabric in the mantle deformed under influence of high stresses, high water contents and low temperatures. In addition, there might also be contribution from the anisotropic body in the upper layer. In the collision domain, spatially averaged fast directions show mostly perpendicular to the plate motion for all deep levels. For shallow level in this region, this trend is mainly governed by the lithospheric deformation process due to the continent-arc collision as also shown by delay time tomographic inversion. For deeper part of the region, the result of tomographic inversion and spatial averaging reveals a high anisotropy followed by rotational pattern of fast directions in the north of Timor. We suggest that this pattern might be related to the induced mantle flow due to lateral tearing of the slab.

Strength variations of the Australian continent: Effects of temperature, strain rate, and rheological changes

The Australian continent is composed of several geologic provinces, showing a general age progression from Archean in the west to Phanerozoic in the east. The lithospheric heterogeneity and complex tectonic history of this region make it a key area for studying the thermal and rheological structure of the geological provinces and testing the influence of different conditions, such as temperature, rheology, and strain rate on the plate strength. In a previous study, temperature and compositional variations of the Australian upper mantle have been determined based on a joint interpretation of the seismic tomography and gravity data. In this study, we further implement a thermal model of the crust, based on available surface heat flow data from regional and global database. The crustal and upper mantle thermal models show different anomalies distribution, indicating a significant variation of the thermal conditions with depth. The new thermal models are used to estimate strength and effective elastic thickness ( Te ) distribution in the lithosphere. For this aim, we assigned the rheology of the crust based on the seismic velocities provided by the AuSREM model and used the strain rate values obtained from a global mantle flow model, constrained by seismic and gravity data. The maximal strength and Te are found in the West Australian Craton, on account of the low temperatures in the lithospheric mantle. We found that locations of the intraplate earthquakes attend to sharp changes in the lithospheric strength. Comparison of the results with those obtained for uniform rheology and strain rate, indicate that in the Officer basin the variations of the crustal rheology enhance the effect of temperature changes, while in the Yilgarn craton they reduce it. On the other hand, the lower values of the strain rate in the cratons than in the Phanerozoic regions influence the strength/ Te in the opposite way with respect to temperatures. Integrated strength (Pa m) and strength (Pa) variations along three cross-sections. Horizontal and vertical axes of the cross-sections display the distance and depth (both in km), respectively. (a) Integrated strength of the lithosphere. The three black lines show the location of the cross-sections; (b) Integrated strength of the crust. White circles show the earthquakes location, occurred from 1990 to present-day, with moment magnitudes (Mw) >2.5, taken from the USGS seismic catalogue ( https://earthquake.usgs.gov/earthquakes/search/ ); (c) Integrated strength of the mantle lithosphere; (d) Percentage of crustal strength; (e) Strength variations along three cross-sections (A-A′, B-B′, and C-C′, respectively). • A crustal themal and rheological model of the Australian continent is set up. • Crustal thermal anomalies correlate with surface heat flow variations. • Upper mantle thermal anomalies correlate with age of the tectonic features. • The strength distribution is strongly influenced by the temperature variations.

Mantle Convection Patterns Reveal the Mechanism of the Red Sea Rifting

AbstractWe present a new model of the stress state and present‐day tectonics of the Red Sea Rift (RSR) based on an instantaneous geodynamic mantle flow model. The initial density and viscosity variations in the mantle are derived from a joint inversion of gravity, residual topography, and tomography, which provides higher resolution than existing models. The calculated mantle flow shows clear distinctions along the rift axis. The tectonics of the southern part of the Red Sea is mainly controlled by the Afar plume and characterized by divergent mantle flow. The passive rifting along the central part of the RSR can be explained either by asthenospheric upwelling due to the Red Sea floor spreading or by the plume, rising from the transition zone and not directly related to the Afar plume. We also observed ridge‐axis‐aligned flow in the uppermost mantle in the northern part of the RSR.

Numerical Modeling of Mantle Flow Beneath Madagascar to Constrain Upper Mantle Rheology Beneath Continental Regions

AbstractOver the past few decades, azimuthal seismic anisotropy measurements have been widely used proxy to study past and present‐day deformation of the lithosphere and to characterize convection in the mantle. Beneath continental regions, distinguishing between shallow and deep sources of anisotropy remains difficult due to poor depth constraints of measurements and a lack of regional‐scale geodynamic modeling. Here, we constrain the sources of seismic anisotropy beneath Madagascar where a complex pattern cannot be explained by a single process such as absolute plate motion, global mantle flow, or geology. We test the hypotheses that either Edge‐Driven Convection (EDC) or mantle flow derived from mantle wind interactions with lithospheric topography is the dominant source of anisotropy beneath Madagascar. We, therefore, simulate two sets of mantle convection models using regional‐scale 3‐D computational modeling. We then calculate Lattice Preferred Orientation that develops along pathlines of the mantle flow models and use them to calculate synthetic splitting parameters. Comparison of predicted with observed seismic anisotropy shows a good fit in northern and southern Madagascar for the EDC model, but the mantle wind case only fits well in northern Madagascar. This result suggests the dominant control of the measured anisotropy may be from EDC, but the role of localized fossil anisotropy in narrow shear zones cannot be ruled out in southern Madagascar. Our results suggest that the asthenosphere beneath northern and southern Madagascar is dominated by dislocation creep. Dislocation creep rheology may be dominant in the upper asthenosphere beneath other regions of continental lithosphere.

Erratum: Present-dayd dynamic topography and lower mantle structure from paleogeographically-constrained mantle flow models

Erratum: Present-dayd dynamic topography and lower mantle structure from paleogeographically-constrained mantle flow models

Present-day dynamic topography and lower-mantle structure from palaeogeographically constrained mantle flow models

Present-day dynamic topography and lower-mantle structure from palaeogeographically constrained mantle flow models

Tectonics and geodynamics of the eastern Tethys and northern Gondwana since the Jurassic

Southeast Asia experienced a complex tectonic and geodynamic history related to the subduction of the eastern Tethyan ocean basins, resulting from the long-term convergence between the Indo-Australian, Eurasian, and Pacific plates since Pangea breakup. The complex collage of continental and island arc terranes can be reconstructed into an estimated ancient arrangement using plate tectonic reconstruction approaches based on a synthesis of continental and marine geological and geophysical data. We use the open-source and cross-platform software GPlates (www.gplates.org) to refine the evolution of the eastern Neo-Tethys since the latest Jurassic rifting episodes along northern Gondwana. We apply the resulting plate motions to drive numerical models of mantle flow in order to predict the evolving mantle structure. New Guinea’s northward motion over subducted slabs, related to the Sepik back-arc basin and the Maramuni subduction system, resulted in long-term flooding of the margin since ~20 Ma, despite falling long-term global sea levels. The Sundaland continental promontory experienced dynamic uplift in the latest Cretaceous to Eocene times due to the accretion of the Woyla Arc at ~80 Ma, leading to slab breakoff and a temporary interruption of subduction. However, renewed subduction along the Sunda margin resulted in renewed dynamic subsidence from ~30 Ma, which was amplified by regional basin rifting events. In addition, the sinking Sunda slab likely triggered a mantle slab avalanche, resulting in a counterintuitive combination of contemporaneous basin inversion and strong dynamic subsidence from ~15 Ma. The evolution of the eastern Tethyan oceanic gateway provides an important framework for understanding the role of plate tectonics in controlling long-term oceanic circulation and climate, as well as shedding light on the complex interplay between deep Earth and surface processes in driving basin formation and evolution. These results provide new avenues for reconciling stratigraphic and tectonic processes, as well as contributing new approaches for basin analysis and hydrocarbon exploration.

The adjoint equations for thermochemical compressible mantle convection: derivation and verification by twin experiments.

The adjoint method is an efficient way to obtain gradient information in a mantle convection model relative to past flow structure, allowing one to retrodict mantle flow from observations of the present-day mantle state. While adjoint equations for isochemical mantle flow have been derived for both incompressible and compressible flows, here we extend the method to thermochemical mantle flow models, and present thermochemical adjoint equations in the elastic-liquid approximation. We verify the method with twin experiments, and retrodict the flow history of a thermochemical reference model (reference twin) assuming for the final state, either a consistent thermochemical interpretation, using the thermochemical adjoint equations, or an inconsistent purely thermal interpretation, using the isochemical adjoint equations. The consistent simulation correctly retrodicts the flow evolution of the reference twin. The inconsistent case, instead, restores a false flow history whereby internal buoyancy forces and convectively maintained topography are overestimated. Because the cost function is reduced in either case, our results suggest that the adjoint method can be used to link assumptions on the role of chemical mantle heterogeneity to geologic inferences of dynamic topography, thus providing additional means to test hypotheses on mantle composition and dynamics.

Tectonic and Geodynamic Evolution of the Northern Australian Margin and New Guinea

Rapid convergence between the Indo-Australian, Southeast Asian, and Pacific plates in the Cenozoic has resulted in a complex tectonic evolution of Australia’s northern margin. A lack of available geologic data leads to large uncertainties, such as the timing of the Sepik collision with the New Guinea margin, currently constrained to sometime between 50 and 30 Ma. Previous work suggested a link between the Sepik collision and a voluminous fast seismic anomaly presently in the mantle beneath Lake Eyre. Following from previous work, this study uses coupled plate reconstruction and numerical geodynamic models to test 50 Ma and 30 Ma collision timings of the Sepik terrane, along with an upper extent back-arc basin, to further refine our understanding of the origin and trajectory of the slab beneath Lake Eyre and address uncertainties in the plate reconstructions. The results of mantle flow models indicate that the Eocene collision timing (~50 Ma) is more likely than an Oligocene collision (~30 Ma). In addition, dynamic topography results support previous suggestions that dynamic subsidence relating to the down-going Sepik slab has influenced the evolution of the Eyre Basin, with up to ~100 m of dynamic subsidence since ~20 Ma. However, further work is required to address numerical issues relating to rapid thermal diffusion of slab material, and to investigate reasonable trench retreat velocities for intermediate (~3000 km) subduction zone lengths. This work highlights the role of numerical experiments in understanding transient plate-mantle processes and their effect on basin evolution.

On the Scales of Dynamic Topography in Whole‐Mantle Convection Models

AbstractMantle convection shapes Earth's surface by generating dynamic topography. Observational constraints and regional convection models suggest that surface topography could be sensitive to mantle flow for wavelengths as short as 1,000 and 250 km, respectively. At these spatial scales, surface processes including sedimentation and relative sea‐level change occur on million‐year timescales. However, time‐dependent global mantle flow models do not predict small‐scale dynamic topography yet. Here we present 2‐D spherical annulus numerical models of mantle convection with large radial and lateral viscosity contrasts. We first identify the range of Rayleigh number, internal heat production rate and yield stress for which models generate plate‐like behavior, surface heat flow, surface velocities, and topography distribution comparable to Earth's. These models produce both whole‐mantle convection and small‐scale convection in the upper mantle, which results in small‐scale (<500 km) to large‐scale (>104 km) dynamic topography, with a spectral power for intermediate scales (500 to 104 km) comparable to estimates of present‐day residual topography. Timescales of convection and the associated dynamic topography vary from five to several hundreds of millions of years. For a Rayleigh number of 107, we investigate how lithosphere yield stress variations (10–50 MPa) and the presence of deep thermochemical heterogeneities favor small‐scale (200–500 km) and intermediate‐scale (500–104 km) dynamic topography by controlling the formation of small‐scale convection and the number and distribution of subduction zones, respectively. The interplay between mantle convection and lithosphere dynamics generates a complex spatial and temporal pattern of dynamic topography consistent with constraints for Earth.

Predicting Rates and Distribution of Carbonate Melting in Oceanic Upper Mantle: Implications for Seismic Structure and Global Carbon Cycling

AbstractCoupling a global mantle flow model with a parameterized carbonate solidus, we examine the global distribution and extent of carbonate melting beneath the ocean basins. We predict carbonate melting in spatially heterogeneous patterns throughout the oceans. The rate of CO2 segregation from the mantle by off‐axis carbonate melting (~1.1 × 1012 mol/year) is comparable to the global ridge flux (~1.2 × 1012 mol/year). As the generation of carbonate melts should be enhanced in regions of mantle upwelling, we compare upwelling patterns with seismic detections of the G‐discontinuity. The upwelling velocities in areas where the G‐discontinuity is detected are not statistically different from areas with no detections—implying that the generation and pooling of carbonate melts are not dominant mechanisms in forming the G‐discontinuity. However, detections of a deeper seismic discontinuity are correlated with enhanced upwelling velocities, suggesting locally higher melt fractions at the transition between carbonate‐only and carbonate‐enhanced silicate melting.