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.     

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.     

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.     

Global 1° × 1° thermal model TC1 for the continental lithosphere: Implications for lithosphere secular evolution

This paper reports a new 1° × 1° global thermal model for the continental lithosphere (TC1). Geotherms for continental terranes of different ages (> 3.6 Ga to present) constrained by reliable data on borehole heat flow measurements (Artemieva, I.M., Mooney, W.D. 2001. Thermal structure and evolution of Precambrian lithosphere: a global study. J. Geophys. Res 106, 16387–16414.), are statistically analyzed as a function of age and are used to estimate lithospheric temperatures in continental regions with no or low-quality heat flow data (ca. 60% of the continents). These data are supplemented by cratonic geotherms based on electromagnetic and xenolith data; the latter indicate the existence of Archean cratons with two characteristic thicknesses, ca. 200 and > 250 km. A map of tectono-thermal ages of lithospheric terranes complied for the continents on a 1° × 1° grid and combined with the statistical age relationship of continental geotherms (z = 0.04  t + 93.6, where z is lithospheric thermal thickness in km and t is age in Ma) formed the basis for a new global thermal model of the continental lithosphere (TC1). The TC1 model is presented by a set of maps, which show significant thermal heterogeneity within continental upper mantle, with the strongest lateral temperature variations (as large as 800 °C) in the shallow mantle. A map of the depth to a 550 °C isotherm (Curie isotherm for magnetite) in continental upper mantle is presented as a proxy to the thickness of the magnetic crust; the same map provides a rough estimate of elastic thickness of old (> 200 Ma) continental lithosphere, in which flexural rigidity is dominated by olivine rheology of the mantle.Statistical analysis of continental geotherms reveals that thick (> 250 km) lithosphere is restricted solely to young Archean terranes (3.0–2.6 Ga), while in old Archean cratons (3.6–3.0 Ga) lithospheric roots do not extend deeper than 200–220 km. It is proposed that the former were formed by tectonic stacking and underplating during paleocollision of continental nuclei; it is likely that such exceptionally thick lithospheric roots have a limited lateral extent and are restricted to paleoterrane boundaries. This conclusion is supported by an analysis of the growth rate of the lithosphere since the Archean, which does not reveal a peak in lithospheric volume at 2.7–2.6 Ga as expected from growth curves for juvenile crust.A pronounced peak in the rate of lithospheric growth (10–18 km³/year) at 2.1–1.7 Ga (as compared to 5–8 km³/year in the Archean) well correlates with a peak in the growth of juvenile crust and with a consequent global extraction of massif-type anorthosites. It is proposed that large-scale variations in lithospheric thickness at cratonic margins and at paleoterrane boundaries controlled anorogenic magmatism. In particular, mid-Proterozoic anorogenic magmatism at the cratonic margins was caused by edge-driven convection triggered by a fast growth of the lithospheric mantle at 2.1–1.7 Ga. Belts of anorogenic magmatism within cratonic interiors can be caused by a deflection of mantle heat by a locally thickened lithosphere at paleosutures and, thus, can be surface manifestations of exceptionally thick lithospheric roots. The present volume of continental lithosphere as estimated from the new global map of lithospheric thermal thickness is 27.8 (± 7.0) × 10⁹ km³ (excluding submerged terranes with continental crust); preserved continental crust comprises ca. 7.7 × 10⁹ km³. About 50% of the present continental lithosphere existed by 1.8 Ga.     

TOPO-OZ: Insights into the various modes of intraplate deformation in the Australian continent

