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.     

Brittle–ductile–brittle deformation during cooling of tonalite (Adamello, Southern Italian Alps)

Late- to post-magmatic deformation in slightly diachronous contiguous intrusions of the north-western Adamello batholith (Southern Alps, Italy) is recorded as, from oldest to youngest: (i) joints, (ii) solid-state ductile shear zones, (iii) faults associated with epidote-K-feldspar veins and (iv) zeolite veins and faults. Structures (ii) to (iv) are localized on the pervasive precursory network of joints (i), which developed during the earliest stages of pluton cooling. High temperature ( 500 °C), ductile overprinting of joints produced lineations, defined by aligned biotite and hornblende, on the joint surfaces and highly localized mylonites. The main phase of faulting, producing cataclasites and pseudotachylytes, occurred at  250 °C and was associated with extensive fluid infiltration. Cataclasites and pseudotachylytes are clustered along different E–W-striking dextral strike-slip fault zones correlated with the activity of the Tonale fault, a major tectonic structure that bounds the Adamello batholith to the north. Ductile deformation and cataclastic/veining episodes occurred at P = 0.25–0.3 GPa during rapid cooling of the batholith to the ambient temperatures ( 250 °C) that preceded the exhumation of the batholith. Timing of the sequence of deformation can be constrained by ³⁹Ar–⁴⁰Ar ages of  30 Ma on pseudotachylytes and various existing mineral ages. In the whole composite Adamello batholith, multiple magma pulses were intruded over the time span 42–30 Ma and each intrusive body shows the same ductile-to-brittle structural sequence localized on the early joint sets. This deformation sequence of the Adamello might be typical of intrusions undergoing cooling at depths close to the brittle–ductile transition.     

Frictional sliding of gabbro gouge under hydrothermal conditions

We investigated the frictional sliding behaviour of gabbro gouge under hydrothermal conditions. Experiments were performed on 1-mm-thick gabbro gouge sandwiched between country rock pieces (with gouge inclined 35° to the sample axis) in a triaxial testing system with argon gas as the confining medium. In the first series, experiments were conducted under effective normal stresses of 200 MPa and 300 MPa respectively, with pore pressure of 10 MPa. For temperature over 400 °C, pore pressure of 30 MPa was also applied to implement supercritical water conditions. At temperatures up to 615 °C, slip rate steps ranging from 0.0488 μm/s to 1.22 μm/s were applied to obtain the rate dependence of friction.At 200 MPa effective normal stress and a pore pressure of 10 MPa, the steady state rate dependence a–b shows velocity-weakening behaviour for temperatures between  200 and  310 °C. The higher temperature limit for velocity-weakening behaviour to occur extends up to  510 °C under supercritical water conditions with a pore pressure of 30 MPa. For the limited sliding distance in our experiments, only velocity-strengthening behaviour occurred at 300 MPa effective normal stress. Considering the limited displacement (< 3.5 mm), velocity-weakening behaviour may not be excluded in the high effective normal stress case for temperature below  510 °C.The coefficient of friction shows an increasing trend with increasing temperature in the low temperature range. The cut-off temperatures for the increasing trend are  250 °C and  440 °C, respectively for the 200 MPa and 300 MPa effective normal stress cases. Above the cut-off temperatures, the coefficient of friction at 1.83 mm permanent displacement varies around an average of 0.73, which is identical to the average for the oven-dried case [He, C., Yao, W., Wang, Z., Zhou, Y., 2006. Strength and stability of frictional sliding of gabbro gouge at elevated temperatures. Tectonophysics 427, 217–229, doi:10.1016/j.tecto.2006.05.023]. Together with the small value of rate dependence (a–b < 0.0073) for the whole temperature range, these results indicate the absence of fluid-assisted creep.With the result of our experiments as a constraint on strength of frictional sliding, comparison between converted strength for strike–slip faults and creep strength of gabbro-like rocks implies fracturing and faulting behaviours in the lower crust of a cool area (Zhangbei) in North China.     

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.     

Glide systems of hematite single crystals in deformation experiments

