Stress fields and fault reactivation angles of the 2000 western Tottori aftershocks and the 2001 northern Hyogo swarm in southwest Japan

We estimated the stress fields of the aftershocks of the 2000 western Tottori earthquake (M_w 6.6) and the northern Hyogo swarm (max M_w 5.2) by a stress tensor inversion of moment tensor solutions reported from the National Research Institute for Earth Science and Disaster Prevention (Japan). The maximum principal stress direction of the western Tottori sequence was estimated as N107°E with a strike–slip regime. In the northern Hyogo swarm, the orientations of the principal stress directions could not be well constrained by the observed data, but after examining the detailed characteristics of the solution, we obtained a most probable solution of N113°E for the σ₁ direction. These solutions are consistent with the maximum horizontal directions roughly estimated from the strike directions of large earthquakes occurring geographically between these two seismic activities. We measured the angle between each fault–slip direction and maximum principal stress direction to investigate the frictional properties of earthquakes. The distribution of the angles was forward modeled to estimate the coefficient of friction and the stress ratio, assuming uniformly distributed fault orientations. For the western Tottori sequence, a homogeneous stress field with a coefficient of friction less than 0.4 was estimated. A high stress level was also suggested because very little change occurred in the stress field during the mainshock. For the northern Hyogo sequence, the coefficient of friction was estimated to be between 0.5 and 1.0.     

Tectonic analysis of an oceanic transform fault zone based on fault-slip data and earthquake focal mechanisms: the Húsavík–Flatey Fault zone, Iceland

The Húsavík–Flatey Fault (HFF) is an oblique dextral transform fault, part of the Tjörnes Fracture Zone (TFZ), that connects the North Volcanic Zone of Iceland and the Kolbeinsey Ridge. We carry out stress inversion to reconstruct the paleostress fields and present-day stress fields along the Húsavík–Flatey Fault, analysing 2700 brittle tectonic data measured on the field and about 700 earthquake focal mechanisms calculated by the Icelandic Meteorological Office. This allows us to discuss the Latest Cenozoic finite deformations (from the tectonic data) as well as the present-day deformations (from the earthquake mechanisms). In both these cases, different tectonic groups are reconstructed and each of them includes several distinct stress states characterised by normal or strike-slip faulting. The stress states of a same tectonic group are related through stress permutations (σ₁−σ₂ and σ₂−σ₃ permutations as well as σ₁−σ₃ reversals). They do not reflect separate tectonic episodes. The tectonic groups derived from the geological data and the earthquake data have striking similarity and are considered to be related. The obliquity of the Húsavík–Flatey Fault implies geometric accommodation in the transform zone, resulting mainly from a dextral transtension along an ENE–WSW trend. This overall mechanism is subject to slip partitioning into two stress states: a Húsavík–Flatey Fault-perpendicular, NE–SW trending extension and a Húsavík–Flatey Fault-parallel, NW–SE trending extension. These three regimes occur in various local tectonic successions and not as a regional definite succession of tectonic events. The largest magnitude earthquakes reveal a regional stress field tightly related to the transform motion, whereas the lowest magnitude earthquakes depend on the local stress fields. The field data also reveal an early extension trending similar to the spreading vector. The focal mechanism data do not reflect this extension, which occurred earlier in the evolution of the HFF and is interpreted as a stage of structural development dominated by the rifting process.     

Cenozoic post-collisional brittle tectonic history and stress reorientation in the High Zagros Belt (Iran, Fars Province)

The Zagros fold-and-thrust belt of SW-Iran is among the youngest continental collision zones on Earth. Collision is thought to have occurred in the late Oligocene–early Miocene, followed by continental shortening. The High Zagros Belt (HZB) presents a Neogene imbricate structure that has affected the thick sedimentary cover of the former Arabian continental passive margin. The HZB of interior Fars marks the innermost part of SE-Zagros, trending NW–SE, that is characterised by higher elevation, lack of seismicity, and no evident active crustal shortening with respect to the outer (SW) parts. This study examines the brittle structures that developed during the mountain building process to decipher the history of polyphase deformation and variations in compressive tectonic fields since the onset of collision. Analytic inversion techniques enabled us to determine and separate different brittle tectonic regimes in terms of stress tensors. Various strike–slip, compressional, and tensional stress regimes are thus identified with different stress fields. Brittle tectonic analyses were carried out to reconstruct possible geometrical relationships between different structures and to establish relative chronologies of corresponding stress fields, considering the folding process. Results indicate that in the studied area, the main fold and thrust structure developed in a general compressional stress regime with an average N032° direction of σ₁ stress axis during the Miocene. Strike–slip structures were generated under three successive strike–slip stress regimes with different σ₁ directions in the early Miocene (N053°), late Miocene–early Pliocene (N026°), and post-Pliocene (N002°), evolving from pre-fold to post-fold faulting. Tensional structures also developed as a function of the evolving stress regimes. Our reconstruction of stress fields suggests an anticlockwise reorientation of the horizontal σ₁ axis since the onset of collision and a significant change in vertical stress from σ₃ to σ₂ since the late stage of folding and thrusting. A late right-lateral reactivation was also observed on some pre-existing belt-parallel brittle structures, especially along the reverse fault systems, consistent with the recent N–S plate convergence. However, this feature was not reflected by large structures in the HZB of interior Fars. The results should not be extrapolated to the entire Zagros belt, where the deformation front has propagated from inner to outer zones during the younger events.     

