Seismic velocities of granulite-facies xenoliths from Central Ireland: Implications for lower crustal composition and anisotropy

V_p and V_s values have been measured experimentally and calculated for granulite-facies lower crustal xenoliths from central Ireland close to the Caledonian Iapetus suture zone. The xenoliths are predominantly foliated and lineated metapelitic (garnet–sillimanite–K-feldspar) granulites. Their metapelitic composition is unusual compared with the mostly mafic composition of lower crustal xenoliths world-wide. Based on thermobarometry, the metapelitic xenoliths were entrained from depths of c. 20–25 ± 3.5 km and rare mafic granulites from depths of 31–33 ± 3.4 km. The xenoliths were emplaced during Lower Carboniferous volcanism and are considered to represent samples of the present day lower crust.V_p values for the metapelitic granulites range between 6.26 and 7.99 km s^{− 1} with a mean value of 7.09 ± 0.4 km s^{− 1}. Psammite and granitic orthogneiss samples have calculated V_p values of 6.51 and 6.23 km s^{− 1}, respectively. V_s values for the metapelites are between 3.86 and 4.34 km s^{− 1}, with a mean value of 4.1 ± 0.15 km s^{− 1}. The psammite and orthogneiss have calculated V_s values of 3.95 and 3.97 km s^{− 1}, respectively.The measured seismic velocities correlate with density and with modal mineralogy, especially the high content of sillimanite and garnet. V_p anisotropy is between 0.15% and 13.97%, and a clear compositional control is evident, mainly in relation to sillimanite abundance. Overall V_s anisotropy ranges from 1% to 11%. Poisson's ratio (σ) lies between 0.25 and 0.35 for the metapelitic granulites, mainly reflecting a high V_p value due to abundant sillimanite in the sample with the highest σ. Anisotropy is probably a function of deformation associated with the closure of the Iapetus ocean in the Silurian as well as later extension in the Devonian. The orientation of the bulk strain ellipsoid in the lower crust is difficult to constrain, but lineation is likely to be NE–SW, given the strike-slip nature of the late Caledonian and subsequent Acadian deformation.When corrected for present-day lower crustal temperature, the experimentally determined V_p values correspond well with velocities from the ICSSP, COOLE I and VARNET seismic refraction lines. Near the xenolith localities, the COOLE I line displays two lower crustal layers with in situ V_p values of 6.85–6.9 and 6.9–8.0 km s^{− 1}, respectively. The upper (lower velocity) layer corresponds well with the metapelitic granulite xenoliths while the lower (higher velocity) layer matches that of the basic granulite xenoliths, though their metamorphic pressures suggest derivation from depths corresponding to the present-day upper mantle.     

Vp/Vs ratio along the Vøring Margin, NE Atlantic, derived from OBS data: implications on lithology and stress field

A total of 13 regional Ocean Bottom Seismograph (OBS) profiles with an accumulated length of 2207 km acquired on the Vøring Margin, NE Atlantic have been travel time modelled with regards to S-waves. The V_p/V_s ratios are found to decrease with depth through the Tertiary layers, which is attributed to increased compaction and consolidation of the rocks. The V_p/V_s ratio in the intra-Campanian to mid-Campanian layer (1.75–1.8) in the central Vøring Basin is significantly lower than for the layers above and beneath, suggesting higher sand/shale ratio. This layer was confirmed by drilling to represent a layer of sandstone. This mid-Cretaceous ‘anomaly’ is also present in the northern Vøring Basin, as well as on the southern Lofoten Margin further north. The V_p/V_s ratio in the extrusive rocks on the Vøring Plateau is estimated to be 1.85, conformable with mafic (basaltic) rocks. Landward of the continent/ocean transition (COT), the V_p/V_s ratio in the layer beneath the volcanics is estimated to be 1.67–1.75. These low values suggest that this layer represents sedimentary rocks, and that the sand/shale ratio might be relatively high here. The V_p/V_s ratio in the crystalline basement is estimated to be 1.67–1.75 in the basin and on the landward part of the Vøring Plateau, indicating the presence of granitic/granodioritic continental crust. In the lower crust, the V_p/V_s ratio in the basin decreases uniformly from southwest to northeast, from 1.85–1.9 to 1.68–1.73, suggesting a gradual change from mafic (gabbroic) to felsic (granodioritic) lower crust. Significant (3–5%) azimuthal S-wave anisotropy is observed for several sedimentary layers, as well as in the lower crust. All these observations can be explained by invoking the presence of liquid-filled microcracks aligned vertically along the direction of the present day maximum compressive stress (NW–SE).     

