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.     

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.     

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.     

Neoproterozoic adakitic plutons in the northern margin of the Yangtze Block, China: Partial melting of a thickened lower crust and implications for secular crustal evolution

Numerous Neoproterozoic felsic and mafic–ultramafic intrusions occur in the Hannan region at the northern margin of the Yangtze Block. Among these, the Wudumen and Erliba plutons consist of granodiorites and have SHRIMP zircon U–Pb ages of  735 Ma. The rocks have high K₂O (0.8–3.6 wt.%) and Na₂O (4.4–6.4 wt.%) and low MgO (0.4–1.7 wt.%). They also have high Sr/Y (32–209) and (La/Yb)_n ratios (4.4–38.6). Their εNd values range from − 0.41 to − 0.92 and zircon initial ¹⁷⁶Hf/¹⁷⁷Hf ratios from 0.282353 to 0.282581. These geochemical features are similar to those of adakitic rocks produced by partial melting of a thickened lower crust. Our new analytical results, combined with the occurrence of voluminous arc-related mafic–ultramafic intrusions emplaced before 740 Ma, lead us to propose that the crustal evolution in the northern margin of the Yangtze Block during Neoproterozoic involved: (1) rapid crustal growth and thickening by underplating of mafic magmas from the mantle which was modified by materials coming from the subducting oceanic slab from  1.0 to  0.74 Ga, and (2) partial melting of the thickened lower crust due to a thermal anomaly induced by upwelling of asthenosphere through an oceanic slab window, producing the  735 Ma adakitic Wudumen and Erliba plutons. Our model suggests that the crustal thickness was more than 50 km at the northern margin of the Yangtze Block at  735 Ma, and rule out the possibility of a mantle plume impact causing the > 735 Ma magmatism in the region.     

Constraints on crustal and mantle structure of the oceanic plate south of Iceland from ocean bottom recorded Rayleigh waves

From April to July 2002 we carried out a deployment of 6 ocean bottom seismometers and 4 ocean bottom hydrophones in the North Atlantic south of Iceland. During the deployment period we recorded clear Rayleigh waves from 2 regional and 14 teleseismic earthquakes. This corresponds to a Rayleigh wave detection rate of nearly 92% for events with M_W ≥ 6.06.0 and epicentral distance less than 110°, close to detection rate estimates based on noise level variability. We measured Rayleigh wave event-station group dispersion and inter-station phase dispersion for one Mid-Atlantic Ridge event. The group dispersion curve is sensitive to the structure of the North-East Atlantic with an average age of  39 Myr. The phase dispersion curve is sensitive to the structure just south of Iceland (average plate age 33 Myr). Both dispersion curves indicate faster velocities than previously postulated for oceanic plate generated at the Reykjanes Ridge. A grid search approach was used to constrain the range of models fitting the data. The high velocity seismic lid just south of Iceland in the model for the phase dispersion path is slower or thinner than in the group dispersion model, which averages over a larger area and a somewhat older plate age, but the velocities in the low velocity half space are similar. We further consider the residual bathymetry in the experimental area. The residual anomaly decreases by 300–400 m from the Reykjanes Ridge to the  30 Myr old plate south of Iceland. This decrease can be explained by the disappearance of a mantle thermal anomaly associated with the Iceland plume. Both the residual bathymetry and the surface wave data are thus consistent with the notion that the southward spreading of the Icelandic plume is channelised underneath the Reykjanes Ridge and does not spread far outside this channel.Based on the experience from the pilot experiment, we estimate that a minimum recording time of 13–15 months in favourable weather conditions (April–September) is required to record enough data to map the spreading plume with surface waves, and to produce a tomographic image to a depth of 1000 km using body waves. This can be achieved by a continuous deployment of at least  20 months, or by two or three deployments during the spring and summer of consecutive years.     

Effect of water and stress on the lattice-preferred orientation of olivine

