Outokumpu revisited: New mineral deposit model for the mantle peridotite-associated Cu–Co–Zn–Ni–Ag–Au sulphide deposits

The metaturbidites of the Palaeoproterozoic Jormua–Outokumpu thrust belt in eastern Finland enclose m- to km-scale ultramafic massifs that are distributed over an area of more than 5000 km². These bodies, which almost entirely consist of highly depleted mantle peridotites (now metaserpentinites and metaperidotites), are intimately associated with massive to semimassive, polymetallic Cu–Co–Zn–Ni–Ag–Au sulphide deposits that sustained mining in the region between 1913 and 1988. Currently, one deposit (Kylylahti) is proceeding into a definitive feasibility study emphasising the renewed economic interest for Outokumpu-type deposits.The origin of these Outokumpu-type Cu–Co–Zn–Ni–Ag–Au deposits is now re-interpreted to be polygenetic. First, their formation requires deposition of a Cu-rich proto-ore within peridotitic sea floor at  1950 Ma. Close modern analogues to the proto-ore setting include, for example, the Logatchev and Rainbow fields at the Mid-Atlantic Ridge, where venting of high-T–low-pH hydrothermal fluid resulted in accumulations of Cu–Zn–Co–Ag–Au sulphides on serpentinised ultramafic seafloor. Second, the Ni-rich composition of Outokumpu sulphide ores calls for a separate source for nickel: Some 40 Ma after the deposition of the Cu-rich proto-ore – concomitant with the obduction of the ultramafic massifs – disseminated Ni sulphides formed through chemical interaction between obducting peridotite massifs and adjacent black schists. This process was related to listwaenite–birbirite type carbonate–silica alteration at margins of the ultramafic massifs. Due to this alteration, silicate nickel was released from the primary Fe–Mg silicates and redeposited as Ni sulphides in the alteration fringes of the massifs.We propose that syntectonic mixing of these two “end-member” sulphides, i.e., the primary Cu-rich proto-ore and the secondary Ni-sulphide disseminations, resulted in the uncommon metal combination of the Outokumpu-type sulphides. Late tectonic solid-state re-mobilisation, related to the duplexing of the ore by isoclinal folding, upgraded the sulphides into economic deposits.     

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.     

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.     

Geologic, fluid inclusion and isotopic characteristics of the Jinding Zn–Pb deposit, western Yunnan, South China: A review

With a reserve of  200 Mt ore grading 6.08% Zn and 1.29% Pb (i.e., a metal reserve of  15 Mt) hosted in Cretaceous and Tertiary terrestrial rocks, the Jinding deposit is the largest Zn–Pb deposit in China, and also the youngest sediment-hosted super giant Zn–Pb deposit in the world. The deposit mainly occurs in the Jinding dome structure as tabular orebodies within breccia-bearing sandstones of the Palaeocene Yunlong Formation (autochthonous) and in the overlying sandstones of the Early Cretaceous Jingxing Formation (allochthonous). The deposit is not stratiform and no exhalative sedimentary rocks have been observed. The occurrence of the orebodies, presence of hangingwall alteration, and replacement and open-space filling textures all indicate an epigenetic origin. Formation of the Jinding Zn–Pb deposit is related to a period of major continental crust movement during the collision of the Indian and Eurasian Plates. The westward thrusts and dome structure were successively developed in the Palaeocene sedimentary rocks in the ore district, and Zn–Pb mineralisation appears to have taken place in the early stage of the doming processes.The study of fluid inclusions in sphalerite and associated gangue minerals (quartz, celestine, calcite and gypsum) shows that homogenisation temperatures ranged from 54 to 309 °C and cluster around 110 to 150 °C, with salinities of 1.6 to 18.0 wt.% NaCl equiv. Inert gas isotope studies from inclusions in ore- and gangue-minerals reveal 2.0 to 15.6% mantle He, 53% mantle Ne and a considerable amount of mantle Xe in the ore-forming fluids. The Pb-isotope composition of ores shows that the metal is mainly of mantle origin, mixed with a lesser amount of crustal lead. The widely variable and negative δ³⁴S values of Jinding sulphides suggest that thermo-chemical or bacterial sulphate reduction produced reduced sulphur for deposition of the Zn–Pb sulphides. The mixing of a mantle-sourced fluid enriched in metals and CO₂ with reduced sulphide-bearing saline formation water in a structural–lithologic trap may have been the key mechanism for the formation of the Jinding deposit.The Jinding deposit differs from known major types of sediment-hosted Zn–Pb deposits in the world, including sandstone-type (SST), Mississippi Valley type (MVT) and sedimentary-exhalative (SEDEX). Although the fine-grained ore texture and high Zn/Pb ratios are similar to those in SEDEX deposits, the Jinding deposit lacks any exhalative sedimentary rocks. Like MVT deposits, Jinding is characterised by simple mineralogy, epigenetic features and involvement of basinal brines in mineralisation, but its host rocks are mainly sandstones and breccia-bearing sandstones. The Jinding deposit is also different from SST deposits with its high Zn/Pb ratios, among other characteristics. Most importantly, the Jinding deposit was formed in an intracontinental terrestrial basin with an active tectonic history in relation to plate collision, and mantle-sourced fluids and metals played a major role in ore formation, which is not the case for SEDEX, MVT, and SST. We propose that Jinding represents a new type of sediment-hosted Zn–Pb deposit, named the ‘Jinding type’.     

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.     

