Crustal velocity inhomogeneities along the Hirapur–Mandla profile, central India and its tectonic implications

Wide-angle seismic and gravity data across the Narmada-Son lineament (NSL) in central India are analyzed to determine crustal structure, velocity inhomogeneities and hence constrain the tectonics of the lineament. We present the 2-D crustal velocity structure from deep wide-angle reflection data by using a ray-trace inverse approach. The main result of the study is the delineation of fault-bounded horst raised to a subsurface depth (1.5 km) and the Moho upwarp beneath the NSL. The crust below the basement consists of three layers with velocities of 6.45–6.7, 6.2–6.5 and 6.7–6.95 km/s and interface depths of about 5.5–8.7, 14–17 and 18–23 km along the profile. The low-velocity (6.2–6.5 km/s) layer goes up to a depth of 5 km and becomes the thickest part (13 km), while the overlying high-velocity (6.45–6.7 km/s) layer becomes the thinnest (3 km) and upper boundary lies at a depth of 1.5 km beneath the NSL. The overall uncertainties of various velocity and boundary nodes are of the order of ±0.12 km/s and ±1.40 km, respectively. The up-lifted crustal block and the up-warping Moho beneath the NSL indicate that the north and south faults bounding the NSL are deeply penetrated through which mafic materials from upper mantle have been intruded into the upper crust. Gravity modeling was also undertaken to assess the seismically derived crustal features and to fill the seismic data gap. The lateral and vertical heterogeneous nature of the structure and velocity inhomogeneities in the crust cause instability to the crustal blocks and played an important role in reactivation of the Narmada south fault during the 1997 Jabalpur earthquake.     

Probing the history of the Mathematician paleoplate using surface waves

This paper investigates the shear velocity structure under the northern East Pacific Rise at the latitude range of 9–18°N, using intermediate-period Rayleigh and Love waves. The selected ocean-bottom seismic records provide source–receiver paths that ideally constrain the lithospheric mantle structure beneath the southern Rivera plate and the Mathematician paleoplate. The Rayleigh wave data infer a relatively thin ( 30 km) lithosphere under the eastern side of the present-day East Pacific Rise. The associated shear velocities are consistent with existing models of oceanic mantle beneath this region, and the estimated plate age of 2–3 million years agrees with results from magnetic dating. The west of the rise axis is characterized by a thicker and faster lithosphere than the eastern flank, and such structural differences suggest the presence of a relatively old Mathematician paleoplate. The discontinuous change in mantle structure across the East Pacific Rise spreading center are observed in both isotropic and anisotropic velocities. The young oceanic lithosphere east of the rise axis shows strong polarization anisotropy, where the dominant orientation of crystallographic axes roughly parallels the spreading direction. However, the western flank of the rise axis is approximately isotropic, and the lack of anisotropy suggests complex deformation mechanisms associated with earlier episodes of ridge segmentation, propagation and dual-spreading on and around the Mathematician paleoplate.     

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.     

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 splitting in the Eastern Alps observed at the TRANSALP network

Data from the temporary TRANSALP seismic network were analysed to investigate the seismic anisotropy in the upper mantle beneath the Eastern Alps. We operated mostly short period and some broadband stations in a dense linear array, which transects the Eastern Alps at 12° E longitude. Recorded SKS and SKKS phases with different backazimuths were used simultaneously to calculate the splitting parameters of delay time and fast axis direction for each seismic station. While we found variations in the delay times between 0.8 and 2.0 s, the determined fast axis directions prove to be rather consistent along the profile (60°–70° N). They coincide well with the trend of the Eastern Alps, thus suggesting orogen-parallel flow in the upper mantle. Our findings support the earlier proposed idea that the Adriatic indenter which was forced northwards into the European plate during the late stage of the Alpine orogeny, triggered an escape movement to the less constrained Pannonian basin to the east.     

Radial anisotropy in the crust and upper mantle beneath the Qinghai-Tibet Plateau and surrounding regions

Relative SV and SH wave speeds are generally attributed to radial seismic anisotropy which can be used as the indicator of crust/mantle deformation styles. Surface wave data were initially collected from events of magnitude Ms  5.0 and shallow or moderate focal depth occurring between 1980 and 2002: 713 of them generated Rayleigh waves and 660 Love waves, which were recorded by 13 broadband digital stations in Eurasia and India. Up to 1525 source-station Rayleigh waveforms and 1464 Love wave trains were earlier analysed by multiple filtering to obtain Love- and Rayleigh wave group velocity curves in the broad period range 10–105 s. We have performed tomographic inversion to obtain period-dependent group velocity and further shear wave velocity at 2° × 2°-sized grid-cells of a mesh covering the model region, after averaging azimuthal effects. Horizontally and vertically varying shear-wave velocities are observed, but the models of isotropic seismic velocity in the crust and upper mantle cannot fit simultaneously the inverted group-velocity dispersion curves due to the discrepancy in the transmission velocities of Love and Rayleigh waves, whose likely origin is the existence of radial anisotropy in the continental crust and topmost mantle. The strength of radial anisotropy computed from the Love–Rayleigh discrepancy and its spatial extent beneath the Qinghai-Tibet Plateau are shown as maps of percentage anisotropy at various depths down to 170 km and cross-sections along five profiles of reference. Areas in which radial anisotropy is in excess of 6% are found in the crust and upper mantle underlying most of the plateau, and even up to 10% in some places. The strength and spatial configuration of radial anisotropy seem to indicate the existence of a regime of horizontal compressive forces in the frame of the convergent Himalayan–Tibetan orogen, the laterally variation of the lithospheric rheology and the differential movement as regards the compressive driving forces.     

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.     

