Statistical occurrence analysis and spatio-temporal distribution of earthquakes in the Apennines (Italy)

We present two examples of statistical analysis of seismicity conducted by integrating geological, geophysical and seismological data with the aim to characterize the active stress field and to define the spatio-temporal distribution of large earthquakes. Moreover, our data will help to improve the knowledge of the “seismogenic behavior” of the areas and to provide useful information for seismic hazard evaluation.The earthquakes are described by two non-parametric statistical procedures integrating also tectonic-physical parameters to study the spatio-temporal variability.The results show that the areas are characterized by: 1) a stress regime with mainly extensional kinematics; 2) tectonic structures mainly oriented with the active stress field (S_hmin = N44° ± 18° in the southern Apennines and S_hmin = N50° ± 17° in the central Apennines); 3) cluster distribution of seismicity and 4) a high probability of earthquake occurrence (M > 5.5) in the next 10 years.     

Seismotectonics of strike–slip earthquakes within the deep crust of southern Italy: Geometry, kinematics, stress field and crustal rheology of the Potenza 1990–1991 seismic sequences (Mmax 5.7)

We present a revision and a seismotectonic interpretation of deep crust strike–slip earthquake sequences that occurred in 1990–1991 in the Southern Apennines (Potenza area). The revision is motivated by: i) the striking similarity to a seismic sequence that occurred in 2002  140 km NNW, in an analogous tectonic context (Molise area), suggesting a common seismotectonic environment of regional importance; ii) the close proximity of such deep strike–slip seismicity with shallow extensional seismicity (Apennine area); and iii) the lack of knowledge about the mechanical properties of the crust that might justify the observed crustal seismicity. A comparison between the revised 1990–1991 earthquakes and the 2002 earthquakes, as well as the integration of seismological data with a rheological analysis offer new constraints on the regional seismotectonic context of crustal seismicity in the Southern Apennines. The seismological revision consists of a relocation of the aftershock sequences based on newly constrained velocity models. New focal mechanisms of the aftershocks are computed and the active state of stress is constrained via the use of a stress inversion technique. The relationships among the observed seismicity, the crustal structure of the Southern Apennines, and the rheological layering are analysed along a crustal section crossing southern Italy, by computing geotherms and two-mechanism (brittle frictional vs. ductile plastic strength) rheological profiles. The 1990–1991 seismicity is concentrated in a well-defined depth range (mostly between 15 and 23 km depths). This depth range corresponds to the upper pat of the middle crust underlying the Apulian sedimentary cover, in the footwall of the easternmost Apennine thrust system. The 3D distribution of the aftershocks, the fault kinematics, and the stress inversion indicate the activation of a right-lateral strike–slip fault striking N100°E under a stress field characterized by a sub-horizontal N142°-trending σ1 and a sub-horizontal N232°-trending σ3, very similar to the known stress field of the Gargano seismic zone in the Apulian foreland. The apparent anomalous depths of the earthquakes (> 15 km) and the confinement within a relatively narrow depth range are explained by the crustal rheology, which consists of a strong brittle layer at mid crustal depths sandwiched between two plastic horizons. This articulated rheological stratification is typical of the central part of the Southern Apennine crust, where the Apulian crust is overthrusted by Apennine units. Both the Potenza 1990–1991 and the Molise 2002 seismic sequences can be interpreted to be due to crustal E–W fault zones within the Apulian crust inherited from previous tectonic phases and overthrusted by Apennine units during the Late Pliocene–Middle Pleistocene. The present strike–slip tectonic regime reactivated these fault zones and caused them to move with an uneven mechanical behaviour; brittle seismogenic faulting is confined to the strong brittle part of the middle crust. This strong brittle layer might also act as a stress guide able to laterally transmit the deviatoric stresses responsible for the strike–slip regime in the Apulian crust and may explain the close proximity (nearly overlapping) of the strike–slip and normal faulting regimes in the Southern Apennines. From a methodological point of view, it seems that rather simple two-mechanism rheological profiles, though affected by uncertainties, are still a useful tool for estimating the rheological properties and likely seismogenic behaviour of the crust.     