The critical resolved shear stresses (CRSSs) of hematite crystals were determined in compression tests for r-twinning, c-twinning and {a}<m>-slip in the temperature range 25 °C to 400 °C, at 400 MPa confining pressure, and a strain rate of 10^− 5 s^− 1 by Hennig-Michaeli, Ch., Siemes, H., 1982. Experimental deformation of hematile crstals betwen 25 °C and 400 °C at 400 MPa confining pressure. In: Schreyer, W. (Ed.) High Pressure Research in Geoscience, Schweizerbart'sche Verlagsbuchhandlung, Stuttgart, p. 133–150. In the present contribution newly performed experiments on hematite single crystals at temperatures up to 800 °C at strain rates of 10^− 5 s^− 1 and 300 MPa confining pressure extends the knowledge about the CRSS of twin and slip modes. Optical observations, neutron diffraction goniometry, SEM forescatter electron images and electron backscatter diffraction are applied in order to identify the glide modes. Both twinning systems and {a}<m>-slip were confirmed by these methods. Besides the known glide systems the existence of the (c)<a>-slip system could be stated. Mechanical data establish that the CRSS of r-twinning decreases from 140 MPa at 25 °C to  5 MPa at 800 °C and for {a}<m>-slip from > 560 MPa at 25 °C to  40 MPa at 700 °C. At room temperature the CRSS for c-twinning is around 90 MPa and at 600 °C  60 MPa. The data indicate that the CRSSs above 200 °C seem to be between the values for r-twinning and {a}<m>-slip. For (c)<a>-slip only the CRSS at 600 °C could be evaluated to  60 MPa. Exact values are difficult to determine because other glide systems are always simultaneously activated.     

Palaeomagnetism of the Loch Doon Granite Complex, Southern Uplands of Scotland: The Late Caledonian palaeomagnetic record and an Early Devonian episode of True Polar Wander

The Southern Uplands terrane is an Ordovician–Silurian back-arc/foreland basin emplaced at the northern margin of the Iapetus Ocean and intruded by granite complexes including Loch Doon (408.3 ± 1.5 Ma) during Early Devonian times. Protracted cooling of this 130 km³ intrusion recorded magnetic remanence comprising a predominant (‘A’) magnetisation linked to initial cooling with dual polarity and mean direction D / I = 237 / 64° (α₉₅ = 4°, palaeopole at 316°E, 21°N). Subsidiary magnetisations include Mesozoic remanence correlating with extensional tectonism in the adjoining Irish Sea Basin (‘B’, D / I = 234/− 59°) and minority populations (‘C’, D / I = 106/− 2° and ‘D’, D / I = 199/1°) recording emplacement of younger ( 395 Ma) granites in adjoining terranes and the Variscan orogenic event. The ‘A’ directions have an arcuate distribution identifying anticlockwise rotation during cooling. A comparable rotation is identified in the Orthotectonic Caledonides to the north and the Paratectonic Caledonides to the south following closure of Iapetus. Continental motion from midsoutherly latitudes ( 40°S) at 408 Ma to equatorial palaeolatitudes by  395 Ma is identified and implies minimum rates of continental movement between 430 and 390 Ma of 30–70 cm/year, more than double maximum rates induced by plate forces and interpreted as a signature of true polar wander. Silurian–Devonian palaeomagnetic data from the British–Scandinavian Caledonides define a 430–385 Ma closed loop comparable to the distributed contemporaneous palaeomagnetic poles from Gondwana. They reconcile pre-430 Ma and post-380 Ma APW from this supercontinent and show that Laurentia–Baltica–Avalonia lay to the west of South America with a relict Rheic Ocean opening to the north which closed to produce Variscan orogeny by a combination of pivotal closure and right lateral transpression.     

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²).     

Ar/Ar thermochronology of the Kampa Dome, southern Tibet: Implications for tectonic evolution of the North Himalayan gneiss domes

