Showing papers in "Geochemistry Geophysics Geosystems in 2011"

Geologically current motion of 56 plates relative to the no-net-rotation reference frame

Abstract: NNR-MORVEL56, which is a set of angular velocities of 56 plates relative to the unique reference frame in which there is no net rotation of the lithosphere, is determined. The relative angular velocities of 25 plates constitute the MORVEL set of geologically current relative plate angular velocities; the relative angular velocities of the other 31 plates are adapted from Bird (2003). NNR-MORVEL, a set of angular velocities of the 25 MORVEL plates relative to the no-net rotation reference frame, is also determined. Incorporating the 31 plates from Bird (2003), which constitute 2.8% of Earth's surface, changes the angular velocities of the MORVEL plates in the no-net-rotation frame only insignificantly, but provides a more complete description of globally distributed deformation and strain rate. NNR-MORVEL56 differs significantly from, and improves upon, NNR-NUVEL1A, our prior set of angular velocities of the plates relative to the no-net-rotation reference frame, partly due to differences in angular velocity at two essential links of the MORVEL plate circuit, Antarctica-Pacific and Nubia-Antarctica, and partly due to differences in the angular velocities of the Philippine Sea, Nazca, and Cocos plates relative to the Pacific plate. For example, the NNR-MORVEL56 Pacific angular velocity differs from the NNR-NUVEL1A angular velocity by a vector of length 0.039 ± 0.011° a−1 (95% confidence limits), resulting in a root-mean-square difference in velocity of 2.8 mm a−1. All 56 plates in NNR-MORVEL56 move significantly relative to the no-net-rotation reference frame with rotation rates ranging from 0.107° a−1 to 51.569° a−1.

Petrolog3: Integrated software for modeling crystallization processes

Abstract: This paper introduces Petrolog3, software for modeling (1) fractional and equilibrium crystallization, (2) reverse fractional crystallization at variable pressure, melt oxidation state and melt H2O contents, and (3) postentrapment reequilibration of melt inclusions in olivine. Petrolog3 offers an algorithm that allows calculations with a potentially unlimited number of (1) mineral-melt equilibrium models for major and trace elements and (2) models describing melt physical parameters such as density and viscosity, melt oxidation state, and solubility of fluid components in silicate melts. The current version of the software incorporates 46 mineral-melt equilibrium models for 8 minerals; 3 models describing distribution of trace elements between minerals and melt; 4 models of melt oxidation state; 1 model for H2O solubility in silicate melts; and 4 models describing melt density and viscosity. The idea behind the program is to provide the community of igneous petrologists and geochemists with a user-friendly interface for using any combinations of available mineral-melt equilibrium models for computer simulation of the crystallization process.

An algorithm for U-Pb isotope dilution data reduction and uncertainty propagation

Abstract: This is the publisher's version, also available electronically from "http://onlinelibrary.wiley.com".

Grain size control of river suspended sediment geochemistry: Clues from Amazon River depth profiles

Abstract: Residual solid products of erosion display a wide range of size, density, shape, mineralogy, and chemical composition and are hydrodynamically sorted in large river channels during their transport. We characterize the chemical and isotopic variability of river sediments of the Amazon Basin, collected at different water depths, as a function of grain size. Absolute chemical concentrations and Sr and Nd isotopic ratios greatly varies along channel depth. The Al/Si ratio, tightly linked to grain size distribution, systematically decreases with depth, mostly reflecting dilution by quartz minerals. A double‐normalization diagram is proposed to correct from dilution effects. Elements define fan‐shaped patterns and can be classified in three different groups with respect to hydrodynamic sorting during transport in the Amazon: (1) " poorly sorted " insoluble elements like Al, Fe, Th, and REEs, (2) " well‐sorted " insoluble elements like Zr and Ti, mainly carried by heavy minerals, and (3) alkali (Na to Cs) and alkali‐earth elements (Mg to Ba), for which a large variety of patterns is observed, related, for alkali, to their variable affinity for phyllosilicates. Sr isotopes show that the Amazon River at the mouth is stratified, the Madeira‐ and Solimoes‐derived sediments being preferentially transported near the channel surface and at depth, respectively. The comparison between the Solimoes and Madeira rivers shows how the interplay between grain sorting, weathering, and crustal composition controls the composition of the suspended river sediments.

Engineering cyber infrastructure for U‐Pb geochronology: Tripoli and U‐Pb_Redux

Abstract: This is the publisher's version, also available electronically from "http://onlinelibrary.wiley.com".

Evaluation of precipitation data sets along the Himalayan front

