Coulomb stress triggering of earthquakes along the Atalanti Fault, central Greece: Two April 1894 M6+ events and stress change patterns

Two M6+ events occurred 15–20 km apart in central Greece on April 20 and April 27, 1894. We identify the April 27, 1894 rupture (2nd in the sequence) with the Atalanti segment of the Atalanti Fault Zone because of unequivocal surface rupturing evidence reported by Skouphos [Skouphos, T., 1894. Die swei grossen Erdbeben in Lokris am 8/20 und 15/27 April 1894. Zeitschrift Ges. Erdkunde zu Berlin, vol. 24, pp. 409–474]. Coulomb stress transfer analysis and macroseismic evidence suggest that the April 20, 1894 event (1st in the sequence) may be associated with the Martinon segment of the same fault zone. Our stress modelling suggests that this segment may have ruptured in an M = 6.4 event producing a 15-km long rupture which transferred 1.14 bar in the epicentral area of the April 27th, 1894 event, thus triggering the second M = 6.6 earthquake along the Atalanti segment and producing a 19-km long rupture. We also examined three alternative fault sources for the first event; however, all these produce smaller stress stresses for triggering the second event. The proposed slip model for the second earthquake is capable of producing coastal subsidence of the order of centimetres to decimetres, which fits the geological data. The 1894 earthquake sequence was followed by a difference in the timing of subsequent M > 5 events in each of the “relaxed” areas (stress shadows; a negative change in Coulomb failure stress > − 0.6 bar), which terminated between 22–37 years (north) and 80 years (south).     

Petrology and mineralogy of granulite-facies mafic xenoliths (Sardinia, Italy): Evidence for KCl metasomatism in the lower crust

Here new mineralogical data is presented on the occurrence of K-feldspar in granulite-facies metagabbronorite xenoliths found in recent alkaline lavas from Western Sardinia, Italy. The xenoliths originated from the underplating of variably evolved subduction-related basaltic liquids, which underwent cooling and recrystallisation in the deep crust (T = 850–900 °C, P = 0.8–1.0 GPa). They consist of orthopyroxene + clinopyroxene + plagioclase porphyroclasts (An = 50–66 mol%) in a granoblastic, recrystallised, quartz-free matrix composed of pyroxene + plagioclase (An = 56–72 mol%) + Fe–Ti oxides ± K-feldspar ± biotite ± fluorapatite ± Ti-biotite. Texturally, the K-feldspar occurs in a variety of different modes. These include: (1) rods, blebs, and irregular patches in a random scattering of plagioclase grains in the form of antiperthite; (2) micro-veins along plagioclase–plagioclase and plagioclase–pyroxene grain rims; (3) myrmekite-like intergrowths with Ca-rich plagioclase along plagioclase–plagioclase grain boundaries; and (4) discrete anhedral grains (sometimes microperthitic). The composition of each type of K-feldspar is characterized by relatively high albite contents (16–33 mol%). An increasing anorthite content in the plagioclase towards the contact with the K-feldspar micro-vein and myrmekite-like intergrowths into the K-feldspar along the plagioclase–K-feldspar grain boundary are also observed. Small amounts of biotite (TiO₂ = 4.7–6.5 wt.%; F = 0.24–1.19 wt.%; Cl = 0.04–0.20 wt.%) in textural equilibrium with the granulite-facies assemblage is present in both K-feldspar-bearing and K-feldspar-free xenoliths. These K-feldspar textures suggest a likely metasomatic origin due to solid-state infiltration of KCl-rich fluids/melts. The presence of such fluids is supported by the fluorapatite in these xenoliths, which is enriched in Cl (Cl = 6–50% of the total F + Cl + OH). These lines of evidence suggest that formation of K-feldspar in the mafic xenoliths reflects metasomatic processes, requiring an external K-rich fluid source, which operated in the lower crust during and after in-situ high-T recrystallisation of relatively dry rocks.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

Ga/Mg ratio as a new geochemical tool to differentiate magmatic from metamorphic blue sapphires

