Seismic velocity model of the crust and uppermost mantle around the Mirnyi kimberlite field in Siberia

We present the first detailed seismic velocity models of the crust and uppermost mantle around the Mirnyi kimberlite field in Yakutia, Siberia. We have digitized vintage seismograms that were acquired in 1981 and 1983 by use of Taiga analogue seismographs along two perpendicular seismic profiles. The 370-km long, northwest striking profile I across the kimberlite pipe was covered by 41 seismographs, which recorded seismic signals from 21 chemical shots along the line, including one off-end shot. The perpendicular, 340-km long profile II across profile I ca. 30 km to the south of the Mirnyi kimberlite field was covered by 45 seismographs, which recorded seismic signals from 22 chemical shots, including four off-end shots. Each shot involved detonation of between 1.5 and 6.0 tons of TNT, distributed in individual charges of 100–200 kg in shallow water (< 2 m deep). The data is of high quality with high signal/noise ratio to the farthest offsets. We present the results from two-dimensional ray tracing, forward modelling.Both velocity models show normal cratonic structure of the ca. 45-km-thick crust with only slight undulation of the Moho. However, relatively small seismic velocity is detected to 25-km depth in a ca. 60-km wide zone around the kimberlite pipe, surrounded by elevated velocity (> 6.3 km/s) in the upper crust. The lower crust has a relatively constant velocity of 6.8–6.9 km/s. It appears relatively unaffected by the presence of the kimberlite field. Extremely large P-wave velocity (> 8.7 km/s) of the sub-Moho mantle is interpreted along profile I, except for a 70-km wide zone with a “normal” Pn velocity of 8.1 km/s below the kimberlite. Profile II mainly shows Pn velocities of 8.0–8.2 km/s, with unusually large velocity (> 8.5 km/s) in two, ca. 100-km wide zones, at its southwestern end, one zone being close to the kimberlite field. The nature of these exceptionally large, sub-Moho mantle velocities is not yet understood. The difference in velocity in the two profile directions indicates anisotropy, but the effect of unusual rock composition, e.g. from a high concentration of garnet, cannot be excluded.     

Structure of the crust and uppermost mantle of the Yanshan Belt and adjacent regions at the northeastern boundary of the North China Craton from Rayleigh Wave Dispersion Analysis

We have used Rayleigh wave dispersion analysis and inversion to produce a high resolution S-wave velocity imaging profile of the crust and uppermost mantle structure beneath the northeastern boundary regions of the North China Craton (NCC). Using waveform data from 45 broadband NCISP stations, Rayleigh wave phase velocities were measured at periods from 10 to 48 s and utilized in subsequent inversions to solve for the S-wave velocity structure from 15 km down to 120 km depth. The inverted lower crust and uppermost mantle velocities, about 3.75 km/s and 4.3 km/s on average, are low compared with the global average. The Moho was constrained in the depth range of 30–40 km, indicating a typical crustal thickness along the profile. However, a thin lithosphere of no more than 100 km was imaged under a large part of the profile, decreasing to only ~ 60 km under the Inner Mongolian Axis (IMA) where an abnormally slow anomaly was observed below 60 km depth. The overall structural features of the study region resemble those of typical continental rift zones and are probably associated with the lithospheric reactivation and tectonic extension widespread in the eastern NCC during Mesozoic–Cenozoic time. Distinctly high velocities, up to ~ 4.6 km/s, were found immediately to the south of the IMA beneath the northern Yanshan Belt (YSB), extending down to > 100-km depth. The anomalous velocities are interpreted as the cratonic lithospheric lid of the region, which may have not been affected by the Mesozoic–Cenozoic deformation process as strongly as other regions in the eastern NCC. Based on our S-wave velocity structural image and other geophysical observations, we propose a possible lithosphere–asthenosphere interaction scenario at the northeastern boundary of the NCC. We speculate that significant undulations of the base of the lithosphere, which might have resulted from the uneven Mesozoic–Cenozoic lithospheric thinning, may induce mantle flows concentrating beneath the weak IMA zone. The relatively thick lithospheric lid in the northern YSB may serve as a tectonic barrier separating the on-craton and off-craton regions into different upper mantle convection systems at the present time.     

Seismic properties of the upper mantle beneath Lanzarote (Canary Islands): Model predictions based on texture measurements by EBSD