The MW = 6.3 2003 Lefkada earthquake (Greece) and induced stress transfer changes

A large earthquake of magnitude M_W = 6.3 occurred on 14 August 2003 NW of the Lefkada Island, which is situated at the Ionian Sea (western Greece). The source parameters of this event are determined using body-wave modeling. The focal depth was found equal to 9 km, the constrained focal mechanism revealed dextral strike–slip motion (φ = 15°, Δ = 80° and λ = 170°), the duration of the source time function was 8 s and the seismic moment 2.9 × 10²⁵ dyn cm. The earthquake occurred close to the northern end of the Kefallinia transform fault, where the 1994 moderate event and its aftershock sequence were also located. The epicentral distribution of the 2003 aftershock sequence revealed the existence of two clusters. The first one is located close to the epicentral area of the mainshock, while the second southern, close to the northwestern coast of the Kefallinia Island. A gap of seismicity is observed between the two clusters. The length of the activated zone is approximately 60 km. The analysis of data revealed that the northern cluster is directly related to the mainshock, while the southern one was triggered by stress transfer caused by the main event.     

Background seismicity in the Central Apennines of Italy: The Abruzzo region case study

We investigate background seismic activity of the Abruzzo region, a 5000 km² area located within the Central Apennines of Italy, where in the past 600 years at least 5 large earthquakes (I = XI–X) have occurred.Between April 2003 and September 2004, a dense temporary seismic network composed of 30 digital three-component seismic stations recorded 850 earthquakes with 0.9 < M_L < 3.7. We present earthquake locations and focal mechanisms obtained by standard procedures and an optimized velocity model computed with a search technique based on genetic algorithms.The seismicity occurs at a low and constant rate of  2.6 e− 04 events/daykm² and is sparsely distributed within the first 15 km of the crust. Minor increases in the seismicity rate are related to the occurrence of small and localised seismic sequences that occur at the tip of major active normal faults along secondary structures.We observe that during the 16 months of study period, the Fucino fault system responsible for the 1915 Fucino earthquake (M_S = 7.0), and the major normal faults of the area, did not produce significant seismic activity.Fault plane solutions evaluated using P-wave polarity data show the predominance of normal faulting mechanisms ( 55%) with NE-trending direction of extension coherent with the regional stress field active in this sector of the Apennines. Around 27% of the focal solutions have pure strike–slip mechanisms and the rest shows transtensional faulting mechanisms that mainly characterise the kinematics of the secondary structures activated by the small sequences.We hypothesize that the largest known NW-trending normal faults are presently locked and we propose that in the case of activation, the secondary structures located at their tips may act as transfer faults accommodating a minor part of the extensional deformation with strike–slip motion.     

Magnitude of σ1, σ2, and σ3 at mid-crustal levels in an orogenic belt: Microboudin method applied to an impure metachert from Turkey

We revised an equation for estimating palaeostress magnitude using the microboudin technique by incorporating the influence of time on the fracture strength of minerals. The equation was used to estimate triaxial palaeostresses from a rare sample of metachert from Turkey that contains microboudinaged, columnar tourmaline grains in a wide range of orientations within the foliation plane. The estimated principal palaeostresses are σ₁ = 605 MPa, σ₂ = 598 MPa, and σ₃ = 597 MPa. As the microboudinage is considered to have occurred immediately before the rock encountered the brittle-plastic transition during exhumation, these stress values correspond to conditions at approximately 18 km depth and 300 °C within a Cretaceous orogenic belt.     

Geometry of the neotectonic stress field in southern Italy: Geological and seismological evidence

The neotectonic regime in southern Italy has been evaluated by making a comparison between all the available structural and seismological data. The area investigated can be subdivided into four distinct zones which are characterized by different stress regimes. In the Southern Apennines the tensile axis of the stress field is oriented approximately NE-SW while the maximum principal stress (σ₁) is subvertical. In Northern Calabria, the tensile axis is ESE-WNW and the σ₁ axis is almost vertical. In the Catanzaro trough both the tensile axis and the σ₁ axis are subhorizontal and act E-W and N-S, respectively. Finally, the Strait of Messina zone is characterized by a tensile axis oriented E-W and by σ₁ being subvertical.     

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.     