As the fastest, lowest, flattest and amongst the most arid of continents, Australia preserves a unique geomorphic record of intraplate tectonic activity, evidencing at least three distinct modes of surface deformation since its rapid northward drift commenced around 43 million years ago. At long wavelengths (several 1000s km) systematic variations in the extent of Neogene marine inundation imply the continent has tilted north–down, southwest–up. At intermediate-wavelengths (several 100s km) several undulations of ~ 100–200 m amplitude have developed on the 1–10 myr timescale. At still shorter wavelengths (several 10s km), fault related motion has produced local relief at rates of up to ~ 100 m/myr over several million years. The long-wavelength, north–down tilting can be related to a dynamic topographic effect associated with Australia's northward drift from the geoid low, dynamic topography low now south of the continent to the geoid high, dynamic topography low centred above the south-east Asian and Melanesian subduction zones. The short wavelength, fault-related deformation is attributed in time to plate-wide increases in compressional stress levels as the result of distant plate boundary interactions and, in space, in part to variations in the thermal structure of the Australian lithosphere. At the intermediate wavelengths, transient, low amplitude undulations can be ascribed to either lithospheric buckling or the development of instabilities in the thermal boundary layer beneath the lithosphere. In the latter case, topographic asymmetries suggest the Australian lithosphere is moving north with respect to the mantle beneath, providing a unique attribution to the progressive alignment of seismic anisotropy and absolute plate motion observed near the base of the Australian lithosphere.     

中国东南地区岩石圈热结构特征及其构造意义

中国东南地区地质演化复杂,中—新生代构造变形强烈,岩石圈深部热力学状态及其对构造活动的影响有待深入。文章结合最新的大地热流数据与地壳结构Crust 1.0模型,利用稳态热传导方程,以岩石捕虏体温压数据和地震学观测为约束,构建了华南地区扬子克拉通、华夏地块以及南海北缘等不同单元的岩石圈热结构。结果表明该区岩石圈热结构存在强烈的不均一性：除了上扬子地区（四川盆地）为“温壳温幔”的热结构,华南其他大部分地区都表现为“热壳热幔”的特征;同一深度下,华夏地块与南海北缘的深部温度显著高于扬子克拉通;热岩石圈厚度从克拉通内部向沿海地区（NWSE）逐渐降低,也即由四川盆地的~200 km减少到华夏地块的~110 km,再到南海的~70 km。此外,我们还发现陆内地震的分布与岩石圈温度密切相关,地震活动集中分布于600℃等温线以内。总体而言,扬子克拉通中西部岩石圈热结构具有冷而厚的特征,而华夏地块和南海北缘受古太平洋平板俯冲和新生代大陆边缘构造—岩浆作用的改造,表现为热且薄的特征,岩石圈的热弱化进而加速了华南大陆边缘的裂解及随后的南海扩张过程。     

Lithospheric structure of the Southwest South China Sea: implications for rifting and extension

ABSTRACT

The South China Sea (SCS) is an excellent site for studying the process of conjugate margin rifting, and the origin and evolution of oceanic basins. Compared with the well-defined northern margin of the SCS, the western and southern segments of the SCS margin have not been researched in significant detail. To investigate the regional structure of the southwestern SCS, a gravity model is constructed, along with the lithospheric thermal structure along a wide-angle seismic profile. The profile extends across the conjugate margins of the Southwest Sub-Basin (SWSB) of the SCS and is based on the latest multiple geophysical measurements (including heat flow and thermo-physical parameters). The results show that the average thicknesses of the crust and thermal lithosphere along the profile are about 15 km and 57 km, respectively. The overall amount of extension of continental crust and lithosphere is more than 200 km. Thermal structure of the lithosphere shows that the continental margins are in a warm thermal state. The southwest SCS is characterized by ultra-wide, thinned continental crust and lithosphere, high Moho heat flow, early syn-rift faulted basins, undeformed late syn-rifting, and high seismic velocities in the lower crust. These various pieces of evidence suggest that the break-up of the mantle lithosphere occurred before that of the continental crust favouring a depth-dependent extension of the southwestern SCS margin.     

Magma genesis beneath Northeast Japan arc: A new perspective on subduction zone magmatism

It is being accepted that earthquakes in subducting slab are caused by dehydration reactions of hydrous minerals. In the context of this “dehydration embrittlement” hypothesis, we propose a new model to explain key features of subduction zone magmatism on the basis of hydrous phase relations in peridotite and basaltic systems determined by thermodynamic calculations and seismic structures of Northeast Japan arc revealed by latest seismic studies. The model predicts that partial melting of basaltic crust in the subducting slab is an inevitable consequence of subduction of hydrated oceanic lithosphere. Aqueous fluids released from the subducting slab also cause partial melting widely in mantle wedge from just above the subducting slab to just below overlying crust at volcanic front. Hydrous minerals in the mantle wedge are stable only in shallow (< 120 km) areas, and are absent in the layer that is dragged into deep mantle by the subducting slab. The position of volcanic front is not restricted by dehydration reactions in the subducting slab but is controlled by dynamics of mantle wedge flow, which governs the thermal structure and partial melting regime in the mantle wedge.     

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.     