Gas–water–rock interactions in sedimentary basins: CO2 sequestration in the Frio Formation, Texas, USA

To investigate the potential for the geologic storage of CO₂ in saline sedimentary aquifers, 1600 ton of CO₂ were injected at 1500 m depth into a 24-m sandstone section of the Frio Formation — a regional reservoir in the US Gulf Coast. Fluid samples obtained from the injection and observation wells before, during and after CO₂ injection show a Na–Ca–Cl type brine with 93,000 mg/L TDS and near saturation of CH₄ at reservoir conditions. As injected CO₂ gas reached the observation well, results showed sharp drops in pH (6.5 to 5.7), pronounced increases in alkalinity (100 to 3000 mg/L as HCO₃) and Fe (30 to 1100 mg/L), and significant shifts in the isotopic compositions of H₂O and DIC. Geochemical modeling indicates that brine pH would have dropped lower, but for buffering by dissolution of calcite and Fe oxyhydroxides. Post-injection results show the brine gradually returning to its pre-injection composition.     

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.     

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.     

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.     

Structural and stratigraphical constraints on the kinematics history of the Southern Tan–Lu Fault Zone during the Mesozoic Anhui Province, China

The Tan–Lu Fault Zone (TLFZ) extends in a NNE–SSW direction for more than 2000 km in Eastern China. It has been considered either as a major sinistral strike-slip fault, as a suture zone or as a normal fault. We have conducted a structural analysis of the southern segment of this fault zone (STLFZ) in the Anhui Province. The ages (Triassic to Palaeocene) of the formations affected by the faults have been re-appraised taking into account recent stratigraphical studies to better constraint the ages of the successive stages of the kinematics of the STLFZ. Subsequently, the kinematics of the faults is presented in terms of strain/stress fields by inversion of the striated fault set data. Finally, the data are discussed in the light of the results obtained by previous workers.We propose the following history of the STLFZ kinematics during the Mesozoic. At the time of collision, a  NNE orientated Tan–Lu margin probably connected two margins located north of the Dabie and Sulu collision belts. During the Middle–Late Triassic, the SCB has been obliquely subducted below the NCB along this margin which has acted as a compressional transfer zone between the Dabie and Sulu continental subduction zones. The STLFZ has been initiated during the Early Jurassic and has acted as a sinistral transform fault during the Jurassic, following which the NCB/SCB collision stopped. A  NW-trending extension related to metamorphic domes was active during the basal Early Cretaceous ( 135–130 Ma); it has been followed by a NW–SE compression and a NE–SW tension during the middle–late Early Cretaceous ( 127 to  105 Ma, possibly  95 Ma); at that time the TLFZ was a sinistral transcurrent fault within the eastern part of the Asian continent. During the Late Cretaceous–Palaeocene, the STLFZ was a normal fault zone under a WNW–ESE tension.     

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.     

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.     

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.     

Fluid inclusion study of the Plomositas–Los Arcos polymetallic epithermal vein tract, Plomosas district, Sinaloa, Mexico