Using ICP-MS–LA analyses, we demonstrate that the use of the Ga/Mg ratio, in conjunction with the Fe concentration, is an efficient tool in discriminating between “metamorphic” and “magmatic” blue sapphires. Magmatic blue sapphires found in alkali basalts (e.g. southeastern Asia, China, Africa) are commonly medium-rich to rich in Fe (with average contents between 2000 and 11000 ppm), high in Ga (> 140 ppm), and low in Mg (generally < 20 ppm) with high Ga/Mg ratios (> 10). Conversely, metamorphic blue sapphires found in basalts (e.g. Pailin pastel) and in metamorphic terrains (e.g. Mogok, Sri Lanka, Ilakaka) are characterized by low average iron contents (< 3000 ppm), low Ga contents (< 75 ppm), and high Mg values (> 60 ppm) with low average Ga/Mg ratios (< 10). Basaltic magmatic sapphires have Fe, Ga and Mg contents similar to those obtained for primary magmatic sapphires found in the Garba Tula syenite. This suggests that these both sets of sapphires have a possible common “syenitic” origin, as previously proposed from other criteria. In addition, plumasite-related sapphires and metamorphic sapphires also exhibit similar composition in trace elements. Based on results from the present study, we suggest that fluid circulations during a metamorphic stage produced metasomatic exchanges between mafic and acidic rocks (plumasite model), thus explaining the high Mg contents and converging Ga/Mg ratios observed in metamorphic sapphires.     

Testing a new hybrid approach to seismic hazard assessment: an application to the Calabrian Arc (Southern Italy)

We tested a new hybrid method for the evaluation of seismic hazard. A recently proposed fault segmentation and earthquake recurrence model of peninsular Italy suggests that the interval for which the local historical catalogue is complete is shorter than the mean recurrence time of individual large faults (1000 years), or at the most comparable. These new findings violate the fundamental assumption of historical probabilistic seismic hazard methods that the historical record is representative of the activity of all the seismogenic sources. The hybrid method we propose uses time-dependent modelling of the major earthquakes and catalogue-based historical probabilistic estimates for all minor events. We assume that the largest earthquakes are characteristic for individual discrete fault segments, model their probability of occurrence by a renewal process and compute the shaking associated with each of them with a simplified procedure. Then we calculate the probability of exceeding a given threshold of peak ground acceleration for specific sites as the aggregate probability of occurrence of large characteristic earthquakes and minor shocks. We apply the method to the Calabrian Arc (Southern Italy) performing the calculations for five major towns. The exposure to seismic hazard of Reggio Calabria, Catanzaro and Vibo Valentia, which locate close to recently activated large faults, decreases with respect to traditional time-independent estimates. On the contrary, an increase of seismic hazard is obtained for Castrovillari, which locates in an area where large faults displaying Holocene activity have been recently recognized but no significant earthquake is reported in the historical catalogue. Cosenza has the highest probability to experience a significant peak ground acceleration with both the new hybrid and the traditional approaches. We wish to stress that the present results should be interpreted only in terms of the differences between the new hybrid and the traditional approaches, not for their absolute values, and that they are not intended to be used for updating or modifying the current national seismic zonation.     

Regional time-predictable modeling in the Hindukush–Pamir–Himalayas region

The seismic characteristic of Hindukush–Pamir–Himalaya (HPH) and its vicinity is very peculiar and has experienced many widely distributed large earthquakes. Recent work on the time-dependent seismicity in the Hindukush–Pamir–Himalayas is mainly based on the so-called “regional time-predictable model”, which is expressed by the relation log T=cM_p+a, where T is the inter-event time between two successive main shocks of a region and M_p is the magnitude of the preceded main shock. Parameter a is a function of the magnitude of the minimum earthquake considered and of the tectonic loading and c is positive (0.3) constant. In 90% of the cases with sufficient data, parameter c was found to be positive, which strongly supports the validity of the model. In the present study, a different approach, which assumes no prior regionalization of the area, is attempted to check the validity of the model. Nine seismic sources were defined within the considered region and the inter-event time of strong shallow main shock were determined and used for each source in an attempt at long-term prediction, which show the clustering and occurrence of at least three earthquakes of magnitude 5.5≤M_s≤7.5 giving two repeat times, satisfying the necessary and sufficient conditions of time-predictable model (TP model). Further, using the global applicability of the regional time- and magnitude-predictable model, the following relations have been obtained: log T_t=0.19 M_min+0.52M_p+0.29 log m₀−10.63 and M_f=1.31M_min−0.60M_p−0.72 log m₀+21.01, where T_t is the inter-event time, measured in years; M_min the surface wave magnitude of the smallest main shock considered; M_p the magnitude of preceding main shock; M_f the magnitude of the following main shock; and m₀ the moment rate in each source per year.

