An estimate of hypocentre location accuracy in a large network: possible implications for tectonic studies in Italy

Data recorded by the Italian Telemetered Seismic Network (ITSN) of the Istituto Nazionale di Geofisica (ING) have been widely used in recent years to image slab structures and to find evidence for active processes along the Italian Peninsula. However, the use of seismic data for geostructural purposes may be affected by the well-known trade-off between earthquake location and seismic-velocity parameters. Furthermore, the confidence ellipse predicted by standard procedures may be inadequate for the representation of the probable error of a computed localization. This paper evaluates the probable errors on the hypocentre determinations of the seismic events recorded by the ITSN, using a Monte Carlo method.
We compute synthetic arrival times using a 1-D velocity model appropriate as an average for the Italian area. The hypocentres used are all those recorded by the ITSN during the period January 1992 to March 1994 (1972 events). Station locations are those of the current ITSN configuration. The synthetic arrival times are perturbed with a Gaussian distribution of errors and input to ING's standard hypocentral location procedure, but using crustal velocities differing by 10 per cent from those used to generate them. Each simulation is repeated at least 30 times. Average absolute shifts of hypocentres are assessed in grid cells of linear dimension 33 km covering the whole Italian region.
For regions within the ITSN, shifts are typically 5–10 km in location and up to 20 km in depth. However, for offshore and coastal regions, they are much greater: 50 km or more in both location and depth (far exceeding the equivalent uncertainties quoted by ING bulletins). Possible consequences of this are highlighted by producing a cross-section of subcrustal hypocentres from the Adriatic to the Tyrrhenian Sea, where the large uncertainty in depth precludes any confident interpretation of dipping tectonic features.     

Paleogeography and structure of the central Mediterranean: Sicily and its offshore area

The geology of the mainland and offshore of Sicily is illustrated by a few geologic sections and seismic profiles across the late Cenozoic orogenic belt of central and western Sicily and across the Sardinia Channel and Sicily Straits. This belt is the result of several tectonic events. Deformation involved mainly the sedimentary cover of the old African continental margin characterized by a broad basinal domain, flanked along its external (southern) margin by a shallow-water carbonate platform attached to Africa in the Triassic. Compressional deformation started in the more internal basinal rock assemblages overlying a thinned crust. The most important structural characteristic of the early phase of thrusting is the duplex pile forming the bulk of the chain in Central Western Sicily. The structure consists of a basal allochthon, made up of Permian to Middle Triassic layers, an intermediate duplex wedge, composed of competent Mesozoic carbonates, and a roof complex, including Upper Mesozoic-Lower Tertiary less competent rocks. Large-scale clockwise rotation of the thrusts accompanied transpressional movements in the hinterland during the Pliocene. Right oblique reverse faults modified the previous tectonic contacts between the allochthons in the hinterland zones. Contemporaneous south-directed imbrications affected the southern external areas, progressively incorporating foreland and piggyback basins. Development of the Gela Thrust System appears to be linked to the transpressional event; its accretion is also related to contemporaneous underthrusting at deeper levels of Mesozoic carbonate substratum. The older buried thrust sheets were pushed up to the surface breaching the deformed Tertiary cover of the Gela TS. Northwards in the belt post-Messinian normal growth faults opened half graben whose sedimentary fill underwent structural inversion. Alternation of extension and compression tectonics characterizes the Sicilian continental margin in the last million years.     

Doubtful Events Reported in Current Italian Earthquake Catalogues: The Medieval Earthquakes of Subiaco (Central Italy)

Within the framework of a study of the seismicity of the Aniene Valley (Central Italy), we analysed the medieval earthquakes of Subiaco (1216, 1227, 1299), the largest events reported for the area. Our main goal was to investigate some doubtful events reported in earthquake catalogues and, as such, currently utilised for seismic hazard estimates. A careful screening of the oldest available sources and their filiation pattern up to the present pointed out the uncertainty on the date and nature of these phenomena. A multidisciplinary approach based on the joint analysis of archaeological, geomorphologic and historical evidence allowed us to propose new interpretations concerning these events and their significance for the assessment of seismic hazard in the Aniene Valley. The main conclusion is that the dates of the 1216 and 1227 events are fairly unsupported. In particular, the 1216 earthquake could be dated back to between AD 1159 and 1181.     

Conceptually-based shot noise modeling of streamflows at short time interval

A conceptual-stochastic approach to short time runoff data modelling is proposed, according to the aim of reproducing the hydrological aspects of the streamflow process and of preserving as much as possible the dynamics of the process itself. This latter task implies preservation of streamflow characteristics at higher scales of aggregation and, within a conceptual framework, involves compatibility with models proposed for the runoff process at those scales. At a daily time scale the watershed response to the effective rainfall is considered as deriving from the response of three linear reservoirs, respectively representing contributions to streamflows of large deep aquifers, with over-year response lag, of aquifers which run dry by the end of the dry season and of subsurface runoff. The surface runoff component is regarded as an uncorrelated point process. Considering the occurrences of effective rainfall events as generated by an independent Poisson process, the output of the linear system represents a conceptually-based multiple shot noise process. Model identification and parameter estimation are supported by information related to the aggregated runoff process, in agreement to the conceptual framework proposed, and this allows parameter parsimony, efficient estimation and effectiveness of the streamflow reproduction. Good performances emerged from the model application and testing made with reference to some daily runoff series from Italian basins.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.