Imaging of Vp, Vs, and Poisson’s ratio anomalies beneath Kyushu, southwest Japan: Implications for volcanism and forearc mantle wedge serpentinization

We determine detailed 3-D V_p and V_s structures of the crust and uppermost mantle beneath the Kyushu Island, southwest Japan, using a large number of arrival times from local earthquakes. From the obtained V_p and V_s models, we further calculate Poisson’s ratio images beneath the study area. By using this large data set, we successfully image the 3-D seismic velocity and Poisson’s ratio structures beneath Kyushu down to a depth of 150 km with a more reliable spatial resolution than previous studies. Our results show very clear low V_p and low V_s anomalies in the crust and uppermost mantle beneath the northern volcanoes, such as Abu, Kujyu and Unzen. Low-velocity anomalies are seen in the mantle beneath most other volcanoes. In contrast, there are no significant low-velocity anomalies in the crust or in the upper mantle between Aso and Kirishima. The subducting Philippine Sea slab is imaged generally as a high-velocity anomaly down to a depth of 150 km with some patches of normal to low seismic wave velocities. The Poisson’s ratio is almost normal beneath most volcanoes. The crustal seismicity is distributed in both the high- and low-velocity zones, but most distinctly in the low Poisson’s ratio zone. A high Poisson’s ratio region is found in the forearc crustal wedge above the slab in the junction area with Shikoku and Honshu; this high Poisson’s ratio could be caused by fluid-filled cracks induced by dehydration from the Philippine Sea slab. The Poisson’s ratio is normal to low in the forearc mantle in middle-south Kyushu. This is consistent with the absence of low-frequency tremors, and may indicate that dehydration from the subducting crust is not vigorous in this region.     

Crustal structure below the Polish Basin: Is it composed of proximal terranes derived from Baltica?

The large-scale seismic refraction and wide-angle reflection experiment POLONAISE'97 together with LT-7 and TTZ profiles carried out with the most modern techniques gave a high resolution of crustal structure of the Trans-European Suture Zone (TESZ) in NW and central Poland. The results of seismic investigations show the presence of relatively low velocity rocks (V_p < 6.1 km/s) down to a depth of 20 km beneath the Polish Basin (PB), and a high velocity lower crust (V_p = 6.8–7.3 km/s). The crustal thickness in the TESZ is intermediate between that of the East European Craton (EEC) to the northeast (40–45 km) and that of the Variscan crust (VB) to the southwest ( 30 km). Velocities in the uppermost mantle are relatively high (V_p = 8.25–8.45 km/s). The crust is three-layered with substantial differences in the velocities and thickness of individual layers. The area of the TESZ in NW and central Poland can be divided into at least two crustal blocks (terranes), called here Pomeranian Unit (PU, in the northwest) and Kuiavian Unit (KU, in the southeast). The postulated boundary between KU and PU is rather sharp at particular levels of the crust. Velocity distribution in the middle and lower crystalline crust in the TESZ area resemble values recognized in the EEC area, the fundamental difference being the much smaller thickness of both these layers. Our hypothesis/speculation is that the attenuated lower and middle crust of the TESZ belong to proximal terranes built of the EEC crust detached in the southeast and re-accreted to the EEC due to the process of anti-clockwise rotation of the Baltica paleocontinent during the Ordovician–Early Silurian.     

Isotope provinces, mechanisms of generation and sources of the continental crust in the Central Asian mobile belt: geological and isotopic evidence