Stress field in the Northern Rhine area, Central Europe, from earthquake fault plane solutions

Fault plane solutions (FPS) from 110 earthquakes in the northern Rhine area with local magnitudes, ranging from 1.0 to 6.1, and occurring between 1976 and 2002 are determined. FPS are retrieved from P-wave first motions using a grid search approach allowing a detailed exploration of the parameter space. The influence of the 1D velocity model on take-off angles and resulting FPS is examined. All events were relocated with a recently developed minimum 1D model of the velocity structure [J. Geophys. Res. (2003)]. Rose diagrams of the orientation of P, T and B axes show a clear preference of trends of P and T axes at N292°E and N27°E, respectively. The majority of B axes trend in northerly directions. Plunges of P and T axes are mostly around 45° while most B axes are subhorizontal. The main direction of the maximum horizontal stress directly inferred from the fault plane solutions is N118°E.To calculate the orientations of the principal stress axes and the shape of the stress tensor, the inversion method of Gephard and Forsyth [J. Geophys. Res. 89 (1984) 9305] was applied to the whole data set and to several subsets of data. The subsets were formed by grouping events from various geological and tectonic areas and by grouping events into different depth ranges. The subset areas include the Lower Rhine Embayment, the Rhenish Massif, the middle Rhine area, the Neuwied Basin and the area known as the Stavelot–Venn Massif. Inversion of the entire data set shows some ambiguity between a strike-slip and extensional stress regime, with a vertical axis for the medium principal stress and a trend of N305°E and N35°E for the σ₁ and σ₃ axis, respectively, as the best fitting tensor. Earthquakes from the Lower Rhine Embayment and, to some degree, from the middle Rhine area indicate an extensional stress regime. In the Lower Rhine Embayment, plunge and trend of the σ₁ axis are 76° and N162°E and for the σ₃ axis 7° and N42°E. The best fitting solution for the area of the Stavelot–Venn Massif is a strike-slip regime with subhorizontal σ₁ and σ₃ axes with a trend of N316°E and N225°E, respectively. Stress orientations found here agree overall with the results from earlier studies based on smaller data sets. The directions of the maximum and minimum horizontal stresses inverted from focal mechanisms agree well with the stress field predicted by the European Stress Map. This confirms earlier interpretations that the stress field of the Rhine Graben system is controlled by plate driving forces acting on the plate boundaries. However, amplitudes of the stresses change on a local scale and with depth. Estimates of the absolute magnitude of principal stresses favor a normal faulting regime in the shallow crust (above 12-km depth) and a strike-slip regime in the lower crust.     

Variation with depth of the stress tensor anisotropy inferred from microfault analysis

In order to study the Pyrenean tectonic phase, a quantitative method of stress analysis using microfault measurements is used on a calcareous plateau located in southern France. The method developed here allows the determination of several tectonic events and the evaluation of (with σ σ₂ σ₃).The Pyrenean compression is seen to occur in two stages, confirming previous geological studies.A recent canyon allows the study of the variation with depth of the R ratio on a vertical cross-section (300 m). With a simple model of gravity and tectonic stresses, the vertical variation of R can be used to estimate quantitatively the Pyrenean tectonic stresses. For realistic values of the parameters in this model, the horizontal tectonic stresses are obtained in the following range: 50–200 bar for the maximum horizontal principal stress, 10–25 bar for the other horizontal stress.These results seem to be consistent with in situ stress measurements, but they are much lower than those predicted by experimental rock mechanics.     

Creep of olivine during hot-pressing

