2D model of the crust and uppermost mantle along rift profile, Siberian craton

A two-dimensional model of the crust and uppermost mantle for the western Siberian craton and the adjoining areas of the Pur-Gedan basin to the north and Baikal Rift zone to the south is determined from travel time data from recordings of 30 chemical explosions and three nuclear explosions along the RIFT deep seismic sounding profile. This velocity model shows strong lateral variations in the crust and sub-Moho structure both within the craton and between the craton and the surrounding region. The Pur-Gedan basin has a 15-km thick, low-velocity sediment layer overlying a 25-km thick, high-velocity crystalline crustal layer. A paleo-rift zone with a graben-like structure in the basement and a high-velocity crustal intrusion or mantle upward exists beneath the southern part of the Pur-Gedan basin. The sedimentary layer is thin or non-existent and there is a velocity reversal in the upper crust beneath the Yenisey Zone. The Siberian craton has nearly uniform crustal thickness of 40–43 km but the average velocity in the lower crust in the north is higher (6.8–6.9 km/s) than in the south (6.6 km/s). The crust beneath the Baikal Rift zone is 35 km thick and has an average crustal velocity similar to that observed beneath the southern part of craton. The uppermost mantle velocity varies from 8.0 to 8.1 km/s beneath the young West Siberian platform and Baikal Rift zone to 8.1–8.5 km/s beneath the Siberian craton. Anomalous high Pn velocities (8.4–8.5 km/s) are observed beneath the western Tunguss basin in the northern part of the craton and beneath the southern part of the Siberian craton, but lower Pn velocities (8.1 km/s) are observed beneath the Low Angara basin in the central part of the craton. At about 100 km depth beneath the craton, there is a velocity inversion with a strong reflecting interface at its base. Some reflectors are also distinguished within the upper mantle at depth between 230 and 350 km.     

Phlogopite in mantle xenoliths and kimberlite from the Grib pipe,Arkhangelsk province,Russia:Evidence for multi-stage mantle metasomatism and origin of phlogopite in kimberlite