Paleomagnetism of the middle Cretaceous Iritono granite in the Abukuma region, northeast Japan

We have studied the paleomagnetism of the middle Cretaceous Iritono granite of the Abukuma massif in northeast Japan together with ⁴⁰Ar–³⁹Ar dating. Paleomagnetic samples were collected from ten sites of the Iritono granite (102 Ma ⁴⁰Ar–³⁹Ar age) and two sites of its associated gabbroic dikes. The samples were carefully subjected to alternating field and thermal demagnetizations and to rock magnetic analyses. Most of natural remanent magnetizations show mixtures of two components: (1) H component, high coercivity (B_c > 50–90 mT) or high blocking temperature (T_b > 350–560 °C) component and (2) L component, relatively low B_c or low T_b component. H component was obtained from all the 12 sites to give a mean direction of shallow inclination and northwesterly declination (I = 29.9°, D = 311.0°, α₉₅ = 2.7°, N = 12). This direction is different from the geocentric axial dipole field at the present latitude (I = 56.5°) and the typical direction of the Cenozoic remagnetization in northeast Japan. Since rock magnetic properties indicate that the H component of the Iritono granite is carried mainly by magnetite inclusions in plagioclase, this component probably retains a primary one. Thus the shallow inclination indicates that the Abukuma massif was located at a low latitude (16.1 ± 1.6°N) about 100 Ma and then drifted northward by about 20° in latitude. The northwesterly deflection is attributed mostly to the counterclockwise rotation of northeast Japan due to Miocene opening of the Japan Sea. According to this model, the low-pressure and high-temperature (low-P/high-T) metamorphism of the Abukuma massif, which has been well known as a typical location, would have not occurred in the present location. On the other hand, the L component is carried mainly by pyrrhotite and its mean direction shows a moderate inclination and a northwesterly declination (I = 42.8°, D = 311.5°, α₉₅ = 3.3°, N = 9). Since this direction is intermediate between the H component and early Cenozoic remagnetization in northeast Japan, some thermal event would have occurred at lower temperature than pyrrhotite Curie point ( 320 °C) during the middle Cretaceous to early Cenozoic time to have resulted in partial remagnetization.     

Lower Eocene carbonate-cemented “chimney” structures (Varna, Bulgaria) — Control of seepage rates on their formation and stable isotopic signature

The genesis of Lower Eocene calcite-cemented columns, “pisoid”-covered structures and horizontal interbeds, clustered in dispersed outcrops in the Pobiti Kamani area (Varna, Bulgaria) is related to fossil processes of hydrocarbon migration. Field observations, petrography and stable isotope geochemistry of the cemented structures and associated early-diagenetic veins, revealed that varying seepage rates of a single, warm hydrocarbon-bearing fluid, probably ascending along active faults, controlled the type of structure formed and its geochemical signature. Slow seepage allowed methane to oxidize within the sediment under ambient seafloor conditions (δ¹⁸O = − 1 ± 0.5‰ V-PDB), explaining columns' depleted δ¹³C ratios of − 43‰. Increasing seepage rates caused methane to emanate into the water column (δ¹³C = − 8‰) and raised precipitation temperatures (δ¹⁸O = − 8‰). Calcite-cemented conduits formed and upward migrating fluids also affected interbed cementation. Even higher-energy fluid flow and temperatures likely controlled the formation of “pisoids”, whereby sediment was whirled up and cemented.     

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.     

Isothermal decompression history in the “Western Granulite” terrain, central Tanzania: Evidence from reaction textures and trapped fluids in metapelites