Structural and thermochronological studies of the Kampa Dome provide constraints on timing and mechanisms of gneiss dome formation in southern Tibet. The core of Kampa Dome contains the Kampa Granite, a Cambrian orthogneiss that was deformed under high temperature (sub-solidus) conditions during Himalayan orogenesis. The Kampa Granite is intruded by syn-tectonic leucogranite dikes and sills of probable Oligocene to Miocene age. Overlying Paleozoic to Mesozoic metasedimentary rocks decrease in peak metamorphic grade from kyanite + staurolite grade at the base of the sequence to unmetamorphosed at the top. The Kampa Shear Zone traverses the Kampa Granite — metasediment contact and contains evidence for high-temperature to low-temperature ductile deformation and brittle faulting. The shear zone is interpreted to represent an exhumed portion of the South Tibetan Detachment System. Biotite and muscovite ⁴⁰Ar/³⁹Ar thermochronology from the metasedimentary sequence yields disturbed spectra with 14.22 ± 0.18 to 15.54 ± 0.39 Ma cooling ages and concordant spectra with 14.64 ± 0.15 to 14.68 ± 0.07 Ma cooling ages. Petrographic investigations suggest disturbed samples are associated with excess argon, intracrystalline deformation, mineral and fluid inclusions and/or chloritization that led to variations in argon systematics. We conclude that the entire metasedimentary sequence cooled rapidly through mica closure temperatures at  14.6 Ma. The Kampa Granite yields the youngest biotite ⁴⁰Ar/³⁹Ar ages of  13.7 Ma immediately below the granite–metasediment contact. We suggest that this age variation reflects either varying mica closure temperatures, re-heating of the Kampa Granite biotites above closure temperatures between 14.6 Ma and 13.7 Ma, or juxtaposition of rocks with different thermal histories. Our data do not corroborate the “inverse” mica cooling gradient observed in adjacent North Himalayan gneiss domes. Instead, we infer that mica cooling occurred in response to exhumation and conduction related to top-to-north normal faulting in the overlying sequence, top-to-south thrusting at depth, and coeval surface denudation.     

Silicic recharge of multiple rhyolite magmas by basaltic intrusion during the 22.6 ka Okareka Eruption Episode, New Zealand

Deposits of the 22.6 ka Okareka Eruption Episode from Tarawera Volcanic Complex record the sequential and simultaneous eruption of three discrete rhyolite magmas following a silicic recharge event related to basaltic intrusion. The episode started with basaltic eruption ( 0.01 km³ magma), and rapidly changed to a plinian eruption involving a moderate temperature (750 °C), cummingtonite-bearing rhyolite magma (T1) with a volume of  0.3 km³. Hybrid basalt/rhyolite clasts demonstrate direct basaltic intrusion that helped trigger the eruption. Crystals, shards and lapilli of two other rhyolite magmas then joined the eruption sequence. They comprise a cooler (720 °C) crystal-rich biotite–hornblende rhyolite magma (T2) ( 0.3 km³), and a hotter (780 °C), crystal-poor, pyroxene–hornblende rhyolite magma (T3) ( 4.5 km³). All mid to late-stage ash units contain various mixtures of T1, T2 and T3 components with a general increase in abundance of T3 and rapid decline of T1 with time. About 4 km³ of T3 magma was extruded as lavas at the end of the episode. Contrasts in melt composition, crystal and volatile contents, and temperatures influenced viscosity and miscibility, and thus limited pre-eruption mixing of the rhyolite magmas. The eruption sequence and the restricted direct basaltic intrusion into only one magma (T1) is consistent with the rhyolites occupying separate melt pods within a large crystal-mush zone. Melt–crystal equilibria and volatile contents in melt inclusions indicate temporary magma storage depths of < 8 km. Each of the magmas display quartz crystals containing melt inclusions that are compositionally highly evolved relative to the accompanying matrix glass, and thus point to a stage of more complete crystallisation. The matrix glass, enriched in Sr and Ti, represents a re-melting event of underlying the crystal pile induced by basaltic intrusion, presumably part of the same event that erupted scoria at the start of the eruption. This recharge rhyolite melt percolated upward and hybridised with the resident melts in each of the three magma pods. The Okareka episode rhyolites contrast with other well-documented rhyolites that are either continuously or discontinuously zoned, or have been homogenised during re-activation to a uniform composition. Rapid basalt dike intrusion to shallow levels appears to have (prematurely?) triggered the Okareka rhyolites into eruption, so that their early ponding in separate melt pods has been recorded before it could be masked by mixing or stratification had amalgamation into a larger body occurred.     