Thermal impact of the break-up of Pangea on the Iberian Peninsula, assessed by thermochronological dating and numerical modelling

Thermochronological studies of Variscan basement in Iberia yield cooling ages typically younger than ~ 200 Ma. In this paper, we explore the regional implications of this recurrent age maximum by examination of low and high temperature thermochronological datasets from all over Iberia. Based on these results, we show that in general the lack of cooling ages older than 200 Ma is the result of several important regional periods of thermal resetting. Resetting took place in areas of extension and burial during the Mesozoic break-up of Pangea. Evidence for large scale magmatism and mineralisation is found in Iberia during the Mesozoic, since at that time Iberia formed part of the Central Atlantic Magmatic Province and a large mineralization province extending from North Africa to Western Europe. Numerical modelling allows us to assess the conditions under which rocks in the upper crust may have been thermally reset and the mechanisms likely involved. Results show that active rifting combined with shallow magmatism, and to a lesser extent deep sedimentary burial, could have led to an increase of the geothermal gradient up to ~ 73 °C/km and the reset of thermochronometers with closure temperatures up to 200 °C. Yet, we suggest that also hydrothermal activity, associated to extensional basins, played an important role to the increase of temperatures of some basement rocks above 300 °C.     

Seismic velocity structure of a large mafic intrusion in the crust of central Denmark from project ESTRID

The origin of regional sedimentary basins is being investigated by the ESTRID project (Explosion Seismic Transects around a Rift In Denmark). This project investigates the mechanisms of the formation of wide, regional basins and their interrelation to previous rifting processes in the Danish–Norwegian Basin in the North Sea region. In May 2004 a 143 km long refraction seismic profile was acquired along the strike direction of a suspected major mafic intrusion in the crust in central Denmark. The data confirms the presence of a body with high seismic velocity (> 6.5 km/s) extending from a depth of  10–12 km depth into the lower crust. There is a remarkable Moho relief between 27 and 34 km depth along this new along-strike profile as based on ray-tracing modelling of PmP reflections. The lack of PmP reflections at a zone of very high velocity in the lowest crust (7.3–7.5 km/s) suggests a possible location of a feeder channel to the batholith. The presence of volcanic rocks of Carboniferous–Permian age above the intrusion (mafic batholith) suggests a similar age of the intrusion. An older obliquely crossing profile and two new fan profiles deployed perpendicular to the main ESTRID profile, show that the batholith is about 30–40 km wide. The existence of this large mafic batholith supports the hypothesis that the origin of the Danish–Norwegian Basin is related to cooling and contraction after intrusion of large amounts of mafic melts into the crust during the late Carboniferous and early Permian. The data and interpretations from project ESTRID will form the basis for subsidence modelling. Tentatively, we interpret the formation of the Danish–Norwegian Basin as a thermal subsidence basin, which developed after widespread rifting of the region.     

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 review of the thermal regime of the Eastern Alps with respect to the effects of paleoclimate and exhumation