The available geological, geochronological and isotopic data on the felsic magmatic and related rocks from South Siberia, Transbaikalia and Mongolia are summarized to improve our understanding of the mechanisms and processes of the Phanerozoic crustal growth in the Central Asian mobile belt (CAMB). The following isotope provinces have been recognised: ‘Precambrian’ (T_DM=3.3–2.9 and 2.5–0.9 Ga) at the microcontinental blocks, ‘Caledonian’ (T_DM=1.1–0.55 Ga), ‘Hercynian’ (T_DM=0.8–0.5 Ma) and ‘Indosinian’ (T_DM=0.3 Ga) that coincide with coeval tectonic zones and formed at 570–475, 420–320 and 310–220 Ma. Continental crust of the microcontinents is underlain by, or intermixed with, ‘juvenile’ crust as evidenced by its isotopic heterogeneity. The continental crust of the Caledonian, Hercynian and Indosinian provinces is isotopically homogeneous and was produced from respective juvenile sources with addition of old crustal material in the island arcs or active continental margin environments. The crustal growth in the CAMB had episodic character and important crust-forming events took place in the Phanerozoic. Formation of the CAMB was connected with break up of the Rodinia supercontinent in consequence of creation of the South-Pacific hot superplume. Intraplate magmatism preceding and accompanying permanently other magmatic activity in the CAMB was caused by influence of the long-term South-Pacific plume or the Asian plume damping since the Devonian.     

Seismological features of island arc crust as inferred from recent seismic expeditions in Japan

Crustal studies within the Japanese islands have provided important constraints on the physical properties and deformation styles of the island arc crust. The upper crust in the Japanese islands has a significant heterogeneity characterized by large velocity variation (5.5–6.1 km/s) and high seismic attenuation (Qp=100–400 for 5–15 Hz). The lateral velocity change sometimes occurs at major tectonic lines. In many cases of recent refraction/wide-angle reflection profiles, a “middle crust” with a velocity of 6.2–6.5 km/s is found in a depth range of 5–15 km. Most shallow microearthquakes are concentrated in the upper/middle crust. The velocity in the lower crust is estimated to be 6.6–7.0 km/s. The lower crust often involves a highly reflective zone with less seismicity, indicating its ductile rheology. The uppermost mantle is characterized by a low Pn velocity of 7.5–7.9 km/s. Several observations on PmP phase indicate that the Moho is not a sharp boundary with a distinct velocity contrast, but forms a transition zone from the upper mantle to the lower crust. Recent seismic reflection experiments revealed ongoing crustal deformations within the Japanese islands. A clear image of crustal delamination obtained for an arc–arc collision zone in central Hokkaido provides an important key for the evolution process from island arc to more felsic continental crust. In northern Honshu, a major fault system with listric geometry, which was formed by Miocene back arc spreading, was successfully mapped down to 12–15 km.     

3-D crustal structure of the extensional Granada Basin in the convergent boundary between the Eurasian and African plates

Three-dimensional P and S wave velocity models of the crust under the Granada Basin in Southern Spain are obtained with a spatial resolution of 5 km in the horizontal direction and 2 to 4 km in depth. We used a total of 15407 P and 13704 S wave high-quality arrival times from 2889 local earthquakes recorded by both permanent seismic networks and portable stations deployed in the area. The computed P and S wave velocities were used to obtain three-dimensional distributions of Poisson's ratio (σ) and the porosity parameter (V_p×V_s). The 3-D velocity images show strong lateral heterogeneities in the region. Significant velocity variations up to ±7% in P and S velocities are revealed in the crust below the Granada Basin. At shallow depth, high-velocity anomalies are generally associated with Mesozoic basement, while the low-velocity anomalies are related to the neogene sedimentary rocks. The south–southeastern part of the Granada Basin exhibits high σ values in the shallowest layers, which may be associated with saturated and unconsolidated sediments. In the same area, V_p×V_s is high outside the basin, indicating low porosity of the mesozoic basement. A low-velocity zone at 18-km depth is found and interpreted as a weak–ductile crust transition that is related to the cut-off depth of the seismic activity. In the lower crust, at 34-km depth, a clear slow V_p and V_s anomalous zone may indicate variations in lithology and/or with the rigidity of the lower crust rocks.     

Magma sources for Mesozoic anorogenic granites of the White Mountain magma series,New England,USA