Timing of thermal stabilization of the Zimbabwe Craton deduced from high-precision Rb–Sr chronology, Great Dyke

Dating low temperature events such as magmatic cooling or (hydro-)thermal surges in Archean and Proterozoic terranes is crucial in defining cratonal thermal stabilization after episodic continental growth during the Archean and Early Proterozoic. Rubidium–Sr chronology is potentially a powerful tool in this regard because of its low closure temperature, i.e., <400 °C in most minerals, but has until now been hampered by its relatively low precision compared to high-temperature chronometers. Consequently, Rb–Sr age investigations have so far failed to provide high-precision age constraints on the cooling of rocks older than 2 Ga. Here, it is demonstrated that internal Rb–Sr microchrons can yield important, high-precision age constraints on the cooling history of Archean intrusions. After careful mineral selection and chemical treatment, a Rb–Sr age of 2543.0 ± 4.4 Ma was obtained from the Archean Great Dyke, Zimbabwe Craton, in contrast to the intrusion age of 2575.8 ± 1 Ma, yielding an ambient average cooling of 5 ± 2 °C/Ma. The non-disturbed magmatic Rb–Sr cooling age of the Great Dyke marks the final stage of Zimbabwe craton stabilization and that the greater craton area did not experience any intensive later reheating event during metamorphic or tectonic events.     

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.     

Development of crystallographic preferred orientations by nucleation and growth of new grains in experimentally deformed quartz single crystals

Deformation experiments have been carried out to investigate the effect of dynamic recrystallisation on crystallographic preferred orientation (CPO) development. Cylindrical samples of natural single crystals of quartz were axially deformed together with 1 vol.% of added water and 20 mg of Mn₂O₃ powder in a Griggs solid medium deformation apparatus in different crystallographic orientations with compression direction: (i) parallel to <c>, (ii) at 45° to <c> and 45° to <a> and (iii) parallel to <a>. The experiments were performed at a temperature of 800 °C, a confining pressure of 1.2 GPa, a strain rate of  10^− 6 s^− 1, to bulk finite strains of  14–36%. The deformed samples were analysed in detail using optical microscopy, electron backscatter diffraction (EBSD) and transmission electron microscopy (TEM). Two different microstructural domains were distinguished in the deformed samples: (i) domains with undulatory extinction and deformation lamellae, and (ii) domains with new recrystallised grains. Within the domains of undulatory extinction, crystal-plastic deformation caused gradual rotations of the crystal lattice up to  30° away from the host orientation. New recrystallised grains show a strong CPO with c-axis maxima at  45° to the compression direction. This is the case in all experiments, irrespective of the initial crystallographic orientation. The results show that c-axes are not continuously rotated towards the new maxima. The new grains thus developed through a mechanism different from subgrain rotation recrystallisation. New grains have a subeuhedral shape and numerous microcavities, voids, fluid channels and fluid inclusions at their grain boundaries. No host control is recorded in misorientation axes across their large angle grain boundaries. New grains might have been created by nucleation from solution in the μm-scale voids and microfractures. The CPO most likely developed due to preferred growth of the freshly precipitated grains with orientations suitable for intracrystalline deformation at the imposed experimental conditions.     

Probing fault zone heterogeneity on the Nojima fault: Constraints from zircon fission-track analysis of borehole samples