We present petrography and mineral chemistry for both phlogopite,from mantle-derived xenoliths(garnet peridotite,eclogite and clinopyroxene-phlogopite rocks)and for megacryst,macrocryst and groundmass flakes from the Grib kimberlite in the Arkhangelsk diamond province of Russia to provide new insights into multi-stage metasomatism in the subcratonic lithospheric mantle(SCLM)and the origin of phlogopite in kimberlite.Based on the analysed xenoliths,phlogopite is characterized by several generations.The first generation(Phil)occurs as coarse,discrete grains within garnet peridotite and eclogite xenoliths and as a rock-forming mineral within clinopyroxene-phlogopite xenoliths.The second phlogopite generation(Phl2)occurs as rims and outer zones that surround the Phil grains and as fine flakes within kimberlite-related veinlets filled with carbonate,serpentine,chlorite and spinel.In garnet peridotite xenoliths,phlogopite occurs as overgrowths surrounding garnet porphyroblasts,within which phlogopite is associated with Cr-spinel and minor carbonate.In eclogite xenoliths,phlogopite occasionally associates with carbonate bearing veinlet networks.Phlogopite,from the kimberlite,occurs as megacrysts,macrocrysts,microcrysts and fine flakes in the groundmass and matrix of kimberlitic pyroclasts.Most phlogopite grains within the kimberlite are characterised by signs of deformation and form partly fragmented grains,which indicates that they are the disintegrated fragments of previously larger grains.Phil,within the garnet peridotite and clinopyroxene-phlogopite xenoliths,is characterised by low Ti and Cr contents(TiO_21 wt.%,Cr_2 O_31 wt.% and Mg# = 100 × Mg/(Mg+ Fe)92)typical of primary peridotite phlogopite in mantle peridotite xenoliths from global kimberlite occurrences.They formed during SCLM metasomatism that led to a transformation from garnet peridotite to clinopyroxene-phlogopite rocks and the crystallisation of phlogopite and high-Cr clinopyroxene megacrysts before the generation of host-kimberlite magmas.One of the possible processes to generate low-Ti-Cr phlogopite is via the replacement of garnet during its interaction with a metasomatic agent enriched in K and H_2O.Rb-Sr isotopic data indicates that the metasomatic agent had a contribution of more radiogenic source than the host-kimberlite magma.Compared with peridotite xenoliths,eclogite xenoliths feature low-Ti phlogopites that are depleted in Cr_2O_3 despite a wider range of TiO_2 concentrations.The presence of phlogopite in eclogite xenoliths indicates that metasomatic processes affected peridotite as well as eclogite within the SCLM beneath the Grib kimberlite.Phl2 has high Ti and Cr concentrations(TiO_22 wt.%,Cr_2O_31 wt.% and Mg# = 100× Mg/(Mg + Fe)92)and compositionally overlaps with phlogopite from polymict brecc:ia xenoliths that occur in global kimberlite formations.These phlogopites are the product of kimberlitic magma and mantle rock interaction at mantle depths where Phl2 overgrew Phil grains or crystallized directly from stalled batches of kimberlitic magmas.Megacrysts,most macrocrysts and microcrysts are disintegrated phlogopite fragments from metasomatised peridotite and eclogite xenoliths.Fine phlogopite flakes within kimberlite groundmass represent mixing of high-Ti-Cr phlogopite antecrysts and high-Ti and low-Cr kimberlitic phlogopite with high Al and Ba contents that may have formed individual grains or overgrown antecrysts.Based on the results of this study,we propose a schematic model of SCLM metasomatism involving phlogopite crystallization,megacryst formation,and genesis of kimberlite magmas as recorded by the Grib pipe.     

Test of the upper mantle low velocity layer in Siberia with surface waves

The existence of the upper mantle low velocity layer (LVL) below 100 km depth in cratonic areas is tested with surface waves dispersion curves. Given the ambient noise we find that a pronounced LVL (80 km thick and 2% velocity reduction or 40 km thick and 5% velocity reduction) can be distinguished from a constant velocity model by comparison of the fundamental mode group velocities, whereas a thin LVL (less than 40 km thick) with small velocity contrast (less than 2%) cannot be resolved. The fundamental modes of Love and Rayleigh waves have similar properties and, in general, the phase velocity differences are smaller than the standard error. Phase velocity alone cannot discriminate between the models, and the group velocity is in general more sensitive to the velocity structure than the phase velocity. The higher modes at short periods could potentially determine a LVL but in reality it is difficult to obtain sufficiently accurate measurements. We invert the synthetic dispersion curves by the non-linear Hedgehog inversion method. A pronounced LVL (more than 40 km thick and with a strong velocity contrast of about 5%) is detectable by the non-linear inversion but for a thin LVL with a strong velocity contrast it is not possible to resolve both velocity and thickness. In the inversions all solutions include a LVL for models with a pronounced LVL, whereas the solution space includes models with and without a LVL for models with a zero or positive gradient velocity–depth structure.We invert also real data with travel path across the Siberian craton with the Hedgehog method. Almost all solutions include a LVL in the depth range of 80–150 km with a velocity contrast up to 2% to the surrounding intervals. Hence, the LVL appears to be a common feature of the Siberian upper mantle, although a constant velocity at the same depth range cannot be totally excluded. Despite low resolution at large depth, a pronounced asthenospheric LVL below a depth of about 225 km is a constant characteristic of the set of solutions.     

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.     

西伯利亚雅库特金伯利岩中四颗古老锆石巨晶的微量元素和Hf同位素特征——对区分金伯利岩中锆石巨晶的启示

金伯利岩中的锆石按照颗粒大小可以分为细粒锆石（一般小于200μm）和巨晶锆石（一般大于500μm）。前人的研究结果显示在金伯利岩中粒径较大的巨晶锆石的U-Pb体系在高温的地幔中一直保持着开放状态,直到寄主金伯利岩浆的喷发才使地幔锆石的U-Pb体系封闭,因此这些巨晶锆石是确定金伯利岩年龄的重要矿物之一。然而,近年来的研究表明,金伯利岩中还存在一些时代远老于金伯利岩年龄的锆石,也具有较大的粒径（以下称古老锆石巨晶）。它们的存在无疑影响了利用锆石U-Pb方法确定金伯利岩年龄的准确性。本文以西伯利亚雅库特（Yakutia）金伯利岩省中的四颗古老锆石巨晶为研究对象,通过形态学、年代学、微量元素和Hf同位素组成,讨论古老锆石巨晶的来源。同时,我们统计和对比了全球多个金伯利岩中能够确定金伯利岩年龄的锆石和古老锆石巨晶的形态学、U、Th含量和微量元素组成、Hf-O同位素等特征。研究结果显示,金伯利岩中的古老锆石巨晶的t_DM年龄和O同位素组成与可以用来确定金伯利岩年龄的锆石巨晶具有明显的差别。这些手段在未来的研究中可以用来区分可确定金伯利岩年龄的锆石巨晶和古老的锆石巨晶。

     

Non-transparent uppermost mantle in the island-arc region of Japan

In Japan, the crust and uppermost mantle seismic character is yet unimaged although many refraction surveys have been recorded. The longest seismic profiles are analyzed. A remarkable feature, a long-duration coda wave after the PmP wave (reflected wave at the Moho boundary), is observed on the record sections. Several possible models are considered to explain the long-duration coda wave. The model with many scatterers located in the uppermost mantle explains the observed data well while the undulating Moho and continuous layering models do not account for some aspects of the observed data. The scatterer distributed uppermost mantle is not consistent with that of continental region which is often characterized as transparent. We estimate the scattering coefficient of the uppermost mantle and crust using simulations. The scattering coefficients obtained for upper crust, lower crust, and uppermost mantle are 0.01, 0.02, and 0.025, respectively. The scattering coefficient of the uppermost mantle is slightly larger than that of lower crust, which is characterized as being reflective. The many scatterers in the uppermost mantle might be related to magmatism in Japan. This will be one of the important observations for understanding formation processes of the Moho boundary and uppermost mantle in the island-arc environment.     

Interaction between protokimberlite melts and mantle lithosphere: Evidence from mantle xenoliths from the Dalnyaya kimberlite pipe,Yakutia (Russia)

The Dalnyaya kimberlite pipe(Yakutia,Russia) contains mantle peridotite xenoliths(mostly Iherzolites and harzburgites) that show both sheared porphyroclastic(deformed) and coarse granular textures,together with ilmenite and clinopyroxene megacrysts.Deformed peridotites contain high-temperature Fe-rich clinopyroxenes,sometimes associated with picroilmenites,which are products of interaction of the lithospheric mantle with protokimberlite related melts.The orthopyroxene-derived geotherm for the lithospheric mantle beneath Dalnyaya is stepped similar to that beneath the Udachnaya pipe.Coarse granular xenoliths fall on a geotherm of 35 mWm-2 whereas deformed varieties yield a 45 mWm-2)geotherm in the 2-7.5 GPa pressure interval.The chemistry of the constituent minerals including garnet,olivine and clinopyroxene shows trends of increasing Fe~#(=Fe/(Fe+Mg))with decreasing pressure.This may suggest that the interaction with fractionating protokimberlite melts occurred at different levels.Two major mantle lithologies are distinguished by the trace element patterns of their constituent minerals,determined by LA-ICP-MS.Orthopyroxenes,some clinopyroxenes and rare garnets are depleted in Ba,Sr,HFSE and MREE and represent relic lithospheric mantle.Re-fertilized garnet and clinopyroxene are more enriched.The distribution of trace elements between garnet and clinopyroxene shows that the garnets dissolved primary orthopyroxene and clinopyroxene.Later high temperature clinopyroxenes related to the protokimberlite melts partially dissolved these garnets.Olivines show decreases in Ni and increases in Al,Ca and Ti from Mg-rich varieties to the more Fe-rich,deformed and refertilized ones.Minerals showing higher Fe~#(0.11-0.15) are found within intergrowths of low-Cr ilmenite-clinopyroxene-garnet related to the crystallization of protokimberlite melts in feeder channels.In P-f(O_2) diagrams,garnets and Cr-rich clinopyroxenes indicate reduced conditions at the base of the lithosphere at-5 log units below a FMQ buffer.However,Cr-poor clinopyroxenes,together with ilmenite and some Fe-Ca-rich garnets,demonstrate a more oxidized trend in the lower part of lithosphere at-2 to 0 log units relative to FMQ.Clinopyroxenes from xenoliths in most cases show conditions transitional between those determined for garnets and megacrystalline Cr-poor suite.The relatively low diamond grade of Dalnyaya kimberlites is explained by a high degree of interaction with the oxidized protokimberlite melts,which is greater at the base of the lithosphere.     