Abstract: [1] Precipitation is one of the main factors which controls surface processes and landscape morphology. Large orogenic belts such as the Himalayas control precipitation distribution as a result of orographic effects due to their prominent relief. However, precipitation is difficult to monitor because mountain regions are largely inaccessible and therefore not sufficiently covered by ground‐based gauge stations. The complexity of orographic effects resulting from the interaction between elevation and climatic processes and the lack of precise meteorological data thus limit our understanding of climatic influence on landscape formation. Therefore, high‐quality precipitation observations with good spatiotemporal coverage are needed. Here we evaluate five gridded precipitation data sets derived from remote sensing and interpolation of rain gauge data with ground‐based precipitation measurements. First, we evaluate the bulk error of each data set, then we evaluate the temporal quality of data within five watersheds, and last we compare the spatial performance along seven swath profiles across strike to the Himalayan range in Nepal. Our evaluation shows that the data sets vary significantly along the orographic front and get more consistent toward the adjacent low‐relief domains, while bulk errors are largest during monsoon season. In particular, where topographic gradients are important, the resolution of gridded data sets cannot incorporate small‐scale spatial changes of precipitation. We show that the data set derived from interpolation of gauge data performs best in the Himalayas. This study gives an overview on the applicability of precipitation data sets within the Himalayan orographic domains where relief has a pronounced impact on precipitation.

Mantle provinces under North America from multifrequency P wave tomography

Abstract: This is a survey of mantle provinces (large-scale seismic anomalies) under North America, from the surface down to 1500–1800 km depth The underlying P velocity model was obtained by multifrequency tomography, a waveform-based method that systematically measures and models the frequency-dependence of teleseismic body waves A novel kind of three-dimensional rendering technique is used to make the considerable structural complexities under North America accessible In the transition zone and below, the North American mantle is dominated by seismically fast provinces, which represent distinct subduction episodes of the Farallon plate I attempt to date and interpret the various slab fragments by reconciling their present positions with paleotrench locations from plate tectonic reconstructions and with major geologic surface episodes Differences in vertical sinking velocity have led to large vertical offsets across adjacent, coeval slabs Some of the mantle provinces have not been discussed much previously, including (1) a seismically slow blanket overlying the oldest Farallon subduction along the eastern continental margin, (2) a transition zone slab coeval with the Laramide orogeny (ca 80–60 Myr), which I discuss in analogy to the “stagnant slab” subduction style commonly found in the western Pacific today, (3) the lower mantle root of present-day Cascadia subduction, which may have started out as intraoceanic subduction,(4) a lower mantle slab under Arizona and New Mexico, the last material to subduct before strike-slip motion developed along the San Andreas boundary, and (5) two narrow plate tears thousands of kilometers long, one of which is the subducted conjugate of the Mendocino Fracture Zone

The nature of the crust beneath the Afar triple junction: Evidence from receiver functions

Abstract: The Afar depression is an ideal locale to study the role of extension and magmatism as rifting progresses to seafloor spreading. Here we present receiver function results from new and legacy experiments. Crustal thickness ranges from ∼45 km beneath the highlands to ∼16 km beneath an incipient oceanic spreading center in northern Afar. The crust beneath Afar has a thickness of 20–26 km outside the currently active rift segments and thins northward. It is bounded by thick crust beneath the highlands of the western plateau (∼40 km) and southeastern plateau (∼35 km). The western plateau shows VP/VS ranging between 1.7–1.9, suggesting a mafic altered crust, likely associated with Cenozoic flood basalts, or current magmatism. The southeastern plateau shows VP/VS more typical of silicic continental crust (∼1.78). For crustal thicknesses 2.0) can only be explained by significant amounts of magmatic intrusions in the lower crust. This suggests that melt emplacement plays an important role in late stage rifting, and melt in the lower crust likely feeds magmatic activity. The crust between the location of the Miocene Red Sea rift axis and the current rift axis is thinner ( 2.0) than beneath the eastern part of Afar (>26 km, VP/VS < 1.9). This suggests that the eastern region contains less partial melt, has undergone less stretching/extension and has preserved a more continental crustal signature than west of the current rift axis. The Red Sea rift axis appears to have migrated eastward through time to accommodate the migration of the Afar triple junction.

Lateral constrictional flow of hot orogenic crust: Insights from the Neoarchean of south India, geological and geophysical implications for orogenic plateaux

Abstract: This study provides an in situ geological perspective on fabrics produced by synconvergence lateral crustal flow of hot orogens. It is based on the example of the Neoarchean orogen of the Dharwar craton (India) and combines structural analysis and ion microprobe U-Pb zircon geochronology. We document a pervasive, three-dimensional flow mode of the lower crust, called lateral constrictional flow (LCF), which combines orogen-normal shortening, lateral constrictional stretching, and transtension. LCF achieves gravity-driven flow, lateral escape, and 3-D mass redistribution in a viscous lower crust submitted to convergence. LCF tends to mechanically and thermally homogenize the lower crust and efficiently compensates topographic relief at a shallow level in the crust. Three type-geodynamic contexts are envisaged for LCF: plateau interiors, inner parts of collisional crustal wedges or plateau edges, and throughout wide ultrahot orogens such as the Neaoarchean orogen of south India. LCF makes the lower crust act as a strain gauge between lateral gravitational collapse or tectonic thickening of the upper crust, thrust stacking in the lowermost crust (collisional crustal wedge case), crustal shortening, and/or lateral flow of the upper mantle. In the case of plateau interiors or ultrahot orogens, LCF of a thick lower crust enables coupling of upper crustal deformation with upper mantle flow through a hot and thin lithosphere being shortened coherently. LCF generates a subhorizontal lamination that should produce the strong seismic reflectivity and lateral anisotropy of the Tibetan Plateau and a variety of hot orogens.