Transient thermal signals such as Pleistocene surface temperature variations or exhumation of great rock volumes are important for the current thermal regime of the Eastern Alpine crust. In this study transient 1-D forward simulations and an analytical approach were used to estimate the order of magnitude of these effects. A comparison with numerical forward simulations and inverse analyses of steady-state heat conduction yields the following main conclusions with respect to the thermal regime of the Eastern Alps along the TRANSALP profile: (1) The change of surface temperatures in the past affects mainly the uppermost part of the Eastern Alpine crust. It results in a maximum thermal signature of more than − 6 K at a depth of 2 km. The deviations from a steady-state temperature gradient and heat flow in the region of the Tauern Window range from 0.3–4 K km^− 1 and 0–6 mW m^− 2, respectively, with maximum values at the surface. (2) Exhumation of the Eastern Alpine lithosphere may result in a thermal signature of up to 4 K at a depth of 1 km. The thermal signature increases further with depth to a maximum of approximately 80 K at a depth of 50 km. As the temperature gradient of the exhumation signal is almost zero at the base of the crust, Moho heat flow appears to be not critically perturbed. (3) The combined effect of exhumation and changing surface temperatures at the Tauern Window amounts to less than 15% of the steady-state temperatures at a depth of  8 km and to less than 10% at the base of Eastern Alpine root. The corresponding perturbation in heat flow is less than 20% at a depth of 4 km, approaching zero below 40 km.     

Mountain building across a lithospheric boundary during arc construction: The Cretaceous Peninsular Ranges batholith in the Sierra San Pedro Martir of Baja California, Mexico

The Jura-Cretaceous Peninsular Ranges batholith (PRB) of Southern and Baja California contains a remarkable example of variation in crustal composition and structure across a batholith-parallel lithospheric-scale discontinuity. This lithospheric boundary between western oceanic-floored and eastern continental-floored crust influenced contractional deformation, arc magmatism, and differential exhumation of western and eastern zones in the batholith during its evolution.In the Sierra San Pedro Martir of Baja California, Mexico, a ca. 20 km wide, doubly vergent fan structure occurs across the PRB basement transition that consists of inward-dipping mylonite thrust sheets on the sides of the fan that gradually transition to a steeply-dipping tectonized zone in the center. A dramatic inverted metamorphic gradient occurs on the western side of this structure where mid-crustal amphibolite metamorphic grade rocks with peak pressures of 5–6 kbar in the center of the fan were thrust over upper-crustal sub-greenschist grade rocks (peak pressures < 2 kbar) in the western zone footwall. An inverted but smaller gradient occurs on the eastern side of the structure where rocks of the fan interior have been thrust eastwards over amphibolite to upper greenschist grade rocks (peak pressures 4–5 kbar).Gradients in cooling ages determined by ⁴⁰Ar/³⁹Ar biotite and K-feldspar and apatite fission track methods coupled with U–Pb zircon ages and Al in hornblende thermobarometry studies on plutons across this zone indicate that structures focused along the transition zone between contrasting lithosphere in the PRB accommodated nearly 15 km of the differential exhumation of western and eastern basement in the orogen. The western zone of the batholith was a major forearc depo-center for thick clastic sequences derived from the uplifting eastern PRB and remained at low average elevation during the Late Cretaceous and Paleogene. In contrast the eastern zone experienced dramatic uplift subsequent to achieving a crustal thickness in excess of 55 km by mid-Cretaceous time. This region had the isostatic potential for 4–5 km surface elevations, and likely formed a topographically high orogenic plateau. Exhumation of the fan structure initiated after 100 Ma and was largely complete by 85 Ma. Eastward-migrating unroofing of the rest of the eastern PRB continued into the Paleogene.A variety of factors were responsible for exhumation in this region. Structural thickening of the eastern zone of the orogen resulted from more than 30 million years of episodic contractional deformation in the fan structure, much of which followed island arc accretion of the western zone along this segment of the batholith. An episode of voluminous magmatism involving the intrusion of the 99–92 Ma La Posta-type magmatic suite across the eastern zone of the PRB triggered exhumation in the fan structure. Denudation in this region appears to have been solely by erosion; no evidence has been found for extensional tectonics during this time. This arc orogen demonstrates the important influence of inherited tectonic boundaries in controlling the spatial distribution of structural thickening and magmatism. It also displays the complex interrelationships among structural thickening, exhumation, and the role of magmatism in triggering exhumation episodes within orogens.     

Crustal structure of the northeastern margin of the Tibetan plateau from the Songpan-Ganzi terrane to the Ordos basin

