Improving automatic location procedure by waveform similarity analysis: An application in the South Western Alps (Italy)

The accuracy of automatic procedures for locating earthquakes is influenced by several factors such as errors in picking seismic phases, network geometry, modeling errors and velocity model uncertainties. The main purpose of this work is to improve the performances of the automatic procedure employed for the “quasi-real-time” location of seismic events in North Western Italy by developing a procedure based on a waveform similarity analysis and by using only one seismic station.To detect “earthquake families” a cross-correlation technique was applied to a data set of seismic waveforms recorded in the period 1985-2002, in a small test area (1600 km²) located in the South Western Alps (Italy). Normalized cross-correlation matrices were calculated using about 2700 seismic events, selected on the basis of the signal to noise ratio, manually picked and located by using the Hypoellipse code. The waveform similarity analysis, based on the bridging technique, allowed grouping about 65% of the selected events into 80 earthquake families (multiplets) located inside the area considered. For each earthquake family a master event is selected, manually re-picked and re-located by using Hypoellipse code. Having chosen a reference station (STV) on the basis of the completeness of the available data set, an automatic procedure has been developed with the aim of cross-correlating new seismic recordings (automatically picked) to the waveforms of the events belonging to the detected families. If the new event is proved to belong to a family (on the basis of the cross-correlation values), its hypocenter co-ordinates are defined by the location of the master event of the associated family. The performance of the proposed procedure is tested and demonstrated using a data set of 104 selected earthquakes recorded in the period January 2003-June 2004 and located in the test area. The automatic procedure is able to locate, associating events with the multiplets detected by the waveform similarity analysis, about 50% of the test events, almost independently of the accuracy of the automatic phase picker and without the biasing of the network geometry and of the velocity model uncertainties.     

Detection of earthquake clusters on the basis of waveform similarity: An application in the monferrato region (Piedmont, Italy)

In August 2000 and July 2001 two seismic sequences, characterized by mainshocks with Ml (local magnitude) respectively 5.1 and 4.8, occurred in the Monferrato region (Italy). The regional seismic network of North-western Italy (RSNI) recorded more than 250 foreshocks and aftershocks. The routine locations, obtained from the Hypoellipse code, show a seismic activity concentrated in a circular area, of a radius of about 15 km, located near Acqui Terme, and randomly distributed in depth. Location errors, due to an uneven azimuthal station distribution of the regional seismic network, prevent recognition of the geometry of the active zone. Waveform analysis revealed the presence of several multiplets. In order to discriminate and successively relocate them, an automatic cross-correlation procedure was applied. Normalized cross-correlation matrixes, for the RSNI stations which recorded almost 90% of considered events, on the basis of a signal to noise ratio analysis, were computed using only S wave time windows. The use of a relocation procedure, based on the double-difference method which incorporates ordinary absolute travel-time measurements and/or cross-correlation differential travel-times, allowed us to overcome the constraints of the uneven distribution of stations giving a quite different frame of seismicity. The improved locations showed that the seismic activity is mainly arranged along a NE-SW oriented volume, at a depth range of 8–20 km, involving the basement crystalline units. This orientation is confirmed by the analysis of the focal mechanisms: most focal solutions show a strike slip component with one of the nodal planes consistent with the main orientation of the seismic events.     

Compositional classification of “kimberlitic” and “non-kimberlitic” ilmenite

Ilmenite is one of the common kimberlitic indicator minerals recovered during diamond exploration, and its distinction from non-kimberlitic rock types is important. This is particularly true for regions where these minerals are present in relatively low abundance, and they are the dominant kimberlitic indicator mineral recovered. Difficulty in visually differentiating kimberlitic from non-kimberlitic ilmenite in exploration concentrates is also an issue, and distinguishing kimberlitic ilmenite from those derive from other similar rocks, such as ultramafic lamprophyres, is practically impossible. Ilmenite is also the indicator mineral whose compositional variety has the most potential to resolve provenance issues related to mineral dispersions with contributions from multiple kimberlite sources.