The densification curves for the hot-pressing of pure olivine powders were obtained as a function of grain size (5 μ–2000 μ), temperature (1000–1600°C), and compacting stress (166–298 bars). This range of variables was found to straddle two fields of hot-pressing behavior, one dominated by power-law creep, one by Coble creep. The time required to density a powder to 99% of the single crystal density could be represented by the shorter of the two times: t₁ = 2.2 · 10³σ^−3.4exp(85,000/RT)t₂ = 1.3 · 10⁴σ^−1.5(G)⁺³exp(85,000/RT) where the compacting stress or pressure, σ, is given in bars and the grain size, G, in centimeters. It was also possible to estimate the parameters appropriate to Coble creep in a solid polycrystalline aggregate from the hot-pressing data; and these were:

Full-size image (1K)

The strain rates computed from this formula are close to those predicted by Stocker and Ashby (1973) and those found by Twiss (1976).     

Neotectonic evolution of the Central Betic Cordilleras (Southern Spain)

Paleostress orientations were calculated from fault-slip data of 36 sites located along a traverse through the Central Betic Cordilleras (southern Spain). Heterogeneous fault sets, which are frequent in the area, have been divided into homogeneous subsets by cross-cutting relationships observed in the field and by a paleostress stratigraphy approach applied on each individual fault population. The state of stress was sorted according to main tectonic events and a new chronology is presented of the Miocene to Recent deformation in the central part of the Betic Cordilleras. The deviatoric stress tensors fall into four distinct groups that are regionally consistent and correlate with three Late Oligocene–Aquitanian to Recent major tectonic events in the Betic Cordilleras. The new chronology of the neotectonic evolution includes, from oldest to youngest, the following main tectonic phases:

(1) Late Oligocene–Aquitanian to Early Tortonian: σ₁ subhorizontal N–S, partly E–W directed, σ₃ subvertical; compressional structures (thrusting of nappes, large-scale folding) and strike-slip faulting in the Alborán Domain and the External Zone of the Betic Cordilleras;

(2) Early Tortonian to Pliocene–Pleistocene: σ₁ subvertical, σ₃ subhorizontal NW–SE, partly N–S directed or E–W-directed (radial extension); large-scale normal faulting in the Central Betic Cordilleras and in the oldest Neogene formations of the Granada Basin related to the gravitational collapse of the Betic Cordilleras and the exhumation of the intensely metamorphosed rock series of the Internal Zones, at the same time formation of the Alborán Basin and intramontane basins such as the Granada Basin;

(3) Pleistocene to Recent: (3a) σ₁ subvertical, σ₃ subhorizontal NE–SW with prominent normal faulting, but coevally; (3b) σ₁ subhorizontal NW directed, σ₃ NE–SW subhorizontal with strike-slip faulting. Extensional structures and strike-slip faulting are related to the ongoing convergence of the Eurasian and African Plates and coeval uplift of the Betic Cordilleras. Reactivation of pre-existing fractures and faults was frequently observed. Phase 3 is interpreted as periodic strike-slip and normal faulting events due to a permutation of the principal stress axes, mainly σ₁ and σ₂.

Comment on “Structural controls on a carbon dioxide-driven mud volcano field in the Northern Apennines (Pieve Santo Stefano, Italy): Relations with pre-existing steep discontinuities and seismicity”

Bonini (2009, Structural controls on a carbon dioxide-driven mud volcano field in the Northern Apennines (Pieve Santo Stefano, Italy): relations with pre-existing steep discontinuities and seismicity. Journal of Structural Geology 31, 44–54) presents a 2D mechanical analysis to infer the failure conditions responsible for the seismicity distribution during an M_w = 4.6 seismic sequence nucleating during 2001 in the Northern Apennines. In my view the mechanical analysis presented in this paper has some weakness or is not well constrained, in particular: 1) the assumption of a dip angle of 50°, is not consistent with the activated structures; 2) the P_f = σ₃ condition, difficult to be attained along a cohesionless fault dipping at 50°; 3) the isotropic stress state, i.e. σ₂ = σ₃, that is not consistent with the active or recent stress field in the area.     

Late Cenozoic stress evolution along the Karasu Valley, SE Turkey

