Evidence for lithospheric detachment in the central Andes from local earthquake tomography

Data from three temporary seismic networks were merged for tomographic inversion. Although the deployments did not coincide in time, spatial overlap was achieved by re-occupying existing sites. Travel times and t operators of about 1600 earthquakes were inverted for 3D models of ν_p, ν_p/ν_s and P-wave attenuation (Q_p^− 1). All three attributes provide a consistent image of the entire subduction zone on a lithospheric scale. The tomographic images reveal low velocities and high attenuation in the crust and mantle underlying the Western Cordillera and most of the Puna plateau, indicative of weak rheology and mostly asthenospheric mantle. In contrast, forearc and eastern foreland are characterized by high Q_p values, corresponding to cold temperatures in accordance with thermal models. In the backarc, between 23°S and 24°S, a high velocity, high Q_p structure beneath the Eastern Cordillera and eastern Puna is interpreted as detaching continental lithosphere that has been thickened in the orogenic process. South of this structure, the mantle is characterized by low velocities, high ν_p/ν_s ratios, and low Q_p values. Here it is believed that lithosphere originally underlying Andean crust has already been removed. This is supported by new estimates of crustal thickness and volcanic activity.     

Intrinsic and scattering attenuation in western India from aftershocks of the 26 January, 2001 Kachchh earthquake

176 vertical-component, short period observations from aftershocks of the M_w 7.7, 26 January, 2001 Kachchh earthquake are used to estimate seismic wave attenuation in western India using uniform and two layer models. The magnitudes (M_w) of the earthquakes are less than 4.5, with depths less than 46 km and hypocentral distances up to 110 km. The studied frequencies are between 1 and 30 Hz. Two seismic wave attenuation factors, intrinsic absorption (Q_i^− 1) and scattering attenuation (Q_s^− 1) are estimated using the Multiple Lapse Time Window method which compares time integrated seismic wave energies with synthetic coda wave envelopes for a multiple isotropic scattering model. We first assume spatial uniformity of Q_i^− 1, Q_s^− 1 and S wave velocity (β). A second approach extends the multiple scattering hypothesis to media consisting of several layers characterized by vertically varying scattering coefficient (g), intrinsic absorption strength (h), density of the media (ρ) and shear wave velocity structure. The predicted coda envelopes are computed using Monte Carlo simulation. Results show that, under the assumption of spatial uniformity, scattering attenuation is greater than intrinsic absorption only for the lowest frequency band (1 to 2 Hz), whereas intrinsic absorption is predominant in the attenuation process at higher frequencies (2 to 30 Hz). The values of Q obtained range from Q_t = 118, Q_i = 246 and Q_s = 227 at 1.5 Hz to Q_t ≈ 4000, Q_i ≈ 4600 and Q_s ≈ 33,300 at 28 Hz center frequencies, being Q_t^− 1 a measure of total attenuation. Results also show that Q_i^− 1, Q_s^− 1 and Q_t^− 1 decrease proportional to f^−ν. Two rates of decay are clearly observed for the low (1 to 6 Hz) and high (6 to 30 Hz) frequency ranges. Values of ν are estimated as 2.07 ± 0.05 and 0.44 ± 0.09 for total attenuation, 1.52 ± 0.21 and 0.48 ± 0.09 for intrinsic absorption and 3.63 ± 0.07 and 0.06 ± 0.08 for scattering attenuation for the low and high frequency ranges, respectively. Despite the lower resolution in deriving the attenuation parameters for a two layered crust, we find that scattering attenuation is comparable to or smaller than the intrinsic absorption in the crust whereas intrinsic absorption dominates in the mantle. Also, for a crustal layer of thickness 42 km, intrinsic absorption and scattering estimates in the crust are lower and greater than those of the mantle, respectively.     

Shear wave velocity and attenuation structure for the shallow crust of the southern Korean peninsula from short period Rayleigh waves

