Petrofabrics and seismic properties of garnet peridotite from the UHP Sulu terrane (China): Implications for olivine deformation mechanism in a cold and dry subducting continental slab

Lattice-preferred orientations (LPO) of olivine, diopside, enstatite and garnet from the Zhimafang garnet peridotite body in the Sulu ultrahigh-pressure (UHP) metamorphic terrane (China) were measured using the electron backscatter diffraction (EBSD) technique. The peridotite was captured from a mantle wedge immediately adjacent the subducted Yangtze slab and then experienced the UHP metamorphism at 750–950 °C and 4–7 GPa. The olivine LPO is characterized by the [001] axis close to the stretching lineation and the (100) plane subparallel to the foliation, indicating the prevailing of (100) [001] slip. Enstatite LPO displays the dominance of (100) [001] slip. Diopside developed complex LPO patterns that are difficult to explain using a single slip system of (100) [001]. Garnet is almost randomly oriented due to its low volume fractions, cubic symmetry and the presence of numerous slip systems. Calculated seismic properties of the peridotite yield a maximum P-wave velocity normal to the foliation and a minimum along the foliation, with anisotropy up to 8% in strongly sheared samples. The S-wave velocity pattern is complex but the fast polarization plane generally normal to the foliation. The inferred shear sense from the olivine LPO is top-to-SE, in contrary to exhumation-induced top-to-NW thrusting recorded in the quartz LPO, implying that the olivine LPO formed at early UHP metamorphic conditions. The olivine crystals have relatively low water contents (141–475 H/10⁶ Si), indicating a fluid-deficient environment for the LPO formation. The present study suggests that a combination of low temperature and UHP plays a much more important role than the water content to promote the activation of (100) [001] slip in olivine.     

Large-strain deformation and strain partitioning in polyphase rocks: Dislocation creep of olivine–magnesiowüstite aggregates

Aggregates composed of olivine and magnesiowüstite have been deformed to large strains at high pressure and temperature to investigate stress and strain partitioning, phase segregation and possible localization of deformation in a polyphase material. Samples with 20 vol.% of natural olivine and 80 vol.% of (Mg_0.7Fe_0.3)O were synthesized and deformed in a gas-medium torsion apparatus at temperatures of 1127 °C and 1250 °C, a confining pressure of 300 MPa and constant angular displacement rates equivalent to constant shear strain rates of 1–3.3 × 10^− 4 s^− 1. The samples deformed homogeneously to total shear strains of up to γ  15. During constant strain rate measurements the flow stress remained approximately stable at 1250 °C while it progressively decreased after the initial yield stress at the lower temperature. Mechanical data, microstructures and textures indicate that both phases were deforming in the dislocation creep regime. The weaker component, magnesiowüstite, controlled the rheological behavior of the bulk material and accommodated most of the strain. Deformation and dynamic recrystallization lead to grain refinement and to textures that were not previously observed in pure magnesiowüstite and may have developed due to the presence of the second phase. At 1127 °C, olivine grains behaved as semi-rigid inclusions rotating in a viscous matrix. At 1250 °C, some olivine grains remained largely undeformed while deformation and recrystallization of other grains oriented for a-slip on (010) resulted in a weak foliation and a texture typical for pure dry olivine aggregates. Both a-slip and c-slip on (010) were activated in olivine even though the nominal stresses were up to 2 orders of magnitude lower than those needed to activate these slip systems in pure olivine at the same conditions.     

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.     

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

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.     

Lattice-preferred orientation of olivine found in diamond-bearing garnet peridotites in Finsch,South Africa and implications for seismic anisotropy

Seismic anisotropy in the upper mantle provides important constraints on mantle dynamics, continental evolution and global tectonics and is believed to be produced by the flow-induced lattice-preferred orientation (LPO) of olivine. Recent experimental studies at high pressure and temperature have suggested that the LPO of olivine is affected by pressure in addition to water and stress. However, there has been no report yet for the pressure-induced LPO of natural olivine because samples from the deep upper mantle are rare and often unsuitable for study due to ambiguous foliation and lineation. Here we show evidence of the pressure-induced LPO of natural olivine in diamond-bearing garnet peridotites from Finsch, South Africa. We found that the [010] axes of olivine are aligned subnormal to foliation and that the [001] axes are aligned subparallel to lineation, which is known as B-type LPO of olivine. The equilibrium pressure of the samples, as estimated using geobarometer, was greater than 4 GPa, indicating that the samples originated from a depth greater than ∼120 km. In addition, FTIR spectroscopy of the olivine showed that the samples are dry, with a water content of less than 90 ± 20 ppm H/Si (5.5 ± 1.2 ppm wt. H₂O). These data suggest that the samples are the first natural examples of olivine displaying B-type LPOs produced due to high pressure under dry condition. Our data indicate that the trench-parallel seismic anisotropy observed in many subduction zones in and below subducting slabs at depths greater than ∼90 km under dry condition may be attributed to the pressure-induced olivine fabrics (B-type LPO) and may be interpreted as the entrainment of the sub-lithospheric mantle in the direction of subduction rather than anomalous trench-parallel flow.     

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.     

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.     

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.     

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.     

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.     

