Vol49_2_2006 853 ANNALS OF GEOPHYSICS, VOL. 49, N. 2/3, April/June 2006 Key words Seismic Bulletin – magnitude – com- pleteness analysis – errors of hypocentral locations – detection capacity 1. Introduction During the first five years of the nineties, the Istituto Nazionale di Geofisica (nowadays INGV) was involved in a European project named «EUROBULL». The aim of the EUROBULL group was to compile a high quality European bulletin to be used both as a reference, to eval- uate bulletins produced by international net- works, and for scientific purposes. Italy, France and Sweden collected data from about fifteen countries and, starting from their own bulletin procedures, processed the da- ta independently from each other. I was directly involved in the software imple- mentation and Giuseppe Smriglio, a very helpful person who was the manager of the Data Centre, let me modify the Italian bulletin procedures adapting them to the European purposes. Moreover, the EUROBULL group produced a series of control routines to test the quality of the new bulletin, so we had the opportunity to apply the same test to the Italian bulletin. This work, that analyzes sixteen years of IN- GV bulletins, recalls the efforts to improve the network and the quality of the data for which Giuseppe was one of the principal promoters. 2. Database Bulletins produced by INGV between 1988 and September 2003, are the database of the present work. Figure 1 shows the area analyzed and the INGV National Centralised Seismic Net- work (RSNC) configuration used in this paper. Figure 2 illustrates the evolution of the RSNC during the analyzed period. Colors of the stacked Evaluation of sixteen years of INGV seismic bulletins Rosalba Di Maro Istituto Nazionale di Geofisica e Vulcanologia, Roma, Italy Abstract This work analyzes earthquake parameters published in the Seismic Bulletin of the Istituto Nazionale di Geofisi- ca e Vulcanologia (INGV). The analysis and interpretation of the digital signals, done by specialist employees daily, produce most of the seismological information that comprises INGV’s earthquake bulletins. After a brief introduction on the criteria we use to obtain seismic parameters, this paper will review the processing procedures employed over a period of sixteen years from 1988 to September 2003.This study also addresses the issue of the comparison between magnitude calculated on signal duration (Md) and on amplitude (ML) and the lack of a cor- rect calibration between them. A completeness analysis of the whole bulletin performed using both the Stepp and Habermann techniques shows the importance of considering changes in the seismicity rate and in the geom- etry of the seismological network. To conclude this excursus, we calculated the errors of hypocentral locations and the detection capacity of individual seismic stations. This final step stresses the increasing improvement of the INGV seismic network over the past 16 years. Mailing address: Dr. Rosalba Di Maro, Istituto Nazio- nale di Geofisica e Vulcanologia, Via di Vigna Murata 605, 00143, Roma, Italy; e-mail: rosalba.dimaro@ingv.it 854 Rosalba Di Maro histograms are related to the years, the different amplitude of the columns is related to the differ- ent percentage working that year (since January 1988, the INGV shift workers draw up a daily re- port on the efficiency of each seismic station). Stations are in descending order based on lati- tude, and the stations reported in the four his- Fig. 1. INGV seismographic network configuration. The star symbols indicate the station locations, the red plus symbols the closed stations. The four colors of the stars are related to four sectors: north, center-north, center-south, south. Fig. 2. Evolution of the seismographic network du- ring the analyzed period. Colors of the stacked histo- grams are related to the years, the different amplitu- de of the columns is associated to the different per- centage working that year. Stations are in descending order based on latitude, and the stations reported in the four histograms are indicated, in the previous fi- gure, by the four sectors respectively. 1 2 855 Evaluation of sixteen years of INGV seismic bulletins tograms are indicated, in the previous figure, by the four sectors respectively. It is easy to identify from the color bars closed stations and new ones, most of which were installed in Central Italy. 3. Local magnitude and duration magnitude Before processing the database, a brief men- tion of the magnitude we used in our determi- nations is useful. For local events we apply the magnitude-duration relation (Console et al., 1988) (3.1) where a = 2.0, b = 0.870 and c = 0.082. τi is the i- th seismic signal duration and ∆i is the linear ( )logM a c bd i i10i )= + -τ ∆ distance of the i-th station compared to the epi- center determination. The (3.1) is derived from the observation that the absolute time (τi + c∆i) is constant, in other words independent from the station. Most of the time, the automatic ac- quisition system calculates τi determining the end of the signal when its amplitude comes back to the noise level. When the signal dura- tion exceeds 8-min (that is the maximum acqui- sition time) the operator evaluates the duration directly from the analog acquisition that has a continuous registration on the paper. The flowchart in fig. 3 represents the algo- rithm that attributes Md value to a regional event. This quantity is calculated as arithmetic average among all the values. The maximum and (/or) the minimum are eventually removed if they are «very distant» from the average. The algorithm that attributes a local magni- tude to an event (fig. 3), is the same as the one for Md. In this case we eliminate values less than 1.5. The Richter formula is used for sta- tions whose distances, from the epicenter, are less than 600 km (3.2) The distance function f (∆) is the same published by Richter in 1935 (Richter, 1935). On the other hand, the deduced attenuation function for Italy closely corresponds to Richter’s one, as describ- ed by Gasperini (2002). At the present time the RSNC is mostly com- posed of short period vertical seismometers (S-13 GEOTECH). For this reason, in eq. (3.2), we use the maximum amplitude picked up automatically on the vertical seismic signal and corrected for the instrumental response, as described below. If Asp is the short period maximum ampli- tude and Hsp the transfer function of the whole instrumental chain, we can obtain the Wood Anderson amplitude (nm) by (3.3) where HWA is the Wood Anderson transfer func- tion. For regional events (600-1600 km) we ap- ply a modified Nuttli formula (Nuttli, 1973) A A H H i i sp i sp WA = ( ) ( ) .logM A A f 2l NS EW i 10i = + + ∆ Fig. 3. Flowchart of the algorithm that attributes the Md value to a regional event. This quantity is cal- culated as the arithmetic average eventually remo- ving the maximum and (/or) the minimum, if they are «very distant» from the average. 856 Rosalba Di Maro (3.4) The original formula was modified by means of regression with Italian data (1988-1993). 4. Data processing 4.1. Discrepancy between the duration and local magnitude In this paper we processed sixteen years of data: from January 1988 to September 2003 (in- clusive). The INGV bulletin reports, for that pe- riod, 42154 locatable events («locatable» means that the event has been located using at least six phases). Among them, 23595 events have both local and duration magnitudes and 9594 have at least five Md values. Figure 4 illustrates the annu- al trend of these quantities. ML and Md are very different from each other most of the time. Figure 5a,b illustrates this dis- crepancy: the first displays the frequency of the events for each class of magnitude Md and ML. It is easy to notice that the ML histogram is shifted more to the left compared to the Md ones. The dif- .log+. ( )log2 34 ) ∆+.4 37-=( )M T A b i i i i 10 10LG c m ferent shape of the two histograms, provided by the same event data-set, shows that the Italian seismicity in the last sixteen years is concentrat- ed between 2.5 and 3.2 considering Md, but under 2.5 if we consider ML. Figure 5b illustrates that Md and ML are not univocally determined. Also the most probable value seems to be the medium one, considering Md in lack of ML creates a non homogeneous data set. This ambiguity is due to the method of deriv- ing Md. The magnitude-duration relation (3.1) Fig. 4. Annual trend of Italian seismicity. Fig. 5a,b. a) Discrepancy between the frequency of the events for class of magnitude Md and ML. b) Three-di- mensional vision of the correspondences between Md and ML. a b 857 Evaluation of sixteen years of INGV seismic bulletins has been obtained analyzing 15 months of data that are probably not enough for a correct deter- mination of parameters. Another problem is re- lated to the new formulation given by (3.1) of the absolute time. To validate the hypothesis (4.1) where τs is the absolute duration; ∆s is the total distance covered by waveforms that travel with velocity v; τi is the i-station signal duration; Ti is the i-station arrival time; T0 is the origin time. We calculated τs and the related standard deviations for the whole data set. Figure 6 illustrates this trend, the biggest magnitude values, that have the biggest number of recordings, are associated to high values of standard deviations. This means that the absolute duration is not independent from the station and (4.1) is not validated. 4.2. Data set completeness To determine the completeness of our data set we apply the Stepp method (Stepp, 1972). This model assumes that: – Earthquake distribution can be modeled by a Poisson distribution whose standard devia- tion is ( ) const v T Ts s i i 0= = + - =τ τ ∆ (4.2) where λ is the average of the occurrences in the period T. – The seismic activity for the whole ob- served period is stationary i.e. λ= cost. Let’s go back to our data set and consider first the local magnitude. Observing the seismic- ity year by year, four steps of local magnitude and duration magnitude have been selected Md ≤ 2.4 ML ≤ 1.8 2.5 ≤ Md ≤ 2.8 1.9 ≤ ML ≤ 2.4 2.9 ≤ Md ≤ 3.4 2.5 ≤ ML ≤ 3.0 Md ≥ 3.5 ML ≥ 3.1. Figures 7a and 8a show the cumulated number of events during the 16-year period, whereas figs. 7b and 8b show, in bi-logarithmic scale, s versus T (starting from 2003). According to Stepp’s theory, when the mean rate of occurrences is constant, the related peri- od of time is ‘complete’ and (4.2) has the linear trend shown in the figs. 7b and 8b, as reference, in green. The four steps of magnitude Md have the same «critical» time Tc under which the da- ta set is considered incomplete. The four steps of magnitude ML have two «critical» times, one, in common with Md corresponding to 1997, and the other for ML ≥ 3.1, corresponding to 1999. The explanation for this «critical» year could derive from fig. 2: numbers of stations were added, during 1997, in the middle-south part of Italy. Consequently, an apparent low seismicity has grown during the last period. The increment in the station density makes the data set for the peri- od before 1997 «incomplete». Moreover we can make another comment: Stepp’s approach to the data-set completeness has the «Poissonianity» of the event distribution underlying the hypothesis. Of course clusters of events as well as aftershocks or foreshocks are not independent, so their spatial, temporary distributions are not «Poissonian». In addition we should mention the «Colfio- rito» sequence that contributed with more than 5000 events during the 1997 and 1998. Of course it is possible to make the database «Poissonian» removing the clusters, but this is not the purpose of this paper. = T σ λ λ Fig. 6. Standard deviations of the absolute time (τs) versus Md. The standard deviation values indicate that τs is not constant. 858 Rosalba Di Maro Habermann (1983) proposed another defini- tion for the minimum magnitude of completeness. Beginning from the Gutenberg-Richter (1944) relationship . (4.3) The determination of the minimum magnitude of log N a bM10 = - completeness MC of a data set is obtained plotting (4.3) where N is the cumulative number of events, M the magnitude and a and b are constants. Fig- ure 9a,b shows the frequency-magnitude distribu- tion of Italian earthquakes during these sixteen years. The straight line represents the best fit of the observations in the area of «complete» rec- ording indicated by MC. We can individuate this a b Fig. 8a,b. a) Cumulated number of events during the 16-year period for four steps of local magnitude. b) Stan- dard deviation, versus T (starting from 2003) in bi-logarithmic scale for four step of local magnitude. Fig. 7a,b. a) Cumulated number of events during the 16-year period for four steps of duration magnitude. b) Standard deviation, versus T (starting from 2003) in bi-logarithmic scale for four steps of duration magnitude. a b 859 Evaluation of sixteen years of INGV seismic bulletins values that for ML coincides with the minimum value 1.5 and for Md is 2.3 ± 0.1. Performing the regression (for both rela- tions the correlation coefficient is approximate- ly 1) we can write (4.3) as follows: . 4.3 Hypocentral locations As mentioned before, in these sixteen years of bulletins, our qualified technicians processed more than 40 000 events, merging data from sev- eral seismic networks to obtain a high quality product. The hypocentral determinations are per- formed, using routines belonging to the Hypo71 family. . .log N M1 03 5 8L10 = - + . .log N M1 54 8 0d10 = - + Fig. 9a,b. a) Minimum magnitude of completeness for Md in accordance with Habermann (1983). b) Mi- nimum magnitude of completeness for ML in accor- dance with Habermann (1983). Fig. 10. The first three events of the 2003 INGV Bulletin. a b 860 Rosalba Di Maro Figure 10 shows an extract of the INGV Bul- letin, in particular the first three events of 2003. Most of the events have the location fixed at a certain depth, sometimes the software itself fixes the depth, and sometimes because the depth value is not consistent the technicians fix them manually. By increasing the number of the stations, the number of events with fixed depth is decreasing as the fig. 11 demonstrates. Figure 12a reports all the depths selected among events with the product of coordinate errors (Lat er- ror ∗ Long error ∗ depth error) below 1000 km3. The interpolation is made on a grid of 30 km side. From the figure the very deep location of the events in the southern part of the Tyrrhenian Sea is clear. The depth of most of the other events is less than 40 km, as fig. 12b shows. Regarding the errors on the hypocentral lo- cations, the routine we use gives latitude, longi- tude and depth errors in kilometres instead of the ellipsoid error. To have an idea on the quali- ty of the locations, let us take a look at the error location distribution. Table I shows (in the sec- ond column) the percentage of the events with a certain error indicated in the first column. Most of these events have a very small error, fig. 13a-d show the distribution of the small- medium location error (below 1000 km3) sepa- rated for class of duration magnitude, it is re- markable that greatest errors are related with ar- eas not included in the RSNC and, of course, Fig. 11. Percentage of the events whose hypocen- tral location has fixed (darker) depth. The percenta- ges are calculated over a three year period, as indica- ted on the abscissa. Fig. 12a,b. a) Hypocentral depths map distribution. The events selected have an error on the coordinates (Lat error ∗ Long error ∗ depth error) below 1000 km3. The interpolation is made on a grid of 30 km per side. b) Hy- pocentral depths distribution in class of 10 km each. The events selected have an error on the coordinates (Lat error ∗ Long error ∗ depth error) below 1000 km3. a b 861 Evaluation of sixteen years of INGV seismic bulletins Table I. Evaluation of the error associated to hypo- central locations. The percentage only considers the events whose depth is not fixed. V (km3) Lat error ∗ Long error ∗ depth error Percentage 0 ≤ V < 100 81.6% 100 ≤ V < 1000 15.6% V ≥ 1000 2.8% Fig. 13a-d. a) Distribution of the small-medium localization error (below 1000 km3 ) for events with Md ≤ 2.4. The interpolation is made on a grid of 50 km per side. b) Distribution of the small-medium localization error (below 1000 km3) for events with 2.5 ≤ Md ≤ 2.8. The interpolation is made on a grid of 50 km per side. c) Distri- bution of the small-medium localization error (below 1000 km3 ) for events with 2.9 ≤ Md ≤ 3.4. The interpolation is made on a grid of 50 km per side. d) Distribution of the small-medium localization error (below 1000 km3) for events with Md≥3.5. The interpolation is made on a grid of 50 km per side. the bigger the event the better the hypocentral location will be. To complete the «panorama», fig. 13e sum- marizes everything about the 2.8% of the largest errors. The radius of the circumferences is in- creasing with Md, the color is darkens with the error values and, in the end, the plus sign indi- cates the events which occurred in the last three years. We can notice that: a b c d 862 Rosalba Di Maro – The hypocentral locations with the big- gest error and magnitude are outside the network. – Most are concentrated in the southern part of the Tyrrhenian Sea where the deepest events occur. – The increasing number of the stations in the last three years improved the quality of lo- cations, only few of them, with low magnitude, are inside the network. 4.4. The detection capacity of the seismographic stations This last paragraph is dedicated to the per- formance of the RSNC stations. All the events recorded by each station have been considered, more precisely, all the arrival times that con- tributed to evaluate the hypocentral location have been considered (not all the phases recorded by the automatic system will contribute to the loca- tion process; it depends on the signal/ noise ratio and on the convergence of the location algorithm solution). The events are collected into four steps of duration magnitude, as usual Md ≤ 2.4 2.5 ≤ Md ≤ 2.8 2.9 ≤ Md ≤ 3.4 Md ≥ 3.5. The area around each station has been divided into twelve parts (30° each), each part has been divided into sectors whose dimensions are re- lated to the distance from the center (i.e., the station). Figure 14 explains the criteria. I selected, to be brief, three stations that, from north to the south are: Trieste (TRI), L’Aquila (AQU) and Sicignano (SGO). The se- lection was on the location basis and also be- cause these stations are currently working. The complete analysis of all the stations will be the topic of a specific paper. The results of this partial investigation are reported in fig. 15a-l: 1) SGO has the best coverage for all the steps of magnitude. Fig. 13e. Largest errors (above 1000 km3) distribu- tion. The radius of the circumferences is related to Md value, the color changes with the error values and, in the end, the plus sign indicates the events which occurred in the last three years (see the legend as reference). Fig. 14. Criteria to determine the detection capacity of the seismographic stations. The area around each station has been divided in twelve parts (30° each), ea- ch part has been divided in to sectors whose dimen- sions are related to the distance from the station: for the first 200 km the sectors are 50 km by 50 km, from 200 km to 600 km the distances are 100 km each. The farthest sector has the distance, from the station, deter- mined by the farthest event detected (marked in the fi- gure by a red or white star). The white sectors, even- tually not closed, indicate a lack of events. 863 Evaluation of sixteen years of INGV seismic bulletins Fig. 15a-h. a) Percentage of events (Md ≤ 2.4) detected by Trieste station. b) Percentage of events (Md ≤ 2.4) detec- ted by L’Aquila station. c) Percentage of events (Md ≤ 2.4) detected by Sicignano station. d) Percentage of events (2.5 ≤ Md ≤ 2.8) detected by Trieste station. e) Percentage of events (2.5 ≤ Md ≤ 2.8) detected by L’Aquila station. f) Percentage of events (2.5 ≤ Md ≤ 2.8) detected by Sicignano station. g) Percentage of events (2.9 ≤ Md ≤ 3.4) detected by Trieste station. h) Percentage of events (2.9 ≤ Md ≤ 3.4) detected by L’Aquila station. The criteria are explained in fig. 14. Green dots are the events recorded by the RSNC network in the range of duration magnitude indicated. a b c d e f g h 864 Rosalba Di Maro 2) The farthest event detected is not in a sector with a high percentage of detection. 3) No privileged direction in the detections is visible. 4) A high percentage of detection does not correspond to those areas with a high seismicity. The «strategy» in the performing of hypocen- tral location is one possible explanation for the last point. In fact, in case of an abundant number of phases, the technicians exclude those arrival times belonging to far stations in order to obtain a robust solution for the location algorithm. This could be the reason for the second point also. 5. Conclusions Sixteen years ago, the possibility to acquire and process digital seismic signals gave a new impulse and speeded up the process of opti- mization of the RSNC. This overview emphasizes the effort to im- prove the software and the hardware related to the correct detection and characterization of the earthquakes. Regarding the parameters of the hypocen- tral location the accuracy is not under discus- sion. The only weak link in the chain is, as in- dicated above, the magnitude and the correla- tion between the local and duration magnitude. Nowadays the introduction of the new digi- tal seismographic network, with new sensors, three-components and broadband, will allow a more correct definition of the local magnitude, so we could calibrate the duration magnitude having much more data to compare. REFERENCES CONSOLE, R., B. DE SIMONI and A. DI SANZA (1988): Rie- same della relazione magnitudo-durata, in Proceedings of VII Yearly Meeting of the Gruppo Nazionale di Geo- Fig. 15i-l. i) Percentage of events (2.9 ≤ Md ≤ 3.4) detected by Sicignano station. j) Percentage of events (Md ≥ ≥ 3.5) detected by Trieste station. k) Percentage of events (Md ≥ 3.5) detected by L’Aquila station. l) Percentage of events (Md ≥ 3.5) detected by Sicignano station. The criteria are explained in fig. 14. Green dots are the events recorded by the RSNC network in the range of duration magnitude indicated. i j k l 865 Evaluation of sixteen years of INGV seismic bulletins fisica della Terra Solida (GNGTS), November 30-De- cember 2, Rome, Italy, 51-62. GASPERINI, P. (2002): Local magnitude revaluation for recent Italian earthquakes (1981-1996), J. Seismol., 6, 503-524. GUTENBERG-RICHTER, N. (1944): Frequency of earthquakes in California, Bull. Seismol. Soc. Am., 35, 185-188. HABERMANN, R.E. (1983): Telesesmic detection in the Aleu- tin Island Arc, J. Geophys. Res., 88, 5056-5064. NUTTLI, O.W. (1973): Seismic wave attenuation and magni- tude relations for Eastern North America, J. Geophys. Res., 78 (5), 876-885. RICHTER, C.F. (1935): An instrumental earthquake magni- tude scale, Bull. Seismol. Soc. Am., 25, 1-31 STEPP, J.C. (1972): Analysis of completenessof the earth- quake sample in the Puget Sound area and its effect on statistical estimates of earthquake hazard, in Proceed- ings of the First Microzonation Conference, Seattle, U.S.A., 897-909.