We analyzed the short period Rayleigh waves from the first crustal-scale seismic refraction experiment in the Korean peninsula, KCRUST2002, to determine the shear wave velocity and attenuation structure of the uppermost 1 km of the crust in different tectonic zones of the Korean peninsula and to examine if this can be related to the surface geology of the study area. The experiment was conducted with two large explosive sources along a 300-km long profile in 2002. The seismic traces, recorded on 170 vertical-component, 2-Hz portable seismometers, show distinct Rayleigh waves in the period range between 0.2 s and 1.2 s, which are easily recognizable up to 30–60 km from the sources. The seismic profiles, which traverse three tectonic regions (Gyeonggi massif, Okcheon fold belt and Yeongnam massif), were divided into five subsections based on tectonic boundaries as well as lithology. Group and phase velocities for the five subsections obtained by a continuous wavelet transform method and a slant stack method, respectively, were inverted for the shear wave models. We obtained shear wave velocity models up to a depth of 1.0 km. Overall, the shear wave velocity of the Okcheon fold belt is lower than that of the Gyeonggi and Yeongnam massifs by  0.4 km/s in the shallowmost 0.2 km and by 0.2 km/s at depths below 0.2 km. Attenuation coefficients, determined from the decay of the fundamental mode Rayleigh waves, were used to obtain the shear wave attenuation structures for three subsections (one for each of the three different tectonic regions). We obtained an average value of Q_β^− 1 in the upper 0.5 km for each subsection. Q_β^− 1 for the Okcheon fold belt ( 0.026) is approximately three times larger than Q_β^− 1 for the massif areas ( 0.008). The low shear wave velocity in the Okcheon fold belt is consistent with the high attenuation in this region.     

The anelasticity of the earth's mantle in the European area: frequency-dependence of Q

The frequency-dependence of Q is estimated for the earth's mantle beneath the Atlantic Ocean and the European continent. Employing an absorption band model, a minimum Q_β of 106 and 324 for the western and eastern part of Europe is found, resp. For both regions the upper cut-off frequency of the absorption band is estimated at 1.6 Hz.     

Compressional wave velocity of granite and amphibolite up to melting temperatures at 1 GPa

Compressional wave velocities (V_P) at above-solidus temperatures and at 1 GPa were obtained for a granite and amphibolite, which are considered to be major constituents of the continental crust. The temperature variation of velocities showed that the V_P values of granite decreased with rising temperature, but substantially increased beyond the melting temperature (850–900 °C). Such an increase may be caused by the α–β transition of quartz. The velocities of amphibolite decreased linearly with increasing temperature and dropped sharply at temperatures above the solidus (700 °C), indicating that partial melting of amphibolite acts to significantly lower the seismic velocities.     

Generalized geophysical model and dynamic properties of the continental crust

Detailed seismic investigations of the continental crust have produced evidence of definite regularities in the general layering of the consolidated crust despite its high degree of inhomogeneity. Three main layers may be resolved in the inner part of a continent: an upper layer with velocities of 5.8–6.4 km/s and a velocity gradient about 0.04–0.05 s⁻¹, an intermediate layer with velocities of 6.2–6.6 km/s and velocity gradient about zero, and a lower layer with velocities of 6.8–7.2 km/s and a high-velocity gradient of 0.05–0.1 s⁻¹. The intermediate layer is characteristically different not only because of its low average velocity gradient, but also because of its more pronounced horizontal layering, inversion zones, and its higher “transparency” and V_p/V_s ratio. The gravity and magnetic data have shown that basement inhomogeneities disappear at the top of the intermediate layer. Also there are few earthquakes in this layer. These pecularities may be interpreted as the result of partial melting (weakening) of rocks and their possible horizontal mobility inside this layer.Thus, dynamic models of tectonic processes must take into consideration the possible existence of a weak zone in the crust.     

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.     

Attenuation of shear-waves in the back-arc region of the Hellenic arc for frequencies from 0.6 to 16 Hz

Q_β for shear-waves is determined for the inner part of the Hellenic arc, the back-arc area, as a function of frequency in the range 0.6–16 Hz. We used 314 digital records from 32 earthquakes with magnitudes (M_w) ranging from 3.9 to 5.1. Epicentral distances ranged from 65 to 515 km. The data were obtained in 1997 during a 6-month operation of a digital portable network in Greece. The Q_β estimates were made for five frequency bands centred at 0.8, 1.5, 3.0, 6.0 and 12.0 Hz and the Q_β values obtained were 47, 79, 143, 271 and 553, respectively. The results show that Q_β for S-waves increases with frequency taking the form Q_β=55f ^0.91 (or Q_β⁻¹0.018f^−0.91). The high attenuation and the strong frequency dependence found, which is close to the frequency dependence of coda Q for Greece, are characteristic of an area with high seismicity, rapid extension, and in agreement with other similar studies in Greece.     

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.     

