The last 40 ka tephrostratigraphic record of Lake Ohrid,Albania and Macedonia: a very distal archive for ash dispersal from Italian volcanoes

A 1075 cm long core (Lz1120) was recovered in the south-eastern part of the Lake Ohrid (Republics of Macedonia and Albania) and sampled for identification of tephra layers. Magnetic susceptibility investigations show rather high magnetic values throughout the core, with peaks unrelated to the occurrence of tephra layers but instead to the relative abundance of detrital magnetic minerals in the sediment. Naked-eye inspection of the core allowed us to identify of two tephra layers, at 896–897 cm and 1070–1075 cm. Laboratory inspection of the grain-size fraction > 125 μm allowed for the identification of a third cryptotephra at 310–315 cm. Major element analyses on glass shards of the tephra layers at 896–897 cm and 1070–1075 cm show a trachytic composition, and indicate a correlation with the regionally dispersed Y-3 and Y-5 tephra layers, dated at ca 30 and 39 cal ka BP. The cryptotephra at 310–315 cm has a mugearitic–benmoreitic composition, and was correlated with the FL eruption of Mt. Etna, dated at 3370 ± 70 cal yr BP. These ages are in agreement with five ¹⁴C AMS measurements carried out on plant remains and macrofossils from the lake sediments at different depths along the core.     

Geothermal energy release at the Solfatara of Pozzuoli (Phlegraean Fields): Phreatic and phreatomagmatic explosion risk implications

The H₂O, CO₂ and H₂S outputs at the Solfatara of Pozzuoli have been measured and a map of the exhaling areas has also been made. The energy released at the surface by the fluids has been estimated to be 10¹⁹ ergs/day.The presence of aquifers at Phlegraean Fields increases the phreatic and phreatomagmatic explosion risk.Our results suggest that even if an uprising magma may interact with water at depth, an explosion could occur only at the shallow levels of a few hundred meters. Since the transfer of energy toward the surface is favoured by the presence of fractures, a detailed analysis of the deep fracture network would help to evaluate the risk levels of the various areas of Phlegraean Fields.     

COMPUTER MODELING OF VERTICAL SEISMIC PROFILING*

There are numerous modeling techniques commonly employed for the computer simulation of seismic wave propagation. The capabilities of these techniques vary according to the theoretical foundations and subsequent approximations upon which the algorithms are based. This paper constitutes a comparative review of seven modeling methods: geometric ray theory, asymptotic ray theory, generalized ray theory, Kirchhoff wave theory, Fourier synthesis, finite differences, and finite elements. These methods can be categorized as ray or wave, acoustic or elastic, and can be contrasted according to their relative abilities to simulate such behavior as wave interference effects, diffractions, and mode conversions. As is implied by their names, geometric ray theory and asymptotic ray theory are both ray methods. The other five methods provide wave theory simulations. Geometric ray theory and Kirchhoff wave theory are normally implemented in acoustic form, while the other methods are readily adapted for computing elastic theory solutions. Generalized ray theory and Fourier synthesis are more limited in the complexity of geological model they can accommodate than are the other techniques. The methods which typically demand the greatest computer resources are the finite-difference and the finite-element techniques. All methods can incorporate at least some multiple events. Diffractions, however, are only inherent in the solutions computed by Kirchhoff wave theory, finite differences and finite elements. Attenuation is readily incorporated in both the Fourier synthesis and the finite-element methods. As an example of the application of seismic modeling, a geological model representative of a typical petroleum exploration target is used to compare vertical seismic profiles calculated by different modeling methods.     

High-resolution garnet chronometry and the rates of metamorphic processes