The influence of water and stress on the lattice-preferred orientation (LPO) of olivine aggregates was investigated through large strain, shear deformation experiments at high pressures and temperatures (P = 0.5–2.1 GPa, T = 1470–1570 K) under both water-poor and water-rich conditions. The specimens are hot-pressed synthetic olivine aggregates or single crystals of olivine. Water was supplied to the sample by decomposition of a mixture of talc and brucite. Deformation experiments were conducted up to γ (shear strain)  6 using the Griggs apparatus where water fugacity was up to  13 GPa at the pressure of 2 GPa. The water content in olivine saturated with water increases with increasing pressure and the solubility of water in olivine at P = 0.5–2 GPa was  400–1200 ppm H/Si. Several new types of LPO in olivine are found depending on water content and stress. Samples deformed in water-poor conditions show a conventional LPO of olivine where the olivine [100] axis is subparallel to the shear direction, the (010) plane subparallel to the shear plane (type-A). However, we identified three new types (type-B, C, and E) of LPO of olivine depending on the water content and stress. The type-B LPO of olivine which was found at relatively high stress and/or under moderate to high water content conditions is characterized by the olivine [001] axis subparallel to the shear direction, the (010) plane subparallel to the shear plane. The type-C LPO which was found at low stress and under water-rich conditions is characterized by the olivine [001] axis subparallel to the shear direction, the (100) plane subparallel to the shear plane. The type-E LPO which was found under low stress and moderate water content is characterized by the olivine [100] axis subparallel to the shear direction, the (001) plane subparallel to the shear plane. Observations by transmission electron microscopy (TEM) and scanning electron microscopy (SEM) show that the dislocations in water-poor samples (type-A) are curved and both b =  [100] and b = [001] dislocations have a similar population. Numerous subgrains are seen in water-poor samples in backscattered electron images. In contrast, water-rich samples (both type-B and type-C) contain mostly b = [001] dislocations and dislocations are straight and sub-grain boundaries are rare compared to those in water-poor samples. These observations suggest that (1) dominant slip systems in olivine change with water fugacity (and stress) and (2) grain boundary migration is enhanced in the presence of water. Seismic anisotropy corresponding to the fabrics under water-rich condition is significantly different from that under water-poor condition. Consequently, the relationship between seismic anisotropy and flow geometry in water-rich regions is expected to be different from that in water-poor regions in which type-A fabric dominates (i.e., the lithosphere). A few cases are discussed including anisotropy in the subduction zone and in the deep upper mantle.     

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.     

Anatolian surface wave evaluated at GEOFON Station ISP Isparta, Turkey

We have studied seismic surface waves of 255 shallow regional earthquakes recently recorded at GEOFON station ISP (Isparta, Turkey) and have selected these 52 recordings with high signal-to-noise ratio for further analysis. An attempt was made by the simultaneous use of the Rayleigh and Love surface wave data to interpret the planar crust and uppermost mantle velocity structure beneath the Anatolian plate using a differential least-square inversion technique. The shear-wave velocities near the surface show a gradational change from approximately 2.2 to 3.6 km s^− 1 in the depth range 0–10 km. The mid-crustal depth range indicating a weakly developed low velocity zone has shear-wave velocities around 3.55 km s^− 1. The Moho discontinuity characterizing the crust–mantle velocity transition appears somewhat gradual between the depth range  25–45 km. The surface waves approaching from the northern Anatolia are estimated to travel a crustal thickness of  33 km whilst those from the southwestern Anatolia and part of east Mediterranean Sea indicate a thicker crust at  37 km. The eastern Anatolia events traveled even thicker crust at  41 km. A low sub-Moho velocity is estimated at  4.27 km s^− 1, although consistent with other similar studies in the region. The current velocities are considerably slower than indicated by the Preliminary Reference Earth Model (PREM) in almost all depth ranges.     

Archean plate tectonics, rise of Proterozoic supercontinentality and onset of regional, episodic stagnant-lid behavior