This study defines the Mio-Pliocene to present-day stress regime acting at the northeastern corner of the eastern Mediterranean region along the Karasu Valley (i.e., the Amanos Range), taking in the Antakya, Osmaniye and Kahramanmaras provinces. The inversion slip vectors measured on fault planes and chronologies between striations indicate that the stress regime varied from transpressional initially to transtensional, having consistent NW- and NE-trending σ_Hmax (σ₁) and σ_Hmin (σ₃) axes, respectively; there are significantly different mean stress-ratio (R_m) values however. The older mean stress state is characterized by N151±11°E-trending σ₁ and N59±12°E-trending σ₃ axes, and by a mean arithmetic R_m value of 0.76, indicating that the regional stress regime is transpressional. The younger stress regime is characterized by N154±8°E-trending σ₁ and N243±8°E-trending σ₃ axes, and by a mean arithmetic R_m value of 0.17, indicating a transtensional character for this regional stress regime. The low R values of the stress deviators related to the recent stress state reflect normal-component slips. The earthquake focal mechanism inversions confirm that the younger stress regime continues into the Recent. The inversion identifies a transtensional stress regime representing strike-slip and an extensional stress state with a consistent NE-trending σ_Hmin (σ₃) axis. These stress states are characterized by N66°E and N249°E-trending σ₃ axes, respectively. Both significant regional stress regimes induce left-lateral displacement along the southern part of the East Anatolian Fault (EAF, or Amanos Fault). The temporal change, probably in Quaternary time, within the regional stress regime—from transpression to transtension—resulted from the coeval influences of subduction processes in the west–southwest (i.e., along the Cyprus arc), continental collision in the east, and westward escape of the Anatolian block.     

Compressional- and transpressional-stress pattern for Pliocene and Quaternary brittle deformation in fore arc and intra-arc zones (Andes of Central and Southern Chile)

Kinematic analysis of fault slip data for stress determination was carried out on Late Miocene to Quaternary rocks from the fore arc and intra-arc regions of the Chilean Andes, between 33° and 46° south latitudes. Studies of Neogene and Quaternary infilling (the Central Depression), as well as plutonic rocks of the North Patagonian Batholith along the Liquiñe–Ofqui Fault Zone, have revealed various compressional and/or transpressional states of stress. In the Pliocene, the maximum compressional stress (σ₁) was generally oriented east–west. During the Quaternary, the deformation was partitioned into two coeval distinctive states of stress. In the fore arc zone, the state of stress was compressional, with σ₁ oriented in a N–S to NNE–SSW direction. In the intra-arc zone the state of stress was transpressional with σ₁ striking NE–SW. Along the coast, in one site (37°30′S) the Quaternary strain deformation is extensional, with an E–W direction, which can be explained by a co-seismic crustal bending readjustment.     

Crustal architecture above the high-pressure belt of the Grenville Province in the Manicouagan area: new structural, petrologic and U–Pb age constraints

Recent U–Pb age determinations and P–T estimates allow us to characterize the different levels of a formerly thickened crust, and provide further constraints on the make up and tectono-thermal evolution of the Grenville Province in the Manicouagan area. An important tectonic element, the Manicouagan Imbricate zone (MIZ), consists of mainly 1.65, 1.48 and 1.17 Ga igneous rocks metamorphosed under 1400–1800 MPa and 800–900 °C at 1.05–1.03 Ga, during the Ottawan episode of the Grenvillian orogenic cycle, coevally with intrusion of gabbro dykes in shear zones. The MIZ has been interpreted as representing thermally weakened deep levels of thickened crust extruded towards the NW over a parautochthonous crustal-scale ramp. Mantle-derived melts are considered as in part responsible for the high metamorphic temperatures that were registered.New data show that mid-crustal levels structurally above the MIZ are represented by the Gabriel Complex of the Berthé terrane, that consists of migmatite with boudins of 1136±15 Ma gabbro and rafts of anatectic metapelite with an inherited monazite age at 1478±30 Ma. These rocks were metamorphosed at about the same time as the MIZ (metamorphic zircon in gabbro: 1046±2 Ma; single grains of monazite in anatectic metapelite: 1053±2 Ma) and under the same T range (800–900 °C) but at lower P conditions (1000–1100 MPa). They are mainly exposed in an antiformal culmination above a high-strain zone, which has tectonic lenses of high P–T rocks from the MIZ and is intruded by synmetamorphic gabbroic rocks. This zone is interpreted as part of the hangingwall of the MIZ during extrusion. A gap of 400 MPa in metamorphic pressures between the tectonic lenses and the country rocks, together with the broad similarity in metamorphic ages, are consistent with rapid tectonic transport of the high P–T rocks over a ramp prior to the incorporation of the mafic lenses in the hangingwall.Between the antiformal culmination of the Gabriel Complex and the MIZ 1.48 Ga old granulites of the Hart Jaune terrane are exposed. They are intruded by unmetamorphosed 1228±3 Ma gabbro sills and 1166±1 Ma anorthosite. Hart Jaune Terrane represents relatively high crustal levels that truncate the MIZ-Gabriel Complex contact and are preserved in a synformal structure.Farther south, the Gabriel Complex is overlain by the Banded Complex, a composite unit including 1403+32/−25 Ma granodiorite and 1238+16/−13−1202+40/−25 Ma granite. This unit has been metamorphosed under relatively low-P (800 MPa) granulite-facies conditions. Metamorphic U–Pb data, limited to zircon lower intercept ages (971±38 Ma and 996±27 Ma) and a titanite (990±5 Ma) age, are interpreted to postdate the metamorphic peak.The general configuration of units along the section is consistent with extrusion of the MIZ during shortening and, finally, normal displacement along discrete shear zones.     