Domains of depleted mantle: New evidence from hafnium and neodymium isotopes

Abstract: Isotope systematics of basalts provide information on the distribution of mantle components and the length scale of mantle heterogeneity. To obtain this information, high data and sampling density are crucial. We present hafnium and neodymium isotope data on more than 400 oceanic volcanics. Over length scales of several hundred to over one thousand kilometers hafnium and neodymium isotopes of mid-ocean ridge basalts are correlated and form an array of parallel trends on a global scale. On a larger scale these domains differ in the amount of highly depleted mantle material with radiogenic hafnium and neodymium isotope ratios. Compared to the Atlantic and Indian Ocean basins the asthenosphere of the Pacific basin seems to have a more uniform and a less radiogenic Hf isotopic composition for a given Nd isotopic composition. The parallel arrays of mid-ocean ridge basalts provide strong constraints on the makeup of the MORB mantle and are evidence for the presence of a highly depleted and highly radiogenic neodymium and hafnium component. This component, because of its highly depleted character, is unrecognized in the strontium-neodymium-lead isotope systems alone. Alternatively, the parallel arrays can have an ancient origin of systematic variations in the degree of depletion. Each array then represents the variations in this fossil melting regime. Individual ocean island basalt suites display different slopes in hafnium-neodymium isotope space, which are also best explained by varying amounts of highly residual mantle rather than isotopic differences in enriched mantle components as previously invoked. The ocean island basalt arrays diverge at the depleted end and project to radiogenic compositions that are similar to those of the asthenosphere through which they travel. This is strong evidence that the plume material interacts with its surrounding mantle as it ascends. The isotopic compositions of the ocean island and ridge basalts suggest that their systematics are influenced by a heretofore unrecognized depleted component.

Physical characteristics of subduction interface type seismogenic zones revisited

Abstract: [1] Based on global earthquake catalogs, the hypocenters, nodal planes, and seismic moments of worldwide subduction plate interface earthquakes were extracted for the period between 1900 and 2007. Assuming that the seismogenic zone coincides with the distribution of 5.5 ≤ M < 7 earthquakes, the subduction interface seismogenic zones were mapped for 80% of the trench systems and characterized with geometrical and mechanical parameters. Using this database, correlations were isolated between significant parameters to identify cause-effect relationships. Empirical laws obtained in previous studies were revisited in light of this more complete, accurate, and uniform description of the subduction interface seismogenic zone. The seismogenic zone was usually found to end in a fore-arc mantle, rather than at a Moho depth. The subduction velocity was the first-order controlling parameter for variations in the physical characteristics of plate interfaces, determining both the geometry and mechanical behavior. As such, the fast subduction zones and cold slabs were associated with large and steep plate interfaces, which, in turn, had large seismic rates. The subduction velocity could not account for the potential earthquake magnitude diversity that was observed along the trenches. Events with Mw ≥ 8.5 preferentially occurred in the vicinity of slab edges, where the upper plate was continental and the back-arc strain was neutral. This observation was interpreted in terms of compressive normal stresses along the plate interface. Large lateral ruptures should be promoted in neutral subduction zones due to moderate compressive stresses along the plate interface that allow the rupture to propagate laterally.

Slab segmentation and late Cenozoic disruption of the Hellenic arc

Abstract: The Hellenic subduction zone displays well-defined temporal and spatial variations in subduction rate and offers an excellent natural laboratory for studying the interaction among slab buoyancy, subduction rate, and tectonic deformation. In space, the active Hellenic subduction front is dextrally offset by 100–120 km across the Kephalonia Transform Zone, coinciding with the junction of a slowly subducting Adriatic continental lithosphere in the north (5–10 mm/yr) and a rapidly subducting Ionian oceanic lithosphere in the south (∼35 mm/yr). Subduction rates can be shown to have decreased from late Eocene time onward, reaching 5–12 mm/yr by late Miocene time, before increasing again along the southern portion of the subduction system. Geodynamic modeling demonstrates that the differing rates of subduction and the resultant trench offset arise naturally from subduction of oceanic (Pindos) lithosphere until late Eocene time, followed by subduction of a broad tract of continental or transitional lithosphere (Hellenic external carbonate platform) and then by Miocene entry of high-density oceanic (Ionian) lithosphere into the southern Hellenic trench. Model results yield an initiation age for the Kephalonia Transform of 6–8 Ma, in good agreement with observations. Consistency between geodynamic model results and geologic observations suggest that the middle Miocene and younger deformation of the Hellenic upper plate, including formation of the Central Hellenic Shear Zone, can be quantitatively understood as the result of spatial variations in the buoyancy of the subducting slab. Using this assumption, we make late Eocene, middle Miocene, and Pliocene reconstructions of the Hellenic system that include quantitative constraints from subduction modeling and geologic constraints on the timing and mode of upper plate deformation.

Structural assessment of Mount Etna volcano from Permanent Scatterers analysis

