Reliability of the automatic procedures for locating earthquakes in southwestern Alps and northern Apennines (Italy)

Reliable automatic procedure for locating earthquake in quasi-real time is strongly needed for seismic warning system, earthquake preparedness, and producing shaking maps. The reliability of an automatic location algorithm is influenced by several factors such as errors in picking seismic phases, network geometry, and velocity model uncertainties. The main purpose of this work is to investigate the performances of different automatic procedures to choose the most suitable one to be applied for the quasi-real-time earthquake locations in northwestern Italy. The reliability of two automatic-picking algorithms (one based on the Characteristic Function (CF) analysis, CF picker, and the other one based on the Akaike’s information criterion (AIC), AIC picker) and two location methods (“Hypoellipse” and “NonLinLoc” codes) is analysed by comparing the automatically determined hypocentral coordinates with reference ones. Reference locations are computed by the “Hypoellipse” code considering manually revised data and tested using quarry blasts. The comparison is made on a dataset composed by 575 seismic events for the period 2000–2007 as recorded by the Regional Seismic network of Northwestern Italy. For P phases, similar results, in terms of both amount of detected picks and magnitude of travel time differences with respect to manual picks, are obtained applying the AIC and the CF picker; on the contrary, for S phases, the AIC picker seems to provide a significant greater number of readings than the CF picker. Furthermore, the “NonLinLoc” software (applied to a 3D velocity model) is proved to be more reliable than the “Hypoellipse” code (applied to layered 1D velocity models), leading to more reliable automatic locations also when outliers (wrong picks) are present.     

A model for the viscosity of rhyolite as a function of H2O-content and pressure: A calibration based on centrifuge piston cylinder experiments

The Newtonian viscosity of synthetic rhyolitic liquids with 0.15-5.24 wt% dissolved water was determined in the interval between 580 and 1640 °C and pressures of 1 atm and 5-25 kbar. Measurements were performed by combining static and accelerated (up to 1000g) falling sphere experiments on water-bearing samples, with high temperature concentric cylinder experiments on 0.15 wt% H₂O melts. These methods allowed viscosity determinations between 10² and 10⁷ Pa s, and cover the complete range of naturally occurring magmatic temperatures, pressures, and H₂O-contents for rhyolites.Our viscosity data, combined with those from previous studies, were modeled by an expression based on the empirical Vogel-Fulcher-Tammann equation, which describes viscosities and derivative properties (glass transition temperature T_g, fragility m, and activation volume of viscous flow Va) of silicic liquids as a function of P-T-X_(H2O). The fitted expressions do not account for composition-dependent parameters other than X_(H2O) and reproduce the entire viscosity database for silicic liquids to within 3.0% average relative error on log η (i.e. std. error of estimate of 0.26 log units).The results yield the expected strong decrease of viscosity with temperature and water content, but show variable pressure dependencies. Viscosity results to be strongly affected by pressure at low pressures; an effect amplified at low temperatures and water contents. Fragility, as a measure for the deviation from Arrhenian behavior, decreases with H₂O-content but is insensitive to pressure. Activation volumes are always largely negative (e.g., less than −10 cm³/mol) and increase strongly with H₂O-content. Variations in melt structure that may account for the observed property variations are discussed.     

Improved 2-D attenuation analysis for Northern Italy using a merged dataset from selected regional seismic networks

A merged, high-quality waveform dataset from different seismic networks has been used to improve our understanding of lateral seismic attenuation for Northern Italy. In a previous study on the same region, Morasca et al. (Bull Seismol Soc Am 98:1936–1946, 2008) were able to resolve only a small area due to limited data coverage. For this reason, the interpretation of the attenuation anomalies was difficult given the complexity of the region and the poor resolution of the available data. In order to better understand the lateral changes in the crustal structure and thickness of this region, we selected 770 earthquakes recorded by 54 stations for a total of almost 16,000 waveforms derived from seismic networks operating totally or partially in Northern Italy. Direct S-wave and coda attenuation images were obtained using an amplitude ratio technique that eliminates source terms from the formulation. Both direct and early-coda amplitudes are used as input for the inversions, and the results are compared. Results were obtained for various frequency bands ranging between 0.3 and 25.0 Hz and in all cases show significant improvement with respect to the previous study since the resolved area has been extended and more crossing paths have been used to image smaller scale anomalies. Quality-factor estimates are consistent with the regional tectonic structure exhibiting a general trend of low attenuation under the Po Plain basin and higher values for the Western Alps and Northern Apennines. The interpretation of the results for the Eastern Alps is not simple, possibly because our resolution for this area is still not adequate to resolve small-scale structures.     