We present a petrophysical analysis of upper mantle xenoliths, collected in the Quaternary alkali basalt fields (Series III and IV) from the island of Lanzarote. The samples consist of eight harzburgite and four dunite nodules, 5 to 15 cm in size, and exhibit a typical protogranular to porphyroclastic texture. An anomalous foliation resulting from strong recovery processes is observed in half of the specimens. The lattice preferred orientations (LPO) of olivine, orthopyroxene and clinopyroxene were measured using electron backscatter diffraction (EBSD). In most samples, olivine exhibits LPOs intermediate between the typical single crystal texture and the [100] fiber texture. Occasionally, the [010] fiber texture was also observed. Simultaneous occurrence of both types of fiber textures suggests the existence of more than one deformation regime, probably dominated by a simple shear component under low strain rate and moderate to high temperature. Orthopyroxene and clinopyroxene display a weaker but significant texture. The LPO data were used to calculate the seismic properties of the xenoliths at PT conditions obtained from geothermobarometry, and were compared to field geophysical data reported from the literature. The velocity of P-waves (7.9 km/s) obtained for a direction corresponding to the existing seismic transect is in good agreement with the most recent geophysical interpretation. Our results are consistent with a roughly W–E oriented fastest P-wave propagation direction in the uppermost mantle beneath the Canary Islands, and with the lithosphere structure proposed by previous authors involving a crust–mantle boundary at around 18 km in depth, overlaid by intermediate material between 11 and 18 km.     

Continuation of the San Andreas fault system into the upper mantle: Evidence from spinel peridotite xenoliths in the Coyote Lake basalt, central California

The Coyote Lake basalt, located near the intersection of the Hayward and Calaveras faults in central California, contains spinel peridotite xenoliths from the mantle beneath the San Andreas fault system. Six upper mantle xenoliths were studied in detail by a combination of petrologic techniques. Temperature estimates, obtained from three two-pyroxene geothermometers and the Al-in-orthopyroxene geothermometer, indicate that the xenoliths equilibrated at 970–1100 °C. A thermal model was used to estimate the corresponding depth of equilibration for these xenoliths, resulting in depths between 38 and 43 km. The lattice preferred orientation of olivine measured in five of the xenolith samples show strong point distributions of olivine crystallographic axes suggesting that fabrics formed under high-temperature conditions. Calculated seismic anisotropy values indicate an average shear wave anisotropy of 6%, higher than the anisotropy calculated from xenoliths from other tectonic environments. Using this value, the anisotropic layer responsible for fault-parallel shear wave splitting in central California is less than 100 km thick. The strong fabric preserved in the xenoliths suggests that a mantle shear zone exists below the Calaveras fault to a depth of at least 40 km, and combining xenolith petrofabrics with shear wave splitting studies helps distinguish between different models for deformation at depth beneath the San Andrea fault system.     

Nature of the mantle roots beneath the North American craton: mantle xenolith evidence from Somerset Island kimberlites

The recently discovered Nikos kimberlite on Somerset Island, in the Canadian Arctic, hosts an unusually well preserved suite of mantle xenoliths dominated by garnet–peridotite (lherzolite, harzburgite, dunite) showing coarse and porphyroclastic textures, with minor garnet–pyroxenite. The whole rock and mineral data for 54 Nikos xenoliths indicate a highly refractory underlying mantle with high olivine forsterite contents (ave. Fo=92.3) and moderate to high olivine abundances (ave. 80 wt.%). These characteristics are similar to those reported for peridotites from the Archean Kaapvaal and Siberian cratons (ave. olivine Fo=92.5), but are clearly distinct from the trend defined by oceanic peridotites and mantle xenoliths in alkaline basalts and kimberlites from post-Archean continental terranes (ave. olivine Fo=91.0). The Nikos xenoliths yield pressures and temperatures of last equilibration between 20 and 55 kb and 650 and 1300°C, and a number of the peridotite nodules appear to have equilibrated in the diamond stability field. The pressure and temperature data define a conductive paleogeotherm corresponding to a surface heat flow of 44 mW/m². Paleogeotherms based on xenolith data from the central Slave province of the Canadian craton require a lower surface heat flow (40 mW/m²) indicating a cooler geothermal regime than that beneath the Canadian Arctic. A large number of kimberlite-hosted peridotites from the Kaapvaal craton in South Africa and parts of the Siberian craton are characterized by high orthopyroxene contents (ave. Kaapvaal 32 wt.%, Siberia 20 wt.%). The calculated modal mineral assemblages for the Nikos peridotites show moderate to low contents of orthopyroxene (ave. 12 wt.%), indicating that the orthopyroxene-rich mineralogy characteristic of the Kaapvaal and Siberian cratons is not a feature of the cratonic upper mantle beneath Somerset Island.     

A composite geologic and seismic profile beneath the southern Rio Grande rift, New Mexico, based on xenolith mineralogy, temperature, and pressure