Seismic reflection imaging of thin, kimberlite dykes and sills: exploration and deposit characterization of the Snap Lake dyke, Canada

Seismic reflection techniques are, for the first time, used to image a thin, diamondiferous, kimberlite dyke from subcrop to depths greater than 1300 m. Exploration for vertical kimberlite pipes generally utilizes potential field techniques that often fail to reveal subhorizontal or shallow-dipping intrusions. In contrast, seismic reflection methods are well suited for imaging targets with this geometry. Therefore, in order to evaluate seismic reflection as a tool for subhorizontal kimberlite dyke/sill exploration and mine planning, a feasibility study and subsequent seismic survey were undertaken on the diamondiferous Snap Lake dyke (Northwest Territories, Canada). A substantial drilling program has mapped the dyke as a gently dipping sheet that averages 2–3 m in thickness. The detailed structural and composition data available at Snap Lake provides a unique opportunity to test reflection techniques on a well-sampled deposit. The feasibility study involved measuring P-velocities and densities of cores drilled from the kimberlite and host rocks. These data were used to model reflection amplitudes, evaluate resolution limitations, and determine the acquisition parameters for the reflection survey. Two 2-D lines were acquired that provide comparative datasets for different sources (explosive and vibroseis) and ground types (land and lake ice). In addition, the exploration-scale survey incorporated high fold (40–260 nominal) and long offsets (3260 m). The explosive-source profile recorded on land yielded a superb image of the dyke from depths of 60 m to more than 1300 m over a lateral distance of 5700 m. The seismic image correlates well with adjacent drill hole data and adds considerable detail to the topography of the kimberlite sheet determined by drilling. The vibroseis source also imaged the dyke, but only when sources and geophones were on land; the dyke was not imaged beneath the ice due to reverberation and attenuation effects. The frequency response and unusually strong reflection amplitudes from the dyke indicate the importance of tuning effects and multiples for this type of target and acquisition environment. Apparent correlations between reflection amplitudes and dyke structure (e.g., thickness, feathering, 3-D geometry) suggest that seismic reflection data may be valuable for guiding drilling programs. The results demonstrate that, in the appropriate situation, seismic methods have great potential for use in kimberlite exploration, subsurface mapping, and detailed imaging for mine development purposes.     