LightSquared effects on estimated C/N 0, pseudoranges and positions

We present experimental results showing the impact of the proposed LightSquared (LS) Long-term Evolution (LTE) signals on both GPS and Galileo civil modulations in the L1/E1 band. The experiments were conducted in radiated mode in a large anechoic chamber. Three Galileo enabled receivers were chosen for the tests, and a state of the art GNSS signal generator was used to simulate both GPS and Galileo signals. The LTE signals were generated by an Agilent Programmable Signal Generator with a license to generate the signals according to the 3GPP LTE FDD standard. The interference impact was measured in terms of a Carrier-to-Noise power spectral density ratio (C/N ₀) degradation, in accordance with the methodology which the LS/GPS Technical Working Group (TWG) established by mandate of the FCC. A model for determining the impact of the LS signal on the considered GNSS signals is provided and is validated against experimental data. It is shown that the Galileo E1 Open Service (OS) signal is marginally more susceptible to this form of interference than the GPS L1 C/A signal due to its greater proximity to the lower edge of the L1 band. The impact of LS interference was further analyzed in terms of pseudorange and position errors. Despite its relevance for most GNSS users, this aspect was not considered by the TWG. Measurement and position domain analysis along with the study of the LS impact on the Galileo OS signals are the major contributions. The analysis confirms the results obtained by the TWG and shows that the receiver front-end plays a major role in protecting GNSS signals against RF interference. While it appears that, for now, the LS network will not be deployed, the approach taken and the results obtained herein can be readily adapted for any future terrestrial mobile network that may take the place of LS.     

Contrasting styles of Mount Vesuvius activity in the period between the Avellino and Pompeii Plinian eruptions,and some implications for assessment of future hazards

Intense explosive activity occurred repeatedly at Vesuvius during the nearly 1,600-year period between the two Plinian eruptions of Avellino (3.5 ka) and Pompeii (79 A.D.). By correlating stratigraphic sections from more than 40 sites around the volcano, we identify the deposits of six main eruptions (AP1-AP6) and of some minor intervening events. Several deposits can be traced up to 20 km from the vent. Their stratigraphic and dispersal features suggest the prevalence of two main contrasting eruptive styles, each involving a complex relationship between magmatic and phreatomagmatic phases. The two main eruption styles are (1) sub-Plinian to phreato-Plinian events (AP1 and AP2 members), where deposits consist of pumice and scoria fall layers alternating with fine-grained, vesiculated, accretionary lapilli-bearing ashes; and (2) mixed, violent Strombolian to Vulcanian events (AP3-AP6 members), which deposited a complex sequence of fallout, massive to thinly stratified, scoria-bearing lapilli layers and fine ash beds. Morphology and density variations of the juvenile fragments confirm the important role played by magma-water interaction in the eruptive dynamics. The mean composition of the ejected material changes with time, and shows a strong correlation with vent position and eruption style. The ranges of intensity and magnitude of these events, derived by estimations of peak column height and volume of the ejecta, are significantly smaller than the values for the better known Plinian and sub-Plinian eruptions of Vesuvius, enlarging the spectrum of the possible eruptive scenarios at Vesuvius, useful in the assessment of its potential hazard.     

Sequential solution to Kepler’s equation

Seven sequential starter values for solving Kepler’s equation are proposed for fast orbit propagation. The proposed methods have constant complexity (not iterative), do not require pre-computed data, and can be implemented in just a few lines of code. The resulting sequential orbit propagation techniques can be done at different levels of accuracy and speed, depending essentially on the value of orbit eccentricity. Accuracy and algorithmic complexity are evaluated for all the proposed approaches and compared with several existing single-point techniques to solve Kepler’s equation. The new methods obtain improved accuracy at lower computational cost as compared to the best existing methods.     