The magma sources for granitic intrusions related to the Mesozoic White Mountain magma series in northern New England, USA, are addressed relying principally upon Nd isotopes. Many of these anorogenic complexes lack significant volumes of exposed mafic lithologies and have been suspected of representing crustal melts. Sm–Nd and Rb–Sr isotope systematics are used to evaluate magma sources for 18 felsic plutons with ages ranging from about 120 to 230 Ma. The possibility of crustal sources is further examined with analyses of representative older crust including Paleozoic granitoids which serve as probes of the lower crust in the region. Multiple samples from two representative intrusions are used to address intrapluton initial isotopic heterogeneities and document significant yet restricted variations (<1 in _Nd). Overall, Mesozoic granite plutons range in _Nd [T] from +4.2 to -2.3, with most +2 to 0, and in initial ⁸⁷Sr/⁸⁶Sr from 0.7031 to 0.709. The isotopic variations are roughly inversely correlated but are not obviously related to geologic, geographic, or age differences. Older igneous and metamorphic crust of the region has much lower Nd isotope ratios with the most radiogenic Paleozoic granitoid at _Nd [180 Ma] of -2.8. These data suggest mid-Proterozoic separation of the crust in central northern New England. Moreover, the bulk of the Mesozoic granites cannot be explained as crustal melts but must have large mantle components. The ranges of Nd and Sr isotopes are attributed to incorporation of crust by magmas derived from midly depleted mantle sources. Crustal input may reflect either magma mixing of crustal and mantle melts or crustal assimilation which is the favored interpretation. The results indicate production of anorogenic granites from mantle-derived mafic magmas.     

Crustal and mantle strengths in continental lithosphere: is the jelly sandwich model obsolete?

The relative importance of the contribution of the lower crust and of the lithospheric mantle to the total strength of the continental lithosphere is assessed systematically for realistic ranges of layer thickness, composition, and temperature. Results are presented as relative strength maps, giving the ratio of the lower crust to upper mantle contribution in terms of crustal thickness and surface heat flow. The lithosphere shows a “jelly sandwich” rheological layering for low surface heat flow, thin to average crustal thickness, and felsic or wet mafic lower crustal compositions. On the other hand, most of the total strength resides in the seismogenic crust in regions of high surface heat flow, crust of any thickness, and dry mafic lower crustal composition.     

青海格尔木市黑海地区三叠系巴颜喀拉山群变质碎屑岩地球化学特征及物源分析

巴颜喀拉地块位于青藏高原东北部,是古特提斯的主体,对其进行深入的研究不仅对确定古特提斯的演化过程与方式具有重要意义,而且对探讨该地区与古板块演化有关的成矿地质作用和矿产分布规律具有重要的现实意义。巴颜喀拉山群位于巴颜喀拉地块,对其物源研究较为薄弱。碎屑岩中的碎屑组分特征和地球化学特征可以反映物源区和沉积盆地的构造环境。本文通过对三叠系巴颜喀拉山群变碎屑岩岩相学特征、沉积结构特征、地球化学判别参数以及地球化学判别图解的研究,结合前人在这一地区所做的如古水流等相关工作,认为巴颜喀拉沉积盆地是在古老基底裂解的基础上形成的,存在裂陷扩张、海水加深的过程,具有大陆岛弧和活动大陆边缘型物源的特征,接受其北部东昆仑造山带碎屑的沉积,沉积物来自亲上地壳长英质源区。     

Crustal structure, fault segmentation, and activity of the Median Tectonic Line in Shikoku, Japan

We conducted a seismic tomographic analysis to estimate the crustal structure beneath the Shikoku and Chugoku regions in Japan. The Philippine Sea slab (PHS slab) subducts continuously in a SE–NW direction beneath this region, and the crustal structure is complex. Furthermore, the Median Tectonic Line (MTL), one of the longest and most active arc-parallel fault systems in Japan (hereafter, the MTL active fault system), is located in this area, and the right-lateral strike–slip movement of this fault system is related to the oblique subduction of the PHS slab. The MTL active fault system has ruptured repeatedly during the last 10 000 years, and has high seismic potential. Our tomographic analysis clarified the heterogeneous crustal structure along the MTL active fault system. This fault system in Shikoku can be divided into two segments, an east segment and a west segment, on the basis of the velocity structure. This segmentation model is consistent with other such models that have been determined from geological and geomorphological data such as fault geometry, slip rate, and faulting history. This consistency suggests that the surface characteristics of the MTL active fault system are related to structural properties of the crust. In particular, a prominent low-velocity (low-V) zone is present in the lower crust beneath the east segment. Our tomographic images show that the lower crust structure beneath the east segment is obviously different from that of the other segment. Furthermore, this low-V zone may indicate the presence of fluid, possibly related to dehydration of the PHS slab. As the presence of fluid in the lower crust affects the activity of the fault, stress accumulation and the fault failure mechanism may differ between the two segments of the MTL active fault system.     

