Site Response for the RSM Seismic Network and Source Parameters in the Central Apennines (Italy)

—?Thirty-three earthquakes which occurred in the Central Apennines (Italy) with Ml ranging from 2.4 to 3.7 have been spectrally analysed using digital recordings from twelve stations of the Rete Sismometrica Marchigiana (RSM) network. Data corrected for geometrical spreading and quality factor Q have been inverted by means of the Generalised Inversion Technique. Site responses have been compared with those obtained by H/V ratio. Site amplifications have been observed both at stations placed on Pleistocene sediments and at one station located at 1800?m altitude. Source parameters have been calculated by fitting the spectra with an automatic procedure adopting the ω² source model. The seismic moments range from 9.23?×?10¹⁹ to 4.28?×?10²¹ dyne-cm with an average M ₀ (S) to M ₀ (P) ratio of 1.13?±?0.38. The stress drops are generally low and they vary between 1.1 and 10.2?bar when estimated by using S source spectra, and between 0.5 and 7.1?bar when the P-source spectra are fitted. For the considered range of seismic moments we observe that the stress drop does not have significant dependence on event size.     

Stress Changes due to Recent Seismic Events in the Central Apennines (Italy)

The Central Apennines, Italy, are characterized by moderate seismic activity on normal faults, oriented in directions parallel to the Apenninic chain. The subject of this study is the Umbria-Marche Apennines, a segment approximately 200-km long, where three main seismic events occurred in the last three decades. The 1979 Norcia earthquake was a M_w = 5.8 event, taking place at the south end of the considered segment. The 1984 Gubbio earthquake was a M_w = 5.6 event which took place at the north end. The 1997-1998 Colfiorito sequence constituted 8 main shocks with magnitudes M_w between 5 and 6 and epicenters located between the Gubbio and the Norcia earthquake areas. A model made of an elastic half-space is considered, in which the seismic sources are represented by rectangular dislocations which have the appropriate values of source parameters, and in which the static stress field produced by each event is calculated. The analysis of the Coulomb stress change (ΔC) as a function of time shows that the coseismic stress transfer and fault interaction played an important role in the region during the past three decades: 7 earthquakes of the 9 considered took place where ΔC>0. Such an interaction has been confirmed by the analysis of the aftershocks in the Colfiorito zone post September 26, 1997: about the 61% of the aftershocks considered took place where ΔC>0. The comparison between the ΔC^s due to the coseismic stress transfer and the rate ΔĊ^t due to the tectonic stress allows us to quantify the time advance of the earthquakes. The ΔC^s pattern shows positive values in two areas that can be regarded as historical seismic gaps.     

Aquifer-induced Seismicity in the Central Apennines (Italy)

—The Gran Sasso chain (Central Apennines, Italy) contains one of the largest aquifers of Central Italy. From 1970–1986 the massif was tunnelled through in order to build up a highway and an international underground laboratory for nuclear physics research. These works have strongly modified the hydrogeological situation of the chain, as shown by the decrease in flow rate that occurred in many springs located at the border of the carbonatic structure, along the boundary between the permeable limestone of the massif and the surrounding aquicludes. The analysis of the seismicity (M≥ 3.0) that occurred in the Gran Sasso area from 1956 to 1995 suggests that after the tunnelling works both the number of earthquakes has increased and epicenters have migrated, gathering at the northwestern border zone. The foremost events which occurred in this zone in recent years took place on May 5, 1992 (M = 3.1), August 25, 1992 (M = 3.9) and March 13, 1994 (M = 3.5). The flow rate data of four springs and water level data of an underground karst pool located at the border of the carbonatic structure of the massif show clear anomalies before the occurrence of the quoted earthquakes. Regardless, these anomalies can be explained by the rapid melting of the thick mantle of snow on the Gran Sasso chain, due to sudden increases of mean temperatures. In this paper we present and discuss the possibility that the quoted earthquakes are induced by the irregular variations of the Gran Sasso aquifer, evidenced by the quoted anomalies in the flow rate and water level.     

Hypocenter location accuracy and seismicity distribution in the Central Apennines (Italy)