S-wave splitting intensity analysis and inversion

Analysing S-wave splitting has become a routine step in processing multicomponent data. Typically, this analysis leads to determining the principal directions of a transversely isotropic medium with a horizontal symmetry axis, which is assumed to be responsible for azimuthal anisotropy, and to the time delays between the fast and slow S-waves. These parameters are commonly estimated layer-by-layer from the top. Errors in layer stripping occurring in shallow layers might propagate to deeper layers. We propose a method for S-wave splitting analysis and compensation that consists of inverting interval values of splitting intensity to obtain a model of anisotropic parameters that vary with time and/or depth. Splitting intensity is a robust attribute with respect to structural variations and is commutative, which means that it can be summed along a ray (or throughout a sensitivity kernel volume) and can be linearly related to anisotropic perturbations at depth. Therefore, it is possible to estimate anisotropic properties within a geological formation (e.g. the reservoir) by analysing the differences of splitting intensity measured at the top and at the bottom of the layer. This allows us to avoid layer stripping, in particular, for shallow layers where anisotropic parameters are difficult to estimate due to poor coverage, and it makes S-wave splitting analysis simpler to apply. We demonstrate this method on synthetic and real data. Because the splitting intensity attribute shows usefulness in S-wave splitting analysis in transversely isotropic media, we extend the splitting intensity theory to lower symmetry classes. It enables the characterization of tilted transversely isotropic and tilted orthorhombic media, opening new opportunities for anisotropic model building.     

Distribution of n-alkanes in lacustrine sediments from subtropical lake in Brazil

The geochemical composition of lake bed sediments of a tropical reservoir in Brazil have been investigated. The C, N, P composition showed almost no variation between the different sampling points. All samples contained inorganic phosphorus (IP) ranging from 52 to 70%. The Redfield ratios show that the lake is without significant anthropogenic inputs. Most of the organic matter is composed by higher plants decomposition revealed in the total organic carbon (TOC):nitrogen (N) ratio ranging from 15.4 to 57.2. Also, TOC:organic phosphorus (OP) ranged from 265.3 to 933.6, suggesting that most of the organic matter has terrestrial origin from wood plants. The ratios suggest that most organic matter is influenced by the terrestrial characteristics of the watershed. Furthermore, the relative abundance of n-alkane homologues in the sediments was investigated. All samples have been analyzed for the n-alkanes C₈ to C₄₀. The sediments were dominated by n-alkanes C₂₅–C₃₈. It is concluded that n-alkane in sediments mainly come from terrestrial plants, however there is a contribution of submerged aquatic plants, especially in those sites in deeper areas of the lake. On the basis of P_aq index the n-alkanes in sediments comes from terrestrial plants, however there is a contribution from emerged/floating plants.The investigated lake seems to be considered as meso to eutrophic. Because of the relatively high primary productivity in the lakes, there is a substantial organic-matter flux to the sediments, which rapidly becomes anoxic. According to the pristine/phytane ratio the sediments represent an anoxic environment with values found between 0.38 and 1.72.     

Anomalously deep earthquakes in northwestern Italy

It is usually assumed that earthquakes in intraplate regions occur in the upper crust, and northwestern Italy is generally assigned to this kind of normal seismicity. In this work, the depth distribution of the events localized in this area by the Istituto Geofisico Geodetico (IGG) seismic network in the period 1991–1997 is analyzed in detail. In particular, the location capability of the network is discussed, adopting as reference quarry blasts (for the epicentral position) and the locations obtained from a dense temporary network (for the depth estimate). Within the so-obtained error limits, the depth distribution of events show a characteristic pattern: while for most of the area covered by the network the well-located seismicity lies within the first 20 km of depth, in a band following the inner arc of the Western Alps, numerous events have anomalously large focal depths, reaching a maximum of 114 km. These depth determinations cannot be attributed to instabilities of the location procedure: different choices of the propagation models used for the hypocentral determination led to very similar depth values, always significantly larger than the standard values for the surrounding areas. A strong correlation has been found between the 3-dimensional distribution of these foci and the P-wave propagation anomalies obtained from tomographic studies, suggesting a direct link between elastic and rheological properties of lower crust and upper mantle in this area.