The 1000-km-long Darlag–Lanzhou–Jingbian seismic refraction profile is located in the NE margin of the Tibetan plateau. This profile crosses the northern Songpan-Ganzi terrane, the Qinling-Qilian fold system, the Haiyuan arcuate tectonic region, and the stable Ordos basin. The P-wave and S-wave velocity structure and Poisson's ratios reveal many significant characteristics in the profile. The crustal thickness increases from northeast to southwest. The average crustal thickness observed increases from 42 km in the Ordos basin to 63 km in the Songpan-Ganzi terrane. The crust becomes obviously thicker south of the Haiyuan fault and beneath the West-Qinlin Shan. The crustal velocities have significant variations along the profile. The average P-wave velocities for the crystalline crust vary between 6.3 and 6.4 km/s. Beneath the Songpan-Ganzi terrane, West-Qinling Shan, and Haiyuan arcuate tectonic region P-wave velocities of 6.3 km/s are 0.15 km/s lower than the worldwide average of 6.45 km/s. North of the Kunlun fault, with exclusion of the Haiyuan arcuate tectonic region, the average P-wave velocity is 6.4 km/s and only 0.5 km/s lower than the worldwide average. A combination of the P-wave velocity and Poisson's ratio suggests that the crust is dominantly felsic in composition with an intermediate composition at the base. A mafic lower crust is absent in the NE margin of the Tibetan plateau from the Songpan-Ganzi terrane to the Ordos basin. There are low velocity zones in the West-Qinling Shan and the Haiyuan arcuate tectonic region. The low velocity zones have low S-wave velocities and high Poisson's ratios, so it is possible these zones are due to partial melting. The crust is divided into two layers, the upper and the lower crust, with crustal thickening mainly in the lower crust as the NE Tibetan plateau is approached. The results in the study show that the thickness of the lower crust increases from 22 to 38 km as the crustal thickness increases from 42 km in the Ordos basin to 63 km in the Songpan-Ganzi terrane south of the Kunlun fault. Both the Conrad discontinuity and Moho in the West-Qinling Shan and in the Haiyuan arcuate tectonic region are laminated interfaces, implying intense tectonic activity. The arcuate faults and large earthquakes in the Haiyuan arcuate tectonic region are the result of interaction between the Tibetan plateau and the Sino–Korean and Gobi Ala Shan platforms.     

Fossil flat-slab subduction beneath the Illinois basin, USA

The Illinois basin is one of several well-studied intracratonic sedimentary basins within the North American craton whose formational mechanisms and subcrustal structure are not well understood. We study the S-velocity structure of the upper mantle beneath the Illinois basin and its surrounding area through seismic tomography. We utilize continental scale waveform data of seismic S and surface waves, enhanced by regional earthquakes located near the Illinois basin. Our 3D tomographic model, IL05, confirms the existence of a slow S-velocity structure in the uppermost mantle beneath the Illinois basin region. This anomalously slow region exists from the base of the crust to depths of  90 km, and is slower than the North American cratonic average by about 200 m/s. This anomalous uppermost mantle beneath the Illinois basin is underlain by a faster lithosphere, typical of the surrounding craton, to depths of  200 km. Excluding the formation of the Reelfoot Rift, this area of North American has been stable for over 1.0 Gy. Thus, we do not expect thermal anomalies from before that time to persist into present day S-velocity anomalies and we consider a delamination origin as an explanation of Illinois basin subsidence unlikely. We cannot rule out that the slow mid-lithosphere beneath the Illinois basin is caused by an uppermost mantle enriched by a deep, but weak plume. We attribute the slow mid-lithosphere to the presence of either oceanic, hydrous crust, or, a relatively cool mantle wedge with preserved hydrous minerals in the Illinois basin's uppermost mantle, related to a fossilized flat subduction zone.     

甘肃龙首山岩带西井镁铁质岩体成因及其构造意义