In this paper, we analyse the seismicity distribution in the CentralApennines (Italy) using the recordings of the Rete SismometricaMarchigiana (RSM). In particular, the selected events are relocated usinga 1-D model calculated by means of an inversion procedure. The robustnessof the 1-D model and the location accuracy are tested. The capability ofthe RSM to well constrain crustal and subcrustal events in the studied areais discussed. We find that in the inner side of the chain the seismicity liesin the upper crustal layers, following the structural trend of the Apenninicbelt. A W-deepening of the events is observed in the Apenninic foredeep,where the seismicity is mainly confined in the lower crust. This evidenceimplies the deepening of the brittle to ductile transition. Some well-locatedsubcrustal events are found. Their locations seem to confirm the W-dippingsubduction of the Adriatic lithosphere beneath the Apennines.     

Fault zone characteristics and scaling properties of the Val d’Agri Fault System (Southern Apennines, Italy)

Systematic in situ analysis of active fault zones in Val d’Agri (southern Italy) suggests that the acquisition of quantitative data on fault-related structural discontinuities is fundamental for (1) discriminating between the latter and older regional features, (2) defining the architecture and related permeability structure of faults, and (3) performing simple statistics in order to evaluate the validity range within which fault characteristics may be considered to be fractals. This type of information can be integrated with regional seismotectonic analysis in order to asses the present-day conditions of deformation characterising the area, and to constrain the possible kinematics of the seismogenic structures controlling earthquake activity in this sector of the southern Apennines. The results of our study are also of interest for modelling earthquake sources, since a knowledge of the permeability structure and scaling properties of a faulted rock volume is potentially relevant for simulating the time and space dependent behaviour of fault zones during a seismic cycle.     

Possible creep-related tilt precursors obtained in the Central Apennines (Italy) and in the Southern Caucasus (Georgia)

Daily averaged tilt component data from two sites of the Central Apennines (Italy) and of the Southern Caucasus (Georgia), respectively, revealed intermediate-term tilts as possible precursors to earthquakes (M=3.0÷4.7) which occurred in the above-mentioned seismic areas within a distance of 50 km from the sites. A good temporal correlation as well as a fair spatial correspondence between these residual tilts (with amplitude and duration of some microradians and months, respectively) and main shocks were pointed out, by removing both secular trends and seasonal thermoelastic effects from the raw tilts. An attempt was made to justify the above-mentioned results, based on the assumption that the observed intermediate-term preseismic tilts are the manifestation of aseismic creep episodes of comparable duration in the fault materials of thrust faults close to the tilt sites. The mechanism refers to a strain field slowly propagating from the preparation focal area to the tilt site, through crustal blocks separated by weak transition zones. This propagation is thought to be the cause of the local aseismic fault slip recorded by the tiltmeters. Previously, both discrete structures and strain propagation effects were revealed in the Central Apennines and are thought also to exist in the Southern Caucasus. As in the past, the rheological properties of fault materials are revealed as viscoelastic ones. In fact, creep equations obtained by applying several viscoelastic models on our data, proved to fit quite well some of the observed tilt precursors, producing viscosity and rigidity values very similar to those reported in literature.Professor Petr Viktorovich Manjgaladze died during the writing of this paper     

Quaternary faults and seismicity in the Umbro-Marchean Apennines (Central Italy): evidence from the 1997 Colfiorito earthquake

Analyses of structural and geomorphological data combined with remote sensing interpretation confirm previous knowledge on the existence of an extensional Quaternary tectonic regime in the Colfiorito area (Umbro-Marchean Central Apennines). This is characterized by a maximum principal axis of finite strain oriented approx. NE–SW, which is the result of a progressive deformation process due to pure and radial extension. Surface geological data, the crustal tectonic setting (reconstructed using a CROP 03 seismic reflection profile), and seismological data relative to the autumn 1997 Colfiorito earthquake sequence constrain the following seismotectonic model. We interpret the seismogenic SW-dipping low-angle normal fault pictured by seismic data as an inverted thrust ramp located in the basement at depth between 5 and 10 km. The surface projection of this seismogenic structure defines a crustal box within which high-angle normal faults are responsible for the deformation of the uppermost crust. The regional patterns of pre-existing basement thrusts therefore control the seismotectonic zoning of the area that cannot be directly related to the high-angle normal fault systems which cut through different crustal boxes; the latter system records, in fact, re-shear along pre-existing normal faults. Moreover, Quaternary slip-rates relative to high-angle normal faults in the Central Apennines are closely related to seismic hazard within each crustal box.     