Intrinsic and scattering attenuation in western India from aftershocks of the 26 January, 2001 Kachchh earthquake

176 vertical-component, short period observations from aftershocks of the M_w 7.7, 26 January, 2001 Kachchh earthquake are used to estimate seismic wave attenuation in western India using uniform and two layer models. The magnitudes (M_w) of the earthquakes are less than 4.5, with depths less than 46 km and hypocentral distances up to 110 km. The studied frequencies are between 1 and 30 Hz. Two seismic wave attenuation factors, intrinsic absorption (Q_i^− 1) and scattering attenuation (Q_s^− 1) are estimated using the Multiple Lapse Time Window method which compares time integrated seismic wave energies with synthetic coda wave envelopes for a multiple isotropic scattering model. We first assume spatial uniformity of Q_i^− 1, Q_s^− 1 and S wave velocity (β). A second approach extends the multiple scattering hypothesis to media consisting of several layers characterized by vertically varying scattering coefficient (g), intrinsic absorption strength (h), density of the media (ρ) and shear wave velocity structure. The predicted coda envelopes are computed using Monte Carlo simulation. Results show that, under the assumption of spatial uniformity, scattering attenuation is greater than intrinsic absorption only for the lowest frequency band (1 to 2 Hz), whereas intrinsic absorption is predominant in the attenuation process at higher frequencies (2 to 30 Hz). The values of Q obtained range from Q_t = 118, Q_i = 246 and Q_s = 227 at 1.5 Hz to Q_t ≈ 4000, Q_i ≈ 4600 and Q_s ≈ 33,300 at 28 Hz center frequencies, being Q_t^− 1 a measure of total attenuation. Results also show that Q_i^− 1, Q_s^− 1 and Q_t^− 1 decrease proportional to f^−ν. Two rates of decay are clearly observed for the low (1 to 6 Hz) and high (6 to 30 Hz) frequency ranges. Values of ν are estimated as 2.07 ± 0.05 and 0.44 ± 0.09 for total attenuation, 1.52 ± 0.21 and 0.48 ± 0.09 for intrinsic absorption and 3.63 ± 0.07 and 0.06 ± 0.08 for scattering attenuation for the low and high frequency ranges, respectively. Despite the lower resolution in deriving the attenuation parameters for a two layered crust, we find that scattering attenuation is comparable to or smaller than the intrinsic absorption in the crust whereas intrinsic absorption dominates in the mantle. Also, for a crustal layer of thickness 42 km, intrinsic absorption and scattering estimates in the crust are lower and greater than those of the mantle, respectively.     

The 2003 earthquake swarm in Anjalankoski, south-eastern Finland

During May 2003 a swarm of 16 earthquakes (M_L = 0.6–2.1) occurred at Anjalankoski, south-eastern Finland. The activity lasted for three weeks, but additional two events were observed at the same location in October 2004. A comparison of the waveforms indicated that the source mechanisms and the hypocentres of the events were nearly identical.A relative earthquake location method was applied to better define the geometry of the cluster and to identify the fault plane associated with the earthquakes. The relocated earthquakes aligned along an ENE–WSW trending zone, with a lateral extent of about 1.0 km by 0.8 km. The relative location and the waveform-modelling of depth sensitive surface wave (Rg) and S-to-P converted body wave (sP) phases indicated that the events were unusually shallow, most likely occurring within the first 2 km of the surface. The revised historical earthquake data confirm that shallow swarm-type seismicity is characteristic to the area.The focal mechanism obtained as a composite solution of the five strongest events corresponds to dip-slip motion along a nearly vertical fault plane (strike = 250°, dip = 80°, rake = 90°). The dip and strike of this nodal plane as well as the relocated hypocentres coincide with an internal intrusion boundary of the Vyborg rapakivi batholith.The events occur under a compressive local stress field, which is explained by large gravitational potential energy differences and ridge-push forces. Pore-pressure changes caused by intrusion of ground water and/or radon gas into the fracture zones are suggested to govern the swarm-type earthquake activity.     

High-temperature flow of wet polycrystalline enstatite

