Crustal accretion and dynamic feedback on mantle melting of a ridge centred plume: The Iceland case

A major consequence of the interaction of a plume with an oceanic ridge is the enhanced melt production and associated crust generation. In the case of Iceland crustal thickness as large as 20 to 40 km has been reported. Crustal seismic velocities are high, and have to be explained by thermal or chemical effects. In the first part of the paper we address the question whether extraction of melt out of the plume beneath a slowly spreading ridge and deposition of extracted basalt volumes at the surface produces a dynamic feedback mechanism on mantle melting. To study this question we solve the convection equations for a ridge centred plume with non-Newtonian rheology including melting, melt extraction associated with deposition of cold crust at the surface of the model, and using a simplified approach for compaction. The assumption of cold crust is justified if the thickness of each deposited basaltic layer is less than roughly 1 km. Depending on the buoyancy flux of the plume, crustal thicknesses between 10 and 40 km are modelled, showing characteristic dipping structures resembling the rift-ward dipping basaltic layers of East- and Western Iceland. Comparing the resulting crustal thickness and magma generation rate with models in which the dynamic effect of crust deposition has been suppressed indicates, that melt generation beneath a slowly spreading ridge is considerably damped by the dynamic feedback mechanism if the plume buoyancy flux exceeds 400 to 600 kg/s. Based on the observed crustal thickness of Iceland our models predict a plume buoyancy flux of 1140 kg/s.In the second part we study the accretion of the Icelandic crust by a thermo-mechanical model in more detail based on the Navier–Stokes-, the heat transport and the mass conservation equations including volumetric sources. Hot (1200 °C) molten crustal material is injected into the newly forming crust with a constant rate at different crustal source regions: a) deep, widespread emplacement of dykes and sills including crustal underplating, b) magma chambers at shallow to mid-crustal level, and c) surface extrusions and intrusions in fissure swarms at shallow depth connected to volcanic centres. We identify the material from the different source regions by a marker approach. Varying the relative dominance of these source regions, characteristic crustal structures evolve, showing shallow dipping upper crustal layers with dip angles between 10 and 15°. The thermal structure of the crust varies between cold crust (shallow-source region dominating) and hot crust (deep-source region dominating). We use observations of maximum depth of seismicity to constrain the depth of the 650 °C isotherm and seismological inferences on the lower crust to constrain temperatures in that region. The best agreement with our models is achieved for crust formation dominated by deep dykes and underplating with a considerable influence of magma chamber accretion.     

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

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

Influence of the precollisional stage on subduction dynamics and the buried crust thermal state: Insights from numerical simulations

At continental subduction initiation, the continental crust buoyancy may induce, first, a convergence slowdown, and second, a compressive stress increase that could lead to the forearc lithosphere rupture. Both processes could influence the slab surface P–T conditions, favoring on one side crust partial melting or on the opposite the formation of ultra-high pressure/low temperature (UHP-LT) mineral. We quantify these two effects by performing numerical simulations of subduction. Water transfers are computed as a function of slab dehydration/overlying mantle hydration reactions, and a strength decrease is imposed for hydrated mantle rocks. The model starts with an old oceanic plate ( 100 Ma) subducting for 145.5 Myr with a 5 cm/yr convergence rate. The arc lithosphere is thermally thinned between 100 km and 310 km away from the trench, due to small-scale convection occuring in the water-saturated mantle wedge. We test the influence of convergence slowdown by carrying on subduction with a decreased convergence rate (≤ 2 cm/yr). Surprisingly, the subduction slowdown yields not only a strong slab warming at great depth (> 80 km), but also a significant cooling of the forearc lithosphere at shallower depth. The convergence slowdown increases the subducted crust temperature at 90 km depth to 705 ± 62 °C, depending on the convergence rate reduction, and might thus favor the oceanic crust partial melting in presence of water. For subduction velocities ≤ 1 cm/yr, slab breakoff is triggered 20–32 Myr after slowdown onset, due to a drastic slab thermal weakening in the vicinity of the interplate plane base. At last, the rupture of the weakened forearc is simulated by imposing in the thinnest part of the overlying lithosphere a dipping weakness plane. For convergence with rates ≥ 1 cm/yr, the thinned forearc first shortens, then starts subducting along the slab surface. The forearc lithosphere subduction stops the slab surface warming by hot asthenosphere corner flow, and decreases in a first stage the slab surface temperature to 630 ± 20 °C at 80 km depth, in agreement with P–T range inferred from natural records of UHP-LT metamorphism. The subducted crust temperature is further reduced to 405 ± 10 °C for the crust directly buried below the subducting forearc. Such a cold thermal state at great depth has never been sampled in collision zones, suggesting that forearc subduction might not be always required to explain UHP-LT metamorphsim.     