Seismic characterization and monitoring of Fucino Basin (Central Italy)

The Fucino basin (Central Italy) is one of the largest intramountain alluvial plain in the Apennines range. It has a tectonic origin related to the presence of important systems of faults located in its northern and eastern edges. Some of these faults are still active and capable of generating strong seismic events. Site effects related to the soft soils filling the basin can be very important. In this paper we show the preliminary results of a seismic network installed in the Fucino area in order to collect information about site amplification effects and geometry of the basin. We analyze ambient seismic vibrations and recordings of about 150 local earthquakes mainly related to the seismic sequence of the April 6th 2009 Mw 6.3 L’Aquila event. Moreover the strongest events of L’Aquila sequence were analyzed at the three permanent strong-motion stations operating in the area. Using standard spectral techniques we investigate the variation of resonance frequencies within the basin. The ground motion recorded in the Fucino plain is mainly characterized by strong energy at low-frequencies (f < 1 Hz) affecting both horizontal and vertical components. This is particularly evident for stations deployed in correspondence of very thick deposits of sedimentary filling, where a significant increase of ground-motion amplitude and duration is likely caused by locally generated surface waves. The amplification at low-frequencies (<1 Hz) on the horizontal components can reach up a factor of 10 in comparison to nearby stiff sites. However, we found evidences of seismic amplification phenomena also for stiff sites surrounding the basin, including stations of the Italian strong motion network. The independent geological information and the shallow shear-velocity profiles available for the basin can be combined with resonance frequencies for deriving representative geological sections to be used as base for future numerical 2D–3D modeling of the basin.     

The Holocene paleoseismicity of the Aremogna-Cinque Miglia Fault (Central Italy)

Geomorphic and trench investigations are used toanalyze the seismic potential of the Aremogna-CinqueMiglia fault, an active N- to NW-trending, W-facingnormal fault located in Central Apennines. Wereconstructed a complex 16 km-long, as much as 6m-high, fault scarp that displaces late Holocenesediments in the Aremogna and Cinque Miglia basins.The complex surface expression of the fault, withdouble sub-parallel scarp sections, a change in strikeof about 40^° and local complexity showingimportant horizontal component, appears to becontrolled by the presence of older tectoniclineaments. We opened two trenches across the faultscarp, used a quarry exposure, and reinterpreted atrench opened by Frezzotti and Giraudi (1989), to findthe geological evidence for three Holocene surfacefaulting earthquakes on the Aremogna-Cinque Migliafault. Based on radiocarbon dating and stratigraphicand climatic considerations timing of the events isconstrained between 800 B.C. and 1030 A.D., between3735 and 2940 B.C., and between 3540 and 5000 B.C.. The most recent event is not reported in the twomillennia-long Italian Catalogues of HistoricalSeismicity. We suggest that the most recent eventcould be one of the Middle Age earthquakes of unknownorigin for which several felt reports exist in Rome.Moreover, we also consider the hypothesis that one ofthe shocks of the ambiguous September 1349 earthquakesequence could be the Aremogna-Cinque Miglia mostrecent event. Anyway, based on historicalconsideration we indicate A.D. 1349 as the youngestpossible age for this event. Finally, we suggest theAremogna-Cinque Miglia fault is part of the easternsecondary Apennines seismogenic belt. The faultparameters we obtain for this fault (i.e., recurrence interval longer than 2000 yr, verticallong-term slip rate of 0.3–0.5 mm/yr and m 6.5–6.8 forthe event) can be used as a first hand reference tocharacterize the seismic behavior of other faultsalong this section of the Apennines.     

Uncertainties in probability of occurrence of strong earthquakes for fault sources in the Central Apennines,Italy