Intense devolatilization and chemical-density differentiation attended late-stage accretion of the primitive Earth; it lessened after crystallization of a magma mush ocean during continued cooling. By 4.3Ga, shallow seas were present, so surface temperatures had fallen far below the 1300, 1120, and 950°C low-pressure solidi of peridotite, basalt, and granite, respectively. At temperatures less than about half their solidi, such materials existed as lithosphere in the near-surface Hadean realm. Stagnant-lid convection probably did not occur because massive heat transfer necessitated vigorous crust–mantle overturn in the early, hot Earth. Instead, bottom-up mantle convection, including voluminous plume ascent, efficiently rid the planet of heat, but lessened over time. Plate thickening and broadening is reflected in the post-Hadean rock record. Stages of geologic evolution included: (a) 4.5–4.4Ga, early, chaotic magma mush ocean overturn and ephemeral lithospheric platelets; (b) 4.4–2.7Ga, growth of oceanic and diminutive continental plates, obliterated by return mantle flow prior to 4.0Ga, but the latter enlarging and gradually accumulating as largely submarine, sutured, sialic crust-capped lithospheric collages; (c) 2.7–1.0Ga, progressive assembly of old shields and younger orogenic belts into supercratonal plates characterized by continental freeboard, sedimentary differentiation, and episodic glaciation during transpolar drift, as well as onset of regionally, temporally limited stagnant-lid convection beneath supercontinents; (d) 1.0Ga-present, modern, laminar-flowing asthenospheric cells capped by giant, stately moving plates. Restriction of komatiitic lavas to the Archean, and of multicycle sediments, most ophiolite complexes ± alkaline igneous rocks, and high-pressure and ultrahigh-pressure metamorphic belts to progressively younger Proterozoic–Phanerozoic orogens reflects increasingly negative buoyancy of the cooler oceanic lithosphere. Attending supercontinent assembly, density instabilities of thickening oceanic plates increasingly began to dominate overturn of the suboceanic mantle as cold, top-down convection. Scales and dynamics of hot asthenospheric upwelling versus lithospheric foundering and asthenospheric return flow (bottom-up versus top-down) changed gradually over geologic time in response to planetary thermal relaxation.     

Cretaceous–Tertiary geology of the Gangdese Arc in the Linzhou area, southern Tibet

Knowledge of the Cretaceous–Tertiary history of upper crustal shortening and magmatism in Tibet is fundamental to placing constraints on when and how the Tibetan plateau formed. In the Lhasa terrane of southern Tibet, the widely exposed angular unconformity beneath uppermost Cretaceous–lower Tertiary volcanic-bearing strata of the Linzizong Formation provides an excellent geologic and time marker to distinguish between deformation that occurred before vs. during the Indo-Asian collision. In the Linzhou area, located  30 km north of the city of Lhasa, a > 3-km-thick section of the Linzizong Formation lies unconformably on Cretaceous and older rocks that were shortened by both northward- and southward-verging structures during the Late Cretaceous. The Linzizong Formation dips northward in the footwall of a north-dipping thrust system that involves Triassic–Jurassic strata and a granite intrusion in the hanging wall. U–Pb zircon geochronologic studies show that the Linzizong Formation ranges in age from 69 Ma to at least 47 Ma and that the hanging wall granite intrusion crystallized at  52 Ma, coeval with dike emplacement into footwall Cretaceous strata. ⁴⁰Ar/³⁹Ar thermochronologic studies suggest slow cooling of the granite between 49 and 42 Ma, followed by an episode of accelerated cooling to upper crustal levels beginning at  42 Ma. The onset of rapid cooling was coeval with the cessation of voluminous arc magmatism in southern Tibet and is interpreted be a consequence of either (1) Tertiary thrusting in this region or (2) regional rock uplift and erosion following removal of overthickened Gangdese arc lower crust and upper mantle or break-off of the Neo-Tethyan oceanic slab.     

Evolution of the passive continental margins of India—a geophysical appraisal