The role of partial melting and extensional strain rates in the development of metamorphic core complexes

Orogenic collapse involves extension and thinning of thick and hot (partially molten) crust, leading to the formation of metamorphic core complexes (MCC) that are commonly cored by migmatite domes. Two-dimensional thermo-mechanical Ellipsis models evaluate the parameters that likely control the formation and evolution of MCC: the nature and geometry of the heterogeneity that localizes MCC, the presence/absence of a partially molten layer in the lower crust, and the rate of extension. When the localizing heterogeneity is a normal fault in the upper crust, the migmatite core remains in the footwall of the fault, resulting in an asymmetric MCC; if the localizing heterogeneity is point like region within the upper crust, the MCC remains symmetric throughout its development. Therefore, asymmetrically located migmatite domes likely reflect the dip of the original normal fault system that generated the MCC. Modeling of a severe viscosity drop owing to the presence of a partially molten layer, compared to a crust with no melt, demonstrates that the presence of melt slightly enhances upward advection of material and heat. Our experiments show that, when associated with boundary-driven extension, far-field horizontal extension provides space for the domes. Therefore, the buoyancy of migmatite cores contributes little to the outer envelope of metamorphic core complexes, although it may play a significant role in the internal dynamics of the partially molten layer. The presence of melt also favors heterogeneous bulk pure shear of the dome as opposed to the bulk simple shear, which dominates in melt-absent experiments. Melt presence affects the shape of P-T-t paths only slightly for material located near the top of the low-viscosity layer but leads to more complex flow paths for material inside the layer. The effect of extension rate is significant: at high extension rate (cm yr^− 1 in the core complex region), partially molten crust crystallizes and cools along a high geothermal gradient (35 to 65 °C km^− 1); material remains partially molten in the dome during ascent. At low strain rate (mm yr^− 1 in the core complex region), the partially molten crust crystallizes at high pressure; this material is subsequently deformed in the solid-state along a cooler geothermal gradient (20 to 35 °C km^− 1) during ascent. Therefore, the models predict distinct crystallization versus exhumation histories of migmatite cores as a function of extensional strain rates. The Shuswap metamorphic core complex (British Columbia, Canada) exemplifies a metamorphic core complex in which an asymmetric, detachment-controlled migmatite dome records rapid exhumation and cooling likely related to faster rates of extension. In contrast the Ruby Mountain-East Humboldt Ranges (Nevada, U.S.A.) exhibits characteristics associated with slower metamorphic core complexes.     

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.     

The state of the upper mantle beneath southern Africa