Abstract: A study of the deformation pattern of Mount Etna volcano based on the results from the Permanent Scatterers (PS) technique is reported. Ground motion data provided by the interferometric synthetic aperture radar (InSAR) PS technique from 1995 to 2000 are compared and validated by GPS data. An analysis of the ascending and descending line of sight (LOS) components of ground velocities has yielded detailed ground deformation maps and cross sections. This analysis allows detection and constraint of discontinuities in the surface velocity field. LOS velocities have then been combined to calculate the vertical and horizontal (E-W) ground velocities. A wide inflation of the edifice has been detected on the western and northern flanks (over an area of about 350 km2). A seaward motion of the eastern and southern flanks has also been measured. PS data allows the geometry and kinematics of the several blocks composing the unstable flanks to be defined even in the highly urbanized areas, and their displacement rates have been measured with millimeter precision. This analysis reveals the extension of some features beyond their field evidences and defines new important features. The results of this work depict a new comprehensive kinematic model of the volcano highlighting the gravitational reorganization of the unbuttressed volcanic pile on its slippery clay basement on the southern flank, but an additional drag force due to a strong subsidence of the continental margin facing the Etna volcano is necessary to explain the PS velocity field observed on the eastern flank.

Grain size control on Sr-Nd isotope provenance studies and impact on paleoclimate reconstructions: An example from deep-sea sediments offshore NW Africa

Abstract: Five sediment cores from offshore NW Africa were analyzed for strontium and neodymium isotope ratios to reconstruct temporal variations in continental weathering regimes. Sediments were taken from three time slices with well-known and distinctive environmental conditions: present-day (dry and warm), similar to 6 ka (wet and warm), and similar to 12 ka (dry and cold). Terrigenous sediment samples were split into two size fractions to distinguish between the two dominant transport mechanisms offshore NW Africa: fluvial (0-10 mu m) and aeolian (10-40 mu m). Sr isotope data record evidence of marked grain size control with higher isotopic ratios in the fine fraction. In contrast, epsilon(Nd) values are largely unaffected by grain size. Minor variability in Nd isotope data at each sampling site indicates near constant sources of terrigenous matter over the last similar to 12 ka. Variations in Sr isotope ratios are interpreted to reflect major changes in the evaporation-precipitation balance. We suggest that the Sr-Nd isotope data record a latitudinal shift of the northern limit of the African rainbelt and associated wind systems causing changes in the humidity and rate of chemical weathering over NW Africa. While hyperarid conditions prevailed similar to 12 ka, more humid conditions and intensified monsoonal rainfall at similar to 6 ka resulted in greater breakdown of easily weathered K-bearing phases and increased Sr-87/Sr-86 in the detritus. In late Holocene times the monsoonal circulation diminished resulting in a return to arid conditions. Our results clearly show that it is of vital importance in paleoenvironmental studies to carry out isotopic analyses on individual sediment fractions that were carried to the studied deposition site by distinct sediment transport mechanisms. If isotopic analyses are carried out on bulk sediments, the observed variability in isotopic values most likely represents changes in the particle size and mixing proportions of the sediment subpopulations.

Structure and serpentinization of the subducting Cocos plate offshore Nicaragua and Costa Rica

Abstract: The Cocos plate experiences extensional faulting as it bends into the Middle American Trench (MAT) west of Nicaragua, which may lead to hydration of the subducting mantle. To estimate the along strike variations of volatile input from the Cocos plate into the subduction zone, we gathered marine seismic refraction data with the R/V Marcus Langseth along a 396 km long trench parallel transect offshore of Nicaragua and Costa Rica. Our inversion of crustal and mantle seismic phases shows two notable features in the deep structure of the Cocos plate: (1) Normal oceanic crust of 6 km thickness from the East Pacific Rise (EPR) lies offshore Nicaragua, but offshore central Costa Rica we find oceanic crust from the northern flank of the Cocos Nazca (CN) spreading center with more complex seismic velocity structure and a thickness of 10 km. We attribute the unusual seismic structure offshore Costa Rica to the midplate volcanism in the vicinity of the Galapagos hot spot. (2) A decrease in Cocos plate mantle seismic velocities from ∼7.9 km/s offshore Nicoya Peninsula to ∼6.9 km/s offshore central Nicaragua correlates well with the northward increase in the degree of crustal faulting outboard of the MAT. The negative seismic velocity anomaly reaches a depth of ∼12 km beneath the Moho offshore Nicaragua, which suggests that larger amounts of water are stored deep in the subducting mantle lithosphere than previously thought. If most of the mantle low velocity zone can be interpreted as serpentinization, the amount of water stored in the Cocos plate offshore central Nicaragua may be about 2.5 times larger than offshore Nicoya Peninsula. Hydration of oceanic lithosphere at deep sea trenches may be the most important mechanism for the transfer of aqueous fluids to volcanic arcs and the deeper mantle.

Compositional trends of Icelandic basalts: Implications for short-length scale lithological heterogeneity in mantle plumes