The Mozambique Belt (MB) of the East Africa Orogen contains large areas of granulite-facies migmatitic gneisses with Archaean and Palaeoproterozoic protolith ages and that were recycled during the Neoproterozoic Pan-African orogeny. The study area is situated along the Great Ruaha River and within the Mikumi National Park in central Tanzania where migmatitic gneisses and mafic to intermediate granulites are interlayered with Neoproterozoic granulite-facies migmatitic metapelites. Mineral textures suggest isothermal decompression, with the peak mineral assemblage comprising Grt–Bt–Ky–Kfs–Pl–Qtz ± Phn ± Ti-Oxide ± melt and amphibolite-facies retrograde assemblage Grt–Bt–Sil–Ms–Kfs–Pl–Qtz ± Fe–Ti-Oxide. The near isothermal retrograde overprint is seen in well-developed formation of pseudomorphs after garnet. The HP granulite-facies assemblages record P–T conditions of 13–14 kbar at 760–800 °C. Retrogression and the release of fluids from crystallizing melts occurred at 7 kbar and 650–700 °C. A fluid inclusion study shows three types of fluid inclusion consisting of nearly pure CO₂, as well as H₂O–NaCl and H₂O–CO₂ mixtures. We suggest that a immiscible CO₂-bearing brine represents the fluid composition during high-grade peak metamorphism, and that the fluid inclusions containing H₂O–NaCl or nearly pure CO₂ represent trapped fluids from in situ crystallised melt. The results suggest strong isothermal decompression, which is probably related to a fast exhumation after crustal thickening in the central part of the Mozambique Belt in Tanzania.     

Proterozoic–Paleozoic development of the basement of the Central Andes (18–26°S) — a mobile belt of the South American craton

The Late Precambrian–Early Paleozoic metamorphic basement forms a volumetrically important part of the Andean crust. We investigated its evolution in order to subdivide the area between 18 and 26°S into crustal domains by means of petrological and age data (Sm–Nd isochrons, K–Ar). The metamorphic crystallization ages and t_DM ages are not consistent with growth of the Pacific margin north of the Argentine Precordillera by accretion of exotic terranes, but favor a model of a mobile belt of the Pampean Cycle. Peak metamorphic conditions in all scattered outcrop areas between 18 and 26°S are similar and reached the upper amphibolite facies conditions indicated by mineral paragensis and the occurrence of migmatite. Sm–Nd mineral isochrons yielded 525±10, 505±6 and 509±1 Ma for the Chilean Coast Range, the Chilean Precordillera and the Argentine Puna, and 442±9 and 412±18 Ma for the Sierras Pampeanas. Conventional K–Ar cooling age data of amphibole and mica cluster around 400 Ma, but are frequently reset by Late Paleozoic and Jurassic magmatism. Final exhumation of the Early Paleozoic orogen is confirmed by Devonian erosional unconformities. Sm–Nd depleted mantle model ages of felsic rocks from the metamorphic basement range from 1.4 to 2.2 Ga, in northern Chile the average is 1.65±0.16 Ga (1σ; n=12), average t_DM of both gneiss and metabasite in NW Argentina is 1.76±0.4 Ga (1σ; n=22), and the isotopic composition excludes major addition of juvenile mantle derived material during the Early Paleozoic metamorphic and magmatic cycle. These new data indicate a largely similar development of the metamorphic basement south of the Arequipa Massif at 18°S and north of the Argentine Precordillera at 28°S. Variations of metamorphic grade and of ages of peak metamorphism are of local importance. The protolith was derived from Early to Middle Proterozoic cratonic areas, similar to the Proterozoic rocks from the Arequipa Massif, which had undergone Grenvillian metamorphism at ca. 1.0 Ga.     

Solution mechanisms of silicate in aqueous fluid and H2O in coexisting silicate melts determined in-situ at high pressure and high temperature