We present a new upper mantle seismic model for southern Africa based on the fitting of a large (3622 waveforms) multi-mode surface wave data set with propagation paths significantly shorter (≤ 6000 km) than those in globally-derived surface wave models. The seismic lithosphere beneath the cratonic region of southern Africa in this model is about 175 ± 25 km thick, consistent with other recent surface wave models, but significantly thinner than indicated by teleseismic body-wave tomography. We determine the in situ geotherm from kimberlite nodules from beneath the same region and find the thermal lithosphere model that best fits the nodule data has a mechanical boundary layer thickness of 186 km and a thermal lithosphere thickness of 204 km, in very good agreement with the seismic measurement. The shear wave velocity determined from analyzes of the kimberlite nodule compositions agree with the seismic shear wave velocity to a depth of 150 km. However, the shear wave velocity decrease at the base of the lid seen in the seismic model does not correspond to a change in mineralogy. Recent experimental studies of the shear wave velocity in olivine as a function of temperature and period of oscillation demonstrate that this wave speed decrease can result from grain boundary relaxation at high temperatures at the period of seismic waves. This decrease in velocity occurs where the mantle temperature is close to the melting temperature (within 100 °C).     

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

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

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

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

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.     

A new tomographic image of the Pyrenean lithosphere from teleseismic data

A new tomographic model of the Pyrenean lithosphere is determined down to 200 km depth from teleseismic P and PKP travel times, with a lateral resolution of 0.25°. Compared to previous models, two important improvements are 1) a larger number of stations with a more even distribution, in particular to the west of the range, and 2) the introduction, before inversion, of crustal corrections inferred from previous refraction and reflection experiments. This last point is crucial because a strong Moho jump (up to 20 km) is present at the North Pyrenean Fault, the former boundary between Eurasian and Iberian plates. The comparison of the models obtained with and without crustal corrections reveals the strong contamination of the models by the crust down to 100 km depth. In the uncorrected model, a large strip with negative P-velocity anomalies, previously interpreted as subduction of lower crust, is observed. It disappears in the corrected model. Moreover, the introduction of crustal corrections allows us to reveal short wavelength heterogeneities which were hidden by the crustal signal.An attempt is made to relate the heterogeneities revealed by the tomographic model with the tectonic history of the Pyrenees, in particular with the Alpine orogeny. The Alpine phase includes an extensive episode with generation of the thin continental crust and possibly the opening of an oceanic sea floor, and then a compressive stage. In our model, no signature of an oceanic subducted slab could be detected all along the range, a result which rules out the opening of a large oceanic floor before the compressive stage. A subduction of continental crust is possible but, due to the transformation of lower crust into eclogite at depth, it can not be detected by seismological methods, whereas it was observed from electrical and gravity data. To the East of the range, large heterogeneities with low velocities are ascribable to the Neogene extension related to the rotation of the Corso–Sardo block and the opening of the Gulf of Lion. A prominent high velocity anomaly extending down to 200 km in eastern-central Pyrenees could possibly be interpreted as a detached piece of the Tethys slab. In north of Iberia outside the range, deep (down to 200 km) low velocity structures oriented N130°E are probably related to Hercynian orogeny.     

Thermal evolution of the Indian shield and its subjacent mantle

The estimate of the paleothermal state (at −3 Ga) of the Indian shield is deduced using pressure and temperature estimates from the mineral assemblages combined with models of the heat transport mechanism and depth distribution of radiogenic heat in the crust. We find that at 3 Ga ago, the mean thermal gradient within the crust, reduced heat flow and mean mantle temperature were 28°C/km 1.08 h.f.u. and 2350°C, respectively.     

The crustal structure of the western Arabian platform from the spectral analysis of long-period P-wave amplitude ratios