A complete understanding of the processes of crustal growth and recycling in the earth remains elusive, in part because data on rock composition at depth is scarce. Seismic velocities can provide additional information about lithospheric composition and structure, however, the relationship between velocity and rock type is not unique. The diverse xenolith suite from the Potrillo volcanic field in the southern Rio Grande rift, together with velocity models derived from reflection and refraction data in the area, offers an opportunity to place constraints on the composition of the crust and upper mantle from the surface to depths of  60 km. In this work, we calculate seismic velocities of crustal and mantle xenoliths using modal mineralogy, mineral compositions, pressure and temperature estimates, and elasticity data. The pressure, temperature, and velocity estimates from xenoliths are then combined with sonic logs and stratigraphy estimated from drill cores and surface geology to produce a geologic and velocity profile through the crust and upper mantle. Lower crustal xenoliths include garnet ± sillimanite granulite, two-pyroxene granulite, charnokite, and anorthosite. Metagabbro and amphibolite account for only a small fraction of the lower crustal xenoliths, suggesting that a basaltic underplate at the crust–mantle boundary is not present beneath the southern Rio Grande rift. Abundant mid-crustal felsic to mafic igneous xenoliths, however, suggest that plutonic rocks are common in the middle crust and were intraplated rather than underplated during the Cenozoic. Calculated velocities for garnet granulite are between  6.9 and 8.0 km/s, depending on garnet content. Granulites are strongly foliated and lineated and should be seismically anisotropic. These results suggest that velocities > 7.0 km/s and a layered structure, which are often attributed to underplated mafic rocks, can also be characteristic of alternating garnet-rich and garnet-poor metasedimentary rocks. Because the lower crust appears to be composed largely of metasedimentary granulite, which requires deep burial of upper crustal materials, we suggest the initial construction of the continental crust beneath the Potrillo volcanic field occurred by thickening of supracrustal material in the absence of large scale magmatic accretion. Mantle xenoliths include spinel lherzolite and harzburgite, dunite, and clinopyroxenite. Calculated P-wave velocities for peridotites range from 7.75 km/s to 7.89 km/s, with an average of 7.82 km/s. This velocity is in good agreement with refraction and reflection studies that report P_n velocities of 7.6–7.8 km/s throughout most of the Rio Grande rift. These calculations suggest that the low P_n velocities compared to average uppermost mantle are the result of relatively high temperatures and low pressures due to thin crust, as well as a fertile, Fe-rich, bulk upper mantle composition. Partial melt or metasomatic hydration of the mantle lithosphere are not needed to produce the observed P_n velocities.     

Monomineral universal clinopyroxene and garnet barometers for peridotitic,eclogitic and basaltic systems

New versions of the universal Jd-Di exchange clinopyroxene barometer for peridotites,pyroxenites and eclogites,and also garnet barometer for eclogites and peridotites were developed.They were checked using large experimental data sets for eciogitic(~530) and peridotitic systems(650).The precision of the universal Cpx barometer for peridotites based on Jd-Di exchange is close to Cr-Tschermalite method produced by Nimis and Taylor(2000).Cpx barometer was transformed by the substitution of major multiplier for K_D by the equations dependent from Al-Na-Fe.Obtained equation in combination with the thermometer of Nimis and Taylor(2000) allow to reconstruct position of the magma feeder systems of the alkali basaltic magma within the mantle diapirs in modern platforms like in Vitim plateau and other Southern Siberia localities and several localities worldwide showing good agreement of pressure ranges for black and green suites.These equations allow construct PTX diagrams for the kimberlite localities in Siberia and worldwide calculating simultaneously the PT parameters for different groups of mantle rocks.They give very good results for the concentrates from kimberlite lamproites and placers with mantle minerals.They are useful for PT estimates for diamond inclusions.The positions of eclogite groups in mantle sections are similar to those determined with new Gar—Cpx barometer produced by C.Beyer et al.(2015).The Fe rich eclogites commonly trace the boundary between the lower upper parts of subcontinental lithospheric mantle(SCLM) at 3-4 CPa marking pyroxenite eclogites layer.Ca-rich eclogites and especially grospydires in SCLM beneath Precambrian kimberlites occurs near pyroxenite layer but in younger mantle sections they became common in the lower parts.The diamondiferous Mg Cr-less group eclogites referring to the ancient island arc complexes are also common in the middle part of mantle sections and near 5-6 GPa.Commonly eclogites in lower apart of mantle sections are remelted and trace the high temperature convective branch.The Mg-and Fe-rich pyroxenites also show the extending in pressure trends which suggest the anatexic melting under the influence of volatiles or under the interaction with plums.     

Pressure–temperature evolution of lawsonite eclogite in Sivrihisar; Tavşanlı Zone–Turkey