Abstract: Lithological variations in the mantle source regions under mid-ocean ridges and ocean islands have been proposed to play a key role in controlling melt generation and basalt composition. Here we combine compositional observations from Icelandic basalts and modelling of melting of a bi-lithologic peridotite-pyroxenite mantle to demonstrate that, while short-lengthscale major element variation is present in the mantle under Iceland, source heterogeneity does not make an important contribution to excess melt production. By identifying the major element characteristics of endmember Icelandic melts, we find enriched melts to be characterised by low SiO2 and CaO, but high FeO. We quantitatively compare endmember compositions to experimental partial melts generated from a range of lithologies, pressures and melt fractions. This comparison indicates that a single source composition cannot account for all the major element variation; depleted Icelandic melts can be produced by depleted peridotite melting, but the major element composition of enriched melts is best matched by melting of mantle sources that have been refertilised by the addition of up to 40% mid-ocean ridge basalt. The enriched source beneath Iceland is more fusible than the source of depleted melts, and as such will be over-represented in accumulated melts compared with its abundance in the source. Modelling of peridotite-pyroxenite melting, combined with our observational constraints on the composition of the Icelandic mantle, indicates that crustal thickness variations in the North Atlantic must be primarily due to mantle temperature and flow field variations.

Pulses of deformation reveal frequently recurring shallow magmatic activity beneath the Main Ethiopian Rift

Abstract: Magmatism strongly influences continental rift development, yet the mechanism, distribution, and timescales on which melt is emplaced and erupted through the shallow crust are not well characterized. The Main Ethiopian Rift (MER) has experienced significant volcanism, and the mantle beneath is characterized by high temperatures and partial melt. Despite its magma-rich geological record, only one eruption has been historically recorded, and no dedicated monitoring networks exist. Consequently, the present-day magmatic processes in the region remain poorly documented, and the associated hazards are neglected. We use satellite-based interferometric synthetic aperture radar observations to demonstrate that significant deformation has occurring at four volcanic edifices in the MER (Alutu, Corbetti, Bora, and Haledebi) from 1993 to 2010. This raises the number of volcanoes known to be deforming in East Africa beyond 12, comparable to many subduction arcs despite the smaller number of recorded eruptions. The largest displacements are at Alutu volcano, the site of a geothermal plant, which showed two pulses of rapid inflation (10–15 cm) in 2004 and 2008 separated by gradual subsidence. Our observations indicate a shallow (<10 km), frequently replenished zone of magma storage associated with volcanic edifices and add to the growing body of observations that indicate shallow magmatic processes operating on a decadal timescale are ubiquitous throughout the East African Rift. In the absence of detailed historical records of volcanic activity, satellite-based observations of monitoring parameters, such as deformation, could play an important role in assessing volcanic hazard.

Constraints on the numerical age of the Paleocene‐Eocene boundary

Abstract: Here we present combined radioisotopic dating (U-Pb zircon) and cyclostratigraphic analysis of the carbon isotope excursion at the Paleocene-Eocene (P-E) boundary in Spitsbergen to determine the numerical age of the boundary. Incorporating the total uncertainty from both radioisotopic and cyclostratigraphic data sets gives an age ranging from 55.728 to 55.964 Ma, within error of a recently proposed astronomical age of ∼55.93 Ma. Combined with the assumption that the Paleocene Epoch spans twenty-five 405 kyr cycles, our new age for the boundary suggests an age of ∼66 Ma for the Cretaceous-Paleogene boundary. Furthermore, our P-E boundary age is consistent with the hypothesis that the onset of the Paleocene-Eocene thermal maximum at the boundary occurred on the falling limb of a 405 kyr cycle, suggesting the event was initiated by a different mechanism to that which triggered the other early Eocene hyperthermals.

Making and breaking an island arc: A new perspective from the Oligocene Kyushu-Palau arc, Philippine Sea

Abstract: The Kyushu-Palau Ridge (KPR) is a 2600 km long remnant island arc that is separated from the active Izu-Bonin-Mariana (IBM) arc by a series of spreading and rift basins. We present 40Ar/39Ar ages and geochemical data for the entire length of the Kyushu-Palau arc as well as for the conjugate arc which is stranded within the IBM fore arc. New 40Ar/39Ar ages indicate that the KPR was active between 25 and 48 Ma, but the majority of the exposed volcanism occurred in the final phase, between 25 and 28 Ma. Rifting of the Kyushu-Palau arc to form the Shikoku and Parece Vela basins occurred simultaneously along the length of the arc (circa 25 Ma), and at a similar distance from the trench. Unlike the IBM, the KPR has only limited systematic along-arc geochemical trends. Two geochemical components within the KPR indicate an origin in the suprasubduction mantle. First, EM-1-like lavas are identified in a restricted section of the arc, suggesting a localized heterogeneity. Second, EM-2-like arc volcanoes formed on juvenile West Philippine Basin crust, potentially reflecting ingress of mantle from the then active EM-2 province which lies in the west. Another geochemical heterogeneity is found at the KPR–Daito Ridge intersection where the arc developed on preexisting Cretaceous Daito Ridge crust. The geochemical characteristics at this intersection likely result from the involvement of sub–Daito Ridge lithospheric mantle. Subduction flux beneath the KPR generally matches post–45 Ma Eocene/Oligocene lavas in the IBM fore arc, involving fluids and melts derived from altered igneous crust.