The heterogeneous upper mantle low velocity zone

The upper mantle low velocity zone (LVZ) is a depth interval with slightly reduced seismic velocity compared to the surrounding depth intervals. The zone is present below a relatively constant depth of 100 km in most continental parts of the world, both in cratonic areas with high average velocity and tectonically active areas with low average velocity. Evidence for the low velocity zone arises from controlled and natural source seismology, including studies of surface waves and of primary and multiple reflections of body waves from the bounding interfaces, calculations of receiver functions, and absolute velocity tomography. The available data indicates a more pronounced reduction in seismic velocity and Q-value for S-waves than P-waves as well as high electrical conductivity in the LVZ. Seismic waves are strongly scattered by the zone, which demonstrates the existence of small-scale heterogeneity. The depth to the base of the LVZ is systematically shallower in cold, stable cratonic areas than in hot, active regions of the world. Because of its global occurrence below a relative constant depth of 100 km, the LVZ cannot be explained by metamorphic or compositional variation and rheological changes. Calculated upper mantle temperatures indicate that the rocks are close to the solidus in an interval with variable thickness below 100 km depth, provided that the rocks contain water and carbon dioxide. The presence of, even small amounts of such fluids in the mantle rocks will lower the solidus by several hundred degrees and introduce a characteristic kink on the solidus curve around 80–100 km depth. The seismic velocities and Q-values are significantly reduced of rocks, which are close to the solidus or contain small amounts of partial melt. Hence, the LVZ may be explained by upper mantle temperatures being close to the solidus in a depth interval below 100 km. Assuming that the rocks contain only limited amounts of fluids, this mechanism may explain the low velocities, Q-values, and resistivity, as well as the intrinsic scattering, and the characteristic variation in thickness of the low velocity zone.     

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

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 < 2-km-thick interface at the top of the subducting plate. We show that another zone of reflectors exists east of Vancouver Island that is approximately 8 km thick, and identified as the D reflectors. These overlie the E reflectors; together the two zones define the landward part of the duplex. The combined zones reach depths as great as 50 km. The duplex structure extends for more than 120 km perpendicular to the margin, has an along-strike extent of 80 km, and at depths between 30 km and 50 km the duplex structure correlates with a region of anomalously deep seismicity, where velocities are less than 7000 m s^− 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.     