The crustal structure of the western Arabian platform is derived using the spectral analysis of long-period P-wave amplitude ratios. The ratio of the vertical to the horizontal component is used to obtain the crustal transfer function based on thickness variations, crustal velocities, densities and the angle of emergence at the lower crust and upper mantle interface. Eleven well-defined earthquakes recorded at the long-period RYD station during the period from 1985 to 1994 were selected for analysis based on the following criteria: focal depths with a range between 7 and 89 km, body-wave magnitudes greater than 4.7, epicentral distances with a range from 8.8° to 26.5°, and back azimuthal coverage from 196° to 340°. Spectral analysis calculations were based on the comparison of the observed spectral ratios with those computed from theoretical P-wave motion obtained using the Thomson–Haskell matrix formulation for horizontally layered crustal models. The selection of the most suitable model was based on the identification of the theoretical model which exhibits the highest cross-correlation coefficient with the observed transfer function ratio. By comparing the spectral peak positions of the observed and theoretical values, the thickness and velocity can be resolved within 3 km and 1 km/s, respectively, of the observed values. The spectral analysis of long-period P-waves can detect a thin layer near the surface of about 1.6 km thick and a velocity contrast of about 10% with that of the underlying layer. A strong velocity gradient of about 0.05 km/s per km was found in the upper crust and 0.02 km/s per km in the lower crust. The derived crustal model is not unique due to the theoretical assumptions (horizontal layering, constant densities and velocities in each layer), quality of the data and complexities of the crustal structure. The crustal model suggests that the crust consists of five distinct layers. The upper crustal layer has a P-wave velocity of about 5.6 km/s and is about 1.6 km thick. The second layer has a velocity of about 6.2 km/s and is 10.2 km thick. The third layer shows a velocity of 6.6 km/s and is 6.8 km thick. The fourth layer has a velocity of about 6.8 km/s and is 12.3 km thick. The lower crustal layer has a velocity of about 7.5 km/s and is 9.3 km thick. The Mohorovicic discontinuity beneath the western Arabian platform indicates a velocity of 8.2 km/s of the upper mantle and 42 km depth.     

Crustal structure at the easternmost termination of the Variscan belt based on CELEBRATION 2000 and ALP 2002 data

The eastern margin of the Variscan belt in Europe comprises plate boundaries between continental blocks and terranes formed during different tectonic events. The crustal structure of that complicated area was studied using the data of the international refraction experiments CELEBRATION 2000 and ALP 2002. The seismic data were acquired along SW–NE oriented refraction and wide-angle reflection profiles CEL10 and ALP04 starting in the Eastern Alps, passing through the Moravo-Silesian zone of the Bohemian Massif and the Fore-Sudetic Monocline, and terminating in the TESZ in Poland. The data were interpreted by seismic tomographic inversion and by 2-D trial-and-error forward modelling of the P waves. Velocity models determine different types of the crust–mantle transition, reflecting variable crustal thickness and delimiting contacts of tectonic units in depth. In the Alpine area, few km thick LVZ with the Vp of 5.1 km s^− 1 dipping to the SW and outcropping at the surface represents the Molasse and Helvetic Flysch sediments overthrust by the Northern Calcareous Alps with higher velocities. In the Bohemian Massif, lower velocities in the range of 5.0–5.6 km s^− 1 down to a depth of 5 km might represent the SE termination of the Elbe Fault Zone. The Fore-Sudetic Monocline and the TESZ are covered by sediments with the velocities in the range of 3.6–5.5 km s^− 1 to the maximum depth of 15 km beneath the Mid-Polish Trough. The Moho in the Eastern Alps is dipping to the SW reaching the depth of 43–45 km. The lower crust at the eastern margin of the Bohemian Massif is characterized by elevated velocities and high Vp gradient, which seems to be a characteristic feature of the Moravo-Silesian. Slightly different properties in the Moravian and Silesian units might be attributed to varying distances of the profile from the Moldanubian Thrust front as well as a different type of contact of the Brunia with the Moldanubian and its northern root sector. The Moho beneath the Fore-Sudetic Monocline is the most pronounced and is interpreted as the first-order discontinuity at a depth of 30 km.     

An explosion seismology investigation of the continental margin west of the Hebrides, Scotland, at 58°N