The passive continental margins of India have evolved as India broke and drifted away from East Antarctica, Madagascar and Seychelles at various geological times. In this study, we have attempted to collate and re-examine gravity and topographic/bathymetry data over India and the adjoining oceans to understand the structure and tectonic evolution of these margins, including processes such as crustal/lithosphere extension, subsidence due to sedimentation, magmatic underplating and so on. The Eastern Continental Margin of India (ECMI) seems to have evolved in a complex rift and shear tectonic settings in its northern and southern segments, respectively, and bears similarities with its conjugate in East Antarctica. Crustal extension rates are uniform along the stretch of the ECMI in spite of the presence or absence of crustal underplated material, variability in lithospheric strength and tectonic style of evolution ranging from rifting to shearing. The Krishna-Godavari basin is underlain by a strong ( 30 km) elastic lithosphere, while the Cauvery basin is underlain by a thin elastic lithosphere ( 3 km). The coupling between the ocean and continent lithosphere along the rifted segment of the ECMI is across a stretched continental crust, while it is direct beneath the Cauvery basin. The Western Continental Margin of India (WCMI) seems to have developed in an oblique rift setting with a strike-slip component. Unlike the ECMI, the WCMI is in striking contrast with its conjugate in the eastern margin of Madagascar in respect of sedimentation processes and alignment of magnetic lineations and fracture zones. The break up between eastern India and East Antarctica seems to have been accommodated along a Proterozoic mobile belt, while that between western India and Madagascar is along a combination of both mobile belt and cratonic blocks.     

A very hydrous mantle under the western Pacific region: Implications for formation of marginal basins and style of Archean plate tectonics

The western Pacific region has been refrigerated by the subducting cold oceanic plates since 450 Ma. However, the region is also characterized by the presence of many oceanic microplates less than 1300 km across, as well as active magmatism; the Philippine Sea plate is representative. We have compiled and examined petrochemical characters of drilled basalts of DSDP from the Philippine Sea plate, and conclude that the source mantle for oceanic basalts is rich in water ca. 0.2 wt.%, and is 50–60 °C lower than that for MORB. The extensive melting is due to the high water content in the source mantle.It is well known that some marginal basins apparently have greater depths than the major oceans. We calculated the age–depth correlation based on a model of transient half-space cooling at given parameters of temperatures of mantle and surface, 1280 and 0 °C, and the thermal diffusivity, 1 mm² s^− 1. The calculation shows the correlation of age-residual depth from a mid-oceanic ridge is 367 for the Philippine Sea, consistent with the bathymetric data. Moreover, the mid-oceanic ridge may be relatively deep because this region is underlain by the cooler mantle.Addition of water to the mantle peridotite lowers the solidus temperature and viscosity. Melting experiments of hydrous peridotite show that addition of 0.2 wt.% H₂O content lowers the solidus temperature by 150 °C. As a result, the mantle under the region may practically correspond to a ca. 90 °C hotter mantle than normal MORB-source mantle in terms of magmatism and rheology. Numerical simulation for a hotter mantle suggests that many small plates should be formed because of extensive heat release by active magmatism, consistent with many microplates in this region. The presence of many oceanic microplates may be analogous to Archean plate tectonics, characterized by a hotter mantle.     

The structure, isostasy and gravity field of the Levant continental margin and the southeast Mediterranean area

A 3-D layered structure of the Levant and the southeastern Mediterranean lithospheric plates was constructed using interpretations of seismic measurements and borehole data. Structural maps of three principal interfaces, elevation, top basement and the Moho, were constructed for the area studied. This area includes the African, Sinai and Arabian plates, the Herodotus and the Levant marine basins and the Nile sedimentary cone. In addition, an isopach map of the Pliocene sediments, as well as the contemporaneous amount of denuded rock units, was prepared to enable setting up the structural map of the base Pliocene sediment. Variable density distributions are suggested for the sedimentary succession in accord with its composition and compaction. The spatial density distribution in the crystalline crust was calculated by weighting the thicknesses of the lower mafic and the upper felsic crustal layers, with densities of 2.9 g/cm³ and 2.77 g/cm³, respectively. Results of the local (Airy) isostatic modeling with compensation on the Moho interface show significant deviations from the local isostasy and require variable density distribution in the upper mantle. Moving the compensation level to the base of the lithosphere ( 100 km depth) and adopting density variations in the mantle lithosphere yielded isostatic compensation (± 200 m) over most of the area studied. The spatial pattern obtained of a density distribution with a range of ± 0.05 g/cm³ is supported by a regional heat flux. Simulations of the flexure (Vening Meinesz) isostasy related to the Pliocene to Recent sedimentary loading and unloading revealed concentric oscillatory negative and positive anomalies mostly related to the Nile sedimentary cone. Such anomalies may explain the rapid subsidence in the Levant Basin and the arching in central Israel, northern Sinai and Egypt during Pliocene–Recent times. Comparison between the observed (Bouguer) gravity and the calculated gravity for the constructed 3-D lithospheric structure, which has variable density distributions, provided a good match and an independent constraint for the large-scale structure suggested and confirmed an oceanic nature for the Levant Basin lithosphere.     