Previous experiments by Raleigh et al. (1971) have shown that at strain rates of 10⁻².sec⁻¹ to 10⁻⁷.sec⁻¹ only slip occurs in dry enstatite at temperatures above 1300°C and 1000°C, respectively.The present experiments have been conducted on polycrystalline enstatite under wet conditions in this regime where enstatite only slips, polygonizes and recrystallizes. Slip occurs throughout the whole regime on the system (100)[001] and at strains greater than 40% the system (010)[001] is observed. Polygonization and intragranular recrystallization begin at about 1300°C and 10⁻⁴.sec⁻¹ and the orientation of these neoblasts is host-controlled. At lower strain rates intergranular neoblasts develop and their fabric is one of [100] maximum parallel with σ₁ and [010] and [001] girdles in the σ₂ = σ₃ plane, similar to those in natural enstatite tectonites.Dislocation substructures of experimentally deformed enstatite have been examined by transmission electron microscopy. The samples were deformed within the field in which slip polygonization and recrystallization are the dominant deformation mechanisms. Samples within this regime have microstructures that are characterized by stacking faults and partial dislocations. Under the conditions of steady-state flow in olivine, these microstructures inhibit the operation of recovery mechanisms in enstatite.Other samples deformed within the polygonization and recrystallization field have microstructures that confirm the optical observations of intragranular and intergranular growth of neoblasts. It is suggested that the former result from strain-induced tilt of subrains, whereas the latter may result from bulge nucleation into adjacent subgrains.Mechanical data from constant strain-rate experiments at steady state, stress relaxation and temperature-differential creep tests are best fit to a power-law creep equation with the stress exponent, n~3 and the apparent activation energy for creep, Q~65 kcal/mole. Extrapolation of this equation to a representative natural geologic strain rate of 10⁻⁴. sec⁻¹, over the temperature interval 1000–2000°C, gives an effective viscosity range of 10²⁰–10¹⁸ poise and stresses in the range of 7-0.1 bar, respectively. Comparison with corrected wet-olivine mechanical data (Carter, 1976) over the same environment indicates that olivine is consistently the weaker of the two minerals and will recrystallize whilst enstatite will only slip and kink, thus accounting for the different habits of olivine and enstatite in ultramafic tectonites.     

Seismic velocities of granulite-facies xenoliths from Central Ireland: Implications for lower crustal composition and anisotropy

V_p and V_s values have been measured experimentally and calculated for granulite-facies lower crustal xenoliths from central Ireland close to the Caledonian Iapetus suture zone. The xenoliths are predominantly foliated and lineated metapelitic (garnet–sillimanite–K-feldspar) granulites. Their metapelitic composition is unusual compared with the mostly mafic composition of lower crustal xenoliths world-wide. Based on thermobarometry, the metapelitic xenoliths were entrained from depths of c. 20–25 ± 3.5 km and rare mafic granulites from depths of 31–33 ± 3.4 km. The xenoliths were emplaced during Lower Carboniferous volcanism and are considered to represent samples of the present day lower crust.V_p values for the metapelitic granulites range between 6.26 and 7.99 km s^{− 1} with a mean value of 7.09 ± 0.4 km s^{− 1}. Psammite and granitic orthogneiss samples have calculated V_p values of 6.51 and 6.23 km s^{− 1}, respectively. V_s values for the metapelites are between 3.86 and 4.34 km s^{− 1}, with a mean value of 4.1 ± 0.15 km s^{− 1}. The psammite and orthogneiss have calculated V_s values of 3.95 and 3.97 km s^{− 1}, respectively.The measured seismic velocities correlate with density and with modal mineralogy, especially the high content of sillimanite and garnet. V_p anisotropy is between 0.15% and 13.97%, and a clear compositional control is evident, mainly in relation to sillimanite abundance. Overall V_s anisotropy ranges from 1% to 11%. Poisson's ratio (σ) lies between 0.25 and 0.35 for the metapelitic granulites, mainly reflecting a high V_p value due to abundant sillimanite in the sample with the highest σ. Anisotropy is probably a function of deformation associated with the closure of the Iapetus ocean in the Silurian as well as later extension in the Devonian. The orientation of the bulk strain ellipsoid in the lower crust is difficult to constrain, but lineation is likely to be NE–SW, given the strike-slip nature of the late Caledonian and subsequent Acadian deformation.When corrected for present-day lower crustal temperature, the experimentally determined V_p values correspond well with velocities from the ICSSP, COOLE I and VARNET seismic refraction lines. Near the xenolith localities, the COOLE I line displays two lower crustal layers with in situ V_p values of 6.85–6.9 and 6.9–8.0 km s^{− 1}, respectively. The upper (lower velocity) layer corresponds well with the metapelitic granulite xenoliths while the lower (higher velocity) layer matches that of the basic granulite xenoliths, though their metamorphic pressures suggest derivation from depths corresponding to the present-day upper mantle.     