Various published data sets from selected kimberlitic (including kimberlites, lamproites, and various ultramafic lamprophyres) and non-kimberlitic rock types have been compiled and evaluated in terms of their major element compositions. Compositional fields and bounding reference lines for ilmenites derived from kimberlites (sensu stricto), ultramafic lamprophyres, and other non-kimberlitic rock types have been defined primarily on MgO–TiO₂ graphs as well as MgO–Cr₂O₃ relationships.     

Covariance-Based Variable Selection for Compositional Data

Omitting variables in compositional data analysis may lead to a substantial change in results from that of multivariate statistical analysis. In particular, this is the case for principal component analysis and the compositional biplot, where both the interpretation of loadings and scores of the remaining subcomposition are affected. A stepwise procedure is introduced that allows for a reduction of the original composition to a subcomposition by avoiding a substantial change of the information, like those carried by the compositional biplot. The subcomposition is easier to handle and interpret. Numerical results give evidence of the usefulness of the procedure.     

The combined inelastic neutron scattering and solid state DFT study of hydrogen atoms dynamics in a highly ordered kaolinite

The dynamics of the hydrogen atoms in the highly ordered kaolinite was studied by vibrational spectroscopy based on inelastic neutron scattering method with the focus on the spectral region of 100–1,250 cm⁻¹. The experimental spectrum was interpreted by means of the solid state density functional theory calculations covering both normal mode analysis and molecular dynamics going beyond the harmonic approximation. The Al–O–H bending modes were found to be spread over the large interval of 100–1,100 cm⁻¹, with the dominant contributions located between 800 and 1,100 cm^-1. The shapes of the individual hydrogen spectra depend on the strengths of the individual interlayer O–H···O hydrogen bonds involving the inner surface hydroxyl groups. The modes assigned to the in-plane movements of the respective hydrogen atoms are well-defined and always appear on the top of the intervals of energy transfer. In contrast, the modes generated by the out-of-plane movements are spread over large intervals of energies spanning down to the region of external (lattice) modes.     

Intercomparison of Boron Isotope and Concentration Measurements. Part II: Evaluation of Results

The Istituto di Geoscienze e Georisorse (IGG), on behalf and with the support of the International Atomic Energy Agency (IAEA), prepared eight geological materials (three natural waters and five rocks and minerals), intended for a blind interlaboratory comparison of measurements of boron isotopic composition and concentration. The materials were distributed to twenty seven laboratories - virtually all those performing geochemical boron isotope analyses in the world -which agreed to participate in the intercomparison exercise. Only fifteen laboratories, however, ultimately submitted the isotopic and/or concentration results they obtained on the intercomparison materials. The results demonstrate that interlaboratory reproducibility is not well reflected by the precision values reported by the individual laboratories and this observation holds true for both boron concentration and isotopic composition. The reasons for the discrepancies include fractionations due to the chemical matrix of materials, relative shift of the zero position on the δ¹¹B scale and a lack of well characterized materials for calibrating absolute boron content measurements. The intercomparison materials are now available at the IAEA (solid materials) and IGG (waters) for future distribution.     

Pre-processing of a sample of multi-scene and multi-date Landsat imagery used to monitor forest cover changes over the tropics

In support to the Remote Sensing Survey of the global Forest Resource Assessment 2010, the TREES-3 project has processed more than 12,000 Landsat TM and ETM+ data subsets systematically distributed over the tropics. The project aims at deriving area estimates of tropical forest cover change for the periods 1990-2000-2005. The paper presents the pre-processing steps applied in an operational and robust manner to this large amount of multi-date and multi-scene imagery: conversion to top-of-atmosphere reflectance, cloud and cloud shadow detection, haze correction and image radiometric normalization. The results show that the haze correction algorithm has improved the visual appearance of the image and significantly corrected the digital numbers for Landsat visible bands, especially the red band. The impact of the normalization procedures (forest normalization and relative normalization) was assessed on 210 image pairs: in all cases the correlation between the spectral values of the same land cover in both images was improved. The developed automatic pre-processing chain provided a consistent multi-temporal data set across the tropics that will constitute the basis for an automatic object-based supervised classification.