The Cretaceous blueschist belt, Tavşanlı Zone, representing the subducted and exhumed northern continental margin of the Anatolide–Tauride platform is exposed in Western Anatolia. The Sivrihisar area east of Tavşanlı is made up of tectonic units consisting of i) metaclastics and conformably overlying massive marbles (coherent blueschist unit), ii) blueschist-eclogite unit, iii) marble–calcschist intercalation and iv) metaperidotite slab. The metaclastics are composed of jadeite–lawsonite–glaucophane and jadeite–glaucophane–chloritoid schists, phengite phyllites, and calcschists with glaucophane–lawsonite metabasite layers. The blueschist-eclogite unit representing strongly sheared, deeply buried and imbricated tectonic slices of accreted uppermost levels of the oceanic crust with minor metamorphosed serpentinite bodies consists of lawsonite-bearing eclogitic metabasites (approximately 90% of the field), lawsonite eclogites, metagabbros, serpentinites, pelagic marbles, omphacite–glaucophane–lawsonite metapelites and metacherts. The mineral assemblage of the lawsonite eclogite (garnet + omphacite > 70%) is omphacite, garnet, lawsonite, glaucophane, phengite and rutile. Lawsonite eclogite lenses are enclosed by garnet–lawsonite blueschist envelopes.Textural evidence from lawsonite eclogites and country rocks reveals that they did not leave the stability field of lawsonite during subduction and exhumation. The widespread preservation of lawsonite in eclogitic metabasites and eclogites can be attributed to rapid subduction and subsequent exhumation in a low geothermal gradient of the oceanic crust material without experiencing a thermal relaxation. Peak P–T conditions of lawsonite eclogites are estimated at 24 ± 1 kbar and 460 ± 25 °C. These P–T conditions indicate a remarkably low geotherm of 6.2 °C/km corresponding to a burial depth of 74 km.     

On the nature of mantle heterogeneities and discontinuities: evidence from a very dense wide-angle shot record

A seismic survey with a receiver spacing of 50 m provided one of the most densely sampled wide-angle seismic reflection images of the lithosphere. This unique data set, recorded by an 18-km-long spread, reveals that at wide-angles the shallow subcrustal mantle features high amplitude reflectivity which contrasts with a lack of reflectivity at latter travel times. This change in the seismic signature is located at approximately 120–150 km depth, which correlates with the depth estimates of the lithosphere–asthenosphere boundary (LAB) of previous DSS studies. This seismic signature can be simulated by two-layer mantle model. Both layers with similar average velocities differ in their degree of heterogeneity. The shallow heterogeneous layer and the deeper and more homogeneous one correlate with the lithosphere and the asthenosphere, respectively. Studies involving surface outcrops of ultramafic massifs and mantle xenoliths infer that the upper mantle is a heterogeneous mixture of ultramafic rocks (lherzolites, harzburgites, pyroxenites, peridotites, dunites, and small amounts of eclogites). Laboratory measurements of physical properties of these mantle rocks indicate that compositional variations alone can account for the wide-angle reflectivity. A temperature increase would homogenize the mixture, decreasing the seismic reflection properties due to melting processes. It is proposed that this would take place below 120–150 km (1200 °C, the LAB).     

Water content in minerals of mantle xenoliths from the Udachnaya pipe kimberlites (Yakutia)

Distribution of water among the main rock-forming nominally anhydrous minerals of mantle xenoliths of peridotitic and eclogitic parageneses from the Udachnaya kimberlite pipe, Yakutia, has been studied by IR spectroscopy. The spectra of all minerals exhibit vibrations attributed to hydroxyl structural defects. The content of H₂O (ppm) in minerals of peridotites is as follows: 23–75 in olivine, 52–317 in orthopyroxene, 29–126 in clinopyroxene, and 0–95 in garnet. In eclogites, garnet contains up to 833 ppm H₂O, and clinopyroxene, up to 1898 ppm (~ 0.19 wt.%). The obtained data and the results of previous studies of minerals of mantle xenoliths show wide variations in H₂O contents both within different kimberlite provinces and within the Udachnaya kimberlite pipe. Judging from the volume ratios of mineral phases in the studied xenoliths, the water content varies over narrow ranges of values, 38–126 ppm. At the same time, the water content in the studied eclogite xenoliths is much higher and varies widely, 391–1112 ppm.     

Lower crustal intrusions beneath the southern Baikal Rift Zone: Evidence from full-waveform modelling of wide-angle seismic data