Crust–upper mantle seismic velocity structure across Southeastern China

One in-line wide-angle seismic profile was conducted in 1990 in the course of the Southeastern China Continental Dynamics project aimed at the study of the contact between the Cathaysia block and the Yangtze block. This 380-km-long profile extended in NW–SE direction from Tunxi, Anhui Province, to Wenzhou, Zhejiang Province. Five in-line shots were fired and recorded at seismic stations with spacing of about 3 km along the recording line. We have used two-dimensional ray tracing to model P- and S-wave arrivals and provide constraints on the velocity structure of the upper crust, middle crust, lower crust, Moho discontinuity, and the top part of the lithospheric mantle. P-wave velocity, S-wave velocity and V_P/V_S ratio are mapped. The crust is 36-km thick on average, albeit it gradually thins from the northwest end to the southeast end (offshore) of the profile. The average crustal velocity is 6.26 km/s for P-waves but 3.6 km/s for S-waves. A relatively narrow low-velocity layer of about 4 km of thickness, with P- and S-wave velocities of 6.2 km/s and 3.5 km/s, respectively, marks the bottom of the middle crust at a depth of 23-km northwest and 17-km southeast. At the crust–mantle transition, the P- and S-wave velocity change quickly from 7.4 to 7.8 km/s (northwest) and 8.0 to 8.2 km/s (southeast) and from 3.9 to 4.2 km/s (northwest) and 3.9 to 4.5 km/s (southeast), respectively. This result implies a lateral contrast in the upper mantle velocity along the 140 km sampled by the profile approximately. The average V_P/V_S ratio ranges from 1.68–1.8 for the upper crust to 1.75 for the middle and 1.75–1.85 for lower crust. With the interpretation of the wide-angle seismic data, Jiangshan–Shaoxin fault is considered as the boundary between the Yangtze and the Cathaysia block.     

Upper crustal seismic structure of the Mazury complex and Mazowsze massif within East European Craton in NE Poland

The large-scale seismic experiment POLONAISE '97 (POlish Lithospheric ONsets—An International Seismic Experiment) was carried out in May 1997 in Poland, Lithuania, and Germany. Its main purpose was to investigate the structure of the crust and the uppermost mantle in the region of the Trans European Suture Zone (TESZ) that lies between the East European Craton (EEC) and the Palaeozoic Platform. This paper covers the interpretation of seismic data along the NW–SE-trending, 180-km-long profile P5 located on the EEC. The recordings were of a high quality with seismic energy clearly visible along the whole profile. We have not found waves refracted below the upper crust in first arrivals. In the NW part of the profile, we have delineated a high-velocity body with the P-wave velocity in the range of 6.5–6.75 km/s in the upper crust. It corresponds to the K trzyn anorthosite massif within the Mazury complex. The Mazowsze massif is rather uniformly characterized by P-wave velocities 5.9–6.05 and 6.2–6.35 km/s in two layers, respectively. Sufficient S-wave data were available to estimate the V_p/V_s ratio (as well as the Poisson ratio), being 1.80 (0.277) in the high-velocity body and 1.67 (0.220) in the upper crust.Apart from the 2-D model along the profile, results of 3-D modelling in the area of the P5 profile are presented. Using off-line recordings, we got P-wave velocity field up to 8 km/s below the P5 profile at the depth of about 40 km as well as horizontal extent of the high-velocity body.     

Seismic attenuation of coda waves in the eastern region of Cuba