A comparative study on the dissolution and solubility of hydroxylapatite and fluorapatite at 25 °C and 45 °C

Dissolution of the synthetic hydroxylapatite (HAP) and fluorapatite (FAP) in pure water was studied at 25 °C and 45 °C in a series of batch experiments. The XRD, FT-IR and SEM analyses indicated that the synthetic, microcrystalline HAP and FAP with apatite structure used in the experiments were found to have no obvious variation after dissolution except that the existence of OH groups in FT-IR spectra for FAP after 2880 h dissolution was observed. During the HAP dissolution (0–4320 h), the aqueous calcium and phosphate concentrations reached the maxima after 120 h and then decreased slowly with time. For the FAP dissolution in pure water, after a transient time of 1440 h (< 60 d), element concentrations and pH became constant suggesting attainment of a steady-state between the solution and solid. During early stages of the FAP dissolution reaction (< 72–120 h), mineral components were released in non-stoichiometric ratios with reacted solution ratios of dissolved Ca:P, Ca:F and P:F being lower than mineral stoichiometric ratios of Ca₅(PO₄)₃F, i.e., 1.67, 5.0 and 3.0, respectively. This indicated that F were preferentially released compared to Ca from the mineral structure. The mean K_sp values were calculated by using PHREEQC for HAP of 10^− 53.28 (10^− 53.02–10^− 53.51) and for FAP of 10^− 55.71 (10^− 55.18–10^− 56.13) at 25 °C, the free energies of formation ΔG_f^o[HAP] and ΔG_f^o[FAP] were calculated to be − 6282.82 kJ/mol and − 6415.87 kJ/mol, respectively.     

Fault slip analysis of Quaternary faults in southeastern Korea

The Quaternary stress field has been reconstructed for southeast Korea using sets of fault data. The subhorizontal direction of the maximum principal stress (σ₁) trended ENE and the direction of the minimum principal stress (σ₃) was nearly vertical. The stress ratio (Φ = (σ₂ − σ₃) / (σ₁ − σ₃)) value was 0.65. Two possible interpretations for the stress field can be made in the framework of eastern Asian tectonics; (1) The σ_Hmax trajectory for southeast Korea fits well with the fan-shaped radial pattern of maximum principal stress induced by the India–Eurasia collision. Thus, we suggest that the main source for this recent stress field in southeast Korea is related to the remote India–Eurasia continental collision. (2) The stress field in Korea shows a pattern similar to that in southwestern Japan. The origin for the E–W trending σ_Hmax in Japan is known to be related to the mantle upwelling in the East China Sea. Thus, it is possible that Quaternary stress field in Korea has evolved synchronously with that in Japan. We suggest further studies (GPS and in situ stress measurement) to test these hypotheses.     

Magnetic fabrics and rock magnetism of Proterozoic dike swarm from the southern São Francisco Craton, Minas Gerais State, Brazil

Magnetic fabric and rock magnetism studies were performed on 32 mafic dikes of a Proterozoic dike swarm from the southern São Francisco Craton (SFC; Minas Gerais State, SE Brazil). Magnetic anisotropies were determined by applying anisotropy of low-field magnetic susceptibility (AMS) and anisotropy of remanent magnetization (ARM). The latter was performed imposing both anhysteretic (total (AAR) and partial pAAR)) and isothermal remanence magnetizations (AIRM). Partial anhysteretic remanence anisotropy was performed based on remanent coercivity spectra from a pilot specimen of each site. In most sites, AMS is dominantly carried by ferromagnetic minerals, however, in some sites, the paramagnetic contribution exceeds 70% of bulk susceptibility. Rock magnetism and thin section analysis allow classifying the dikes as non-hydrothermalized and hydrothermalized. Magnetic measurement shows that the mean magnetic susceptibility is usually lower than 5×10⁻³ (SI). Ti-poor titanomagnetites up to pure magnetite pseudo-single-domain (PSD) grain sizes carry the majority of magnetic fabrics for non-hydrothermalized dikes whereas coarse to fine grained Ti-poor titanomagnetites carry the majority of magnetic fabrics for hydrothermalized dikes.Three primary AMS fabrics are recognized which are coaxial with ARM fabric, except for two dikes, from both non-hydrothermalized and hydrothermalized dikes. Normal AMS fabric surprisingly is not dominant (31%). The parallelism between AMS, pAAR_0–30, pAAR_30–60 and pAAR_60–90 fabrics in the hydrothermalized dikes indicates that magnetic grains formed due to late-stage crystallization or to remobilization of iron oxides due to hydrothermal alteration after dike emplacement have acquired a mimetic fabric coaxial with the primary fabric given by coarse-grained early crystallized Ti-poor titanomagnetites. This fabric is interpreted as magma flow in which the analysis of K_max inclination permitted the inference that the dikes were fed by horizontal or subhorizontal fluxes (K_max<30°). Intermediate AMS fabric is the most important (41%) in the investigated swarm. It is interpreted as due to vertical compaction of a static magma column with the minimum stress along the dike strike. ARM determinations for these sites also remained intermediate except for two dikes. In one of them, AIRM fabric resulted in normal AMS fabric while for the other AAR fabric resulted in inverse AMS fabric. A combination of AMS and ARM fabrics suggest that magmatic fabric for both dikes were overprinted by some late local event, probably related to Brasiliano orogenic processes after dike emplacement. InverseInverse AMS fabric is a minority (four dikes). ARM determinations also remained inverse suggesting a primary origin for inverse AMS fabric.     