Combined Sr isotope and trace element analysis of melt inclusions at sub-ng levels using micro-milling, TIMS and ICPMS

This paper describes a technique, which allows precise and accurate Sr isotope measurement combined with trace element analysis of individual melt inclusions, of sample sizes  1 ng of Sr. The technique involves sampling by micro-milling, chemical dissolution, micro Sr column chemistry, TIMS, and ICPMS analyses. A 10% aliquot of each sample solution is used for trace element analysis by double focusing magnetic sector field ICPMS, while Sr is chemically separated from the remaining 90% and used for ⁸⁷Sr/⁸⁶Sr determinations by TIMS.During the development of the technique outlined above, we documented in detail the potential sources of blank contributions and their magnitude. The average size and Sr isotope composition of our laboratory total procedural blank during this study was 5.4 pg ± 0.3 pg Sr (n = 21) with an ⁸⁷Sr/⁸⁶Sr of 0.7111 ± 0.0002 (2SE, n = 3). The total procedural Rb blank was 1.9 ± 0.7 pg (n = 21). The total procedural blank was found to have minimal effect (< 150 ppm shift) on the ⁸⁷Sr/⁸⁶Sr of sample material containing down to  250 pg Sr. Applying a blank correction allows ‘in house’ standards of this size to be corrected back to within 175 ppm of their accepted values. By applying blank corrections we can confidently measure the Sr isotope composition on sample sizes down to  25 pg Sr to an accuracy better than 400 ppm.The utility of the technique is illustrated by application to a suite of melt inclusions from NW Iceland and their host olivines. It is shown that the effect of a small amount of entrainment of the host olivine during sampling of 50 μm melt inclusions has a negligible effect on the measured Sr isotope and trace element composition. Furthermore, where melt inclusions are < 50 μm it is possible to obtain Sr isotope and trace element data on multiple melt inclusions hosted in a single olivine. This provides similar information to that of the single melt inclusions.     

Natural deformation related to serpentinisation of an ultramafic inclusion within a continental shear zone: The key role of fluids