Constraining models of the tectonic setting of the giant Olympic Dam iron oxide–copper–gold deposit, South Australia, using deep seismic reflection data

In the Gawler Craton, the completeness of cover concealing the crystalline basement in the region of the giant Olympic Dam Cu–Au deposit has impeded any sufficient understanding of the crustal architecture and tectonic setting of its IOCG mineral-system. To circumvent this problem, deep seismic reflection data were recently acquired from  250 line-km of two intersecting traverses, centered on the Olympic Dam deposit. The data were recorded to 18 s TWT ( 55 km). The crust consists of Neoproterozoic cover, in places more than 5 km thick, over crystalline basement with the Moho at depths of 13–14 s TWT ( 40–42 km). The Olympic Dam deposit lies on the boundary between two distinct pieces of crust, one interpreted as the Archean–Paleoproterozoic core to the craton, the other as a Meso–Neoproterozoic mobile belt. The host to the deposit, a member of the  1590 Ma Hiltaba Suite of granites, is situated above a zone of reduced impedance contrast in the lower crust, which we interpret to be source-region for its  1000 °C magma. The crystalline basement is dominated by thrusts. This contrasts with widely held models for the tectonic setting of Olympic Dam, which predict extension associated with heat from the mantle producing the high temperatures required to generate the Hiltaba Suite granites implicated in mineralization. We use the seismic data to test four hypotheses for this heat-source: mantle underplating, a mantle-plume, lithospheric extension, and radioactive heating in the lower crust. We reject the first three hypotheses. The data cannot be used to reject or confirm the fourth hypothesis.     

Holocene climate variability and cultural evolution in the Near East from the Dead Sea sedimentary record

A comprehensive record of lake level changes in the Dead Sea has been reconstructed using multiple, well dated sediment cores recovered from the Dead Sea shore. Interpreting the lake level changes as monitors of precipitation in the Dead Sea drainage area and the regional eastern Mediterranean palaeoclimate, we document the presence of two major wet phases ( 10–8.6 and  5.6–3.5 cal kyr BP) and multiple abrupt arid events during the Holocene. The arid events in the Holocene Dead Sea appear to coincide with major breaks in the Near East cultural evolution (at  8.6, 8.2, 4.2, 3.5 cal kyr BP). Wetter periods are marked by the enlargement of smaller settlements and growth of farming communities in desert regions, suggesting a parallelism between climate and Near East cultural development.     

Cretaceous and Cenozoic cooling history across the ultrahigh pressure Tongbai–Dabie belt, central China, from apatite fission-track thermochronology