The structure of H₂O-saturated silicate melts, coexisting silicate-saturated aqueous solutions, and supercritical silicate liquids in the system Na₂O·4SiO₂–H₂O has been characterized with the sample at high temperature and pressure in a hydrothermal diamond anvil cell (HDAC). Structural information was obtained with confocal microRaman and with FTIR microscopy. Fluids and melts were examined along pressure-temperature trajectories defined by the isochores of H₂O at nominal densities, ρ_fluid, (from EOS of pure H₂O) of 0.90 and 0.78 g/cm³. With ρ_fluid = 0.78 g/cm³, water-saturated melt and silicate-saturated aqueous fluid coexist to the highest temperature (800 °C) and pressure (677 MPa), whereas with ρ_fluid = 0.90 g/cm³, a homogeneous single-phase liquid phase exists through the temperature and pressure range (25–800 °C, 0.1–1033 MPa). Less than 5 vol% quartz precipitates near 650 °C in both experimental series, thus driving Na/Si-ratios of melt + fluid phase assemblages to higher values than that of the Na₂O·4SiO₂ starting material.Molecular H₂O (H₂O°) and structurally bonded OH groups were observed in coexisting melts and fluids as well as in supercritical liquids. Their OH/(H₂O)-ratio is positively correlated with temperature. The OH/(H₂O)° in melts is greater than in coexisting fluids. Structural units of Q³, Q², Q¹, and Q⁰ type are observed in all phases under all conditions. An expression of the form, 12Q³ + 13H₂O2Q² + 6Q¹ + 4Q⁰, describes the equilibrium among those structural units. This equilibrium shifts to the right with increasing pressure and temperature with a ΔH of the reaction near 425 kJ/mol.     

Late Paleozoic remagnetization of the Trenton Formation in Ordovician petroleum reservoirs of southwestern Ontario

A paleomagnetic study of subsurface core samples from dolomitized carbonates of two producing reservoirs in the Upper Ordovician Trenton Formation, collected from four wells in southwestern Ontario yielded a paleomagnetic direction of D = 152.3°, I = − 12.3° (N = 49, α₉₅ = 8.7). This characteristic remanent magnetization (ChRM) direction was azimuth-corrected by aligning the viscous remanence magnetization (VRM) with the present Earth's magnetic field direction. A drilling-induced magnetization (VRMdi) was present in less than half the specimens sampled in this study. In addition, where the VRM correction could not be made, a paleolatitudinal arc calculated from the inclination-only mean of I = − 9.0° (N = 34, α₉₅ = 3.0°) intersected the apparent polar wander path in the Late Permian–Early Triassic. These paleodirections are similar to the paleomagnetic directions observed in Ordovician Trenton carbonates from the Michigan Basin and New York State, U.S.A., suggesting a related regional late Paleozoic remagnetization.     

Evaluation of earthquake potential along the Northern Anatolian Fault Zone in the Marmara Sea using comparisons of GPS strain and seismotectonic parameters

Seismotectonic parameters including the Gutenberg-Richter b-value and multifractal dimensions D₂ and D₁₅ of seismicity patterns (both spatial and temporal) were compared to GPS-derived maximum shear and dilatation strains measured in the Marmara Sea region of western Turkey along the Northern Anatolian Fault Zone (NAFZ). Comparisons of seismotectonic parameters and GPS-derived maximum shear and dilatation strain along the NAFZ in the vicinity of the 1999 M7.4 Izmit earthquake reveal a positive correlation (r = 0.5, p = 0.05) between average dilatation and the Gutenberg-Richter b-value. Significant negative correlation (r = − 0.56, p = 0.03 and r = − 0.56, p = 0.02) was also observed between the spatial fractal dimension D₂ and GPS-derived maximum geodetic and shear strain. This relationship suggests that, as maximum geodetic and shear strains increase, seismicity becomes increasingly clustered.Anomalous interrelationships are observed in the Marmara Sea region prior to the Izmit event along a bend in the NAFZ near the eastern end of the Marmara Sea known as the Northern Boundary Fault (NBF). An asperity is located near the northwest end of the NBF. Along the 50-km length of the NBF, GPS strains become slightly compressive. The correlation between b-value and GPS-derived dilatation suggests that regions in compression have increased probability of larger magnitude rupture. The NBF appears to serve as an impediment to the transfer of strain from east to west along the NAFZ. Recurrence times for large earthquakes along the NBF are larger than in surrounding areas. Temporal clustering of seismicity in the vicinity of the NBF may represent foreshocks of an impending rupture.     