These relations may be used for seismic hazard assessment in the region. Based on these relations and taking into account the time of occurrence and the magnitude of the last main shock in each seismogenic source, time-dependent conditional probabilities for the occurrence of the next large (M_s≥5.5) shallow main shocks during the next 20 years as well as the magnitudes of the expected main shocks are determined.     

Petrogenesis and tectonic setting of the peralkaline Pine Canyon caldera, Trans-Pecos Texas, USA

The Pine Canyon caldera is a small (6–7 km diameter) ash-flow caldera that erupted peralkaline quartz trachyte, rhyolite, and high-silica rhyolite lavas and ash-flow tuffs about 33–32 Ma. The Pine Canyon caldera is located in Big Bend National Park, Texas, USA, in the southern part of the Trans-Pecos Magmatic Province (TPMP). The eruptive products of the Pine Canyon caldera are assigned to the South Rim Formation, which represents the silicic end member of a bimodal suite (with a “Daly Gap” between 57 and 62 wt.% SiO₂); the mafic end member consists primarily of alkali basalt to mugearite lavas of the 34–30 Ma Bee Mountain Basalt. Approximately 60–70% crystallization of plagioclase, clinopyroxene, olivine, magnetite, and apatite from alkali basalt coupled with assimilation of shale wall rock (M_a/M_c = 0.3–0.4) produced the quartz trachyte magma. Variation within the quartz trachyte–rhyolite suite was the result of 70% fractional crystallization of an assemblage dominated by alkali feldspar with subordinate clinopyroxene, fayalite, ilmenite, and apatite. High-silica rhyolite is not cogenetic with the quartz trachyte–rhyolite suite, and can be best explained as the result of  5% partial melting of a mafic granulite in the deep crust under the fluxing influence of fluorine. Variation within the high-silica rhyolite is most likely due to fractional crystallization of alkali feldspar, quartz, magnetite, biotite, and monazite. Lavas and tuffs of the South Rim Formation form A-type rhyolite suites, and are broadly similar to rock series described in anorogenic settings both in terms of petrology and petrogenesis. The Pine Canyon caldera is interpreted to have developed in a post-orogenic tectonic setting, or an early stage of continental rifting, and represents the earliest evidence for continental extension in the TPMP.     

Fluid-controlled crustal metasomatism within a high-pressure subducted mélange (Mt. Hochwart, Eastern Italian Alps)

The Nonsberg–Ultental Region of northern Italy contains a Palaeozoic mélange that was partially subducted during the Variscan orogeny. This mélange is constituted mainly by metapelites characterized by shale-type REE-patterns, displaying partial melting which began under high-pressure conditions. The resulting migmatites enclose minor slivers of mantle-wedge peridotites that have been incorporated into the mélange during subduction. Peridotites display important large ion lithophile elements (LILE) enrichment consequent to amphibole recrystallization contemporaneously with metapelite migmatization at P ≈ 2.7 GPa and T ≈ 850 °C in the garnet–peridotite field. Crustal and mantle (ultramafic) rocks of the mélange display the same Sm–Nd ages of about 330 ± 6 Ma, which dates both the metamorphic peak and the migmatization event. The zircon U–Pb age of the metasomatic amphibolitic contact between garnet peridotite and migmatite is identical (333.3 ± 2.4 Ma) within analytical errors. Therefore, metasomatism, migmatization and peak metamorphism are constrained to the same event. The presence of Cl-rich apatite and ferrokinoshitalite in the contact amphibolite, together with the trace-element patterns of peridotites, suggest that metasomatism was driven by Cl- and LILE-rich fluids derived from ocean water transported into the subduction zone by sediments and crustal rocks. These fluids interacted with the crust, prompting partial melting under water oversaturated conditions and partitioning LILE from the crust itself. Peridotites, which were well below their wet solidus temperature, could not melt but they recrystallized in the crustal mélange under garnet-facies conditions. Crustal fluids caused extensive hydration and LILE-enrichment in peridotites and severe Sm–Nd isotope disequilibrium between minerals, especially in the recrystallized peridotites. The proposed scenario suggests massive entrapment of crustal aqueous fluids at high-pressure conditions within subduction zones.     

Diamonds and their mineral inclusions from the A154 South pipe, Diavik Diamond Mine, Northwest territories, Canada

Mineral inclusions recovered from 100 diamonds from the A154 South kimberlite (Diavik Diamond Mines, Central Slave Craton, Canada) indicate largely peridotitic diamond sources (83%), with a minor (12%) eclogitic component. Inclusions of ferropericlase (4%) and diamond in diamond (1%) represent “undetermined” parageneses.