The crystalline terrane of the Tongbai–Dabie region, central China, comprising the Earth's largest ultrahigh-pressure (UHP) exposure was formed during Triassic collision between the Sino–Korean and Yangtze cratons. New apatite fission-track (AFT) data presented here from the UHP terrane, extends over a significantly greater area than reported in previous studies, and includes the (eastern) Dabie, the Hong'an (northwestern Dabie) and Tongbai regions. The new data yield ages ranging from 44 ± 3 to 142 ± 36 Ma and mean track lengths between 10 and 14.4 μm. Thermal history models based on the AFT data taken together with published ⁴⁰Ar/³⁹Ar, K–Ar, apatite and zircon (U–Th)/He and U–Pb data, exhibit a three-stage cooling pattern that is similar across the study region, commencing with an Early Cretaceous rapid cooling event, followed by a period of relative thermal stability during which rocks remained at temperatures within the AFT partial annealing zone (60–110 °C) and ending with a possible renewed phase of accelerated cooling during Pliocene to Recent time. The first cooling phase followed large-scale transtensional deformation between 140 and 110 Ma and is related to Early Cretaceous eastward tectonic escape and Pacific back arc extension. Between this phase and the subsequent slow cooling phase, a transition period from 120 to 80 Ma (to 70 to 45 Ma along the Tan–Lu fault) was characterised by a relatively low cooling rate (3–5 °C/Ma). This transition is likely related to a tectonic response associated with the mid-Cretaceous subduction of the Izanagi–Pacific plate as well as lithospheric extension and thinning in eastern Asia. The present regional AFT age pattern is therefore basically controlled by the Early Cretaceous rapid cooling event, but finally shaped through active Cenozoic faulting. Following the transition phase the subsequent slow cooling phase pattern implies a net reduction in horizontal compressional stress corresponding to increased extension rates along the continental margin due to the decrease in plate convergence. Modelling of the AFT data suggests a possible Pliocene–Recent cooling episode, which may be supported by increased rates of sedimentation observed in adjacent basins. This cooling phase may be interpreted as a response to the far-field effects of the frontal India–Eurasia collision to the west. Approximate estimates suggest that the total amount of post 120 Ma denudation across the UHP orogen ranged from 2.4 to 13.2 km for different tectonic blocks and ranged from 0.8 to 9.7 km during the Cretaceous to between 1.7 and 3.8 km during the Cenozoic.     

SHRIMP zircon U–Pb geochronological and whole-rock geochemical evidence for an early Neoproterozoic Sibaoan magmatic arc along the southeastern margin of the Yangtze Block

It has been generally accepted that the South China Block was formed through amalgamation of the Yangtze and Cathaysia Blocks during the Proterozoic Sibaoan orogenesis, but the timing and kinematics of the Sibao orogeny are still not well constrained. We report here SHRIMP U–Pb zircon geochronological and geochemical data for the Taohong and Xiqiu tonalite–granodiorite stocks from northeastern Zhejiang, southeastern margin of the Yangtze Block. Our data demonstrate that these rocks, dated at 913 ± 15 Ma and 905 ± 14 Ma, are typical amphibole-rich calc-alkaline granitoids formed in an active continental margin. Combined with previously reported isotopic dates for the  1.0 Ga ophiolites and  0.97 Ga adakitic rocks from northeastern Jiangxi, the timing of the Sibao orogenesis is thus believed to be between  1.0 and  0.9 Ga in its eastern segment. It is noted that the Sibao orogeny in South China is in general contemporaneous with some other early Neoproterozoic (1.0–0.9 Ga) orogenic belts such as the Eastern Ghats Belt of India and the Rayner Province in East Antarctica, indicating that the assembly of Rodinia was not finally completed until  0.9 Ga.     

Kinematics of the Iberia–Maghreb plate contact from seismic moment tensors and GPS observations

The Iberian Peninsula and the Maghreb experience moderate earthquake activity and oblique,  NW–SE convergence between Africa and Eurasia at a rate of  5 mm/yr. Coeval extension in the Alboran Basin and a N35°E trending band of active, left-lateral shear deformation in the Alboran–Betic region are not straightforward to understand in the context of regional shortening, and evidence complexity of deformation at the plate contact. We estimate 86 seismic moment tensors (M_W 3.3 to 6.9) from time domain inversion of near-regional waveforms in an intermediate period band. Those and previous moment tensors are used to describe regional faulting style and calculate average stress tensors. The solutions associated to the Trans-Alboran shear zone show predominantly strike-slip faulting, and indicate a clockwise rotation of the largest principal stress orientation compared to the regional convergence direction (σ₁ at N350°E). At the N-Algerian and SW-Iberian margins, reverse faulting solutions dominate, corresponding to N350°E and N310°E compression, respectively. Over most of the Betic range and intraplate Iberia, we observe predominately normal faulting, and WSW–ENE extension (σ₃ at N240°E). From GPS observations we estimate that more than 3 mm/yr of African (Nubian)–Eurasian plate convergence are currently accommodated at the N-Algerian margin,  2 mm/yr in the Moroccan Atlas, and  2 mm/yr at the SW-Iberian margin. 2 mm/yr is a reasonable estimate for convergence within the Alboran region, while Alboran extension can be quantified as  2.5 mm/yr along the stretching direction (N240°E). Superposition of both motions explains the observed left-lateral transtensional regime in the Trans-Alboran shear zone. Two potential driving mechanisms of differential motion of the Alboran–Betic–Gibraltar domain may coexist in the region: a secondary stress source other than plate convergence, related to regional-scale dynamic processes in the upper mantle of the Alboran region, as well as drag from the continental-scale motion of the Nubian plate along the southern limit of the region. In the Atlantic Ocean, the  3.5 mm/yr, westward motion of the Gibraltar Arc relative to intraplate Iberia can be accommodated at the transpressive SW-Iberian margin, while available GPS observations do not support an active subduction process in this area.     