Palaeomagnetic results from the Early Permian Copacabana Group, southern Peru: Implication for Pangaea palaeogeography

Samples collected from folded carbonate rocks of the Early Permian Copacabana Group exposed in the Peruvian Subandean Zone have been subjected to detailed palaeomagnetic analysis. Thermal demagnetisation of most samples yield stable high unblocking temperature directions dominantly carried by titanomagnetite minerals. This remanence, identified in 32 samples (43 specimens), is exclusively of reverse polarity consistent with the Permian–Carboniferous Reversal Superchron (PCRS). The overall directions pass the fold test at the 99% confidence level and are considered as being a pre-folding remanence acquired in Early Permian times. The Copacabana Group yields an overall mean direction of D = 166°, I = +49° (α₉₅ = 4.5°, k = 131.5, N = 9 sites) in stratigraphic coordinates and a corresponding palaeosouth pole position situated at λ = 68°S,  = 321°E (A₉₅ = 5.2°, K = 100). Combining this pole with the coeval high quality data from South America, Africa and Australia results in a mean pole for Gondwana situated at λ = 34.4°S,  = 065.6°E (A₉₅ = 4.9°, K = 73.6, N = 13 studies) in African coordinates. This pole position supports a Pangaea B palaeogeography in Early Permian times. In contrast, the combined pole for Gondwana diverges from the coeval Laurasian mean pole when assuming the Pangaea A-type configuration. Poor quality of the Gondwana dataset and inclination shallowing in sediments seem to play no role in the misfit between the Permian–Triassic poles from Gondwana and Laurasia in Pangaea A reconstruction.     

Magnetic studies of magma-supply and sea-floor metamorphism: Troodos ophiolite dikes

Dikes of the eastern Troodos ophiolite of Cyprus intruded at slow ocean-spreading axes with dips ranging up to 15° from vertical and with bimodal strikes (now NE–SW and N–S due to post-88 Ma sinistral microplate rotation). Varied dike orientations may represent local stress fields during dike-crack propagation but do not influence the spatial-distributions or orientation-distributions of dikes' magnetic fabrics, nor of their palaeomagnetic signals. Anisotropy of magnetic susceptibility (AMS) integrates mineral orientation-distributions from each of 1289 specimens sampled from dikes at 356 sites over 400 km² in the eastern Troodos ophiolite of Cyprus. In 90% of dikes, AMS fabrics define a foliation (k_MAX–k_INT) parallel to dike walls and a lineation (k_MAX) that varies regionally and systematically. Magma-flow alignment of accessory magnetite controls the AMS with a subordinate contribution from the mafic silicate matrix that is reduced in anisotropy by sea-floor metamorphism. Titanomagnetite has less influence on anisotropy. Occasionally, intermediate and minimum susceptibility axes are switched so as to be incompatible with the kinematically reasonable flow plane but maximum susceptibility (k_MAX) still defines the magmatic flow axis. Such blended subfabrics of kinematically compatible mafic-silicate and misaligned multidomain magnetite subfabrics; are rare. Areas of steep magma flow (k_MAX plunge ≥ 70°) and of shallow magma-flow alternate in a systematic and gradual spatial pattern. Foci of steep flow were spaced 4 km parallel to the spreading axes and 6 km perpendicular to the spreading axes. Ridge-parallel separation of steep flow suggest the spacing of magma-feeders to the dikes whereas ridge-perpendicular spacing of 6 km at a spreading rate of 50 mm/a implies the magma sources may have been active for 240 Ka. The magma feeders feeding dikes may have been ≤ 2 km in diameter. Stable paleomagnetic vectors, in some cases verified by reversal tests, are retained by magnetite and titanomagnetite. In all specimens, the stable components were isolated by three cycles of low-temperature demagnetization (LTD) followed by ≥ 10 steps of incremental thermal demagnetization (TD). 47% of primary A-components [338.2 /+ 57.2 n = 207, α₉₅ = 3.9; mean T_UB = 397 ± 8 °C] are overprinted by a B-component [341.4 /+ 63.5, n = 96, α₉₅ = 8.7; mean T_UB = 182 ± 11 °C]. A- and B-components are ubiquitous and shared equally by the N–S and NE–SW striking dikes. A-component unblocking temperatures (T_UB) are zoned subparallel to the fossil spreading axis. Their spatial pattern is consistent with chemical remagnetization at some certain off-axis distance determined by sea-floor spreading. A-components indicate less microplate rotation and more northerly palaeolatitudes that are consistent with metamorphic remagnetization after some spreading from the ridge-axis. Thus, their magnetizations are younger than those of the overlying volcanic sequence for which ChRMs are commonly reported as 274 /+ 33 (88 Ma).     