Compared to inclusions in diamonds from the Kaapvaal Craton, overall higher CaO contents (2.6 to 6.0 wt.%) of harzburgitic garnets and lower Mg-numbers (90.6 to 93.6) of olivines indicate diamond formation in a chemically less depleted environment. Peridotitic diamonds at A154 South formed in an exceptionally Zn-rich environment, with olivine inclusions containing more than twice the value (of  52 ppm) established for normal mantle olivine. Harzburgitic garnet inclusions generally have sinusoidal rare earth element (REE_N) patterns, enriched in LREE and depleted in HREE. A single analyzed lherzolitic garnet is re-enriched in middle to heavy REE resulting in a “normal” REE_N pattern. Two of the harzburgitic garnets have “transitional” REE_N patterns, broadly similar to that of the lherzolitic garnet. Eclogitic garnet inclusions have normal REE_N patterns similar to eclogitic garnets worldwide but at lower REE concentrations.

Carbon isotopic values (δ¹³C) range from − 10.5‰ to + 0.7‰, with 94% of diamonds falling between − 6.3‰ and − 4.0‰. Nitrogen concentrations range from below detection (< 10 ppm) to 3800 ppm and aggregation states cover the entire spectrum from poorly aggregated (Type IaA) to fully aggregated (Type IaB). Diamonds without evidence of previous plastic deformation (which may have accelerated nitrogen aggregation) typically have < 25% of their nitrogen in the fully aggregated B-centres. Assuming diamond formation beneath the Central Slave to have occurred in the Archean [Westerlund, K.J., Shirey, S.B., Richardson, S.H., Gurney, J.J., Harris, J.W., 2003b. Re–Os systematics of diamond inclusion sulfides from the Panda kimberlite, Slave craton. VIIIth International Kimberlite Conference, Victoria, Canada, Extended Abstracts, 5p.], such low aggregation states indicate mantle residence at fairly low temperatures (< 1100 °C). Geothermometry based on non-touching inclusion pairs, however, indicates diamond formation at temperatures around 1200 °C. To reconcile inclusion and nitrogen based temperature estimates, cooling by about 100–200 °C shortly after diamond formation is required.     

Provenance of the Neoproterozoic Maricá Formation (Sul-rio-grandense Shield, Southern Brazil): Petrographic and Sm–Nd isotopic constraints

This paper reports the integrated application of petrographic and Sm–Nd isotopic analyses for studying the provenance of the Neoproterozoic Maricá Formation, southern Brazil. This unit encompasses sedimentary rocks of fluvial and marine affiliations. In the lower fluvial succession, sandstones plot in the “craton interior” and “transitional continental” fields of the QFL diagram. Chemical weathering probably caused the decrease of the ¹⁴⁷Sm/¹⁴⁴Nd ratios to 0.0826 and 0.0960, consequently lowering originally > 2.0 Ga T_DM ages to 1.76 and 1.81 Ga. ¹⁴³Nd/¹⁴⁴Nd ratios are also low (0.511521 to 0.511633), corresponding to negative ε_Nd present-day values (− 21.8 and − 19.6). In the intermediate marine succession, sandstones plot in the “dissected arc” field, reflecting the input of andesitic clasts. Siltstones and shales reveal low ¹⁴³Nd/¹⁴⁴Nd ratios (0.511429 to 0.511710), ε_Nd values of − 18.1 and − 23.6, and T_DM ages of 2.16 and 2.37 Ga. Sandstones of the upper fluvial succession have “dissected arc” and “recycled orogen” provenance. ¹⁴³Nd/¹⁴⁴Nd isotopic ratios are also relatively low, from 0.511487 to 0.511560, corresponding to ε_Nd values of − 22.4 and − 21.0 and T_DM of 2.07 Ga. A uniform granite–gneissic basement block of Paleoproterozoic age, with subordinate volcanic rocks, is suggested as the main sediment source of the Maricá Formation.     

Decelerating precursory seismicity in Vrancea

Preshock seismic excitation followed by seismic quiescence has been observed in the seismogenic region of strong shallow mainshocks. The strain released by such preshocks is decelerating with the time to the mainshock and is fitted by a power-law with a power value larger than unit. This model is tested in the present work for the intermediate-depth earthquakes of the Vrancea region, generated in an isolated seismogenic zone proper for such testing. A backward application of this “decelerating preshock strain” model for the case of 4.3.1977 (M = 7.5) earthquake, for which reliable data are available, shows a good fit of the power-law pattern to the seismic activity preceding the main shock. The occurrence rate of recent intermediate-depth shocks in Vrancea indicates that this region is currently in a state of decelerating seismic deformation, which may lead to the generation of a strong intermediate-depth mainshock there at about the beginning of the third decade of the present century. The respective uncertainties are unknown due to lack of previous relative studies.     