Mantle flow-induced crustal thinning in the area between the easternmost part of the Anatolian plate and the Arabian Foreland (E Turkey) deduced from the geological and geophysical data

Eastern Anatolia consisting of an amalgamation of fragments of oceanic and continental lithosphere is a current active intercontinental contractional zone that is still being squeezed and shortened between the Arabian and Eurasian plates. This collisional and contractional zone is being accompanied by the tectonic escape of most of the Anatolian plate to the west by major strike-slip faulting on the right-lateral North Anatolian Transform Fault Zone (NATFZ) and left-lateral East Anatolian Transform Fault Zone (EATFZ) which meet at Karlıova forming an east-pointing cusp. The present-day crust in the area between the easternmost part of the Anatolian plate and the Arabian Foreland gets thinner from north (ca 44 km) to south (ca 36 km) relative to its eastern (EAHP) and western sides (central Anatolian region). This thinner crustal area is characterized by shallow CPD (12–16 km), very low Pn velocities (< 7.8 km/s) and high Sn attenuation which indicate partially molten to eroded mantle lid or occurrence of asthenospheric mantle beneath the crust. Northernmost margin of the Arabian Foreland in the south of the Bitlis–Pötürge metamorphic gap area is represented by moderate CPD (16–18 km) relative to its eastern and western sides, and low Pn velocities (8 km/s). We infer from the geophysical data that the lithospheric mantle gets thinner towards the Bitlis–Pötürge metamorphic gap area in the northern margin of the Arabian Foreland which has been most probably caused by mechanical removal of the lithospheric mantle during mantle invasion to the north following the slab breakoff beneath the Bitlis–Pötürge Suture Zone. Mantle flow-driven rapid extrusion and counterclockwise rotation of the Anatolian plate gave rise to stretching and hence crustal thinning in the area between the easternmost part of the Anatolian plate and the Arabian Foreland which is currently dominated by wrench tectonics.     

Insights into Wasatch fault vertical slip rates using the age of sediments in Timpanogos Cave, Utah

Timpanogos Cave, located near the Wasatch fault, is about 357 m above the American Fork River. Fluvial cave sediments and an interbedded carbonate flowstone yield a paleomagnetic and U–Th depositional age of 350 to 780 ka. Fault vertical slip rates, inferred from calculated river downcutting rates, range between 1.02 and 0.46 mm yr^− 1. These slip rates are in the range of the 0–12 Ma Wasatch Range exhumation rate ( 0.5–0.7 mm yr^− 1), suggesting that the long-term vertical slip rate remained stable through mid-Pleistocene time. However, the late Pleistocene (0–250 ka) decelerated slip rate ( 0.2–0.3 mm yr^− 1) and the accelerated Holocene slip rate ( 1.2 mm yr^− 1) are consistent with episodic fault activity. Assuming that the late Pleistocene vertical slip rate represents an episodic slowing of fault movement and the long-term (0–12 Ma) average vertical slip rate, including the late Pleistocene and Holocene, should be  0.6 mm yr^− 1, there is a net late Pleistocene vertical slip deficit of  50–75 m. The Holocene and late Pleistocene slip rates may be typical for episodes of accelerated and slowed fault movement, respectively. The calculated late Pleistocene slip deficit may mean that the current accelerated Wasatch fault slip rate will extend well into the future.