The role of fluids in the deformation of continental serpentinites is investigated from structural, microstructural and petrographic analyses applied to a km-scale porphyroclast mantled in a viscous matrix of amphibolites. The clast is sited within a shear zone of the Palaeozoic Maures massif, France. Syntectonic fluid–rock interactions occurred from km to mm scales, at first on the clast borders (along the main rheological boundaries) then within the clast. They are accommodated macroscopically by slickenfibers faults and microscopically by shear microcracks within crack-seal veins, typifying an intermediate, brittle–ductile behaviour of serpentinites. Three main stages of deformation–serpentinisation processes occurred in relation with the left-lateral movement of the hosted shear zone. They developed under metamorphic conditions evolving from amphibolites to green-schists facies conditions ( 400 MPa/550 °C to  200 MPa/< 300 °C), as inferred from the surrounding sheared amphibolites. Deformation and serpentinisation increase through time although fluid pressure decreases. If the shape of the inclusion and its orientation relative to the shear zone mainly controlled the deformation pattern though time (P then R' shears), fluid pressure is required for starting deformation–serpentinisation processes along inherited anisotropy planes. Whatever the origin of fluids, they play a key role all along the deformation processes by influencing stress states within the shear zone at the onset of deformation and by changing at various scales and through time behaviour of the rock, depending of the intensity of serpentinisation and the rate of deformation.     

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.     

Komatiites from west Shandong, North China craton: Implications for plume tectonics

The Mengyin komatiites are located at the base of the Taishan Complex ( 2.7 Ga), in the western Shandong greenstone belt of the eastern block of North China craton. The Mengyin komatiites have high-magnesian, low-titanium, and aluminum-undepleted chemical characteristics. Although the ultramafic rocks underwent amphibolite to greenschist facies metamorphism, they preserve remnant igneous spinifex textures. Some komatiite samples even preserve rare fresh olivine. According to the mass equilibrium between the olivine and melt, the potential eruption temperature of the Mengyin komatiites at 1 atm is about 1270 °C which might indicate a plume tectonic setting. It is inferred that komatiites in Mengyin, west Shandong may be part of a plume-related global magmatic event related to the break-up of a supercontinent that included the eastern block of the North China craton.     

Linking source and sink: Evaluating the balance between onshore erosion and offshore sediment accumulation since Gondwana break-up, South Africa

The geomorphic origin and evolution of the tectonically unique interior highland of southern Africa, the Kalahari Plateau, and its flanking low-lying coastal planes, remain largely unresolved because of a lack of regional quantitative analyses of its uplift and erosion history. Here we focus on the southern Cape, South Africa and link onshore denudation, based on new apatite fission track thermochronology results, to offshore sediment accumulation, using abundant well data and a seismic reflection profile. We attempt to relate source and sink in order to resolve some first order issues concerning timing of the exhumation and development of the topographic features of southern Africa. The volume of sediment accumulated off South Africa's south coast is calculated using 173 wells and a seismic reflection profile. A total, uncompacted, sediment volume of 268,500 km³ accumulated off South Africa's south coasts since  136 Ma, in the Outeniqua and Southern Outeniqua Basins. Accumulation volumes and rates were highest in the early Cretaceous (48,800 × 10⁴ km³ at  8150 km³/Ma from  136 to 130 Ma, and 57,500 × 10⁴ km³ at 5750 km³/Ma from  130 to 120 Ma) and mid–late Cretaceous (83,700 × 10⁴ km³ at 3200 km³/Ma from  93 to 67 Ma). Volumes and accumulation rates were lowest for the early–mid-Cretaceous (47,400 × 10⁴ km³ at 1750 km³/Ma from  120 to 93 Ma) and the Cenozoic (31,200 × 10⁴ km³ at 450 km³/Ma from  67 to 0 Ma). Although our analysis shows that the accumulated volume of offshore sediments does not match the calculated volume of onshore erosion, as quantified through apatite fission track thermochronology (e.g. Tinker, J.H., de Wit, M.J., Brown, R., 2008. Mesozoic exhumation of the 439 southern Cape, South Africa, quantified using apatite fission track thermochronology. Tectonophysics, doi: 10.1016/j.tecto.2007.10.009), the timing of increased sediment accumulation closely matches the timing of increased onshore denudation. This suggests that the greatest volumes of material were transported from source to sink during two distinct Cretaceous episodes, and that the processes driving onshore denudation decreased by an order of magnitude during the Cenozoic.     

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.     

Deformation-induced polymorphic transformation: experimental deformation of kyanite, andalusite, and sillimanite

Torsion experiments were performed on the Al₂SiO₅ polymorphs in the sillimanite stability field to determine basic rheological characteristics and the effect of deformation on polymorphic transformation. The experiments resulted in extensive transformation of andalusite and kyanite to sillimanite. No transformation occurred during the hot-press (no deformation) stage of sample preparation, which was carried out at similar P–T conditions and duration as the torsion experiments. Experiments were conducted on fine-grained (< 15 µm) aggregates of natural andalusite, kyanite and sillimanite at 1250 °C, 300 MPa, and a constant shear strain rate of 2 × 10^− 4/s to a maximum shear strain of 400%. Electron back-scattered diffraction (EBSD) analysis of the experiments revealed development of lattice-preferred orientations, with alignment of sillimanite and andalusite [001] slightly oblique to the shear plane. The kyanite experiment could not be analyzed using EBSD because of near complete transformation to sillimanite. Very little strain ( 30%) is required to produce widespread transformation in kyanite and andalusite. Polymorphic transformation in andalusite and kyanite experiments occurred primarily along 500 µm wide shear bands oriented slightly oblique and antithetic to the shear plane and dominated by sub-µm (100–150 nm) fibrolitic sillimanite. Shear bands are observed across the entire strain field preserved in the torsion samples. Scanning transmission electron microscope imaging shows evidence for transformation away from shear bands; e.g. fibrolitic rims on relict andalusite or kyanite. Relict grains typically have an asymmetry that is consistent with shear direction. These experimental results show that sillimanite is by far the weakest of the polymorphs, but no distinction can yet be made on the relative strengths of kyanite and andalusite. These observations also suggest that attaining high bulk strain energy in strong materials such as the Al₂SiO₅ polymorphs is not necessary for triggering transformation. Strain energy is concentrated along grain boundaries, and transformation occurs by a dynamic recrystallization type process. These experiments also illustrate the importance of grain-size sensitive creep at high strains in a system with simultaneous reaction and deformation.     

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.