The most important deposit in the Plomosas–Rosario district, Sinaloa, is the vein tract named Plomositas–Plomosas–Los Arcos. These are NNW–SSE striking veins hosted in rocks of the Lower Volcanic Supergroup (LVS), and also in rocks at the bottom of the Upper Volcanic Supergroup (UVS). Both supergroups belong to the Sierra Madre Occidental. These veins evolved from an early intermediate sulfidation stage (1), rich in base metal sulfides, to a low sulfidation stage (2), rich in silver sulfides and sulfosalts. There is also a 45 m-wide stockwork with native silver and gold. Stage 1 is found in the deeper portion of the veins whereas stage 2 is found in the most shallow portion of the deposit. These stages record fluid inclusion salinities ranging from 7 to 12 wt.% NaCl equiv., and from 0.2 to 3.5 wt.% NaCl equiv., respectively. Homogenization temperatures range from 120 °C for surface samples to 200 °C at a depth of 320 m. The low homogenization temperatures recorded, and the dispersion of veins within host rocks as veinlets, suggest that this deposit formed at shallow depths and was probably blind.     

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.     

Dating pseudotachylyte of the Asuke Shear Zone using zircon fission-track and U–Pb methods

Zircon fission-track (FT) and U–Pb analyses were performed on zircon extracted from a pseudotachylyte zone and surrounding rocks of the Asuke Shear Zone (ASZ), Aichi Prefecture, Japan. The U–Pb ages of all four samples are  67–76 Ma, which is interpreted as the formation age of Ryoke granitic rocks along the ASZ. The mean zircon FT age of host rock is 73 ± 7 (2σ) Ma, suggesting a time of initial cooling through the zircon closure temperature. The pseudotachylyte zone however, yielded a zircon FT age of 53 ± 9 (2σ) Ma, statistically different from the age of the host rock. Zircon FTs showed reduced mean lengths and intermediate ages for samples adjacent to the pseudotachylyte zone. Coupled with the new zircon U–Pb ages and previous heat conduction modeling, the present FT data are best interpreted as reflecting paleothermal effects of the frictional heating of the fault. The age for the pseudotachylyte coincides with the change in direction of rotation of the Pacific plate from NW to N which can be considered to initialize the NNE–SSW trending sinistral–extensional ASZ before the Miocene clockwise rotation of SW Japan. The present study demonstrates that a history of fault motions in seismically active regions can be reconstructed by dating pseudotachylytes using zircon FT thermochronology.     

Analysis of crustal deformation in Luzon, Philippines using geodetic observations and earthquake focal mechanisms

We utilize regional GPS velocities from Luzon, Philippines, with focal mechanism data from the Harvard Centroid Moment Tensor (CMT) Catalog, to constrain tectonic deformation in the complex plate boundary zone between the Philippine Sea Plate and Eurasia (the Sundaland block). Processed satellite imagery and digital elevation models are used with existing gravity anomaly, seismicity, and geologic maps to define a suite of six elastic blocks. Geodetic and focal mechanism data are inverted simultaneously to estimate plate rotations and fault-locking parameters for each of the tectonic blocks and faults comprising Luzon. Major tectonic structures that were found to absorb the plate convergence include the Manila Trench (20–100 mm yr^− 1) and East Luzon Trough ( 9–15 mm yr^− 1)/Philippine Trench ( 29–34 mm yr^− 1), which accommodate eastward and westward subduction beneath Luzon, respectively; the left-lateral strike-slip Philippine Fault ( 20–40 mm yr^− 1), and its northward extensions, the Northern Cordillera Fault ( 17–37 mm yr^− 1 transtension), and the Digdig Fault ( 17–27 mm yr^− 1 transpression). The Macolod Corridor, a zone of active volcanism, crustal thinning, extension, and extensive normal and strike-slip faulting in southwestern Luzon, is associated with left-lateral, transtensional slip of  5–10 mm yr^− 1. The Marikina Fault, which separates the Central Luzon block from the Southwestern Luzon block, reveals  10–12 mm yr^− 1 of left-lateral transpression. Our analysis suggests that much of the Philippine Fault and associated splays are locked to partly coupled, while the Manila and Philippine trenches appear to be poorly coupled. Luzon is best characterized as a tectonically active plate boundary zone, comprising six mobile elastic tectonic blocks between two active subduction zones. The Philippine Fault and associated intra-arc faults accommodate much of the trench-parallel component of relative plate motion.     

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.     

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

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.