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.     

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.     

Characterizing hydrological processes within kettle holes using stable water isotopes in the Uckermark of northern Brandenburg,Germany

Understanding the hydrologic connectivity between kettle holes and shallow groundwater, particularly in reaction to the highly variable local meteorological conditions, is of paramount importance for tracing water in a hydro(geo)logically complex landscape and thus for integrated water resource management. This article is aimed at identifying the dominant hydrological processes affecting the kettle holes' water balance and their interactions with the shallow groundwater domain in the Uckermark region, located in the north-east of Germany. For this reason, based on the stable isotopes of oxygen (δ¹⁸O ) and hydrogen (δ²H ), an isotopic mass balance model was employed to compute the evaporative loss of water from the kettle holes from February to August 2017. Results demonstrated that shallow groundwater inflow may play the pivotal role in the processes taking part in the hydrology of the kettle holes in the Uckermark region. Based on the calculated evaporation/inflow (E/I) ratios, most of the kettle holes (86.7%) were ascertained to have a partially open, flow-through-dominated system. Moreover, we identified an inverse correlation between E/I ratios and the altitudes of the kettle holes. The same holds for electrical conductivity (EC) and the altitudes of the kettle holes. In accordance with the findings obtained from this study, a conceptual model explaining the interaction between the shallow groundwater and the kettle holes of Uckermark was developed. The model exhibited that across the highest altitudes, the recharge kettle holes are dominant, where a lower ratio of E/I and a lower EC was detected. By contrast, the lowest topographical depressions represent the discharge kettle holes, where a higher ratio of E/I and EC could be identified. The kettle holes existing in between were categorized as flow-through kettle holes through which the recharge takes place from one side and discharge from the other side.     

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.     

Easy method of interpreting dyke anomalies

Summary An entirely new procedure is proposed in this paper for interpreting anomalies of dykes. This reduces the ambiguity in magnetic interpretation. Measurements of the first vertical derivative simplifies the problem of interpretation and also supplies additional information concerning the parameters of a dyke. The first vertical derivative profile cuts the distance axis at two points whose separation is related to the depth and direction of magnetisation of a dyke.