Constraints from Moho geometry and crustal thickness on the geodynamic origin of the Vrancea Seismogenic Zone (Romania)

Reprocessing of industry deep seismic reflection data (Ramnicu Sarat and Braila profiles) from the SE Carpathian foreland of Romania provides important new constraints on geodynamic models for the origin of the intermediate depth Vrancea Seismogenic Zone (VSZ). Mantle (70–200 km) earthquakes of the VSZ are characterized by high magnitudes (greater than 6.5), frequent occurrence rates (approximately 25 years), and confinement in a very narrow (30 × 70 × 200 km³) near vertical zone atypical for a Wadati–Benioff plane, located in front of the orogen. These two deep (20 s) seismic reflection profiles (70 km length across the foreland) reveal (1) a high-amplitude, gently east-dipping reflection across most of the section from what we interpret to be the Moho at  15 s (40–42 km) on the Ramnicu Sarat line to  16 s (47–48 km) on the Braila line, (2) a thick sedimentary cover increasing in thickness from east (1 s;  800 m) to west (7.5 s; 14 km), (3) an eastward increase in crustal thickness from 38 km (near VSZ) to  45 km, (4) seismic and topographic evidence for a newly imaged, possibly seismically active basement fault with a surface offset of 30 m observed on the Ramnicu Sarat line, (5) a lack of notable west-dipping structures in the crust and across the Moho, and (6) variable displacements on Peceneaga–Camena Fault of  5 km at Moho and  200 m at the basement–sedimentary cover contact.These observations appear to argue against recent models for west-dipping subduction of oceanic lithosphere at or in the vicinity of the Vrancea Seismogenic Zone given the lack of west-dipping fabrics in the lower crust and across the crust–mantle boundary. Consequently, one possible explanation for the geodynamic origin of VSZ could be partial delamination of the continental lithosphere in an intra-plate setting along a sub-horizontal lithospheric interface in the Carpathian hinterland that likely involves remnant lithospheric coupling between the crust and uppermost mantle in the foreland.     

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.     

Seismic properties of the upper crust in the central Friuli area (northeastern Italy) based on petrophysical data

The compressional and shear wave velocities have been measured at room temperature and pressure up to 450 MPa on 5 sedimentary rock samples, representative of the most common lithologies of the upper crust in the central Friuli area (northeastern Italy). At 400 MPa confining pressure the Triassic dolomitic rock shows the highest velocities (Vp  7 km/s, Vs  3.6 km/s), the Jurassic and Triassic limestones samples intermediate velocities (Vp  6.3 /s, Vs  3.5 km/s) and the Cenozoic and Paleozoic sandstones the lowest velocities (Vp  6.15 km/s, Vs  3.35 km/s). The Paleozoic sandstone sample is characterized by the strongest anisotropy (10%) and significant birefringence (0.2 km/s) is found only on the Cenozoic sandstone sample. We elaborated the synthetic profiles of seismic velocities, density, elastic parameters and reflection coefficient, related to 4 one-dimensional geological models extended up to 22 km depth. The synthetic profiles evidence high rheological contrasts between Triassic dolomitic rocks and the soft sandstones and the Jurassic limestones. The Vp profiles obtained from laboratory measurements match very well the in-situ Vp profile measured by sonic log for the limestones and dolomitic rocks, supporting our one-dimensional modelling of the calcareous-carbonatic stratigraphic series. The Vp and Vs values of the synthetic profiles are compared with the corresponding ones obtained from the 3-D tomographic inversion of local earthquakes. The laboratory Vp are generally higher than the tomographic ones with major discrepancies for the dolomitic lithology. The comparison with the depth location of seismicity reveals that the seismic energy is mainly released in correspondence of high-contrast rheological boundaries.     

3D seismic analysis of the Coast Shear Zone in SE Alaska and Western British Columbia: Broadside analysis of ACCRETE wide-angle data