Fission-track (FT) thermochronologic analysis was performed on zircon separates from rocks in and around the Nojima fault, which was activated during the 1995 Kobe earthquake. Samples were collected from the University Group 500 m (UG-500) borehole and nearby outcrops. FT lengths in zircons from localities > 25 m away from the fault plane as well as one 0.1 m away from the fault in the footwall are characterized by concordant mean values of  10–11 μm and unimodal distributions with negative skewness, which showed no signs of appreciable reduction in FT length. In contrast, those adjacent (< 3 m) to the fault at depths on the hanging wall side showed significantly reduced mean track lengths of  6–8 μm and distributions having a peak around 6–7 μm with rather positive skewness. The former pattern is interpreted to reflect cooling through the zircon partial annealing zone (ZPAZ), without later, partial thermal overprints. The latter indicates substantial track shortening due probably to secondary heating by a thermal event(s) that locally perturbed the geothermal structure. Modeled zircon FT length and age data of partially annealed samples from the UG-500 borehole revealed a cooling episode in the ZPAZ that started at  4 Ma within  3 m from the fault plane, whereas those from the Geological Survey of Japan 750 m borehole record cooling started at  31–38 Ma within  25 m from the fault. On the basis of one-dimensional heat conduction modeling as well as the consistency between the degree of FT annealing and the degree of deformation/alteration of borehole rocks, these cooling ages in both boreholes are interpreted as consequences of ancient thermal overprints by heat transfer or dispersion via fluids in the fault zone. Together with the zircon FT data of a pseudotachylyte layer recently analyzed, it is suggested that the present Nojima fault system was reactivated in the Middle Quaternary from an ancient fault initiated at  56 Ma at mid-crustal depths. Also shown is a temporal/spatial variation in terms of the thermal anomalies recorded in the fault rocks, implying heterogeneity of hot fluid flows in the fault zone.     

Diagenetic mineralization in Pennsylvanian coals from Indiana, USA: C/C and O/O implications for cleat origin and coalbed methane generation

Cleats and fractures in southwestern Indiana coal seams are often filled with authigenic kaolinite and/or calcite. Carbon- and oxygen-stable isotope ratios of kaolinite, calcite, and coalbed CO₂ were evaluated in combination with measured values and published estimates of δ¹⁸O of coalbed paleowaters that had been present at the time of mineralization. δ¹⁸O_mineral and δ¹⁸O_water values jointly constrain the paleotemperature of mineralization. The isotopic evidence and the thermal and tectonic history of this part of the Illinois Basin led to the conclusion that maximum burial and heat-sterilization of coal seams approximately 272 Ma ago was followed by advective heat redistribution and concurrent precipitation of kaolinite in cleats at a burial depth of < 1600 m at  78 ± 5 °C. Post-Paleozoic uplift, the development of a second generation of cleats, and subsequent precipitation of calcite occurred at shallower burial depth between  500 to  1300 m at a lower temperature of 43 ± 6 °C. The available paleowater in coalbeds was likely ocean water and/or tropical meteoric water with a δ¹⁸O_water  − 1.25‰ versus VSMOW. Inoculation of coalbeds with methanogenic CO₂-reducing microbes occurred at an even later time, because modern microbially influenced ¹³C-enriched coalbed CO₂ (i.e., the isotopically fractionated residue of microbial CO₂ reduction) is out of isotopic equilibrium with ¹³C-depleted calcite in cleats.     

Finite-fault slip models for the 15 April 1979 (Mw 7.1) Montenegro earthquake and its strongest aftershock of 24 May 1979 (Mw 6.2)