On the location of plumes and lateral movement of thermochemical structures with high bulk modulus in the 3-D compressible mantle

Abstract: The two large low shear velocity provinces (LLSVPs) at the base of the lower mantle are prominent features in all shear wave tomography models. Various lines of evidence suggest that the LLSVPs are thermochemical and are stable on the order of hundreds of million years. Hot spots and large igneous province eruption sites tend to cluster around the edges of LLSVPs. With 3-D global spherical dynamic models, we investigate the location of plumes and lateral movement of chemical structures, which are composed of dense, high bulk modulus material. With reasonable values of bulk modulus and density anomalies, we find that the anomalous material forms dome-like structures with steep edges, which can survive for billions of years before being entrained. We find that more plumes occur near the edges, rather than on top, of the chemical domes. Moreover, plumes near the edges of domes have higher temperatures than those atop the domes. We find that the location of the downwelling region (subduction) controls the direction and speed of the lateral movement of domes. Domes tend to move away from subduction zones. The domes could remain relatively stationary when distant from subduction but would migrate rapidly when a new subduction zone initiates above. Generally, we find that a segment of a dome edge can be stationary for 200 million years, while other segments have rapid lateral movement. In the presence of time-dependent subduction, the computations suggest that maintaining the lateral fixity of the LLSVPs at the core-mantle boundary for longer than hundreds of million years is a challenge.

Fluidity: A fully unstructured anisotropic adaptive mesh computational modeling framework for geodynamics

Abstract: We present a new computational modeling framework, Fluidity, for application to a range of two- and three-dimensional geodynamic problems, with the focus here on mantle convection. The approach centers upon a finite element discretization on unstructured simplex meshes, which represent complex geometries in a straightforward manner. Throughout a simulation, the mesh is dynamically adapted to optimize the representation of evolving solution structures. The adaptive algorithm makes use of anisotropic measures of solution complexity, to vary resolution and allow long, thin elements to align with features such as boundary layers. The modeling framework presented differs from the majority of current mantle convection codes, which are typically based upon fixed structured grids. This necessitates a thorough and detailed validation, which is a focus of this paper. Benchmark comparisons are undertaken with a range of two- and three-dimensional, isoviscous and variable viscosity cases. In addition, model predictions are compared to experimental results. Such comparisons highlight not only the robustness and accuracy of Fluidity but also the advantages of anisotropic adaptive unstructured meshes, significantly reducing computational requirements when compared to a fixed mesh simulation.

Imaging the seismic lithosphere-asthenosphere boundary of the oceanic plate

Abstract: [1] The seismic lithosphere-asthenosphere boundary (LAB) or G discontinuity, a seismologically characterized abrupt drop in wave speed in the uppermost mantle, is one of the key issues in current geodynamics. Although plate tectonics started as a theory for the ocean, reports on LAB for normal oceanic regions are scarce due to paucity of seismic data, and whether or not the oceanic LAB grows with age is the key issue to be resolved. We conduct a systematic survey for the oceanic LAB using S-to-p converted seismic waves along three margins of oceanic plates whose crustal age ranges from ∼10 Myr to ∼130 Myr, and we observe laterally continuous oceanic LAB images. The thickness of the oceanic plate estimated from LABs increases with the plate age, though scattered, suggesting that the evolution of oceanic lithosphere is predominantly governed by temperature and that the oceanic seismic LAB represents a boundary that grows with age.

Geochemical characteristics and origin of the HIMU reservoir: A possible mantle plume source in the lower mantle

Abstract: Combined Pb-Sr-Nd-Hf-Os isotopes, together with major and trace element compositions, were determined from clinopyroxene and olivine phenocrysts, along with whole rocks, for ocean island basalts with high μ (μ = 238U/204Pb) (HIMU) and enriched mantle isotopic characteristics from Cook-Austral Islands Clinopyroxene and olivine separates record reliable isotopic information of the sources because of minimized in situ radiogenic ingrowth and their lower susceptibility to crustal contamination Coherent isotopic systematics in multi-isotope spaces defined by the HIMU samples are best explained by recent mixing of melts derived from the HIMU reservoir and the local shallow mantle The isotopic compositions of the HIMU reservoir are constrained to be low ɛNd (≤+4), low ɛHf (≤+3), and moderately radiogenic 187Os/188Os (014–015) in association with radiogenic Pb isotopes (206Pb/204Pb ≥ 215) Since ancient oceanic crust would have had exceptionally radiogenic 187Os/188Os, moderately high 187Os/188Os precludes recycled oceanic crust as the only contributor to the HIMU reservoir Instead, mantle metasomatized with partial melts from subducted oceanic crust is a likely candidate for the HIMU reservoir Moreover, partial melting of oceanic crust in equilibrium with Mg perovskite would fractionate U/Pb, Sm/Nd, and Lu/Hf, which are in accordance with the time-integrated U/Pb, Sm/Nd, and Lu/Hf deduced from Pb, Nd, and Hf isotopic compositions of the HIMU reservoir, respectively, with a formation age of 2–3 Ga We thus propose that the HIMU reservoir was formed by hybridization of a subducted oceanic crust-derived melt with the ambient mantle and then stored for several billion years in the lower mantle