Perspectives on basaltic magma crystallization and differentiation: Lava-lake blocks erupted at Mauna Loa volcano summit, Hawaii

Explosive eruptions at Mauna Loa summit ejected coarse-grained blocks (free of lava coatings) from Moku'aweoweo caldera. Most are gabbronorites and gabbros that have 0–26 vol.% olivine and 1–29 vol.% oikocrystic orthopyroxene. Some blocks are ferrogabbros and diorites with micrographic matrices, and diorite veins (≤ 2 cm) cross-cut some gabbronorites and gabbros. One block is an open-textured dunite.

The MgO of the gabbronorites and gabbros ranges  7–21 wt.%. Those with MgO > 10 wt.% have some incompatible-element abundances (Zr, Y, REE; positive Eu anomalies) lower than those in Mauna Loa lavas of comparable MgO; gabbros (MgO < 10 wt.%) generally overlap lava compositions. Olivines range Fo_83–58, clinopyroxenes have Mg#s  83–62, and orthopyroxene Mg#s are 84–63 — all evolved beyond the mineral-Mg#s of Mauna Loa lavas. Plagioclase is An_75–50. Ferrogabbro and diorite blocks have  3–5 wt.% MgO (TiO₂ 3.2–5.4%; K₂O 0.8–1.3%; La 16–27 ppm), and a diorite vein is the most evolved (SiO₂ 59%, K₂O 1.5%, La 38 ppm). They have clinopyroxene Mg#s 67–46, and plagioclase An_57–40. The open-textured dunite has olivine  Fo_83.5. Seven isotope ratios are ⁸⁷Sr/⁸⁶Sr 0.70394–0.70374 and ¹⁴³Nd/¹⁴⁴Nd 0.51293–0.51286, and identify the suite as belonging to the Mauna Loa system.

Gabbronorites and gabbros originated in solidification zones of Moku'aweoweo lava lakes where they acquired orthocumulate textures and incompatible-element depletions. These features suggest deeper and slower cooling lakes than the lava lake paradigm, Kilauea Iki, which is basalt and picrite. Clinopyroxene geobarometry suggests crystallization at < 1 kbar P. Highly evolved mineral Mg#s, < 75, are largely explained by cumulus phases exposed to evolving intercumulus liquids causing compositional ‘shifts.’ Ferrogabbro and diorite represent segregation veins from differentiated intercumulus liquids filter pressed into rigid zones of cooling lakes. Clinopyroxene geobarometry suggests < 300 bar P. Open-textured dunite represents olivine-melt mush, precursor to vertical olivine-rich bodies (as in Kilauea Iki). Its Fo_83.5 identifies the most primitive lake magma as  8.3 wt.% MgO. Mass balancing and MELTS show that such a magma could have yielded both ferrogabbro and diorite by ≥ 50% fractional crystallization, but under different fO₂: < FMQ (250 bar) led to diorite, and FMQ (250 bar) yielded ferrogabbro. These segregation veins, documented as similar to those of Kilauea, testify to appreciable volumes of ‘rhyolitic’ liquid forming in oceanic environments. Namely, SiO₂-rich veins are intrinsic to all shields that reached caldera stage to accommodate various-sized cooling, differentiating lava lakes.     

Magnitude-frequency distribution in the European-Mediterranean earthquake regions

The European-Mediterranean earthquake catalogue from 1901 to 1985, which comprises uniformly determined magnitudes M_S and m_B(h ≥ 60 km) of 13300 events, was used in the study of cumulative magnitude-frequency relationships N_c(M) compiled for 75 earthquake regions and 25 larger provinces. In the whole magnitude range observed, the Gutenberg-Richter formula log N_c(M) = a − bM very rarely fits the cumulative (log N_c, M) distributions. The b-values of log-linear segments of N_c(M) vary regionally from b = 0.7 to b = 1.3; averaging of all values leads to (shallow events, M_S and ).