Paleo-deviatoric stress magnitudes from calcite twins and related structural permeability evolution in minor faults: Example from the toarcian shale of the French Causses Basin, Aveyron, France

The relationship between deformation and so-called fluid paleotransfers in minor faults has been analysed in an argillaceous formation located in the Causses Basin in France. The fluid paleotransfers are related to the fault activity to a large extent. We attempt to estimate the intensity of paleo-deviatoric stress magnitudes under which the fault activity may have occurred and consequently, the change in the structural fault permeability. The paleo-deviatoric stress magnitudes were calculated with the inverse method of Etchecopar applied to calcite twinning. The measured crystals are contained within the core zone of minor faults and this study is based on a previous complete microtectonic and microstructural analysis of the faults. In this paper, analysis of calcite twinning has been applied for the first time to vein fillings associated small faults in a context of relatively weak deformation, a condition ensured by the tectonic and structural analysis. Calculation and discussion of the paleo-deviatoric stress tensors in relation to the evolution of the structural fault permeability and to the hydraulic behaviour of the faults are the aim of this paper. The analysed faults, created and active during the same tectonic event, were permeable under a (σ₁–σ₃) mean value of 40–50 MPa. On the other hand, the reactivation of faults during a second tectonic event implies mean (σ₁–σ₃) value higher than 40–50 MPa, especially for the faults that are poorly oriented with respect to the principal tectonic stress directions. The core zone of these faults remained sealed and impermeable or became permeable by development of microcracks inside the pre-existing fillings.     

Iron isotopes in acid mine waters and iron-rich solids from the Tinto–Odiel Basin (Iberian Pyrite Belt, Southwest Spain)

The isotopic composition of Fe was determined in water, Fe-oxides and sulfides from the Tinto and Odiel Basins (South West Spain). As a consequence of sulfide oxidation in mine tailings both rivers are acidic (1.45 < pH < 3.85) and display high concentrations of dissolved Fe (up to 420 mmol l^− 1) and sulphates (up to 1190 mmol l^− 1).The δ⁵⁶Fe of pyrite-rich samples from the Rio Tinto and from the Tharsis mine ranged from − 0.56 ± 0.08‰ to + 0.25 ± 0.1‰. δ⁵⁶Fe values for Fe-oxides precipitates that currently form in the riverbed varied from − 1.98 ± 0.10‰ to 1.57 ± 0.08‰. Comparatively narrower ranges of values (− 0.18 ± 0.08‰ and + 0.21 ± 0.14‰) were observed in their fossil analogues from the Pliocene–Pleistocene and in samples from the Gossan (the oxidized layer that formed through exposure to oxygen of the massive sulfide deposits) (− 0.36 ± 0.12‰ to 0.82 ± 0.07‰). In water, δ⁵⁶Fe values ranged from − 1.76 ± 0.10‰ to + 0.43 ± 0.05‰.At the source of the Tinto River, fractionation between aqueous Fe(III) and pyrite from the tailings was less than would be expected from a simple pyrite oxidation process. Similarly, the isotopic composition of Gossan oxides and that of pyrite was different from what would be expected from pyrite oxidation. In rivers, the precipitation of Fe-oxides (mainly jarosite and schwertmannite and lesser amounts of goethite) from water containing mainly (more than 99%) Fe(III) with concentrations up to 372 mmol l^− 1 causes variable fractionation between the solid and the aqueous phase (− 0.98‰ < Δ⁵⁶Fe_{solid–water} < 2.25‰). The significant magnitude of the positive fractionation factor observed in several Fe(III) dominated water may be related to the precipitation of Fe(III) sulphates containing phases.