Comparison of crustal velocity profiles determined by seismic refraction and teleseismic methods

Recently, two diverse seismic techniques were applied independently to the study of the crustal structure of the Cumberland Plateau, eastern Tennessee. One involved a reinterpretation of a refraction experiment performed in 1965 by the U.S. Geological Survey, consisting of two 400 km long, reversed refraction lines. The other entailed the inversion of broadband teleseismic P waveforms recorded at a single three-component broadband station, RSCP, located at the intersection of the two refraction profiles. A comparison of the two sets of velocity profiles revealed many similarities and some significant differences. Both sets of velocity models consist of three major crustal layers: (1) an upper crust (V_p = 6.1–6.4 km/s) down to about 17 km, (2) a mid-crust (V_p = 6.7–6.9 km/s) between 17 and 40 km depth, (3) a lower crust (V_p = 7.2–7.4 km/s) from 40 to 51 km depth. The refraction models have linear transition zones up to 11 km thick at the base of each layer, whereas the teleseismic models have more irregular transition zones at the base of the mid- and lower crust. The differences in the results of these studies are attributed to the differing frequency bandwidths of the data sets; the predominant sensitivity of the teleseismic data to shear velocities, compared to compressional velocities for the refraction data; and the different analysis procedures involved in each method. Nevertheless, the similarities indicate that the teleseismic waveform method with broadband data is capable of retreiving comparable crustal information as the Cumberland Plateau refraction survey. In addition, it provides the kind of complementary information required to constrain the composition of the continental lower crust and uppermost mantle.     

Lithosphere structure of the Ukrainian Shield and Pripyat Trough in the region of EUROBRIDGE-97 (Ukraine and Belarus) from gravity modelling