Most distributions pertain to the Mediterranean area (b = 0.86 from the graph for shallow events) and many of them indicate the existence of characteristic earthquakes in accordance with the theoretical single-fault model. Other observed shapes of N_c(M) can be explained by the superposition of populations of different M_max values or by the presence of swarm-type activity. The observed N_c(M) distributions depend very much on the delineation of earthquake regions i.e. on the number and dimension of seismoactive faults in the investigated region.

A premonitory enhancement of medium earthquake activity (M = 4.5–5.5) can be observed only very rarely.     

Spinel–sapphirine–quartz bearing composite inclusion within garnet from an ultrahigh-temperature pelitic granulite: Implications for metamorphic history and P–T path

We report here a multiphase mineral inclusion composed of quartz, plagioclase, K-feldspar, sapphirine, spinel, orthopyroxene, and biotite, in porphyroblastic garnet within a pelitic granulite from Rajapalaiyam in the Madurai Granulite Block, southern India. In this unique textural association, hitherto unreported in previous studies, sapphirine shows four occurrences: (1) as anhedral mineral between spinel and quartz (Spr-1), (2) subhedral to euhedral needles mantled by quartz (Spr-2), (3) subhedral to anhedral mineral in orthopyroxene, and (4) isolated inclusion with quartz (Spr-4). Spr-1, Spr-2, and Spr-4 show direct grain contact with quartz, providing evidence for ultrahigh-temperature (UHT) metamorphism at temperatures exceeding 1000 °C. Associated orthopyroxene shows high Mg/(Fe + Mg) ratio ( 0.75) and Al₂O₃ content (up to 9.6 wt.%), also suggesting T > 1050 °C and P > 10 kbar during peak metamorphism.

Coarse spinel (Spl-1) with irregular grain morphology and adjacent quartz grains are separated by thin films of Spr-1 and K-feldspar, suggesting that Spl-1 and quartz were in equilibrium before the stability of Spr-1 + quartz. This texture implies that the P–T conditions of the rock shifted from the stability field of spinel + quartz to sapphirine + quartz. Petrogenetic grid considerations based on available data from the FMAS system favour exhumation along a counterclockwise P–T trajectory. The irregular shape of the inclusion and chemistry of the inclusion minerals are markedly different from the matrix phases suggesting the possibility that the inclusion minerals could have equilibrated from cordierite-bearing silicate-melt pockets during the garnet growth at extreme UHT conditions.     

An Anisotropic Modelling for the Determination of Regional Attenuation Coefficients

The anisotropic attenuation of macroseismic intensity for a seismogenetic zone is dealt with using a new modelling of intensity distribution. The analysis, carried out starting from the intensity maps of the earthquakes of different seismogenetic zones of Central and Southern Italy, allows the determination of the attenuation coefficients for each seismogenetic zone by an anisotropic attenuation law. The obtained results show the reliability of the proposed modelling within seismic hazard evaluation studies.     

Late Quaternary basin evolution of the Gulf of Corinth: Sequence stratigraphy, sedimentation, fault–slip and subsidence rates

The Gulf of Corinth is a graben, which has undergone extension during the Late Quaternary. The subsidence rate is rapid in the currently marine part whereas uplift now affects a large part of the initially subsiding area in the North Peloponnese. In this paper, we document the rates of subsidence/uplift and extension based on new subsurface data, including seismic data and long piston coring in the deepest part of the Gulf. Continuous seismic profiling data (air gun) have shown that four (at least) major oblique prograding sequences can be traced below the northern margin of the central Gulf of Corinth. These sequences have been developed successively during low sea level stands, suggesting continuous and gradual subsidence of the northern margin by 300 m during the Late Quaternary (last 250 ka). Subsidence rates of 0.7–1.0 m kyr^− 1 were calculated from the relative depth of successive topset to foreset transitions. The differential total vertical displacement between the northern and the southern margins of the Corinth graben is estimated at about 2.0–2.3 m kyr^− 1.