西井岩体位于北祁连造山带以北,阿拉善地块西南缘的龙首山隆起带。以往的研究多以沿龙首山断裂分布的镁铁-超镁铁质岩带作为和金川岩体相关的岩浆事件进行,而本次选择西井镁铁质岩体进行了精确的地质年代学和地球化学研究,确定了西井岩体岩性主要为橄榄辉石岩和辉长岩,成岩时代为 (421.0±9.0) Ma,可以和北祁连高压变质带榴辉岩年龄相对应;ε_Nd(t)为4.06~5.52,(⁸⁷Sr/⁸⁶Sr)_i为0.704 548~0.707 575,具有地幔岩石圈特征;微量元素及其同位素计算表明西井岩体经历了约10%的下地壳物质混染。据此得出西井岩体及其龙首山岩带早志留世镁铁质侵入岩体成因模式为:祁连洋壳连续俯冲过程中洋壳与陆壳分离,热的软流圈物质持续冲击地幔岩石圈的底部;由于热传导效应,大地热流沿着地幔岩石圈上升,使得80 km深度的湿的橄榄岩层发生熔融,产生玄武质岩浆作用,玄武质岩浆上升过程中与下地壳物质发生约10%混染,形成西井岩体及其龙首山镁铁超镁铁质岩带。     

The structure of the lower crust at the Argentine continental margin, South Atlantic at 44°S

It is well established that the Argentine passive margin is of the rifted volcanic margin type. This classification is based primarily on the presence of a buried volcanic wedge beneath the continental slope, manifested by seismic data as a seaward dipping reflector sequence (SDRS). Here, we investigate the deep structure of the Argentine volcanic margin at 44°S over 200 km from the shelf to the deep oceanic Argentine Basin. We use wide-angle reflection/refraction seismic data to perform a joint travel time inversion for refracted and reflected travel times. The resulting P-wave velocity-depth model confirms the typical volcanic margin structure. An underplated body is resolved as distinctive high seismic velocity (v_p up to 7.5 km/s) feature in the lower crust in the prolongation of a seaward dipping reflector sequence. A remarkable result is that a second, isolated body of high seismic velocity (v_p up to 7.3 km/s) exists landward of the first high-velocity feature. The centres of both bodies are 60 km apart. The high-velocity lower-crustal bodies likely were emplaced during transient magmatic–volcanic events accompanying the late rifting and initial drifting stages. The lateral variability of the lower crust may be an expression of a multiple rifting process in the sense that the South Atlantic rift evolved by instantaneous breakup of longer continental margin segments. These segments are confined by transfer zones that acted as rift propagation barriers. A lower-crustal reflector was detected at 3 to 5 km above the modern Moho and probably represents the lower boundary of stretched continental crust. With this finding we suggest that the continent–ocean boundary is situated 70 km more seaward than in previous interpretations.     

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.     

冈底斯弧的岩浆作用：从新特提斯俯冲到印度亚洲碰撞

位于青藏高原南部的冈底斯岩浆弧形成于中生代新特提斯大洋岩石圈的长期俯冲过程中,而且在印度与亚洲大陆碰撞过程中叠加了强烈的新生代岩浆作用,是世界上典型的复合型大陆岩浆弧,已经成为研究汇聚板块边缘岩浆作用和大陆地壳生长与再造的天然实验室。基于对现有研究成果的总结,我们将冈底斯岩浆弧的岩浆构造演化划分为5个阶段：第1阶段发生在晚白垩世之前,以新特提斯洋岩石圈长期正常俯冲和钙碱性弧岩浆岩的发育为特征;第2阶段发生在晚白垩世时期,以活动的新特提斯洋中脊发生俯冲和强烈的岩浆作用与显著的新生地壳生长为特征;第3阶段发生在晚白垩世晚期,以残余的新特提斯大洋岩石圈俯冲和正常弧型岩浆作用为特征;第4阶段发生在古新世至中始新世,以印度与亚洲大陆碰撞、俯冲的新特提斯洋岩石圈回转和断离,及其诱发的幔源岩浆作用、新生和古老地壳的强烈再造为特征;第5阶段为发生在晚渐新世到中中新世的后碰撞阶段,深俯冲印度岩石圈的回转和断离,或加厚岩石圈地幔的对流移去导致了加厚下地壳的部分熔融和埃达克质岩石的广泛发育,同时伴随幔源钾质超钾质岩浆作用。冈底斯弧岩浆作用与岩浆成分的系统时空变化很好地记录了从新特提斯洋俯冲到印度亚洲大陆碰撞的完整构造演化过程。