During summer 1975, a line of large shots was fired across the continental margin between the Rockall Trough and the Hebridean shelf along 58°N. Arrivals were observed at temporary seismic stations set up across Scotland and in northwestern Ireland. A clear P₂ phase was observed to cross the margin and a converted phase P₁ also seen on the records is interpreted as travelling in the sub-sedimentary oceanic crust of Rockall Trough and in the upper continental crust beneath the shelf.The continental crust beneath the Hebridean shelf is estimated to be 27 ± 2 km thick, with P_g = 6.22 ± 0.03 km/s and P_n = 8.01 ± 0.04 km/s as determined by time-term analysis. P_g delays on the outer shelf are interpreted in terms of a seaward thickening wedge of Mesozoic sediments which pre-date the split. P_n beneath the Rockall Trough was poorly determined at 8.20 ± 0.17 km/s and the Moho is estimated to be 18 ± 2 km deep at 58°N. This and other seismic and gravity work indicates a northward thickening of the crust along the Rockall Trough, accounting for the northward decrease in the height of the slope.Our results, and those of gravity interpretations, indicate a relatively abrupt transition between continental and oceanic crust, possibly correlating with the lack of major shelf subsidence. This is attributed to a relatively cool origin for this margin. The main thinning of the continental crust beneath the slope is attributed to outslip of continental crustal material into and beneath the newly forming oceanic crust during the first few million years after the split, possibly enhanced by pre-split stretching.     

Heat flow and lithospheric thermal regime in the Northeast German Basin

New values of surface heat flow are reported for 13 deep borehole locations in the Northeast German Basin (NEGB) ranging from 68 to 91 mW m^− 2 with a mean of 77 ± 3 mW m^− 2. The values are derived from continuous temperature logs, measured thermal conductivity, and log-derived radiogenic heat production. The heat-flow values are supposed free of effects from surface palaeoclimatic temperature variations, from regional as well as local fluid flow and from thermal refraction in the vicinity of salt structures and thus represent unperturbed crustal heat flow. Two-D numerical lithospheric thermal models are developed for a 500 km section along the DEKORP-BASIN 9601 deep seismic line across the basin with a north-eastward extension across the Tornquist Zone. A detailed conceptual model of crustal structure and composition, thermal conductivity, and heat production distribution is developed. Different boundary conditions for the thickness of thermal lithosphere were used to fit surface heat flow. The best fit is achieved with a thickness of thermal lithosphere of about 75 km beneath the NEGB. This estimate is corroborated by seismological studies and somewhat less than typical for stabilized Phanerozoic lithosphere. Modelled Moho temperatures in the basin are about 800 °C; heat flow from the mantle is about 35 to 40 mW m^− 2. In the southernmost part of the section, beneath the Harz Mountains, higher Moho temperatures up to 900 to 1000 °C are shown. While the relatively high level of surface heat flow in the NEGB obviously is of longer wave length and related to lithosphere thickness, changes in crustal structure and composition are responsible for short-wave-length anomalies.     

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

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

Tomographic image of crustal structures across the Chelungpu fault: Is the seismogenic layer structure- or depth-dependent?

In 1999, a large earthquake (M_w = 7.6) occurred along the Chelungpu fault in the fold-and-thrust belt of western Taiwan. To shed more light on the subsurface structures and the seismogenic layers, three-dimensional velocity structures were inverted by using the travel times of both P- and S-waves from 2391 aftershocks recorded by the Central Weather Bureau during the 15 months that followed. From tomography, a typical image of the large-scale thrusting structures in the upper crust across the Chelungpu fault was obtained. In general, high velocities beneath the Western Foothills and Central Ranges are separated from low velocities beneath the Coastal Plain by an east-dipping boundary that is roughly consistent with the Chelungpu fault on the surface. The contrast in velocity on either side of the Chelungpu fault is indicative of about a 7- to 9-km vertical offset in the upper crust. The relocated hypocenter for the Chi-Chi earthquake shifts by 2.2 km toward the northwest, and its focal depth decreases by 0.7 km. A plot of focal depths versus rock velocities where the aftershocks occurred shows earthquakes are more inclined to occur in rock with a velocity of around 5.6 km/s. This strongly suggests the seismogenic layer in the fold-and-thrust belt of Taiwan is more structure-dependent than depth-dependent.     