The importance of the effective intermediate principal stress (σ′2) to fault slip patterns

In normal faulting regimes, the magnitudes and orientations of the maximum and minimum principal compressive stresses may be known with some confidence. However, the magnitude of the intermediate principal compressive stress is generally much more difficult to constrain and is often not considered to be an important factor. In this paper, we show that the slip characteristics of faults and fractures with complex or nonoptimal geometry are highly sensitive to variation or uncertainty in the ambient effective intermediate principal stress (σ′₂). Optimally oriented faults and fractures may be less sensitive to such variations or uncertainties. Slip tendency (T_s) analysis provides a basis for quantifying the effects of uncertainty in the magnitudes and orientations of all principal stresses and in any stress regime, thereby focusing efforts on the most important components of the system. We also show, for a normal faulting stress regime, that the proportion of potential surfaces experiencing high slip tendency (e.g., T_s ≥ 0.6) decreases from a maximum of about 38% where σ′₂ = σ′₃, to a minimum of approximately 14% where σ′₂ is halfway between σ′₃ and σ′₁, and increases to another high of approximately 29% where σ′₂ = σ′₁. This analysis illustrates the influence of the magnitude of σ′₂ on rock mass strength, an observation previously documented by experimental rock deformation studies. Because of the link between fault and fracture slip characteristics and transmissivity in critically stressed rock, this analysis can provide new insights into stress-controlled fault transmissivity.     

Paleomagnetism of the Fort Knox Stock, Alaska, and rotation of the Yukon–Tanana terrane after 92.5 Ma

The 92.5 Ma Fort Knox granodiorite stock, near the western end of the Fairbanks Belt in the Yukon–Tanana terrane (YTT) of central Alaska, hosts a world-class gold mine. The stock has been analysed paleomagnetically using thermal and alternating-field step demagnetization and isothermal remanence methods. This pluton retains a primary thermoremanent magnetization at 18 sites (232 specimens) that resides mainly in single-to pseudosingle-domain magnetite with a direction of D = 228.8°, I = 84.3° (N = 18, k = 130, α₉₅ = 3.0°), giving a paleopole at 56.5°N, 197.1°E (dp = 5.9°, dm = 5.8°). The pluton's host rock, the Fairbanks schist, does not retain a stable coherent remanence. Relative to the North American craton, the stock's paleoinclination indicates that the Fairbanks Belt has undergone nonsignificant poleward (northwesterly) translation of 25 ± 750 km only. Analysed in concert with the few available paleoinclinations available for the YTT in Yukon, the paleoinclination suggests further that the YTT has undergone only  250 to 450 km of dextral displacement along the Tintina fault in the past  100 Ma and, therefore, is parautocthonous since the mid-Cretaceous. The stock's paleodeclination records 121 ± 35° of counterclockwise rotation relative to the North American craton. Consideration of models published for Alaska's tectonic evolution suggests that this paleodeclination discordance is caused by rotations associated with the opening of the Canada Basin, with dextral displacement on the Tintina fault, and with development of the western Alaskan orocline. Thus the paleomagnetic results for the Fort Knox stock support a thin-skin tectonic model for the accretion of the YTT and Intermontane Belt terranes to the northern Cordillera.     

Solid inclusions of magmatic halite and sylvite in felsic granitoids, Sierra Norte, Córdoba, Argentina