The Cenozoic Baikal Rift Zone (BRZ) is situated in south-central Siberia in the suture between the Precambrian Siberian Platform and the Amurian plate. This more than 2000-km long rift zone is composed of several individual basement depressions and half-grabens with the deep Lake Baikal at its centre. The BEST (Baikal Explosion Seismic Transect) project acquired a 360-km long, deep seismic, refraction/wide-angle reflection profile in 2002 across southern Lake Baikal. The data from this project is used for identification of large-scale crustal structures and modelling of the seismic velocities of the crust and uppermost mantle. Previous interpretation and velocity modelling of P-wave arrivals in the BEST data has revealed a multi layered crust with smooth variation in Moho depth between the Siberian Platform (41 km) and the Sayan-Baikal fold belt (46 km). The lower crust exhibits normal seismic velocities around the rift structure, except for beneath the rift axis where a distinct 50–80-km wide high-velocity anomaly (7.4–7.6 ± 0.2 km/s) is observed. Reverberant or “ringing” reflections with strong amplitude and low frequency originate from this zone, whereas the lower crust is non-reflective outside the rift zone. Synthetic full-waveform reflectivity modelling of the high-velocity anomaly suggests the presence of a layered sequence with a typical layer thickness of 300–500 m coinciding with the velocity anomaly. The P-wave velocity of the individual layers is modelled to range between 7.4 km/s and 7.9 km/s. We interpret this feature as resulting from mafic to ultra-mafic intrusions in the form of sills. Petrological interpretation of the velocity values suggests that the intrusions are sorted by fractional crystallization into plagioclase-rich low-velocity layers and pyroxene- and olivine-rich high-velocity layers. The mafic intrusions were probably intruded into the ductile lower crust during the main rift phase in the Late Pliocene. As such, the intrusive material has thickened the lower crust during rifting, which may explain the lack of Moho uplift across southern BRZ.     

Three-dimensional P- and S-wave velocity structures beneath the Japan Islands obtained by high-density seismic stations by seismic tomography

We construct fine-scale 3D P- and S-wave velocity structures of the crust and upper mantle beneath the whole Japan Islands with a unified resolution, where the Pacific (PAC) and Philippine Sea (PHS) plates subduct beneath the Eurasian (EUR) plate. We can detect the low-velocity (low-V) oceanic crust of the PAC and PHS plates at their uppermost part beneath almost all the Japan Islands. The depth limit of the imaged oceanic crust varies with the regions. High-V_P/V_S zones are widely distributed in the lower crust especially beneath the volcanic front, and the high strain rate zones are located at the edge of the extremely high-V_P/V_S zone; however, V_P/V_S at the top of the mantle wedge is not so high. Beneath northern Japan, we can image the high-V subducting PAC plate using the tomographic method without any assumption of velocity discontinuities. We also imaged the heterogeneous structure in the PAC plate, such as the low-V zone considered as the old seamount or the highly seismic zone within the double seismic zone where the seismic fault ruptured by the earthquake connects the upper and lower layer of the double seismic zone. Beneath central Japan, thrust-type small repeating earthquakes occur at the boundary between the EUR and PHS plates and are located at the upper part of the low-V layer that is considered to be the oceanic crust of the PHS plate. In addition to the low-V oceanic crust, the subducting high-V PAC plate is clearly imaged to depths of approximately 250 km and the subducting high-V PHS zone to depths of approximately 180 km is considered to be the PHS plate. Beneath southwestern Japan, the iso-depth lines of the Moho discontinuity in the PHS plate derived by the receiver function method divide the upper low-V layer and lower high-V layer of our model at depths of 30–50 km. Beneath Kyushu, the steeply subducting PHS plate is clearly imaged to depths of approximately 250 km with high velocities. The high-V_P/V_S zone is considered as the lower crust of the EUR plate or the oceanic crust of the PHS plate at depths of 25–35 km and the partially serpentinized mantle wedge of the EUR plate at depths of 30–45 km beneath southwestern Japan. The deep low-frequency nonvolcanic tremors occur at all parts of the high-V_P/V_S zone—within the zone, the seaward side, and the landward side where the PHS plate encounters the mantle wedge of the EUR plate. We prove that we can objectively obtain the fine-scale 3D structure with simple constraints such as only 1D initial velocity model with no velocity discontinuity.     

Imaging global chemical and thermal heterogeneity in the subcontinental lithospheric mantle with garnets and xenoliths: Geophysical implications