Geochemical and oxygen isotope signatures of mantle corundum megacrysts from the Mbuji-Mayi kimberlite,Democratic Republic of Congo,and the Changle alkali basalt,China

Oxygen isotope signatures of ruby and sapphire megacrysts, combined with trace-element analysis, from the Mbuji-Mayi kimberlite, Democratic Republic of Congo, and the Changle alkali basalt, China, provide clues to specify their origin in the deep Earth. At Mbuji-Mayi, pink sapphires have δ¹⁸O values in the range 4.3 to 5.4‰ (N = 10) with a mean of 4.9 ± 0.4‰, and rubies from 5.5 to 5.6‰ (N = 3). The Ga/Mg ratio of pink sapphires is between 1.9 and 3.9, and in rubies, between 0.6 and 2.6. The blue or yellow sapphires from Changle have δ¹⁸O values from 4.6 to 5.2 ‰, with a mean of 4.9 ± 0.2‰ (N = 9). The Ga/Mg ratio is between 5.7 and 11.3. The homogenous isotopic composition of ruby suggests a derivation from upper mantle xenoliths (garnet lherzolite, pyroxenite) or metagabbros and/or lower crustal garnet clinopyroxenite eclogite-type xenoliths included in kimberlites. Data from the pink sapphires from Mbuji-Mayi suggest a mantle origin, but different probable protoliths: either subducted oceanic protolith transformed into eclogite with δ¹⁸O values buffered to the mantle value, or clinopyroxenite protoliths in peridotite. The Changle sapphires have a mantle O-isotope signature. They probably formed in syenitic magmas produced by low degree partial melting of a spinel lherzolite source. The kimberlite and the alkali basalt acted as gem conveyors from the upper mantle up to the surface.     

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

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 (P_n and S_n). 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 (P_cc) 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 P_n 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.     

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.     

Seismic wave velocity and anisotropy of serpentinized peridotite in the Oman ophiolite

Shallow seismic measurements in harzburgite from the Oman ophiolite performed in a zone where the maximum horizontal anisotropy is expected (vertical foliation and horizontal lineation) point to a dominant dependence of seismic properties on fracturing.

Optical microscopy studies show that microcracks are guided by the serpentine (lizardite) penetrative network oriented subparallel to the harzburgite foliation and subperpendicular to the mineral lineation, and that serpentine (lizardite) vein filling has a maximum concentration of (001) planes parallel to the veins walls. The calculated elastic properties of the oriented alteration veins filled with serpentine in an anisotropic matrix formed by oriented crystals of olivine and orthopyroxene are compared with seismic velocities measured on hand specimens.