Using the characteristic earthquake model, we calculate the probability of occurrence of earthquakes M _w > 5.5 for individual fault sources in the Central Apennines for the 30-year period (2007–2037). We show the effect of time-dependent and time-independent occurrence (Brownian passage time (BPT) and Poisson) models together with uncertain slip rates and uncertain maximum magnitudes and, hence, uncertain recurrence times. In order to reduce the large prior geological slip rate uncertainty distribution for most faults, we obtain a posterior slip rate uncertainty distribution using a likelihood function obtained from regional historical seismicity. We assess the uncertainty of maximum magnitude by assuming that the uncertainty in fault width and length are described by a normal distribution with standard deviation equal to ±20% of the mean values. We then estimate the uncertainties of the 30-year probability of occurrence of a characteristic event using a Monte Carlo procedure. Uncertainty on each parameter is represented by the 16th and the 84th percentiles of simulated values. These percentiles bound the range that has a 68% probability of including the real value of the parameter. We do these both for the Poisson case and for the BPT case by varying the aperiodicity parameter (α value) using the values 0.3, 0.5, and 0.7. The Bayesian posterior slip rate uncertainties typically differ by a factor of about 2 from the 16th to the 84th percentile. Occurrence probabilities for the next 30 years at the 84th percentile typically range from 1% to 2% for faults where the Poisson model dominates and from 2% to 21% where one of the BPT models dominates. The uncertainty in occurrence probability under the time-dependent hypothesis is very large, when measured by the ratio of the 84th to the 16th percentile, frequently being as much as two orders of magnitude. On the other hand, when measured by standard deviation, these standard deviations range from 2% to 6% for those faults whose elapsed time since previous event is large, but always 2% or less for faults with relatively recent previous occurrence, because the probability of occurrence is always small.     

Seismic expression of active extensional faults in northern Umbria (Central Italy)

This paper deals with the geometry and kinematics of the active normal faults in northern Umbria, and their relationship with the seismicity observed in the area. In particular, we illustrate the contribution of seismic reflection data (a network of seismic profiles, NNW–SSE and WSW–ENE trending) in constraining at depth the geometry of the different active fault systems and their reciprocal spatial relationships. The main normal fault in the area is the Alto Tiberina fault, NNW trending and ENE dipping, producing a displacement of about 5 km, and generating a continental basin (Val Tiberina basin), infilled by up to 1500 m with Upper Pliocene–Quaternary deposits. The fault has a staircase trajectory, and can be traced on the seismic profiles to a depth of about 13 km. A set of WSW-dipping, antithetic faults can be recognised on the profiles, the most important of which is the Gubbio fault, bordering an extensional Quaternary basin and interpreted as an active fault based on geological, geomorphologic and seismological evidence. The epicentral distribution of the main historical earthquakes is strictly parallel to the general trend of the normal faults. The focal mechanisms of the major earthquakes show a strong similarity with the attitude of the extensional faults, mapped at the surface and recognised on the seismic profiles. These observations demonstrate the connection between seismicity in the area and the activity of the normal faults. Moreover, the distribution of the instrumental seismicity suggests the activity of the Alto Tiberina fault as the basal detachment for the extensional tectonics of the area. Finally, the action of the Alto Tiberina fault was simulated using two dimensional finite element modelling: a close correspondence between the concentration of shear stresses in the model and the distribution of the present earthquakes was obtained.     

Defining a model of 3D seismogenic sources for Seismic Hazard Assessment applications: The case of central Apennines (Italy)

Geology-based methods for Probabilistic Seismic Hazard Assessment (PSHA) have been developing in Italy. These methods require information on the geometric, kinematic and energetic parameters of the major seismogenic faults. In this paper, we define a model of 3D seismogenic sources in the central Apennines of Italy. Our approach is mainly structural-seismotectonic: we integrate surface geology data (trace of active faults, i.e. 2D features) with seismicity and subsurface geological–geophysical data (3D approach). A fundamental step is to fix constraints on the thickness of the seismogenic layer and deep geometry of faults: we use constraints from the depth distribution of aftershock zones and background seismicity; we also use information on the structural style of the extensional deformation at crustal scale (mainly from seismic reflection data), as well as on the strength and behaviour (brittle versus plastic) of the crust by rheological profiling. Geological observations allow us to define a segmentation model consisting of major fault structures separated by first-order (kilometric scale) structural-geometric complexities considered as likely barriers to the propagation of major earthquake ruptures. Once defined the 3D fault features and the segmentation model, the step onward is the computation of the maximum magnitude of the expected earthquake (M _max). We compare three different estimates of M _max: (1) from association of past earthquakes to faults; (2) from 3D fault geometry and (3) from geometrical estimate corrected by earthquake scaling laws. By integrating all the data, we define a model of seismogenic sources (seismogenic boxes), which can be directly used for regional-scale PSHA. Preliminary applications of PSHA indicate that the 3D approach may allow to hazard scenarios more realistic than those previously proposed.     

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.     