Suites of mantle-derived xenoliths in volcanic rocks provide estimates of the geothermal gradient and composition of the subcontinental lithospheric mantle (SCLM) at the time of the volcanic eruption. The development of single-grain thermometry and barometry, applied to xenocryst minerals in volcanic rocks, has greatly expanded the number of localities for which such data can be obtained and made it feasible to map the geology of the SCLM on a broader scale, both vertically and laterally. From garnet xenocrysts, it is possible to derive profiles showing mean values of olivine composition, bulk-rock composition, density and seismic velocities, as well as geotherm parameters and constraints on the thickness of the SCLM. Geochemical profiles, coupled with Re–Os dating of peridotites and their enclosed sulfide minerals, show that Archean or Proterozoic SCLM is preserved at shallow levels beneath many areas of younger tectonothermal age; this implies rapid vertical variations in Vs and Vp with depth, which may affect seismic interpretations. Data from several hundred localities worldwide define a secular evolution in the composition of the SCLM, related to the tectonothermal age of the overlying crust. Archean SCLM is typically strongly depleted in basaltic components, highly magnesian and thick (160–250 km), and has low geotherms; Phanerozoic SCLM is typically fertile (rich in basaltic components), Fe-rich, thin (50–100 km) and has a range of high geotherms; Proterozoic SCLM (much of which may be reworked Archean mantle) tends to be intermediate in all respects. The correlated variations in SCLM fertility, lithospheric thickness and geotherm reinforce the effects of each on seismic velocity, and produce more rapid lateral variations in seismic response than would result from thermal effects alone. These correlations are the key to using seismic tomography images to map the lateral extent of different types of SCLM.     

Anatolian surface wave evaluated at GEOFON Station ISP Isparta, Turkey

We have studied seismic surface waves of 255 shallow regional earthquakes recently recorded at GEOFON station ISP (Isparta, Turkey) and have selected these 52 recordings with high signal-to-noise ratio for further analysis. An attempt was made by the simultaneous use of the Rayleigh and Love surface wave data to interpret the planar crust and uppermost mantle velocity structure beneath the Anatolian plate using a differential least-square inversion technique. The shear-wave velocities near the surface show a gradational change from approximately 2.2 to 3.6 km s^− 1 in the depth range 0–10 km. The mid-crustal depth range indicating a weakly developed low velocity zone has shear-wave velocities around 3.55 km s^− 1. The Moho discontinuity characterizing the crust–mantle velocity transition appears somewhat gradual between the depth range  25–45 km. The surface waves approaching from the northern Anatolia are estimated to travel a crustal thickness of  33 km whilst those from the southwestern Anatolia and part of east Mediterranean Sea indicate a thicker crust at  37 km. The eastern Anatolia events traveled even thicker crust at  41 km. A low sub-Moho velocity is estimated at  4.27 km s^− 1, although consistent with other similar studies in the region. The current velocities are considerably slower than indicated by the Preliminary Reference Earth Model (PREM) in almost all depth ranges.     

Petrofabrics and seismic properties of garnet peridotite from the UHP Sulu terrane (China): Implications for olivine deformation mechanism in a cold and dry subducting continental slab

Lattice-preferred orientations (LPO) of olivine, diopside, enstatite and garnet from the Zhimafang garnet peridotite body in the Sulu ultrahigh-pressure (UHP) metamorphic terrane (China) were measured using the electron backscatter diffraction (EBSD) technique. The peridotite was captured from a mantle wedge immediately adjacent the subducted Yangtze slab and then experienced the UHP metamorphism at 750–950 °C and 4–7 GPa. The olivine LPO is characterized by the [001] axis close to the stretching lineation and the (100) plane subparallel to the foliation, indicating the prevailing of (100) [001] slip. Enstatite LPO displays the dominance of (100) [001] slip. Diopside developed complex LPO patterns that are difficult to explain using a single slip system of (100) [001]. Garnet is almost randomly oriented due to its low volume fractions, cubic symmetry and the presence of numerous slip systems. Calculated seismic properties of the peridotite yield a maximum P-wave velocity normal to the foliation and a minimum along the foliation, with anisotropy up to 8% in strongly sheared samples. The S-wave velocity pattern is complex but the fast polarization plane generally normal to the foliation. The inferred shear sense from the olivine LPO is top-to-SE, in contrary to exhumation-induced top-to-NW thrusting recorded in the quartz LPO, implying that the olivine LPO formed at early UHP metamorphic conditions. The olivine crystals have relatively low water contents (141–475 H/10⁶ Si), indicating a fluid-deficient environment for the LPO formation. The present study suggests that a combination of low temperature and UHP plays a much more important role than the water content to promote the activation of (100) [001] slip in olivine.     

Pressure dependence and anisotropy of P-wave velocities in ultrahigh-pressure metamorphic rocks from the Dabie–Sulu orogenic belt (China): Implications for seismic properties of subducted slabs and origin of mantle reflections

