Preliminary analysis of seismic anisotropy in the uppermost mantle beneath NW Pacific reveled by the Normal Oceanic Mantle project

NASA Astrophysics Data System (ADS)

Takeo, A.; Isse, T.; Nishida, K.; Kawakatsu, H.; Shiobara, H.; Sugioka, H.; Ito, A.; Utada, H.

2013-12-01

Seismic structure including anisotropy in the oceanic uppermost mantle is essential for understanding deformation related to plate tectonics. Recent reports of a sharp discontinuity between the high velocity LID and the low velocity zone (LVZ) especially emphasize the importance of observation in oceanic basins apart from ridges and hotspots for determining the structure including LID and LVZ. In this study, we analyzed records of four broadband ocean bottom seismometers (BBOBSs) deployed in the northwest of Shatsky Rise by the pilot observation of the Normal Oceanic Mantle (NOMan) project in 2010-2011. We first measured average phase velocities of surface waves at periods of 5-30 s by the ambient-noise cross correlation method. Based on the method of Takeo et al. (in prep. GJI), we analyzed fundamental- and first higher- mode Rayleigh waves and fundamental-mode Love wave simultaneously by waveform fitting after the correction of clock delay. At periods of 25-100 s, we measured phase velocities of fundamental-mode surface waves by the array analysis of teleseismic waveforms. We then determined one-dimensional radially anisotropic structure beneath the array by the method of Takeo et al. (2013, JGR). The obtained structure shows transition from LID to LVZ at depths of 50-80km, which is marginally consistent with the depth of ~80 km estimated by a receiver function analysis at WP2 station situated at east of the studies area (Kawakatsu et al., 2009). The velocity gradient in the LID is almost zero and inconsistent with the simple cooling model of homogeneous oceanic plate. The average intensity of S-wave radial anisotropy at depths of ~10-220 km is ~3% (VSH>VSV). We further estimated S-wave azimuthal anisotropy at depths of ~30-100 km by analyzing teleseismic fundamental-mode Rayleigh waves at periods of 25-50 s. The intensity of anisotropy is 2-3%. The fastest direction is about N35W, close to that of Sn-wave velocity around WP2 station obtained by a refraction survey (Shinohara et al., 2008), and indicates the presence of past mantle flow almost perpendicular to the ancient mid ocean ridge or the presence of current mantle flow parallel to the plate motion at depths of 30-100 km. We will further analyze new records after the recovery of 13 BBOBSs in August 2013 and will present more detailed structure around Shatsky Rise. BBOBS stations of pilot observation of NOMan project (white crosses), WP2 station (circle), isochrons (white lines). Black bars show the fastest directions of Rayleigh wave at periods of 25-50 s and the fastest direction of Sn-wave velocity (Shinohara et al. 2008).

Probability density functions for radial anisotropy: implications for the upper 1200 km of the mantle

NASA Astrophysics Data System (ADS)

Beghein, Caroline; Trampert, Jeannot

2004-01-01

The presence of radial anisotropy in the upper mantle, transition zone and top of the lower mantle is investigated by applying a model space search technique to Rayleigh and Love wave phase velocity models. Probability density functions are obtained independently for S-wave anisotropy, P-wave anisotropy, intermediate parameter η, Vp, Vs and density anomalies. The likelihoods for P-wave and S-wave anisotropy beneath continents cannot be explained by a dry olivine-rich upper mantle at depths larger than 220 km. Indeed, while shear-wave anisotropy tends to disappear below 220 km depth in continental areas, P-wave anisotropy is still present but its sign changes compared to the uppermost mantle. This could be due to an increase with depth of the amount of pyroxene relative to olivine in these regions, although the presence of water, partial melt or a change in the deformation mechanism cannot be ruled out as yet. A similar observation is made for old oceans, but not for young ones where VSH> VSV appears likely down to 670 km depth and VPH> VPV down to 400 km depth. The change of sign in P-wave anisotropy seems to be qualitatively correlated with the presence of the Lehmann discontinuity, generally observed beneath continents and some oceans but not beneath ridges. Parameter η shows a similar age-related depth pattern as shear-wave anisotropy in the uppermost mantle and it undergoes the same change of sign as P-wave anisotropy at 220 km depth. The ratio between dln Vs and dln Vp suggests that a chemical component is needed to explain the anomalies in most places at depths greater than 220 km. More tests are needed to infer the robustness of the results for density, but they do not affect the results for anisotropy.

Differentiating flow, melt, or fossil seismic anisotropy beneath Ethiopia

NASA Astrophysics Data System (ADS)

Hammond, J. O. S.; Kendall, J.-M.; Wookey, J.; Stuart, G. W.; Keir, D.; Ayele, A.

2014-05-01

Ethiopia is a region where continental rifting gives way to oceanic spreading. Yet the role that pre-existing lithospheric structure, melt, mantle flow, or active upwellings may play in this process is debated. Measurements of seismic anisotropy are often used to attempt to understand the contribution that these mechanisms may play. In this study, we use new data in Afar, Ethiopia along with legacy data across Ethiopia, Djibouti, and Yemen to obtain estimates of mantle anisotropy using SKS-wave splitting. We show that two layers of anisotropy exist, and we directly invert for these. We show that fossil anisotropy with fast directions oriented northeast-southwest may be preserved in the lithosphere away from the rift. Beneath the Main Ethiopian Rift and parts of Afar, anisotropy due to shear segregated melt along sharp changes in lithospheric thickness dominates the shear-wave splitting signal in the mantle. Beneath Afar, away from regions with significant lithospheric topography, melt pockets associated with the crustal and uppermost mantle magma storage dominate the signal in localized regions. In general, little anisotropy is seen in the uppermost mantle beneath Afar suggesting melt retains no preferential alignment. These results show the important role melt plays in weakening the lithosphere and imply that as rifting evolves passive upwelling sustains extension. A dominant northeast-southwest anisotropic fast direction is observed in a deeper layer across all of Ethiopia. This suggests that a conduit like plume is lacking beneath Afar today, rather a broad flow from the southwest dominates flow in the upper mantle.

Permeability of serpentinites and implication for the oceanic mantle hydration along the outer rise faults

NASA Astrophysics Data System (ADS)

Hatakeyama, K.; Katayama, I.

2016-12-01

Recent geophysical surveys indicate that hydration (serpentinization) of oceanic mantle is related to outer-rise faulting prior to subduction (e.g., Fujie et al., 2013, Shilington et al., 2015). The serpentinization of oceanic mantle influences the generation of intermediate-depth earthquakes (e.g., Seno and Yamanaka, 1996) and the subduction water flux (e.g., Hacker, 2008). Since the chemical reactions that produce serpentinite are geologically rapid at low temperatures (Martin and Fyfe, 1970), the rate of water delivery to the reaction front likely controls the extent of serpentinization (Macdonald and Fyfe, 1985). Because the water through existing serpentinite is supplied to reaction front, permeability of serpentinite has important role of the extent of serpentinization along the outer-rise fault. In this study, we measured permeability of low-temperature serpentinites composed of lizardite and chrysotile, and calculated the extent of serpentinization along an outer-rise fault from Darcy's law. Our experimental results show that the permeability of serpnetinites decreases with increasing confining pressure, and reaches to 10-19 m2 to 10-21 m2 at confining pressure of 100 MPa. In extrapolating our experimental results to pressure of oceanic mantle, permeability of serpentinite can be as low as 10-22 m2 at the top of oceanic mantle (7 km depth beneath seafloor). If we assume that the time scale of water supply to the reaction front of 1.0 My, the lateral extent of serpentinization is approximately 9 km along the outer-rise fault in the uppermost oceanic mantle. Based on these estimate, we calculated the global water flux carried by serpentinized oceanic mantle to be 3.4×1012 kg/year, which is markedly higher than the water flux of hydrated oceanic crust (1.3×1012 kg/year). Since the subduction water flux is much greater than the output flux through magmatic degassing, the amount of present-day ocean might be decreasing, and this may result in the disappearance of the Earth's oceans in the future.

Possible emplacement of crustal rocks into the forearc mantle of the Cascadia Subduction Zone

USGS Publications Warehouse

Calvert, A.J.; Fisher, M.A.; Ramachandran, K.; Trehu, A.M.

2003-01-01

Seismic reflection profiles shot across the Cascadia forearc show that a 5-15 km thick band of reflections, previously interpreted as a lower crustal shear zone above the subducting Juan de Fuca plate, extends into the upper mantle of the North American plate, reaching depths of at least 50 km. In the extreme western corner of the mantle wedge, these reflectors occur in rocks with P wave velocities of 6750-7000 ms-1. Elsewhere, the forearc mantle, which is probably partially serpentinized, exhibits velocities of approximately 7500 ms-1. The rocks with velocities of 6750-7000 ms-1 are anomalous with respect to the surrounding mantle, and may represent either: (1) locally high mantle serpentinization, (2) oceanic crust trapped by backstepping of the subduction zone, or (3) rocks from the lower continental crust that have been transported into the uppermost mantle by subduction erosion. The association of subparallel seismic reflectors with these anomalously low velocities favours the tectonic emplacement of crustal rocks. Copyright 2003 by the American Geophysical Union.

Seismic velocity structure of the incoming Pacific Plate subducting into the central part of the Japan Trench revealed by traveltime tomography using OBS data

NASA Astrophysics Data System (ADS)

Obana, K.; Fujie, G.; Kodaira, S.; Takahashi, T.; Yamamoto, Y.; Miura, S.; Shinohara, M.

2016-12-01

Subduction of oceanic plates plays an important role in the water transportation from the earth surface into the deep mantle. Recent active seismic survey studies succeed to image that the seismic velocities within the oceanic crust and the uppermost mantle in the outer rise region decreases toward the trench axis. These velocity changes are considered as an indication of the hydration and alteration of the incoming oceanic plates prior to the subduction. However, the area with sufficient resolution of the active seismic studies is often limited at depths corresponding to the oceanic crust and several km beneath the oceanic Moho. In this study, we have examined the seismic velocity structure of the incoming/subducting Pacific Plate beneath the trench axis and outer trench-slope of the central part of the Japan Trench. The seismicity in the Pacific Plate, including several M7-class intra-plate earthquakes, has been active since the 2011 Tohoku-Oki earthquake in the study area. These activities were observed by the ocean bottom seismographs (OBS) deployed repeatedly. The data obtained from these OBS observations allow us to resolve the seismic velocity structures at greater depths compared to the active seismic surveys. We conducted 3-D traveltime tomography by using double-difference tomography method (Zhang and Thurber, 2003). The results show that the seismic velocities within the oceanic mantle decreased toward the trench axis. The velocity reduction begins at about 80 km seaward of the trench axis and extended to a depth of at least 30 km beneath the trench axis area. If the observed P-wave velocity reduction from 8.4 km/s to 7.7 km/s at a depth of 15 km below the oceanic Moho is caused by the serpentinization of the oceanic mantle (Carlson and Miller, 2003), roughly 2.5 weight per cent of water is expected in the low velocity anomalies in the oceanic mantle.

Consequences of an unstable chemical stratification on mantle dynamics

NASA Astrophysics Data System (ADS)

Plesa, Ana-Catalina; Tosi, Nicola; Breuer, Doris

2013-04-01

Early in the history of terrestrial planets, the fractional crystallization of primordial magma oceans may have led to the formation of large scale chemical heterogeneities. These may have been preserved over the entire planetary evolution as suggested for Mars by the isotopic analysis of the so-called SNC meteorites. The fractional crystallization of a magma ocean leads to a chemical stratification characterized by a progressive enrichment in heavy elements from the core-mantle boundary to the surface. This results in an unstable configuration that causes the overturn of the mantle and the subsequent formation of a stable chemical layering. Assuming scaling parameters appropriate for Mars, we first performed simulations of 2D thermo-chemical convection in Cartesian geometry with the numerical code YACC [1]. We investigated systems heated either solely from below or from within by varying systematically the buoyancy ratio B, which measures the relative importance of chemical to thermal buoyancy, and the mantle rheology, by considering systems with constant, strongly temperature-dependent and plastic viscosity. We ran a large set of simulations spanning a wide parameter space in order to understand the basic physics governing the magma ocean cumulate overturn and its consequence on mantle dynamics. Moreover, we derived scaling laws that relate the time over which chemical heterogeneities can be preserved (mixing time) and the critical yield stress (maximal yield stress that allows the lithosphere to undergo brittle failure) to the buoyancy ratio. We have found that the mixing time increases exponentially with B, while the critical yield stress shows a linear dependence. We investigated then Mars' early thermo-chemical evolution using the code GAIA in a 2D cylindrical geometry [2] and assuming a detailed magma ocean crystallization sequence as obtained from geochemical modeling [3]. We used an initial composition profile adapted from [3], accounted for an exothermic phase transition between lower and upper mantle and assumed all radiogenic heat sources to be enriched during the freezing-phase of the magma ocean in the uppermost 50 km [4]. A stagnant lid forms rapidly because of the strong temperature dependence of the viscosity. This prevents the uppermost dense cumulates to sink, even when allowing for a plastic yielding mechanism. Below this dense stagnant lid, the mantle chemical gradient settles to a stable configuration. The convection pattern is dominated by small-scale structures, which are difficult to reconcile with the large-scale volcanic features observed over Mars' surface. Assuming that the stagnant lid will break, a stable density gradient is obtained, with the densest material and the entire amount of heat sources lying above the core-mantle-boundary. This leads to a strong overheating of the lowermost mantle, whose temperature increases to values that exceed the liquidus. Therefore a fractionated global and deep magma ocean is difficult to reconcile with observations. Different scenarios assuming, for instance, a hemispherical or shallow magma ocean will have to be considered. References [1] N. Tosi, D.A. Yuen and O. Dadek; EPSL (2010) (Yet Another Convection Code, https://code.google.com/p/yacc-convection/) [2] C. Huettig and K. Stemmer; PEPI (2008) [3] L.T. Elkins-Tanton, E.M. Parmentier and P.C. Hess; Meteoritic and Planetary Science (2003) [4] L.T. Elkins-Tanton, S.E. Zaranek, E.M. Parmentier and P.C. Hess; EPSL (2005)

Seismic structure of the crust and uppermost mantle of north America and adjacent oceanic basins: A synthesis

USGS Publications Warehouse

Chulick, G.S.; Mooney, W.D.

2002-01-01

We present a new set of contour maps of the seismic structure of North America and the surrounding ocean basins. These maps include the crustal thickness, whole-crustal average P-wave and S-wave velocity, and seismic velocity of the uppermost mantle, that is, Pn and Sn. We found the following: (1) The average thickness of the crust under North America is 36.7 km (standard deviation [s.d.] ??8.4 km), which is 2.5 km thinner than the world average of 39.2 km (s.d. ?? 8.5) for continental crust; (2) Histograms of whole-crustal P- and S-wave velocities for the North American crust are bimodal, with the lower peak occurring for crust without a high-velocity (6.9-7.3 km/sec) lower crustal layer; (3) Regions with anomalously high average crustal P-wave velocities correlate with Precambrian and Paleozoic orogens; low average crustal velocities are correlated with modern extensional regimes; (4) The average Pn velocity beneath North America is 8.03 km/sec (s.d. ?? 0.19 km/sec); (5) the well-known thin crust beneath the western United States extends into northwest Canada; (6) the average P-wave velocity of layer 3 of oceanic crust is 6.61 km/ sec (s.d. ?? 0.47 km/sec). However, the average crustal P-wave velocity under the eastern Pacific seafloor is higher than the western Atlantic seafloor due to the thicker sediment layer on the older Atlantic seafloor.

Salt partitioning between water and high-pressure ices. Implication for the dynamics and habitability of icy moons and water-rich planetary bodies

NASA Astrophysics Data System (ADS)

Journaux, Baptiste; Daniel, Isabelle; Petitgirard, Sylvain; Cardon, Hervé; Perrillat, Jean-Philippe; Caracas, Razvan; Mezouar, Mohamed

2017-04-01

Water-rich planetary bodies including large icy moons and ocean exoplanets may host a deep liquid water ocean underlying a high-pressure icy mantle. The latter is often considered as a limitation to the habitability of the uppermost ocean because it would limit the availability of nutrients resulting from the hydrothermal alteration of the silicate mantle located beneath the deep ice layer. To assess the effects of salts on the physical properties of high-pressure ices and therefore the possible chemical exchanges and habitability inside H2O-rich planetary bodies, we measured partitioning coefficients and densities in the H2O-RbI system up to 450 K and 4 GPa; RbI standing as an experimentally amenable analog of NaCl in the H2O-salt solutions. We measured the partitioning coefficient of RbI between the aqueous fluid and ices VI and VII, using in-situ Synchrotron X-ray Fluorescence (XRF). With in-situ X-ray diffraction, we measured the unit-cell parameters and the densities of the high-pressure ice phases in equilibrium with the aqueous fluid, at pressures and temperatures relevant to the interior of planetary bodies. We conclude that RbI is strongly incompatible towards ice VI with a partitioning coefficient Kd(VI-L) = 5.0 (± 2.1) ṡ10-3 and moderately incompatible towards ice VII, Kd(VII-L) = 0.12 (± 0.05). RbI significantly increases the unit-cell volume of ice VI and VII by ca. 1%. This implies that RbI-poor ice VI is buoyant compared to H2O ice VI while RbI-enriched ice VII is denser than H2O ice VII. These new experimental results might profoundly impact the internal dynamics of water-rich planetary bodies. For instance, an icy mantle at moderate conditions of pressure and temperature will consist of buoyant ice VI with low concentration of salt, and would likely induce an upwelling current of solutes towards the above liquid ocean. In contrast, a deep and/or thick icy mantle of ice VII will be enriched in salt and hence would form a stable chemical boundary layer on top of the silicate mantle. Such a contrasted dynamics in the aqueous-ice VI-ice VII system would greatly influence the migration of nutrients towards the uppermost liquid ocean, thus controlling the habitability of moderate to large H2O-rich planetary bodies in our solar system (e.g., Ganymede, Titan, Calisto) and beyond.

Origin of dipping structures in fast-spreading oceanic lower crust offshore Alaska imaged by multichannel seismic data

NASA Astrophysics Data System (ADS)

Bécel, Anne; Shillington, Donna J.; Nedimović, Mladen R.; Webb, Spahr C.; Kuehn, Harold

2015-08-01

Multi-channel seismic (MCS) reflection profiles across the Pacific Plate south of the Alaska Peninsula reveal the internal structure of mature oceanic crust (48-56 Ma) formed at fast to intermediate spreading rates during and after a major plate re-organization. Oceanic crust formed at fast spreading rates (half spreading rate ∼ 74 mm /yr) has smoother basement topography, thinner sediment cover with less faulting, and an igneous section that is at least 1 km thicker than crust formed at intermediate spreading rates (half spreading rate ∼ 28- 34 mm /yr). MCS data across fast-spreading oceanic crust formed during plate re-organization contain abundant bright reflections, mostly confined to the lower crust above a highly reflective Moho transition zone, which has a reflection coefficient (RC) of ∼0.1. The lower crustal events dip predominantly toward the paleo-ridge axis at ∼10-30°. Reflections are also imaged in the uppermost mantle, which primarily dip away from the ridge at ∼10-25°, the opposite direction to those observed in the lower crust. Dipping events in both the lower crust and upper mantle are absent on profiles acquired across the oceanic crust formed at intermediate spreading rates emplaced after plate re-organization, where a Moho reflection is weak or absent. Our preferred interpretation is that the imaged lower crustal dipping reflections within the fast spread crust arise from shear zones that form near the spreading center in the region characterized by interstitial melt. The abundance and reflection amplitude strength of these events (RC ∼ 0.15) can be explained by a combination of solidified melt that was segregated within the shear structures, mylonitization of the shear zones, and crystal alignment, all of which can result in anisotropy and constructive signal interference. Formation of shear zones with this geometry requires differential motion between the crust and upper mantle, where the upper mantle moves away from the ridge faster than the crust. Active asthenospheric upwelling is one possible explanation for these conditions. The other possible interpretation is that lower crustal reflections are caused by magmatic (mafic/ultramafic) layering associated with accretion from a central mid-crustal magma chamber. Considering that the lower crustal dipping events have only been imaged in regions that have experienced plate re-organizations associated with ridge jumps or rift propagation, we speculate that locally enhanced mantle flow associated with these settings may lead to differential motion between the crust and the uppermost mantle, and therefore to shearing in the ductile lower crust or, alternatively, that plate reorganization could produce magmatic pulses which may lead to mafic/ultramafic banding.

Mantle Serpentinization near the Central Mariana Trench Constrained by Ocean Bottom Surface Wave Observations

NASA Astrophysics Data System (ADS)

Cai, C.; Wiens, D. A.; Lizarralde, D.; Eimer, M. O.; Shen, W.

2017-12-01

We investigate the crustal and uppermost mantle seismic structure across the Mariana trench by jointly inverting Rayleigh wave phase and group velocities from ambient noise and longer period phase velocities from Helmholtz tomography of teleseismic waveforms. We use data from a temporary deployment in 2012-2013, consisting of 7 island-based stations and 20 broadband ocean bottom seismographs, as well as data from the USGS Northern Mariana Islands Seismograph Network. To avoid any potential bias from the starting model, we use a Bayesian Monte-Carlo algorithm to invert for the azimuthally-averaged SV-wave velocity at each node. This method also allows us to apply prior constraints on crustal thickness and other parameters in a systematic way, and to derive formal estimates of velocity uncertainty. The results show the development of a low velocity zone within the incoming plate beginning about 80 km seaward of the trench axis, consistent with the onset of bending faults from bathymetry and earthquake locations. The maximum depth of the velocity anomaly increases towards the trench, and extends to about 30 km below the seafloor. The low velocities persist after the plate is subducted, as a 20-30 km thick low velocity layer with a somewhat smaller velocity reduction is imaged along the top of the slab beneath the forearc. An extremely low velocity zone is observed beneath the serpentine seamounts in the outer forearc, consistent with 40% serpentinization in the forearc mantle wedge. Azimuthal anisotropy results show trench parallel fast axis within the incoming plate at uppermost mantle depth (2%-4% anisotropy). All these observations suggest the velocity reduction in the incoming plate prior to subduction results from both serpentinized normal faults and water-filled cracks. Water is expelled from the cracks early in subduction, causing a modest increase in the velocity of the subducting mantle, and moves upward and causes serpentinization of the outer forearc. Assuming the velocity anomaly remaining in the subducting plate mantle is caused by serpentinization, calculations suggest the top 20 km of the slab mantle retains 10-15% serpentinization beyond the outer forearc. The amount of water carried into the deep mantle by this layer ( 54 Tg/Myr/m) is two to three times greater than previous estimates for the entire slab.

CaSiO3 perovskite in diamond indicates the recycling of oceanic crust into the lower mantle

NASA Astrophysics Data System (ADS)

Nestola, F.; Korolev, N.; Kopylova, M.; Rotiroti, N.; Pearson, D. G.; Pamato, M. G.; Alvaro, M.; Peruzzo, L.; Gurney, J. J.; Moore, A. E.; Davidson, J.

2018-03-01

Laboratory experiments and seismology data have created a clear theoretical picture of the most abundant minerals that comprise the deeper parts of the Earth’s mantle. Discoveries of some of these minerals in ‘super-deep’ diamonds—formed between two hundred and about one thousand kilometres into the lower mantle—have confirmed part of this picture. A notable exception is the high-pressure perovskite-structured polymorph of calcium silicate (CaSiO3). This mineral—expected to be the fourth most abundant in the Earth—has not previously been found in nature. Being the dominant host for calcium and, owing to its accommodating crystal structure, the major sink for heat-producing elements (potassium, uranium and thorium) in the transition zone and lower mantle, it is critical to establish its presence. Here we report the discovery of the perovskite-structured polymorph of CaSiO3 in a diamond from South African Cullinan kimberlite. The mineral is intergrown with about six per cent calcium titanate (CaTiO3). The titanium-rich composition of this inclusion indicates a bulk composition consistent with derivation from basaltic oceanic crust subducted to pressures equivalent to those present at the depths of the uppermost lower mantle. The relatively ‘heavy’ carbon isotopic composition of the surrounding diamond, together with the pristine high-pressure CaSiO3 structure, provides evidence for the recycling of oceanic crust and surficial carbon to lower-mantle depths.

Next-generation seismic experiments: wide-angle, multi-azimuth, three-dimensional, full-waveform inversion

NASA Astrophysics Data System (ADS)

Bell, Rebecca; Morgan, Joanna; Warner, Michael

2016-04-01

There are many outstanding plate-tectonic scale questions that require us to know information about sub-surface physical properties, for example ascertaining the geometry and location of magma chambers and estimating the effective stress along plate boundary faults. These important scientific targets are often too deep, impractical and expensive for extensive academic drilling. Full-waveform inversion (FWI) is an advanced seismic imaging technique that has recently become feasible in three dimensions, and has been widely adopted by the oil and gas industry to image reservoir-scale targets at shallow-to-moderate depths. In this presentation we explore the potential for 3-D FWI, when combined with appropriate marine seismic acquisition, to recover high-resolution high-fidelity P-wave velocity models for sub-sedimentary targets within the crystalline crust and uppermost mantle. Using existing geological and geophysical models, we construct P-wave velocity models over three potential sub-sedimentary targets: the Soufrière Hills Volcano on Montserrat and its associated crustal magmatic system, the downgoing oceanic plate beneath the Nankai subduction margin, and the oceanic crust-uppermost mantle beneath the East Pacific Rise mid-ocean ridge. We use these models to generate realistic multi-azimuth 3-D synthetic seismic data, and attempt to invert these data to recover the original models. We explore the resolution and accuracy, sensitivity to noise and acquisition geometry, ability to invert elastic data using acoustic inversion codes, and the trade-off between low frequencies and starting velocity model accuracy. We will show that FWI applied to multi-azimuth, refracted, wide-angle, low-frequency data can resolve features in the deep crust and uppermost mantle on scales that are significantly better than can be achieved by any other geophysical technique, and that these results can be obtained using relatively small numbers (60-90) of ocean-bottom receivers combined with large numbers of air-gun shots. We demonstrate that multi-azimuth 3-D FWI is robust in the presence of noise, that acoustic FWI can invert elastic data successfully, and that the typical errors to be expected in starting models derived using travel times will not be problematic for FWI given appropriately designed acquisition. In this presentation we will also discuss a recent field-example of the use of FWI to image the Endeavour spreading centre in the northeastern Pacific. FWI is a rapidly maturing technology; its transfer from the petroleum sector to tackle a broader range of targets now appears entirely achievable.

Seismic anisotropy and mantle creep in young orogens

USGS Publications Warehouse

Meissner, R.; Mooney, W.D.; Artemieva, I.

2002-01-01

Seismic anisotropy provides evidence for the physical state and tectonic evolution of the lithosphere. We discuss the origin of anisotropy at various depths, and relate it to tectonic stress, geotherms and rheology. The anisotropy of the uppermost mantle is controlled by the orthorhombic mineral olivine, and may result from ductile deformation, dynamic recrystallization or annealing. Anisotropy beneath young orogens has been measured for the seismic phase Pn that propagates in the uppermost mantle. This anisotropy is interpreted as being caused by deformation during the most recent thermotectonic event, and thus provides information on the process of mountain building. Whereas tectonic stress and many structural features in the upper crust are usually orientated perpendicular to the structural axis of mountain belts, Pn anisotropy is aligned parallel to the structural axis. We interpret this to indicate mountain-parallel ductile (i.e. creeping) deformation in the uppermost mantle that is a consequence of mountain-perpendicular compressive stresses. The preferred orientation of the fast axes of some anisotropic minerals, such as olivine, is known to be in the creep direction, a consequence of the anisotropy of strength and viscosity of orientated minerals. In order to explain the anisotropy of the mantle beneath young orogens we extend the concept of crustal 'escape' (or 'extrusion') tectonics to the uppermost mantle. We present rheological model calculations to support this hypothesis. Mountain-perpendicular horizontal stress (determined in the upper crust) and mountain-parallel seismic anisotropy (in the uppermost mantle) require a zone of ductile decoupling in the middle or lower crust of young mountain belts. Examples for stress and mountain-parallel Pn anisotropy are given for Tibet, the Alpine chains, and young mountain ranges in the Americas. Finally, we suggest a simple model for initiating mountain parallel creep.

3D Density Structure of Oceanic Lithosphere Affected by A Plume: A Case Study from the Greater Jan Mayen-East Greenland Region (NE Atlantic)

NASA Astrophysics Data System (ADS)

Tan, P.; Sippel, J.; Breivik, A. J.; Scheck-Wenderoth, M.; Meeßen, C.

2017-12-01

Unraveling the density structure of the oceanic lithosphere north of Iceland is key for understanding the effects of the Iceland Plume on the mid-ocean ridges of the greater Jan Mayen-East Greenland Region. We use a data-integrative approach for 3D gravity modeling to develop new insights into the crust and upper mantle density structure of this region. First, we obtain the 3D density structure of the sediments and crust from interpretations of regional reflection and refraction seismic lines. Then, the temperature and density structure of the mantle between 50 and 250 km are derived from a published shear-wave velocity (Vs) tomography model. To assess the density configuration between the Moho and 50 km depth, we follow a combined forward and inverse 3D gravity modeling approach. The Vs tomography and derived density of the deeper mantle (>50 km depth) reveal that the low-density anomaly related to the Iceland plume gets weaker with increasing distance from the plume, i.e. from the strongly influenced Middle Kolbeinsey Ridge (MKR) to the Mohn's Ridge. The West Jan Mayen Fracture Zone is identified as a main mantle density contrast, indicative of differences in the thermal evolution of the ridge systems it separates. Beneath the MKR region, the low-density anomaly at depths of >50 km continues upwards into the uppermost mantle, where its lateral dimensions narrow considerably. This elongated density anomaly is consistent with a basement high and indicates a channelization of the Iceland plume effects. The NE-SW elongated mantle anomaly does not, however, coincide with the topographical NNE-SSW striking ridge axis. Thus, the modelled plume-affected oceanic lithosphere reveals discrepancies with the half-space cooling model. We discuss the 3D density model in terms of such spatial relations between deeper mantle anomalies and the shallow crustal structure.

Geophysical investigations of the area between the Mid-Atlantic Ridge and the Barents Sea: From water to the lithosphere-asthenosphere system

NASA Astrophysics Data System (ADS)

Grad, Marek; Mjelde, Rolf; Krysiński, Lech; Czuba, Wojciech; Libak, Audun; Guterch, Aleksander

2015-03-01

As a part of the large international panel "IPY Plate Tectonics and Polar Gateways" within the "4th International Polar Year" framework, extensive geophysical studies were performed in the area of southern Svalbard, between the Mid-Atlantic Ridge and the Barents Sea. Seismic investigations were performed along three refraction and wide-angle reflection seismic lines. Integrated with gravity data the seismic data were used to determine the structure of the oceanic crust, the transition between continent and ocean (COT), and the continental structures down to the lithosphere-asthenosphere system (LAB). We demonstrate how modeling of multiple water waves can be used to determine the sound velocity in oceanic water along a seismic refraction profile. Our 2D seismic and density models documents 4-9 km thick oceanic crust formed at the Knipovich Ridge, a distinct and narrow continent-ocean transition (COT), the Caledonian suture zone between Laurentia and Barentsia, and 30-35 km thick continental crust beneath the Barents Sea. The COT west of southern Spitsbergen expresses significant excess density (more than 0.1 g/cm3 in average), which is characteristic for mafic/ultramafic and high-grade metamorphic rocks. The results of the gravity modeling show relatively weak correlation of the density with seismic velocity in the upper mantle, which suggests that the horizontal differences between oceanic and continental mantle are dominated by mineralogical changes, although a thermal effect is also present. The seismic velocity change with depth suggests lherzolite composition of the uppermost oceanic mantle, and dunite composition beneath the continental crust.

Icelandic-type crust

USGS Publications Warehouse

Foulger, G.R.; Du, Z.; Julian, B.R.

2003-01-01

Numerous seismic studies, in particular using receiver functions and explosion seismology, have provided a detailed picture of the structure and thickness of the crust beneath the Iceland transverse ridge. We review the results and propose a structural model that is consistent with all the observations. The upper crust is typically 7 ?? 1 km thick, heterogeneous and has high velocity gradients. The lower crust is typically 15-30 ?? 5 km thick and begins where the velocity gradient decreases radically. This generally occurs at the V p ??? 6.5 km s-1 level. A low-velocity zone ??? 10 000 km2 in area and up to ??? 15 km thick occupies the lower crust beneath central Iceland, and may represent a submerged, trapped oceanic microplate. The crust-mantle boundary is a transition zone ???5 ?? 3 km thick throughout which V p increases progressively from ???7.2 to ???8.0 km s-1. It may be gradational or a zone of alternating high- and low-velocity layers. There is no seismic evidence for melt or exceptionally high temperatures in or near this zone. Isostasy indicates that the density contrast between the lower crust and the mantle is only ???90 kg m-3 compared with ???300 kg m-3 for normal oceanic crust, indicating compositional anomalies that are as yet not understood. The seismological crust is ???30 km thick beneath the Greenland-Iceland and Iceland-Faeroe ridges, and eastern Iceland, ???20 km beneath western Iceland, and ???40 km thick beneath central Iceland. This pattern is not what is predicted for an eastward-migrating plume. Low attenuation and normal V p/V s ratios in the lower crust beneath central and southwestern Iceland, and normal uppermost mantle velocities in general, suggest that the crust and uppermost mantle are subsolidus and cooler than at equivalent depths beneath the East Pacific Rise. Seismic data from Iceland have historically been interpreted both in terms of thin-hot and thick-cold crust models, both of which have been cited as supporting the plume hypothesis. This suggests that the plume model for Iceland is an a priori assumption rather than a hypothesis subject to testing. The long-extinct Ontong-Java Plateau, northwest India and Parana??, Brazil large igneous provinces, beneath which mantle plumes are not expected are all underlain by mantle low-velocity bodies similar to that beneath Iceland. A plume interpretation for the mantle anomaly beneath Iceland is thus not required.

Subduction factory 1. Theoretical mineralogy, densities, seismic wave speeds, and H2O contents

NASA Astrophysics Data System (ADS)

Hacker, Bradley R.; Abers, Geoffrey A.; Peacock, Simon M.

2003-01-01

We present a new compilation of physical properties of minerals relevant to subduction zones and new phase diagrams for mid-ocean ridge basalt, lherzolite, depleted lherzolite, harzburgite, and serpentinite. We use these data to calculate H2O content, density and seismic wave speeds of subduction zone rocks. These calculations provide a new basis for evaluating the subduction factory, including (1) the presence of hydrous phases and the distribution of H2O within a subduction zone; (2) the densification of the subducting slab and resultant effects on measured gravity and slab shape; and (3) the variations in seismic wave speeds resulting from thermal and metamorphic processes at depth. In considering specific examples, we find that for ocean basins worldwide the lower oceanic crust is partially hydrated (<1.3 wt % H2O), and the uppermost mantle ranges from unhydrated to ˜20% serpentinized (˜2.4 wt % H2O). Anhydrous eclogite cannot be distinguished from harzburgite on the basis of wave speeds, but its ˜6% greater density may render it detectable through gravity measurements. Subducted hydrous crust in cold slabs can persist to several gigapascals at seismic velocities that are several percent slower than the surrounding mantle. Seismic velocities and VP/VS ratios indicate that mantle wedges locally reach 60-80% hydration.

Geophysical and petrological modeling of the lower crust and uppermost mantle in the Variscan and Proterozoic surroundings of the Trans-European Suture Zone in Central Europe

NASA Astrophysics Data System (ADS)

Puziewicz, Jacek; Polkowski, Marcin; Grad, Marek

2017-04-01

High-quality seismic data on the lower crust and uppermost lithospheric mantle in the Central European part of the Trans European Suture Zone, together with thermal and gravimetric data, enables the quantitative modeling of the rocks occurring in those parts of the lithosphere, including their mineral compositions and the chemical composition of individual minerals. The P3 seismic profile is located at the SW margin of the East European Craton. The lower crust is dominated by gabbronoritic intrusions (plagioclase An45Ab55, clinopyroxene Di80Hed20, orthopyroxene En74Fs26), and the uppermost mantle is harzburgitic (olivine and orthopyroxene Mg# 0.91). The lower crust and upper mantle of the P1 seismic profile belong to the Trans European Suture Zone, albeit the upper crust is of Variscan affinity. The P1 lower crust has gabbronoritic composition which is layered from plagioclase-rich compositions on the top to the orthopyroxene-rich ones at the bottom (plagioclase An45Ab55, clinopyroxene Di80Hed20, orthopyroxene En85Fs15), and is lithologically different Proterozoic and Variscan surroundings. The 100 × 200 km eclogite slice (garnet Alm48Gr25Py27, clinopyroxene Di51Hed10Jd39), with a thickness of 5-10 km, occurs in the uppermost mantle sampled by the P1 profile. The Niedźwiedź Massif is located at the NE margin of the Bohemian Massif, which shows an exposed Variscan basement. The lower crust beneath the Niedźwiedź Massif consists of gabbroic rock of variable proportions of plagioclase (An45Ab55) and clinopyroxene (Di80Hed20), whereas the uppermost mantle is supposedly spinel harzburgite (olivine, ortho- and clinopyroxene Mg# 0.90). Our models show that the lowermost crust and uppermost mantle of the East European Craton do not continue to the SW into the Trans European Suture Zone in its Central European section in Poland.

Heterogeneous seismic anisotropy in the transition zone and uppermost lower mantle: evidence from South America, Izu-Bonin and Japan

NASA Astrophysics Data System (ADS)

Lynner, Colton; Long, Maureen D.

2015-06-01

Measurements of seismic anisotropy are commonly used to constrain deformation in the upper mantle. Observations of anisotropy at mid-mantle depths are, however, relatively sparse. In this study we probe the anisotropic structure of the mid-mantle (transition zone and uppermost lower mantle) beneath the Japan, Izu-Bonin, and South America subduction systems. We present source-side shear wave splitting measurements for direct teleseismic S phases from earthquakes deeper than 300 km that have been corrected for the effects of upper mantle anisotropy beneath the receiver. In each region, we observe consistent splitting with delay times as large as 1 s, indicating the presence of anisotropy at mid-mantle depths. Clear splitting of phases originating from depths as great as ˜600 km argues for a contribution from anisotropy in the uppermost lower mantle as well as the transition zone. Beneath Japan, fast splitting directions are perpendicular or oblique to the slab strike and do not appear to depend on the propagation direction of the waves. Beneath South America and Izu-Bonin, splitting directions vary from trench-parallel to trench-perpendicular and have an azimuthal dependence, indicating lateral heterogeneity. Our results provide evidence for the presence of laterally variable anisotropy and are indicative of variable deformation and dynamics at mid-mantle depths in the vicinity of subducting slabs.

Crustal accretion along the global mid-ocean ridge system based on basaltic glass and olivine-hosted melt inclusion compositions

NASA Astrophysics Data System (ADS)

Wanless, V. D.; Behn, M. D.

2015-12-01

The depth and distribution of crystallization at mid-ocean ridges controls the overall architecture of the oceanic crust, influences hydrothermal circulation, and determines geothermal gradients in the crust and uppermost mantle. Despite this, there is no overall consensus on how crystallization is distributed within the crust/upper mantle or how this varies with spreading rate. Here, we examine crustal accretion at mid-ocean ridges by combining crystallization pressures calculated from major element barometers on mid-ocean ridge basalt (MORB) glasses with vapor-saturation pressures from melt inclusions to produce a detailed map of crystallization depths and distributions along the global ridge system. We calculate pressures of crystallization from >11,500 MORB glasses from the global ridge system using two established major element barometers (1,2). Additionally, we use vapor-saturation pressures from >400 olivine-hosted melt inclusions from five ridges with variable spreading rates to constrain pressures and distributions of crystallization along the global ridge system. We show that (i) crystallization depths from MORB glasses increase and become less focused with decreasing spreading rate, (ii) maximum glass pressures are greater than the maximum melt inclusion pressure, which indicates that the melt inclusions do not record the deepest crystallization at mid-ocean ridges, and (iii) crystallization occurs in the lower crust/upper mantle at all ridges, indicating accretion is distributed throughout the crust at all spreading rates, including those with a steady-state magma lens. Finally, we suggest that the remarkably similar maximum vapor-saturation pressures (~ 3000 bars) in melt inclusion from all spreading rates reflects the CO2 content of the depleted upper mantle feeding the global mid-ocean ridge system. (1) Michael, P. & W. Cornell (1998), Journal of Geophysical Research, 103(B8), 18325-18356; (2) Herzberg, C. (2004), Journal of Petrology, 45(12), 2389.

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

NASA Astrophysics Data System (ADS)

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

2018-01-01

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

Understanding the geodynamic setting of São Miguel, Azores: A peculiar bit of mantle in the Central Atlantic

NASA Astrophysics Data System (ADS)

Wilson, M.; Houlie, N.; Khan, A.; Lithgow-Bertelloni, C. R.

2012-12-01

The Azores Plateau and Archipelago in the Central Atlantic Ocean has traditionally been considered as the surface expression of a deep mantle plume or hotspot that has interacted with a mid-ocean ridge. It is geodynamically associated with the triple junction between the North American, African and Eurasian plates. (Yang et al., 2006) used finite frequency seismic tomography to demonstrate the presence of a zone of low P-wave velocities (peak magnitude -1.5%) in the uppermost 200km of the mantle beneath the plateau. The tomographic model is consistent with SW deflection of a mantle plume by regional upper mantle shear flow driven by absolute plate motions. The volcanic island of Sao Miguel is located within the Terceira Rift, believed to represent the boundary between the African and Eurasian plates; magmatic activity has been characterised by abundant basaltic eruptions in the past 30,000 years. The basalts are distinctive within the spectrum of global ocean island basalts for their wide range in isotopic composition, particularly in 87Sr/86Sr. Their Sr-Nd-Pb isotopic compositions show systematic variations from west to east across the island which can be interpreted in terms of melting of a two-component mantle source. The low melting point (enriched) component in the source has been attributed to recycled ancient (~3 Ga) oceanic crust(Elliott et al., 2007). Using the thermo-barometry approach of (Lee et al., 2009) we demonstrate that the pressure and temperature of magma generation below Sao Miguel increase from west (2 GPa, 1425 °C) to east (3.8 GPa, 1575 °C), consistent with partial melting along a mantle geotherm with a potential temperature of ~ 1500 °C. This is consistent with the magnitude of the thermal anomaly beneath the Azores Plateau (ΔT ~ 150-200 °C) inferred on the basis of the seismic tomography study. The site of primary magma generation extends from the base of the local lithosphere (~ 50 km) to ~ 125 km depth. To understand the geodynamic setting of the Sao Miguel magmatism we combine GPS data and mantle convection models with our interpretation of the geochemistry of the basalts. We demonstrate strong south-westerly and downward flow in the asthenospheric mantle above the Transition Zone (410 km seismic discontinuity), consistent with a zone of upper mantle shearing below the base of the lithosphere. The maximum flow velocity is broadly consistent with the depth of magma generation. The advection of the mantle with respect to the oceanic plate "moves" an isotopically distinct mantle source component beneath the active volcanoes of Sao Miguel and carries its previous melting residues to the south-west. We discuss the nature of this mantle source and its contribution to the mantle velocity anomalies determined by seismic tomography. This study opens-up new perspectives for seismic tomography and potentially new connections between the fields of geophysics and geochemistry in oceanic domains.

Terrestrial aftermath of the Moon-forming impact.

PubMed

Sleep, Norman H; Zahnle, Kevin J; Lupu, Roxana E

2014-09-13

Much of the Earth's mantle was melted in the Moon-forming impact. Gases that were not partially soluble in the melt, such as water and CO2, formed a thick, deep atmosphere surrounding the post-impact Earth. This atmosphere was opaque to thermal radiation, allowing heat to escape to space only at the runaway greenhouse threshold of approximately 100 W m(-2). The duration of this runaway greenhouse stage was limited to approximately 10 Myr by the internal energy and tidal heating, ending with a partially crystalline uppermost mantle and a solid deep mantle. At this point, the crust was able to cool efficiently and solidified at the surface. After the condensation of the water ocean, approximately 100 bar of CO2 remained in the atmosphere, creating a solar-heated greenhouse, while the surface cooled to approximately 500 K. Almost all this CO2 had to be sequestered by subduction into the mantle by 3.8 Ga, when the geological record indicates the presence of life and hence a habitable environment. The deep CO2 sequestration into the mantle could be explained by a rapid subduction of the old oceanic crust, such that the top of the crust would remain cold and retain its CO2. Kinematically, these episodes would be required to have both fast subduction (and hence seafloor spreading) and old crust. Hadean oceanic crust that formed from hot mantle would have been thicker than modern crust, and therefore only old crust underlain by cool mantle lithosphere could subduct. Once subduction started, the basaltic crust would turn into dense eclogite, increasing the rate of subduction. The rapid subduction would stop when the young partially frozen crust from the rapidly spreading ridge entered the subduction zone. © 2014 The Author(s) Published by the Royal Society. All rights reserved.

Cascadia subducting plate fluids channelled to fore-arc mantle corner: ETS and silica deposition

USGS Publications Warehouse

Hyndman, Roy D.; McCrory, Patricia A.; Wech, Aaron; Kao, Han; Ague, Jay

2015-01-01

In this study we first summarize the constraints that on the Cascadia subduction thrust, there is a 70 km gap downdip between the megathrust seismogenic zone and the Episodic Tremor and Slip (ETS) that lies further landward; there is not a continuous transition from unstable to conditionally stable sliding. Seismic rupture occurs mainly offshore for this hot subduction zone. ETS lies onshore. We then suggest what does control the downdip position of ETS. We conclude that fluids from dehydration of the downgoing plate, focused to rise above the fore-arc mantle corner, are responsible for ETS. There is a remarkable correspondence between the position of ETS and this corner along the whole margin. Hydrated mineral assemblages in the subducting oceanic crust and uppermost mantle are dehydrated with downdip increasing temperature, and seismic tomography data indicate that these fluids have strongly serpentinized the overlying fore-arc mantle. Laboratory data show that such fore-arc mantle serpentinite has low permeability and likely blocks vertical expulsion and restricts flow updip within the underlying permeable oceanic crust and subduction shear zone. At the fore-arc mantle corner these fluids are released upward into the more permeable overlying fore-arc crust. An indication of this fluid flux comes from low Poisson's Ratios (and Vp/Vs) found above the corner that may be explained by a concentration of quartz which has exceptionally low Poisson's Ratio. The rising fluids should be silica saturated and precipitate quartz with decreasing temperature and pressure as they rise above the corner.

Sn-wave velocity structure of the uppermost mantle beneath the Australian continent

NASA Astrophysics Data System (ADS)

Wei, Zhi; Kennett, Brian L. N.; Sun, Weijia

2018-06-01

We have extracted a data set of more than 5000 Sn traveltimes for source-station pairs within continental Australia, with 3-D source relocation using Pn arrivals to improve data consistency. We conduct tomographic inversion for S-wave-speed structure down to 100 km using the Fast Marching Tomography (FMTOMO) method for the whole Australian continent. We obtain a 3-D model with potential resolution of 3.0° × 3.0°. The new S-wave-speed model provides strong constraints on structure in a zone that was previously poorly characterized. The S velocities in the uppermost mantle are rather fast, with patterns of variation generally corresponding to those for Pn. We find strong heterogeneities of Swave speed in the uppermost mantle across the entire continent of Australia with a close relation to crustal geological features. For instance, the cratons in the western Australia usually have high S velocities (>4.70 km s-1), while the volcanic regions on the eastern margin of Australia are characterized by low S velocities (<4.40 km s-1). Exploiting an equivalent Pn inversion, we also determine the Vp/Vs ratios across the whole continent. We find that most of the uppermost mantle has Vp/Vs between 1.65 and 1.85, but with patches in central Australia and in the east with much higher Vp/Vs ratios. Distinctive local anomalies on the eastern margin may indicate the positions of remnants of mantle plumes.

Fertile lithospheric mantle beneath the northwestern North China and its implication for the subduction of the Paleo-Asian Ocean

NASA Astrophysics Data System (ADS)

Dai, H. K.; Zheng, J.; Su, Y. P.; Xiong, Q.; Pan, S. K.

2017-12-01

The nature of the sub-continental lithospheric mantle (SCLM) beneath the western North China Craton (NCC) is poorly known, which hinders understanding the cratonic response to the southward subduction of the Paleo-Asian Ocean. Mineral chemical data of spinel lherzolite xenoliths from newly discovered Cenozoic Langshan basalts in the northwestern part of the craton have been integrated with data from other localities across the western NCC, to put constrains on the SCLM nature and to explore the reworking processes involved. Compositions of mineral cores (i.e., Mg# in olivine = 88 91) and P-T estimates ( 1.2 GPa, 950 oC) suggest the Langshan xenoliths/xenocrysts represent fragments of the uppermost SCLM and experienced <15% melt extraction. These characteristics are similar to those of mantle xenoliths from other locaties (Siziwangqi and Hannuoba) along the northern margin of the western NCC. Disequilibrium characteristics are observed in xenoliths/xenocrysts in this study, including pyroxene spongy coronae and compositionally zoned olivine. They are interpreted to be induced by partial melting and by ironic diffusion with silicate melts in the mantle respectively, shortly before the eruption of host basalt. Metasomatism is recorded in clinopyroxene cores by concomitant enrichments in light rare earth elements and high field strength elements and was likely related to the migration of silicate melts derived from a mantle modified by slab melts during the Paleozoic time. The SCLM along the northern margin of the western NCC is fertile in nature constrained by mantle xenoliths from several localities (Langshan in this study, Siziwangqi and Hannuoba in references). Considering 1) the coexistence of fertile lithospheric mantle (similar to the Phanerozoic SCLM of the eastern NCC) and the overlying ancient continental crust, and 2) the sharp decrease in lithospheric thickness from the inner part to the northern margin of the western NCC, the SCLM beneath the northwestern part should have been strongly rejuvenated or replaced by fertile and non-cratonic mantle. Combined with other geological evidence on the northwestern margin, the mantle replacement and metasomatism were likely triggered by southward subduction of the Paleo-Asian Ocean.

A resolved mantle anomaly as the cause of the Australian-Antarctic Discordance

NASA Astrophysics Data System (ADS)

Ritzwoller, M. H.; Shapiro, N. M.; Leahy, G. M.

2003-12-01

We present evidence for the existence of an Australian-Antarctic Mantle Anomaly (AAMA), which trends northwest-southeast (NW-SE) through the Australian-Antarctic Discordance (AAD) on the Southeast Indian Ridge (SEIR), is confined to the upper 120 km of the mantle beneath the AAD, and dips shallowly to the west so that it extends to a depth of about 150 km west of the AAD. Average temperatures within the AAMA are depressed about 100°C relative to surrounding lithosphere and suggest very rapid cooling of newly formed lithosphere at the AAD to an effective thermal age between 20 and 30 Ma. A convective down welling beneath the AAD is not consistent with the confinement of the AAMA in the uppermost mantle. In substantial agreement with the model of [1998], we argue that the AAMA is the suspended remnant of a slab that subducted at the Gondwanaland-Pacific convergent margin more than 100 Myr ago, foundered in the deeper mantle, and then ascended into the shallow mantle within the past 30 Myr, cutting any ties to deeper roots. The stability of the AAMA and its poor correlation with residual topography and gravity imply that it is approximately neutrally buoyant. The thermally induced density anomaly can be balanced by bulk iron depletion of less than 0.8%, consistent with the warmer conditions of formation for the Pacific than Indian lithosphere. We hypothesize that the low temperatures in the AAMA inhibit crustal formation and the AAD depth anomaly is formed at the intersection of the SEIR and the AAMA. The northward migration of the SEIR overriding the cold NW-SE trending AAMA therefore presents a simple kinematic explanation for both the V-shaped residual depth anomaly in the southeast Indian Ocean and the western migration of the AAD along the SEIR. Neither explanation requires the Pacific asthenospheric mantle to push westward and displace Indian asthenosphere. The AAMA may also act as a barrier to large-scale flows in the shallow asthenosphere and may therefore define a boundary for mantle convection and between the Indian and Pacific isotopic provinces. The westward dip of the AAMA would also favor along-axis flow from the Indian Ocean asthenosphere to the AAD that may contribute to the penetration of Indian Ocean mid-ocean ridge basalts into the AAD.

Analysis of converted S-waves and gravity anomaly along the Aegir Ridge: implications for crustal lithology

NASA Astrophysics Data System (ADS)

Rai, A. K.; Breivik, A. J.; Mjelde, R.; Hanan, B. B.; Ito, G.; Sayit, K.; Howell, S.; Vogt, P. R.; Pedersen, R.

2012-12-01

The Aegir Ridge is an extinct spreading ridge in North-East Atlantic ocean. A thinner than normal crust around the Aegir Ridge appears as a hole in the extensively magmatic surroundings. Its proximity to the Iceland hot-spot makes it particularly important for understanding the changing dynamics of hotspot-ridge interaction. An integrated seismic and dredging experiment was conduced during the summer of 2010 with the primary aim to understand the nature of magmatism along the ridge shortly before cessation of seafloor spreading through variations of sub-seafloor lithological properties. Here, we present results of analysis of converted shear-waves recorded on OBS-sesimic data, and ship-gravity data. The shear-wave study enables us to quantify the variation of Vp/Vs in the sediments, crust and the upper-most mantle. We also inverted the gravity data to determine the sub-seafloor density distribution. The P- to S- converted shear-waves were identified on 20 OBSs along a profile with a total length of 550 km parallel to the ridge-axis. The sedimentary section on top of the crystalline crust is well illuminated in the streamer data. The forward modelling of the OBS data reveals that the Vp/Vs ratio in sediments are as high as 4.8, decreasing rapidly to a value of 3.00, primarily due to compaction of sediments with depth. Identification of sufficient PnS and PSn phases enable us to model the crustal and upper-most mantle Vp/Vs. The upper crystalline crust requires a Vp/Vs value of 1.99 and 1.89 for the southern and the northern profiles respectively, to fit the observations. The lower crust and upper-most part of the mantle have a Vp/Vs of ~1.82 and 1.795 respectively. Slightly lower Vp and moderate increase in Vp/Vs in parts of the crust and upper mantle presumably indicate presence of faulting, fracturing in the crust and moderate degree of serpentinization of the upper mantle. A sub-seafloor density model is derived by non-linear inversion of the gravity anomaly. The distribution of sediments appear to control the short-wavelength features of the gravity data, whereas density variations are required in the upper mantle to optimally fit the overall gravity anomaly. Our results suggest certain degree of temperature and/or compositional heterogeneities towards the southern ends of Aegir Ridge, near the Iceland-Faroes Ridge.

S-wave attenuation structure beneath the northern Izu-Bonin arc

NASA Astrophysics Data System (ADS)

Takahashi, Tsutomu; Obana, Koichiro; Kodaira, Shuichi

2016-04-01

To understand temperature structure or magma distribution in the crust and uppermost mantle, it is essential to know their attenuation structure. This study estimated the 3-D S-wave attenuation structure in the crust and uppermost mantle at the northern Izu-Bonin arc, taking into account the apparent attenuation due to multiple forward scattering. In the uppermost mantle, two areas of high seismic attenuation (high Q -1) imaged beneath the volcanic front were mostly colocated with low-velocity anomalies. This coincidence suggests that these high- Q -1 areas in low-velocity zones are the most likely candidates for high-temperature regions beneath volcanoes. The distribution of random inhomogeneities indicated the presence of three anomalies beneath the volcanic front: Two were in high- Q -1 areas but the third was in a moderate- Q -1 area, indicating a low correlation between random inhomogeneities and Q -1. All three anomalies of random inhomogeneities were rich in short-wavelength spectra. The most probable interpretation of such spectra is the presence of volcanic rock, which would be related to accumulated magma intrusion during episodes of volcanic activity. Therefore, the different distributions of Q -1 and random inhomogeneities imply that the positions of hot regions in the uppermost mantle beneath this arc have changed temporally; therefore, they may provide important constraints on the evolutionary processes of arc crust and volcanoes.

The P wavespeed structure in the mantle to 800 km depth below the Philippines region: geodynamic implications

NASA Astrophysics Data System (ADS)

Wright, C.

2009-03-01

P waves from earthquakes south of Taiwan, recorded by the BATS seismic array and CWB seismic network, were used define the P wavespeed structure between depths of 100 and 800 km below the Philippines region. The presence of a low wavespeed zone in the upper mantle is inferred, although the details are unclear. Wavespeeds in the uppermost mantle are low, as expected for seismic energy propagating within an oceanic plate. The estimated depths of the 410- and 660-km discontinuities are 325 and 676 km respectively. The unusually shallow depth of the upper discontinuity below and to the east of Luzon is inferred by clearly resolving the travel-time branch produced by refraction through the transition zone. A possible explanation for the northern part of the region covered is that seismic energy reaches its maximum depth within or close to the cool, subducted oceanic South China Sea slab where subduction has been slow and relatively recent. Further south, however, the presence of a broken remnant of the South China Sea slab, formed during a period of shallower subduction, is suggested at depths below 300 km to explain the broad extent of the elevated 410-km discontinuity. The 660-km discontinuity is slightly deeper than usual, implying that low temperatures persist to lower mantle depths. The wavespeed gradients within the transition zone between depths of 450 and 610 km are higher than those predicted by both the pyrolite and piclogite models of the mantle, possibly due to the presence of water in the transition zone.

Petrochemical and petrophysical characterization of the lower crust and the Moho beneath the West African Craton, based on Xenoliths from Kimberlites

NASA Technical Reports Server (NTRS)

Haggerty, Stephen E.; Toft, Paul B.

1988-01-01

Additional evidence to the composition of the lower crust and uppermost mantle was presented in the form of xenolith data. Xenoliths from the 2.7-Ga West African Craton indicate that the Moho beneath this shield is a chemically and physically gradational boundary, with intercalations of garnet granulite and garnet eclogite. Inclusions in diamonds indicate a depleted upper mantle source, and zenolith barometry and thermometry data suggest a high mantle geotherm with a kink near the Moho. Metallic iron in the xenoliths indicates that the uppermost mantle has a significant magnetization, and that the depth to the Curie isotherm, which is usually considered to be at or above the Moho, may be deeper than the Moho.

Physical and chemical consequences of crustal melting in fossil mature intra-oceanic arcs

NASA Astrophysics Data System (ADS)

Berger, J.; Burg, J.-P.

2012-04-01

Seismic velocity models of active intra-oceanic arcs show roots with densities and P-wave velocities intermediate to classical lower oceanic crust (density; ~3.0, Vp: ~7.0 km/s) and uppermost harzburgitic mantle (density: 3.2-3.3, Vp: 7.9-8.0 km/s). Most studies on active and fossil exhumed island arcs interpret the petrological nature of this root as ultramafic cumulates crystallized from primitive melts and/or as pyroxenites formed via basalt-peridotite reactions. Igneous cumulates and pyroxenites have densities close to or above that of uppermost mantle rocks; they can consequently undergo gravity-driven delamination, a process thought to drive the bulk composition of the arc toward an andesitic, continental crust-like composition. Dehydration and melting reactions are reported from exposed arc roots (Jijal complex in Kohistan; Amalaoulaou arc in Mali; Fiordland arc in New-Zealand). Intense influx of mantle-derived basaltic magmas at high pressure in a thickening island arc can enable lower crustal rocks to locally cross the dehydration-melting solidus of hydrous subalkaline basalts. Thermodynamic modeling using Perple_X, geochemical analysis and compilation of experimental and field data have been combined to constrain processes, conditions and consequences of intra-arc melting. The position of the solidus in a P-T grid is strongly dependent of the bulk water content: at 1 GPa, it is as low as 750 °C for water saturated hornblende-gabbros (>1 wt% H2O) and 830°C for gabbros with 0.1 wt% H2O. Incipient melting (F <10 %) near the solidus produces trondhjemitic melt and garnet granulites residue. The latter has composition very close to that of igneous precursors but is characterized by contrasted physical properties (density: 3.2-3.3, Vp: 6.9-7.4 km/s). Higher partial melting degrees (F: 10-20 %) lead to the formation of anorthositic melts in equilibrium with garnet-clinopyroxene-rutile residues (density: up to 3.45, Vp: up to 7.7 km/s). These melts are rich in LILE (Rb, Ba, Sr) and LREE but strongly depleted in HREE and Y, while the residues are moderately enriched in Ti, Zr, Nb, HREE and Y but depleted in LREE relative to their igneous precursors. Compared to depleted mantle values, the residues also have low Rb/Sr but high Sm/Nd and Lu/Hf ratios. Partial melting in the lowermost oceanic arc crust thus produces the conditions to trigger gravity-driven delamination of the root and could lead to introduction of fertile arc garnet pyroxenites within the upper mantle. However, in Kohistan and at Amalaoulaou, the dense garnet-clinopyroxene residues are dispersed in the arc roots; they are intermingled with hornblendite and pyroxenite bodies. The small density contrast between garnet granulites and the harzburgitic mantle, and the low volumes of garnet-clinopyroxene residues preclude massive delamination of the partial melting residues. Further numerical modeling of physical modifications induced by dehydration-melting together with igneous mineral segregation in arc roots will help constraining fundamental parameters (mantle and arc crust rheology and density, composition, P-T conditions, volume and rate of incoming basaltic fluxes…) that control the stability of the lowermost arc crust.

Local thickening of the Cascadia forearc crust and the origin of seismic reflectors in the uppermost mantle

USGS Publications Warehouse

Calvert, A.J.; Ramachandran, K.; Kao, H.; Fisher, M.A.

2006-01-01

Seismic reflection profiles from three different surveys of the Cascadia forearc are interpreted using P wave velocities and relocated hypocentres, which were both derived from the first arrival travel time inversion of wide-angle seismic data and local earthquakes. The subduction decollement, which is characterized beneath the continental shelf by a reflection of 0.5 s duration, can be traced landward into a large duplex structure in the lower forearc crust near southern Vancouver Island. Beneath Vancouver Island, the roof thrust of the duplex is revealed by a 5–12 km thick zone, identified previously as the E reflectors, and the floor thrust is defined by a short duration reflection from a − 1. We suggest that these relatively low velocities indicate the presence of either crustal rocks from the oceanic plate that have been underplated to the continent or crustal rocks from the forearc that have been transported downward by subduction erosion. The absence of seismicity from within the E reflectors implies that they are significantly weaker than the overlying crust, and the reflectors may be a zone of active ductile shear. In contrast, seismicity in parts of the D reflectors can be interpreted to mean that ductile shearing no longer occurs in the landward part of the duplex. Merging of the D and E reflectors at 42–46 km depth creates reflectivity in the uppermost mantle with a vertical thickness of at least 15 km. We suggest that pervasive reflectivity in the upper mantle elsewhere beneath Puget Sound and the Strait of Georgia arises from similar shear zones.

CRUST 5.1: A global crustal model at 5° x 5°

USGS Publications Warehouse

Mooney, Walter D.; Laske, Gabi; Masters, T. Guy

1998-01-01

We present a new global model for the Earth's crust based on seismic refraction data published in the period 1948–1995 and a detailed compilation of ice and sediment thickness. An extensive compilation of seismic refraction measurements has been used to determine the crustal structure on continents and their margins. Oceanic crust is modeled with both a standard model for normal oceanic crust, and variants for nonstandard regions, such as oceanic plateaus. Our model (CRUST 5.1) consists of 2592 5° × 5° tiles in which the crust and uppermost mantle are described by eight layers: (1) ice, (2) water, (3) soft sediments, (4) hard sediments, (5) crystalline upper, (6) middle, (7) lower crust, and (8) uppermost mantle. Topography and bathymetry are adopted from a standard database (ETOPO-5). Compressional wave velocity in each layer is based on field measurements, and shear wave velocity and density are estimated using recently published empirical Vp- Vs and Vp-density relationships. The crustal model differs from previous models in that (1) the thickness and seismic/density structure of sedimentary basins is accounted for more completely, (2) the velocity structure of unmeasured regions is estimated using statistical averages that are based on a significantly larger database of crustal structure, (3) the compressional wave, shear wave, and density structure have been explicitly specified using newly available constraints from field and laboratory studies. Thus this global crustal model is based on substantially more data than previous models and differs from them in many important respects. A new map of the thickness of the Earth's crust is presented, and we illustrate the application of this model by using it to provide the crustal correction for surface wave phase velocity maps. Love waves at 40 s are dominantly sensitive to crustal structure, and there is a very close correspondence between observed phase velocities at this period and those predicted by CRUST 5.1. We find that the application of crustal corrections to long-period (167 s) Rayleigh waves significantly increases the variance in the phase velocity maps and strengthens the upper mantle velocity anomalies beneath stable continental regions. A simple calculation of crustal isostacy indicates significant lateral variations in upper mantle density. The model CRUST 5.1 provides a complete description of the physical properties of the Earth's crust at a scale of 5° × 5° and can be used for a wide range of seismological and nonseismological problems.

Formation and evolution of a metasomatized lithospheric root at the motionless Antarctic plate: the case of East Island, Crozet Archipelago (Indian Ocean)

NASA Astrophysics Data System (ADS)

Meyzen, Christine; Marzoli, Andrea; Bellieni, Giuliano; Levresse, Gilles

2016-04-01

Sitting atop the nearly stagnant Antarctic plate (ca. 6.46 mm/yr), the Crozet archipelago midway between Madagascar and Antarctica constitutes a region of unusually shallow (1543-1756 m below sea level) and thickened oceanic crust (10-16.5 km), high geoid height, and deep low-velocity zone, which may reflect the surface expression of a mantle plume. Here, we present new major and trace element data for Quaternary sub-aerial alkali basalts from East Island, the easterly and oldest island (ca. 9 Ma) of the Crozet archipelago. Crystallization at uppermost mantle depth and phenocryst accumulation have strongly affected their parental magma compositions. Their trace element patterns show a large negative K anomaly relative to Ta-La, moderate depletions in Rb and Ba with respect to Th-U, and heavy rare earth element (HREE) depletions relative to light REE. These characteristics allow limits to be placed upon the composition and mineralogy of their mantle source. The average trace element spectrum of East Island basalts can be matched by melting of about 2 % of a garnet-phlogopite-bearing peridotite source. The stability field of phlogopite restricts melting depth to lithospheric levels. The modelled source composition requires a multistage evolution, where the mantle has been depleted by melt extraction before having been metasomatized by alkali-rich plume melts. The depleted mantle component may be sourced by residual mantle plume remnants stagnated at the melting locus due to a weak lateral flow velocity inside the melting regime, whose accumulation progressively edifies a depleted lithospheric root above the plume core. Low-degree alkali-rich melts are likely derived from the plume source. Such a mantle source evolution may be general to both terrestrial and extraterrestrial environments where the lateral component velocity of the mantle flow field is extremely slow.

Seismic Structure of India from Regional Waveform Matching

NASA Astrophysics Data System (ADS)

Gaur, V.; Maggi, A.; Priestley, K.; Rai, S.

2003-12-01

We use a neighborhood adaptive grid search procedure and reflectivity synthetics to model regional distance range (500-2000~km) seismograms recorded in India and to determine the variation in the crust and uppermost mantle structure across the subcontinent. The portions of the regional waveform which are most influenced by the crust and uppermost mantle structure are the 10-100~s period Pnl and fundamental mode surface waves. We use the adaptive grid search algorithm to match both portions of the seismogram simultaneously. This procedure results in a family of 1-D path average crust and upper mantle velocity and attenuation models whose propagation characteristics closely match those of the real Earth. Our data set currently consist of ˜20 seismograms whose propagation paths are primarily confined to the Ganges Basin in north India and the East Dharwar Craton of south India. The East Dharwar Craton has a simple and uniform structure consisting of a 36+/-2 km thick two layer crust, and an uppermost mantle with a sub-Moho velocity of 4.5~km/s. The structure of northern India is more complicated, with pronounced low velocities in the upper crustal layer due to the large sediment thicknesses in the Ganges basin.

Crust-mantle coupling mechanism in Cameroon, West Africa, revealed by 3D S-wave velocity and azimuthal anisotropy

NASA Astrophysics Data System (ADS)

Ojo, Adebayo Oluwaseun; Ni, Sidao; Chen, Haopeng; Xie, Jun

2018-01-01

To understand the depth variation of deformation beneath Cameroon, West Africa, we developed a new 3D model of S-wave isotropic velocity and azimuthal anisotropy from joint analysis of ambient seismic noise and earthquake surface wave dispersion. We found that the Cameroon Volcanic Line (CVL) is well delineated by slow phase velocities in contrast with the neighboring Congo Craton, in agreement with previous studies. Apart from the Congo Craton and the Oubanguides Belt, the uppermost mantle revealed a relatively slow velocity indicating a thinned or thermally altered lithosphere. The direction of fast axis in the upper crust is mostly NE-SW, but trending approximately N-S around Mt. Oku and the southern CVL. The observed crustal azimuthal anisotropy is attributed to alignment of cracks and crustal deformation related to magmatic activities. A widespread zone of weak-to-zero azimuthal anisotropy in the mid-lower crust shows evidence for vertical mantle flow or isotropic mid-lower crust. In the uppermost mantle, the fast axis direction changed from NE-SW to NW-SE around Mt. Oku and northern Cameroon. This suggests a layered mechanism of deformation and revealed that the mantle lithosphere has been deformed. NE-SW fast azimuths are observed beneath the Congo Craton and are consistent with the absolute motion of the African plate, suggesting a mantle origin for the observed azimuthal anisotropy. Our tomographically derived fast directions are consistent with the local SKS splitting results in some locations and depths, enabling us to constrain the origin of the observed splitting. The different feature of azimuthal anisotropy in the upper crust and the uppermost mantle implies decoupling between deformation of crust and mantle in Cameroon.

Xenolith constraints on seismic velocities in the upper mantle beneath southern Africa

NASA Astrophysics Data System (ADS)

James, D. E.; Boyd, F. R.; Schutt, D.; Bell, D. R.; Carlson, R. W.

2004-01-01

We impose geologic constraints on seismic three-dimensional (3-D) images of the upper mantle beneath southern Africa by calculating seismic velocities and rock densities from approximately 120 geothermobarometrically calibrated mantle xenoliths from the Archean Kaapvaal craton and adjacent Proterozoic mobile belts. Velocity and density estimates are based on the elastic and thermal moduli of constituent minerals under equilibrium P-T conditions at the mantle source. The largest sources of error in the velocity estimates derive from inaccurate thermo-barometry and, to a lesser extent, from uncertainties in the elastic constants of the constituent minerals. Results are consistent with tomographic evidence that cratonic mantle is higher in velocity by 0.5-1.5% and lower in density by about 1% relative to off-craton Proterozoic samples at comparable depths. Seismic velocity variations between cratonic and noncratonic xenoliths are controlled dominantly by differences in calculated temperatures, with compositional effects secondary. Different temperature profiles between cratonic and noncratonic regions have a relatively minor influence on density, where composition remains the dominant control. Low-T cratonic xenoliths exhibit a positive velocity-depth curve, rising from about 8.13 km/s at uppermost mantle depths to about 8.25 km/s at 180-km depth. S velocities decrease slightly over the same depth interval, from about 4.7 km/s in the uppermost mantle to 4.65 km/s at 180-km depth. P and S velocities for high-T lherzolites are highly scattered, ranging from highs close to those of the low-T xenoliths to lows of 8.05 km/s and 4.5 km/s at depths in excess of 200 km. These low velocities, while not asthenospheric, are inconsistent with seismic tomographic images that indicate high velocity root material extending to depths of at least 250 km. One plausible explanation is that high temperatures determined for the high-T xenoliths are a nonequilibrium consequence of relatively recent thermal perturbation and compositional modification associated with emplacement of kimberlitic fluids into the deep tectospheric root. Seismic velocities and densities for cratonic xenoliths differ significantly from those predicted for both primitive mantle peridotite and mantle eclogite. A model primitive mantle under cratonic P-T conditions exhibits velocities about 1% lower for P and about 1.5% lower for S, a consequence of a more fertile composition and different modal composition. Primitive mantle is also about 2% more dense at 150-km depth than low-T garnet lherzolite at cratonic P-T conditions. Similar calculations based on an oceanic geotherm are consistent with the isopycnic hypothesis of comparable density columns beneath oceanic and cratonic regions. Calculations for a hypothetical "cratonic" eclogite (50:50 garnet/omphacite) with an assumed cratonic geotherm produce extremely high VP and VS (8.68 km/s and 4.84 km/s, respectively, at 150 km depth) as well as high density (˜3.54 gm/cc). The very high velocity of eclogite should render it seismically conspicuous in the cratonic mantle if present as large volume blocks or slabs. We discuss how the seismic velocity data we have compiled in this paper from both xenoliths and generic petrologic models of the upper mantle differ from commonly used standard earth models IASPEI and PREM.

Effect of bend faulting on the hydration state of oceanic crust: Electromagnetic constraints from the Middle America Trench

NASA Astrophysics Data System (ADS)

Naif, S.; Key, K.; Constable, S.; Evans, R. L.

2017-12-01

In Northern Central America, the portion of the incoming Cocos oceanic plate formed at the East Pacific Rise has a seafloor spreading fabric that is oriented nearly parallel to the trench axis, whereby flexural bending at the outer rise reactivates a dense network of dormant abyssal hill faults. If bending-induced normal faults behave as fluid pathways they may promote extensive mantle hydration and significantly raise the flux of fluids into the subduction system. Multi-channel seismic reflection data imaged bend faults that extend several kilometers beneath the Moho offshore Nicaragua, coincident with seismic refraction data showing significant P-wave velocity reductions in both the crust and uppermost mantle. Ignoring the effect of fracture porosity, the observed mantle velocity reduction is equivalent to an upper bound of 15-20% serpentinization (or 2.0-2.5 wt% H2O). Yet the impact of bend faulting on porosity structure and crustal hydration are not well known. Here, we present results on the electrical resistivity structure of the incoming Cocos plate offshore Nicaragua, the first controlled-source electromagnetic (CSEM) experiment at a subduction zone. The CSEM data imaged several sub-vertical conductive channels extending beneath fault scarps to 5.5 km below seafloor, providing independent evidence for fluid infiltration into the oceanic crust via bending faults. We applied Archie's Law to estimate porosity from the resistivity observations: the dike and gabbro layers increase from 2.7% and 0.7% porosity at 100 km to 4.8% and 1.7% within 20 km of the trench, respectively. In contrast the resistivity, and hence porosity, remain relatively unchanged at sub-Moho depths. Therefore, either the faults do not provide an additional flux of free water to the mantle or, in light of the reduced seismic velocities, the volumetric expansion resulting from mantle serpentinization rapidly consumes any fault-generated porosity. Since our crustal porosity estimates seaward of the outer rise are in very good agreement with drilling observations, we conclude that bending faults effectively double the subducted free water budget of the intrusive oceanic crust.

Chemical heterogeneities in the interior of terrestrial bodies

NASA Astrophysics Data System (ADS)

Plesa, Ana-Catalina; Maurice, Maxime; Tosi, Nicola; Breuer, Doris

2016-04-01

Mantle chemical heterogeneities that can strongly influence the interior dynamics have been inferred for all terrestrial bodies of the Solar System and range from local to global scale. Seismic data for the Earth, differences in surface mineral compositions observed in data sets from space missions, and isotopic variations identified in laboratory analyses of meteorites or samples indicate chemically heterogeneous systems. One way to generate large scale geochemical heterogeneities is through the fractional crystallization of a liquid magma ocean. The large amount of energy available in the early stages of planetary evolution can cause melting of a significant part or perhaps even the entire mantle of a terrestrial body resulting in a liquid magma ocean. Assuming fractional crystallization, magma ocean solidification proceeds from the core-mantle boundary to the surface where dense cumulates tend to form due to iron enrichment in the evolving liquid. This process leads to a gravitationally unstable mantle, which is prone to overturn. Following cumulate overturn, a stable stratification may be reached that prevents efficient material transport. As a consequence, mantle reservoirs may be kept separate, possibly for the entire thermo-chemical evolution of a terrestrial body. Scenarios assuming fractional crystallization of a liquid magma ocean have been suggested to explain lavas with distinct composition on Mercury's surface [1], the generation of the Moon's mare basalts by sampling a reservoir consisting of overturned ilmenite-bearing cumulates [2], and the preservation of Mars' geochemical reservoirs as inferred by isotopic analysis of the SNC meteorites [3]. However, recent studies have shown that the style of the overturn as well as the subsequent density stratification are of extreme importance for the subsequent thermo-chemical evolution of a planetary body and may have a major impact on the later surface tectonics and volcanic history. The rapid formation of a stagnant lid that traps the uppermost dense cumulates close to the surface and prevents them from sinking into the mantle or the difficulty to initiate thermal convection because of the stable compositional gradient established after the overturn are difficult to reconcile with observations [4, 5]. More recent results show that the crystallization achieved upon solidification of a liquid magma ocean is considerably more complex than previously assumed. In fact, the onset of solid-state convection prior to complete crystallization of the mantle can efficiently reduce mantle chemical heterogeneities [5]. We thus conclude that mantle mixing may partly or even completely erase the effects of fractional crystallization well before complete solidification. Nevertheless, the subsequent differentiation caused by partial melting, may introduce additional heterogeneities between residual and primitive mantle that could explain compositional differences observed over the surface of terrestrial bodies [6]. References: [1] Charlier et al., 2013, EPSL; [2] Elkins-Tanton et al., 2011, EPSL; [3] Elkins-Tanton et al., 2005, JGR; [4] Tosi et al., 2013, JGR; [5] Plesa et al., 2014, EPSL; [5] Maurice et al, 2015, EGU; [6] Plesa & Breuer, 2014, PSS.

Seismic properties of the crust and uppermost mantle of North America

NASA Technical Reports Server (NTRS)

Braile, L. W.; Hinze, W. J.; Vonfrese, R. R. B.; Keller, G. R.

1983-01-01

Seismic refraction profiles for the North American continent were compiled. The crustal models compiled data on the upper mantle seismic velocity (P sub n), the crustal thickness (H sub c) and the average seismic velocity of the crystalline crust (V sub p). Compressional wave parameters were compared with shear wave data derived from surface wave dispersion models and indicate an average value for Poisson's ratio of 0.252 for the crust and of 0.273 for the uppermost mantle. Contour maps illustrate lateral variations in crustal thickness, upper mantle velocity and average seismic velocity of the crystalline crust. The distribution of seismic parameters are compared with a smoothed free air anomaly map of North America and indicate that a complidated mechanism of isostatic compensation exists for the North American continent. Several features on the seismic contour maps also correlate with regional magnetic anomalies.

High-resolution 3D seismic model of the crustal and uppermost mantle structure in Poland

NASA Astrophysics Data System (ADS)

Grad, Marek; Polkowski, Marcin; Ostaficzuk, Stanisław R.

2016-01-01

In the area of Poland a contact between the Precambrian and Phanerozoic Europe and the Carpathians has a complicated structure and a complex P-wave velocity of the sedimentary cover, crystalline crust, Moho depth and the uppermost mantle. The geometry of the uppermost several kilometers of sediments is relatively well recognized from over 100,000 boreholes. The vertical seismic profiling (VSP) from 1188 boreholes provided detailed velocity data for regional tectonic units and for stratigraphic successions from Permian to the Tertiary and Quaternary deposits. These data, however, do not provide information about the velocity and basement depth in the central part of the Trans-European suture zone (TESZ) and in the Carpathians. So, the data set is supplemented by 2D velocity models from 32 deep seismic sounding refraction profiles which also provide information about the crust and uppermost mantle. Together with the results of other methods: vertical seismic profiling, magnetotelluric, allow for the creation of a detailed, high-resolution 3D model for the entire Earth's crust and the uppermost mantle down to a depth of 60 km. The thinnest sedimentary cover in the Mazury-Belarus anteclise is only 0.3 to 1 km thick, which increases to 7 to 8 km along the East European Craton (EEC) margin, and 9 to 12 km in the TESZ. The Variscan domain is characterized by a 1-4 km thick sedimentary cover, while the Carpathians are characterized by very thick sedimentary layers, up to about 20 km. The crystalline crust is differentiated and has a layered structure. The crust beneath the West European Platform (WEP; Variscan domain) is characterized by P-wave velocities of 5.8-6.6 km/s. The upper and middle crusts beneath the EEC are characterized by velocities of 6.1-6.6 km/s, and are underlain by a high velocity lower crust with a velocity of about 7 km/s. A general decrease in velocity is observed from the older to the younger tectonic domains. The TESZ is associated with a steep dip in the Moho depth, from 30-35 km in the Paleozoic Platform to 42-52 km in the Precambrian craton. The new model confirms the Moho depth derived from previous compilations. In the TESZ the lower crust has a very high seismic velocity (> 7.0 km/s) which correlates to the high P-wave velocity (about 8.4 km/s) in the uppermost mantle beneath the Polish Basin. The Cratonic area is generally characterized by high P-wave velocities (> 8.2 km/s), while the Phanerozoic area is characterized by velocities of ~ 8.0 km/s. In the TESZ very high velocities of 8.3-8.4 km/s are observed, and the southwestern limitation of this area coincides with a high velocity lower crust, and could be continued to the NW toward the Elbe line. The influence of the structure for teleseismic tomography time residuals of seismic waves traveling through the 3D seismic model was analyzed. Lithological candidates for the crust and uppermost mantle of the EEC and WEP were suggested by comparison to laboratory data. The presented 3D seismic model may make more reliable studies on global dynamics, and geotectonic correlations, particularly for sedimentary basins in the Polish Lowlands, the napped flysch sediment series in the Carpathians, the basement shape, the southwestern edge of the EEC, a high-velocity lower crust and the high-velocity uppermost mantle in the TESZ. Finally, the new 3D velocity model of the crust shows a heterogeneous structure and offers a starting point for the numerical modeling of deeper structures by allowing for a correction of the crustal effects in studies of the mantle heterogeneities.

A Detailed 3D Seismic Velocity Structure of the Subducting Pacific Slab Beneath Hokkaido, Tohoku and Kanto, Japan, by Double-Difference Tomography

NASA Astrophysics Data System (ADS)

Tsuji, Y.; Nakajima, J.; Kita, S.; Okada, T.; Matsuzawa, T.; Hasegawa, A.

2007-12-01

Three-dimensional heterogeneous structure beneath northeastern (NE) Japan has been investigated by previous studies and an inclined seismic low-velocity zone is imaged in the mantle wedge sub-parallel to the down-dip direction of the subducting slab (Zhao et al., 1992, Nakajima et al., 2001). However, the heterogeneous structure within the slab has not been well studied even though it is very important to understand the whole process of water transportation from the slab to the surface. Here we show a detailed 3D seismic velocity structure within the subducted Pacific slab around Japan and propose a water-transportation path from the slab to the mantle wedge. In this study, we estimated 3D velocity structure within the Pacific slab by the double-difference tomography (Zhang and Thurber, 2003). We divided the study area, from Hokkaido to Kanto, into 6 areas due to the limitation of memory and computation time. In each area, arrival-time data of 7,500-17,000 events recorded at 70-170 stations were used in the analysis. The total number of absolute travel-time data was about 140,000-312,000 for P wave and 123,000-268,000 for S wave, and differential data were about 736,000-1,920,000 for P wave and 644,000-1,488,000 for S wave. Horizontal and vertical grid separations are 10-25 km and 6.5 km, respectively. RMS residuals of travel times for P wave decreased from 0.23s to 0.09s and for S wave from 0.35s to 0.13s. The obtained results are as follows: (1) a remarkable low-Vs zone exists in the uppermost part of the subducting slab, (2) it extends down to a depth of about 80 km, (3) the termination of this low-Vs zone almost corresponds to the "seismic belt" recently detected in the upper plane of the double seismic zone (Kita et al.,2006; Hasegawa et al., 2007), (4) at depths deeper than 80 km, a low-Vs and high-Vp/Vs zone is apparently distributed in the mantle wedge, immediately above the slab crust. We consider that these features reflect water-transportation processes from the slab to the mantle wedge. A low- Vs zone in the uppermost part of the subducting slab corresponds to the hydrous oceanic crust since its absolute velocity is about 4.0 km/s, comparable to that expected for the oceanic crust (Hacker et al., 2003). Dehydration reactions occur in the oceanic crust as temperature and pressure increase, and a relatively large amount of water is released at depths of about 80-100 km. The water generated by dehydration reactions could migrate upward and react peridotite at the base of the mantle wedge, forming a thin-serpentine layer there. Then, the layer is dragged by the subducting slab to deeper depths (e.g. Iwamori, 1998). Such water-transportation processes from the slab to the mantle wedge are partly constrained by a recent receiver function analysis (Kawakatsu and Watada, 2007). We further found an along-arc variation of the termination depth of the low-velocity oceanic crust, suggesting the along-arc variation in the amount of fluids released from the slab.

Crustal magmatism and lithospheric geothermal state of western North America and their implications for a magnetic mantle

NASA Astrophysics Data System (ADS)

Wang, Jian; Li, Chun-Feng

2015-01-01

The western North American lithosphere experienced extensive magmatism and large-scale crustal deformation due to the interactions between the Farallon and North American plates. To further understand such subduction-related dynamic processes, we characterize crustal structure, magmatism and lithospheric thermal state of western North America based on various data processing and interpretation of gravimetric, magnetic and surface heat flow data. A fractal exponent of 2.5 for the 3D magnetization model is used in the Curie-point depth inversion. Curie depths are mostly small to the north of the Yellowstone-Snake River Plain hotspot track, including the Steens Mountain and McDermitt caldera that are the incipient eruption locations of the Columbia River Basalts and Yellowstone hotspot track. To the south of the Yellowstone hotspot track, larger Curie depths are found in the Great Basin. The distinct Curie depths across the Yellowstone-Snake River Plain hotspot track can be attributed to subduction-related magmatism induced by edge flow around fractured slabs. Curie depths confirm that the Great Valley ophiolite is underlain by the Sierra Nevada batholith, which can extend further west to the California Coast Range. The Curie depths, thermal lithospheric thickness and surface heat flow together define the western edge of the North American craton near the Roberts Mountains Thrust (RMT). To the east of the RMT, large Curie depths, large thermal lithospheric thickness, and low thermal gradient are found. From the differences between Curie-point and Moho depth, we argue that the uppermost mantle in the oceanic region is serpentinized. The low temperature gradients beneath the eastern Great Basin, Montana and Wyoming permit magnetic uppermost mantle, either by serpentinization/metasomatism or in-situ magnetization, which can contribute to long-wavelength and low-amplitude magnetic anomalies and thereby large Curie-point depths.

Electrical conductivity imaging in the western Pacific subduction zone

NASA Astrophysics Data System (ADS)

Utada, Hisashi; Baba, Kiyoshi; Shimizu, Hisayoshi

2010-05-01

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

Three-dimensional thermal structure and seismogenesis in the Tohoku and Hokkaido subduction system

NASA Astrophysics Data System (ADS)

van Keken, P. E.; Kita, S.; Nakajima, J.; Bengtson, A. K.; Hacker, B. R.; Abers, G. A.

2010-12-01

The Northern Japan arc is characterized by fast subduction of old oceanic lithosphere. The high density instrumentation and high seismicity make this an ideal natural laboratory to study the interplay between subduction zone dynamics, dehydration, migration of fluids, and seismogenesis. In this study we use high resolution finite element models to predict the thermal structure of the subduction slab below Tohoku (Northern Honshu) and Hokkaido. These models allow us to predict the pressure, temperature and mineralogy of the subducted crust and mantle. We use these models to predict the (p,T) conditions of earthquakes that are relocated with a precision of around 1 km by double difference techniques. Below Northern Hokkaido and Tohoku we find that the earthquake activity is strong in crust and the uppermost mantle for temperatures < 450 C. Above this temperature earthquakes occur more sporadically and have significantly reduced integrated seismic moment. The strongest 3D variations in this arc occur below southern Hokkaido. This 200 km wide region is characterized by a change in trench geometry, anomalously low heatflow and an anomalous velocity structure in the mantle wedge. Tomographic imaging suggest that continental crust is subducted to significant depth, thereby insulating the subducting slab from the hot mantle wedge at least at intermediate depths. The thermal insulation is also suggested by the deepening of the earthquakes in the slab (Kita et al., EPSL, 2010). This region may be characterized by active crustal erosion which would lead to a further blanketing of the crust by a sedimentary layer. Further modifications in thermal structure are possible due to the 3D wedge flow that is generated by the along-arc variations in trench geometry. We quantitatively verify the relative importance of these processes using 2D and 3D dynamical models. Without the seismically imaged crustal structure the earthquake temperatures are significantly elevated compared to the Tohoku and (northern) Hokkaido sections. If we take the modified crustal structure into account we find a (p,T) pattern that is quite similar to that in the other sections, suggesting that the processes that lead to earthquakes in crust and uppermost mantle of the downgoing slab are similar across the northern Japan arc.

The role of harzburgite layers in the morphology of subducting plates and the behavior of oceanic crustal layers

NASA Astrophysics Data System (ADS)

Yoshida, Masaki

2014-05-01

Previous numerical studies of mantle convection focusing on subduction dynamics have indicated that the viscosity contrast between the subducting plate and the surrounding mantle have a primary effect on the behavior of subducting plates. The seismically observed plate stagnation at the base of the mantle transition zone (MTZ) under the Western Pacific and Eastern Eurasia is considered to mainly result from a viscosity increase at the ringwoodite to perovskite + magnesiowüstite (Rw→Pv+Mw) phase decomposition boundary, i.e., the boundary between the upper and lower mantle. The harzburgite layer, which is sandwiched between basaltic crust and depleted peridotite (lherzolite) layers, is a key component of highly viscous, cold oceanic plates. However, the possible sensitivity of the effective viscosity of harzburgite layers in the morphology of subducting plates that are flattened in the MTZ and/or penetrated in the lower mantle has not been examined systematically in previous three-dimensional (3D) numerical modeling studies that consider the viscosity increase at the boundary between the upper and lower mantle. In this study, in order to investigate the role of harzburgite layers in the morphology of subducting plates and the behavior of oceanic crustal layers, I performed a series of numerical simulations of mantle convection with semi-dynamic plate subduction in 3D regional spherical-shell geometry. The results show that a buckled crustal layer is observed under the "heel" of the stagnant slab that begins to penetrate into the lower mantle, regardless of the magnitude of the viscosity contrast between the harzburgite layer and the underlying mantle, when the factor of viscosity increase at the boundary of the upper and lower mantle is larger than 60-100. As the viscosity contrast between the harzburgite layer and the underlying mantle increases, the curvature of buckling is larger. When the viscosity increase at the boundary of the upper and lower mantle and the viscosity contrast between the harzburgite layer and the underlying mantle are larger, the volumes of crustal and harzburgite materials trapped in the mantle transition zone (MTZ) are also larger, although almost all of the materials penetrate into the lower mantle. These materials are trapped in the MTZ for over tens of millions of years. The bending of crustal layers numerically observed in the present study is consistent with seismological evidence that there is a piece of subducted oceanic crust in the uppermost lower mantle beneath the subducting slab under the Mariana trench [Niu et al., 2003, JGR]. The results of the present study suggest that when the viscosity increase at the boundary of the upper and lower mantle is larger than 60-100, a seismically observed stagnant slab is reproduced. This result is consistent with the previous independent geodynamic studies. For instance, a 2D geodynamic model with lateral viscosity variations suggested that it would need to be substantially greater than 30, say, around 100, to explain the positive geoid anomaly in the subduction zones where the subducting slab reaches the boundary between the upper and lower mantle such as that of the western Pacific [Tosi et al., 2009, GJI]. References: [1] Tajima, F. Yoshida, M. and Ohtani, E., Conjecture with water and rheological control for subducting slab in the mantle transition zone, Geoscience Frontiers, doi:10.1016/j.gsf.2013.12.005, 2014. [2] Yoshida, M. The role of harzburgite layers in the morphology of subducting plates and the behavior of oceanic crustal layers, Geophys. Res. Lett., 40(20), 5387-5392, doi:10.1002/2013GL057578, 2013. [3] Yoshida, M. and Tajima, F., On the possibility of a folded crustal layer stored in the hydrous mantle transition zone, Phys. Earth Planet. Inter., 219, 34-48, doi:10.1016/j.pepi.2013.03.004, 2013.

Molybdenum mobility and isotopic fractionation during subduction at the Mariana arc

NASA Astrophysics Data System (ADS)

Freymuth, Heye; Vils, Flurin; Willbold, Matthias; Taylor, Rex N.; Elliott, Tim

2015-12-01

The fate of crustal material recycled into the convecting mantle by plate tectonics is important for understanding the chemical and physical evolution of the planet. Marked isotopic variability of Mo at the Earth's surface offers the promise of providing distinctive signatures of such recycled material. However, characterisation of the behaviour of Mo during subduction is needed to assess the potential of Mo isotope ratios as tracers for global geochemical cycles. Here we present Mo isotope data for input and output components of the archetypical Mariana arc: Mariana arc lavas, sediments from ODP Sites 800, 801 and 802 near the Mariana trench and the altered mafic, oceanic crust (AOC), from ODP Site 801, together with samples of the deeper oceanic crust from ODP Site 1256. We also report new high precision Pb isotope data for the Mariana arc lavas and a dataset of Pb isotope ratios from sediments from ODP Sites 800, 801 and 802. The Mariana arc lavas are enriched in Mo compared to elements of similar incompatibility during upper mantle melting, and have distinct, isotopically heavy Mo (high 98Mo/95Mo) relative to the upper mantle, by up to 0.3 parts per thousand. In contrast, the various subducting sediment lithologies dominantly host isotopically light Mo. Coupled Pb and Mo enrichment in the Mariana arc lavas suggests a common source for these elements and we further use Pb isotopes to identify the origin of the isotopically heavy Mo. We infer that an aqueous fluid component with elevated [Mo], [Pb], high 98Mo/95Mo and unradiogenic Pb is derived from the subducting, mafic oceanic crust. Although the top few hundred metres of the subducting, mafic crust have a high 98Mo/95Mo, as a result of seawater alteration, tightly defined Pb isotope arrays of the Mariana arc lavas extrapolate to a fluid component akin to fresh Pacific mid-ocean ridge basalts. This argues against a flux dominantly derived from the highly altered, uppermost mafic crust or indeed from an Indian-like mantle wedge. Thus we infer that the Pb and Mo budgets of the fluid component are dominated by contributions from the deeper, less altered (cooler) portion of the subducting Pacific crust. The high 98Mo/95Mo of this flux is likely caused by isotopic fractionation during dehydration and fluid flow in the slab. As a result, the residual mafic crust becomes isotopically lighter than the upper mantle from which it was derived. Our results suggest that the continental crust produced by arc magmatism should have an isotopically heavy Mo composition compared to the mantle, whilst a contribution of deep recycled oceanic crust to the sources of some ocean island basalts might be evident from an isotopically light Mo signature.

Gravity and geoid anomalies of the Philippine Sea: Evidence on the depth of compensation for the negative residual water depth anomaly

NASA Technical Reports Server (NTRS)

Bowin, C.

1982-01-01

A negative free-air gravity anomaly which occurs in the central part of the Philippine Sea was examined to determine the distribution and nature of possible regional mass excesses or deficiencies. Geoid anomalies from GEOS-3 observation were positive. A negative residual geoid anomaly consistent with the area of negative free-air gravity anomalies were found. Theoretical gravity-topography and geoid-topography admittance functions indicated that high density mantle at about 60 km dept could account for the magnitudes of the gravity and residual geoid anomaly and the 1 km residual water depth anomaly in the Philippine Sea. The negative residual depth anomaly may be compensated for by excess density in the uppermost mantle, but the residual geoid and regional free-air gravity anomalies and a slow surface wave velocity structure might result from low-density warm upper mantle material lying beneath the zone of high-density uppermost mantle. From a horizontal disk approximation, the depth of the low-density warm mantle was estimated to be on the order of 200 km.

Flow, melt and fossil seismic anisotropy beneath Ethiopia

NASA Astrophysics Data System (ADS)

Hammond, James; Kendall, J.-Michael; Wookey, James; Stuart, Graham; Keir, Derek; Ayele, Atalay

2014-05-01

Ethiopia is a region where continental rifting gives way to oceanic spreading. Yet the role that pre-existing lithospheric structure, melt, mantle flow or active upwellings may play in this process is debated. Measurements of seismic anisotropy are often used to attempt to understand the contribution that these mechanisms may play. In this study we use new data in Afar, Ethiopia along with legacy data across Ethiopia, Djibouti and Yemen to obtain estimates of mantle anisotropy using SKS-wave splitting. We show that two layers of anisotropy exist, and use shear-wave splitting tomography to invert for these. We show that fossil anisotropy with fast directions oriented northeast-southwest may be preserved in the lithosphere away from the rift. Beneath the Main Ethiopian Rift and parts of Afar, anisotropy due aligned melt due to sharp changes in lithospheric thickness dominate the shear-wave splitting signal in the mantle. Beneath Afar, away from lithospheric topography, melt pockets associated with the crustal magma storage dominate the signal and little anisotropy is seen in the uppermost mantle suggesting melt retains no preferential alignment, possibly due to a lack of mantle lithosphere. These results show the important role melt plays in weakening the lithosphere and imply that as rifting evolves passive upwelling sustains extension. A dominant northeast-southwest anisotropic fast direction is observed in a deeper layer across all of Ethiopia. This suggests that a conduit like plume is absent beneath Afar today, rather a broad flow from the southwest dominates in the upper mantle.

Characterising the range of seismogenic behaviour on detachment faults - the case of 13o20'N, Mid Atlantic Ridge.

NASA Astrophysics Data System (ADS)

Craig, T. J.; Parnell-Turner, R.

2017-12-01

Extension at slow- and intermediate-spreading mid-ocean ridges is commonly accommodated through slip on long-lived detachment faults. These curved, convex-upward faults consist of a steeply-dipping section thought to be rooted in the lower crust or upper mantle which rotates to progressively shallower dip-angles at shallower depths, resulting in a domed, sub-horizontal oceanic core complex at the seabed. Although it is accepted that detachment faults can accumulate kilometre-scale offsets over millions of years, the mechanism of slip, and their capacity to sustain the shear stresses necessary to produce large earthquakes, remains debated. In this presentation we will show a comprehensive seismological study of an active oceanic detachment fault system on the Mid-Atlantic Ridge near 13o20'N, combining the results from a local ocean-bottom seismograph deployment with waveform inversion of a series of larger, teleseismically-observed earthquakes. The coincidence of these two datasets provides a more complete characterisation of rupture on the fault, from its initial beginnings within the uppermost mantle to its exposure at the surface. Our results demonstrate that although slip on the steeply-dipping portion of detachment fault is accommodated by failure in numerous microearthquakes, the shallower-dipping section of the fault within the upper few kilometres is relatively strong, and is capable of producing large-magnitude earthquakes. Slip on the shallow portion of active detachment faults at relatively low angles may therefore account for many more large-magnitude earthquakes at mid-ocean ridges than previously thought, and suggests that the lithospheric strength at slow-spreading mid-ocean ridges may be concentrated at shallow depths.

Seismic Structure of the Oceanic Plate Entering the Central Part of the Japan Trench Obtained from Ocean-Bottom Seismic Data

NASA Astrophysics Data System (ADS)

Ohira, A.; Kodaira, S.; Fujie, G.; No, T.; Nakamura, Y.; Miura, S.

2017-12-01

In trench-outer rise regions, the normal faults develop due to the bending of the incoming plate, which cause numerous normal-faulting earthquakes and systematic structural variations toward trenches. In addition to the effects on the bend-related normal fault, structural variations which are interpreted to be attributed to pseudofaults, a fracture zone, and petit-spot volcanic activities are observed in the oceanic plate entering the central part of the Japan Trench, off Miyagi. In May-June 2017, to understand detail structural variations and systematic structural changes of the oceanic plate toward the trench, we conducted an active-source seismic survey off Miyagi using R/V Kaimei, a new research vessel of JAMSTEC. Along a 100 km-long seismic profile which is approximately perpendicular to the trench axis, we deployed 40 ocean-bottom seismometers at intervals of 2 km and fired a large airgun array (total volume 10,600 cubic inches) with 100 m shooting intervals. Multi-channel seismic reflection data were also collected along the profile. On OBS records we observed refractions from the sedimentary layer and the oceanic crust (Pg), wide-angle reflections from the crust-mantle boundary (PmP), and refractions from the uppermost mantle (Pn). Pg is typically observed clearly at near offsets (approximately 20 km) but it highly attenuates at far offsets (> 20 km). A triplication of Pg-PmP-Pn with strong amplitudes is observed at ranges from 30 km to 60 km offsets. Pn is typically weak and its apparent velocity is approximately 8 km/sec. High attenuation of Pg and weak Pn may indicate the complex crustal structure related to petit-spot volcanic activities and/or a fracture zone, which are recognized in bathymetry data around the profile.

Electrically Anisotropic 35 Ma Pacific Lithosphere

NASA Astrophysics Data System (ADS)

Chesley, C. J.; Key, K.; Constable, S.; Behrens, J.; MacGregor, L.

2017-12-01

Geophysical studies of anisotropy in the oceanic lithosphere and asthenosphere can yield crucial insights into the processes of plate formation and evolution as the plate cools and thickens. While most previous studies have employed seismic methods to investigate anisotropy, here we examine the electrical conductivity anisotropy as constrained by controlled-source electromagnetic (CSEM) data collected during the Anisotropy and Physics of the Pacific Lithosphere Experiment (APPLE). Unlike passive magnetotelluric data, which are not particularly sensitive to the resistive part of the lithosphere or its anisotropy, CSEM data are highly sensitive to anisotropy in both the resistive crust and uppermost mantle. The APPLE data include a 30 km radius circular deep-tow of a Horizontal Electric Dipole (HED) transmitter around orthogonal pairs of HED receivers. The circular tow was optimized to measure azimuthal anisotropy, while radially oriented data at ranges from 14 to 70 km provided constraints on depth dependence of bulk conductivity. We inverted these data with a nonlinear anisotropic inversion that allows for laterally transverse isotropy, with the vertical plane of isotropy aligned orthogonal to the paleo-spreading direction. Our best model shows at least an order of magnitude resistivity difference between the paleo-spreading and paleo-ridge strike directions in both the crust and upper mantle. In the crust, conductivity is higher in the paleo-ridge and vertical directions. The opposite is true in the upper mantle, where conductivity is ten times higher in the paleo-spreading direction. Since the study area is centered on 35 Ma lithosphere, it is unlikely that melt plays a role in the observed anisotropy. Instead we propose that the crustal anisotropy is due to conductive clay minerals in normal faults promoted by hydration during paleo-extension close to the mid-ocean ridge. The upper mantle anisotropy potentially results from a crystal preferred orientation of olivine induced by shear deformation. These findings offer clues about the processes associated with oceanic spreading and may be of import to ophiolite studies.

Resolving Discrepancies Between Observed and Predicted Dynamic Topography on Earth

NASA Astrophysics Data System (ADS)

Richards, F. D.; Hoggard, M.; White, N. J.

2017-12-01

Compilations of well-resolved oceanic residual depth measurements suggest that present-day dynamic topography differs from that predicted by geodynamic simulations in two significant respects. At short wavelengths (λ ≤ 5,000 km), much larger amplitude variations are observed, whereas at long wavelengths (λ > 5,000 km), observed dynamic topography is substantially smaller. Explaining the cause of this discrepancy with a view to reconciling these different approaches is central to constraining the structure and dynamics of the deep Earth. Here, we first convert shear wave velocity to temperature using an experimentally-derived anelasticity model. This relationship is calibrated using a pressure and temperature-dependent plate model that satisfies age-depth subsidence, heat flow measurements, and seismological constraints on the depth to the lithosphere-asthenosphere boundary. In this way, we show that, at short-wavelengths, observed dynamic topography is consistent with ±150 ºC asthenospheric temperature anomalies. These inferred thermal buoyancy variations are independently verified by temperature measurements derived from geochemical analyses of mid-ocean ridge basalts. Viscosity profiles derived from the anelasticity model suggest that the asthenosphere has an average viscosity that is two orders of magnitude lower than that of the underlying upper mantle. The base of this low-viscosity layer coincides with a peak in azimuthal anisotropy observed in recent seismic experiments. This agreement implies that lateral asthenospheric flow is rapid with respect to the underlying upper mantle. We conclude that improved density and viscosity models of the uppermost mantle, which combine a more comprehensive physical description of the lithosphere-asthenosphere system with recent seismic tomographic models, can help to resolve spectral discrepancies between observed and predicted dynamic topography. Finally, we explore possible solutions to the long-wavelength discrepancy that exploit the velocity to density conversion described above combined with radial variation of mantle viscosity.

Depth-varying seismogenesis on an oceanic detachment fault at 13°20‧N on the Mid-Atlantic Ridge

NASA Astrophysics Data System (ADS)

Craig, Timothy J.; Parnell-Turner, Ross

2017-12-01

Extension at slow- and intermediate-spreading mid-ocean ridges is commonly accommodated through slip on long-lived faults called oceanic detachments. These curved, convex-upward faults consist of a steeply-dipping section thought to be rooted in the lower crust or upper mantle which rotates to progressively shallower dip-angles at shallower depths. The commonly-observed result is a domed, sub-horizontal oceanic core complex at the seabed. Although it is accepted that detachment faults can accumulate kilometre-scale offsets over millions of years, the mechanism of slip, and their capacity to sustain the shear stresses necessary to produce large earthquakes, remains subject to debate. Here we present a comprehensive seismological study of an active oceanic detachment fault system on the Mid-Atlantic Ridge near 13°20‧N, combining the results from a local ocean-bottom seismograph deployment with waveform inversion of a series of larger teleseismically-observed earthquakes. The unique coincidence of these two datasets provides a comprehensive definition of rupture on the fault, from the uppermost mantle to the seabed. Our results demonstrate that although slip on the deep, steeply-dipping portion of detachment faults is accommodated by failure in numerous microearthquakes, the shallow, gently-dipping section of the fault within the upper few kilometres is relatively strong, and is capable of producing large-magnitude earthquakes. This result brings into question the current paradigm that the shallow sections of oceanic detachment faults are dominated by low-friction mineralogies and therefore slip aseismically, but is consistent with observations from continental detachment faults. Slip on the shallow portion of active detachment faults at relatively low angles may therefore account for many more large-magnitude earthquakes at mid-ocean ridges than previously thought, and suggests that the lithospheric strength at slow-spreading mid-ocean ridges may be concentrated at shallow depths.

Seismic structure of the crust and uppermost mantle of South America and surrounding oceanic basins

USGS Publications Warehouse

Chulick, Gary S.; Detweiler, Shane; Mooney, Walter D.

2013-01-01

We present a new set of contour maps of the seismic structure of South America and the surrounding ocean basins. These maps include new data, helping to constrain crustal thickness, whole-crustal average P-wave and S-wave velocity, and the seismic velocity of the uppermost mantle (Pn and Sn). We find that: (1) The weighted average thickness of the crust under South America is 38.17 km (standard deviation, s.d. ±8.7 km), which is ∼1 km thinner than the global average of 39.2 km (s.d. ±8.5 km) for continental crust. (2) Histograms of whole-crustal P-wave velocities for the South American crust are bi-modal, with the lower peak occurring for crust that appears to be missing a high-velocity (6.9–7.3 km/s) lower crustal layer. (3) The average P-wave velocity of the crystalline crust (Pcc) is 6.47 km/s (s.d. ±0.25 km/s). This is essentially identical to the global average of 6.45 km/s. (4) The average Pn velocity beneath South America is 8.00 km/s (s.d. ±0.23 km/s), slightly lower than the global average of 8.07 km/s. (5) A region across northern Chile and northeast Argentina has anomalously low P- and S-wave velocities in the crust. Geographically, this corresponds to the shallowly-subducted portion of the Nazca plate (the Pampean flat slab first described by Isacks et al., 1968), which is also a region of crustal extension. (6) The thick crust of the Brazilian craton appears to extend into Venezuela and Colombia. (7) The crust in the Amazon basin and along the western edge of the Brazilian craton may be thinned by extension. (8) The average crustal P-wave velocity under the eastern Pacific seafloor is higher than under the western Atlantic seafloor, most likely due to the thicker sediment layer on the older Atlantic seafloor.

A record of spontaneous subduction initiation in the Izu-Bonin-Mariana arc

NASA Astrophysics Data System (ADS)

Arculus, Richard J.; Ishizuka, Osamu; Bogus, Kara A.; Gurnis, Michael; Hickey-Vargas, Rosemary; Aljahdali, Mohammed H.; Bandini-Maeder, Alexandre N.; Barth, Andrew P.; Brandl, Philipp A.; Drab, Laureen; Do Monte Guerra, Rodrigo; Hamada, Morihisa; Jiang, Fuqing; Kanayama, Kyoko; Kender, Sev; Kusano, Yuki; Li, He; Loudin, Lorne C.; Maffione, Marco; Marsaglia, Kathleen M.; McCarthy, Anders; Meffre, Sebastién; Morris, Antony; Neuhaus, Martin; Savov, Ivan P.; Sena, Clara; Tepley, Frank J., III; van der Land, Cees; Yogodzinski, Gene M.; Zhang, Zhaohui

2015-09-01

The initiation of tectonic plate subduction into the mantle is poorly understood. If subduction is induced by the push of a distant mid-ocean ridge or subducted slab pull, we expect compression and uplift of the overriding plate. In contrast, spontaneous subduction initiation, driven by subsidence of dense lithosphere along faults adjacent to buoyant lithosphere, would result in extension and magmatism. The rock record of subduction initiation is typically obscured by younger deposits, so evaluating these possibilities has proved elusive. Here we analyse the geochemical characteristics of igneous basement rocks and overlying sediments, sampled from the Amami Sankaku Basin in the northwest Philippine Sea. The uppermost basement rocks are areally widespread and supplied via dykes. They are similar in composition and age--as constrained by the biostratigraphy of the overlying sediments--to the 52-48-million-year-old basalts in the adjacent Izu-Bonin-Mariana fore-arc. The geochemical characteristics of the basement lavas indicate that a component of subducted lithosphere was involved in their genesis, and the lavas were derived from mantle source rocks that were more melt-depleted than those tapped at mid-ocean ridges. We propose that the basement lavas formed during the inception of Izu-Bonin-Mariana subduction in a mode consistent with the spontaneous initiation of subduction.

Geochemical constraints on the spatial distribution of recycled oceanic crust in the mantle source of late Cenozoic basalts, Vietnam

NASA Astrophysics Data System (ADS)

Hoang, Thi Hong Anh; Choi, Sung Hi; Yu, Yongjae; Pham, Trung Hieu; Nguyen, Kim Hoang; Ryu, Jong-Sik

2018-01-01

This study presents a comprehensive analysis of the major and trace element, mineral, and Sr, Nd, Pb and Mg isotopic compositions of late Cenozoic intraplate basaltic rocks from central and southern Vietnam. The Sr, Nd, and Pb isotopic compositions of these basalts define a tight linear array between Indian mid-ocean-ridge basalt (MORB)-like mantle and enriched mantle type 2 (EM2) components. These basaltic rocks contain low concentrations of CaO (6.4-9.7 wt%) and have high Fe/Mn ratios (> 60) and FeO/CaO-3MgO/SiO2 values (> 0.54), similar to partial melts derived from pyroxenite/eclogite sources. This similarity is also supported by the composition of olivine within these samples, which contains low concentration of Ca and high concentrations of Ni, and shows high Fe/Mn ratios. The basaltic rocks have elevated Dy/Yb ratios that fall within the range of melts derived from garnet lherzolite material, although their Yb contents are much higher than those of modeled melts derived from only garnet lherzolite material and instead plot near the modeled composition of eclogite-derived melts. The Vietnamese basaltic rocks have lighter δ26Mg values (- 0.38 ± 0.06‰) than is expected for the normal mantle (- 0.25 ± 0.07‰), and these values decrease with decreasing Hf/Hf* and Ti/Ti* ratios, indicating that these basalts were derived from a source containing carbonate material. On primitive mantle-normalized multi-element variation diagrams, the central Vietnamese basalts are characterized by positive Sr, Eu, and Ba anomalies. These basalts also plot within the pelagic sediment field in Pbsbnd Pb isotopic space. This suggests that the mantle source of the basalts contained both garnet peridotite and recycled oceanic crust. A systematic analysis of variations in geochemical composition in basalts from southern to central Vietnam indicates that the recycled oceanic crust (possibly the paleo-Pacific slab) source material contains varying proportions of gabbro, basalt, and sediment. The basalts from south-central Vietnam (12°N-14°N) may be dominated by the lowest portion of the residual slab that contains rutile-bearing plagioclase-rich gabbroic eclogite, whereas the uppermost portion of the recycled slab, including sediment and basaltic material with small amounts of gabbro, may be a major constituent of the source for the basalts within the central region of Vietnam (14°N-16°N). Finally, the southern region (10°N-12°N) contains basalts sourced mainly from recycled upper oceanic crust that is basalt-rich and contains little or no sediment.

The bright spot in the West Carpathian upper mantle: a trace of the Tertiary plate collision-and a caveat for a seismologist

NASA Astrophysics Data System (ADS)

Środa, Piotr

2010-07-01

The 2-D full waveform modelling of the mantle arrivals from the CELEBRATION 2000 profiles crossing the Carpathian orogen suggests two possible tectonic models for the collision of ALCAPA (Alpine-Carpathian-Pannonian) and the European Plate in the West Carpathians in southern Poland and Slovakia. Due to an oblique (NE-SW) convergence of plates, the character of the collision may change along the zone of contact of the plates: in the western part of the area an earlier collision might have caused substantial crustal shortening and formation of a crocodile-type structure, with the delaminated lower crust of ~100km length acting as a north-dipping reflecting discontinuity in the uppermost mantle. In the eastern part, a less advanced collision only involved the verticalization of the subducted slab remnant after a slab break-off. The lower crustal remnant of ~10km size in the uppermost mantle acts as a pseudo-diffractor generating observable mantle arrivals. Due to the similarity of synthetic data generated by both models, the question of the non-uniqueness of seismic data interpretation, that may lead to disparate tectonic inferences, is also discussed.

Cordilleran Longevity, Elevation and Heat Driven by Lithospheric Mantle Removal

NASA Astrophysics Data System (ADS)

Mackay-Hill, A.; Currie, C. A.; Audet, P.; Schaeffer, A. J.

2017-12-01

Cordilleran evolution is controlled by subduction zone back-arc processes that generate and maintain high topography due to elevated uppermost mantle temperatures. In the northern Canadian Cordillera (NCC), the persisting high mean elevation long after subduction has stopped (>50 Ma) requires a sustained source of heat either from small-scale mantle convection or lithospheric mantle removal; however direct structural constraints of these processes are sparse. We image the crust and uppermost mantle beneath the NCC using scattered teleseismic waves recorded on an array of broadband seismograph stations. We resolve two sharp and flat seismic discontinuities: a downward velocity increase at 35 km that we interpret as the Moho; and a deeper discontinuity with opposite velocity contrast at 50 km depth. Based on petrologic estimates, we interpret the deeper interface as the lithosphere-asthenosphere boundary (LAB), which implies an extremely thin ( 15 km) lithospheric mantle. We calculate the temperature at the Moho and the LAB in the range 800-900C and 1200-1300C, respectively. Below the LAB, we find west-dipping features far below the LAB beneath the eastern NCC that we associate with laminar downwelling of Cordilleran lithosphere. Whether these structures are fossilized or active, they suggest that lithospheric mantle removal near the Cordillera-Craton boundary may have provided the source of heat and elevation and therefore played a role in the longevity and stability of the Cordillera.

Seismic and thermal structure of the crust and uppermost mantle beneath Antarctica from inversion of multiple seismic datasets

NASA Astrophysics Data System (ADS)

Wiens, D.; Shen, W.; Anandakrishnan, S.; Aster, R. C.; Gerstoft, P.; Bromirski, P. D.; Dalziel, I.; Hansen, S. E.; Heeszel, D.; Huerta, A. D.; Nyblade, A.; Stephen, R. A.; Wilson, T. J.; Winberry, J. P.; Stern, T. A.

2017-12-01

Since the last decade of the 20th century, over 200 broadband seismic stations have been deployed across Antarctica (e.g., temporary networks such as TAMSEIS, AGAP/GAMSEIS, POLENET/ANET, TAMNNET and RIS/DRIS by U.S. geoscientists as well as stations deployed by Japan, Britain, China, Norway, and other countries). In this presentation, we discuss our recent efforts to build reference crustal and uppermost mantle shear velocity (Vs) and thermal models for continental Antarctica based on those seismic arrays. By combing the high resolution Rayleigh wave dispersion maps derived from both ambient noise and teleseismic earthquakes, together with P receiver function waveforms, we develop a 3-D Vs model for the crust and uppermost mantle beneath Central and West Antarctica to a depth of 200 km. Additionally, using this 3-D seismic model to constrain the crustal structure, we re-invert for the upper mantle thermal structure using the surface wave data within a thermodynamic framework and construct a 3-D thermal model for the Antarctic lithosphere. The final product, a high resolution thermal model together with associated uncertainty estimates from the Monte Carlo inversion, allows us to derive lithospheric thickness and surface heat flux maps for much of the continent. West Antarctica shows a much thinner lithosphere ( 50-90 km) than East Antarctica ( 130-230 km), with a sharp transition along the Transantarctic Mountains (TAM). A variety of geological features, including a slower/hotter but highly heterogeneous West Antarctica and a much faster/colder East Antarctic craton, are present in the 3-D seismic/thermal models. Notably, slow seismic velocities observed in the uppermost mantle beneath the southern TAM are interpreted as a signature of lithospheric foundering and replacement with hot asthenosphere. The high resolution image of these features from the 3-D models helps further investigation of the dynamic state of Antarctica's lithosphere and underlying asthenosphere and provides key constraints on the interaction between the solid Earth and the West Antarctic Ice Sheet.

Evolution of Meso-Cenozoic lithospheric thermal-rheological structure in the Jiyang sub-basin, Bohai Bay Basin, eastern North China Craton

NASA Astrophysics Data System (ADS)

Xu, Wei; Qiu, Nansheng; Wang, Ye; Chang, Jian

2018-01-01

The Meso-Cenozoic lithospheric thermal-rheological structure and lithospheric strength evolution of the Jiyang sub-basin were modeled using thermal history, crustal structure, and rheological parameter data. Results indicate that the thermal-rheological structure of the Jiyang sub-basin has exhibited obvious rheological stratification and changes over time. During the Early Mesozoic, the uppermost portion of the upper crust, middle crust, and the top part of the upper mantle had a thick brittle layer. During the early Early Cretaceous, the top of the middle crust's brittle layer thinned because of lithosphere thinning and temperature increase, and the uppermost portion of the upper mantle was almost occupied by a ductile layer. During the late Early Cretaceous, the brittle layer of the middle crust and the upper mantle changed to a ductile one. Then, the uppermost portion of the middle crust changed to a thin brittle layer in the late Cretaceous. During the early Paleogene, the thin brittle layer of the middle crust became even thinner and shallower under the condition of crustal extension. Currently, with the decrease in lithospheric temperature, the top of the upper crust, middle crust, and the uppermost portion of the upper mantle are of a brittle layer. The total lithospheric strength and the effective elastic thickness ( T e) in Meso-Cenozoic indicate that the Jiyang sub-basin experienced two weakened stages: during the late Early Cretaceous and the early Paleogene. The total lithospheric strength (approximately 4-5 × 1013 N m-1) and T e (approximately 50-60 km) during the Early Mesozoic was larger than that after the Late Jurassic (2-7 × 1012 N m-1 and 19-39 km, respectively). The results also reflect the subduction, and rollback of Pacific plate is the geodynamic mechanism of the destruction of the eastern North China Craton.

Mantle transition zone beneath the central Tien Shan: Lithospheric delamination and mantle plumes

NASA Astrophysics Data System (ADS)

Kosarev, Grigoriy; Oreshin, Sergey; Vinnik, Lev; Makeyeva, Larissa

2018-01-01

We investigate structure of the mantle transition zone (MTZ) under the central Tien Shan in central Asia by using recordings of seismograph stations in Kyrgyzstan, Kazakhstan and adjacent northern China. We apply P-wave receiver functions techniques and evaluate the differential time between the arrivals of seismic phases that are formed by P to SV mode conversion at the 410-km and 660-km seismic boundaries. The differential time is sensitive to the thickness of the MTZ and insensitive to volumetric velocity anomalies above the 410-km boundary. Under part of the southern central Tien Shan with the lowest S wave velocity in the uppermost mantle and the largest thickness of the crust, the thickness of the MTZ increases by 15-20 km relative to the ambient mantle and the reference model IASP91. The increased thickness is a likely effect of low (about - 150 K) temperature. This anomaly is indicative of delamination and sinking of the mantle lithosphere. The low temperature in the MTZ might also be a relic of subduction of the oceanic lithosphere in the Paleozoic, but this scenario requires strong coupling and coherence between structures in the MTZ and in the lithosphere during plate motions in the last 300 Myr. Our data reveal a reduction of thickness of the MTZ of 10-15 km under the Fergana basin, in the neighborhood of the region of small-scale basaltic volcanism at the time near the Cretaceous-Paleogene boundary. The reduced thickness of the MTZ is the effect of a depressed 410-km discontinuity, similar to that found in many hotspots. This depression suggests a positive temperature anomaly of about 100-150 K, consistent with the presence of a thermal mantle plume. A similar depression on the 410-km discontinuity is found underneath the Tarim basin.

S-Wave Velocity Models Under the Saudi Arabian Portable Broadband Deployment: Evidence for Lithospheric Erosion Beneath the Arabian Shield

NASA Astrophysics Data System (ADS)

Julià, J.; Ammon, C. J.; Herrmann, R. B.

2002-12-01

Models of crustal evolution strongly rely on our knowledge on the mineralogical composition of subsurface rocks, as well as pressure and temperature conditions. Direct sampling of subsurface rocks is often not possible, so that constraints have to be placed from indirect estimates of rock properties. Detailed seismic imaging of subsurface rocks has the potential for providing such constraints, and probe the extent at depth of surface geologic observations. In this study, we provide detailed S-wave velocity profiles for the crust and uppermost mantle beneath the Saudi Arabian Portable Broadband Deployment stations. Seismic velocities have been estimated from the joint inversion of receiver functions and fundamental mode group velocities. Receiver functions are sensitive to S-wave velocity contrasts and vertical travel times, and surface-wave dispersion is sensitive to vertical S-wave velocity averages, so that their combination bridge resolution gaps associated with each individual data set. Our resulting models correlate well with surface geology observations in the Arabian Shield and characterize its terranes at depth: the Asir terrane consists of a 10-km thick upper crust of 3.3~km/s overlying a lower crust with shear-wave velocities of 3.7-3.8 km/s; the Afif terrane is made of a 20-km thick upper crust with average velocity of 3.6 km/s and a lower crust with a shear-velocity of about 3.8~km/s; the Nabitah mobile belt has a gradational, 15-km thick upper crust up to 3.6 km/s overlying a gradational lower crust with velocities up to 4.0 km/s. The crust-mantle transition is sharper in terranes of continental affinity and more gradational beneath terranes of oceanic affinity. In the uppermost mantle, our models suggest a thin lid between up to 50-60 km depth overlying a low velocity zone beneath station TAIF, located close to a region of upwelling mantle material. Temperatures in the lid are estimated to be about 1000 C, which are in good agreement with independent xenolith data, and suggest that the lithosphere could be eroded to a thickness as little as 50~km under this station.

Faulting induced by precipitation of water at grain boundaries in hot subducting oceanic crust.

PubMed

Zhang, Junfeng; Green, Harry W; Bozhilov, Krassimir; Jin, Zhenmin

2004-04-08

Dehydration embrittlement has been proposed to explain both intermediate- and deep-focus earthquakes in subduction zones. Because such earthquakes primarily occur at shallow depths or within the core of the subducting plate, dehydration at relatively low temperatures has been emphasized. However, recent careful relocation of subduction-zone earthquakes shows that at depths of 100-250 km, earthquakes continue in the uppermost part of the slab (probably the former oceanic crust that has been converted to eclogite) where temperatures are higher. Here we show that at such pressures and temperatures, eclogite lacking hydrous phases but with significant hydroxyl incorporated as defects in pyroxene and garnet develops a faulting instability associated with precipitation of water at grain boundaries and the production of very small amounts of melt. This new faulting mechanism satisfactorily explains high-temperature earthquakes in subducting oceanic crust and could potentially be involved in much deeper earthquakes in connection with similar precipitation of water in the mantle transition zone (400-700 km depth). Of potential importance for all proposed high-pressure earthquake mechanisms is the very small amount of fluid required to trigger this instability.

Evidences of a Stalled-slab Beneath the Coast Ranges, California, From Seismicity and Converted Phases

NASA Astrophysics Data System (ADS)

Cao, A.; Liu, K. H.; Gao, S. S.

2001-12-01

In spite of numerous geophysical studies, the existence and geometry of a stalled slab beneath the Coast Ranges remains vague. In this study we use the distribution of mantle earthquakes and P-to-S converted phases from tilt interfaces to address the problem. Based on the CNSS catalog, in the period between 01/1960 and 04/2001, there were about 450 earthquakes occurred at depth larger than 35 km in the vicinity of the Coast Ranges. When plotted along east-west cross-sections, those earthquakes show a clear slab-like image, similar to the upper part of classic Benioff zones along subducting oceanic slabs. One of such cross-sections, which has a width of 20 km and a latitude of 39N, is located in the so-called 'slabless window' suggested by several previous geologic and geophysic studies, implying the existence of a stalled-slab along the cross-section. The mantle earthquakes can be explained as the result of stress concentration caused by heterogeneities in elastic properties associated with the cold slab, and of changes in mineralogical phases in the upper-most mantle in and around the slab. The existence of the slab is supported by clear azimuthal variations of the amplitude and arrival time of P-to-S converted phases from a tilt interface at about 70 km depth recorded by a broadband seismic station in the area. Our analysis shows that the converted phase is probably from a subducted oceanic lithosphere dipping to the east. The strike of the slab is approximately parallel to the Coast Ranges.

Seismicity and Structure of the Incoming Pacific Plate Subducting into the Japan Trench off Miyagi

NASA Astrophysics Data System (ADS)

Obana, K.; Fujie, G.; Kodaira, S.; Takahashi, T.; Yamamoto, Y.; Sato, T.; Yamashita, M.; Nakamura, Y.; Miura, S.

2015-12-01

Stresses within the oceanic plate in trench axis and outer-rise region have been characterized by shallow extension and deep compression due to the bending of the plate subducting into the trench. The stress state within the incoming/subducting oceanic plate is an important factor not only for the occurrence of shallow intraplate normal-faulting earthquakes in the trench-outer rise region but also the hydration of the oceanic plate through the shallow normal faults cutting the oceanic lithosphere. We investigate seismic velocity structure and stress state within the incoming/subducting Pacific Plate in the Japan Trench based on the OBS aftershock observations for the December 2012 intraplate doublet, which consists of a deep reverse faulting (Mw 7.2) and a shallow normal faulting (Mw 7.2) earthquake, in the Japan Trench off Miyagi. Hypocenter locations and seismic velocity structures were estimated from the arrival time data of about 3000 earthquakes by using double-difference tomography method (Zhang and Thurber, 2003). Also, focal mechanisms were estimated from first motion polarities by using the program HASH by Hardebeck and Shearer (2002). The results show that the earthquakes occurred mainly within the oceanic crust and the uppermost mantle. The deepest event was located at a depth of about 60 km. Focal mechanisms of the earthquakes shallower than a depth of 40 km indicate normal-faulting with T-axis normal to the trench. On the other hand, first motion polarities of the events at depths between 50 and 60 km can be explained a reverse faulting. The results suggest that the neutral plane of the stress between shallow extension and deep compression locates at 40 to 50 km deep. Seismic velocity structures indicate velocity decrease in the oceanic mantle toward the trench. Although the velocity decrease varies with locations, the results suggest the bending-related structure change could extend to at least about 15 km below the oceanic Moho in some locations.

Crust and uppermost mantle structure of the Kyushu-Palau Ridge, remnant arc on the Philippine Sea plate

NASA Astrophysics Data System (ADS)

Nishizawa, Azusa; Kaneda, Kentaro; Oikawa, Mitsuhiro

2016-02-01

We acquired 27 wide-angle seismic profiles to investigate variation in crustal structure along the Kyushu-Palau Ridge (KPR), a 2600-km-long remnant island arc in the center of the Philippine Sea plate; 26 lines were shot across the strike of the KPR at 13°-31°N, and one was shot along the northernmost KPR. The derived P-wave velocity (Vp) models show that the KPR has a crustal thickness of 8-23 km, which is thicker than the neighboring backarc basin oceanic crusts of the West Philippine Basin to the west and the Shikoku and Parece Vela Basins to the east. While the KPR crust consists mainly of lower crusts with a Vp of 6.8-7.2 km/s, the thicker crust contains a thick middle crust with Vp of 6.0-6.8 km/s. In general, the KPR crust is thicker in the north than in the south. The uppermost mantle velocities just below the KPR bathymetric highs are lower than 8.0 km/s and are commonly associated with a slightly high Vp of 7.2 km/s at the base of the crust. Large amplitude reflection signals are sometimes observed at far offsets on several lines suggesting the existence of several reflectors at depths of 23-40 km in the mantle beneath the KPR. The characteristics of these reflections are similar to these observed beneath the Izu-Ogasawara (Bonin) island arc, the tectonically conjugate arc of the KPR before backarc basin spreading. Very thin crust and high Pn velocities characterize the transition between the KPR and the eastern basins, which is probably a relic of the initial stage of the rifting. West of the KPR, the crust varies in structure from north to south as a result of the different tectonic settings in which it evolved.

Flow in the Deep Mantle from Seisimc Anisotropy: Progress and Prospects

NASA Astrophysics Data System (ADS)

Long, M. D.

2017-12-01

Observations of seismic anisotropy, or the directional dependence of seismic wavespeeds, provide one some of the most direct constraints on the pattern of flow in the Earth's mantle. In particular, as our understanding of crystallographic preferred orientation (CPO) of olivine aggregates under a range of deformation conditions has improved, our ability to exploit observations of upper mantle anisotropy has led to fundamental discoveries about the patterns of flow in the upper mantle and the drivers of that flow. It has been a challenge, however, to develop a similar framework for understanding flow in the deep mantle (transition zone, uppermost lower mantle, and lowermost mantle), even though there is convincing observational evidence for seismic anisotropy at these depths. Recent progress on the observational front has allowed for an increasingly detailed view of mid-mantle anisotropy (transition zone and uppermost lower mantle), particularly in subduction systems, which may eventually lead to a better understanding of mid-mantle deformation and the dynamics of slab interaction with the surrounding mid-mantle. New approaches to the observation and modeling of lowermost mantle anisotropy, in combination with constraints from mineral physics, are progressing towards interpretive frameworks that allow for the discrimination of different mantle flow geometries in different regions of D". In particular, observational strategies that involve the use of multiple types of body wave phases sampled over a range of propagation azimuths enable detailed forward modeling approaches that can discriminate between different mechanisms for D" anisotropy (e.g., CPO of post-perovskite, bridgmanite, or ferropericlase, or shape preferred orientation of partial melt) and identify plausible anisotropic orientations. We have recently begun to move towards a full waveform modeling approach in this work, which allows for a more accurate simulation for seismic wave propagation. Ongoing improvements in seismic observational strategies, experimental and computational mineral physics, and geodynamic modeling approaches are leading to new avenues for understanding flow in the deep mantle through the study of seismic anisotropy.

Spatiotemporal distribution of low-frequency earthquakes in Southwest Japan: Evidence for fluid migration and magmatic activity

NASA Astrophysics Data System (ADS)

Yu, Zhiteng; Zhao, Dapeng; Niu, Xiongwei; Li, Jiabiao

2018-01-01

Low-frequency earthquakes (LFEs) in the lower crust and uppermost mantle are widely observed in Southwest Japan, and they occur not only along the subducting Philippine Sea (PHS) slab interface but also beneath active arc volcanoes. The volcanic LFEs are still not well understood because of their limited quantities and less reliable hypocenter locations. In this work, seismic tomography is used to determine detailed three-dimensional (3-D) P- and S-wave velocity (Vp and Vs) models of the crust and upper mantle beneath Southwest Japan, and then the obtained 3-D Vp and Vs models are used to relocate the volcanic LFEs precisely. The results show that the volcanic LFEs can be classified into two types: pipe-like and swarm-like LFEs, and both of them are located in or around zones of low-velocity and high-Poisson's ratio anomalies in the crust and uppermost mantle beneath the active volcanoes. The pipe-like LFEs may be related to the fluid migration from the lower crust or the uppermost mantle, whereas the swarm-like LFEs may be related to local magmatic activities or small magma chambers. The number of LFEs sometimes increases sharply before or after a nearby large crustal earthquake which may cause cracks and fluid migration. The spatiotemporal distribution of the LFEs may indicate the track of migrating fluids. As compared with the tectonic LFEs along the PHS slab interface, the volcanic LFEs are more sensitive to fluid migration and local magmatic activities. High pore pressures play an important role in triggering both types of LFEs in Southwest Japan.

Along-axis crustal structure of the Porcupine Basin from seismic refraction data modelling

NASA Astrophysics Data System (ADS)

Prada, Manel; Watremez, Louise; Chen, Chen; O'Reilly, Brian; Minshull, Tim; Reston, Tim; Wagner, Gerlind; Gaws, Viola; Klaschen, Dirk; Shannon, Patrick

2016-04-01

The Porcupine Basin is a tongue-shaped offshore basin SW of Ireland that formed during the opening of the North Atlantic Ocean. Its history of development involved several rifting and subsidence phases during the Late Paleozoic and Cenozoic, with a particular major rift phase occurring in Late Jurassic-Early Cretaceous times. Previous work, focused on subsidence analysis, showed that stretching factors (β) in the northern part of the basin are < 1.5 and increase significantly southwards, where they were estimated to be > 6. However, recent studies based on seismic reflection and refraction profiles concluded that β in places along the basin axis were significantly higher, and suggested the presence of major crustal faulting and uppermost mantle serpentinization in the basin. Constraining β and the processes related to the formation of the basin will provide insights into aspects such as the tectonic response to lithospheric extension and the thermal evolution of the basin. Here we present the tomography results of five wide-angle seismic (WAS) profiles acquired across and along the basin axis. We used a travel time inversion method to model the WAS data and obtain P-wave velocity (Vp) models of the crust and uppermost mantle, together with the geometry of the main geological interfaces along each of these lines. Coincident seismic reflection profiles to each WAS line were also used to integrate the tectonic structure with the Vp model. These results improved constrains on the location of the base of the crust and allow to estimate maximum β (βmax) along each profile. The analysis shows that βmax values in the northern part of the basin are 5-6 times larger than estimates based on subsidence analysis. Towards the south, βmax increases up to 10, but then rapidly decreases to 3.3 southwards. These values are well within the range of crustal extension at which the crust becomes entirely brittle at magma-poor margins allowing the formation of major crustal faulting and serpentinization of the mantle. In agreement with this observation, Vp values of the mantle are lower than those expected for a non-altered mantle (i.e. ~8 km/s) supporting mantle serpentinization. The outcome of this study reveals the complexity of the crustal structure of the Porcupine Basin and demonstrates the importance and value of this type of analysis in understanding rift systems. This project is funded by the Irish Shelf Petroleum Studies Group (ISPSG) of the Irish Petroleum Infrastructure Programme Group 4.

Sea level fall during glaciation stabilized atmospheric CO2 by enhanced volcanic degassing

PubMed Central

Hasenclever, Jörg; Knorr, Gregor; Rüpke, Lars H.; Köhler, Peter; Morgan, Jason; Garofalo, Kristin; Barker, Stephen; Lohmann, Gerrit; Hall, Ian R.

2017-01-01

Paleo-climate records and geodynamic modelling indicate the existence of complex interactions between glacial sea level changes, volcanic degassing and atmospheric CO2, which may have modulated the climate system’s descent into the last ice age. Between ∼85 and 70 kyr ago, during an interval of decreasing axial tilt, the orbital component in global temperature records gradually declined, while atmospheric CO2, instead of continuing its long-term correlation with Antarctic temperature, remained relatively stable. Here, based on novel global geodynamic models and the joint interpretation of paleo-proxy data as well as biogeochemical simulations, we show that a sea level fall in this interval caused enhanced pressure-release melting in the uppermost mantle, which may have induced a surge in magma and CO2 fluxes from mid-ocean ridges and oceanic hotspot volcanoes. Our results reveal a hitherto unrecognized negative feedback between glaciation and atmospheric CO2 predominantly controlled by marine volcanism on multi-millennial timescales of ∼5,000–15,000 years. PMID:28681844

Evolution and hydration of the Juan de Fuca crust and uppermost mantle: a plate-scale seismic investigation from ridge to trench

NASA Astrophysics Data System (ADS)

Carbotte, S. M.; Canales, J.; Carton, H. D.; Nedimovic, M. R.; Han, S.; Marjanovic, M.; Gibson, J. C.; Janiszewski, H. A.; Horning, G.; Delescluse, M.; Watremez, L.; Farkas, A.; Biescas Gorriz, B.; Bornstein, G.; Childress, L. B.; Parker, B.

2012-12-01

The evolution of oceanic lithosphere involves incorporation of water into the physical and chemical structure of the crust and shallow mantle through fluid circulation, which initiates at the mid-ocean ridge and continues on the ridge flanks long after crustal formation. At subduction zones, water stored and transported with the descending plate is gradually released at depth, strongly influencing subduction zone processes. Cascadia is a young-lithosphere end member of the global subduction system where relatively little hydration of the downgoing Juan de Fuca (JdF) plate is expected due to its young age and presumed warm thermal state. However, numerous observations support the abundant presence of water within the subduction zone, suggesting that the JdF plate is significantly hydrated prior to subduction. Knowledge of the state of hydration of the JdF plate is limited, with few constraints on crustal and upper mantle structure. During the Cascadia Ridge-to-Trench experiment conducted in June-July 2012 over 4000 km of active source seismic data were acquired as part of a study of the evolution and state of hydration of the crust and shallow mantle of the JdF plate prior to subduction at the Cascadia margin. Coincident long-streamer (8 km) multi-channel seismic (MCS) and wide-angle ocean bottom seismometer (OBS) data were acquired in a two-ship program with the R/V Langseth (MGL1211), and R/V Oceanus (OC1206A). Our survey included two ridge-perpendicular transects across the full width of the JdF plate, a long trench-parallel line ~10 km seaward of the Cascadia deformation front, as well as three fan lines to study mantle anisotropy. The plate transects were chosen to provide reference sections of JdF plate evolution over the maximum range of JdF plate ages (8-9 Ma), offshore two contrasting regions of the Cascadia Subduction zone, and provide the first continuous ridge-to-trench images acquired at any oceanic plate. The trench-parallel line was designed to characterize variations in plate structure and hydration linked to JdF plate segmentation for over 450 km along the margin. Shipboard brute stacks of the MCS data reveal evidence for reactivation of abyssal hill faulting in the plate interior far from the trench. Ridgeward-dipping lower crustal reflectors are observed, similar to those observed in mature Pacific crust elsewhere, as well as conjugate reflectivity near the deformation front along the Oregon transect. Bright intracrustal reflectivity is also observed along the trench-parallel transect with marked changes in reflectivity along the Oregon and Washington margins. Initial inspection of the OBS record sections indicate good quality data with the expected oceanic crustal and upper mantle P-wave arrivals: Ps and Pg refractions through sedimentary and igneous layers, respectively, PmP wide-angle reflections from the crust-mantle transition zone, and Pn upper mantle refractions. The Pg-PmP-Pn triplication is typically observed at 40-50 km source-receiver offsets. Pn characteristics show evidence for upper mantle azimuthal anisotropic propagation: along the plate transects Pn is typically weaker and difficult to observe beyond ~80 km offsets, while along the trench-parallel transect Pn arrivals have higher amplitude and are easily observed up to source-receiver offsets of 160-180 km. An overview on the Cascadia Ridge to Trench data acquisition program and preliminary results will be presented.

A Review of Geophysical Constraints on the Deep Structure of the Tibetan Plateau, the Himalaya and the Karakoram, and their Tectonic Implications

NASA Astrophysics Data System (ADS)

Molnar, P.

1988-09-01

The Tibetan Plateau, the Himalaya and the Karakoram are the most spectacular consequences of the collision of the Indian subcontinent with the rest of Eurasia in Cainozoic time. Accordingly, the deep structures beneath them provide constraints on both the tectonic history of the region and on the dynamic processes that have created these structures. The dispersion of seismic surface waves requires that the crust beneath Tibet be thick: nowhere less than 50 km, at least 65 km, in most areas, but less than 80 km in all areas that have been studied. Wide-angle reflections of P-waves from explosive sources in southern Tibet corroborate the existence of a thick crust but also imply the existence of marked lateral variations in that thickness, or in the velocity structure of the crust. Thus isostatic compensation occurs largely by an Airy-type mechanism, unlike that, for instance, of the Basin and Range Province of western North America where a hot upper mantle buoys up a thin crust. The P-wave and S-wave velocities in the uppermost mantle of most of Tibet are relatively high and typical of those of Precambrian shields and stable platforms: Vp = 8.1 km s-1 or higher, and Vs≈ 4.7 km s-1. Travel times and waveforms of S-waves passing through the uppermost mantle of much of Tibet, however, require a much lower average velocity in the uppermost mantle than that of the Indian, or other, shields. They indicate a thick low-velocity zone in the upper mantle beneath Tibet, reminiscent of tectonically active regions. These data rule out a shield structure beneath northern Tibet and suggest that if such a structure does underlie part of the plateau, it does so only beneath the southern part. Lateral variations in the upper-mantle structure of Tibet are apparent from differences in travel times of S-waves from earthquakes in different parts of Tibet, in the attenuation of short-period phases, Pn and Sn, that propagate through the uppermost mantle of Tibet, and in surface-wave dispersion for different paths. The notably lower velocities and the greater attenuation in the mantle of north--central Tibet than elsewhere imply higher temperatures there and are consistent with the occurrence of active and young volcanism in roughly the same area. Surface-wave dispersion across north--central Tibet also requires a thinner crust in that area than in most of the plateau. Consequently the relatively uniform height of the plateau implies that isostatic compensation in the north--central part of Tibet occurs partly because the density of the relatively hot material in the upper mantle is lower than that elsewhere beneath Tibet, the mechanism envisioned by Pratt. Several seismological studies provide evidence consistent with a continuity of the Indian Shield, and its cold thick lithosphere, beneath the Himalaya. Fault-plane solutions and focal depths of the majority of moderate earthquakes in the Himalaya are consistent with their occurring on the top surface of the gently flexed, intact Indian plate that has underthrust the Lesser Himalaya roughly 80-100 km or more. P-waves from explosions in southern Tibet and recorded in Nepal can be interpreted as wide-angle reflections from this fault zone. P-wave delays across the Tarbela network in Pakistan from distant earthquakes indicate a gentle dip of the Moho beneath the array without pronounced later variations in upper-mantle structure. High Pn and Sn velocities beneath the Himalaya and normal to early S-wave arrival times from Himalayan earthquakes recorded at teleseismic distances are consistent with Himalaya being underlain by the same structure that underlies India. Results from explosion seismology indicate an increase in crustal thickness from the Indo--Gangetic Plain across the Himalaya to southern Tibet, but Hirn, Lepine, Sapin and their co-workers inferred that the depth of the Moho does not increase smoothly northward, as it would if the Indian Shield had been underthrust coherently beneath the Himalaya. They interpreted wide-angle reflections as evidence for steps in the Moho displaced from one another on southward-dipping faults. Although I cannot disprove this interpretation, I think that one can recognize a sequence of signals on their wide-angle reflection profiles that could be wide-angle reflections from a northward-dipping Moho. Gravity anomalies across the Himalaya show that both the Indo--Gangetic Plain and the Himalaya are not in local isostatic equilibrium. A mass deficit beneath the plain is apparently caused by the flexure of the Indian Shield and by the low density of the sedimentary rock in the basin formed by the flexure. The mass excess in the Himalaya seems to be partly supported by the strength of the Indian plate, for which the flexural rigidity is particularly large. An increase in the Bouguer gravity gradient from about 1 mGal km-1 (1 mGal = 10-3 cm s-2) over the Indo--Gangetic Plain to 2 mGal km-1 over the Himalaya implies a marked steepening of the Moho, and therefore a greater flexure of the Indian plate, beneath the Himalaya. This implies a northward decrease in the flexural rigidity of the part of the Indian plate underlying the range. Nevertheless, calculations of deflections of elastic plates with different flexural rigidities and flexed by the weight of the Himalaya show larger deflections and yield more negative gravity anomalies than are observed. Thus, some other force, besides the flexural strength of the plate, must contribute to the support of the range. A bending moment applied to the end of the Indian plate could flex the plate up beneath the range and provide the needed support. The source of this moment might be gravity acting on the mantle portion of the subducting Indian continental lithosphere with much or all of the crust detached from it. Seismological studies of the Karakoram are consistent with its being underlain by particularly cold material in the upper mantle. Intermediate-depth earthquakes occur between depths of 70 and 100 km but apparently do not define a zone of subducted oceanic lithosphere. Rayleigh-wave phase velocities are particularly high for paths across this area and imply high shear wave velocities in the upper mantle. Isostatic gravity anomalies indicate a marked low of 70 mGal over the Karakoram, which could result from a slightly thickened crust pulled down by the sinking of cold material beneath it. Geophysical constraints on the structure of Tibet, the Himalaya and the Karakoram are consistent with a dynamic uppermost mantle that includes first, the plunging of cold material into the asthenosphere beneath southern Tibet and the Karakoram, as the Indian plate slides beneath the Himalaya, and second, an upwelling of hot material beneath north--central Tibet. The structure is too poorly resolved to require such dynamic flow, but the existence for both a hot uppermost mantle beneath north--central Tibet and a relatively cold uppermost mantle beneath southern Tibet and the Karakoram seem to be required. Both group and phase velocities of Rayleigh waves and Love waves are delayed along paths crossing Tibet. The low velocities require a crustal thickness in excess of 50 km, and for most regions in excess of 60 km. Crustal thicknesses in excess of 80 km can be ruled out for all paths studied, and for most of Tibet, a crustal thickness of 65-70 km seems required. Clear evidence for lateral heterogeneity beneath Tibet is provided not only by body waves (discussed below) but also by surface waves (Brandon & Romanowicz 1986), which show an area of lower uppermost shear-wave velocity and thinner crust in north--central Tibet than elsewhere in the plateau. These variations might explain the differences in group velocities measured by different workers, and the different structures that they deduced, but if so, they also render the regionalization of surface-wave dispersion into arbitrary tectonic provinces risky. Although Rayleigh-wave phase velocities can resolve large differences in upper-mantle velocities for regions the size of Tibet, constraints on these velocities are best derived from body waves. Thus, with the exceptions of Brandon & Romanowicz's (1986) detailed investigation of north--central Tibet, the study of southernmost Tibet by Jobert et al. (1985) and that of Romanowicz (1982) for the northeasternmost part of the plateau, I do not think that surface waves have placed an important bound on the velocity in the upper mantle beneath Tibet. The seismic data are broadly consistent with partial melting of the uppermost mantle of north--central Tibet, where recent volcanism has been observed. Correspondingly, there is no suggestion of such low velocities, and such high temperatures, in the mantle elsewhere beneath Tibet, for which late-Cainozoic volcanism has not been reported. The results are also consistent with a slightly thinner crust in north--central Tibet than farther south, suggesting that both Airy and Pratt isostasy share compensation for north--central Tibet's great height. Finally, the average shear-wave velocity in the upper mantle of southern Tibet seems to be higher than that in northern Tibet, but neither is the degree of difference well determined, nor is the location of the transition from one to the other well mapped.

Identifying Moho depths and velocity anomalies in the uppermost mantle of the Mississippi Embayment from Pn tomography and anisotropy studies

NASA Astrophysics Data System (ADS)

Basu, U.; Powell, C. A.

2017-12-01

Lateral depth variations of the Mohorovicic discontinuity, Pn velocities, and anisotropy features at uppermost mantle depths below the central U.S. are determined using Pn tomography. Excellent raypath coverage throughout the northern Mississippi Embayment (ME) is obtained using the NELE (Northern Embayment Lithosphere Experiment) and US TA (Transportable Array) datasets. High Pn velocities are present below the northern portion of the Reelfoot Rift and the New Madrid seismic zone. Prominent regions of low velocity are present to the east and north of the ME, in agreement with recent teleseismic tomography studies indicating the presence of low P- and S-wave velocities in the uppermost mantle. A prominent region of low velocity coincides with the southwestern portion of the Illinois Basin. Higher velocities are located west of the Illinois Basin and west of the Ozark Plateau. Crustal thicknesses obtained from the Pn station delays indicate thinner crust in the southern Coastal Plain and ME and thicker crust north of the ME. Strong Pn anisotropy and rotation of the fast directions are associated with the northern ME. Fast directions differ from present absolute plate motion directions and from fast directions determined from SKS splitting, suggesting the presence of multiple anisotropic layers. Parameter errors estimated using the bootstrap method are all less than 0.1 km/s for velocity and magnitude of the anisotropy.

Kinematic signature of India/Australia plates break-up

NASA Astrophysics Data System (ADS)

Iaffaldano, G.; Bunge, H.

2008-12-01

The paradigm of Plate Tectonics states that the uppermost layer of the Earth is made of a number of quasi- rigid blocks moving at different rates in different directions, while most of the deformation is focused along their boundaries. Perhaps one of the most interesting and intriguing processes in Plate Tectonics is the generation of new plate boundaries. The principle of inertia implies that any such event would invariably trigger changes in plate motions, because the budget of mantle basal-drag and plate-boundary forces would be repartitioned. A recent episode is thought to have occurred in the Indian Ocean, where a variety of evidences - including localized seismicity along the Nienty East Ridge, compression-generated unconformities of ocean-floor sediments, and identified paleomagnetic isochrones - suggest the genesis of a boundary separating the India and Australia plates. Here we use global numerical models of the coupled mantle/lithosphere system to show for the first time that an event of separation between India and Australia, having occurred sometime between 11 and 8 Myrs ago, has left a distinct signature in the observed record of plate motions. Specifically, while motions of India and Australia relative to fixed Eurasia are almost indistinguishable prior to 11 Myrs ago, their convergence to Eurasia since then differs significantly, by as much as 2 cm/yr. Finally, we speculate about possible causes for the separation between India and Australia plates.

Kinematic signature of India/Australia plates break-up

NASA Astrophysics Data System (ADS)

Iaffaldano, G.; Bunge, H.-P.

2009-04-01

The paradigm of Plate Tectonics states that the uppermost layer of the Earth is made of a number of quasi-rigid blocks moving at different rates in different directions, while most of the deformation is focused along their boundaries. Perhaps one of the most interesting and intriguing processes in Plate Tectonics is the generation of new plate boundaries. The principle of inertia implies that any such event would invariably trigger changes in plate motions, because the budget of mantle basal-drag and plate-boundary forces would be repartitioned. A recent episode is thought to have occurred in the Indian Ocean, where a variety of evidences - including localized seismicity along the Nienty East Ridge, compression-generated unconformities of ocean-floor sediments, and identified paleomagnetic isochrones - suggest the genesis of a boundary separating the India and Australia plates. Here we use global numerical models of the coupled mantle/lithosphere system to show for the first time that an event of separation between India and Australia, having occurred sometime between 11 and 8 Myrs ago, has left a distinct signature in the observed record of plate motions. Specifically, while motions of India and Australia relative to fixed Eurasia are almost indistinguishable prior to 11 Myrs ago, their convergence to Eurasia since then differs significantly, by as much as 2 cm/yr. Finally, we speculate about possible causes for the separation between India and Australia plates.

High resolution image of uppermost mantle beneath NE Iran continental collision zone

NASA Astrophysics Data System (ADS)

Motaghi, K.; Tatar, M.; Shomali, Z. H.; Kaviani, A.; Priestley, K.

2012-10-01

We invert 3775 relative P wave arrival times using the ACH damped least square method of Aki et al. (1977) to study upper mantle structure beneath the NE Iran continental collision zone. The data for this study were recorded by 17 three component broad-band stations operated from August 2006 to February 2008 along a profile from the center of Iranian Plateau, near Yazd, to the northeastern part of Iran on the Turan Platform just north of the Kopeh Dagh Mountains. The results confirm the previously known low velocity upper mantle beneath Central Iran. Our tomographic model reveals a deep high velocity anomaly. The surficial expressions of this anomaly are between the Ashkabad and Doruneh Faults, where the resolution and ray coverage are good. A transition zone in uppermost mantle is recognized under the Binalud foreland that we interpreted as suture zone between Iran and Turan platform. Our results indicate that Atrak Valley which is the boundary between the Binalud and Kopeh Dagh Mountains can be considered as the northeastern suture of the Iranian Plateau where Eurasia and Turan Platform under-thrust beneath the Binalud range and Central Iran.

Pre-subduction metasomatic enrichment of the oceanic lithosphere induced by plate flexure

NASA Astrophysics Data System (ADS)

Pilet, S.; Abe, N.; Rochat, L.; Kaczmarek, M.-A.; Hirano, N.; Machida, S.; Buchs, D. M.; Baumgartner, P. O.; Müntener, O.

2016-12-01

Oceanic lithospheric mantle is generally interpreted as depleted mantle residue after mid-ocean ridge basalt extraction. Several models have suggested that metasomatic processes can refertilize portions of the lithospheric mantle before subduction. Here, we report mantle xenocrysts and xenoliths in petit-spot lavas that provide direct evidence that the lower oceanic lithosphere is affected by metasomatic processes. We find a chemical similarity between clinopyroxene observed in petit-spot mantle xenoliths and clinopyroxene from melt-metasomatized garnet or spinel peridotites, which are sampled by kimberlites and intracontinental basalts respectively. We suggest that extensional stresses in oceanic lithosphere, such as plate bending in front of subduction zones, allow low-degree melts from the seismic low-velocity zone to percolate, interact and weaken the oceanic lithospheric mantle. Thus, metasomatism is not limited to mantle upwelling zones such as mid-ocean ridges or mantle plumes, but could be initiated by tectonic processes. Since plate flexure is a global mechanism in subduction zones, a significant portion of oceanic lithospheric mantle is likely to be metasomatized. Recycling of metasomatic domains into the convecting mantle is fundamental to understanding the generation of small-scale mantle isotopic and volatile heterogeneities sampled by oceanic island and mid-ocean ridge basalts.

Magmatic arc structure around Mount Rainier, WA, from the joint inversion of receiver functions and surface wave dispersion

NASA Astrophysics Data System (ADS)

Obrebski, Mathias; Abers, Geoffrey A.; Foster, Anna

2015-01-01

The deep magmatic processes in volcanic arcs are often poorly understood. We analyze the shear wave velocity (VS) distribution in the crust and uppermost mantle below Mount Rainier, in the Cascades arc, resolving the main velocity contrasts based on converted phases within P coda via source normalization or receiver function (RF) analysis. To alleviate the trade-off between depth and velocity, we use long period phase velocities (25-100 s) obtained from earthquake surface waves, and at shorter period (7-21 s) we use seismic noise cross correlograms. We use a transdimensional Bayesian scheme to explore the model space (VS in each layer, number of interfaces and their respective depths, level of noise on data). We apply this tool to 15 broadband stations from permanent and Earthscope temporary stations. Most results fall into two groups with distinctive properties. Stations east of the arc (Group I) have comparatively slower middle-to-lower crust (VS = 3.4-3.8 km/s at 25 km depth), a sharp Moho and faster uppermost mantle (VS = 4.2-4.4 km/s). Stations in the arc (Group II) have a faster lower crust (VS = 3.7-4 km/s) overlying a slower uppermost mantle (VS = 4.0-4.3 km/s), yielding a weak Moho. Lower crustal velocities east of the arc (Group I) most likely represent ancient subduction mélanges mapped nearby. The lower crust for Group II ranges from intermediate to felsic. We propose that intermediate-felsic to felsic rocks represent the prearc basement, while intermediate composition indicates the mushy andesitic crustal magmatic system plus solidified intrusion along the volcanic conduits. We interpret the slow upper mantle as partial melt.

Major Element Geochemistry of Peridotites from Santa Elena Ophiolite Complex, NW Costa Rica and Their Tectonic Implications

NASA Astrophysics Data System (ADS)

Wright, S.; Snow, J. E.; Gazel, E.; Sisson, V.

2010-12-01

The Santa Elena Ophiolite Complex (SEOC) is located on the west coast of Northern Costa Rica, near the Nicaraguan border. It consists primarily of preserved oceanic crustal rocks and underlying upper mantle thrust onto an accretionary complex. The petrogenesis and tectonic origin of this complex have widely been interpreted to be either a preserved mantle portion of the Caribbean Large Igneous Province (CLIP) as it drifted between North and South America from the Galapagos hotpot into the present day Caribbean Ocean around 80 Ma or as the mantle section to the nearby Nicoya complex. Previous structural work suggests that SEOC is a supra-subduction complex, not related to the CLIP or Nicoya. Our preliminary results agree. Mantle peridotites collected from the Santa Elena Ophiolite Complex consist primarily of spinel lherzolite (61 %) with minor amounts of harzburgite and dunite (22 % and 16 % respectively). Spinel Cr# [molar Cr / (Cr+Al)*100] is widely accepted to constrain mantle partial melting and lithospheric melt stagnation. Cr# of spinels within Santa Elena lherzolites fall between 12 and 35, suggesting an extent of 3 % to 13 % partial melting. Cr# of harzburgites range from 35 to 39, suggesting 13 % to 14 % partial melting. This range of partial melting suggests only modest depletion of this exposed portion of the ancient uppermost mantle. TiO2 concentrations of the lherzolite and harzburgite range from 0.004% to 0.128%, with the exception of one sample, SE10 - 17 (0.258%), and fall within the normal melting trend for mantle peridotites. The presence of dunite indicates that melt flow and associated melt - rock reaction with the surrounding peridotite took place within this portion of the mantle. A Cr# of 84.5 from one of these dunite samples indicate that significant melt rock reaction with refractory melts took place. Such results are rarely found in mid-ocean ridge abyssal peridotite settings, and are currently found primarily in forearc tectonic settings. However, due to the overall "normal" TiO2 concentrations in all but one spinel peridotite requires that if melt flow did occur, that the melt be nearly depleted in titanium. The relatively low Cr#'s and TiO2 concentrations of spinel in these peridotites that suggest low degrees of partial melting along with the paleo presence of melt flow and melt-rock reaction by low titanium melts, such as boninites, point toward a young fore-arc model for the tectonic origin of this ophiolite body rather than a preserved mantle portion of the CLIP. Additionally, two lines of evidence suggest SEOC was emplaced prior to the collision of the CLIP with North and South America. The SEOC is 1) capped by a Campanian (83.5 - 70.6 Ma) rudist limestone and 2) lies uncomformably atop Cenomanian (93.6 - 99.6 Ma) radiolarite beds. This suggests that the mantle portion of the SEOC was emplaced and exposed at the Caribbean ocean floor prior to the Late Cretaceous (Campanian), but no earlier than the Cenomanian. This combined tectonic and geochemical evidence suggests SEOC may be a portion of the proto-arc that existed between the Americas in the Cretaceous prior to assault by the CLIP.

Melt-rock interactions and fabric development of peridotites from North Pond in the Kane area, Mid-Atlantic Ridge: Implications of microstructural and petrological analyses of peridotite samples from IODP Hole U1382A

NASA Astrophysics Data System (ADS)

Harigane, Yumiko; Abe, Natsue; Michibayashi, Katsuyoshi; Kimura, Jun-Ichi; Chang, Qing

2016-06-01

North Pond is an isolated sedimentary pond on the western flank of the Kane area along the Mid-Atlantic Ridge. Drill-hole U1382A of IODP Expedition 336 recovered peridotite and gabbro samples from a sedimentary breccia layer in the pond, from which we collected six fresh peridotite samples. The peridotite samples came from the southern slope of the North Pond where an oceanic core complex is currently exposed. The samples were classified as spinel harzburgite, plagioclase-bearing harzburgite, and a vein-bearing peridotite that contains tiny gabbroic veins. No obvious macroscopic shear deformation related to the formation of a detachment fault was observed. The spinel harzburgite with a protogranular texture was classified as refractory peridotite. The degree of partial melting of the spinel harzburgite is estimated to be ˜17%, and melt depletion would have occurred at high temperatures in the uppermost mantle beneath the spreading axis. The progressive melt-rock interactions between the depleted spinel harzburgite and the percolating melts of Normal-Mid Ocean Ridge Basalt (N-MORB) produced the plagioclase-bearing harzburgite and the vein-bearing peridotite at relatively low temperatures. This implies that the subsequent refertilization occurred in an extinct spreading segment of the North Pond after spreading at the axis. Olivine fabrics in the spinel and plagioclase-bearing harzburgites are of types AG, A, and D, suggesting the remnants of a mantle flow regime beneath the spreading axis. The initial olivine fabrics appear to have been preserved despite the later melt-rock interactions. The peridotite samples noted above preserve evidence of mantle flow and melt-rock interactions beneath a spreading ridge that formed at ˜8 Ma.

Direct evidence of hydration into mantle during shearing below a trasform fault: Prince Edward transform fault, Southwest Indian Ridge

NASA Astrophysics Data System (ADS)

Michibayashi, K.; Kakihata, Y.; Dick, H. J.

2017-12-01

Southwest Indian Ridge (SWIR) is located to the southwest of Rodriguez Triple Junction, where three Indian ocean ridges meet (Zhou & Dick, 2013, Nature). SWIR is one of the slowest spreading ocean ridges in the world. In this study, we studied microstructural development of 21 peridotite samples obtained from Prince Edward transform fault of SWIR by PROTEA5 cruise in 1983. The peridotites consist dominantly of olivine, orthopyroxene and clinopyroxene with minor amounts of amphibole and plagioclase as well as secondary minerals such as serpentine and magnetite. The peridotites were classified into four groups based on their microstructures: 3 ultramylonites mostly consisting of extremely fine crystals (3-5µm), 13 heterogeneous tectonites consisting of coarse-grained crystals and fine-grained matrix, 1 cataclasite and 4 intensely serpentinized peridotites. Olivine Mg# is 0.90-0.91 and spinel Cr# is 0.1-0.35. Amphibole crystals have chemical compositions of tremolite and magnesio-hornblende and they were intensely deformed within the ultramylonites and the heterogeneous tectonites, indicating that they have occurred before or during intense shearing in mantle. Moreover, extremely fine grain sizes of olivine and microboudin textures in both pyroxene and spinel crystals suggest that these peridotites have been sheared under high stress conditions. Furthermore, olivine crystal-fabrics within the amphibole bearing peridotites have B and E types that could be developed under hydrous conditions, whereas olivine fabrics within the other peridotites have A and D types that could be developed under anhydrous conditions (Karato et al., 2008, Annu. Rev. Earth Planet. Sci.). Consequently, the petrophysical characteristics of peridotites in this study indicate that the uppermost mantle below the Prince Edward transform fault has been locally but intensely hydrated during shearing due to transform movement.

Dependence of Subduction Zone seismicity on Strain-Rate-Dependent Critical Homologous Temperature

NASA Astrophysics Data System (ADS)

Davis, P. M.

2016-12-01

Earthquakes mainly occur in crust or mantle that is below a critical temperature for the tectonic strain-rate, such that stress builds up to the breaking point before it can relax due to creep. Then long-range stress correlation gives rise to power law seismicity with large events. The limiting temperature depends on pressure, which is taken into account by finding a critical homologous temperature THc=T/TM above which earthquakes are rarely observed. We find that THc for ocean plates is ˜0.55. For California earthquakes, it is also close to 0.55. The uppermost mantle layer of oceanic plates of thickness ˜50 km is composed of harzburgite and depleted peridotite from which basalt has been removed to form ocean crust. Thus it has a higher melting temperature than the peridotite of the surrounding mantle, or the lower halves of plates. Thicknesses of seismicity in deep subduction zones, determined from 2D polynomial fits to a relocated catalog, are ˜50 km, which suggests that the earthquake channel is confined to this layer. We construct models to find homologous temperatures in slabs, and find that seismicity thicknesses are also, on average, confined to TH ≤ 0.55 ± 0.05. The associated rheology is compared with that obtained from flexure models of ocean lithosphere. The brittle-ductile transition occurs where viscosity drops from high values in the cold cores of slabs to values of 1022 to $1023 Pa s, i.e., where creep strain-rates become comparable to tectonic rates. The cutoff for deep earthquakes is not sharp. However they appear unlikely to occur if homologous temperature is high TH>0.55. Exceptions to the rule are anomalously deep earthquakes such as those beneath the Iceland and the Hawaiian hotspots, and the Newport Inglewood Fault. These are smaller events with short-range stress correlation, and can be explained if strain-rates are 2 to 3 orders of magnitude higher than those associated with earthquakes located where TH ≤0.55. We conclude that the brittle-ductile transition corresponds to the transition from long-range to short-range stress correlation.

Pn anisotropy beneath the South Island of New Zealand and implications for distributed deformation in continental lithosphere

NASA Astrophysics Data System (ADS)

Collins, John A.; Molnar, Peter

2014-10-01

Pn travel times from regional earthquakes recorded both by stations on New Zealand and by ocean bottom seismographs deployed offshore indicate anisotropy in the uppermost mantle beneath the region. The largest anisotropy of ~8% (±2%, 1σ) lies beneath the deforming part of the South Island to just off its West Coast, a zone roughly 100-200 km wide. The fastest propagation is aligned N60°E (±3°), essentially parallel to the largely strike-slip relative plate motion since 20 Ma, also ~ N60°E. The magnitude of anisotropy decreases abruptly northwest and southeast of this zone, and on the southeast side of the island, the orientation of fastest propagation is between N32°W and N-S. The ~ N60°E orientation of fast propagation is consistent with finite strain within the uppermost part of the mantle lithosphere if the measured 850 km of displacement of the Pacific plate past the Australia plate is spread over a region with a width of 100-200 km. The agreement of this orientation of fast propagation with the orientation or relative plate motion suggests the possibility of but does not require some dynamic recrystallization in rock as cold as 500-800°C, where Peierls creep seems to be the likely deformation mechanism. Such a strain distribution matches deformation of a thin viscous sheet that obeys a constitutive relationship of the form ɛ>˙ ~ τn, where ɛ>˙ is the average strain rate and τ is the operative deviatoric stress, with an average value of n ≈ 3-10. Presumably, the NW-SE fast propagation in the region southeast of the island results from strain that precedes the Cenozoic deformation that has shaped the island.

Upper-mantle seismic structure in a region of incipient continental breakup: northern Ethiopian rift

NASA Astrophysics Data System (ADS)

Bastow, Ian D.; Stuart, Graham W.; Kendall, J.-Michael; Ebinger, Cynthia J.

2005-08-01

The northern Ethiopian rift forms the third arm of the Red Sea, Gulf of Aden triple junction, and marks the transition from continental rifting in the East African rift to incipient oceanic spreading in Afar. We determine the P- and S-wave velocity structure beneath the northern Ethiopian rift using independent tomographic inversion of P- and S-wave relative arrival-time residuals from teleseismic earthquakes recorded by the Ethiopia Afar Geoscientific Lithospheric Experiment (EAGLE) passive experiment using the regularised non-linear least-squares inversion method of VanDecar. Our 79 broad-band instruments covered an area 250 × 350 km centred on the Boset magmatic segment ~70 km SE of Addis Ababa in the centre of the northern Ethiopian rift. The study area encompasses several rift segments showing increasing degrees of extension and magmatic intrusion moving from south to north into the Afar depression. Analysis of relative arrival-time residuals shows that the rift flanks are asymmetric with arrivals associated with the southeastern Somalian Plate faster (~0.65 s for the P waves; ~2 s for the S waves) than the northwestern Nubian Plate. Our tomographic inversions image a 75 km wide tabular low-velocity zone (δVP~-1.5 per cent, δVS~-4 per cent) beneath the less-evolved southern part of the rift in the uppermost 200-250 km of the mantle. At depths of >100 km, north of 8.5°N, this low-velocity anomaly broadens laterally and appears to be connected to deeper low-velocity structures under the Afar depression. An off-rift low-velocity structure extending perpendicular to the rift axis correlates with the eastern limit of the E-W trending reactivated Precambrian Ambo-Guder fault zone that is delineated by Quaternary eruptive centres. Along axis, the low-velocity upwelling beneath the rift is segmented, with low-velocity material in the uppermost 100 km often offset to the side of the rift with the highest rift flank topography. Our observations from this magmatic rift zone, which is transitional between continental and oceanic rifting, do not support detachment fault models of lithospheric extension but instead point to strain accommodation via magma assisted rifting.

Crustal structure of the Sunda-Banda arc transition: results from marine geophysical investigations offshore eastern Indonesia

NASA Astrophysics Data System (ADS)

Planert, L.; Shulgin, A.; Kopp, H.; Lueschen, E.; Mueller, C.; Flueh, E.; Djajadihardja, Y.; Engels, M.

2009-04-01

The Sunda-Banda arc transition, the easternmost portion of the Indonesian convergent margin, presents a probably unique natural laboratory to study lower plate variability and related upper plate deformation in the so-called ‘subduction factory' for a deeper understanding of forearc evolution. In neighboring margin segments, we can observe strong changes of the incoming plate (transition from an oceanic to a continental lower plate, increasing plate age to the East, presence/absence of an oceanic plateau, variability in plate roughness) as well as a wide range of corresponding forearc structures, including large sedimentary basins and an accretionary prism/outer arc high of variable size and shape. During RV Sonne cruise SO190 in 2006 (SINDBAD: Seismic and Geoacoustic Investigations along the Sunda-Banda Arc Transition), we acquired a combination of seismic wide-angle OBH/OBS refraction, multichannel streamer and gravity data in order to study the seismic velocity structure of the subducting crust and the overriding island arc along a number of trench normal corridors located between 113°E and 121°E. Additionally, a number of trench parallel profiles were conducted which mainly focus on the internal structure of the large sedimentary basins and which were also intended for further clarifying the type of underlying forearc crust and mantle respectively. We used a tomographic approach for refracted and reflected phases to obtain seismic velocity models which again were used for prestack depth-migration of the MCS data. In turn, we incorporated the highly resolved sedimentary portions as a priori structure in our tomography. The results show the seismic velocity structure of the incoming plate, starting 100 km seaward of the trench, and the adjoining forearc down to depths of 20-28 km, i.e. well into the upper mantle, and at the same time fit the gravity data very well, using simple velocity-density relations. In the Argo abyssal plain, the models show 8.0-8.5 km thick oceanic crust. The velocities in the crust and uppermost mantle are reduced within distances of ~50 km seaward of the trench, which coincides with the onset of normal faulting on the incoming oceanic plate. Anomalously low mantle velocities of 7.5 km/s directly beneath the Moho are possibly due to the intrusion of seawater and subsequent serpentinisation of mantle peridotite. Landward of the trench in the outer arc high, velocities do not exceed 5.5 km/s down to the top of the subducting slab, which can be traced over ~70 km length beneath the forearc down to ~13 km depth. The plate boundary is of irregular shape, obviously imprinted by the complex deformation of the oceanic basement prior to subduction, which is further amplified as response to thrusting/downbending of the dissected oceanic blocks. Offshore Lombok island, our models reveal the geometry of the Lombok basin as well as the forearc Moho in ~16 km depth. Reduced upper mantle velocities suggest a hydrated shallow mantle wedge for this corridor. Further east offshore Sumba island, where the Java trench terminates and the transition to the collisional regime further east occurs, our models show a subducting oceanic plate of similar thickness and structure. But different to the situation offshore Lombok, we find no evidence for a shallow mantle wedge beneath the forearc; crustal-type velocities are found down to depths of ~20 km. The different forearc regime is most likely related to the collision with the Sumba block. Our results give a detailed view into the complex structure in both the deeper and shallower portions of this convergent margin.

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

PubMed Central

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

2015-01-01

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

Seismic anisotropy of 70 Ma Pacific-plate upper mantle

NASA Astrophysics Data System (ADS)

Mark, H. F.; Lizarralde, D.; Collins, J. A.; Miller, N. C.; Hirth, G.; Gaherty, J. B.; Evans, R. L.

2017-12-01

We present a new measurement of seismic anisotropy and velocity gradients in the Pacific-plate upper mantle based on data from the NoMelt experiment. The seismic velocity structure of oceanic lithosphere reflects the processes involved in its formation at mid-ocean ridges and subsequent evolution off-axis. Increasing mantle depletion with depth due to melt extraction predicts negative velocity gradients, as does cooling with age. Alignment of olivine by corner flow predicts azimuthal anisotropy. Some models predict the strength of anisotropy should decrease with depth. Measurements of uppermost mantle velocities have not fully verified these predictions. Observations of direct Pn phases demonstrate that positive velocity gradients exist; and anisotropy measurements, while consistent with strain-induced olivine alignment, vary widely and generally suggest weaker fabric development than is observed in ophiolite samples. These discrepancies raise questions about the extent to which mantle structure evolves through time due to processes such as cracking and alteration, and hinder the use of seismic measurements to make more detailed inferences on aspects of lithospheric formation processes. We have measured anisotropy and vertical velocity gradients to 10 km below the Moho on 70 Ma lithosphere between the Clarion and Clipperton fracture zones. The lithosphere at the study site has not been obviously affected by tectonic or magmatic events since its formation. We find 6.2% anisotropy at the Moho with a mean velocity of 8.14 km/s and the fast direction parallel to paleospreading. Velocity gradients are estimated at 0.02 km/s/km in the fast direction and near 0 km/s/km in the slow direction. The gradient estimates can be explained by aligned microcracks oriented perpendicular to spreading that close with depth. Cracks are expected to close by 10 km below the Moho. At that depth the strength of anisotropy increases to 9%, close to the strength estimated from ophiolite fabrics. These results are consistent with observed olivine fabrics and the predicted effects of lithospheric formation processes, and suggest that lithospheric evolution is modest even at 70 Ma, involving microcracks oriented by a stress field consistent with thermal contraction.

3D Crust and Uppermost Mantle Structure beneath Tian Shan Region from ambient noise and earthquake surface waves

NASA Astrophysics Data System (ADS)

Xiao, X.; Wen, L.

2017-12-01

As a typical active intracontinental mountain range in Central Asia, Tian Shan Mt serves as the prototype in studying geodynamic processes and mechanism of intracontinental mountain building. We study 3D crust and the uppermost mantle structure beneath Tian Shan region using ambient noise and earthquake surface waves. Our dataset includes vertical component records of 62 permanent broadband seismic stations operated by the Earthquake Administration of China. Firstly, we calculate two-year stacked Cross-Correlation Functions (CCFs) of ambient noise records between the stations. The CCFs are treated as the Empirical Green's Functions (EGFs) of each station pair, from which we measured phase velocities of fundamental-mode Rayleigh wave in the period of 3-40 s using a frequency-time analysis method. Secondly, we collect surface wave data from tele-seismic events with Mw > 5.5 and depth shallower than 200 km and measure phase velocities of the fundamental-mode of Rayleigh wave in the period of 30-150 s using a two-station method. Finally, we combine the phase velocity measurements from ambient noise and earthquake surface waves, obtain lateral isotropic phase velocity maps at different periods based on tomography and invert a 3D Vsv model of crust and uppermost mantle down to about 150 km using a Monte Carlo Inversion method. We will discuss our inversion results in detail, as well as their implications to the tectonics in the region.

Crustal and Uppermost Mantle Structure across the Tibet-Qinling Transition Zone in NE Tibet: Implications for Material Extrusion of the Tibetan Plateau

NASA Astrophysics Data System (ADS)

Ye, Z.; Li, J.; Gao*, R.; Song, X.; Li, Q.; Li, Y.; Huang, X.; Xiong, X.; Li, W.; WANG, Y.

2017-12-01

Based on a dense linear seismic array traversing across the eastern margin of the Tibetan plateau to the Qinling belt, we conducted joint inversion of receiver functions and surface wave dispersions under constraints of P-wave velocity and derived a crustal and uppermost mantle Vs profile simultaneously with a Vp/Vs ratio profile. Our observations indicate that the Qinling belt, which shows ratio Vp/Vs<1.8 indicative of intermediate-to-felsic components in the lower crust, is currently not acting as a channel accommodating extrusion of the mid-lower crustal flow; and extrusion of Tibet's ductile mantle flow through the Qinling belt as a channel would only be feasible in the sub-lithosphere depth (asthenosphere). Our results suggest that ductile material extrusion of the mid-lower crustal flow accompanied with fault-related tectonics and gravitational buoyancy resulted from lithospheric detachment (triggered by the asthenospheric flow) may jointly work on the plateau uplift and expansion in this Tibet-Qinling transition zone. Corresponding Author: R.Gao, ruigao126@126.com

Gravity Spectra from the Density Distribution of Earth's Uppermost 435 km

NASA Astrophysics Data System (ADS)

Sebera, Josef; Haagmans, Roger; Floberghagen, Rune; Ebbing, Jörg

2018-03-01

The Earth masses reside in a near-hydrostatic equilibrium, while the deviations are, for example, manifested in the geoid, which is nowadays well determined by satellite gravimetry. Recent progress in estimating the density distribution of the Earth allows us to examine individual Earth layers and to directly see how the sum approaches the observed anomalous gravitational field. This study evaluates contributions from the crust and the upper mantle taken from the LITHO1.0 model and quantifies the gravitational spectra of the density structure to the depth of 435 km. This is done without isostatic adjustments to see what can be revealed with models like LITHO1.0 alone. At the resolution of 290 km (spherical harmonic degree 70), the crustal contribution starts to dominate over the upper mantle and at about 150 km (degree 130) the upper mantle contribution is nearly negligible. At the spatial resolution <150 km, the spectra behavior is driven by the crust, the mantle lid and the asthenosphere. The LITHO1.0 model was furthermore referenced by adding deeper Earth layers from ak135, and the gravity signal of the merged model was then compared with the observed satellite-only model GOCO05s. The largest differences are found over the tectonothermal cold and old (such as cratonic), and over warm and young areas (such as oceanic ridges). The misfit encountered comes from the mantle lid where a velocity-density relation helped to reduce the RMS error by 40%. Global residuals are also provided in terms of the gravitational gradients as they provide better spatial localization than gravity, and there is strong observational support from ESA's satellite gradiometry mission GOCE down to the spatial resolution of 80-90 km.

Variable styles of rifting expressed in crustal structure across three rift segments of the Gulf of California

NASA Astrophysics Data System (ADS)

Lizarralde, D. D.; Axen, G. J.; Brown, H. E.; Fletcher, J. M.; Fernandez, A. G.; Harding, A. J.; Holbrook, W. S.; Kent, G. M.; Paramo, P.; Sutherland, F. H.; Umhoefer, P. J.

2007-05-01

We present a summary of results from a crustal-scale seismic experiment conducted in the southern Gulf of California. This experiment, the PESCADOR experiment, imaged crustal structure across three rift segments, the Alarcon, Guaymas, and San José del Cabo to Puerto Vallarta (Cabo-PV) segments, using seismic refraction/wide-angle reflection data acquired with airgun sources and recorded by closely spaced (10-15 km) ocean-bottom seismometers (OBSs). The imaged crustal structure reveals a surprisingly large variation in rifting style and magmatism between these segments: the Alarcon segment is a wide rift with apparently little syn-rift magmatism; the Guaymas segment is a narrow, magmatically robust rift; and the Cabo-PV segment is a narrow, magmatically "normal" rift. Our explanation for the observed variability is non-traditional in that we do not invoke mantle temperature, the factor commonly invoked to explain end-member volcanic and non-volcanic rifted margins, as the source of the considerable, though non-end-member variability we observe. Instead, we invoke mantle depletion related to pre-rift arc volcanism to account for observed wide, magma-poor rifting and mantle fertility and possibly the influence of sediments to account for robust rift and post-rift magmatism. These factors may commonly vary over small lateral spatial scales in regions that have transitioned from convergent to extensional tectonics, as is the case for the Gulf of California and many other rifts. Our hypothesis suggests that substantial lateral variability may exist within the uppermost mantle beneath the Gulf of California today, and it is hoped that ongoing efforts to image upper mantle structure here will provide tests for this hypothesis.

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

USGS Publications Warehouse

Pollitz, Fred F.; Thatcher, Wayne R.

2010-01-01

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

Seismic imaging of the upper mantle beneath the northern Central Andean Plateau: Implications for surface topography

NASA Astrophysics Data System (ADS)

Ward, K. M.; Zandt, G.; Beck, S. L.; Wagner, L. S.

2015-12-01

Extending over 1,800 km along the active South American Cordilleran margin, the Central Andean Plateau (CAP) as defined by the 3 km elevation contour is second only to the Tibetan Plateau in geographic extent. The uplift history of the 4 km high Plateau remains uncertain with paleoelevation studies along the CAP suggesting a complex, non-uniform uplift history. As part of the Central Andean Uplift and the Geodynamics of High Topography (CAUGHT) project, we use surface waves measured from ambient noise and two-plane wave tomography to image the S-wave velocity structure of the crust and upper mantle to investigate the upper mantle component of plateau uplift. We observe three main features in our S-wave velocity model including (1), a high velocity slab (2), a low velocity anomaly above the slab where the slab changes dip from near horizontal to a normal dip, and (3), a high-velocity feature in the mantle above the slab that extends along the length of the Altiplano from the base of the Moho to a depth of ~120 km with the highest velocities observed under Lake Titicaca. A strong spatial correlation exists between the lateral extent of this high-velocity feature beneath the Altiplano and the lower elevations of the Altiplano basin suggesting a potential relationship. Non-uniqueness in our seismic models preclude uniquely constraining this feature as an uppermost mantle feature bellow the Moho or as a connected eastward dipping feature extending up to 300 km in the mantle as seen in deeper mantle tomography studies. Determining if the high velocity feature represents a small lithospheric root or a delaminating lithospheric root extending ~300 km into the mantle requires more integration of observations, but either interpretation shows a strong geodynamic connection with the uppermost mantle and the current topography of the northern CAP.

The North American upper mantle: density, composition, and evolution

USGS Publications Warehouse

Mooney, Walter D.; Kaban, Mikhail K.

2010-01-01

The upper mantle of North America has been well studied using various seismic methods. Here we investigate the density structure of the North American (NA) upper mantle based on the integrative use of the gravity field and seismic data. The basis of our study is the removal of the gravitational effect of the crust to determine the mantle gravity anomalies. The effect of the crust is removed in three steps by subtracting the gravitational contributions of (1) topography and bathymetry, (2) low-density sedimentary accumulations, and (3) the three-dimensional density structure of the crystalline crust as determined by seismic observations. Information regarding sedimentary accumulations, including thickness and density, are taken from published maps and summaries of borehole measurements of densities; the seismic structure of the crust is based on a recent compilation, with layer densities estimated from P-wave velocities. The resultant mantle gravity anomaly map shows a pronounced negative anomaly (−50 to −400 mGal) beneath western North America and the adjacent oceanic region and positive anomalies (+50 to +350 mGal) east of the NA Cordillera. This pattern reflects the well-known division of North America into the stable eastern region and the tectonically active western region. The close correlation of large-scale features of the mantle anomaly map with those of the topographic map indicates that a significant amount of the topographic uplift in western NA is due to buoyancy in the hot upper mantle, a conclusion supported by previous investigations. To separate the contributions of mantle temperature anomalies from mantle compositional anomalies, we apply an additional correction to the mantle anomaly map for the thermal structure of the uppermost mantle. The thermal model is based on the conversion of seismic shear-wave velocities to temperature and is consistent with mantle temperatures that are independently estimated from heat flow and heat production data. The thermally corrected mantle density map reveals density anomalies that are chiefly due to compositional variations. These compositional density anomalies cause gravitational anomalies that reach ~250 mGal. A pronounced negative anomaly (−50 to −200 mGal) is found over the Canadian shield, which is consistent with chemical depletion and a corresponding low density of the lithospheric mantle, also referred to as the mantle tectosphere. The strongest positive anomaly is coincident with the Gulf of Mexico and indicates a positive density anomaly in the upper mantle, possibly an eclogite layer that has caused subsidence in the Gulf. Two linear positive anomalies are also seen south of 40°N: one with a NE-SW trend in the eastern United States, roughly coincident with the Grenville-Appalachians, and a second with a NW-SE trend beneath the states of Texas, New Mexico, and Colorado. These anomalies are interpreted as being due to (1) the presence of remnants of an oceanic slab in the upper mantle beneath the Grenville-Appalachian suture and (2) mantle thickening caused by a period of shallow, flat subduction during the Laramie orogeny, respectively. Based on these geophysical results, the evolution of the NA upper mantle is depicted in a series of maps and cartoons that display the primary processes that have formed and modified the NA crust and lithospheric upper mantle.

Uppermost mantle structure beneath eastern China and its surroundings from Pn and Sn tomography

NASA Astrophysics Data System (ADS)

Sun, Weijia; Kennett, B. L. N.

2016-04-01

The Pn and Sn residuals from regional events provide strong constraints on the structure and lithological characteristics of the uppermost mantle beneath eastern China and its surroundings. With the dense Chinese Digital Seismic Network in eastern China, separate Pn and Sn tomographic inversions have been exploited to obtain P and S velocities at a resolution of 2° × 2° or better. The patterns of P velocities are quite consistent with the S velocities at depth of 50 and 60 km, but the amplitude of P wave speed anomalies are a little larger than those of S wave speed. The low P wave speed, high S wave speed, and low Vp/Vs ratio beneath the northern part of Ordos Basin are related to upwelling hot material. Abrupt changes in material properties are indicated from the rapid variations in the Vp/Vs ratio.

Evidence of phase nucleation during olivine diffusion creep: a new perspective for mantle strain localisation

NASA Astrophysics Data System (ADS)

Précigout, Jacques; Stünitz, Holger

2017-04-01

Mantle strain localisation is of great importance for lithosphere dynamics, but the cause for this phenomenon remains very elusive, particularly in conditions of the strong and ductile uppermost mantle. In these latter, grain size reduction leading to diffusion creep in olivine is believed to be one of the best candidates to account for strain localisation. However, the mechanisms of grain size reduction in this regime are still poorly understood. Here we show the results of Griggs-type experiments that document grain size reduction and material weakening during wet olivine diffusion creep at 900 °C and 1.2 GPa. While occurring for both, mono-phase and two-phase aggregates, grain size reduction is coeval with strain localisation and local phase mixing in olivine-pyroxene aggregates. Based on evidence of fluid inclusions and cracks filled with a fine-grained phase mixture, we conclude that grain size reduces as a result of fluid-assisted nucleation. Cavitation induced by grain boundary sliding (creep cavitation) can be inferred, and may play a critical role for olivine grain size reduction. Amongst their implications for rock rheology in general, our findings highlight a key process for strain localisation in the ductile uppermost mantle. This study has been published under the reference: "Précigout, J., and Stünitz, H. (2016) Evidence of phase nucleation during olivine diffusion creep: a new perspective for mantle strain localisation. Earth and Planetary Science Letters 455: 94-105, doi:101016/j.epsl.2016.09.029".

Evidence of dehydration in peridotites from Eifel Volcanic Field and estimates of the rate of magma ascent

NASA Astrophysics Data System (ADS)

Denis, Carole M. M.; Demouchy, Sylvie; Shaw, Cliff S. J.

2013-05-01

We report major element compositions and water contents in upper mantle minerals from peridotites transported by silica undersaturated, mafic alkaline lavas from three volcanoes Rockeskyllerkopf, Dreiser Weiher, and Meerfelder Maar of the Eifel Volcanic Field (West Germany). The hydrogen concentrations (expressed in ppm wt. H2O) obtained from unpolarized and polarized Fourier transform infrared (FTIR) spectroscopy give water contents for olivine, enstatite and diopside of ~ 6 ppm wt. H2O, ~ 200 ppm wt. H2O and ~ 285 ppm wt. H2O, respectively. The hydrogen concentration in individual olivine grains is strongly heterogeneous whereas that in pyroxenes is homogeneous. Profiles measured across crystallographically oriented single-crystals of olivine using polarized infrared radiation reveal hydrogen depleted rims that are interpreted to be due to partial dehydration by ionic diffusion during the ascent of the xenolith to the surface. Using experimentally obtained diffusion coefficients for hydrogen in olivine at high temperature and high pressure, we estimate that the duration of the dehydration for the spinel-bearing xenoliths is limited to a few hours yielding rates of magma ascent from 3 ms- 1 to 12 ms- 1. Our study suggests that the water contents of the upper mantle based solely on measurements of mantle-derived olivine, when concentration is not homogeneous, underestimate the true water content of the equilibrated uppermost mantle and that pyroxenes are a better proxy to constrain uppermost mantle water contents.

Role of upper-most crustal composition in the evolution of the Precambrian ocean-atmosphere system

NASA Astrophysics Data System (ADS)

Large, R. R.; Mukherjee, I.; Zhukova, I.; Corkrey, R.; Stepanov, A.; Danyushevsky, L. V.

2018-04-01

Recent research has emphasized the potential relationships between supercontinent cycles, mountain building, nutrient flux, ocean-atmosphere chemistry and the origin of life. The composition of the Upper-Most Continental Crust (UMCC) also figures prominently in these relationships, and yet little detailed data on each component of this complex relationship has been available for assessment. Here we provide a new set of data on the trace element concentrations, including the Rare Earth Elements (REE), in the matrix of 52 marine black shale formations spread globally through the Archean and Proterozoic. The data support previous studies on the temporal geochemistry of shales, but with some important differences. Results indicate a change in provenance of the black shales (upper-most crustal composition), from more mafic in the Archean prior to 2700 Ma, to more felsic from 2700 to 2200 Ma, followed by a return to mafic compositions from 2200 to 1850 Ma. Around 1850 to 1800 Ma there is a rapid change to uniform felsic compositions, which remained for a billion years to 800 Ma. The shale matrix geochemistry supports the assertion that the average upper-most continental source rocks for the shales changed from a mix of felsic, mafic and ultramafic prior to 2700 Ma to more felsic after 1850 Ma, with an extended transition period between. The return to more mafic UMCC from 2200 to 1850 Ma is supported by the frequency of Large Igneous Provinces (LIPs) and banded iron formations, which suggest a peak in major mantle-connected plume events and associated Fe-rich hydrothermal activity over this period. Support for the change to felsic UMCC around 1850 Ma is provided by previous geological data which shows that felsic magmas, including, A-type granites and K-Th-U-rich granites intruded vast areas of the continental crust, peaking around 1850 Ma and declining to 1000 Ma. The implications of this change in UMCC are far reaching and may go some way to explain the distinct features of the Boring Billion (1800-800 Ma). Firstly, because mafic-ultramafic rocks contain significantly higher levels of the bio-essential nutrient elements (e.g. Fe, P, Ni, Cr, Co, Cu, Se, Mn, Zn) compared with felsic rocks, the flux of macro- and micro-nutrients to the ocean would have decreased significantly post 1850 Ma. This would have contributed to a drop in productivity and a drop in atmosphere O2 as suggested by the marine pyrite proxy. In addition, a change from mafic to felsic dominant composition of the UMCC post 1850 Ma, would have led to a decrease in the erosive flux of Ca and Mg to the ocean, affecting the oceanic carbonate equilibrium and likely contributing to a rise in atmosphere CO2. On this basis, we speculate that the commencement of the middle Proterozoic, commonly known as the Boring Billion period from 1800 to 800 Ma, marks the start of an extended time in Earth's evolution when the UMCC became dominated by felsic rocks, particularly K-U-Th-anorogenic granites. This led to a period of anomalously low concentrations of bio-essential trace elements, but elevated REE, U, Th, Pb, Tl, Rb/Al and K/Na in the oceans.

Gravity field over northern Eurasia and variations in the strength of the upper mantle

NASA Technical Reports Server (NTRS)

Kogan, Mikhail G.; Mcnutt, Marcia K.

1993-01-01

The correlation of long-wavelength gravity anomalies in northern Eurasia with seismic velocity anomalies in the upper mantle reverses in sign between western and eastern Eurasia. The difference between western and eastern Eurasia can be explained by the presence of a low-viscosity zone in the uppermost mantle beneath eastern Eurasia that is absent to the west. The location of the lateral change in viscosity corresponds with the geologic boundary between the older shields and platforms of the Baltics, Russia, and Siberia and the younger, geologically active mountain belts of eastern Asia. This relation provides evidence that differences in the strength of the upper mantle control the locus of intracontinental deformation.

Lithospheric evolution of the Northern Arabian Shield: Chemical and isotopic evidence from basalts, xenoliths and granites

NASA Technical Reports Server (NTRS)

Stein, M.

1988-01-01

The evolution of the upper-mantle and the lower-crust (the conteinental lithosphere), is the area of Israel and Sinai was studied, using the chemical composition and the Nd-Sr isotopic systematics from mantle and crustal nodules, their host basalts, and granites. The magmatism and the metasomatism making the lithosphere are related to uprise of mantle diapirs in the uppermost mantle of the area. These diapirs heated the base of the lithosphere, eroded, and replaced it with new hot material. It caused a domal uplift of the lithosphere (and the crust). The doming resulted in tensional stresses that in turn might develop transport channels for the basalt.

Surface wave tomography applied to the North American upper mantle

NASA Astrophysics Data System (ADS)

van der Lee, Suzan; Frederiksen, Andrew

Tomographic techniques that invert seismic surface waves for 3-D Earth structure differ in their definitions of data and the forward problem as well as in the parameterization of the tomographic model. However, all such techniques have in common that the tomographic inverse problem involves solving a large and mixed-determined set of linear equations. Consequently these inverse problems have multiple solutions and inherently undefinable accuracy. Smoother and rougher tomographic models are found with rougher (confined to great circle path) and smoother (finite-width) sensitivity kernels, respectively. A powerful, well-tested method of surface wave tomography (Partitioned Waveform Inversion) is based on inverting the waveforms of wave trains comprising regional S and surface waves from at least hundreds of seismograms for 3-D variations in S wave velocity. We apply this method to nearly 1400 seismograms recorded by digital broadband seismic stations in North America. The new 3-D S-velocity model, NA04, is consistent with previous findings that are based on separate, overlapping data sets. The merging of US and Canadian data sets, adding Canadian recordings of Mexican earthquakes, and combining fundamental-mode with higher-mode waveforms provides superior resolution, in particular in the US-Canada border region and the deep upper mantle. NA04 shows that 1) the Atlantic upper mantle is seismically faster than the Pacific upper mantle, 2) the uppermost mantle beneath Precambrian North America could be one and a half times as rigid as the upper mantle beneath Meso- and Cenozoic North America, with the upper mantle beneath Paleozoic North America being intermediate in seismic rigidity, 3) upper-mantle structure varies laterally within these geologic-age domains, and 4) the distribution of high-velocity anomalies in the deep upper mantle aligns with lower mantle images of the subducted Farallon and Kula plates and indicate that trailing fragments of these subducted oceanic plates still reside in the transition zone. The thickness of the high-velocity layer beneath Precambrian North America is estimated to be 250±70 km thick. On a smaller scale NA04 shows 1) high-velocities associated with subduction of the Pacific plate beneath the Aleutian arc, 2) the absence of expected high velocities in the upper mantle beneath the Wyoming craton, 3) a V-shaped dent below 150 km in the high-velocity cratonic lithosphere beneath New England, 4) the cratonic lithosphere beneath Precambrian North America being confined southwest of Baffin Bay, west of the Appalachians, north of the Ouachitas, east of the Rocky Mountains, and south of the Arctic Ocean, 5) the cratonic lithosphere beneath the Canadian shield having higher S-velocities than that beneath Precambrian basement that is covered with Phanerozoic sediments, 6) the lowest S velocities are concentrated beneath the Gulf of California, northern Mexico, and the Basin and Range Province.

Two lithospheric profiles across southern California derived from gravity and seismic data

USGS Publications Warehouse

Romanyuk, T.; Mooney, W.D.; Detweiler, S.

2007-01-01

We present two detailed 2-D density transects for the crust and uppermost mantle across southern California using a linear gravity inversion technique. This technique parameterizes the crust and upper mantle as a set of blocks that are based on published geologic and seismic models. Each block can have a range of densities that are constrained where possible by borehole measurements, seismic velocities, and petrologic data. To further constrain the models, it is assumed that the lithosphere is close to isostatic equilibrium at both ends of the profiles, in the deep ocean and east of the Mojave Desert. We calculate the lithostatic pressure variations field for the whole cross section to rule out the geophysically insignificant solutions. In the linear equation, ?? = a + bV (V, seismic P-wave velocity; ??, density), which approximates the mantle density-velocity (??-V) relationship, different coefficients for b were evaluated. Lower coefficients (b 0.3) imply that other effects, such as composition and/or metamorphic changes, play an important role in the mantle. Density models were constructed with the coefficient b ranging from 0 to 0.6. The results indicate that a high b value in the mantle ??-V relationship is associated with less dense crust in the Mojave block and more dense crust in the Catalina schist block. In the less dense Mojave block, the average density of the whole crust is ???2.75 g/cm3, while that of the lower crust is ???2.72 g/cm3. These densities imply a high silica content in the crust, and a minor fraction of basic rock in the lower crust, or perhaps the absence of a basaltic layer altogether. By comparison, the average density of a typical continental stable platform is ???2.85 g/cm3. Models with higher b coefficients (0.5-0.6) are characterized by a large isostatic imbalance. On the other hand, lower b values (0-0.2) require a consolidated whole crust density in the Mojave Desert of ???2.78 g/cm3, and a lower crust density of ???2.89 g/cm3 with mostly basaltic composition. This contradicts the observed, lower Vp/Vs-ratio in the Mojave Desert associated with mostly felsic and low-density crust. Models with lower b coefficients (0.1-0.2) are characterized by an absence of local Airy compensation beneath the San Gabriel Mountains at the LARSE-1 profile. These, and other non-gravity arguments, suggest optimal solutions to the mantle ??-V relation of b ??? 0.2-0.4. This, in turn, means that both thermal and petrological effects occur inside the downwelling of the uppermost mantle high velocity body located beneath the Transverse Ranges. During the development of this mantle downwelling, the basaltic layer of the Mojave block was likely eroded and pulled down into the high velocity body. Those basaltic fragments may have been transformed into eclogites, and this metamorphic change implies a higher b-coefficient density-velocity relationship than would be expected for a purely thermal process.

Radial Anisotropy in the Mantle Transition Zone and Its Implications

NASA Astrophysics Data System (ADS)

Chang, S. J.; Ferreira, A. M.

2016-12-01

Seismic anisotropy is a useful tool to investigate mantle flow, mantle convection, and the presence of melts in mantle, since it provides information on the direction of mantle flow or the orientation of melts by combining it with laboratory results in mineral physics. Although the uppermost and lowermost mantle with strong anisotropy have been well studied, anisotropic properties of the mantle transition zone is still enigmatic. We use a recent global radially anisotropic model, SGLOBE-rani, to examine the patterns of radial anisotropy in the mantle transition zone. Strong faster SV velocity anomalies are found in the upper transition zone beneath subduction zones in the western Pacific, which decrease with depth, thereby nearly isotropic in the lower transition zone. This may imply that the origin for the anisotropy is the lattice-preferred orientation of wadsleyite, the dominant anisotropic mineral in the upper transition zone. The water content in the upper transition zone may be inferred from radial anisotropy because of the report that anisotropic intensity depends on the water content in wadsleyite.

How Much Ocean Is Left Between Libya and Crete

NASA Astrophysics Data System (ADS)

Makris, J.; Yegorova, T.

The intense deformation of the Hellenides is due to crustal shortening and the collision between the European and African Plates. This processes creates the Mediterranean accretionary wedge known as Mediterranean Ridge, which is composed of thick sedi- mentary sequences exceeding 10 km in thickness. The stage of this collision has been under dispute for many years. We performed wide aperture seismic soundings between Crete and Libya along 5 seismic lines. The results were used to constrain gravity mod- elling and develop density models in 2D and 3D between Libya and the Cretan Sea. We identified the limits of the European continental crust extending south of Crete for more than 100 km and building the backstop of the sediment accumulation . The African continental crust extends to the north for about 80 to 100 km, so that the remaining space floored by the oceanic Thethian basement is at its narrowest point not more than 100 to 120 km wide. By modelling in 3D the gravity field of the sedi- ments, crust and uppermost mantle we identified significant variations of the density distribution of the upper mantle. The young intensely deforming area of the Aegean domain is floored by low density upper mantle due to the mobilization of magma and the activation of the thermal regime. The subducted cold oceanic slab sinks below the Cretan crust in NE orientation and is decupled from the continental crust between central Crete and the southeastern edge of the Peloponnese. The deformation of the sediments controlled by the compressional processes have their maximum accumu- lation at the limits of the backstop. Here the transition of the deep trough to the flat and nearly undeformed sedimentary sequence is very abrupt and the transition oc- curs along vertical displacements of 6 to 8 km near vertical throw. Near the southern transition of the oceanic crust to the African continental domain obducted ophiolites extend over large areas explaining gravity highs and also observed intense magnetic anomalies. The computed gravity field fits the observed one in all its low frequency spectrum. We avoided modelling the high frequency part of the field since the seismic 1 information was not dense enough to justify the effort. 2

Dependence of the brittle ductile transition on strain-rate-dependent critical homologous temperature

NASA Astrophysics Data System (ADS)

Davis, Paul M.

2017-05-01

Earthquakes mainly occur in crust or mantle that is below a critical temperature for the tectonic strain-rate, \\dot{e}_t, such that stress builds up to the breaking point before it can relax due to creep. Then long-range stress correlation gives rise to power law seismicity including large events. The limiting temperature depends on pressure, which is taken into account by finding a critical homologous temperature THc = T/TM above which earthquakes are rarely observed (where T, TM are temperature and average melting temperature of constituent minerals). We find that THc for ocean plates is ∼0.55. For California earthquakes, it is also close to 0.55. The uppermost mantle layer of oceanic plates of thickness ∼50 km is composed of harzburgite and depleted peridotite from which basalt has been removed to form ocean crust. Thus it has a higher melting temperature than the peridotite of the surrounding mantle, or the lower halves of plates. Thicknesses of seismicity in deep subduction zones, determined from 2-D polynomial fits to a relocated catalogue, are ∼50 km, which suggests that the earthquake channel is confined to this layer. We construct models to find homologous temperatures in slabs, and find that seismicity thicknesses are also, on average, confined to TH ≤ 0.55 ± 0.05. The associated rheology is compared with that obtained from flexure models of ocean lithosphere. The brittle-ductile transition occurs where viscosity drops from high values in the cold cores of slabs to values of 1022-1023 Pa s, that is, where creep strain-rates become comparable to tectonic rates. The cut-off for deep earthquakes is not sharp. However they appear unlikely to occur if homologous temperature is high TH > 0.55. Exceptions to the rule are anomalously deep earthquakes such as those beneath the Iceland and the Hawaiian hotspots, and the Newport Inglewood Fault. These are smaller events with short-range stress correlation, and can be explained if strain-rates are two to three orders of magnitude higher than those associated with earthquakes located where TH ≤ 0.55. We conclude that the brittle-ductile transition corresponds to the transition from long-range (regional) to short-range (localized on asperities) stress correlation.

Ambient Noise Tomography and Microseism Directionalities across the Juan de Fuca Plate

NASA Astrophysics Data System (ADS)

Tian, Ye

Ambient noise tomography has been well developed over the past decade and proven to be effective in studying the crust and upper mantle structure beneath the Earth’s continents. With new seismic array deployments beginning in the oceans, the application of the tomographic methods based on ambient noise observed at ocean bottom seismometers (OBSs) has become an important topic for research. In this thesis, I investigate the application of ambient noise tomography to oceanic bottom seismic data recorded by the Cascadia Initiative experiment across the Juan de Fuca plate. With higher local noise levels recorded by OBSs, I find that traditional data processing procedures used in ambient noise tomography produce measurable Rayleigh wave Green’s functions between deep ocean stations, whereas the shallow water stations are severely contaminated by both tilt noise and compliance noise and require new methods of processing. Because the local noise level varies across the study region, four semi-independent studies are conducted to both utilize the quieter deep-water stations and to address the problem posed by noisy shallow water stations. First, I construct an age-dependent shear wave speed model of the crust and uppermost mantle with 18 deep-water stations near the Juan de Fuca Ridge. The model possess a shallow low shear velocity zone near the ridge and has its sedimentary thickness, lithospheric thickness, and mantle shear wave speeds increase systematically with age Second, I investigate the locations and mechanisms of microseism generation using ambient noise cross-correlations constructed between 61 OBSs and 42 continental stations near the western US coast and find that the primary and secondary microseisms are generated at different locations and possibly have different physical mechanisms. Third, I show that tilt and compliance noise on the vertical components of the OBSs can be reduced substantially using the horizontal components and the differential pressure gauge records. Removal of these types of noise improves the signal-to-noise ratio of ambient noise cross-correlations significantly at beyond 10 sec period. Lastly, I present a new single-station method to estimate the microseism Rayleigh wave strength and directionality based on the horizontal-to-vertical transfer function. The high spatial and temporal resolution of this method may open up the microseism Rayleigh waves for a wider range of studies.

Difference of the seismic structure between the Hyuga-nada and the Nankai seismogenic segments

NASA Astrophysics Data System (ADS)

Yamamoto, Y.; Obana, K.; Takahashi, T.; Nakanishi, A.; Kodaira, S.; Kaneda, Y.

2010-12-01

In the Nankai Trough, three major seismogenic zones of megathrust earthquake exist (Tokai, Tonankai and Nankai earthquake regions). The Hyuga-nada region was distinguished from these seismogenic zones because of the lack of megathrust earthquake. In the Hyuga-nada region, interplate earthquakes of M~7 occur repeatedly at intervals of about 20 years whereas no megathrust (M > 8) earthquakes had been recognized up to now. However, recent studies show the possibility of simultaneous rupture of the Tokai, Tonankai, Nankai and Hyuga-nada segments was also pointed out [e.g., Hori et al., 2009 AOGS]. To understand the possibility of seismic linkage of Nankai and Hyuga-nada segments, Japan Agency for Marine-Earth Science and Technology has been carried out a wide-angle active source survey and local seismic observation in the western end of the Nankai seismogenic zone, as a part of Research concerning Interaction Between the Tokai, Tonankai and Nankai Earthquakes' funded by Ministry of Education, Culture, Sports, Science and Technology, Japan. Nakanishi et al [2009, AGU] showed that subducting Philippine Sea Plate can be divided into three zones and there is the zone of the thin oceanic crust of the subducting Philippine Sea Plate between Nankai segment and Kyushu-Palau Ridge segment by analyzing of the active source survey. Deep structure of the subducting slab is also important to consider the possibility of the seismic linkage and the location of the boundary among three zones described above. To obtain the deep seismic image, we performed a three-dimensional seismic tomography using the local seismic data recorded on 158 ocean bottom seismographs and 105 land seismic stations. From these data, we could detect 1141 earthquakes in the Hyuga-nada region. From the result of hypocenter relocation, microseismicity near the trough axis is active on the western part of the ‘thin oceanic crust’, whereas inactive on the eastern part. Besides, velocity structure of the uppermost part of the subducting slab mantle shows spatial heterogeneities. In the thin oceanic crust zone, high velocity slab mantle is imaged from near the trough to coastline. On the other hands, there is low velocity zone in the slab mantle near the trough axis in the Kyusyu-Palau Ridge segment. This low velocity zone may be related to the location of the eastern end of subducted Kyusyu-Palau Ridge.

Shear Wave Structure in the Lithosphere of Texas from Ambient Noise Tomography

NASA Astrophysics Data System (ADS)

Yao, Y.; Li, A.

2014-12-01

Texas contains several distinct tectonic provinces, the Laurentia craton, the Ouachita belt, and the Gulf coastal plain. Although numerous geophysical experiments have been conducted in Texas for petroleum exploration, the lithosphere structure of Texas has not been well studied. We present here the Texas-wide shear wave structure using seismic ambient noise data recorded at 87 stations from the Transportable Array of the USArray between March 2010 and February 2011. Rayleigh wave phase velocities between pairs of stations are obtained by cross-correlating long ambient noise sequences and are used to develop phase velocity maps from 6 to 40 s. These measured phase velocities are used to construct 1-D and 3-D shear wave velocity models, which consist of four crust layers and one upper mantle layer. Shear wave velocity maps reveal a close correlation with major geological features. From the surface to 25 km depth, Positive anomalies coincide with the Laurentia craton, and negative anomalies coincide with the continental margin. The boundary of positive-negative anomaly perfectly matches the Ouachita belt. The Llano Uplift is imaged as the highest velocity through the mid-crust because the igneous rock forming the uplift has faster seismic velocity than the normal continental crust. Similarly, three small high-velocity areas exist beneath the Waco Uplift, Devils River Uplift, and Benton Uplift, even though surface geological traces are absent in these areas. The lowest velocity at the shallow crust appears in northeastern and southeastern Texas separated by the San Marcos Arch, correlating with thick sediment layers. An exceptional low velocity is imaged in southernmost Texas in the lower crust and upper mantle, probably caused by subducted wet oceanic crust before the rifting in the Gulf of Mexico. In the uppermost mantle, positive shear wave anomalies extend southeastward from the Ouachita belt to the Gulf coast, likely evidencing the subducted oceanic lithosphere during the Ouachita orogeny. This observation need be further tested using long period surface wave dispersions from earthquakes, which help to improve model resolution in the upper mantle.

Deep mantle cycling of oceanic crust: evidence from diamonds and their mineral inclusions.

PubMed

Walter, M J; Kohn, S C; Araujo, D; Bulanova, G P; Smith, C B; Gaillou, E; Wang, J; Steele, A; Shirey, S B

2011-10-07

A primary consequence of plate tectonics is that basaltic oceanic crust subducts with lithospheric slabs into the mantle. Seismological studies extend this process to the lower mantle, and geochemical observations indicate return of oceanic crust to the upper mantle in plumes. There has been no direct petrologic evidence, however, of the return of subducted oceanic crustal components from the lower mantle. We analyzed superdeep diamonds from Juina-5 kimberlite, Brazil, which host inclusions with compositions comprising the entire phase assemblage expected to crystallize from basalt under lower-mantle conditions. The inclusion mineralogies require exhumation from the lower to upper mantle. Because the diamond hosts have carbon isotope signatures consistent with surface-derived carbon, we conclude that the deep carbon cycle extends into the lower mantle.

How did the Lunar Magma Ocean crystallize?

NASA Astrophysics Data System (ADS)

Davenport, J.; Neal, C. R.

2012-12-01

It is generally accepted that the lunar crust and at least the uppermost (500 km) mantle was formed by crystallization of a magma ocean. How the magma ocean cooled and crystallized is still under debate. Parameters such as bulk composition, lunar magma ocean (LMO) crystallization method (fractional vs. equilibrium), depth of the LMO, and time for LMO solidification (effects of tidal heating mechanisms, insulating crustal lid, etc.) are still under debate. Neal (2001, JGR 106, 27865-27885) argues for the presence of garnet in the deep lunar mantle via compositional differences between low- and high-Ti mare basalts and volcanic glasses. Neal (2001) suggests that these compositional differences are due to the presence of garnet in the source regions of certain volcanic glass bead groups. As Neal (2001, JGR 106, 27865-27885) points out, determining if there is garnet in the lunar mantle is important in determining if the LMO was a "whole-Moon" event or if it was limited to certain areas. In the latter case, garnet would have been preserved in the lunar mantle and would have been used in the source material for some of the volcanic glasses. High-pressure experimental work concludes that with the right T-P conditions (2.5-4.5 GPa and 1675-1800° C) there could be a garnet-bearing pyroxene rich protolith at ~500 km depth. This also has significant implications for the bulk Al2O3 composition of the initial bulk Moon. If the LMO was not global, the volcanic glass beads that show evidence of garnet in their sources were formed from the deep, primitive lunar mantle, it begs the questions how was the non-LMO regions of the Moon formed and what was it's bulk composition? To try to answer these questions, it is necessary to thoroughly model the evolution of the LMO and then use that work to model the sources and formation of mare basalts, the volcanic glass beads, and other regions in question. To begin to answer these questions, we developed a scenario we have termed reverse LMO modeling. Geochemical compositions such as KREEP, ur-KREEP and FAN will be run backwards through various LMO models that have been proposed in the literature. The concentration of the initial bulk Moon, according to the concentrations of the particular type of rock being used, can be modeled by taking this from 0 percent liquid (PCL; a completely solidified Moon) to 100 PCL. Using the KREEP composition reported by Warren and Wasson (1979, Rev. of Geophysics and Space Physics 17, 73-88), Warren (1988, Proc. 18th LPSC, 233-241) and Warren (1989, LPI Tech. Report 89, 149-153), the Mg numbers (Mg#) for the bulk initial Moon were calculated yielding 0.87, 0.76, and 0.86 respectively. The major element compositions of calculated bulk Moon compositions have elevated Al2O3, FeO, and TiO2, consistent with the presence of garnet in the lunar mantle as well as generating high-Ti basalts. Using these data we can model the petrogenesis of the low- and high-Ti mare basalt and volcanic glass source regions. Furthermore, using remote sensing and the calculated source data we can compare the modeled concentrations of these rocks to where these ranges of concentrations fall on the Moon's surfaces, so that we can constrain the areas where the presence of a magma ocean on the Moon was possible.

Understanding the nature of mantle upwelling beneath East-Africa

NASA Astrophysics Data System (ADS)

Civiero, Chiara; Hammond, James; Goes, Saskia; Ahmed, Abdulhakim; Ayele, Atalay; Doubre, Cecile; Goitom, Berhe; Keir, Derek; Kendall, Mike; Leroy, Sylvie; Ogubazghi, Ghebrebrhan; Rumpker, Georg; Stuart, Graham

2014-05-01

The concept of hot upwelling material - otherwise known as mantle plumes - has long been accepted as a possible mechanism to explain hotspots occurring at Earth's surface and it is recognized as a way of removing heat from the deep Earth. Nevertheless, this theory remains controversial since no one has definitively imaged a plume and over the last decades several other potential mechanisms that do not require a deep mantle source have been invoked to explain this phenomenon, for example small-scale convection at rifted margins, meteorite impacts or lithospheric delamination. One of the best locations to study the potential connection between hotspot volcanism at the surface and deep mantle plumes on land is the East African Rift (EAR). We image seismic velocity structure of the mantle below EAR with higher resolution than has been available to date by including seismic data recorded by stations from many regional networks ranging from Saudi Arabia to Tanzania. We use relative travel-time tomography to produce P- velocity models from the surface down into the lower mantle incorporating 9250 ray-paths in our model from 495 events and 402 stations. We add smaller earthquakes (4.5 < mb < 5.5) from poorly sampled regions in order to have a more uniform data coverage. The tomographic results allow us to image structures of ~ 100-km length scales to ~ 1000 km depth beneath the northern East-Africa rift (Ethiopia, Eritrea, Djibouti, Yemen) with good resolution also in the transition zone and uppermost lower mantle. Our observations provide evidence that the shallow mantle slow seismic velocities continue trough the transition zone and into the lower mantle. In particular, the relatively slow velocity anomaly beneath the Afar Depression extends up to depths of at least 1000 km depth while another low-velocity anomaly beneath the Main Ethiopian Rift seems to be present in the upper mantle only. These features in the lower mantle are isolated with a diameter of about 400 km indicating deep multiple sources of upwelling that converge in broader low-velocity bodies along the rift axis at shallow depths. Moreover, our preliminary models show that the low-velocity feature in the transition zone and uppermost lower mantle beneath Afar trends to the northeast beneath the Red Sea and Saudi Arabia as opposed to being linked to the African Superplume towards the southwest.

The uppermost mantle shear wave velocity structure of eastern Africa from Rayleigh wave tomography: constraints on rift evolution

NASA Astrophysics Data System (ADS)

O'Donnell, J. P.; Adams, A.; Nyblade, A. A.; Mulibo, G. D.; Tugume, F.

2013-08-01

An expanded model of the 3-D shear wave velocity structure of the uppermost mantle beneath eastern Africa has been developed using earthquakes recorded by the AfricaArray East African Seismic Experiment in conjunction with data from permanent stations and previously deployed temporary stations. The combined data set comprises 331 earthquakes recorded on a total of 95 seismic stations spanning Kenya, Uganda, Tanzania, Zambia and Malawi. In this study, data from 149 earthquakes were used to determine fundamental-mode Rayleigh wave phase velocities at periods ranging from 20 to 182 s using the two-plane wave method, and then combined with the similarly processed published measurements and inverted for a 3-D shear wave velocity model of the uppermost mantle. New features in the model include (1) a low-velocity region in western Zambia, (2) a high-velocity region in eastern Zambia, (3) a low-velocity region in eastern Tanzania and (4) low-velocity regions beneath the Lake Malawi rift. When considered in conjunction with mapped seismicity, these results support a secondary western rift branch striking southwestwards from Lake Tanganyika, likely exploiting the relatively weak lithosphere of the southern Kibaran Belt between the Bangweulu Block and the Congo Craton. We estimate a lithospheric thickness of ˜150-200 km for the substantial fast shear wave anomaly imaged in eastern Zambia, which may be a southward subsurface extension of the Bangweulu Block. The low-velocity region in eastern Tanzania suggests that the eastern rift branch trends southeastwards offshore eastern Tanzania coincident with the purported location of the northern margin of the proposed Ruvuma microplate. Pronounced velocity lows along the Lake Malawi rift are found beneath the northern and southern ends of the lake, but not beneath the central portion of the lake.

Crust and uppermost-mantle structure of Greenland and the Northwest Atlantic from Rayleigh wave group velocity tomography

NASA Astrophysics Data System (ADS)

Darbyshire, Fiona A.; Dahl-Jensen, Trine; Larsen, Tine B.; Voss, Peter H.; Joyal, Guillaume

2018-03-01

The Greenland landmass preserves ˜4 billion years of tectonic history, but much of the continent is inaccessible to geological study due to the extensive inland ice cap. We map out, for the first time, the 3-D crustal structure of Greenland and the NW Atlantic ocean, using Rayleigh wave anisotropic group velocity tomography, in the period range 10-80 s, from regional earthquakes and the ongoing GLATIS/GLISN seismograph networks. 1-D inversion gives a pseudo-3-D model of shear wave velocity structure to depths of ˜100 km with a horizontal resolution of ˜200 km. Crustal thickness across mainland Greenland ranges from ˜25 km to over 50 km, and the velocity structure shows considerable heterogeneity. The large sedimentary basins on the continental shelf are clearly visible as low velocities in the upper ˜5-15 km. Within the upper continental basement, velocities are systematically lower in northern Greenland than in the south, and exhibit a broadly NW-SE trend. The thinning of the crust at the continental margins is also clearly imaged. Upper-mantle velocities show a clear distinction between typical fast cratonic lithosphere (Vs ≥4.6 km s-1) beneath Greenland and its NE margin and anomalously slow oceanic mantle (Vs ˜4.3-4.4 km s-1) beneath the NW Atlantic. We do not observe any sign of pervasive lithospheric modification across Greenland in the regions associated with the presumed Iceland hotspot track, though the average crustal velocity in this region is higher than that of areas to the north and south. Crustal anisotropy beneath Greenland is strong and complex, likely reflecting numerous episodes of tectonic deformation. Beneath the North Atlantic and Baffin Bay, the dominant anisotropy directions are perpendicular to the active and extinct spreading centres. Anisotropy in the subcontinental lithosphere is weaker than that of the crust, but still significant, consistent with cratonic lithosphere worldwide.

A Preliminary Teleseismic Investigation of the Crust and Mantle Lithosphere Obtained from BISN in the Western Canadian Arctic

NASA Astrophysics Data System (ADS)

Schaeffer, A. J.; Snyder, D. B.; Cairns, S.; Elliot, B.; Audet, P.; Esteve, C.; Murray-Bergquist, L.; Falck, H.

2016-12-01

The tectonic evolution of the Beaufort Sea continental margin has contributed to the maturation of these rocks into a major petroleum reservoir. Recent shallow offshore seismic reflection data suggest that Banks Island represents thin crust along a rifted margin established during the opening of the Arctic Ocean. In this case, rifting of the margin caused Banks Island to subside and accumulate sediments rich in petroleum source material. The cooling history and further subsidence of these sediments is important for understanding the thermal maturation of petroleum products. Recently published surface-wave velocity models of North America indicate seismic velocities at 100-150 km depths similar to those beneath Canada's diamond mines in the central Slave craton north of Yellowknife. These results imply that Banks Island is part of the Canadian Shield and that any kimberlites found thereon might contain diamonds. However, the fast velocities are inconsistent with this being a tectonically disrupted and thinned lithosphere along the Arctic margin of the Canada Basin. The problem is therefore to reconcile mantle structure typical of the Canadian Shield with crust typical of a rifted passive margin. Possibly related seismicity beneath the Mackenize River Delta and offshore in the Beaufort Sea has been observed for decades but its origin remains unknown, although has been suggested as due to incipient subduction of oceanic lithosphere beneath the North American craton. Resolving these questions requires high-resolution 3-D seismic models obtained from an array of broadband seismograph stations. Here we present preliminary results on the structure of the crust and uppermost mantle underlying the western Canadian Arctic. These results are generated using new data from the Banks Island Seismograph Network (BISN), three stations installed over the summer of 2014 and 2015; augmented with several USArray Transportable Array stations and older POLARIS and CNSN stations on neighbouring Arctic Islands.

Scales of Heterogeneities in the Continental Crust and Upper Mantle

NASA Astrophysics Data System (ADS)

Tittgemeyer, M.; Wenzel, F.; Ryberg, T.; Fuchs, K.

1999-09-01

A seismological characterization of crust and upper mantle can refer to large-scale averages of seismic velocities or to fluctuations of elastic parameters. Large is understood here relative to the wavelength used to probe the earth.¶In this paper we try to characterize crust and upper mantle by the fluctuations in media properties rather than by their average velocities. As such it becomes evident that different scales of heterogeneities prevail in different layers of crust and mantle. Although we cannot provide final models and an explanation of why these different scales exist, we believe that scales of inhomogeneities carry significant information regarding the tectonic processes that have affected the lower crust, the lithospheric and the sublithospheric upper mantle.¶We focus on four different types of small-scale inhomogeneities (1) the characteristics of the lower crust, (2) velocity fluctuations in the uppermost mantle, (3) scattering in the lowermost lithosphere and on (4) heterogeneities in the mantle transition zone.

Implications of spinel compositions for the petrotectonic history of abyssal peridotite from Southwest Indian Ridge (SWIR)

NASA Astrophysics Data System (ADS)

Chen, T.; Jin, Z.; Wang, Y.; Tao, C.

2012-12-01

Abyssal peridotites generate at mid-ocean ridges. Lherzolite and harzburgite are the main rock types of peridotites in the uppermost mantle. The lherzolite subtype, less depleted and less common in ophiolites, characterizes mantle diapirs and slow-spreading ridges. Along the Earth's mid-ocean ridges, abyssal peridotites undergo hydration reactions to become serpentinite minerals, especially in slow to ultraslow spreading mid-ocean ridges. Spinel is common in small quantities in peridotites, and its compositions have often been used as petrogenetic indicators [1]. The Southwest Indian Ridge (SWIR) is one of the two ultraslow spreading ridges in the world. The studied serpentinized peridotite sample was collected by the 21st Voyage of the Chinese oceanic research ship Dayang Yihao (aka Ocean No. 1) from a hydrothermal field (63.5°E, 28.0°S, and 3660 m deep) in SWIR. The studied spinels in serpentinized lherzolite have four zones with different compositions: relic, unaltered core is magmatic Al-spinels; micro- to nano- sized ferrichromite zoned particles; narrow and discontinuous magnetite rim; and chlorite aureoles. The values Cr# of the primary Al-spinels indicate the range of melting for abyssal peridotites from SWIR extends from ~4% to ~7% [2]. The alteration rims of ferrichromite have a chemical composition characterized by Fe enrichment and Cr# increase indicating chromite altered under greenschist-amphibolite facies. Magnetites formed in syn- and post- serpentinization. Chlorite (clinochlore) formed at the boundary and crack of spinel indicating it had undergone with low-temperature MgO- and SiO2-rich hydrothermal fluids [3]. It suggests that serpentinized lherzolite from SWIR had undergone poly-stage hydration reactions with a wide range of temperature. Acknowledgments: EMPA experiment was carried out by Xihao Zhu and Shu Zheng in The Second Institute of Oceanography and China University of Geosciences, respectively. The work was supported by NSFC. References [1] Dick, H.J.B., et al., Contrib Mineral Petr., 86:54, 1984. [2] Hellebrand et al., Nature, 410: 677, 2001. [3] Hamdy, M.M. , et al., JGMR, 3(9): 232, 2011.

The role of solid-solid phase transitions in mantle convection

NASA Astrophysics Data System (ADS)

Faccenda, Manuele; Dal Zilio, Luca

2017-01-01

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

Constraints on the Amount of deeply subducted Water from numerical Models in comparison with natural Samples

NASA Astrophysics Data System (ADS)

Konrad-Schmolke, M.; Halama, R.

2014-12-01

The subduction of hydrated slab mantle to beyond-arc depths is the most important and yet weakly constrained factor in the quantification of the Earth's deep geologic water cycle. During subduction of hydrated oceanic lithosphere, dehydration reactions in the downgoing plate lead to a partitioning of water between upper and lower plate. Water retained in the slab is recycled into the mantle where it controls its rheology and thus plate tectonic velocities. Hence, quantification of the water partitioning in subduction zones is crucial for the understanding of mass transfer between the Earth's surface and the mantle. Combined thermomechanical and thermodynamic models yield quantitative constraints on the water cycle in subduction zones, but unless model results can be linked to natural observations, the reliability of such models remains speculative. We present combined thermomechanical, thermodynamic and geochemical models of active and paleo-subduction zones, whose results can be tested with independent geochemical features in natural rocks. In active subduction zones, evidence for the validity of our model comes from the agreement between modeled and observed across-arc trends of boron concentrations and isotopic compositions in arc volcanic rocks. In the Kamchatkan subduction zone, for example, the model successfully predicts complex geochemical patterns and the spatial distribution of arc volcanoes. In paleo-subduction zones (e.g. Western Gneiss Region and Western Alps), constraints on the water budget and dehydration behavior of the subducting slab come from trace element zoning patterns in ultra-high pressure (UHP) garnets. Distinct enrichments of Cr, Ni and REE in the UHP zones of the garnets can be reconciled by our models that predict intense rehydration and trace element re-enrichment of the eclogites at UHP conditions by fluids released from the underlying slab mantle. Models of present-day subduction zones indicate the presence of 2.5-6 wt.% of water within the uppermost 15 km of the subducted slab mantle. Depending on hydration depth, between 25 and 90% of this water is recycled into the deeper mantle. The Lower Devonian example from the Western Gneiss Region indicates that subduction of water into the Earth's deeper mantle is an active process at least since the middle Paleozoic.

Oceanic mantle rocks reveal evidence for an ancient, 1.2-1.3 Ga global melting event

NASA Astrophysics Data System (ADS)

Dijkstra, A. H.; Sergeev, D.; McTaminey, L.; Dale, C. W.; Meisel, T. C.

2011-12-01

It is now increasingly being recognized that many oceanic peridotites are refertilized harzburgites, and that the refertilization often masks an extremely refractory character of the original mantle rock 'protolith'. Oceanic peridotites are, when the effects of melt refertilization are undone, often too refractory to be simple mantle melting residues after the extraction of mid-ocean ridge basalts at a spreading center. Rhenium-osmium isotope analysis is a powerful method to look through the effects of refertilization and to obtain constraints on the age of the melting that produced the refractory mantle protolith. Rhenium-depletion model ages of such anomalously refractory oceanic mantle rocks - found as abyssal peridotites or as mantle xenoliths on ocean islands - are typically >1 Ga, i.e., much older than the ridge system at which they were emplaced. In my contribution I will show results from two case studies of refertilized anciently depleted mantle rocks (Macquarie Island 'abyssal' peridotites and Lanzarote mantle xenoliths). Interestingly, very refractory oceanic mantle rocks from sites all around the world show recurring evidence for a Mesoproterozoic (~1.2-1.3 Ga) melting event [1]. Therefore, oceanic mantle rocks seem to preserve evidence for ancient melting events of global significance. Alternatively, such mantle rocks may be samples of rafts of ancient continental lithospheric mantle. Laser-ablation osmium isotope 'dating' of large populations of individual osmium-bearing alloys from mantle rocks is the key to better constrain the nature and significance of these ancient depletion events. Osmium-bearing alloys form when mantle rocks are melted to high-degrees. We have now extracted over >250 detrital osmium alloys from placer gold occurrences in the river Rhine. These alloys are derived from outcrops of ophiolitic mantle rocks in the Alps, which include blocks of mantle rocks emplaced within the Tethys Ocean, and ultramafic lenses of unknown (Precambrian?) age in the pre-Alpine Massifs. Populations of model ages of these Rhine alloys show prominent peaks at 0.5 and 1.2-1.3 Ga. The 1.2-1.3 Mesoproterozoic age peak recorded by the Rhine Os alloy population does also occur in Os alloy age distributions of other ophiolites worldwide, generally as a subsidiary peak [2]. In summary, osmium isotope model ages from mantle rocks and mantle-derived individual osmium alloys worldwide collectively point to a Mesoproterozoic, 1.2-1.3 Ga high-degree mantle melting event of global significance. This event may be related to a slab-avalanche or whole-mantle overturn event in Mesoproterozoic times. [1] Dijkstra et al. (2010) J. Petrology 51, 469-493 [2] Pearson et al. (2007) Nature 449, 202-205

Mid-ocean ridges produced thicker crust in the Jurassic than in Recent times

NASA Astrophysics Data System (ADS)

Van Avendonk, H. J.; Harding, J.; Davis, J. K.; Lawver, L. A.

2016-12-01

We present a compilation of published marine seismic refraction data to show that oceanic crust was 1.7 km thicker on average in the mid-Jurassic (170 Ma) than along the present-day mid-ocean ridge system. Plate reconstructions in a fixed hotspot framework show that the thickness of oceanic crust does not correlate with proximity to mantle hotspots, so it is likely that mid-plate volcanism is not the cause of this global trend. We propose that more melt was extracted from the upper mantle beneath mid-ocean ridges in the Jurassic than in recent times. Numerical studies show that temperature increase of 1 degree C in the mantle can lead to approximately 50-70 m thicker crust, so the upper mantle may have cooled 15-20 degrees C/100 Myr since 170 Ma. This average temperature decrease is larger than the secular cooling rate of the Earth's mantle, which is roughly 10 degrees C/100 Myr since the Archean. Apparently, the present-day configuration and dynamics of continental and oceanic plates removes heat more efficiently from the Earth's mantle than in its earlier history. The increase of ocean crustal thickness with plate age is also stronger in the Indian and Atlantic oceans than in the Pacific Ocean basin. This confirms that thermal insulation by the supercontinent Pangaea raised the temperature of the underlying asthenospheric mantle, which in turn led to more magmatic output at the Jurassic mid-ocean ridges of the Indian and Atlantic oceans.

Osmium isotopes and mantle convection.

PubMed

Hauri, Erik H

2002-11-15

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

The lateral variation of P n velocity gradient under Eurasia

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yang, Xiaoning

We report that mantle lid P wave velocity gradient, or P n velocity gradient, reflects the depth and lateral variations of thermal and rheological state of the uppermost mantle. Mapping the P n velocity gradient and its lateral variation helps us gain insight into the temperature, composition, and dynamics of the uppermost mantle. In addition, because P n velocity gradient has profound influence on P n propagation behavior, an accurate mapping of P n velocity gradient also improves the modeling and prediction of P n travel times and amplitudes. In this study, I used measured P n travel times tomore » derive path-specific P n velocity gradients. I then inverted these velocity gradients for two-dimensional (2-D) P n velocity-gradient models for Eurasia based on the assumption that a path-specific Pn velocity gradient is the mean of laterally varying P n velocity gradients along the P n path. Result from a Monte Carlo simulation indicates that the assumption is appropriate. The 2-D velocity-gradient models show that most of Eurasia has positive velocity gradients. High velocity gradients exist mainly in tectonically active regions. Most tectonically stable regions show low and more uniform velocity gradients. In conclusion, strong velocity-gradient variations occur largely along convergent plate boundaries, particularly under overriding plates.« less

The lateral variation of P n velocity gradient under Eurasia

DOE PAGES

Yang, Xiaoning

2017-05-03

We report that mantle lid P wave velocity gradient, or P n velocity gradient, reflects the depth and lateral variations of thermal and rheological state of the uppermost mantle. Mapping the P n velocity gradient and its lateral variation helps us gain insight into the temperature, composition, and dynamics of the uppermost mantle. In addition, because P n velocity gradient has profound influence on P n propagation behavior, an accurate mapping of P n velocity gradient also improves the modeling and prediction of P n travel times and amplitudes. In this study, I used measured P n travel times tomore » derive path-specific P n velocity gradients. I then inverted these velocity gradients for two-dimensional (2-D) P n velocity-gradient models for Eurasia based on the assumption that a path-specific Pn velocity gradient is the mean of laterally varying P n velocity gradients along the P n path. Result from a Monte Carlo simulation indicates that the assumption is appropriate. The 2-D velocity-gradient models show that most of Eurasia has positive velocity gradients. High velocity gradients exist mainly in tectonically active regions. Most tectonically stable regions show low and more uniform velocity gradients. In conclusion, strong velocity-gradient variations occur largely along convergent plate boundaries, particularly under overriding plates.« less

Deep structure of the Alborz Mountains by joint inversion of P receiver functions and dispersion curves

NASA Astrophysics Data System (ADS)

Rastgoo, Mehdi; Rahimi, Habib; Motaghi, Khalil; Shabanian, Esmaeil; Romanelli, Fabio; Panza, Giuliano F.

2018-04-01

The Alborz Mountains represent a tectonically and seismically active convergent boundary in the Arabia - Eurasia collision zone, in western Asia. The orogenic belt has undergone a long-lasted tectono-magmatic history since the Cretaceous. The relationship between shallow and deep structures in this complex tectonic domain is not straightforward. We present a 2D velocity model constructed by the assemblage of 1D shear wave velocity (Vs) models from 26 seismic stations, mainly distributed along the southern flank of the Alborz Mountains. The shear wave velocity structure has been estimated beneath each station using joint inversion of P-waves receiver functions and Rayleigh wave dispersion curves. A substantiation of the Vs inversion results sits on the modeling of Bouguer gravity anomaly data. Our velocity and density models show low velocity/density anomalies in uppermost mantle of western and central Alborz at a depth range of ∼50-100 km. In deeper parts of the uppermost mantle (depth range of 100-150 km), a high velocity/density anomaly is located beneath most of the Mountain range. The spatial pattern of these low and high velocity/density structures in the upper mantle is interpreted as the result of post collisional delamination of lower part of the western and central Alborz lithosphere.

Can Nitrogen be a Candidate for the Fe-Core Formation?

NASA Astrophysics Data System (ADS)

Dobrzhinetskaya, L.; Wirth, R.; Yang, J.; Weber, P.; Hutcheon, I.; Green, H. W.

2008-12-01

Among the light elements that have been added to mineral physics experiments concerning the Fe-rich core of the Earth, nitrogen is less favorable. In general, this is because metal-nitrides are thought to be rare within Earth. However this may not be because they are rare, but because nitrogen is difficult to detect by conventional electron microprobe analysis unless one is specifically looking for it. Theoretically, metal-nitrides could be equally considered as potential candidates for the light element in the core, not only because nitrogen forms strong metallic bonds, but also because metal-nitrides are common constituents of many iron meteorites. Some Fe-nitrides are found to be stable at extreme pressures and temperatures corresponding to Earth's core in both diamond anvil cell and shock experiments (Adler and Williams, 2005; Sekine et al., 2007). We have discovered a metal-nitride phase, TiN (osbornite) within a mantle mineralogical assemblage, opening a new opportunity to understand the history of Earth's core formation. The TiN was found in the mantle section of an unmetamorphosed Tibetan ophiolite, a fragment of former mid-ocean spreading center, which now marks the tectonic boundary between Asia and India. The osbornite occurs as inclusions in coesite pseudomorphic after stishovite, in association with FeTi alloy, native Fe, TiO2 II, cubic BN and diamond included in Os-Ir alloy, all from a massive chromitite ore body enclosed within harzburgite (Yang et al., 2007; Dobrzhinetskaya et al., 2007). The chromite also exhibits coesite and diopside exsolution lamellae (Yamamoto et al., 2007) that might suggest the calcium-ferrite polymorph of chromite as a precursor decompressed during upwelling. Measurements of δ15N with a Cameca 50 NanoSIMS using the same Focused Ion Beam foils prepared and used for earlier TEM studies suggest that the Tibetan osbornite is characterized by negative δ15N (-10 ‰). The δ15N results from the Tibetan osbornite are somewhat more negative than the most commonly measured value for Earth's uppermost mantle (δ15N = -3 to -5 ‰), and they are clearly different from the δ15N of shallow reservoirs. The latter include atmosphere, ocean, and crust having values of delta δ15N -- 0 - +5 ‰ for the atmosphere and ocean and +5 - +12 ‰.) for the crustal rocks and sediments. We conclude that the Tibetan osbornite contains mantle N, perhaps from an old and/or deep mantle reservoir. Apropos of the suggestion of N in the core, we point out that most iron meteorites have extremely negative δ15N values of -60 ‰ or more, hence it is conceivable that part of the N signal in our materials comes from a leaky core.

Ferropericlase inclusions in ultradeep diamonds from Sao Luiz (Brazil): high Li abundances and diverse Li-isotope and trace element compositions suggest an origin from a subduction mélange

NASA Astrophysics Data System (ADS)

Seitz, Hans-Michael; Brey, Gerhard P.; Harris, Jeffrey W.; Durali-Müller, Soodabeh; Ludwig, Thomas; Höfer, Heidi E.

2018-05-01

The most remarkable feature of the inclusion suite in ultradeep alluvial and kimberlitic diamonds from Sao Luiz (Juina area in Brazil) is the enormous range in Mg# [100xMg/(Mg + Fe)] of the ferropericlases (fper). The Mg-richer ferropericlases are from the boundary to the lower mantle or from the lower mantle itself when they coexist with ringwoodite or Mg- perovskite (bridgmanite). This, however, is not an explanation for the more Fe-rich members and a lowermost mantle or a "D" layer origin has been proposed for them. Such a suggested ultra-deep origin separates the Fe-rich fper-bearing diamonds from the rest of the Sao Luiz ultradeep diamond inclusion suite, which also contains Ca-rich phases. These are now thought to have an origin in the uppermost lower mantle and in the transition zone and to belong either to a peridotitic or mafic (subducted oceanic crust) protolith lithology. We analysed a new set of more Fe-rich ferropericlase inclusions from 10 Sao Luiz ultradeep alluvial diamonds for their Li isotope composition by solution MC-ICP-MS (multi collector inductively coupled plasma mass spectrometry), their major and minor elements by EPMA (electron probe micro-analyser) and their Li-contents by SIMS (secondary ion mass spectrometry), with the aim to understand the origin of the ferropericlase protoliths. Our new data confirm the wide range of ferropericlase Mg# that were reported before and augment the known lack of correlation between major and minor elements. Four pooled ferropericlase inclusions from four diamonds provided sufficient material to determine for the first time their Li isotope composition, which ranges from δ7Li + 9.6 ‰ to -3.9 ‰. This wide Li isotopic range encompasses that of serpentinized ocean floor peridotites including rodingites and ophicarbonates, fresh and altered MORB (mid ocean ridge basalt), seafloor sediments and of eclogites. This large range in Li isotopic composition, up to 5 times higher than `primitive upper mantle' Li-abundances, and an extremely large and incoherent range in Mg# and Cr, Ni, Mn, Na contents in the ferropericlase inclusions suggests that their protoliths were members of the above lithologies. This mélange of altered rocks originally contained a variety of carbonates (calcite, magnesite, dolomite, siderite) and brucite as the secondary products in veins and as patches and Ca-rich members like rodingites and ophicarbonates. Dehydration and redox reactions during or after deep subduction into the transition zone and the upper parts of the lower mantle led to the formation of diamond and ferropericlase inclusions with variable compositions and a predominance of the Ca-rich, high-pressure silicate inclusions. We suggest that the latter originated from peridotites, mafic rocks and sedimentary rocks as redox products between calcite and SiO2.

Crustal and uppermost mantle S-wave velocity structure beneath the Japanese islands from seismic ambient noise tomography

NASA Astrophysics Data System (ADS)

Guo, Zhi; Gao, Xing; Shi, Heng; Wang, Weiming

2013-04-01

In this study, the crustal and uppermost mantle shear wave velocities beneath the Japanese islands have been determined by inversion from seismic ambient noise tomography using data recorded at 75 Full Range Seismograph Network of Japan broad-band seismic stations, which are uniformly distributed across the Japanese islands. By cross-correlating 2 yr of vertical component seismic ambient noise recordings, we are able to extract Rayleigh wave empirical Green's functions, which are subsequently used to measure phase velocity dispersion in the period band of 6-50 s. The dispersion data are then inverted to yield 2-D tomographic phase velocity maps and 3-D shear wave velocity models. Our results show that the velocity variations at short periods (˜10 s), or in the uppermost crust, correlate well with the major known surface geological and tectonic features. In particular, the distribution of low-velocity anomalies shows good spatial correlation with active faults, volcanoes and terrains of sediment exposure, whereas the high-velocity anomalies are mainly associated with the mountain ranges. We also observe that large upper crustal earthquakes (5.0 ≤ M ≤ 8.0, depth ≤ 25 km) mainly occurred in low-velocity anomalies or along the boundary between low- and high-velocity anomalies, suggesting that large upper crustal earthquakes do not strike randomly or uniformly; rather they are inclined to nucleate within or adjacent to low-velocity areas.

Sensitivities of seismic velocities to temperature, pressure and composition in the lower mantle

NASA Astrophysics Data System (ADS)

Trampert, Jeannot; Vacher, Pierre; Vlaar, Nico

2001-08-01

We calculated temperature, pressure and compositional sensitivities of seismic velocities in the lower mantle using latest mineral physics data. The compositional variable refers to the volume proportion of perovskite in a simplified perovskite-magnesiowüstite mantle assemblage. The novelty of our approach is the exploration of a reasonable range of input parameters which enter the lower mantle extrapolations. This leads to realistic error bars on the sensitivities. Temperature variations can be inferred throughout the lower mantle within a good degree of precision. Contrary to the uppermost mantle, modest compositional changes in the lower mantle can be detected by seismic tomography, with a larger uncertainty though. A likely trade-off between temperature and composition will be largely determined by uncertainties in tomography itself. Given current sources of uncertainties on recent data, anelastic contributions to the temperature sensitivities (calculated using Karato's approach) appear less significant than previously thought. Recent seismological determinations of the ratio of relative S to P velocity heterogeneity can be entirely explain by thermal effects, although isolated spots beneath Africa and the Central Pacific in the lowermost mantle may ask for a compositional origin.

Upper mantle heterogeneity: Comparisons of regions south of Australia with Philippine Basin

NASA Technical Reports Server (NTRS)

1982-01-01

The nature of mass anomalies that occur beneath the regions of negative residual depth anomalies were identified. Residual geoid anomalies with negative residual depth anomalies are identified in the Philippine Basin (negative) and in the region south of Australia (positive and negative). In the latter region the geoid anomalies are eastward and the depth anomaly is northeast. It is suggested that the negative depth anomaly and the compensating mass excess in the uppermost mantle developed in the Eocene as the lithosphere of the west Philippine basin formed. Heating of the deeper upper mantle which causes slow surface wave velocities and negative gravity and geoid anomalies may be a younger phenomenon which is still in progress.

Clinopyroxenite dykes within a banded unit in the basal mantle section of the northern part of the Oman ophiolite: A record of the latest deep-seated magmatism

NASA Astrophysics Data System (ADS)

Ishimaru, Satoko; Arai, Shoji; Tamura, Akihiro

2017-11-01

We found clinopyroxenite dykes in a banded harzburgite block within the Sumeini area in the uppermost part of the metamorphic sole of the northern part of the Oman ophiolite. The dykes clearly cut the deformational structure of the harzburgite and contain its fragments, indicating dyke formation during obduction of the ophiolite. The Mg# [= Mg / (Mg + total Fe)] of clinopyroxenes in the dykes ranges from 0.81 to 0.91, and increases up to 0.93 proximal to harzburgite fragments. Mantle minerals in the harzburgite fragments were modified chemically through interaction with the magma that formed the dyke, yielding lower clinopyroxene and spinel Mg#, and spinels with higher TiO2 contents than those in the unaltered harzburgite. These geochemical features indicate that the clinopyroxenite dykes are cumulates derived from a relatively deep-seated primitive magma enriched in light rare earth elements (LREE) with an ocean island basalt (OIB)-like affinity, geochemically similar to the V3 lavas of an off-ridge origin. Combining these data with geological observations suggests that the clinopyroxenite dykes represent root system of the V3 lavas. Our analyses of the clinopyroxenite dykes testify to the external nature of the V3 magmas, which was added to the sliced oceanic lithosphere from the outside. It is likely that the V3 magma underwent deep-seated crystallization of clinopyroxene and had limited interaction with mantle peridotite en route to the surface. The mode of occurrence of the Sumeini clinopyroxenites (i.e., emplaced into a banded harzburgite block surrounded by garnet amphibolite) is consistent with the generation of OIB-like magmas (V3 lava) beneath the Oman ophiolite resulting from the break-off of the "subducting slab" and subsequent infiltration of hot asthenospheric mantle. This view is consistent with the limited distribution of V3-related rocks in the Oman ophiolite. The production of such OIB-like magmas during ophiolite obduction is not a rare event, especially during the subduction of young and hot oceanic lithosphere. Footnote: Mg#, Mg/(Mg + Fe total) atomic ratio; Cr#, Cr/(Cr + Al) atomic ratio; YAl, Al/(Al + Cr + Fe3 +) atomic ratio; YCr, Al/(Al + Cr + Fe3 +) atomic ratio; YFe, Al/(Al + Cr + Fe3 +) atomic ratio. RSD: relative standard deviation. Ref. values: Reference values of JEOL_Kfs from a JEOL database.

Decrease in oceanic crustal thickness since the breakup of Pangaea

NASA Astrophysics Data System (ADS)

van Avendonk, Harm J. A.; Davis, Joshua K.; Harding, Jennifer L.; Lawver, Lawrence A.

2017-01-01

Earth's mantle has cooled by 6-11 °C every 100 million years since the Archaean, 2.5 billion years ago. In more recent times, the surface heat loss that led to this temperature drop may have been enhanced by plate-tectonic processes, such as continental breakup, the continuous creation of oceanic lithosphere at mid-ocean ridges and subduction at deep-sea trenches. Here we use a compilation of marine seismic refraction data from ocean basins globally to analyse changes in the thickness of oceanic crust over time. We find that oceanic crust formed in the mid-Jurassic, about 170 million years ago, is 1.7 km thicker on average than crust produced along the present-day mid-ocean ridge system. If a higher mantle temperature is the cause of thicker Jurassic ocean crust, the upper mantle may have cooled by 15-20 °C per 100 million years over this time period. The difference between this and the long-term mantle cooling rate indeed suggests that modern plate tectonics coincide with greater mantle heat loss. We also find that the increase of ocean crustal thickness with plate age is stronger in the Indian and Atlantic oceans compared with the Pacific Ocean. This observation supports the idea that upper mantle temperature in the Jurassic was higher in the wake of the fragmented supercontinent Pangaea due to the effect of continental insulation.

Joint Inversion of Crustal and Uppermost Mantle Structure in Western China

DTIC Science & Technology

2013-11-02

western China from this project. Approved for public release; distribution is unlimited. 21 References Ammon, C. J., G. E. Randall , and G. Zandt, On...the non-uniqueness of receiver function inversions, J. Geophys. Res., 95(B10), pp. 15303-15318, 1990. Benson , G., et al., Processing seismic

Evolution of the lithosphere in the Indian Ocean from combined earthquake and ambient noise tomography

NASA Astrophysics Data System (ADS)

Ma, Zhitu; Dalton, Colleen A.

2017-01-01

Rayleigh wave dispersion extracted from ambient seismic noise has been widely used to image crustal and uppermost mantle structure. Applications of this approach in continental settings are abundant, but there have been relatively few studies within ocean basins. By cross-correlating ambient noise recorded at broadband seismic stations around the Indian Ocean, we demonstrate the feasibility of extracting high-quality, long-period (10-30 mHz) Rayleigh waves that traverse the entire ocean basin. High-quality Rayleigh wave cross-correlation functions can be obtained from stacking waveforms over less than 2 years at land stations and less than 4 years at island stations. We show that adding the dispersion information extracted from ambient noise to a global earthquake data set can improve the resolution of phase velocity maps by about 20% in the northern Indian Ocean, where the station distribution is the best. We find that a plate cooling model with a potential temperature of 1450°C and plate thickness of 125 km can fit both the seismic observations and seafloor topography. The Seychelles-Mascarene Plateau is characterized by anomalously slow velocity at 30 mHz. The inclusion of ambient noise data in the tomographic inversion shifts the slow velocity anomaly into better agreement with the topographic relief, allowing us to estimate its crustal thickness and confirm that the plateau's elevation is supported by thick crust. The 10 and 20 mHz phase velocity maps show a strong asymmetry across the Central Indian Ridge that is best explained by eastward asthenospheric flow emanating from nearby hot spots.

How broad and deep is the region of chemical alteration of oceanic plates at trenches?

NASA Astrophysics Data System (ADS)

Ranero, C. R.; Grevemeyer, I.; Barckhausen, U.

2017-12-01

Different lines of evidence indicate that oceanic plates are affected by pervasive bending-related deformation approaching ocean trenches. Results from active-seismic work support that deformation provides paths for exchange between hydrosphere and lithosphere, possibly causing chemical alteration of the incoming lithosphere. Much work focused on the potential transformation of peridotite to serpentine in the uppermost mantle of incoming plates, but there is no consensus on the region where it may occur or the intensity of alteration, let alone on limiting factors for the process. Teleseismic (large-great) earthquakes with normal-fault mechanism in the outer rise region have been often called to speculate on the depth of penetration of plate hydration. However, large-great outer-rise earthquakes may be related to stress changes due to slab pull after decoupling along the inter-plate boundary, and not necessarily controlled by bending stresses only. If so, the majority of the time the depth of water percolation may be related to local bending stresses expressed by micro-earthquakes rather than large events. Seismic images and multibeam bathymetry from lithosphere of similar thermal thickness from different trenches display a remarkable variability of the intensity of bending-related deformation along the subduction zones where plate age does not change significantly indicating that the intensity of deformation (not the depth) and perhaps hydration is very variable in space and not controlled by plate age. Seismic images showing hundreds of kilometers perpendicular to the trench into the incoming plate show that the bending-related deformation reaches mantle under the outer rise, well before the lithosphere plunges into the trench and develops the marked bend-faulting fabric observable in bathymetric maps. Thus, alteration occurs in a hundreds-of-km wide area, with deformation intensity related to local characteristics, and deformation depth to plate age.

Absolute Plate Motion Control Since the Triassic from the Cocos Slab and its Associated Subduction Record in Mexico

NASA Astrophysics Data System (ADS)

Boschman, L.; Van Hinsbergen, D. J. J.; Langereis, C. G.; Molina-Garza, R. S.; Kimbrough, D. L.; Spakman, W.

2017-12-01

A positive wave speed anomaly interpreted as the Cocos slab stretches from the uppermost mantle at the Middle America trench in the west, to the lowermost mantle below the Atlantic in the east. The length and continuity of this slab indicates long-lived, uninterrupted eastward subduction of the attached Cocos Plate and its predecessor, the Farallon Plate. The geological record of Mexico contains Triassic to present day evidence of subduction, of which the post-Late Cretaceous phase is of continental margin-style. Interpretations of the pre-Upper Cretaceous subduction-related rock assemblages are under debate, and vary from far-travelled exotic intra-oceanic island arc character to in-situ extended continental margin origin. We present new paleomagnetic data that show that Triassic, Jurassic and Cretaceous subduction-related rocks from the Vizcaíno Peninsula and the Guerrero terrane have a paleolatitudinal plate motion history that is equal to that of the North American continent. This suggests that these rock assemblages were part of the overriding plate and were perhaps only separated from the North American continent by temporal fore- or back-arc spreading. The entire Triassic-present day subduction record, and hence, reconstructed trench location, can therefore be linked to the Cocos slab, which provides control on longitudinal plate motion of North America since the time of Pangea. Compared to the latest state of the art mantle frames, in which longitudes are essentially unconstrained for pre-Cretaceous times, our reconstructed absolute position of North America requires a significant westward longitudinal shift for Mesozoic times.

Unveiling the lithospheric structure of the US Interior using the USArray Transportable Array

NASA Astrophysics Data System (ADS)

Moschetti, M. P.; Ritzwoller, M. H.; Lin, F.; Shen, W.; Yang, Y.

2009-12-01

We present current results from ambient noise tomography (ANT) and earthquake surface wave tomography applied to the USARRAY Transportable Array (TA) for the western and central US. We have processed ambient seismic noise data since October 2004 to produce cumulative Rayleigh and Love wave dispersion maps (from about 6 to 40 sec period) within the footprint of the TA. The high spatial density of these instruments results in dispersion maps with a resolution of about the average inter-station distance (70 km) and far exceeds previous surface wave tomographic results for the US interior. The dispersion maps from ANT are complemented by Rayleigh wave phase speed maps from teleseismic earthquake tomography (25 - 100 sec period). The development of a new method of surface wave tomography, termed Eikonal tomography, that models wavefront complexity and off great-circle propagation allows for the robust estimation of phase velocity azimuthal anisotropy. Eikonal tomography has been applied to ambient seismic noise and earthquake measurements and provides a means to compare and vet results in the period band of overlap (25 - 40 sec). In addition, the recent application of this method to Love waves from teleseismic earthquakes provides dispersion measurements up to 50 sec period. These longer period Love wave dispersion measurements may improve the characterization of anisotropy in the uppermost mantle. In addition to the current dispersion maps, we present regional-scale 3-D models of isotropic and anisotropic shear-velocities for the crust and uppermost mantle beneath the western US. Because dispersion measurements from ambient seismic noise include short period (<20 sec) information, they provide a strong constraint on the shear-velocity structure of the crust and uppermost mantle. A radially anisotropic shear-velocity model of the crust and uppermost mantle is constructed by simultaneously inverting Rayleigh and Love wave dispersion measurements from ANT and from earthquake tomography. Models with isotropic and radially anisotropic mantle shear-velocities do not fit the Rayleigh and Love wave measurements simultaneously across large regions of the western US, and the models present a Rayleigh-Love misfit discrepancy at the periods most sensitive to crustal velocity structures. However, by introducing positive radial anisotropy (Vsh>Vsv) to the middle and lower crust, this misfit discrepancy is resolved. Higher amplitude crustal radial anisotropy is observed in the predominant extensional provinces of the western US and is thought to result from the alignment of anisotropic crustal minerals during extension and deformation. Several regions of the western US remain poorly fit by the 3-D radially anisotropic shear-velocity model. These include the Olympic Peninsula, Mendocino Triple Junction, southern Cascadia backarc, Yakima Fold Belt, Wasatch Front, Salton Trough and Great Valley. We investigate various additional model parametrizations and the effect of breaking the constraint on the monotonic increase of crustal velocities with depth to resolve crustal shear-velocity structure in these regions. These techniques will readily be applied to data from the US Interior as the TA moves to the east.

Age-Orientation Relationships of Northern Hemisphere Martian Gullies and "Pasted-on" Mantling Unit: Implications for Near-Surface Water Migration in Mars' Recent History

NASA Technical Reports Server (NTRS)

Bridges, N. T.; Lackner, C. N.

2005-01-01

The finding of abundant, apparently young, Martian gullies with morphologies indicative of formation by flowing fluid was surprising in that volumes of near-surface liquid water in sufficient quantities to modify the surface geology were not thought possible under current conditions. Original hypotheses on origin of gullies were mostly centered on groundwater seepage and surface runoff and melting of near-surface ground ice. More recently, melting of snow deposited in periods of higher obliquity has been proposed as a possible origin of the gullies. Tied to this hypothesis is the supposition that the "pasted-on" mantling unit observed in association with many gullies is composed of remnant snowpack. The mantling unit has distinct rounded edge on its upper boundary and exhibits features suggestive of flow noted that the uppermost part of the mantle marks where gullies begin, suggesting that the source of water for the gullies was within the mantle. The mantle is found preferentially on cold, pole-facing slopes and, where mantled and non-mantled slopes are found together, gullies are observed incised into the latter. In other cases, the mantling material lacks gullies.

Origin of the Indian Ocean-type isotopic signature in basalts from Philippine Sea plate spreading centers: An assessment of local versus large-scale processes

NASA Astrophysics Data System (ADS)

Hickey-Vargas, Rosemary

1998-09-01

Basalts erupted from spreading centers on the Philippine Sea plate between 50 Ma and the present have the distinctive isotopic characteristics of Indian Ocean mid-ocean ridge basalt (MORB), such as high 208Pb/204Pb and low 143Nd/144Nd for a given 206Pb/204Pb compared with Pacific and Atlantic Ocean MORB. This feature may indicate that the upper mantle of the Philippine Sea plate originated as part of the existing Indian Ocean upper mantle domain, or, alternatively, that local processes duplicated these isotopic characteristics within the sub-Philippine Sea plate upper mantle. Synthesis of new and published isotopic data for Philippine Sea plate basin basalts and island arc volcanic rocks, radiometric ages, and tectonic reconstructions of the plate indicates that local processes, such as contamination of the upper mantle by subducted materials or by western Pacific mantle plumes, did not produce the Indian Ocean-type signature in Philippine Sea plate MORB. It is more likely that the plate originated over a rapidly growing Indian Ocean upper mantle domain that had spread into the area between Australia/New Guinea and southeast Asia before 50 Ma.

Upper mantle electrical resistivity structure beneath back-arc spreading centers

NASA Astrophysics Data System (ADS)

Seama, N.; Shibata, Y.; Kimura, M.; Shindo, H.; Matsuno, T.; Nogi, Y.; Okino, K.

2011-12-01

We compare four electrical resistivity structure images of the upper mantle across back-arc spreading centers (Mariana Trough at 18 N and 13 N, and the Eastern Lau at 19.7 S and 21.3 S) to provide geophysical constraints on issues of mantle dynamics beneath the back-arc spreading system related to the subducting slab. The central Mariana Trough at 18 N has the full spreading rate of 25 km/Myr, and shows characteristic slow-spreading features; existence of median valley neovolcanic zone and "Bull's eyes" mantle Bouguer anomaly (MBA) along the axes. On the other hand, the southern Mariana Trough at 13 N shows an EPR type axial relief in morphology and lower MBA than that in the central Mariana Trough (Kitada et al., 2006), suggesting abundance of magma supply, even though the full spreading rate is 35 km/Myr that is categorized as a slow spreading ridge. At the Eastern Lau spreading center, crustal thickness and morphology vary systematically with arc proximity and shows the opposed trends against spreading rate: The full spreading rate increases from 65 km/Myr at 21.3 S to 85 km/Myr at 19.7 S, while the crustal thicknesses decrease together with morphology transitions from shallow peaked volcanic highs to a deeper flat axis (Martinez et al., 2006). Matsuno et al. (2010) provides a resistivity structure image of the upper mantle across the central Mariana subduction system, which contains several key features: There is an uppermost resistive layer with a thickness of 80-100 km beneath the central Mariana Trough, suggesting dry residual from the plate accretion process. But there is no evidence for a conductive feature beneath the back-arc spreading center at 18 N, and this feature is clearly independent from the conductive region beneath the volcanic arc below 60 km depth that reflects melting and hydration driven by water release from the subducting slab. The resultant upper mantle resistivity structure well support that the melt supply is not abundant, resulting in characteristic slow-spreading features at the surface. We have conducted marine magnetotelluric (MT) surveys at the southern Mariana in 2010 and at the Eastern Lau in 2009-2010. We obtained 10 ocean bottom electro-magnetometer (OBEM) data from a 130 km length MT transect across the southern Mariana spreading axis at 13 N, while we obtained 2 OBEM data and 11 ocean bottom magnetometer data from two 160 km length MT transects across the Eastern Lau spreading axes at 19.7 S and 21.3 S. After calculation of MT response functions and their correction for topographic distortion, two-dimensional electrical resistivity structures will be derived using an inversion algorithm. At this meeting, first we will show the resistivity structure images of the upper mantle beneath these spreading axes. Then, these structure images will be compared to identify differences in the mantle dynamics and the melt supply beneath the back-arc spreading system related to the subducting slab.

Seismic evidence for broad attenuation anomalies in the asthenosphere beneath the Pacific Ocean

NASA Astrophysics Data System (ADS)

Adenis, Alice; Debayle, Eric; Ricard, Yanick

2017-06-01

We present QADR17, a global model of Rayleigh-wave attenuation based on a massive surface wave data set (372 629 frequency-dependent attenuation curves in the period range 50-260 s). We correct for focusing-defocusing effects and geometrical spreading, and perform a stringent selection to only keep robust observations. Then, data with close epicentres recorded at the same station are clustered, as they sample the same Earth's structure. After this pre-selection, our data set consists of about 35 000 curves that constrain the Rayleigh-wave intrinsic attenuation in the upper mantle. The logarithms of the attenuation along the individual rays are then inverted to obtain global maps of the logarithm of the local attenuation. After a first inversion, outliers are rejected and a second inversion yields a variance reduction of about 45 per cent. Our attenuation maps present strong agreement with surface tectonics at periods lower than 200 s, with low attenuation under continents and high attenuation under oceans. Over oceans, attenuation decreases with increasing crustal ages, but at periods sensitive to the uppermost 150 km, mid-ocean ridges are not characterized by a very localized anomaly, in contrast to what is commonly observed for seismic velocity models. Attenuation is rather well correlated with hotspots, especially in the Pacific ocean, where a strong attenuating anomaly is observed in the long wavelength component of our signal at periods sampling the oceanic asthenosphere. We suggest that this anomaly results from the horizontal spreading of several thermal plumes within the asthenosphere. Strong velocity reductions associated with high attenuation anomalies of moderate amplitudes beneath the East Pacific Rise, the Red Sea and the eastern part of Asia may require additional mechanisms, such as partial melting.

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

NASA Astrophysics Data System (ADS)

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

2014-12-01

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

Platinum-group elements, S, Se and Cu in highly depleted abyssal peridotites from the Mid-Atlantic Ocean Ridge (ODP Hole 1274A): Influence of hydrothermal and magmatic processes

NASA Astrophysics Data System (ADS)

Marchesi, Claudio; Garrido, Carlos J.; Harvey, Jason; González-Jiménez, José María; Hidas, Károly; Lorand, Jean-Pierre; Gervilla, Fernando

2013-11-01

Highly depleted harzburgites and dunites were recovered from ODP Hole 1274A, near the intersection between the Mid-Atlantic Ocean Ridge and the 15°20'N Fracture Zone. In addition to high degrees of partial melting, these peridotites underwent multiple episodes of melt-rock reaction and intense serpentinization and seawater alteration close to the seafloor. Low concentrations of Se, Cu and platinum-group elements (PGE) in harzburgites drilled at around 35-85 m below seafloor are consistent with the consumption of mantle sulfides after high degrees (>15-20 %) of partial melting and redistribution of chalcophile and siderophile elements into PGE-rich residual microphases. Higher concentrations of Cu, Se, Ru, Rh and Pd in harzburgites from the uppermost and lowest cores testify to late reaction with a sulfide melt. Dunites were formed by percolation of silica- and sulfur-undersaturated melts into low-Se harzburgites. Platinum-group and chalcophile elements were not mobilized during dunite formation and mostly preserve the signature of precursor harzburgites, except for higher Ru and lower Pt contents caused by precipitation and removal of platinum-group minerals. During serpentinization at low temperature (<250 °C) and reducing conditions, mantle sulfides experienced desulfurization to S-poor sulfides (mainly heazlewoodite) and awaruite. Contrary to Se and Cu, sulfur does not record the magmatic evolution of peridotites but was mostly added in hydrothermal sulfides and sulfate from seawater. Platinum-group elements were unaffected by post-magmatic low-temperature processes, except Pt and Pd that may have been slightly remobilized during oxidative seawater alteration.

Deep seismic reflection images of the Wharton Basin oceanic crust and uppermost mantle offshore Northern Sumatra: Relation with active and past deformation

NASA Astrophysics Data System (ADS)

Carton, Hélène; Singh, Satish C.; Hananto, Nugroho D.; Martin, James; Djajadihardja, Yusuf S.; Udrekh; Franke, Dieter; Gaedicke, Christoph

2014-01-01

present deep seismic reflection images along two profiles collected in 2006 in the Wharton Basin offshore Northern Sumatra. The main profile is located subparallel to the Sumatran trench at a distance of 32-66 km. Faulting of the entire sedimentary section (strike-slip deformation sometimes accompanied by a dip-slip component) is imaged over two fracture zones of the extinct Wharton Spreading Center that prior studies have shown to be reactivated as left-lateral faults. The western fracture zone is associated with a wide region of strong basement topography, a difference in crustal thickness of 1.5 km, and an age offset of 9 Ma. The epicenters of the 11 April 2012 Mw 8.6 great strike-slip earthquake, its Mw 7.2 foreshock, and Mw 8.2 aftershock align along this major structure > 100 km south of the profile intersection. Our high-quality long-offset seismic reflection data also reveal bright dipping reflections extending down to a maximum of 24 km into the oceanic mantle ( 37 km below sea level). Apparent dips are mostly 25-35°, corresponding to 30-55° along either N-S to NNE-SSW or E-W to WNW-ESE directions, which encompass the directions of plate fabric and nodal planes of the Mw 8.6 event. We suggest that these enigmatic reflections arise from presently inactive dip-slip fault planes reaching for the deepest ones to the base of the brittle layer. Possible origins include extension related to plate bending or an episode of now inactive thrust-type deformation reactivating paleonormal faults, similar to that taking place in the Central Indian Basin.

Magnetotelluric investigations of the lithosphere beneath the central Rae craton, mainland Nunavut, Canada

NASA Astrophysics Data System (ADS)

Spratt, Jessica E.; Skulski, Thomas; Craven, James A.; Jones, Alan G.; Snyder, David B.; Kiyan, Duygu

2014-03-01

New magnetotelluric soundings at 64 locations throughout the central Rae craton on mainland Nunavut constrain 2-D resistivity models of the crust and lithospheric mantle beneath three regional transects. Responses determined from colocated broadband and long-period magnetotelluric recording instruments enabled resistivity imaging to depths of > 300 km. Strike analysis and distortion decomposition on all data reveal a regional trend of 45-53°, but locally the geoelectric strike angle varies laterally and with depth. The 2-D models reveal a resistive upper crust to depths of 15-35 km that is underlain by a conductive layer that appears to be discontinuous at or near major mapped geological boundaries. Surface projections of the conductive layer coincide with areas of high grade, Archean metasedimentary rocks. Tectonic burial of these rocks and thickening of the crust occurred during the Paleoproterozoic Arrowsmith (2.3 Ga) and Trans-Hudson orogenies (1.85 Ga). Overall, the uppermost mantle of the Rae craton shows resistivity values that range from 3000 Ω m in the northeast (beneath Baffin Island and the Melville Peninsula) to 10,000 Ω m beneath the central Rae craton, to >50,000 Ω m in the south near the Hearne Domain. Near-vertical zones of reduced resistivity are identified within the uppermost mantle lithosphere that may be related to areas affected by mantle melt or metasomatism associated with emplacement of Hudsonian granites. A regional decrease in resistivities to values of 500 Ω m at depths of 180-220 km, increasing to 300 km near the southern margin of the Rae craton, is interpreted as the lithosphere-asthenosphere boundary.

Upper mantle P velocity structure beneath the Baikal Rift from modeling regional seismic data

NASA Astrophysics Data System (ADS)

Brazier, Richard A.; Nyblade, Andrew A.

2003-02-01

Uppermost mantle P wave velocity structure beneath the Baikal rift and southern margin of the Siberian Platform has been investigated by using a grid search method to model Pnl waveforms from two moderate earthquakes recorded by station TLY at the southwestern end of Lake Baikal. The results yielded a limited number of successful models which indicate the presence of upper mantle P wave velocities beneath the rift axis and the margin of the platform that are 2-5% lower than expected. The magnitude of the velocity anomalies and their location support the presence of a thermal anomaly that extends laterally beyond the rift proper, possibly created by small-scale convection or a plume-like, thermal upwelling.

Crustal and uppermost mantle structures of the South China from joint analysis of receiver functions and Rayleigh wave dispersions

NASA Astrophysics Data System (ADS)

Guo, Zhi; Gao, Xing; Li, Tong; Wang, Wei

2018-05-01

We use P-wave receiver function H-k stacking and joint inversion of receiver functions and Rayleigh wave dispersions to investigate crustal and uppermost mantle structure beneath the South China. The obtained results reveal prominent crustal structure variations in the study area, Moho depth increases from ∼30 km in the Cathaysia Block to more than ∼60 km in the eastern Tibetan Plateau. A Moho undulation and Vp/Vs ratio variations can be observed from the Cathaysia Block to Yangtze Craton. These observations consistent with the crustal structures predict by the flat slab subduction model. We interpret these lateral crustal structure variations reflect the tectonic evolution of the Yangtze Craton and Cathaysia Block prior the Mesozoic and the post-orogenic magmatism due to the breaking up of the subducted flat slab and subsequent slab rollback in the South China. The observed variations of the crustal structures not only reveal the lateral crustal inhomogeneity, but also provide constraints on the geodynamic evolution of the South China.

Oceanic crust recycling and the formation of lower mantle heterogeneity

NASA Astrophysics Data System (ADS)

van Keken, Peter E.; Ritsema, Jeroen; Haugland, Sam; Goes, Saskia; Kaneshima, Satoshi

2016-04-01

The Earth's lower mantle is heterogeneous at multiple scales as demonstrated for example by the degree-2 distribution of LLSVPs seen in global tomography and widespread distribution of small scale heterogeneity as seen in seismic scattering. The origin of this heterogeneity is generally attributed to leftovers from Earth's formation, the recycling of oceanic crust, or a combination thereof. Here we will explore the consequences of long-term oceanic crust extraction and recycling by plate tectonics. We use geodynamical models of mantle convection that simulate plates in an energetically consistent manner. The recycling of oceanic crust over the age of the Earth produces persistent lower mantle heterogeneity while the upper mantle tends to be significantly more homogeneous. We quantitatively compare the predicted heterogeneity to that of the present day Earth by tomographic filtering of the geodynamical models and comparison with S40RTS. We also predict the scattering characteristics from S-P conversions and compare these to global scattering observations. The geophysical comparison shows that lower mantle heterogeneity is likely dominated by long-term oceanic crust recycling. The models also demonstrate reasonable agreement with the geochemically observed spread between HIMU-EM1-DMM in ocean island basalts as well as the long-term gradual depletion of the upper mantle as observed in Lu-Hf systematics.

Shear wave velocity structure of the Anatolian Plate and surrounding regions using Ambient Noise Tomography

NASA Astrophysics Data System (ADS)

Delph, J. R.; Beck, S. L.; Zandt, G.; Biryol, C. B.; Ward, K. M.

2013-12-01

The Anatolian Plate consists of various lithospheric terranes amalgamated during the closure of the Tethys Ocean, and is currently extruding to the west in response to a combination of the collision of the Arabian plate in the east and the roll back of the Aegean subduction zone in the west. We used Ambient Noise Tomography (ANT) at periods <= 40s to investigate the crust and uppermost mantle structure of the Anatolian Plate. We computed a total of 13,779 unique cross-correlations using one sample-per-second vertical component broadband seismic data from 215 stations from 8 different networks over a period of 7 years to compute fundamental-mode Rayleigh wave dispersion curves following the method of Benson et al. (2007). We then inverted the dispersion data to calculate phase velocity maps for 11 periods from 8 s - 40 s throughout Anatolia and the Aegean regions (Barmin et al. 2001). Using smoothed Moho values derived from Vanacore et al. (2013) in our starting models, we inverted our dispersion curves using a linear least-squares iterative inversion scheme (Herrmann & Ammon 2004) to produce a 3-D shear-wave velocity model of the crust and uppermost mantle throughout Anatolia and the Aegean. We find a good correlation between our seismic shear wave velocities and paleostructures (suture zones) and modern deformation (basin formation and fault deformation). The most prominent crustal velocity contrasts occur across intercontinental sutures zones, resulting from the juxtaposition of the compositionally different basements of the amalgamated terranes. At shallow depths, seismic velocity contrasts correspond closely with surficial features. The Thrace, Cankiri and Tuz Golu basins, and accretionary complexes related to the closure of the Neotethys are characterized by slow shear wave velocities, while the Menderes and Kirsehir Massifs, Pontides, and Istanbul Zone are characterized by fast velocities. We find that the East Anatolia Plateau has slow shear-wave velocities, as expected due to high heat flow and active volcanism. The Tuz Golu fault has a visible seismic signal down to ~15 km below sea level, and the eastern Inner-Tauride Suture corresponding to the Central Anatolian Fault Zone may extend into the mantle. The Isparta Angle separates the actively extending portion of western Anatolia from the plateau regions in the east, and the largest anomaly (slow velocities) extending into the upper mantle is observed under the western flank of the Isparta Angle, corresponding to the Fethiye-Burdur fault zone. We attribute these slow shear-wave velocities to the effects of complex deformations within the crust as a result of the interactions of the African and Anatolian Plates. In the upper mantle, slow shear-wave velocities are consistent with a slab tear along a STEP fault corresponding to the extensions of the Pliny and Strabo Transform faults, allowing asthenosphere to rise to very shallow depths. The upper mantle beneath the Taurides exhibits very slow shear-wave velocities, in agreement with possible delamination or slab-breakoff (Cosentino et al. 2012) causing rapid uplift in the last 8 million years.

Can the composition and structure of the lower ocean crust and upper mantle be known without deep ocean drilling?

NASA Astrophysics Data System (ADS)

Dick, H.; Natland, J.

2003-04-01

No. With few exceptions, lower ocean crust sampled by dredge or submersible in tectonic windows such as Atlantis Bank in the Indian Ocean or the MARK area on the Mid-Atlantic Ridge are not representative of the ocean crust. They represent tectonic mixing of rocks from the mantle and crust on large faults that also localize late magmatic intrusion. Where this can be sorted out, the in-situ crustal sections may generally represent a sub-horizontal cross-section through the lower crust and mantle and not a vertical one. The gabbroic rocks exposed represent largely high-level intrusions, highly hybridized by late melt flow along deep faults, or highly evolved gabbro at the distal ends of larger intrusions emplaced into the mantle near transforms. Oceanic gabbros have average compositions that lie outside the range of primary MORB compositions, and rarely are equivalent to spatially associated MORB either as a parent to, or as a residue of their crystallization. Oceanic gabbros sampled from these complexes generally are very coarse-grained, and are unlike those seen in nearly all ophiolites and layered intrusions. In addition, there are few exposures of gabbro and lower ocean crust and mantle in Pacific tectonic windows, though there the possibility of more representative sections is greater due to their exposure in propagating rifts. Limited samples of the mantle from near the midpoints of ocean ridge segments at slow-spreading rifts are from anomalous crustal environments such as ultra-slow spreading ridges or failed rifts. These include abundant dunites, as opposed to samples from fracture zones, which contain only about 1% dunite. While this indicates focused mantle flow towards the midpoint of a ridge, it also shows that fracture zone peridotites are not fully representative of the oceanic upper mantle. Major classes of rocks common in ophiolites, such as fine to medium grained layered primitive olivine gabbros, troctolites, wherlites and dunites, sheeted dikes, and epidosites are rarely or even not exposed. Models of lower ocean crust stratigraphy drawn from deep sea sampling, certainly from slow spreading ridges, do not match those for major intact ophiolites. Thus the ophiolite hypothesis remains unconfirmed for the lower ocean crust and shallow mantle, and it is nearly impossible to accurately identify the ocean ridge environment of any one ophiolite. The one deep drill hole that exists in lower ocean crust, 1.5 km Hole 735B, has a bulk composition too fractionated to mass balance MORB back to a primary mantle melt composition. Thus, a large mass of primitive cumulates is missing and could be situated in the crust below the base of the hole or in the underlying mantle. This is an unresolved question that is critical to understanding the evolution of the most common magma on earth: MORB. Since lower ocean crust and mantle represent a major portion of the crust and the exchange of mass, heat and volatiles from the earth's interior to its exterior this leaves a major hole in our understanding of the global geochemical and tectonic cycle which can only be filled by deep drilling.

Metasomatic Enrichment of Oceanic Lithospheric Mantle Documented by Petit-Spot Xenoliths

NASA Astrophysics Data System (ADS)

Pilet, S.; Abe, N.; Rochat, L.; Hirano, N.; Machida, S.; Kaczmarek, M. A.; Muntener, O.

2015-12-01

Oceanic lithosphere is generally interpreted as mantle residue after MORB extraction. It has been proposed, however, that metasomatism could take place at the interface between the low-velocity zone and the cooling and thickening oceanic lithosphere or by the percolation of low-degree melts produced in periphery of Mid Ocean Ridges. This later process is observed in slow spreading ridges and ophiolites where shallow oceanic lithospheric mantle could be metasomatized/refertilized during incomplete MORB melt extraction. Nevertheless, direct evidence for metasomatic refertilization of the deep part of the oceanic lithospheric mantle is still missing. Xenoliths and xenocrysts sampled by petit-spot volcanoes interpreted as low-degree melts extracted from the base of the lithosphere in response to plate flexure, provide important new information about the nature and the processes associated with the evolution of oceanic lithospheric mantle. Here, we report, first, the presence of a garnet xenocryst in petit-spot lavas from Japan characterized by low-Cr, low-Ti content and mostly flat MREE-HREE pattern. This garnet is interpreted as formed during subsolidus cooling of pyroxenitic or gabbroic cumulates formed at ~1 GPa during the incomplete melt extraction at the periphery of the Pacific mid-ocean ridge. It is the first time that such processes are documented in fast spreading context. Second, we report petit-spot mantle xenoliths with cpx trace element "signatures" characterized by high U, Th, relative depletion in Nb, Pb, Ti and high but variable LREE/HREE ratio suggesting equilibration depth closed to the Gt/Sp transition zone. Such "signatures" are unknown from oceanic settings and show unexpected similarity to melt-metasomatized gt-peridotites sampled by kimberlites. This similarity suggests that metasomatic processes are not restricted to continental setting, but could correspond to a global mechanism at the lithosphere-asthenosphere boundary. As plate flexure represents a global mechanism in subduction zone, a portion of oceanic lithospheric mantle is likely to be metasomatized; recycling of these enriched domains into the convecting mantle is fundamental to understand the generation of small scale mantle isotopic and volatile heterogeneities sampled by OIBs and MORBs.

Effect of water on the fluorine and chlorine partitioning behavior between olivine and silicate melt.

PubMed

Joachim, Bastian; Stechern, André; Ludwig, Thomas; Konzett, Jürgen; Pawley, Alison; Ruzié-Hamilton, Lorraine; Clay, Patricia L; Burgess, Ray; Ballentine, Christopher J

2017-01-01

Halogens show a range from moderate (F) to highly (Cl, Br, I) volatile and incompatible behavior, which makes them excellent tracers for volatile transport processes in the Earth's mantle. Experimentally determined fluorine and chlorine partitioning data between mantle minerals and silicate melt enable us to estimate Mid Ocean Ridge Basalt (MORB) and Ocean Island Basalt (OIB) source region concentrations for these elements. This study investigates the effect of varying small amounts of water on the fluorine and chlorine partitioning behavior at 1280 °C and 0.3 GPa between olivine and silicate melt in the Fe-free CMAS+F-Cl-Br-I-H 2 O model system. Results show that, within the uncertainty of the analyses, water has no effect on the chlorine partitioning behavior for bulk water contents ranging from 0.03 (2) wt% H 2 O (D Cl ol/melt = 1.6 ± 0.9 × 10 -4 ) to 0.33 (6) wt% H 2 O (D Cl ol/melt = 2.2 ± 1.1 × 10 -4 ). Consequently, with the effect of pressure being negligible in the uppermost mantle (Joachim et al. Chem Geol 416:65-78, 2015), temperature is the only parameter that needs to be considered for the determination of chlorine partition coefficients between olivine and melt at least in the simplified iron-free CMAS+F-Cl-Br-I-H 2 O system. In contrast, the fluorine partition coefficient increases linearly in this range and may be described at 1280 °C and 0.3 GPa with ( R 2  = 0.99): [Formula: see text]. The observed fluorine partitioning behavior supports the theory suggested by Crépisson et al. (Earth Planet Sci Lett 390:287-295, 2014) that fluorine and water are incorporated as clumped OH/F defects in the olivine structure. Results of this study further suggest that fluorine concentration estimates in OIB source regions are at least 10% lower than previously expected (Joachim et al. Chem Geol 416:65-78, 2015), implying that consideration of the effect of water on the fluorine partitioning behavior between Earth's mantle minerals and silicate melt is vital for a correct estimation of fluorine abundances in OIB source regions. Estimates for MORB source fluorine concentrations as well as chlorine abundances in both mantle source regions are within uncertainty not affected by the presence of water.

Normal Mode Derived Models of the Physical Properties of Earth's Outer Core

NASA Astrophysics Data System (ADS)

Irving, J. C. E.; Cottaar, S.; Lekic, V.; Wu, W.

2017-12-01

Earth's outer core, the largest reservoir of metal in our planet, is comprised of an iron alloy of an uncertain composition. Its dynamical behaviour is responsible for the generation of Earth's magnetic field, with convection driven both by thermal and chemical buoyancy fluxes. Existing models of the seismic velocity and density of the outer core exhibit some variation, and there are only a small number of models which aim to represent the outer core's density.It is therefore important that we develop a better understanding of the physical properties of the outer core. Though most of the outer core is likely to be well mixed, it is possible that the uppermost outer core is stably stratified: it may be enriched in light elements released during the growth of the solid, iron enriched, inner core; by elements dissolved from the mantle into the outer core; or by exsolution of compounds previously dissolved in the liquid metal which will eventually be swept into the mantle. The stratified layer may host MAC or Rossby waves and it could impede communication between the chemically differentiated mantle and outer core, including screening out some of the geodynamo's signal. We use normal mode center frequencies to estimate the physical properties of the outer core in a Bayesian framework. We estimate the mineral physical parameters needed to best produce velocity and density models of the outer core which are consistent with the normal mode observations. We require that our models satisfy realistic physical constraints. We create models of the outer core with and without a distinct uppermost layer and assess the importance of this region.Our normal mode-derived models are compared with observations of body waves which travel through the outer core. In particular, we consider SmKS waves which are especially sensitive to the uppermost outer core and are therefore an important way to understand the robustness of our models.

The harzburgites-lherzolite cycle: depletion and refertilization processes

NASA Astrophysics Data System (ADS)

Dijkstra, A. H.

2011-12-01

Lherzolites or clinopyroxene-rich harzburgites sampled at the ocean floor are now generally interpreted as refractory harzburgites refertilized by melt-rock reaction or melt impregnation at the spreading center, rather than as relatively undepleted bulk upper mantle. The key evidence for a melt refertilization origin is often textural. Critically, the refertilization can mask the underlying very refractory character: oceanic peridotites prior to melt refertilization at the ridge are often too refractory to be simple mantle residues of bulk upper mantle that was melted at the ridge. This suggests that the upper mantle contains large domains that record prior melting histories. This is supported by ancient rhenium-depletion ages that are common in oceanic peridotites. In this presentation, I will discuss some key examples (e.g., Macquarie Island [1], Pindos, Totalp, Lanzarote) of refertilized oceanic peridotites, which all have recorded previous, ancient depletions. I will show the textural and geochemical evidence for melt refertilization. It has often been assumed that melt refertilization occurs by interaction with mantle melts. However, there is now evidence for melt refertilization through a reaction with eclogite-derived melts, probably at the base of the melting column underneath the ridge system. These eclogitic mantle heterogeneities themselves do not normally survive the melting underneath the spreading center, but their isotopic signature can be recognized in the reacted peridotites. In summary, we have moved away from the idea that oceanic mantle rocks are simple melting residues of homogeneous bulk upper mantle. The picture that emerges is a rich and complex one, suggesting that oceanic mantle rocks record dynamic histories of melting and refertilization. In particular, the melting event in refertilized peridotites can be much older than the age of the ridge system at which they are sampled. Many oceanic peridotites contain evidence for a Mesoproterozoic melting event of perhaps global significance. Regardless of the nature of these melting events, it is now clear that in their complex overprinting history, oceanic peridotites more and more resemble polygenetic metamorphic rocks.

Ferric Iron Production in Magma Oceans and Evolution of Mantle Oxidation State

NASA Astrophysics Data System (ADS)

Schaefer, L.; Elkins-Tanton, L. T.; Pahlevan, K.

2018-05-01

Self-oxidation of the magma ocean by ferric iron production at high pressure may explain the mantle oxidation state of the Earth. Partitioning during fractional crystallization can further increase the mantle oxygen fugacity during solidification.

Geophysical study of the structure and processes of the continental convergence zones: Alpine-Himalayan Belt

NASA Technical Reports Server (NTRS)

Toksoez, M. N.

1981-01-01

The seismic wave velocity structure in the crust and upper mantle region beneath the Tibetan plateau was studied in detail. Also, a preliminary study of the uppermost mantle P wave velocity beneath Iran and Turkey was carried out, and the results are compared with those for the Tibetan plateau. These two studies compose the bulk of the efforts on the observational aspects of continental collision zones in addition to satellite derived data. On the theoretical aspects the thermal evolution of converging plate boundaries was explored using a finite difference scheme.

Crustal and Uppermost Mantle Structure beneath the Western United States from USArray Regional Phase Analysis

NASA Astrophysics Data System (ADS)

Buehler, Janine Sylvia

The aim of this dissertation is to improve our understanding of the crust and uppermost mantle structure in the western United States, profiting from the wealth of regional phase data recorded at USArray stations. USArray, a transportable seismic array of ˜400 seismometers, has greatly increased seismic data coverage across the United States in the past few years, and allows imaging of the lithosphere of the North American continent with better resolution and new methods. The regional phases are often challenging to analyze, especially in a tectonically-active region like the western United States, because of their sensitivities to the heterogeneities of the crust and uppermost mantle. However, knowledge of the seismic structure of the lithosphere is not only essential in order to accurately image the velocity structure at greater depths, but also for constraining geodynamic models that reconstruct the tectonic evolution of the continent, and hence the information that is carried by the regional phases is very valuable. The data set used in this study consists mostly of the regional seismic phases Pn and Sn, which propagate horizontally along the Moho in the mantle lid and constrain the seismic velocity structure at a confined depth. We applied traditional tomographic methods that profit from the improved ray coverage through USArray, but also employed array-based techniques that take advantage of the regular station spacing of the transportable array and don't depend on regularization. In addition, we used stacking methods to image the propagation efficiency of the Sn phase, which is often highly attenuated in tectonically active regions, on a regional scale. The results complement other seismic studies that average over greater depth intervals, such as surface- and body-wave tomographies and anisotropy analysis from shear-wave splitting, to provide information on temperature, composition, and tectonic processes at depth. Comparisons between Pn azimuthal anisotropy and fast polarization direction from shear wave splitting suggest significant vertical changes in anisotropy in several regions of the upper mantle beneath the western United States. Sn can in theory further constrain the nature of anisotropy in the mantle lid. However, we have so far been unable to resolve shear-wave splitting directly in the Sn waveforms as the phase is often attenuated and difficult to detect. Still, we obtained evidence for Sn propagation in several regions of the western United States such as the central Great Basin and the northeastern part of the Colorado Plateau. We found that there are enough quality Sn picks for joint Pn-Sn tomography and identified prominent Vp/Vs anomalies, such as large high Vp/Vs regions --- typically associated with partial melt --- below the Snake River Plain and the Colorado Plateau.

Deformation, Fluid Flow and Mantle Serpentinization at Oceanic Transform Faults

NASA Astrophysics Data System (ADS)

Rupke, L.; Hasenclever, J.

2017-12-01

Oceanic transform faults (OTF) and fracture zones have long been hypothesized to be sites of enhanced fluid flow and biogeochemical exchange. In this context, the serpentine forming interaction between seawater and cold lithospheric mantle rocks is particularly interesting. The transformation of peridotite to serpentinite not only leads to hydration of oceanic plates and is thereby an important agent of the geological water cycle, it is also a mechanism of abiotic hydrogen and methane formation, which can support archeal and bacterial communities at the seafloor. Inferring the likely amount of mantle undergoing serpentinization reactions therefore allows estimating the amount of biomass that may be autotrophically produced at and around oceanic transform faults and mid-ocean ridges Here we present results of 3-D geodynamic model simulations that explore the interrelations between deformation, fluid flow, and mantle serpentinization at oceanic transform faults. We investigate how slip rate and fault offset affect the predicted patterns of mantle serpentinization around oceanic transform faults. Global rates of mantle serpentinization and associated H2 production are calculated by integrating the modeling results with plate boundary data. The global additional OTF-related production of H2 is found to be between 6.1 and 10.7 x 1011 mol per year, which is comparable to the predicted background mid-ocean ridge rate of 4.1 - 15.0 x 1011 mol H2/yr. This points to oceanic transform faults as potential sites of intense fluid-rock interaction, where chemosynthetic life could be sustained by serpentinization reactions.

A nearly water-saturated mantle transition zone inferred from mineral viscosity.

PubMed

Fei, Hongzhan; Yamazaki, Daisuke; Sakurai, Moe; Miyajima, Nobuyoshi; Ohfuji, Hiroaki; Katsura, Tomoo; Yamamoto, Takafumi

2017-06-01

An open question for solid-earth scientists is the amount of water in Earth's interior. The uppermost mantle and lower mantle contain little water because their dominant minerals, olivine and bridgmanite, have limited water storage capacity. In contrast, the mantle transition zone (MTZ) at a depth of 410 to 660 km is considered to be a potential water reservoir because its dominant minerals, wadsleyite and ringwoodite, can contain large amounts of water [up to 3 weight % (wt %)]. However, the actual amount of water in the MTZ is unknown. Given that water incorporated into mantle minerals can lower their viscosity, we evaluate the water content of the MTZ by measuring dislocation mobility, a property that is inversely proportional to viscosity, as a function of temperature and water content in ringwoodite and bridgmanite. We find that dislocation mobility in bridgmanite is faster by two orders of magnitude than in anhydrous ringwoodite but 1.5 orders of magnitude slower than in water-saturated ringwoodite. To fit the observed mantle viscosity profiles, ringwoodite in the MTZ should contain 1 to 2 wt % water. The MTZ should thus be nearly water-saturated globally.

Big mantle wedge, anisotropy, slabs and earthquakes beneath the Japan Sea

NASA Astrophysics Data System (ADS)

Zhao, Dapeng

2017-09-01

The Japan Sea is a part of the western Pacific trench-arc-backarc system and has a complex bathymetry and intense seismic activities in the crust and upper mantle. Local seismic tomography revealed strong lateral heterogeneities in the crust and uppermost mantle beneath the eastern margin of the Japan Sea, which was determined using P and S wave arrival times of suboceanic earthquakes relocated precisely with sP depth phases. Ambient-noise tomography revealed a thin crust and a thin lithosphere beneath the Japan Sea and significant low-velocity (low-V) anomalies in the shallow mantle beneath the western and eastern margins of the Japan Sea. Observations with ocean-bottom seismometers and electromagnetometers revealed low-V and high-conductivity anomalies at depths of 200-300 km in the big mantle wedge (BMW) above the subducting Pacific slab, and the anomalies are connected with the low-V zone in the normal mantle wedge beneath NE Japan, suggesting that both shallow and deep slab dehydrations occur and contribute to the arc and back-arc magmatism. The Pacific slab has a simple geometry beneath the Japan Sea, and earthquakes occur actively in the slab down to a depth of ∼600 km beneath the NE Asian margin. Teleseismic P and S wave tomography has revealed that the Philippine Sea plate has subducted aseismically down to the mantle transition zone (MTZ, 410-660 km) depths beneath the southern Japan Sea and the Tsushima Strait, and a slab window is revealed within the aseismic Philippine Sea slab. Seismic anisotropy tomography revealed a NW-SE fast-velocity direction in the BMW, which reflects corner flows induced by the fast deep subduction of the Pacific slab. Large deep earthquakes (M > 7.0; depth > 500 km) occur frequently beneath the Japan Sea western margin, which may be related to the formation of the Changbai and Ulleung intraplate volcanoes. A metastable olivine wedge is revealed within the cold core of the Pacific slab at the MTZ depth, which may be related to the deep seismicity. However, many of these results are still preliminary, due to the lack of seismic stations in the Japan Sea. The key to resolving these critical geoscientific issues is seismic instrumentation in the Japan Sea, for which international cooperation of geoscience communities in the East Asian countries is necessary.

Effects of Fertile Mantle Compositional Variation and Spreading Rate Variation on the Working of Global Ocean Ridges

NASA Astrophysics Data System (ADS)

Niu, Y.; O'Hara, M. J.

2014-12-01

Mantle temperature variation, plate spreading rate variation and mantle compositional variation have been considered to be the three fundamental variables that govern the working of global ocean ridges [1]. An analysis demonstrates that mantle compositional variation exerts the primary control on ocean ridge processes; it determines (1) variation in both composition and mode of mantle mineralogy, (2) variation of mantle density, (3) variation of ridge axial depth, (4) source-inherited MORB compositional variation, (4) density-controlled variation in the amplitude of mantle upwelling, (5) apparent variation in the extent of melting, and (6) the correlated variation of MORB chemistry with ridge axial depth [2]. The above interpretations are reinforced by the updated MORB database [3]. The new database also confirms spreading rate control on the extent of melting as shown previously [4]. Mantle temperature variation could play a part, but its overstated role [3,5] results from a basic error (1) in treating ridge axial depth variation as evidence of mantle temperature variation by ignoring the intrinsic control of mantle composition, (2) in treating "mantle plume" influenced ridges (e.g., Iceland) as normal ridges of plate spreading origin, and (3) in treating low Vs at greater depths (> 300 km vs. < 200 km beneath ridges) beneath these "mantle plume" influenced ridges as evidence for hot ridge mantle. In order to understand the working of global ocean ridges, we must avoid plume-influenced ridges (e.g., in the vicinity of Iceland) and remove/average out data from such ridges. As a result, the correlations (e.g., between ridge axial depth, mantle low Vs anomaly, and some geochemical parameters) required for the interpretation of mantle temperature control all disappear. There is thus no evidence for large mantle temperature variation away from ridges influenced by "mantle plumes". References: [1] Niu et al., 2001, Earth Planet Sci. Lett., 186, 383-399; [2] Niu & O'Hara, 2008, J. Petrol., 49, 633-664; [3] Gale et al., 2014, J. Petrol, 55, 1051-1082; [4] Niu & Hékinian, 1997, Nature, 385, 326-329; [5] Dalton et al., 2014, Science, 334, 80-83; [6]Niu & Hékinian, 2004, In Oceanic Hotspots, Springer-Verlag, 285-307.

Complex Anisotropic Structure of the Mantle Wedge Beneath Kamchatka Volcanoes

NASA Astrophysics Data System (ADS)

Levin, V.; Park, J.; Gordeev, E.; Droznin, D.

2002-12-01

A wedge of mantle material above the subducting lithospheric plate at a convergent margin is among the most dynamic environments of the Earth's interior. Deformation and transport of solid and volatile phases within this region control the fundamental process of elemental exchange between the surficial layers and the interior of the planet. A helpful property in the study of material deformation and transport within the upper mantle is seismic anisotropy, which may reflect both microscopic effects of preferentialy aligned crystals of olivine and orthopyroxene and macroscopic effects of systematic cracks, melt lenses, layering etc. Through the mapping of anisotropic properties within the mantle wedge we can establish patterns of deformation. Volatile content affects olivine alignment, so regions of anomalous volatile content may be evident. Indicators of seismic anisotropy commonly employed in upper mantle studies include shear wave birefringence and mode-conversion between compressional and shear body waves. When combined together, these techniques offer complementary constraints on the location and intensity of anisotropic properties. The eastern coast of southern Kamchatka overlies a vigorous convergent margin where the Pacific plate descends at a rate of almost 80 mm/yr towards the northwest. We extracted seismic anisotropy indicators from two data sets sensitive to the anisotropic properties of the uppermost mantle. Firstly, we evaluated teleseismic receiver functions for a number of sites, and found ample evidence for anisotropicaly-influenced P-to-S mode conversion. Secondly, we measured splitting in S waves of earthquakes with sources within the downgoing slab. The first set of observations provides constraints on the depth ranges where strong changes in anisotropic properties take place. The local splitting data provides constraints on the cumulative strength of anisotropic properties along specific pathways through the mantle wedge and possibly parts of the slab. To explain the vertical stratification of anisotropy implied from receiver functions, and the strong lateral dependence of shear-wave splitting observations, we cannot rely on simple models of mantle wedge behaviour e.g., olivine-crystal alignment through subduction-driven corner flow. Diverse mechanisms can contribute to the observed pattern of anisotropic properties, with volatiles likely being a key influence. For instance, we find evidence in favor of a slow-symmetry-axis anisotropy within the uppermost 10-20 km of the mantle wedge, implying either excessive hydration of the mantle or else a presence of systematically aligned volatile-filled cracks or lenses. Also, shear-wave splitting is weak beneath the Avachinsky-Koryaksky volcanic center, suggesting either vertical flow or the influence of volatiles and/or thermally-enhanced diffusion creep.

Extracting Rayleigh wave dispersion from ambient noise across the Indian Ocean

NASA Astrophysics Data System (ADS)

Ma, Z.; Dalton, C. A.

2016-12-01

Rayleigh wave dispersion extracted from ambient seismic noise has been widely used to image crustal and uppermost mantle structure. Applications of this approach in continental settings are abundant, but there have been relatively few studies within ocean basins. In this presentation, we will first demonstrate the feasibility of extracting high quality Rayleigh wave dispersion information from ambient noise across the entire Indian Ocean basin. Phase arrival times measured from ambient noise are largely consistent with the ones predicted from 2-D phase velocity maps that were determined from earthquake data alone. Secondly, we show that adding dispersion information extracted from ambient noise to existing earthquake data can indeed improve the resolution of phase velocity maps by about 20% in the western Indian Ocean region where the station distribution is the densest. High quality Rayleigh wave dispersion information can be obtained from stacking waveforms over less than two years at land stations and less than four years at island stations. After removing the age dependent average velocities, the 2-D phase velocity maps show slow anomalies associated with the Seychelles-Mascarene plateau. Forward modeling suggests that the crust is about 15-25 km thick in this area, which agrees with previous estimates obtained from gravity data. We also observe that the slow anomaly related to the Central Indian Ridge is asymmetric. The center of this slow anomaly lies to the west side of ridge, which is opposite to the ridge migration direction. This asymmetry probably reflects the interactions between the ridge and nearby hotspots.

Seismically imaging the Afar plume

NASA Astrophysics Data System (ADS)

Hammond, J. O.; Kendall, J. M.; Bastow, I. D.; Stuart, G. W.; Keir, D.; Ayele, A.; Ogubazghi, G.; Ebinger, C. J.; Belachew, M.

2011-12-01

Plume related flood basalt volcanism in Ethiopia has long been cited to have instigated continental breakup in northeast Africa. However, to date seismic images of the mantle beneath the region have not produced conclusive evidence of a plume-like structure. As a result the nature and even existence of a plume in the region and its role in rift initiation and continental rupture are debated. Previous seismic studies using regional deployments of sensors in East-Africa show that low seismic velocities underlie northeast Africa, but their resolution is limited to the top 200-300km of the Earth. Thus, the connection between the low velocities in the uppermost mantle and those imaged in global studies in the lower mantle is unclear. We have combined new data from Afar, Ethiopia with 6 other regional experiments and global network stations across Ethiopia, Eritrea, Djibouti and Yemen, to produce high-resolution models of upper mantle P- and S- wave velocities to the base of the transition zone. Relative travel time tomographic inversions show that the top 100km is dominated by focussed low velocity zones, likely associated with melt in the lithosphere/uppermost asthenosphere. Below these depths a broad SW-NE oriented sheet like upwelling extends down to the top of the transition zone. Within the transition zone two focussed sharp-sided low velocity regions exist: one beneath the Western Ethiopian plateau outside the rift valley, and the other beneath the Afar depression. The nature of the transition zone anomalies suggests that small upwellings may rise from a broader low velocity plume-like feature in the lower mantle. This interpretation is supported by numerical and analogue experiments that suggest the 660km phase change and viscosity jump may impede flow from the lower to upper mantle creating a thermal boundary layer at the base of the transition zone. This allows smaller, secondary upwellings to initiate and rise to the surface. Our images of secondary upwellings suggest that there is no evidence for a plume in the classical sense (i.e. a narrow conduit). Instead, we propose that secondary upwellings rise from the base of the transition zone and connect in the upper mantle. This coupled with measurements of seismic anisotropy suggest that mantle material flows northeast towards Arabia, and may be responsible for the dramatic dynamic topography observed in northeast Africa and western Arabia.

Upper mantle beneath foothills of the western Himalaya: subducted lithospheric slab or a keel of the Indian shield?

NASA Astrophysics Data System (ADS)

Vinnik, L.; Singh, A.; Kiselev, S.; Kumar, M. Ravi

2007-12-01

The fate of the mantle lithosphere of the Indian Plate in the India-Eurasia collision zone is not well understood. Tomographic studies reveal high P velocity in the uppermost mantle to the south of the western Himalaya, and these high velocities are sometimes interpreted as an image of subducting Indian lithosphere. We suggest that these high velocities are unrelated to the ongoing subduction but correspond to a near-horizontal mantle keel of the Indian shield. In the south of the Indian shield upper-mantle velocities are anomalously low, and relatively high velocities may signify a recovery of the normal shield structure in the north. Our analysis is based on the recordings of seismograph station NIL in the foothills of the western Himalaya. The T component of the P receiver functions is weak relative to the Q component, which is indicative of a subhorizontally layered structure. Joint inversion of the P and S receiver functions favours high uppermost mantle velocities, typical of the lithosphere of Archean cratons. The arrival of the Ps converted phase from 410 km discontinuity at NIL is 2.2 s earlier than in IASP91 global model. This can be an effect of remnants of Tethys subduction in the mantle transition zone and of high velocities in the keel of the Indian shield. Joint inversion of SKS particle motions and P receiver functions reveals a change in the fast direction of seismic azimuthal anisotropy from 60° at 80-160 km depths to 150° at 160-220 km. The fast direction in the lower layer is parallel to the trend of the Himalaya. The change of deformation regimes at a depth of 160 km suggests that this is the base of the lithosphere of the Indian shield. A similar boundary was found with similar techniques in central Europe and the Tien Shan region, but the base of the lithosphere in these regions is relatively shallow, in agreement with the higher upper-mantle temperatures. The ongoing continental collision is expressed in crustal structure: the crust beneath NIL is very thick (58 +/- 2 km), and the S velocity in the intermediate and lower crust is around 4.0 km s-1. This anomalously large velocity and thickness can be explained by scraping off the lower crust, when the Indian lithosphere underthrusts the Himalaya.

The nature of the crust under Cayman Trough from gravity

USGS Publications Warehouse

ten Brink, Uri S.; Coleman, D.F.; Dillon, William P.

2002-01-01

Considerable crustal thickness variations are inferred along Cayman Trough, a slow-spreading ocean basin in the Caribbean Sea, from modeling of the gravity field. The crust to a distance of 50 km from the spreading center is only 2–3 km thick in agreement with dredge and dive results. Crustal thickness increases to ∼5.5 km at distances between 100 and 430 km west of the spreading center and to 3.5–6 km at distances between 60 and 370 km east of the spreading center. The increase in thickness is interpreted to represent serpentinization of the uppermost mantle lithosphere, rather than a true increase in the volume of accreted ocean crust. Serpentinized peridotite rocks have indeed been dredged from the base of escarpments of oceanic crust rocks in Cayman Trough. Laboratory-measured density and P-wave speed of peridotite with 40–50% serpentine are similar to the observed speed in published refraction results and to the inferred density from the model. Crustal thickness gradually increases to 7–8 km at the far ends of the trough partially in areas where sea floor magnetic anomalies were identified. Basement depth becomes gradually shallower starting 250 km west of the rise and 340 km east of the rise, in contrast to the predicted trend of increasing depth to basement from cooling models of the oceanic lithosphere. The gradual increase in apparent crustal thickness and the shallowing trend of basement depth are interpreted to indicate that the deep distal parts of Cayman Trough are underlain by highly attenuated crust, not by a continuously accreted oceanic crust.

Towards understanding carbon recycling at subduction zones - lessons from Central America

NASA Astrophysics Data System (ADS)

Hilton, D. R.; Barry, P. H.; Fischer, T. P.

2010-12-01

Subduction zones provide the essential pathways for input of carbon from Earth’s external reservoirs (crust, sediments, oceans) to the mantle. However, carbon input to the deep interior is interrupted by outputs via the fore-arc, volcanic front, and back-arc regions. Coupled CO2 and He isotope data for geothermal fluids from throughout Central American (CA) are used to derive estimates of the output carbon flux for comparison with inputs estimated for the subducting Cocos Plate. The carbon flux carried by the incoming sediments is ~1.6 × 109 gCkm-1yr-1[1], as is the ratio of input carbon derived from pelagic limestone (L) and organic sediment (S), i.e., L/S ~10.7. Additionally, the upper 7 km of oceanic (crustal) basement supplies ~9.1 × 108 gCkm-1yr-1[2]: this flux is dominated by L-derived CO2. In terms of output, measured carbon concentrations coupled with flow rates for submarine cold seeps sites at the Costa Rica outer forearc yield CO2 and CH4 fluxes of ~ 6.1 × 103 and 8.0 × 105 (gCkm-1yr-1), respectively [3]. On the Nicoya Peninsula, the Costa Rica Pacific coastline (including the Oso Peninsula) and the Talamanca Mountain Range, coupled CO2-He studies allow recognition of a deep input (3He/4He up to 4RA) and resolution of CO2 into L- and S-components. There is an increase in the L/S ratio arc-ward with the lowest values lying close to diatomaceous ooze in the uppermost sequence of subducting sediment package. This observation is consistent with under-plating and removal of the uppermost organic-rich sediment from deeper subduction. As the input carbon fluxes of the individual sedimentary layers are well constrained [1], we can limit the potential steady-state flux of carbon loss at the subaerial fore-arc to ~ 6 × 107 gCkm-1yr-1, equivalent to ~88% of the input flux of the diatomaceous ooze, or < 4% of the total incoming sedimentary carbon. The greatest loss of slab-derived carbon occurs at the volcanic front. Estimates of the output CO2 flux along the CA front - 2-5 (× 108 gCkm-1yr-1) [4-5] together with identification of a slab origin (~90%) of the CO2, gives output estimates between 12% (Costa Rica) and 29% (El Salvador) of the sedimentary input [6]. The low L/S ratio found along the entire strike of the volcanic front precludes a significant C-contribution from oceanic basement of the subducting slab. Finally, arc-like L/S ratios behind the volcanic front in Honduras [6] indicates the back-arc inventory is composed of either entrained or ancient CO2 but not slab carbon released beyond the region of arc magma generation. Thus, at the CA subduction zone, significant carbon influx to the mantle can occur due to limited fore-arc and back-arc losses and modest C-outputs via the volcanic front. These observations are compared with other subduction zones where sediment lithologies, thermal conditions and water budgets differ, to address the question of understanding intrinsic and extrinsic controls on the mass balance of the mantle carbon reservoir. [1] Li and Bebout, JGR, 2005; [2] Hilton et al., Rev. Min. Geochem., 2002; [3] Furi et al., G-cubed, 2010; [4] Rodriguez et al., JVGR, 2004; [5] Zimmer et al., G-cubed, 2004; [6] De Leeuw et al., EPSL, 2007.

Synthetic Pn and Sn phases and the frequency dependence of Q of oceanic lithosphere

NASA Astrophysics Data System (ADS)

Sereno, Thomas J., Jr.; Orcutt, John A.

1987-04-01

The oceanic lithosphere is an extremely efficient waveguide for high-frequency seismic energy. In particular, the propagation of the regional to teleseismic oceanic Pn and Sn phases is largely controlled by properties of the oceanic plates. The shallow velocity gradient in the sub-Moho lithosphere results in a nearly linear travel time curve for these oceanic phases and an onset velocity near the material velocity of the uppermost mantle. The confinement of Pn/Sn to the lithosphere imposes a constraint on the maximum range that a normally refracted wave can be observed. The rapid disappearance of Sn and the discontinuous drop in Pn/Sn group velocity beyond a critical distance, dependent upon the local thickness of the lithosphere, are interpreted as a shadowing effect of the low Q asthenosphere. Wave number integration was used to compute complete synthetic seismograms for a model of oceanic lithosphere. The results were compared to data collected during the 1983 Ngendei Seismic Experiment in the southwest Pacific. The Pn/Sn coda is successfully modeled as a sum of leaky organ-pipe modes in the sediment layer and oceanic water column. While scattering is present to some degree, it is not required to explain the long duration and complicated nature of the Pn/Sn wave trains. The presence of extremely high frequencies in Pn/Sn phases and the greater efficiency of Sn than Pn propagation are interpreted in terms of an absorption band rheology. A shorter high-frequency relaxation time for P waves than for S waves results in a rheology with the property that Qα > Qβ at low frequency while Qβ > Qα at high frequency, consistent with the teleseismic Pn/Sn observations. The absorption band model is to viewed as only an approximation to the true frequency dependence of Q in the oceanic lithosphere for which analytic expressions for the material dispersion have been developed.

Crust And Upper Mantle Structure Of The Bengal Basin And Bay Of Bengal From Surface Wave Group Velocity Dispersion Studies

NASA Astrophysics Data System (ADS)

Dhali, K. K.; Majhi, S.; Mitra, S.; Priestley, K.

2007-12-01

Fundamental mode Rayleigh and Love wave group velocity dispersion for paths crossing the Bay of Bengal have been calculated for earthquakes in the Indo-Burman arc and the Andaman-Sumatra subduction zone recorded at seismographs in the eastern part of Peninsula India and Sri Lanka. The ray-path coverage in this study provides a better spatial sampling than any previous studies of the region. The individual dispersion curves range from 12 to 70~s and have been clustered in four spatial groups to form average dispersion curves representative of the Bengal basin, northern, central and southern Bay of Bengal. These average dispersion curves for Rayleigh and Love waves are jointly inverted to obtain shear wave velocity structure of the lithosphere. The higher frequencies/shorter periods (12--30~s) used in the inversion constrains the sediment shear wave speed and thickness while the longer periods provide information of the upper mantle structure. The results show a remarkable increase in the sediments thickness along the Bengal Fan from south to north ranging from 6 km, around the southern tip of India, to 23 km beneath the Bengal basin. The shear wave velocity models reveal a sediment saturation beyond 7-10 km of burial leading to metamorphism and eventual increase in velocity to continent like material with depth. The average crustal thickness (loose sediments overlying consolidated sediments followed by metasediments and oceanic crust) is anomalously continental (~20-36 km) rather than being simply oceanic crust overlain by sediments. The average shear wave velocity is about 3.5-3.8 km/s which is more representative of continental crusts. Finally the low velocity zone in the uppermost mantle is possibly an effect of the expected increase in temperature due to blanketing of the fan sediments over the Bay of Bengal crust. The misfits to parts of the dispersion data using a 1D isotropic model provides an indication of the presence of polarization anisotropy in the lithosphere and sets a good starting point for modeling the anisotropic structure.

Abundant carbon in the mantle beneath Hawai`i

NASA Astrophysics Data System (ADS)

Anderson, Kyle R.; Poland, Michael P.

2017-09-01

Estimates of carbon concentrations in Earth’s mantle vary over more than an order of magnitude, hindering our ability to understand mantle structure and mineralogy, partial melting, and the carbon cycle. CO2 concentrations in mantle-derived magmas supplying hotspot ocean island volcanoes yield our most direct constraints on mantle carbon, but are extensively modified by degassing during ascent. Here we show that undegassed magmatic and mantle carbon concentrations may be estimated in a Bayesian framework using diverse geologic information at an ocean island volcano. Our CO2 concentration estimates do not rely upon complex degassing models, geochemical tracer elements, assumed magma supply rates, or rare undegassed rock samples. Rather, we couple volcanic CO2 emission rates with probabilistic magma supply rates, which are obtained indirectly from magma storage and eruption rates. We estimate that the CO2 content of mantle-derived magma supplying Hawai`i’s active volcanoes is 0.97-0.19+0.25 wt%--roughly 40% higher than previously believed--and is supplied from a mantle source region with a carbon concentration of 263-62+81 ppm. Our results suggest that mantle plumes and ocean island basalts are carbon-rich. Our data also shed light on helium isotope abundances, CO2/Nb ratios, and may imply higher CO2 emission rates from ocean island volcanoes.

Processes of lithosphere evolution: New evidence on the structure of the continental crust and uppermost mantle

USGS Publications Warehouse

Artemieva, I.M.; Mooney, W.D.; Perchuc, E.; Thybo, H.

2002-01-01

We discuss the structure of the continental lithosphere, its physical properties, and the mechanisms that formed and modified it since the early Archean. The structure of the upper mantle and the crust is derived primarily from global and regional seismic tomography studies of Eurasia and from global and regional data on seismic anisotropy. These data as documented in the papers of this special issue of Tectonophysics are used to illustrate the role of different tectonic processes in the lithospheric evolution since Archean to present. These include, but are not limited to, cratonization, terrane accretion and collision, continental rifting (both passive and active), subduction, and lithospheric basal erosion due to a relative motion of cratonic keels and the convective mantle. ?? 2002 Elsevier Science B.V. All rights reserved.

Initiation and evolution of the Oligo-Miocene rift basins of southwestern Europe: Contribution of deep seismic reflection profiling

NASA Astrophysics Data System (ADS)

Bois, C.

1993-11-01

Southwestern European Oligo-Miocene rift basins have recently been investigated by deep seismic reflection profiling. The study of these data, together with other geophysical and geological data, shows that the rifts, which run from the Rhinegraben to the western Mediterranean, do not form a single clearcut system. The N-trending rifts (Rhinegraben, Bresse and Limagne) were developed on a cold and rigid lithosphere affected by the Alpine collision. The NE-trending rifts (southeastern France, Gulf of Lions and Valencia Trough) were formed slightly later in a backarc basin associated with an active segment of the European-Iberian plate that was heated, affected by widespread calcalkaline volcanism and probably weakened. All the southwestern European rifts and basins together may, however, be related to a common heritage represented by the boundary between the European-Iberian and African-Apulian plates that was created in the Jurassic with the initiation of the Tethys Ocean. The present features of the southwestern European Oligo-Miocène rift basins may be explained by a combination of three geodynamic mechanisms: mechanical stretching of the lithosphere, active mantle uplifting, and subordinate lithospheric flexuring. All the rifts were probably initiated by passive stretching. A systematic discrepancy between stretching derived from fault analysis and attenuation of the crust has been observed in all the rifts. This suggests that these rifts were subsequently reworked by one or several active mantle upwelling events associated with late shoulder uplift, asthenosphere upwelling and anomalous P-wave velocities in the lowermost crust and the uppermost mantle. Crustal attenuation may have been achieved by mantle intrusion, metamorphism of the deep crust and/or its delamination. Some of the rifts were affected by lithospheric flexuring. Combinations, in various proportions, of a small number of geodynamic mechanisms probably controlled many basins in the world. This explains the unique characteristics of each basin, difficulties in basin classification and the frequent failure of single-mechanism models to explain the geological observations.

Improved Nazca slab structure from teleseismic P-wave tomography along the Andean margin

NASA Astrophysics Data System (ADS)

Portner, D. E.; Beck, S. L.; Scire, A. C.; Zandt, G.

2017-12-01

South America marks the longest continuous ocean-continent subduction zone. As such, there is significant along-strike variability in the subducting Nazca slab structure and the tectonics of the South American margin. Most notably two gaps in the otherwise continuous volcanic arc are correlated with regions of flat slab subduction, indicating that the structure of the Nazca slab plays a controlling role in South American tectonics. Traditionally in subduction zones, our knowledge of slab structure is defined by Wadati-Benioff zone earthquakes. While this method allows for the determination of large-scale variations in Nazca slab structure such as regions of flat slab subduction, a scarcity of intermediate-depth earthquakes hinders our ability to observe the smaller-scale structural variations in the slab that may be critical to our understanding of the geologic record. We use an updated, larger dataset for finite-frequency teleseismic P-wave tomography including relative arrival times from >700 seismic stations along the Andean margin to image the detailed Nazca slab structure throughout the upper mantle and uppermost lower mantle between latitudes 5°S and 45°S. Our results show prominent variations in slab character along the margin. Slab dip varies significantly, from sub-vertical inboard of the Peruvian flat slab segment to 30° dip south of the Pampean flat slab, while the slab's velocity anomaly amplitude changes dramatically near the Pampean flat slab region. High slab velocities north of the Pampean region relative to the south indicate variable slab thermal structures that correspond roughly with the locations of deep (>500 km depth) earthquakes that also occur exclusively north of the Pampean region. Additionally, a wider regional footprint increases our sampling of the upper-lower mantle boundary, improving constraints on the slab's interaction with the 660 km discontinuity along strike. We see that the Nazca slab appears to penetrate into the lower mantle along the majority of the margin.

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

NASA Astrophysics Data System (ADS)

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

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

Long wavelength magnetic anomalies over continental rifts in cratonic region

NASA Astrophysics Data System (ADS)

Friedman, S. A.; Persaud, P.; Ferre, E. C.; Martín-Hernández, F.; Feinberg, J. M.

2017-12-01

New collections of unaltered mantle xenoliths shed light on potential upper mantle contributions to long wavelength magnetic anomalies (LWMA) in continental rifts in cratonic / shield areas. The new material originates from the East African Rift (Tanzania), the Rio Grande Rift (U.S.A.), the Rhine Rift (Germany), and the West Antarctic Rift (Antarctica). The xenoliths sample the uppermost (<80 km depth) lithospheric mantle in these regions in the spinel-peridotite and plagioclase-peridotite stability fields. The most common lithology by far (95% of samples) is a spinel-lherzolite indicating relatively low oxygen fugacities (FMQ -1). Chrome spinel in these peridotites is non-magnetic (Al + Mg > 0.2 or Fe < 0.3) and primary magnetite (Fe3O4) occurs only in trace amounts, typically yielding low natural remanent magnetizations (NRM < 10-2 A/m). The low Koenigsberger ratios (Qn < 1) of these materials, combined with high geotherms (>60ºC/km) that are characteristic of rifted regions preclude any contribution to LWMA at depths >10 km. Hence, only upper basalts and hypovolcanic mafic sills would constitute potential magnetic sources. In contrast, the margins of these rifted regions consist of refractory cratonic domains, often characterized by oxidized sublithospheric mantle that host significant concentrations of primary magnetite. The higher NRMs of these peridotites (up to 15 A/m, Qn > 2.5) combined with much lower geotherms (as low as 15ºC/km) allows for a 5 to 10 km layer of uppermost mantle to potentially contribute to LWMA. Assuming that Qn values in rift margins are also <1, the new data presented here suggests that relatively young rifts would display a central negative magnetic anomaly surrounded by two broad positive anomalies. The latitude of the rift is predicted to exert a primary control on the magnitude of such anomalies, while the steepness of the magnetic gradient across the rift would primarily reflect thermal equilibration over time.

Imaging Canary Island hotspot material beneath the lithosphere of Morocco and southern Spain

NASA Astrophysics Data System (ADS)

Miller, Meghan S.; O'Driscoll, Leland J.; Butcher, Amber J.; Thomas, Christine

2015-12-01

The westernmost Mediterranean has developed into its present day tectonic configuration as a result of complex interactions between late stage subduction of the Neo-Tethys Ocean, continental collision of Africa and Eurasia, and the Canary Island mantle plume. This study utilizes S receiver functions (SRFs) from over 360 broadband seismic stations to seismically image the lithosphere and uppermost mantle from southern Spain through Morocco and the Canary Islands. The lithospheric thickness ranges from ∼65 km beneath the Atlas Mountains and the active volcanic islands to over ∼210 km beneath the cratonic lithosphere in southern Morocco. The common conversion point (CCP) volume of the SRFs indicates that thinned lithosphere extends from beneath the Canary Islands offshore southwestern Morocco, to beneath the continental lithosphere of the Atlas Mountains, and then thickens abruptly at the West African craton. Beneath thin lithosphere between the Canary hot spot and southern Spain, including below the Atlas Mountains and the Alboran Sea, there are distinct pockets of low velocity material, as inferred from high amplitude positive, sub-lithospheric conversions in the SRFs. These regions of low seismic velocity at the base of the lithosphere extend beneath the areas of Pliocene-Quaternary magmatism, which has been linked to a Canary hotspot source via geochemical signatures. However, we find that this volume of low velocity material is discontinuous along strike and occurs only in areas of recent volcanism and where asthenospheric mantle flow is identified with shear wave splitting analyses. We propose that the low velocity structure beneath the lithosphere is material flowing sub-horizontally northeastwards beneath Morocco from the tilted Canary Island plume, and the small, localized volcanoes are the result of small-scale upwellings from this material.

Seismic and Petrological Constraints on Deep Crustal Evolution in North America: Where and What are 7.x Layers?

NASA Astrophysics Data System (ADS)

Mahan, K. H.; Schulte-Pelkum, V.; Shen, W.; Ritzwoller, M. H.

2012-12-01

Continental crust worldwide has been found to have areas with a lowermost layer characterized by unusually high seismic P velocities of over 7 km/s, often called 7.x layers. Such layers are commonly ascribed to underplating - in some cases by underthrusting, but in most cases by magmatic processes. In North America, high-velocity lower crust underlies upper crust of Archean, Proterozoic, and younger ages. Its presence reflects the tectonic and magmatic processes associated with continental rifting, collision, subduction, and other evolutionary (e.g. thermal) trends, and its occurrence also provides clues on the nature of the underlying mantle. Detection of a lower crustal high-velocity layer stems mostly from seismic refraction and wide-angle reflection experiments, and information on its geographical extent is very spotty. Similarly sparse are age determinations and knowledge of the tectonic processes responsible for construction of these layers. Despite glimpses of 7.x layers on many profiles across the continental U.S. and Canada, there is no systematic geographical and age information on this fundamental process of crustal growth, and many of the existing observations contradict current hypotheses on underplating. We compare compositional and physical property data of lower crustal and uppermost mantle xenoliths from Montana, Wyoming, and other localities with maps of lower crustal and uppermost mantle seismic velocities obtained from joint inversions of receiver functions with surface waves, and to mapped distinct high-velocity lower crustal layers in receiver functions in areas covered by the EarthScope Transportable Array. Xenolith observations from Montana indicate that portions of metasomatized uppermost mantle exist in that area that may be difficult to distinguish from mafic lower crust based on seismic velocities alone, raising the interesting question of whether a 7.x layer may be below rather than above the seismic Moho in some cases. The persistence of high-velocity, presumably strong lower crust under the Laramide-affected Wyoming craton and the Colorado Plateau suggest that crustal strength may influence surface deformation. The Rocky Mountain Front and Rio Grande rift largely separate fast lower crust to the East from slower lower crust to the West, cutting across NE-SW trends inherited from continental assembly and suggesting that the velocity distribution may be dominated by thermal effects; however, recent volcanics do not correlate well geographically with lower crustal velocity.

The origin of volatiles in the Earth's mantle

NASA Astrophysics Data System (ADS)

Hier-Majumder, Saswata; Hirschmann, Marc M.

2017-08-01

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