Fault growth by segment linkage in seismically active settings: Examples from the Southern Apennines, Italy, and the Coast Ranges, California

Most models for fault growth and scaling are based on analysis of faults which display dip-slip (i.e. reverse, normal) and strike-slip kinematics; by contrast, little information is derived from faults displaying oblique-slip kinematics. Observations on mesoscopic transpressional faults from the Salinian Block of California and transtensional faults from the Southern Apennines of Italy reveal a complex kinematic history of fault propagation. Faults initially nucleate as isolate segments, which are later kinematically and mechanically linked via development of diffuse deformation zones and/or localised oblique connecting splays. The geometry of observed mesoscopic faults is similar to that of the host, larger structures, thus suggesting that the produced fault patterns are scale independent. Moreover, the overprinting relationships among minor fault-related fabrics permit to define a relative chronology within fault arrays, thus enabling a general sequence of structural stages to be correctly established. Based on minor fabrics and their overprinting relationships, a kinematic deformation model of fault growth by segment linkage is presented, which may have a wide applicability in the field of seismic hazard evaluation.     

Methods for determining the Pleistocene–Holocene component of displacement on active faults reactivating pre-Quaternary structures: examples from the Central Apennines (Italy)

In the tectonically active fault system of the Central Apennines Ridge (Italy) several normal/transtensive faults mapped as Quaternary structures show evidence of pre-Neogene thrusting activity. Therefore, determining the amount of fault-slip during Pleistocene–Holocene times is crucial for seismic hazard assessment. Three principal lines of evidence have been used in this study for estimating the pre-Quaternary activity of extensional fault systems in the Central Apennines. The first is the geological and stratigraphic record (i.e. thickness and facies variations) for pre-Quaternary successions. The second is the widespread structural analysis evidence of extensional faults involved in Neogene thrust tectonics; whether considering a passive role (e.g. fault carried and tilted above a thrust-sheet) or inversion (e.g. positive inversion and related minor structures). The third is the geomorphological evidence, particularly erosion surface analysis, which permits the understanding of the relief generated by tectonics after the formation of post-thrusting erosional surfaces. Preliminary results from some faults which belong to the Sibillini Mts. and the Norcia extensional fault zones show clear evidence of Quaternary reactivations.     

Detachments and normal faulting in the Marche fold-and-thrust belt (central Apennines, Italy): inferences on fluid migration paths