Laboratory ultrasonic data indicate that open microcracks are closed at 100 MPa pressure, e.g. (J. Geophys. Res. 65, (1960) 1083) and (Proc. ODP Sci. Results Leg 118, (1990) 227). Above this pressure, laboratory measurements and modeling show that P-compressional and S-shear wave velocities are mainly controlled by the mineral preferred orientation. Veins sealed with serpentine are effective in slightly lowering P and S velocities and increasing anisotropy. The penetrative lizardite network does not affect directly the geometry of seismic anisotropy, but contributes indirectly in the fact that this network controls the microcrack orientations.

Comparison between seismic measurements of peridotite and gabbro in the same conditions suggest that P- and S-waves anisotropies are a possible discriminating factor between the two lithologies in the suboceanic lithosphere.     

Lithological interpretation of crustal composition in the Fennoscandian Shield with seismic velocity data

In this study, we report the results of an investigation of lithological interpretation of the crust in the central Fennoscandian Shield (in Finland) using seismic wide-angle velocity models and laboratory measurements on P- and S-wave velocities of different rock types. The velocities adopted from wide-angle velocity models were compared with laboratory velocities of different rock types corrected for the crustal PT conditions in the study area. The wide-angle velocity models indicate that the P-wave velocity does not only increase step-wise at boundaries of major crustal layers, but there is also gradual increase of velocity within the layers. On the other hand, the laboratory measurements of velocities indicate that no single rock type is able to provide the gradual downward increasing trends. Thus, there must be gradual vertical changes in rock composition. The downward increase of velocities indicates that the composition of the crust becomes gradually more mafic with increasing depth. We have calculated vertical velocity profiles for a range of possible crustal lithological compositions. The Finnish crustal velocity profiles require a more mafic composition than an average global continental model would suggest. For instance, on the SVEKA'81 transect, the calculated models suggest that the crustal velocity profiles can be simulated with rock type mixtures where the upper crust consists of felsic gneisses and granitic–granodioritic rocks with a minor contribution of amphibolite and diabase. In the middle crust, the amphibolite proportion increases. The lower crust consists of tonalitic gneiss, mafic garnet granulite, hornblendite, pyroxenite and minor mafic eclogite. Assuming that these rock types are present in sufficiently extensive and thick layers, they would also have sufficiently high acoustic reflection coefficients for generating the generally well-developed reflectivity in the crust in the central part of the shield. Density profiles calculated from the lithological models suggest that there is practically no density contrast at Moho in areas of the high-velocity lower crust. Comparison of reflectors from FIRE-1 and FIRE-3 transects and the velocity model from SVEKA'81 wide-angle transect indicated that the reflectors correlate with velocity layering, but the three-dimensional structures of the crust complicate such comparisons.     

Significance of inhomogeneous compression of mineral substance for explaining crystalline crust and uppermost mantle velocity profiles

Consideration of the Voronezh Crystalline Massif shows that three-layered seismic models and velocity profiles with a characteristic distribution of layer velocities and their gradients in each layer, reflect a crystalline crust and uppermost mantle structure formed through regional metamorphism and magmatism processes.Velocity and seismic parameter vertical gradients from the lithosphere profile are several times higher than those for homogeneous compression of minerals and rocks.For the identified layers the contribution of different factors for each seismic parameter gradient value is estimated. Anomalously high rates of increase with depth for V_p, V_s and ø are shown to be caused by the combined effects of structure-phase transformations and substance composition variations in lithospheric complex mineral parageneses.     

The evolution of refertilized lithospheric mantle beneath the northeastern Siberian craton: Links between mantle metasomatism,thermal state and diamond potential