The compressional wave velocities (Vp), pressure derivatives (Vp′) and anisotropy (A) of three types of eclogites and country rocks from the Dabie–Sulu ultrahigh-pressure (UHP) metamorphic belt, China, have been measured under confining pressures up to 800 MPa. Type-1 eclogites, which are coarse-grained and subjected to almost no retrograde metamorphism, experienced recovery-accommodated dislocation creep at peak metamorphic conditions (in the diamond stability field). Type-2 eclogites are fine-grained reworked Type-1 materials that experienced recrystallization-accommodated dislocation creep under quartz/coesite boundary conditions during the early stage of exhumation. Type-3 eclogites are retrogressed samples that were overprinted by significant amphibolite facies metamorphism during a late stage of exhumation within the crust. Type-1 eclogites are richer in Al₂O₃ and MgO but poorer in SiO₂ and Na₂O than Type-2 and Type-3 eclogites. Anisotropy of Type-1 and Type-2 eclogites is generally low (<4%) because volumetrically important garnet is elastically quasi-isotropic, while Type-3 eclogites can exhibit high anisotropy (>10%) due to the presence of strongly anisotropic retrograde minerals such as amphibole, plagioclase and mica. The transition of the pressure dependence of velocity from the poroelastic to elastic regimes occurs at a critical pressure (P_c), which depends mainly on the density and distribution of microcracks and in turn on the exhumation history of rocks. The Vp–pressure relationship can be expressed by Vp=a(lnP)²+blnP+c (P≤P_c) and Vp=V₀+DP (P≥P_c), where P is the confining pressure, a and b are constants describing the closure of microcracks below P_c, c is the velocity when P is equal to one (MPa), V₀ is the projected velocity of a crack-free sample at room pressure, and D is the intrinsic pressure derivative above P_c. When data are curve-fit, pressure derivatives and anisotropy as functions of pressure are determined. The average Vp of the eclogites in the linear regime is 8.42+1.41×10⁻⁴P for Type-1, 7.80+1.58×10⁻⁴P for Type-2, and 7.33+2.04×10⁻⁴P for Type-3, where Vp is in km/s and P in MPa. The decrease in V₀ and increase in D from Type-1 to Type-3 eclogites are attributed to a decrease in garnet content and an increase in retrograde minerals. The NE–SW trending, NW-dipping, slab-like high Vp anomaly (8.72 km/s at a depth of 71 km) which extends from the Moho to at least 110 km beneath the Dabie–Sulu region, can be interpreted as the remnant of a subducted slab which is dominated by Type-1 eclogites and has frozen in the upper mantle since about 200–220 Ma. Such relic crustal materials, subducted and preserved as eclogite layers intercalated with felsic gneiss, garnet–jadeite quartzite, marble and serpentinized peridotite, could be responsible for regionally observed seismic reflectors in the upper mantle.     

Upper mantle structure of the South American continent and neighboring oceans from surface wave tomography

We present a new three-dimensional SV-wave velocity model for the upper mantle beneath South America and the surrounding oceans, built from the waveform inversion of 5850 Rayleigh wave seismograms. The dense path coverage and the use of higher modes to supplement the fundamental mode of surface waves allow us to constrain seismic heterogeneities with horizontal wavelengths of a few hundred kilometres in the uppermost 400 km of the mantle.The large scale features of our tomographic model confirm previous results from global and regional tomographic studies (e.g. the depth extent of the high velocity cratonic roots down to about 200–250 km).Several new features are highlighted in our model. Down to 100 km depth, the high velocity lid beneath the Amazonian craton is separated in two parts associated with the Guyana and Guapore shields, suggesting that the rifting episode responsible for the formation of the Amazon basin has involved a significant part of the lithosphere. Along the Andean subduction belt, the structure of the high velocity anomaly associated with the sudbduction of the Nazca plate beneath the South American plate reflects the along-strike variation in dip of the subducting plate. Slow velocities are observed down to about 100 km and 150 km at the intersection of the Carnegie and Chile ridges with the continent and are likely to represent the thermal anomalies associated with the subducted ridges. These lowered velocities might correspond to zones of weakness in the subducted plate and may have led to the formation of “slab windows” developed through unzipping of the subducted ridges; these windows might accommodate a transfer of asthenospheric mantle from the Pacific to the Atlantic ocean. From 150 to 250 km depth, the subducting Nazca plate is associated with high seismic velocities between 5°S and 37°S. We find high seismic velocities beneath the Paraná basin down to about 200 km depth, underlain by a low velocity anomaly in the depth range 200–400 km located beneath the Ponta Grossa arc at the southern tip of the basin. This high velocity anomaly is located southward of a narrow S-wave low velocity structure observed between 200 and 500–600 km depth in body wave studies, but irresolvable with our long period datasets. Both anomalies point to a model in which several, possibly diachronous, plumes have risen to the surface to generate the Paraná large igneous province (LIP).     