Discretization errors and free surface stabilization in the finite difference and marker-in-cell method for applied geodynamics: A numerical study

Abstract: The finite difference–marker-in-cell (FD-MIC) method is a popular method in thermomechanical modeling in geodynamics. Although no systematic study has investigated the numerical properties of the method, numerous applications have shown its robustness and flexibility for the study of large viscous deformations. The model setups used in geodynamics often involve large smooth variations of viscosity (e.g., temperature-dependent viscosity) as well large discontinuous variations in material properties (e.g., material interfaces). Establishing the numerical properties of the FD-MIC and showing that the scheme is convergent adds relevance to the applications studies that employ this method. In this study, we numerically investigate the discretization errors and order of accuracy of the velocity and pressure solution obtained from the FD-MIC scheme using two-dimensional analytic solutions. We show that, depending on which type of boundary condition is used, the FD-MIC scheme is a second-order accurate in space as long as the viscosity field is constant or smooth (i.e., continuous). With the introduction of a discontinuous viscosity field characterized by a viscosity jump (η*) within the control volume, the scheme becomes first-order accurate. We observed that the transition from second-order to first-order accuracy will occur with only a small increase in the viscosity contrast (η* ≈ 5). We have employed two methods for projecting the material properties from the Lagrangian markers onto the Eulerian nodes. The methods are based on the size of the interpolation volume (4-cell, 1-cell). The use of a more local interpolation scheme (1-cell) decreases the absolute velocity and pressure discretization errors. We also introduce a stabilization algorithm that damps the potential oscillations that may arise from quasi free surface calculations in numerical codes that employ the strong form of the Stokes equations. This correction term is of particular interest for topographic modeling, since the surface of the Earth is generally represented by a free surface. Including the stabilization enables physically meaningful solutions to be obtained from our simulations, even in cases where the time step value exceeds the isostatic relaxation time. We show that including the stabilization algorithm in our FD stencil does not affect the convergence properties of our scheme. In order to verify our approach, we performed time-dependent simulations of free surface Rayleigh-Taylor instability.

Rheology of magmas with bimodal crystal size and shape distributions: Insights from analog experiments

Abstract: [1] Magmas in volcanic conduits commonly contain microlites in association with preexisting phenocrysts, as often indicated by volcanic rock textures. In this study, we present two different experiments that investigate the flow behavior of these bidisperse systems. In the first experiments, rotational rheometric methods are used to determine the rheology of monodisperse and polydisperse suspensions consisting of smaller, prolate particles (microlites) and larger, equant particles (phenocrysts) in a bubble-free Newtonian liquid (silicate melt). Our data show that increasing the relative proportion of prolate microlites to equant phenocrysts in a magma at constant total particle content can increase the relative viscosity by up to three orders of magnitude. Consequently, the rheological effect of particles in magmas cannot be modeled by assuming a monodisperse population of particles. We propose a new model that uses interpolated parameters based on the relative proportions of small and large particles and produces a considerably improved fit to the data than earlier models. In a second series of experiments we investigate the textures produced by shearing bimodal suspensions in gradually solidifying epoxy resin in a concentric cylinder setup. The resulting textures show the prolate particles are aligned with the flow lines and spherical particles are found in well-organized strings, with sphere-depleted shear bands in high-shear regions. These observations may explain the measured variation in the shear thinning and yield stress behavior with increasing solid fraction and particle aspect ratio. The implications for magma flow are discussed, and rheological results and textural observations are compared with observations on natural samples.

Polarized Earth's ambient microseismic noise

Abstract: We quantify, analyze, and characterize the frequency-dependent microseismic noise recorded by worldwide distributed seismic stations. Microseismic noise is generated through the interaction of ocean waves. It is the strongest ambient noise, and it is observed everywhere on Earth. We introduce a new approach which permits us to detect polarized signals in the time-frequency domain and which we use to characterize the microseismic noise. We analyze 7 years of continuous seismograms from the global GEOSCOPE network. Microseisms are dominated by Rayleigh waves, and we therefore focus on elliptically polarized signals. The polarized signals are detected in the time-frequency domain through a degree of polarization measure. We design polarization spectra and show that microseismic noise is more strongly polarized than noise in other frequency bands. This property is used to measure the directions of the polarized noise at individual stations as a function of time and frequency. Seasonal variations are found for the back azimuths and for the number of polarized signals at many stations. We show that the back azimuth directions are robust measurements that point toward the source areas computed from ocean wave models.

The 3‐D geometry of detachment faulting at mid‐ocean ridges