Although the diamond potential of cratons is linked mainly to thick and depleted Archean lithospheric keels, there are examples of craton-edge locations and circum-cratonic Proterozoic terranes underlain by diamondiferous mantle. Here, we use the results of comprehensive major and trace-element studies of detrital garnets from diamond-rich Late Triassic (Carnian) sedimentary rocks in the northeastern Siberia to constrain the thermal and chemical state of the pre-Triassic mantle and its ability to sustain the diamond storage. The studied detrital mantle-derived garnets are dominated by low- to medium-Cr lherzolitic (～45%) and low-Cr megacrystic (～39%) chemistries, with a significant proportion of eclogitic garnets (～11%), and only subordinate contribution from harzburgitic garnets (～5%) with variable Cr₂O₃ contents (1.2–8.4 wt.%). Low-Cr megacrysts display uniform, “normal” rare-earth element (REE) patterns with no Eu/Eu* anomalies, systematic Zr and Ti enrichment (mainly within 2.5–5), which are evidence of their crystallization from deep metasomatic melts. Lherzolitic (G9) garnets exhibit normal or humped to MREE-depleted sinusoidal REE patterns and elevated Nd/Y (up to 0.33–0.41) and Zr/Y ratios (up to 7.62). Rare low- to high-Cr harzburgitic (G10) garnets have primarily “depleted”, sinusoidal REE-patterns, low Ti, Y and HREE, but vary significantly in Zr-Hf, Ti and MREE-HREE contents, Nd/Y (within 0.1–2.4) and Zr/Y (1.53–19.9) ratios. The observed trends of chemical enrichment from the most depleted, harzburgitic garnets towards lherzolitic (including high-Ti high-Cr G11-type) garnets and megacrysts result from either voluminous high-temperature metasomatism by plume-derived silicate melts or recurrent mobilization of less voluminous kimberlitic or related carbonated mantle melts, rather than the initially primitive, fertile nature of the Proterozoic SCLM. Calculated Ni-in-garnet temperatures (primarily within ～1150–1250 °C) indicate their derivation from at least ～220 km thick Cr-undersaturated lithosphere at the relevant Devonian to Triassic thermal flow of ～45 mW/m² or cooler. We suggest the existence of rare harzburgitic domains in the primarily lherzolitic diamond-facies SCLM beneath the northeastern Siberian craton at least by Triassic, whereas the abundance of eclogitic garnets, predominance of E-type inclusions in placer diamonds and specific morphologies argue for diamondiferous eclogites occurring within a ～50–65 kbar diamond window of the Olenek province by the same time.     

Must magmatic intrusion in the lower crust produce reflectivity?

The Færoe–Iceland Ridge (FIR) provides a laboratory in which to investigate the reflectivity and velocity structure of thick crust generated above a mantle plume in order to constrain models of underplating and the origins of lower-crustal layering in an environment dominated by young igneous processes. Over 600 km of common midpoint (cmp) data were collected along and across the FIR using a large airgun array with a 240-channel streamer. The interpretation of these data has been integrated with a velocity model of the crust and upper mantle along the FIR obtained from wide-angle seismic arrivals into ocean bottom and land seismometers. Due to the intermediate water depths and the presence of basalt near the water bottom, specialized processing steps were required for the cmp data. A wave equation-based multiple attenuation scheme was applied to the prestack data, which used a forward model of the multiple series to predict and attenuate multiple energy. Array simulations were applied in the shot and receiver domains in order to minimize spatial aliasing and reduce low apparent-velocity noise. Most of the sections over the central (oceanic) portion of the FIR show no pronounced reflectivity, although occasional Moho and/or lower-crustal reflections are observed. We believe that the poor reflectivity results largely from a lack of physical property contrasts rather than being an effect of acquisition or processing, although we also conclude that residual energy from strong multiple reflection remains in the final sections. Amplitude decay and reflection strength vary along the FIR, but there is good signal-to-noise ratio to travel times of at least 9 s (i.e., into the lower crust), implying that the reduced reflectivity beneath the main part of the FIR is not an artifact of signal penetration loss. We conclude that the addition of melt to the lower crust along the trace of the plume apparently did not produce strong physical property contrasts in the lower crust, where little reflectivity is apparent. Perhaps this was because the entire crust was hot at the time of formation. In contrast, igneous intrusion into preexisting continental crust (at the Færoe Islands end of the FIR) and into older igneous crust (at the Iceland end of the FIR) produces significant lower-crustal reflectivity. Strong lower-crustal reflectivity elsewhere beneath the northwestern European continental margins may have a similar intrusive origin.