This study provides evidence for the existence of halite and sylvite solid inclusions in igneous quartz and feldspars, the first to be reported in intrusive rocks, and to partially constrain the physicochemical environment that lets halides crystallize under magmatic conditions.Halite and sylvite solid inclusions were found included in quartz and feldspars from a micrographic–granophyric assemblage in a miarolitic aplite and, rarer, in alkali-feldspar from a miarolitic monzogranite. Monzogranite and aplite represent I-type, K-enriched postcollisional rocks of the Late Cambrian–Early Ordovician Sierra Norte–Ambargasta batholith in the Eastern Sierras Pampeanas. Both granitoids fall among the most evolved felsic rocks of the batholith, with aplite approaching haplogranitic compositions. Halite is far more common than sylvite and the presence and distribution of one or both halides are erratic within the felsic intrusive bodies. Halides occur as small skeletal grains, commonly in cross-shaped aggregates of less than 50 μm. No K or Na was found at the detection limits of EDS in either halite or sylvite respectively. Textural relationships suggest that the alkali-chlorides separated from the melt near the minima along the quartz–feldspar cotectics of P_H2O > 160 < 200 MPa in a silica-, and potassium-rich magmatic system at approximately 750–700 °C, prior to the H₂O-vapor saturated miarole-forming stage.Computed ratios for the magmatic volatile phase (MVP) coexisting with melt at the early stage of aplite crystallization are: NaCl/HCl = 0.11–0.97 and KCl/HCl = 0.24–1.62, being the highest range of values (0.79–0.97 and 1.45–1.62, respectively) found in those alkali-chloride-bearing samples. Maximum HCl/ΣCl_(MVP) (0.28 to 0.31) indicates higher total Cl concentration in the MVP of alkali-chloride-bearing aplites, which is much higher in the halite-free aplite samples (HCl/ΣCl_(MVP) = 0.59 to 0.74). One miarolitic monzogranite sample, where halite solid inclusions are present, also yielded the highest ratios for NaCl/HCl_(MVP) (0.91) and KCl/HCl_(MVP) (1.46), and the HCl/ΣCl_(MVP) is 0.30. A high HCl concentration in the fluid phase is suggested by the log f(HF)/f(H₂O) = − 4.75 to − 4.95, log f(HCl)/f(H₂O) = − 3.73 to − 3.86, and log f(HF)/f(HCl) = − 0.88 to − 1.22, computed at 750 °C after biotite composition. The Cl concentrations at 800 °C, computed with a D^v/l_Cl = 0.84 + 26.6P (P at 200 MPa), yielded values within the range of  70 to 700 ppm Cl in the melt and  4000 to 40 000 ppm Cl in the coexisting MVP. The preferential partitioning of Cl in the vapor phase is controlled by the D^v/l_Cl; however, the low concentration of Cl in the melt suggests that high concentrations of Cl are not necessary to saturate the melt in NaCl or KCl.Cl-saturation of the melt and coexisting MVP might have been produced by a drop in Cl solubility due to the near-haplogranitic composition of the granitoids after extreme fractionation, probably enhanced by fluctuating reductions of the emplacement pressure in the brittle monzogranite host. Liquid immiscibility, based in the differential viscosity and density among alkali-chloride saturated hydrosaline melt, aluminosilicate felsic melt, and H₂O-rich volatiles is likely to have crystallized halite and sylvite from exsolved hydrosaline melt. High degrees of undercooling might have been important at the time of alkali-chloride exsolution. The effectiveness of alkali-chloride separation from the melt at magmatic temperatures is in line with the interpretation of “halite subtraction” as a necessary process to understand the origin of the “halite trend” in highly saline fluid inclusions from porphyry copper and other hydrothermal mineralizations, despite the absence of the latter in the Cerro Baritina aplites, where this process preceded the exsolution of halite-undersaturated fluids.Pervasive alteration of the monzogranite country rock as alkali-metasomatic mineral assemblages, the mineral chemistry of some species, and the association of weak molybdenite mineralization are compatible with the activity of alkaline hypersaline fluids, most likely exsolved during the earliest stages of aplite consolidation.