Abstract: Seismic images of Cretaceous slow spreading crust from the eastern Central Atlantic provide new constraints on the process of seafloor spreading and the importance of detachment faulting. The seismic depth sections image detachment faults that appear to have exhumed footwall massifs of similar geometry to massifs of plutonic and mantle rocks mapped at the present Mid-Atlantic Ridge. The detachments are consistent with the structure and microearthquakes of the Mid-Atlantic Ridge at 26°N, with the footwall rotation inferred from paleomagnetic data and with numerical modeling of oceanic detachments. Other seismically imaged detachments have similar dimensions and geometry, but are covered by a layer of small fault blocks. The detachment types differ in whether or not the fault locks up in the subsurface, probably controlled by the fault strength, the elastic thickness, and whether the exhumed footwall is partly covered by basalts. Toward the segment middle, decreasing mantle serpentinization, decreasing elastic thickness, and thicker median valley basalts all increase the likelihood that the fault locks up, and a new fault propagates upwards from the still active root zone, transferring a slice of the hanging wall to the footwall, to be rafted with the footwall out of the median valley. As a result an oceanic detachment fault, exhuming the footwall at a segment end to form an oceanic core complex, may disappear laterally beneath rafted blocks; detachment faulting may be more widespread at slow spreading ridges than interpreted from seafloor mapping.

Three‐dimensional Moho topography in Italy: New constraints from receiver functions and controlled source seismology

Abstract: In complex tectonics regions, seismological, geophysical, and geodynamic modeling require accurate definition of the Moho geometry. Various active and passive seismic experiments performed in the central Mediterranean region revealed local information on the Moho depth, in some cases used to produce interpolated maps. In this paper, we present a new and original map of the 3-D Moho geometry obtained by integrating selected high-quality controlled source seismic and teleseismic receiver function data. The very small cell size makes the retrieved model suitable for detailed regional studies, crustal corrections in teleseismic tomography, advanced 3-D ray tracing in regional earthquake location, and local earthquake tomography. Our results show the geometry of three different Moho interfaces: the European, Adriatic-Ionian, and Tyrrhenian. The three distinct Moho are fashioned following the Alpine and Apennines subduction, collision, and back-arc spreading and show medium- to high-frequency topographic undulations reflecting the complexity of the geodynamic evolution.

Continents, supercontinents, mantle thermal mixing, and mantle thermal isolation: Theory, numerical simulations, and laboratory experiments

Abstract: Super-continental insulation refers to an increase in mantle temperature below a supercontinent due to the heat transfer inefficiency of thick, stagnant continental lithosphere relative to thinner, subducting oceanic lithosphere. We use thermal network theory, numerical simulations, and laboratory experiments to provide tighter physical insight into this process. We isolate two end-member dynamic regimes. In the thermally well mixed regime the insulating effect of continental lithosphere can not cause a localized increase in mantle temperature due to the efficiency of lateral mixing in the mantle. In this regime the potential temperature of the entire mantle is higher than it would be without continents, the magnitude depending on the relative thickness of continental and oceanic lithosphere (i.e., the insulating effects of continental lithosphere are communicated to the entire mantle). Thermal mixing can be short circuited if subduction zones surround a supercontinent or if the convective flow pattern of the mantle becomes spatially fixed relative to a stationary supercontinent. This causes a transition to the thermal isolation regime: The potential temperature increases below a supercontinent whereas the potential temperature below oceanic domains drops such that the average temperature of the whole mantle remains constant. Transition into this regime would thus involve an increase in the suboceanic viscosity, due to local cooling, and consequently a decrease in the rate of oceanic lithosphere overturn. Transition out of this regime can involve the unleashing of flow driven by a large lateral temperature gradient, which will enhance global convective motions. Our analysis highlights that transitions between the two states, in either direction, will effect not only the mantle below a supercontinent but also the mantle below oceanic regions. This provides a larger set of predictions that can be compared to the geologic record to help determine if a hypothesized super-continental thermal effect did or did not occur on our planet.

Slumping and mass transport deposition in the Nankai fore arc: Evidence from IODP drilling and 3‐D reflection seismic data

Abstract: Multiple lines of evidence exist for a range of sediment mass movement processes within the shallow megasplay fault zone (MSFZ) area and the adjacent slope basin in the outer fore arc of the Nankai subduction zone, Japan. Diagnostic features observed in three-dimensional reflection seismic data and in cores of the Integrated Ocean Drilling Program (IODP) document a multifarious mass movement history spanning ∼2.87 million years. Various modes and scales of sediment remobilization can be related to the different morphotectonic settings in which they occurred. From this evidence, we decipher the tectonic control on slumping and mass transport deposition in the Nankai fore arc. Three periods of intensified mass wasting coincided with pulses of enhanced activity on the splay fault: (1) an initial phase of juvenile out-of-sequence thrusting ∼1.95 to 1.7 Ma, (2) a reactivation phase between ∼1.55 and 1.24 Ma, and (3) at about 1 Ma, during a phase of uplift of the fore-arc high and motion along the MSFZ. We suggest that slope oversteepening, extensional stress regimes, and lateral transmission of fluid overpressures may have preconditioned the slope sediments to fail. Individual mass-wasting events may have been triggered by dynamic loading from earthquake waves and/or transient pulses of pore pressure along the splay fault. Overall, our results provide insights into the complicated interplay between tectonic and submarine mass movement processes. We demonstrate that detailed knowledge about the spatial and temporal distribution of submarine mass movements can be integrated into a holistic reconstruction of tectonostratigraphic evolution of accretionary margins.