The southern segment of the seismic profile EUROBRIDGE—EUROBRIDGE-97 (EB'97)—located in Belarus and Ukraine, crosses the suture zone between two main segments of the East European Craton—Fennoscandia and Sarmatia—as well as Sarmatia itself. At the initial stage of our study, a 3-D density model has been constructed for the crust of the study region, including the major part of the Osnitsa–Mikashevichi Igneous Belt (OMIB) superimposed by sediments of the Pripyat Trough (PT), and three domains in the Ukrainian Shield—the Volhyn Domain (VD) with the anorthosite–rapakivi Korosten Pluton (KP), the Podolian Domain (PD), and the Ros–Tikich Domain (RTD). The model comprises three layers—sediments with maximum thickness (6 km) in the PT and two heterogeneous layers in the crystalline crust separated at a depth of 15 km. 3-D calculations show the main features of the observed gravity field are caused by density heterogeneities in the upper crust. Allocation of density domains deeper than 15 km is influenced by Moho topography. Fitting the densities here reveals an increase (up to 2960 kg m⁻³) in the modelled bodies accompanied by a Moho deepening to 50 km. In contrast, a Moho uplift to a level of 35–37 km below the KP and major part of the PT is associated with domains of reduced densities. An important role for the deep Odessa–Gomel tectonic zone, dividing the crust into two regions one of basically Archean consolidation in the west (PD and RTD) and one of Proterozoic crust in the east (Kirovograd Domain)—was confirmed.2-D density modelling on the EB'97 profile shows that in the upper crust three main domains of different Precambrian evolution—the OMIB (with the superimposed PT), the VD with the KP, and the PD—can be distinguished. Deeper, in the middle and lower crust, layered structures having no connection to the surface geology are dominant features of the models. Least thickness of the crust was obtained below the KP. Greatest crustal thickness (more than 50 km) was found below the PD, characterised also by maximum deviation of velocity/density relation in the rocks from a standard one. The velocity and density models along the EB'97 profile have been interpreted together with inferred V_p/V_s ratios to estimate crustal composition in terms of SiO₂ content. In the course of the modelling, the status of the PD as a centre of Archean granulitic consolidation has been confirmed. The crustal structure of the anorthosite–rapakivi KP is complex. For the first time, a complicated structure for the lower crust and lower crust–upper mantle transition zone beneath the KP has been determined. The peculiarities of the crustal structure of the KP are quite well explained in terms of formation of rapakivi–anorthosite massifs as originating from melt chambers in the upper mantle and lower crust. An important role for the South Pripyat Fault (SPF), repeatedly activated during Proterozoic–Palaeozoic times, has been ascertained. At the subplatform stage of crustal evolution the SPF was, probably, a magma channel facilitating the granitic intrusions of the KP. In the Palaeozoic the fault was reactivated during rifting in the PT.     

Crustal thickness and VP/VS variations in the Grenville orogen (Ontario, Canada) from analysis of teleseismic receiver functions

We have developed a simple semblance-weighted stacking technique to estimate crustal thickness and average V_P/V_S ratio using teleseismic receiver functions. We have applied our method to data from 32 broadband seismograph stations that cover a 700 × 400 km² region of the Grenville orogen, a 1.2–0.98 Ga Himalayan-scale collisional belt in eastern North America. Our seismograph network partly overlaps with Lithoprobe and other crustal refraction surveys. In 8 out of 9 cases where a crustal-refraction profile passes within 30 km of a seismograph station, the two independent crustal thickness estimates agree to within 7%. Our regional crustal-thickness model, constructed using both teleseismic and refraction observations, ranges between 34.0 and 52.4 km. Crustal-thickness trends show a strong correlation with geological belts, but do not correlate with surface topography and are far in excess of relief required to maintain local isostatic equilibrium. The thickest crust (52.4 ± 1.7 km) was found at a station located within the 1.1 Ga mid-continent (failed) rift. The Central Gneiss Belt, which contains rocks exhumed from deep levels of the crust, is characterized by V_P/V_S ranging from 1.78 to 1.85. In other parts of the Grenville orogen, V_P/V_S is found to be generally less than 1.80. The thinnest crust (34.5–37.0 km) occurs northeast of the 0.7 Ga Ottawa–Bonnechere graben and correlates with areas of high intraplate seismicity.     

2-D seismic tomographic and ray tracing modelling of the crustal structure across the Sudetes Mountains basing on SUDETES 2003 experiment data

The seismic data obtained during SUDETES 2003 experiment are analysed, and detailed crustal structure for profiles S02, S03 and S06 is presented using three different 2-D techniques: (1) “smooth” tomography of refracted waves travel times, (2) ray tracing of reflected and refracted waves, and (3) joint velocity and depth of reflector tomographic inversion. In spite of different interpretation techniques used, the models of the crustal structure show common characteristic features. The low velocity (V_p < 4 km/s) sedimentary layer was documented in the northeastern part of the study area. The topmost basement has in general a velocity of 5.8–6.0 km/s, and velocities at ca. 20 km depth are 6.15–6.25 km/s. The strong reflecting boundaries were found at 20–23 and 25–28 km depth with a velocity contrast about 0.4 km/s, and the highest velocities in the lowermost crust are 6.8–7.2 km/s. In general, the crust of the Bohemian Massif is slightly thicker (33–35 km) than in the northern part of the area. Velocities beneath Moho are relatively low, of 7.95 km/s. On the basis of well recorded reflected waves, mantle reflectors were discovered in the depth interval ca. 40–70 km. Apart of new results for the geology and tectonics of the area, some conclusion could be made about different techniques used. In the 2-D case the “clasical” ray tracing method with using all correlated phases gives the most adequate model of the structure, because of full, manual control of the model creation. The “smooth” first arrival travel times tomography, although very fast, is not satisfactory enough to describe the complex structure. So, the best candidate in 3-D case seems to be travel time tomography for both refracted and reflected waves in multi-layers models.     

The lower crust of the Gangdese magmatic arc,southern Tibet,implication for the growth of continental crust

Magmatic arcs are thought to be the primary sites of modern-day continental crustal growth, and arc crustal sections provide an exceptional opportunity to directly observe the geological processes that occur there, yet few deeply exposed arc sections are available for direct study. The Gangdese magmatic arc, southern Tibet, formed during the Mesozoic subduction of Neo-Tethyan oceanic lithosphere and Cenozoic collision between the Indian and Asian continents, and represent juvenile continental crust. However, the petrological components and compositions of the lower crust of the Gangdese arc remain unknown. Based on detailed geological mapping, we conducted a systemic geochemical, geochronological and zircon Hf isotopic study of well-exposed high-grade metamorphic and migmatitic rocks from the lower crust of the eastern Gangdese arc. The results obtained show that Late Cretaceous garnet amphibolites, dioritic and granitic gneisses, and Paleocene–Eocene garnet amphibolites and granitic gneisses are the main components of the Gangdese lower arc crust. These meta-intrusive rocks witnessed a long period of magmatic, and metamorphic and anatectic processes from the Middle Jurassic to the Late Eocene, and have chemical compositions that range from ultramafic to felsic, with an average SiO₂ content of 57.61 wt% and Mg# value of 0.49. These new data indicate firstly that the Gangdese lower arc crust has an overall intermediate composition and typical feature of juvenile crusts, and therefore supports the recent proposition that continental lower crusts are relatively felsic in composition, instead of mafic. We consider that the downward transport of felsic intrusives and associated sedimentary rocks into the deep crustal levels and subsequent partial melting resulted in componential and compositional changes of the Gangdese arc lower crust over time. This is a potential key mechanism in transforming primary lower arc crust to mature continental lower crust for the magmatic arcs with a complete growth history.     

Upper mantle structure of the Northern Eurasia from peaceful nuclear explosion data

Several long-range seismic profiles were carried out in Russia with Peaceful Nuclear Explosions (PNE). The data from 25 PNEs recorded along these profiles were used to compile a 3-D upper mantle velocity model for the central part of the Northern Eurasia. 2-D crust and upper mantle models were also constructed for all profiles using a common methodology for wavefield interpretation. Five basic boundaries were traced over the study area: N1 boundary (velocity level, V = 8.35 km/s; depth interval, D = 60–130 km), N2 (V = 8.4 km/s; D = 100–140 km), L (V = 8.5 km/s; D = 180–240 km) and H (V = 8.6 km/s; D = 300–330 km) and structural maps were compiled for each boundary. Together these boundaries describe a 3-D upper mantle model for northern Eurasia. A map characterised the velocity distribution in the uppermost mantle down to a depth of 60 km is also presented. Mostly horizontal inhomogeneity is observed in the uppermost mantle, and the velocities range from the average 8.0–8.1 km/s to 8.3–8.4 km/s in some blocks of the Siberian Craton. At a depth of 100–200 km, the local high velocity blocks disappear and only three large anomalies are observed: lower velocities in West Siberia and higher velocities in the East-European platform and in the central part of the Siberian Craton. In contrast, the depths to the H boundary are greater beneath the craton and lower beneath in the West Siberian Platform. A correlation between tectonics, geophysical fields and crustal structure is observed. In general, the old and cold cratons have higher velocities in the mantle than the young platforms with higher heat flows.Structural peculiarities of the upper mantle are difficult to describe in form of classical lithosphere–asthenosphere system. The asthenosphere cannot be traced from the seismic data; in contrary the lithosphere is suggested to be rheologically stratified. All the lithospheric boundaries are not simple discontinuities, they are heterogeneous (thin layering) zones which generate multiphase reflections. Many of them may be a result of fluids concentrated at some critical P–T conditions which produce rheologically weak zones. The most visible rheological variations are observed at depths of around 100 and 250 km.     

The Anfeg post-collisional Pan-African high-K calc-alkaline batholith (Central Hoggar, Algeria), result of the LATEA microcontinent metacratonization

The Anfeg batholith (or composite laccolith) occupies a large surface (2000 km²) at the northern tip of the Laouni terrane, just south of Tamanrasset in Hoggar. It is granodioritic to granitic in composition and comprises abundant enclaves that are either mafic microgranular enclaves (MME) or gneissic xenoliths. It intruded an Eburnian (≈2 Ga) high-grade basement belonging to the LATEA metacraton at approximately 608 Ma (recalculated from the U–Pb dating of [Tectonics 5 (1986) 955]) and cooled at approximately 4 kbar, with a temperature of about 750 °C. This emplacement occurred mainly along subhorizontal thrust planes related to Pan-African subvertical mega-shear zones close to the attachment zone of a strike-slip partitioned transpression system. Although affected by some LILE mobility, the Anfeg batholith can be ascribed to a high-K calc-alkaline suite but characterized by low heavy REE contents and high LREE/HREE ratios. The MME belong to the Anfeg magmatic trend while some xenoliths belong to Neoproterozoic island arc rocks.The Anfeg batholith defines a Nd–Sr isotopic initial ratios trend (_Nd/(⁸⁷Sr/⁸⁶Sr)_i from −2.8/0.7068 to −11.8/0.7111) pointing to a mixing between a depleted mantle and an old Rb-depleted granulitic lower crust. Both sources have been identified within LATEA and elsewhere in the Tuareg shield (_Nd/⁸⁷Sr/⁸⁶Sr)_i of +6.2/0.7028 for the depleted mantle, −22/0.708 for the old lower crust.The model proposed relates the above geochemical features to a lithospheric delamination along the subvertical mega-shear zones that dissected the rigid LATEA former passive margin without major crustal thickening (metacratonization) during the general northward tectonic escape of the Tuareg terranes, a consequence of the collision with the West African craton. This delamination allowed the uprise of the asthenosphere. In turn, this induced the melting of the asthenosphere by adiabatic pressure release and of the old felsic and mafic lower crust due to the high heat flow. A gradient in the mantle/crust ratio within the source of the Pan-African magmatism is observed in LATEA from the northeast (Egéré-Aleksod terrane) where rare plutons are rooted within the Archaean/Eburnian basement to the southwest (Laouni terrane) where abundant batholiths, including Anfeg, have a mixed signature. Some mantle melts with only slight crustal contamination (Laouni troctolitic layered intrusions) are even present. This suggests that the southern boundary of LATEA microcontinent is not far south of the Tuareg shield.     

Crustal structure in the Changbaishan volcanic area, China, determined by modeling receiver functions

During 1998–1999, we installed a temporary broadband seismic network in the Changbaishan volcanic region, NE China. We estimated crustal structure using teleseismic seismograms collected at the network. We detected a near surface region of strong anisotropy directly under the main volcanic edifice of the volcanic area. We modeled 109 receiver functions from 19 broadband stations using three techniques. First we used a “slant-stacking” method to model the principal crustal P reverberation phases to estimate crustal thickness and the average crustal P to S speed ratio (v_p/v_s), assuming an average P-wave velocity in the crust. We then estimated crustal S-wave velocity (v_s) and v_p/v_s profiles by modeling stacked receiver functions using a direct search. Finally, we inverted several receiver functions recorded at stations closest to the main volcanic edifice using least squares to estimate v_s velocity profiles, assuming a v_p/v_s value. The results from the three estimation techniques were consistent, and generally we found that the receiver functions constrained estimates of changes in wave speeds better than absolute values. We resolved that the crust is 30–39 km thick under the volcanic region and 28–32 km thick away from the volcanic region, with a midcrust velocity transition at about 10–15 km depth. We estimated that the average crust P-wave velocity is about 6.0–6.2 km/s surrounding the main volcanic region, while it is slightly lower in the vicinity of the main volcanic edifice. The estimates of v_p/v_s were more ambiguous, but we inferred that the bulk crustal Poisson's ratio (which is related to v_p/v_s) ranges between 0.20 and 0.30, with a suggestion that the Poisson's ratio is lower under the central volcanic region compared to the surrounding areas. We resolved low S-wave velocities (down to about 3 km/s) in the middle crust in the region of the main volcanic edifice. The low velocity anomaly extends from about 5–10 to 15–25 km below the surface, probably indicating a region of elevated temperatures. We were unable to determine if partial melt is present with the data we considered in this paper.     

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.