We revisit the April 1979 Montenegro earthquake sequence to invert for finite-fault slip models for the mainshock of 15 April 1979 (M_w 7.1) and of the strongest aftershock of 24 May 1979 (M_w 6.2) using P, SH and SV waveforms, retrieved from IRIS data center. We also used body waveform modelling inversion to confirm the focal mechanism of the mainshock as a pure thrust mechanism and rule out the existence of considerable strike slip component in the motion. The mainshock occurred along a shallow (depth 7 km), low angle (14°) thrust fault, parallel to the coastline and dipping to the NE. Our preferred slip distribution model for the mainshock indicates that rupture initiated from SE and propagated towards NW, with a speed of 2.0 km/s. Moment was released in a main slip patch, confined in an area of L  50 km × W  23 km. The maximum slip ( 2.7 m) occurred  30 km to the NW of the hypocenter (location of rupture initiation). The average slip is 49 cm and the total moment release over the fault is 4.38e19 Nm. The slip model adequately fits the distribution of the M_w ≥ 4.3 aftershocks, as most of them are located in the regions of the fault plane that did not slip during the mainshock. The 24 May 1979 (M_w 6.2) strongest aftershock occurred  40 km NW of the mainshock. Our preferred slip model for this event showed a characteristic two-lobe pattern, where each lobe is  7.5 × 7.5 km². Rupture initiated in the NW lobe, where the slip obtained its maximum value of 45 cm, very close to the hypocenter, and propagated towards the south-eastern lobe where it reached another maximum value — for this lobe — of 30 cm, approximately 10 km away from the hypocenter. To indirectly validate our slip models we produced synthetic PGV maps (Shake maps) and we compared our predictions with observations of ground shaking from strong motion records. All comparisons were made for rock soil conditions and in general our slip models adequately fit the observations especially at the closest stations where the shaking was considerably stronger. Through the search of the parameter space for our inversions we obtained an optimum location for the mainshock at 42.04°N and 19.21° E and we also observed that better fit to the observations was obtained when the fault was modeled as a blind thrust fault.     

A late Quaternary paleotemperature record from Hanging Lake, northern Yukon Territory, eastern Beringia

The late Quaternary paleoclimate of eastern Beringia has primarily been studied by drawing qualitative inferences from vegetation shifts. To quantitatively reconstruct summer temperatures, we analyzed lake sediments for fossil chironomids, and additionally we analyzed the sediments for fossil pollen and organic carbon content. A comparison with the δ¹⁸O record from Greenland indicates that the general climatic development of the region throughout the last glaciation–Holocene transition differed from that of the North Atlantic region. Between  17 and 15 ka, mean July air temperature was on average 5°C colder than modern, albeit a period of near-modern temperature at  16.5 ka. Total pollen accumulation rates ranged between  180 and 1200 grains cm^− 2 yr^− 1. At  15 ka, approximately coeval with the Bølling interstadial, temperatures again reached modern values. At  14 ka, nearly 1000 yr after warming began, Betula pollen percentages increased substantially and mark the transition to shrub-dominated pollen contributors. Chironomid-based inferences suggest no evidence of the Younger Dryas stade and only subtle evidence of an early Holocene thermal maximum, as temperatures from  15 ka to the late Holocene were relatively stable. The most recognizable climatic oscillation of the Holocene occurred from  4.5 to 2 ka.     

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.     

Extensional Late Paleozoic deformation on the western margin of Pangea, Patlanoaya area, Acatlán Complex, southern Mexico