Garnets in an amphibolite-facies metasediment from Sulitjelma, North Norway yield precise and concordant SmNd, UPb and RbSr ages that relate directly to the pressure (P) and temperature (T) conditions of mineral growth. Differential mineral reaction between graphitic and non-graphitic layers within this sample preserves a record of theP-T and time (t) history experienced during Barrovian regional metamorphism. Garnets in graphitic layers grew during prograde metamorphism at462 ± 16°C and5.2 ± 0.5 kbar under conditions of lowaH₂O, and yield indistinguishable¹⁴⁷Sm¹⁴³Nd and²³⁸U²⁰⁶Pb ages of434.1 ± 1.2 Ma and433.9 ± 1.0 Ma, respectively. In contrast, garnet growth in adjacent graphite-free layers did not occur untilP-T conditions of540 ± 18°C and8.0 ± 1.0 kbar were attained, with continued growth in response to minor heating and decompression with final matrix equilibration at544 ± 16°C and7.0 ± 1.0 kbar. The inclusion-free garnet rims in this assemblage record indistinguishable¹⁴⁷Sm¹⁴³Nd and²³⁸U²⁰⁶Pb ages of424.6 ± 1.2 Ma and423.4± 1.7 Ma, respectively. These results provide precise estimates for average heating and burial rates during prograde metamorphism of 8.6_−4.4^+7.5°C Ma⁻¹ and 0.8_−0.5^+0.9 km Ma⁻¹, respectively. Rb and Sr exchange between coexisting silicates in the graphite-free assemblage continued for some 37 Ma after the “peak” of metamorphism, and require an average cooling rate of about 4.0°C Ma⁻¹ during uplift. These results illustrate a clear relationship between reaction history and the timing of mineral growth and provide definitive constraints on the rates of thermal and tectonic processes accompanying regional metamorphism.     

Impact of user parameters in inversion velocity analysis

Migration velocity analysis is a method devoted to the evaluation of both reflectivity and background velocity models, associated with the high and low wavenumber components of the model, respectively. Inversion velocity analysis is one of its improved versions, leading to more stable background velocity updates. Still, the impact of the user parameters should be understood for an optimal update of the background velocity. We show that a sign reversal of the background velocity gradient could occur when the selected surface offset range or the space lag range is too small. We derive the theoretical limits and check their consistency through simulations in a simple model with a single interface. These guidelines determine the necessary ranges of surface offsets and space lags for a proper update of the background velocity model. We discuss their applicability on the Marmousi model. Artefacts in the retrieved background velocity model are observed when the guidelines are not satisfied.     

Statistical assessment of seismicity patterns in Italy: Are they precursors of subsequent events?

The variations of seismicity rate in Central Apenninesprior to the sequence started in September, 1997 (at00:33 UTC, M _L5.6) has been analysedby statistical methods, with the purpose of pointingup eventual periods of quiescence. The analysis wascarried out on the instrumental catalogue of theIstituto Nazionale di Geofisica (ING), covering theperiod from January 1975 to March 1998. In apreliminary phase, the catalogue was declustered usingthe Reasenberg algorithm. After that, eventualmagnitude shifts due to variations in the modalitiesof observation have been individuated and corrected.The subsequent analysis, carried out making use of theZmap software package, has put in evidence thatthe sequence of September 1997 was preceded bya 2.5 year period characterised by absence of eventsof magnitude larger than 3.2, in an area approximately20 × 40 km wide, including the epicentre of themain shock. The statistical methodology shows thatonly 1/10³ of the space-time volumes analysed inthis study, exhibited quiescence of the same level.The study of seismicity rate change correlated toprevious main shocks in a larger area of CentralApennines shows that none of them were preceded by aseismic quiescence, specially close to the epicentreof the main shock, and lasting until the time ofoccurrence of the main shock as in the 1997 case.Actually, we found other patterns of precursoryquiescence with different time or space distribution.We conclude that precursory quiescence is a realfeature of Central Apennines seismicity, but it isdifficult to define a simple hypothesis, which appliesto the generality of cases and can be tested beforeimplementation in a system of earthquake riskmitigation.     

Observations of vertical ground accelerations exceeding gravity during the 1997 Umbria-Marche (central Italy) earthquakes

We found extensive evidence that the vertical ground accelerations produced during the largest shock (M = 6.0) of the 1997 Umbria-Marche earthquake sequence exceeded 1g in two areas close to the heavily-damaged villages of Annifo and Colle Croce. This evidence comes from the striking observation of thousands of freshly fractured and broken rocks and stones in these areas. Some of the broken stones lie isolated on soft detritic soil while others had been previously piled up, probably a long time agoto clear the fields for farming. The freshness of the cuts and fractures and the consistency of the observations for thousands of rocks and stones in these areas indicate that these rocks were thrown upwards during the earthquake, with breakage occurring at the time of impact. Ground motion calculations consistent with the static deformation inferred from GPS and interferometry data, show that the broken stones and rocks are found in the zone where the strongest shaking took place during the earthquake and that most of the shaking there was vertical.