Sequence stratigraphic interpretation of seismic profiles from the basin suggests that the upper sediments (0.6 s twtt thick) in the depocenter were accumulated during the last 250 ka at a mean rate of 2.2–2.4 m kyr^− 1. Long piston coring in the central Gulf of Corinth basin enabled the recovery of lacustrine sediments, buried beneath 12–13.5 m of Holocene marine sediments. The lacustrine sequence consists of varve-like muddy layers interbedded with silty and fine sand turbidites. AMS dating determined the age of the marine–lacustrine interface (reflector Z) at about 13 ka BP. Maximum sedimentation rates of 2.4–2.9 m kyr^− 1 were calculated for the Holocene marine and the last glacial, lacustrine sequences, thus verifying the respective rates obtained by the sequence stratigraphic interpretation. Recent accumulation rates obtained by the ²¹⁰Pb-radiometric method on short sediment box cores coincide with the above sedimentation rates. Vertical fault slip rates were measured by using fault offsets of correlated reflector Z. The maximum subsidence rate of the depocenter (3.6 m kyr^− 1) exceeds the maximum sedimentation rate by 1.8 m kyr^− 1, which, consequently, corresponds to the rate of deepening of the basin's floor. The above rates indicate that the 2.2 km maximum sediment thickness as well as the 870 m maximum depth of the basin may have formed during the last 1 Ma, assuming uniform mean sedimentation rate throughout the evolution of the basin.     

Age and emplacement of late-Variscan granites of the western Bohemian Massif with main focus on the Hauzenberg granitoids (European Variscides, Germany)

The Variscan Hauzenberg pluton consists of granite and granodiorite that intruded late- to postkinematically into HT-metamorphic rocks of the Moldanubian unit at the southwestern margin of the Bohemian Massif (Passauer Wald). U–Pb dating of zircon single-grains and monazite fractions, separated from medium- to coarse-grained biotite-muscovite granite (Hauzenberg granite II), yielded concordant ages of 320 ± 3 and 329 ± 7 Ma, interpreted as emplacement age. Zircons extracted from the younger Hauzenberg granodiorite yielded a ²⁰⁷Pb–²⁰⁶Pb mean age of 318.6 ± 4.1 Ma. The Hauzenberg granite I has not been dated. The pressure during solidification of the Hauzenberg granite II was estimated at 4.6 ± 0.6 kbar using phengite barometry on magmatic muscovite, corresponding to an emplacement depth of 16-18 km. The new data are compatible with pre-existing cooling ages of biotite and muscovite which indicate the Hauzenberg pluton to have cooled below T = 250–400 °C in Upper Carboniferous times. A compilation of age data from magmatic and metamorphic rocks of the western margin of the Bohemian Massif suggests a west- to northwestward shift of magmatism and HT/LP metamorphism with time. Both processes started at > 325 Ma within the South Bohemian Pluton and magmatism ceased at ca. 310 Ma in the Bavarian Oberpfalz. The slight different timing of HT metamorphism in northern Austria and the Bavarian Forest is interpreted as being the result of partial delamination of mantle lithosphere or removal of the thermal boundary layer.     

Coseismic reactivation of the Samambaia fault, Brazil

Northeastern Brazil is, within the present knowledge of historical and instrumental seismicity, one the most seismic active areas in intraplate South America. Seismic activity in the region has occurred mainly around the Potiguar basin. This seismicity includes earthquake swarms characterized by instrumentally-recorded events ≤ 5.2 m_b and paleoseismic events ≥ 7.0. Our study concentrates in the João Câmara (JC) epicentral area, where an earthquake swarm composed of more than 40,000 aftershocks occurred mainly from 1986 to 1990 along the Samambaia fault; 14 of which had m_b > 4.0 and two of which had 5.1 and 5.0 m_b. We describe and compare this aftershock sequence with the present-day stress field and the tectonic fabric in an attempt to understand fault geometry and local control of seismogenic faulting. Earthquake data indicate that seismicity decreased steadily from 1986 to 1998. We selected 2,746 epicenters, which provided a high-quality and precise dataset. It indicates that the fault trends 37° azimuth, dips 76°–80° to NW, and forms an alignment  27 km long that cuts across the NNE–SSW-trending ductile Precambrian fabric. The depth of these events ranged from  1 km to  9 km. The fault forms an echelon array of three main left-bend segments: one in the northern and two in the southern part of the fault. A low-seismicity zone, which marks a contractional bend, occurs between the northern and southern segments. Focal mechanisms indicate that the area is under an E–W-oriented compression, which led to strike–slip shear along the Samambaia fault with a small normal component. The fault is at 53° to the maximum compression and is severely misoriented for reactivation under the present-day stress field. The seismicity, however, spatially coincides with a brittle fabric composed of quartz veins and silicified-fault zones. We conclude that the Samambaia fault is a discontinuous and reactivated structure marked at the surface by a well-defined brittle fabric, which is associated with silica-rich fluids.     

PGE and Os-isotopic variations in lavas from Kohala Volcano, Hawaii: Constraints on PGE behavior and melt/crust interaction