New mapping in the northern part of the Paleozoic Acatlán Complex (Patlanoaya area) records several ductile shear zones and brittle faults with normal kinematics (previously thought to be thrusts). These movement zones separate a variety of units that pass structurally upwards from: (i) blueschist-eclogitic metamorphic rocks (Piaxtla Suite) and mylonitic megacrystic granites (Columpio del Diablo granite ≡ Ordovician granites elsewhere in the complex); (ii) a gently E-dipping, listric, normal shear zone with top to the east kinematic indicators that formed under upper greenschist to lower amphibolite conditions; (iii) the Middle–Late Ordovician Las Minas quartzite (upper greenschist facies psammites with minor interbedded pelites intruded by mafic dikes and a leucogranite dike from the Columpio del Diablo granite) unconformably overlain by the Otate meta-arenite (lower greenschist facies psammites and pelites): roughly temporal equivalents are the Middle–Late Ordovician Mal Paso and Ojo de Agua units (interbedded metasandstone and slate, and metapelite and mafic minor intrusions, respectively) — some of these units are intruded by the massive, 461 ± 2 Ma, Palo Liso megacrystic granite: decussate, contact metamorphic muscovite yielded a ⁴⁰Ar/³⁹Ar plateau age of 440 ± 4 Ma; (iv) a steeply-moderately, E-dipping normal fault; (v) latest Devonian–Middle Permian sedimentary rocks (Patlanoaya Group: here elevated from formation status). The upward decrease in metamorphic grade is paralleled by a decrease in the number of penetrative fabrics, which varies from (i) three in the Piaxtla Suite, through (ii) two in the Las Minas unit (E-trending sheath folds deformed by NE-trending, subhorizontal folds with top to the southeast asymmetry, both associated with a solution cleavage), (iii) one in the Otate, Mal Paso, and Ojo de Agua units (steeply SE-dipping, NE–SW plunging, open-close folds), to (iv) none in the Patlanoaya Group. ⁴⁰Ar/³⁹Ar analyses of muscovite from the earliest cleavage in the Las Minas unit yielded a plateau age of 347 ± 3 Ma and show low temperature ages of  260 Ma. Post-dating all of these structures and the Patlanoaya Group are NE-plunging, subvertical folds and kink bands. An E–W, vertical normal fault juxtaposes the low-grade rocks against the Anacahuite amphibolite that is cut by megacrystic granite sheets, both of which were deformed by two penetrative fabrics. Amphibole from this unit has yielded a ⁴⁰Ar/³⁹Ar plateau age of 299 ± 6 Ma, which records cooling through  490 °C and is probably related to a Permo-Carboniferous reheating event during exhumation. The extensional deformation is inferred to have started in the latest Devonian ( 360 Ma) during deposition of the basal Patlanoaya Group, lasting through the rapid exhumation of the Piaxtla Suite at  350–340 Ma synchronous with cleavage development in the Las Minas unit, deposition of the Patlanoaya Group with active fault-related exhumation suggested by Mississippian and Early Permian conglomerates ( 340 and 300 Ma, respectively), and continuing at least into the Middle Permian (≡ 260 Ma muscovite ages). The continuity of Mid-Continent Mississippian fauna from the USA to southern Mexico suggests that this extensional deformation occurred on the western margin of Pangea after closure of the Rheic Ocean.     

Oscillations in the southern extent of the Indo-Pacific Warm Pool during the mid-Holocene

The Indo-Pacific Warm Pool (IPWP) is thought to play a key role in the propagation and amplification of climate changes through its influence on the global distribution of heat and water vapour. However, little is known about past changes in the size and position of the IPWP. In this study, we use a total of 48 modern and fossil coral records from the Mentawai Islands (Sumatra, Indonesia) and Muschu/Koil Islands (Papua New Guinea) to reconstruct oscillations in the extent of the IPWP since the mid-Holocene. We show that reliable estimates of mean sea surface temperature (SST) can be obtained from fossil corals by using low-resolution Sr/Ca analysis of a suite of corals to overcome the large uncertainties associated with mean Sr/Ca-SST estimates from individual coral colonies. The coral records indicate that the southeastern and southwestern margins of the IPWP were cooler than at present between 5500 and 4300 years BP (1.2 °C ± 0.3 °C) and were similarly cool before 6800 years BP. This mid-Holocene cooling was punctuated by an abrupt, short-lived shift to mean SSTs that were warmer than at present between 6600 and 6300 years BP (1.3 °C ± 0.3 °C), while similarly warm conditions may have also existed after 4300 years BP. We suggest that mid-Holocene cooling at our study sites was related to contractions of the southeastern and southwestern margins of the IPWP, associated with the more northerly position of the Inter-tropical Convergence Zone (ITCZ) that accompanied mid-Holocene strengthening of the Asian summer monsoon. Conversely, intervals of abrupt warming appear to correspond with widespread episodes of monsoon weakening and accompanying southward migrations of the ITCZ that caused the IPWP to expand beyond our coral sites. Intervals of a strengthened Asian monsoon and cooling in the southwestern IPWP during the mid-Holocene appear to correspond with a more positive Indian Ocean Dipole (IOD)-like mean configuration across the tropical Indian Ocean, suggesting that the Asian monsoon–IOD interaction that exists at interannual time scales also persists over centennial to millennial scales. Associated mean changes in the Pacific ENSO modes may have also occurred during the mid-Holocene. The dynamic and inter-connected behaviour of the IPWP with tropical climate systems during the mid-Holocene highlights the fundamental importance of the warm pool region for understanding climate change throughout the tropics and beyond.     

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.