Cuba's seismic attenuation had never been studied in detail. In this paper we present the results of the research on the seismic attenuation of Cuba's eastern zone based upon the information collected by the seismological Cuban network from 1998 to 2003. 581 earthquakes were selected from the Cuban catalogue to make this study. All of them, recorded by at least three seismic stations, had their epicenters located in the eastern Cuban region (19.3–22 N, 79–73 W), epicentral distances between 15 km and 213 km, their coda duration magnitudes ranging from 2 and 4.1 and their focal depths reaching up to 30 km. The seismic wave attenuation was studied using coda waves. The single scattering method proposed by Sato in 1977 was applied, the attenuation and frequency dependency for different paths and the correlation of the results with the geotectonics of the region are presented in this paper.The mean Q_c value calculated was Q_c = (64 ± 2)f^{0.84 ± 0.01}. The relatively low Q₀ and the high frequency dependency agree with the values of a region characterized by a high tectonic activity. The Q_c values of seven subregions of eastern Cuba were calculated and correlated with the geology and tectonics of the area.     

Explosion seismic P and S velocity and attenuation constraints on the lower crust of the North–Central Tibetan Plateau, and comparison with the Tethyan Himalayas: Implications on composition, mineralogy, temperature, and tectonic evolution

P and S velocity and attenuation estimates in the lower crust are obtained from a set of wide angle reflection–refraction profiles in the region of active tectonics at the NE edge of the Tibetan Plateau and discussed together with respect to similar data at its Himalaya–south Tibet edge.The quality factor is estimated in the lower half of the crust by accounting for the differential effect on amplitude–frequency observed between waves of different penetrations, and both in P and S modes. Attenuation values allow to exclude a significant proportion of partial melt and to estimate the homologous temperature, ratio of in situ to solidus absolute temperatures. The latter depend on the physical conditions being of dry, wet or dehydration melting, which are found different among the regions of the northern Bayan Har and northern Qang Tang boundaries between blocks, as well as the Tethyan–Himalayas, south of the Indus–Tsangpo suture. Their in situ temperatures differ also as estimated from their different V_p for a similar felsic composition.Joint measurement of several parameters, V_p, V_s, Q_p and Q_s reveals the composition, mineralogy, temperature and hydration conditions of the lower half of the thickened crust of Tibet that may be discussed in terms of evolution. The material presently in the thickened crust, even its lower part, has a felsic composition, upper to middle crustal lithology, and the temperature conditions estimated suggest that basic material that could have underlain it could be eclogitized and not appear anymore above the seismic Moho.Under northern Qang Tang, the felsic material in the lower half of the crust appears as hot and dry. Its burial may have occurred earlier or may have been moderate in the postcollisional phase. This is consistent with a model of indentation of the Qang Tang crust by an originally thinner Bayan Har crust to bring part of its crust to greater depth, suggested from imaging the crustal architecture. Under northern Bayan Har, the material in the lower half of the crust appears as felsic, at low temperature and not dry conditions. This is evidence that it has been transported from a shallower depth, and this recently enough not to be yet dehydrated and temperature equilibrated in a conductive geotherm. It supports a model of recent overriding of the middle crust of the north Kun Lun block to the north independently suggested from the image of crustal architecture. The Tethyan Himalayas case appears bracketed by these two cases in northern Tibet for V_p and temperature conditions, but shows highest attenuation in the lower crust that is colder but less dry than under northern Qang Tang.     

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.     

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.     

Poisson's ratio in the oceanic crust — a review

Laboratory samples from the upper oceanic crust (tholeiitic basalt flows) that have not been significantly weathered, hydrothermally altered or fractured have a typical Poisson's ratio of 0.30 ( ) and a compressional velocity of 6.0 km s⁻¹; from the middle crust (dolerite sheeted dykes) a ratio of 0.28 ( ) and a velocity of 6.7 km s⁻¹; from the lower crust (gabbro) a ratio of 0.31 ( ) and a velocity of 7.1 km s⁻¹; and from the uppermost mantle a ratio of 0.24 ( ) and a velocity of 8.4 km s⁻¹. These sample values are representative of the large scale insitu values for the middle and lower crust and for the upper mantle. The upper crust is modified by several processes that decrease the velocity and generally increase Poisson's ratio: (1) the formation of an irregular layer of low temperature weathering generally less than 50 m thick; (2) large scale porosity in the form of drained pillows and lava tubes, of talus and rubble and of large open fractures; (3) where there was a high sedimentation rate over the ridge that formed the crust, hydrothermal alteration and intercalation of basalt and sediments. The Poisson's ratios of both high velocity sediments and of crystalline continental crustal rocks generally are significantly lower than the ratios of oceanic crustal rocks of similar compressional wave velocity. Thus, the use of shear wave velocities should permit the separation of these different formations which frequently cannot be distinguished on the basis of compressional wave seismic refraction data alone.     

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.     