The multidisciplinary ACCRETE project addresses the question of continental assemblage in southeast Alaska and western British Columbia by terrane accretion and magmatic addition. The previous studies of this project yielded important information for understanding the structures across the Coast Shear Zone (CSZ) and the formation of the CSZ and the Coast Mountains Batholith (CMB). The present study extends these interpretations into pseudo-3-D by using two additional wide-angle ACCRETE seismic lines. By analyzing the broadside wide-angle data using a series of laterally homogeneous 2-D models, we derive a lower-resolution 3-D velocity model of the outboard terranes and constrain variations in crustal thickness across and along the CSZ. Models of the broadside data confirms major structural and compositional trends extend along strike to the northwest. The key features are: a) a steep Moho ramp only  15-km wide is coincident with the CSZ and divides thin (25 ± 1 km) crust to the west below the west-vergent thrust belt (WTB) from thicker ( 31 ± 1 km) crust to the east below the CMB, (b) low-velocity mantle (7.7--7.9 km/s) extends beneath the entire study region indicating high crustal and upper-mantle temperatures below the WTB and CMB, and (c) the Alexander terrane is characterized by strong mid-crustal reflectivity and high lower crustal velocities that are consistent with gabbroic composition. This study extends the earlier interpretation and implies that the ramp is indeed likely associated with transpressional tectonics and magmatic crustal addition east of the CSZ.     

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.     

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.     

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.     

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.     

Evolution of mineralizing fluids in the Zn–Pb–Cu(–Ag ± Au) skarn and epithermal deposits of the world-class San Martín district, Zacatecas, Mexico

Mineralizing fluids at the San Martín skarn show an evolution characterized by prograde and retrograde associations. The prograde mineral associations consist of (1) a massive garnet zone, (2) a tremolite ± garnet zone, and (3) a late association of quartz, sphalerite, calcite and fluorite lining the vugs in the garnet zone. The fluids of the prograde associations exhibit decreasing temperatures of homogenization (Th) and variable salinities. The fluids of the massive garnet zone have salinities of 36 wt.% NaCl equiv. and Th of 645 to 570 °C, corresponding to pressures of 1055 bar. At the tremolite ± garnet zone, Th range from 438 to 354 °C. In the late association at the endoskarn, the following evolution can be drawn: (a) salinities of 50 to 42 wt.% NaCl equiv., and Th of 455 to 346 °C in quartz, (b) salinities of 46 wt.% NaCl equiv., and Th of 415 to 410 °C in sphalerite, (c) salinities of 50 to 37 wt.% NaCl equiv., and Th of 479 to 310 °C in calcite, (d) salinities of 33 to 28 wt.% NaCl equiv. and of 24 to 22 wt.% KCl in fluorite, and (e) two types of fluids with salinities of 2 and 39 wt.% NaCl equiv. and Th 344 and 300 °C, respectively, in later saccharoidal quartz segregations. The retrograde mineral associations comprise pervasive propylitic alteration to carbonization, and mantos with sulfides. Fluids in epidote have salinities of 7.6 wt.% NaCl equiv. and Th of 287 to 252 °C, and in calcite have salinities of 9.2 to 1 wt.% NaCl equiv. and Th of 188 to 112 °C. Fluids in the sulfide assemblages in the mantos have salinities of 8 to 3 wt.% NaCl equiv. and Th 300 °C, with corresponding pressures of 94 bar. Fluids in late epithermal veins close to the intrusive body have salinities of 10 to 5 wt.% NaCl equiv. and Th of 275 to 200 °C, and distal veins show salinities of 2 to 1 wt.% NaCl equiv. and Th of 160 °C.     

Experimental frictional heating of coal gouge at seismic slip rates: Evidence for devolatilization and thermal pressurization of gouge fluids

High velocity (1 m/s) friction experiments on bituminous coal gouge display several earthquake-related phenomena, including devolatilization by frictional heating, gas pressurization, and slip weakening. Stage I is characterized by sample shortening and reduction in the coefficient of friction (μ) from  1 to 0.6. Stage II is characterized by high frequency ( 5 Hz) oscillations in stress and strain records and by gas emissions. Stage III is marked by rapid weakening (μ  0.1 to 0.35) and sample shortening, together with continued gas emissions. Stage IV produces stable stress records and continued weakness (μ  0.2), but without gas emission. Stage I shortening is due to compaction of the gouge and the weakening is attributed to mechanical or thermal effects. Stage II behavior is interpreted as due to coal gasification and fluctuations in fluid pressure, resulting in high frequency stick-slip type behavior. Dramatic reduction in shear stress in stage III is attributed to gas pressurization by pore collapse and corresponds to a frictional instability, analogous to nucleation of an earthquake. Microstructural observations indicate the deformation was brittle during stages I and II but ductile during stages III and IV. Time dependent finite element frictional heat models indicate the center of the samples became hot ( 900 °C) during stage II, whereas the edge of samples remained relatively cold (< 300 °C). Vitrinite reflectance of coal samples shows an increase in reflectance from  0.5 to  0.8% over the displacement interval 20–40 m (20–40 s), indicating that the reflectance responds to frictional heating on a short time scale. The energy expended per unit area in these low stress, large displacement experiments is similar to that of higher stress ( 50 MPa), short displacement ( 1 m) earthquakes ( 10⁷ J/m²).     

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.     

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.