Lateral and vertical variability in crustal velocity structure in a small area of the North Atlantic

The detailed velocity structure of a 70 by 35 km area of 6–10 Ma old crust on the flank of the mid Atlantic Ridge at 24°N was studied using 358 explosive charges and several hundred 16.4-1 airgun shots fired into an array of eight ocean bottom hydrophones. Inversion of the first arrival refracted travel times shows that the crust comprises a normal oceanic section about 5 km thick with a steep velocity gradient in the upper crust increasing from about 3.5 km/sec at the seafloor overlying a typical oceanic layer 3 and a probably anisotropic mantle. Delay time function mapping using two datasets containing arrivals from layer 3 and from the mantle show that lateral variability is generally low over most of the survey area, with a small region of high delay times in the northwest corner caused by the presence of abnormal crust probably associated with a minor fracture zone. We find that the topography of the base of layer 2 is similar to that of the top, indicating that the normal faulting which occurs along the margins of the median valley extends down at least into layer 3. Our observations from mantle arrivals are consistent with a much flatter Moho which constrains possible models of crustal formation at the spreading centre.     

The post-collisional stabilization of thickened crust in the Southern Urals

The Urals are characterized by a depression of the Moho to a depth of 57 km. This structure is interpreted as a relic orogenic root, which has been conserved because no significant post-collisional processes occurred. However, there is evidence that voluminous post-collisional magmatism affected the lower crust. In this paper, we use thermal finite element models to quantify the influence of the post-collisional magmatism on the stabilization of the root. We show that at least 70% of the heat producing elements migrated in granitic melts from the lower crust to the upper crust. As a result the crustal heat flow reduced and the lithosphere could stabilize at a thickness of 180 km. Furthermore, we propose that a granulite metamorphic event during the thermal relaxation of the collision zone prevented the 57 km thick crust from delamination. These results strongly indicate that post-collisional processes were necessary for the stabilization of the Uralian crust and lithosphere.     

Ultrahigh pressure macro diamonds from Copeton (New South Wales, Australia), based on Raman spectroscopy of inclusions

Mining of Cenozoic alluvial deposits at Copeton and Bingara (Eastern Australia) has produced two million macrodiamonds (0.25 ct median size). Raman spectroscopy is used to identify included minerals within uncut Copeton diamonds, with sealed chamber remnant pressures of 31.7 to 35.6 kbar for coesite, 13.6 and 22.7 kbar for clinopyroxene, and 7.6 kbar for grossular garnet. Assuming elastic behaviour, these values generate inclusion entrapment PT loci which intersect, restricting diamond formation conditions: from 250 °C, 43 kbar to 800 °C, 52 kbar. Larger than error (± 100 °C and ± 4 kbar), this range shows a systematic variation in inclusion composition with diamond zoning and N properties. Published research shows 1) Copeton and Bingara diamonds are unique, and 2) modern alluvium in the Bingara district carries mantle-formed garnet, captured by post-tectonic alkali basalt from an extensive diamondiferous ultrahigh pressure (UHP) terrane that stalled at depth because it is dominated by mafic eclogite. The combined Raman and geological results indicate two sets of subduction UHP diamond formation conditions/protolith are required, firstly cooler oceanic slab and secondly including higher temperature continental crust. The Copeton and Bingara stones are UHP macrodiamonds, and Carboniferous ⁴⁰Ar/³⁹Ar age dates on clinopyroxene inclusions should be interpreted as ages of crystallisation, representing the termination of subduction. The characteristic features of ruptured inclusions and etched percussion marks on Copeton and Bingara diamond indicate volcanic delivery to the earth's surface. Alluvial deposits elsewhere in Eastern Australia may carry similar diamond along with diamond of different origin.