Seismic velocity discontinuities in the crust and uppermost mantle beneath the Kanto district, central Japan, identified from receiver function imaging and repeating earthquake activity

We constructed vertical cross-sections of depth-converted receiver function images to estimate the seismic velocity structure of the crust and uppermost mantle beneath the Kanto district, central Japan. Repeating earthquake data for the plate boundary were also used to estimate geometries of the subducting Philippine Sea plate and the subducting Pacific plate. As a result, we present images of some major seismic discontinuities. The upper boundary of the Pacific plate dips to the northwest in northern Kanto and to the west–southwest in southern Kanto with some undulations. On the other hand, the upper boundary of the Philippine Sea plate as a whole dips to the northwest. However, it is concave to the northeast in the southern Boso peninsula. We suggest that the low-velocity mantle wedge may be indicated on the top of both subducting plates. Plate thickness gradually decreases to the northeast. The northeastern end of the Philippine Sea plate is interpreted to be at depths of 45–90 km. The Moho discontinuity in the overriding plate is deeper than 25 km in the northern Kanto. It contacts the subducting Philippine Sea plate in the southwestern part near 35.8°N.     

Strain localisation in the subcontinental mantle — a ductile alternative to the brittle mantle

It is now admitted that the high strength of the subcontinental uppermost mantle controls the first order strength of the lithosphere. An incipient narrow continental rift therefore requires an important weakening in the subcontinental mantle to promote lithosphere-scale strain localisation and subsequent continental break-up. Based on the classical rheological layering of the continental lithosphere, the origin of a lithospheric mantle shear/fault zone has been attributed to the existence of a brittle uppermost mantle. However, the lack of mantle earthquakes and the absence of field occurrences in the mantle fault zone led to the idea of a ductile-related weakening mechanism, instead of brittle-related, for the incipient mantle strain localisation. In order to provide evidence for this mechanism, we investigated the microstructures and lattice preferred orientations of mantle rocks in a kilometre-scale ductile strain gradient in the Ronda Peridotites (Betics cordillera, Spain). Two main features were shown: 1) grain size reduction by dynamic recrystallisation is found to be the only relevant weakening mechanism responsible for strain localisation and 2), with increasing strain, grain size reduction is coeval with both the scattering of orthopyroxene neoblasts and the decrease of the olivine fabric strength (LPO). These features allow us to propose that grain boundary sliding (GBS) partly accommodates dynamic recrystallisation and subsequent grain size reduction.A new GBS-related experimental deformation mechanism, called dry-GBS creep, has been shown to accommodate grain size reduction during dynamic recrystallisation and to induce significant weakening at low temperatures (T < 800 °C). The present microstructural study demonstrates the occurrence of the grain size sensitive dry-GBS creep in natural continental peridotites and allows us to propose a new rheological model for the subcontinental mantle. During dynamic recrystallisation, the accommodation of grain size reduction by three competing deformation mechanisms, i.e., dislocation, diffusion and dry-GBS creeps, involves a grain size reduction controlled by the sole dislocation creep at high temperatures (> 800 °C), whereas dislocation creep and dry-GBS creep, are the accommodating mechanisms at low temperatures (< 800 °C). Consequently, weakening is very limited if the grain size reduction occurs at temperatures higher than 800 °C, whereas a large weakening is expected in lower temperatures. This large weakening related to GBS creep would occur at depths lower than 60 km and therefore provides an explanation for ductile strain localisation in the uppermost continental mantle, thus providing an alternative to the brittle mantle.     

Tomographic imaging of hydrated crust and mantle in the subducting Pacific slab beneath Hokkaido, Japan: Evidence for dehydration embrittlement as a cause of intraslab earthquakes

We estimate detailed three-dimensional seismic velocity structures in the subducting Pacific slab beneath Hokkaido, Japan, using a large number of arrival-time data from 6902 local earthquakes. A remarkable low-velocity layer with a thickness of ~ 10 km is imaged at the uppermost part of the slab and is interpreted as hydrated oceanic crust. The layer gradually disappears at depths of 70–80 km, suggesting the breakdown of hydrous minerals there. We find prominent low-velocity anomalies along the lower plane of the double seismic zone and above the aftershock area of the 1993 Kushiro-oki earthquake (M7.8). Since seismic velocities of unmetamorphosed peridotite are much higher than the observations, hydrous minerals are expected to exist in the lower plane as well as the hypocentral area of the 1993 earthquake. On the other hand, regions between the upper and lower planes, where seismic activity is not so high compared to the both planes, show relatively high velocities comparable to those of unmetamorphosed peridotite. Our observations suggest that intermediate-depth earthquakes occur mainly in regions with hydrous minerals, which support dehydration embrittlement hypothesis as a cause of earthquake in the subducting slab.