We have examined Re, Platinum-Group Element (PGE) and Os-isotope variations in suites of variably fractionated lavas from Kohala Volcano, Hawaii, in order to evaluate the effects of melt/crust interaction on the mantle isotopic signature of these lavas. This study reveals that the behavior of Os and other PGEs changes during magma differentiation. The concentrations of all PGEs strongly decrease with increasing fractionation for melts with MgO < 8 wt.%. Fractionation trends indicate significantly higher bulk partition coefficients for PGEs in lavas with less than 8 wt.% MgO (D_PGE = 35–60) when compared to values for more primitive lavas with MgO > 8 wt.% (D_PGE ≤ 6). This sudden change in PGE behavior most likely reflects the onset of sulfur saturation and sulfide fractionation in Hawaiian magmas at about 8 wt.% MgO.

The Os-rich primitive lavas (≥ 8 wt.% MgO, > 0.1 ppb Os) display a narrow range of ¹⁸⁷Os/¹⁸⁸Os values (0.130–0.133), which are similar to values in high-MgO lavas from Mauna Kea and Haleakala Volcanoes and likely represent the mantle signature of Kohala lavas. However, Os-isotopic ratios become more radiogenic with decreasing MgO and Os content in evolved lavas, ranging from 0.130 to 0.196 in the shield-stage Pololu basalts and from 0.131 to 0.223 in the post-shield Hawi lavas. This reflects assimilation of local oceanic crust material during fractional crystallization of the magma at shallow level (AFC processes). AFC modeling suggests that assimilation of up to 10% upper oceanic crust could produce the most radiogenic Os-isotope ratios recorded in the Pololu lavas. This amount of upper crust assimilation has a negligible effect on the Sr and Nd-isotopic compositions of Kohala lavas. Thus, these isotopic compositions likely represent the composition of the mantle source of Kohala lavas.     

Mineralogical and chemical transformation of Oldoinyo Lengai natrocarbonatites, Tanzania

The minerals of Oldoinyo Lengai natrocarbonatite lavas are unstable under atmospheric conditions. Subsolidus mineral assemblages in natrocarbonatites were studied in 105 samples from contemporary eruptions ranging from present day to about 100 years old. The subsolidus minerals in natrocarbonatites were formed (i) along cracks on the lava surface from hot gases escaping during cooling, (ii) as atmospheric alteration by solution of water-soluble minerals, in particular halides and gregoryite, and by hydration of nyerereite under the influence of meteoric water and (iii) by reaction with fumarole gases. After solidification, the lavas were cut by a network of thin cracks, the edges of which are covered by polymineralic encrustations. Samples collected 2–24 h after eruption contain nahcolite, trona, sylvite, and halite with accessory kalicinite and villiaumite. Atmospheric humidity results immediately (≥ 2 h after eruption) in alteration of black lavas that is marked by the appearance of white powdery thermonatrite with nahcolite on the lava surface. Subsequent reaction (weeks, months, years) of natrocarbonatite with meteoric water and the atmosphere results in the formation of pirssonite, gaylussite, shortite, trona, thermonatrite, nahcolite and calcite. Generally, the first important step is the formation of pirssonite and the end-members are calcite carbonate rocks or loose aggregates. Fumarolic activity is common for the active northern crater of the volcano. Reaction of hot (54–141 °C) fumarolic gases with natrocarbonatite leads to the formation of sulphur, gypsum, calcite, anhydrite, monohydrocalcite, barite and celestine. Changes in mineralogy of the natrocarbonatite lead to substantial chemical transformation. The most obvious chemical changes in this process are the loss of Na, K, Cl and S, combined with an increase in H₂O, Ca, Sr, Ba, F and Mn. The oxygen and carbon isotopic composition of altered natrocarbonatites shows a significant shift from the primary “Lengai Box” to high values of δ¹⁸O and δ¹³C. Calcite exhibits δ¹³C values between − 2‰ and − 4‰ PDB and δ¹⁸O values of + 23‰ to + 26‰ SMOW. The observed assemblages of secondary minerals formed by reaction with atmosphere and meteoric water, the changes in chemical composition of the natrocarbonatite and field observations suggest that alteration of natrocarbonatite is an open-system low-temperature process. It takes place at temperatures between 8 and 43 °C with the addition of H₂O to the system and the removal of Na, K, Cl and S from the carbonatites. Low-temperature thermodynamic models developed for alkali carbonate systems can be used for the interpretation of Oldoinyo Lengai subsolidus mineralization.