Lithosphere structure of European sedimentary basins from regional three-dimensional gravity modelling

A three-dimensional (3D) density model, approximated by two regional layers—the sedimentary cover and the crystalline crust (offshore, a sea-water layer was added), has been constructed in 1° averaging for the whole European continent. The crustal model is based on simplified velocity model represented by structure maps for main seismic horizons—the “seismic” basement and the Moho boundary. Laterally varying average density is assumed inside the model layers. Residual gravity anomalies, obtained by subtraction of the crustal gravity effect from the observed field, characterize the density heterogeneities in the upper mantle. Mantle anomalies are shown to correlate with the upper mantle velocity inhomogeneities revealed from seismic tomography data and geothermal data. Considering the type of mantle anomaly, specific features of the evolution and type of isostatic compensation, the sedimentary basins in Europe may be related into some groups: deep sedimentary basins located in the East European Platform and its northern and eastern margins (Peri-Caspian, Dnieper–Donets, Barents Sea Basins, Fore–Ural Trough) with no significant mantle anomalies; basins located on the activated thin crust of Variscan Western Europe and Mediterranean area with negative mantle anomalies of −150 to −200×10⁻⁵ ms⁻² amplitude and the basins associated with suture zones at the western and southern margins of the East European Platform (Polish Trough, South Caspian Basin) characterized by positive mantle anomalies of 50–150×10⁻⁵ ms⁻² magnitude. An analysis of the main features of the lithosphere structure of the basins in Europe and type of the compensation has been carried out.     

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.     

Zircon U–Pb ages, Hf and O isotopes constrain the crustal architecture of the ultrahigh-pressure Dabie orogen in China

The crustal structure of the Dabie orogen was reconstructed by a combined study of U–Pb ages, Hf and O isotope compositions of zircons from granitic gneiss from North Dabie, the largest lithotectonic unit in the orogen. The results were deciphered from metamorphic history to protolith origin with respect to continental subduction and exhumation. Zircon U–Pb dating provides consistent ages of 751 ± 7 Ma for protolith crystallization, and two group ages of 213 ± 4 to 245 ± 17 Ma and 126 ± 4 to 131 ± 36 Ma for regional metamorphism. Majority of zircon Hf isotope analyses displays negative ε_Hf(t) values of − 5.1 to − 2.9 with crust Hf model ages of 1.84 to 1.99 Ga, indicating protolith origin from reworking of middle Paleoproterozoic crust. The remaining analyses exhibit positive ε_Hf(t) values of 5.3 to 14.5 with mantle Hf model ages of 0.74 to 1.11 Ga, suggesting prompt reworking of Late Mesoproterozoic to Early Neoproterozoic juvenile crust. Zircon O isotope analyses yield δ¹⁸O values of − 3.26 to 2.79‰, indicating differential involvement of meteoric water in protolith magma by remelting of hydrothermally altered low δ¹⁸O rocks. North Dabie shares the same age of Neoproterozoic low δ¹⁸O protolith with Central Dabie experiencing the Triassic UHP metamorphism, but it was significantly reworked at Early Cretaceous in association with contemporaneous magma emplacement. The Rodinia breakup at about 750 Ma would lead to not only the reworking of juvenile crust in an active rift zone for bimodal protolith of Central Dabie, but also reworking of ancient crust in an arc-continent collision zone for the North Dabie protolith. The spatial difference in the metamorphic age (Triassic vs. Cretaceous) between the northern and southern parts of North Dabie suggests intra-crustal detachment during the continental subduction. Furthermore, the Dabie orogen would have a three-layer structure prior to the Early Cretaceous magmatism: Central Dabie in the upper, North Dabie in the middle, and the source region of Cretaceous magmas in the lower.