In the outermost domains of the central Apennines fold-and-thrust belt, the structural architecture of the late Miocene–early Pliocene contractional edifice was controlled by competence contrasts in the Calcareous–Marly sequences of Mesozoic–Tertiary age, and by a different state of lithification of the rock units at the onset of deformation. Field data on relative chronology of outcrop-scale structures (cleavage, veins, faults, folds) are presented for the three largest thrust-ramp anticlines of the Marche fold-and-thrust belt: Monte Gorzano, Acquasanta and Montagna dei Fiori-Montagnone. The data show that the timing and geometry of deformation structures differ for: (1) the lower Calcareous interval of late Triassic–early Cretaceous age (LCI) bounded on top by the intermediate detachment (ID) of the Fucoidi Marls; (2) the upper Calcareous–Marly interval (UCMI) of late Cretaceous–Oligocene age; (3) the uppermost detachment zone (UDZ) of lower–middle Miocene age; (4) the topmost Messinian Flysch sequence (FS). In the UDZ early episodes of deformation are manifested by compaction of a poorly lithified sequence followed by pervasive development of layer-parallel pressure-solution cleavage. Reverse faults ramp obliquely across the stratigraphic sequence, and are coated by multiple overgrowths of calcite fibers. These structures are deformed by large, eastward-verging asymmetric folds with N–S axial trends, and are cut by late generations of reverse faults. Normal faults started to develop in the fold backlimbs during the final stages of shortening, in middle–late Pliocene times. These early normal faults were reactivated during episodes of late Pliocene–Pleistocene extensional downfaulting, and are now superposed on the compressional edifice. The UDZ is interpreted to have temporarily sealed the upward escape of fluids during the initial episodes of shortening. Pervasive interlayer flow in the poorly lithified sequence was responsible for development of broken beds and scaly fabrics, similar to those observed in accretionary prisms. Only in the latest stages of deformation did propagation of discrete faults provide an interconnected pathway for fluid migration, until the final offset of the UDZ. The structural relationships suggest that fluids trapped within the fold cores and sealed by the UDZ were finally driven upwards due to progressive disruption of the thrust belt by late normal faults of late Pliocene to Pleistocene and Holocene age. Large-scale fluid migration along structurally-controlled pathways was enhanced by the strong components of uplift consequent to the final stages of deformation in the Marche fold-and-thrust belt, and was eventually associated with episodes of normal seismic faulting.     

Monte Carlo Inversion of DInSAR Data for Dislocation Modeling: Application to the 1997 Umbria-Marche Seismic Sequence (Central Italy)

— The study of surface deformation due to seismic activity is often made using dislocations with uniform slip and simple geometries. A better modeling of coseismic and postseismic surface displacements can be obtained by using dislocations with variable slip and nonregular shapes. This is consistent with the asperity model of fault surfaces, assuming a friction distribution on faults made of locked zones with much higher friction than surrounding zones. In this paper we consider the 1997–1998 Colfiorito seismic sequence. The coseismic surface displacements in the Colfiorito zone are used in order to infer the slip distribution on the fault surface at different stages of the sequence. The displacement field has been modeled varying the slip distribution on the fault, and comparing the deformation observed by SAR and GPS techniques with model results. The slip distribution is calculated by Monte Carlo simulations on a normal fault with the dip angle equal to 40°. A good approximation is obtained by using square asperity units of 1.5×1.5 km². In the first stage, we employed a simplified model with uniform slip, in which each asperity unit is allowed to slip a constant amount or not to slip at all, and in the second stage, we evaluate the slip distribution in the dislocation area determined by the Monte Carlo inversion: in this case we allow unit cells to undergo different values of slip in order to refine the initial dislocation model. The results show that the 1997 seismic events of the sequence can be modeled by irregular dislocations, obtaining a good fit to the DInSAR and GPS observations. The model also confirms the results of previous studies by a different methodology, defining the distribution of asperities on the fault plane using the fault geometry, the geodetic data and the seismic moment of the 1997–1998 Colfiorito seismic sequence. Furthermore, the analysis of 1997 aftershocks in the seismogenic region shows a strong correlation between most events and the asperity distribution, which can be considered as an independent test of the validity of the model.     

Seismic Source Classification by Means of a Sonogram- Correlation Approach: Application to Data of the RSM Seismic Network (Central Italy)

—?In this paper, we apply a new sonogram correlation approach for seismic source classification. In order to test the technique, we use the waveforms recorded by 6 among 16 of the stations of the Rete Sismometrica Marchigiana (RSM) network, relevant to the vertical component of 365 earthquakes that occurred in the analyzed area from April 1996 until September 1999. All the sonograms of each station are cross-correlated in a single station procedure and the calculated degrees of resemblance are used to classify the events by means of a bridging technique. The results obtained with the single-station approach are integrated into a network classification. We identified three families in the investigated area. Once the efficiency of the method is tested, we analyzed the families obtained for the selected subset of events. Considering the temporal and the spatial distribution of the events forming each family, we classified them as quarry blasts. This classification, even if relevant to a subset of earthquakes occurring in Central Italy, modifies the present knowledge of the seismicity distribution in the area and allows estimation of the errors affecting the hypocenter locations provided by the RSM network. An extension of this new technique to the whole Central Italy earthquake catalog promises an interesting review of the seismicity in the area.