Vol49,4_5,2006 1041 ANNALS OF GEOPHYSICS, VOL. 49, N. 4/5, August/October 2006 Key words seismic array – microearthquakes – Central Apennines 1. Introduction Over the past 40 years, much attention has been given to the use of instrumental arrays to detect plane-wave signals in seismological time series and to measure wave slowness vectors (Rost and Thomas, 2002; Schweitzer et al., 2002). Compared to classical seismic networks, the main advantage of seismic arrays consists in their ability to detect weak or emergent signals, and to allow for an effective noise reduction through multichannel waveform stacking. More- over, once combined with data from adjacent networks, the wave azimuths determined at one or several arrays provide strong constraints on the location of low-magnitude events detected by a small number of stations (e.g., Bratt and Bache, 1985). Such reduction in the detection threshold represents a critical step toward a successful as- sessment of the location, geometry and kinemat- ics of active fault systems. Following these arguments, the University of Salerno, the Italian National Institute for Nu- clear Physics and National Institute for Geo- Performances of the UNDERground SEISmic array for the analysis of seismicity in Central Italy Gilberto Saccorotti (1), Bellina Di Lieto (2), Fabrizio Tronca (1), Costantino Fischione (1), Roberto Scarpa (2) and Rosanna Muscente (3) (1) Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Vesuviano, Napoli, Italy (2) Dipartimento di Fisica «E.R. Caianello», Università degli Studi di Salerno, Baronissi (SA), Italy (3) Parco Scientifico e Tecnologico d’Abruzzo, L’Aquila, Italy Abstract This paper presents the first results from the operation of a dense seismic array deployed in the underground Physics Laboratories at Gran Sasso (Central Italy). The array consists of 13 short-period, three-component seis- mometers with an aperture of about 550 m and average sensor spacing of 90 m. The reduced sensor spacing, joined to the spatially-white character of the background noise allows for quick and reliable detection of coher- ent wavefront arrivals even under very poor SNR conditions. We apply high-resolution frequency-slowness and polarization analyses to a set of 27 earthquakes recorded between November, 2002, and September, 2003, at epi- central distances spanning the 20-140 km interval. We locate these events using inversion of P- and S-wave back- azimuths and S-P delay times, and compare the results with data from the Centralized National Seismic Network catalog. For the case of S-wave, the discrepancies among the two set of locations never exceed 10 km; the largest errors are instead observed for the case of P-waves. This observation may be due to the fact that the small array aperture does not allow for robust assessment of waves propagating at high apparent velocities. This informa- tion is discussed with special reference to the directions of future studies aimed at elucidating the location of seismogenetic structures in Central Italy from extended analysis of the micro-seismicity. Mailing address: Dr. Gilberto Saccorotti, Istituto Na- zionale di Geofisica e Vulcanologia, Osservatorio Vesuvia- no, Via Diocleziano 328, 80124 Napoli, Italy; e-mail: gil- berto@ov.ingv.it 1042 Gilberto Saccorotti, Bellina Di Lieto, Fabrizio Tronca, Costantino Fischione, Roberto Scarpa and Rosanna Muscente physics and Volcanology (hereinafter referred to as INFN and INGV, respectively) conducted a joint effort for the deployment of a small- aperture, dense array of short period seismome- ters in the Gran Sasso (Central Italy) under- ground Physics Laboratories. Designed and tested during the late 90’s, the underground ar- ray (UNDERSEIS) became fully operative by early May 2002, providing low-noise, high-dy- namic continuous recordings of ground shaking signals (Scarpa et al., 2004). In Italy, multi- channel deployments have already been suc- cessfully tested for their detection capabilities on local seismicity (e.g., Braun et al., 2004); however, UNDERSEIS represents the first ex- ample of a seismic antenna deployed in an un- derground setting. This paper presents the first results from the operation of the UNDERSEIS array. We first describe the array configuration and setting, and discuss its resolving capabili- ties in detecting and measuring plane-wave sig- nals. We then describe the techniques routinely used for retrieving the propagation and polar- ization parameters from the multichannel recordings. Examples of application are pre- sented for a set of low magnitude (Md < 3) earth- quakes recorded at local and regional distances. S-P delay times and backazimuths of these events are then used to obtain epicentral esti- mates, the robustness of which is demonstrated by comparison with the catalogue locations from the INGV’s Centralised National Seismic Network (CNSN). These results will be dis- cussed further with special reference to the di- rections for future works aimed at a) improving the array performance, and b) investigating the details of wave propagation and the location of low-magnitude events in the investigated area. 2. The UNDERSEIS arrey Located in Central Italy, some 120 km ENE of Rome, the Gran Sasso is a limestone massif Fig. 1. Map of Central Italy with location of the Gran Sasso INFN laboratories (square) and epicenters of the analysed earthquakes (black dots). Sites mentioned in the text are labeled. The inset at the upper rigth shows the configuration of the UNDERSEIS seismic array. Performances of the UNDERground SEISmic array for the analysis of seismicity in Central Italy constituting the highest peak (2914 m) of the Apennines Chain (fig. 1). Since 1986, subnu- clear particle physics laboratories have been hosted in underground tunnel systems located at a mean elevation of 1500 m, being therefore naturally shielded by a 1400-m-thick limestone cap. The area is affected by major seismicity as- sociated with the mainly distensive tectonics af- fecting the Apennines since the Late Pliocene. The major historical event was the 1915, Ms = 6.8 Avezzano earthquake, which caused about 32 000 casualties (Amoruso et al., 1998). In addition, low-to moderate-size earthquakes occur quite frequently, the most recent in 1992, 1994 and 1996 (De Luca et al., 2000; Chiarab- ba et al., 2005). In the same area a temporary, large-scale seismic deployment was installed aimed at defining the location and geometry of seismogenic structures from microseismicity data (Chiarabba et al., 2004). The low noise associated with the under- ground setting (De Luca et al., 1998), and the quite considerable seismicity make the Gran Sas- so labs an ideal site for high resolution seismic observations. The design of the UNDERSEIS hardware and software components began on the late 90’s; the array became fully operative by May, 2002. UNDERSEIS intersects a main seis- mogenetic fault where the occurrence of slow earthquakes has recently been detected through two wide-band geodetic laser interferometers (Crescentini et al., 1999; Amoruso et al., 2002). In its present configuration, the array consists of 13 elements, each equipped with a MARK Prod- uct L4C-3D, 1 Hz 3-component seismometer, whose signals are digitized at a sampling interval of 0.01 s and a dynamic range of 24 bits. Further details on the design and working principles of the array are reported in Scarpa et al. (2004). In its present configuration, UNDERSEIS depicts an average sensor spacing of about 90 m and the widest distance among sensors of about 550 m (fig. 1). Location of the array elements was constrained by the geometry of the tunnel system, thus hindering the achievement of uni- form station spacing and relative azimuths dis- tribution which are warranted for multichannel analyses (fig. 2). Recently, however, new gal- leries have been made available, and installa- tion of additional array elements is already in progress with the specific goal of improving ar- ray aperture and geometry. 3. Array performance In this section we investigate the array’s ca- pabilities and limitations, in turn deriving use- ful hints for the selection of parameters to be adopted for the subsequent frequency-slowness measurements. 3.1. Array response patterns The resolving capabilities of a multichannel system are conveniently represented by the «Beam Pattern» function (Capon, 1969) 1043 Fig. 2. Histograms of station spacing (left) and relative orientation (right) for all the indepedent station pairs of the array. 1044 Gilberto Saccorotti, Bellina Di Lieto, Fabrizio Tronca, Costantino Fischione, Roberto Scarpa and Rosanna Muscente (3.1) where ⎢K⎪= 2 π/ λ is the wave number correspon- ding to wavelength λ, N is the number of stations, the x r i s are the station coordinates, and x r 0 is the vector position of the reference sensor. This func- tion describes the shape of the broad-band wave number spectrum for a vertically-incident plane wave. In practice, it represents the ’pulse re- sponse’ of the array in the wavenumber domain, and is particularly useful for investigating the ef- fects of spatial aliasing which manifest as spuri- ous peaks located at wave-numbers other than the origin. The UNDERSEIS wavenumber spectrum is shown in fig. 3. As a consequence of the ENE- WSW elongation of the array, the main spectral peak depicts a marked NNW-SSE trend, this im- plying low resolving capabilites once measuring waves propagating along that particular direction. Aliasing peaks appear at wavenumbers of abouts 20 cycles/km, corresponding to wave- lengths in the order of 300 m. On the other side, B N e 1 S i x x i N K i 0= $− −2 2 ^ ^h h/ the lowest bounds for body waves apparent veloc- ities correspond to the seismic velocities of the shallowest crust, here taken equal to 4 km/s (P- waves) and 2.3 km/s (S-waves). Consequently, the Nyquist frequencies for the two wavetypes are in the order of ≈ 13 Hz and ≈ 6 Hz, respectively. 3.2. Signal coherency All multichannel processing methods are based upon the common waveform model of the signal: in other words, it is assumed that the sig- nal recorded by each array element is the same except for a phase factor associated with the propagation across the deployment. In order to fullfil such assumption, it is crucial for the signal to maintain significant coherency throughout the different sensors. Using a number of P- and S- wave arrivals recorded at different SNRs, we cal- culated the array-averaged signal coherency as a function of frequency for both type of waves us- ing the formula (Saccorotti and Del Pezzo, 2000): (3.2) where f is the frequency, Np = N (N−1)/2 is the number of independent station pairs in the ar- ray, and czi,j is the Fisher’s Z-transform of the correlation coefficient among the band-pass filtered signals from stations i and j. Signal coherency is expected to vary from event to event as a function of a number of factors, including SNR and the frequency content of the signal and noise. However, we observed that in general P-waves maintain significant correlations over the 1-10 Hz frequency range, while S-wave are significantly coherent only over the 1-6 Hz frequency range. Combined with the Nyquist fre- quency limits for spatial aliasing mentioned above, these considerations provided constraints on the frequency bands to adopt for the analysis of earthquake wavefields described throughout the next sections. Similar correlation analyses were performed over different noise windows collected at day- and nigh-times. For the whole set of inter- station distances, significant noise coherency was only observed in association with marine micro- seisms, which typically manifest at frequencies C f a N cz i jtanh 1 , , p i j i j N !=^ eh o/ Fig. 3. Broad-band beam pattern spectrum (see eq. (3.1) in text) representing the response of the array to a monochromatic, vertically-inciding plane wave. Kx and Ky are the two components of the wavenumber vector with respect to a Cartesian system in which the x- and y-axes are oriented W-E and S-N, respectively. Performances of the UNDERground SEISmic array for the analysis of seismicity in Central Italy lower than ≈ 1 Hz (fig. 4). At higher frequencies, the correlation plots reported in fig. 4 fully sup- port the assumption of spatially-uncorrelated seis- mic noise required by the multichannel process- ing schemes presented in the following. 4. Methods of analysis We analysed UNDERSEIS data using up- to-date techniques aimed at fully exploiting the resolving capabilities of a multi-channel system toward measurement of the propagation and po- larization parameters of the incident wavefields. 4.1. Polarization For polarization analysis, we adopt the tech- nique proposed by Jurkevics (1988), which con- sists of a multichannel extension of the original covariance method of Kanasewich (1981). For a given signal window, the polarization parameters (azimuth and incidence angles, degree of lineari- ty) are derived from the eigenvectors and relation- ships among eigenvalues of an array-averaged co- variance matrix of the three-components of ground motion. In theory, the array-averaging procedure allows the contribution of noise (here assumed to be spatially uncorrelated) to be re- duced by a factor to the covariance estimates. For each window of analysis, signals from indi- vidual stations should be aligned according to the slowness measured for that window. In our situa- tion, however, alignment is not necessary consid- ering that the inter-station time delays are much smaller than the dominant period of the signal. 4.2. Frequency-slowness analysis For extracting the propagation parameters of UNDERSEIS data we use MUSIC, the MUlti- ple SIgnal Classification technique (Schmidt, N 1045 Fig. 4. Noise correlation as a function of frequency for different distance ranges. The correlation functions have been obtained from averaging the frequency-dependent correlation estimates associated with station pairs whose spacing d is included into the bounds reported at the top of each plot. 1046 Gilberto Saccorotti, Bellina Di Lieto, Fabrizio Tronca, Costantino Fischione, Roberto Scarpa and Rosanna Muscente 1986; Goldstein and Archuleta, 1987, 1991). Once compared to more classical methods (e.g., Beamforming or Capon’s (1969) High Resolu- tion), MUSIC has demonstrated superior per- formances for both the cases of low SNR signals and more plane waves simultaneously imping- ing on the array (e.g., Goldstein and Archuleta, 1987, 1991; Chiou and Bolt, 1993). For a given window of analysis, the MUSIC slowness spec- trum is expressed as (Chiou and Bolt, 1993) (4.1) where M is the number of signals impinging at the array, whose number is estimated applying Akaike’s information Criterium (AIC; Wang and Kaveh, 1985) to the ordered set of eigenval- ues of the spatial Cross-Spectral Matrix (CSM); Vi is the i-th eigenvector of the CSM, and A(S r )H is the Hermitian of a column vector containing the expected inter-station phase delays for a monochromatic plane-wave of frequency ω0 propagating with slowness S r . Note that the ele- Q S SA V1 1 H i i M 2 1 $ = − = ^ ^ h h/ ments of vector A correspond to the terms in the summation of eq. (3.1) above. Maximizing eq. (3.2) means finding those slowness vector(s) for which the array response has a maximum pro- jection onto the signal subspace, which is spanned by the i=1, ..., M eigenvectors Vi. Peaks in the Q(S r ) function above are thus associated with the horizontal slownesses of individual plane-wave components crossing the array. The amplitude of these peaks gives a measure of the extent to which the observed inter-station phase delays fit the model of a plane-wave propagating at the particular slowness S r , and is completely independe of the actual amplitude of the signal. 4.3. Uncertainties estimation We obtain an empirical estimate of the un- certainties associated with the measurement of propagation parameters by considering the range of horizontal slownesses and propagation azimuths associated with points of the slowness spectrum whose power is greater than 90% of Fig. 5. Sketch of the procedures used for retrieving uncertainties in the estimates of the propagation (left) and polarization (right) parameters. On the left, a sample slowness spectrum is displayed, with evidenced the con- tour line encompassing that region for which the power is greater than 90% of the main peak’s power. This re- gion is used to define errors in azimuth and ray parameter estimates. On the right plot, the shaded area marks the projection onto the horizontal plane of the polarization ellipsoid. El are the horizontal projections of the three eigenvectors; ϕ is the estimated polarization azimuths, whose uncertainties are defined through the angles ∆ϕ. Performances of the UNDERground SEISmic array for the analysis of seismicity in Central Italy 1047 Fig. 6a,b. a) Raw velocity seismograms and amplitude spectra for the P- and S-wave arrivals as recorded by the vertical-and rotated transverse-component, respectively, for a Md = 2.4 earthquake on August 26, 2003, at 08:43:56 GMT. The epicentral range is about 25 km. Seismograms and spectra are normalized to the largest-amplitude trace. Thin lines in the spectral plots refer to noise spectra calculated over 2-s-long windows preceeding the onset of the event. b) The same as in a), but for a Md = 3.3 earthquake on December 1, 2002, at 23:29:23 GMT. The epicentral range is about 124 km. the peak’s power (fig. 5, left). For the polariza- tion azimuth, angular uncertainties are derived from the ratio of the horizontal projection of the principal and secondary eigenvectors (fig. 5, right). This procedure represents a first attempt toward accounting for the large azimuthal un- certainties which arise for small incidence an- gles, even when the motion is very rectilinear. A similar procedure is followed for obtaining the uncertainties on the incidence angle. 5. Analysis of local and regional seismicity Figure 6a,b shows sample seismograms and amplitude spectra for the vertical-and transverse- a b 10481048 component of ground velocity associated with two earthquakes recorded at local and regional distances. For both cases, the figure also reports the amplitude spectra associated with noise win- dows preceeding the onset of the events. Despite the low magnitude of these earthquakes, seismo- grams depict an excellent Signal-to-Noise Ratio, and significant wave coherence is observed among the different array elements. In the follow- ing, we proceed by describing application of the techniques illustrated above to an heterogenous set of earthquakes as recorded by UNDERSEIS. The first step in calculating MUSIC slowness spectra is to obtain robust estimates of the nar- row-band spatial covariance of the signal. There- fore, we first band-pass-filter the array recordings over a narrow frequency band using a zero-phase- shift, 2-pole Butterworth filter. The filtered sig- nals are then Hilbert-transformed, and converted to Analytic Signals which eventually serve to ob- tain the complex-valued estimates of the spatial covariance. Eigenvalues and eigenvectors of the covariance matrix are derived using Singular Val- ue Decomposition. In order to account for the variable frequency content associated with differ- ent wavetypes, we conduct the analysis over con- secutive, 2-Hz-wide frequency bands spanning the 1-10 Hz and 1-6 Hz frequency range for P- and S-waves, respectively. Slowness spectra eval- uated over individual frequency bands are then Fig. 7a-e. a) Vertical-component seismograms for an event occurred at 10:32:30 GMT on August 31, 2003. b) Time interval used for P-wave slowness analysis. c) Multichannel signal Coherency (eq. (3.2) in text) over the frequency band selected for the analysis. The different lines refer to the five consecutive time position of the win- dow of analysis. d) The ordered set of normalized eigenvalues derived from diagonalization of the spatial co- variance matrix. e) Time-Frequency stacked slowness spectrum. Propagation azimuth φ and ray parameter ⎢P⎪ associated with the peak value are reported at the bottom left of the plot. The arrow points to the direction of wave propagation. a b c d e Gilberto Saccorotti, Bellina Di Lieto, Fabrizio Tronca, Costantino Fischione, Roberto Scarpa and Rosanna Muscente 1049 Performances of the UNDERground SEISmic array for the analysis of seismicity in Central Italy 1049 Performances of the UNDERground SEISmic array for the analysis of seismicity in Central Italy stacked to obtain a broad-band spectrum in which the contribution of body-waves is enhanced and the effects of the frequency-dependent peaks due to either spatial aliasing or surface-wave contam- inations is instead attenuated (Spudich and Op- penheimer, 1986). This procedure is repeated over five consecutive time windows encompass- ing the selected wave arrival, and a final slowness spectrum is derived from the stack of these time- varying estimates. We first apply slowness analy- sis to P-wave arrivals recorded at the vertical components. Then, we rotate the horizontal com- ponents into the transverse direction to the back- azimuth obtained from the P-wave, and evaluate S-wave slowness spectra over the rotated compo- nents. Examples of application of MUSIC slow- ness measurements are displayed in figs. 7a-e and 8a-e, which report P- and S-waves slow- ness spectra obtained for a Md = 2.4 earthquake recorded on August 31, 2003, at an epicentral distance of about 25 km. Figure 9a-f depicts the particle motion patterns for the P- and S-wave arrivals of the same event. In addition to the slowness spectra and co- herence plots associated with the selected time windows, figs. 7a-e and 8a-e also report the or- dered set of eigenvalues derived from diagonal- ization of the spatial covariance matrix. The ab- solute dominance of the main eigenvalue over the remaining ones indicates that a single signal is impinging on the array, in turn confirming the marked separation between the signal and noise subspaces previously postulated on the basis of noise correlation analyses. The performances of UNDERSEIS toward detection and location of local seismicity be- come even more evident for the case of low Fig. 8a-e. The same as in fig. 7a-e, but for the S-wave arrival measured on the transverse components. a b c ed 1050 Gilberto Saccorotti, Bellina Di Lieto, Fabrizio Tronca, Costantino Fischione, Roberto Scarpa and Rosanna Muscente SNR signals. An example is given by a Md = 2 earthquake which occurred on December 27, 2002, and located by the RSNC at coordinates 41.938°N, 14.007°E, ESE of the town of Sul- mona. This event was recorded by UNDER- SEIS with a very low amplitude, thus making impossible any visual estimate of the P-wave onset (fig. 10). However, the multichannel analyses provide consistent slowness estimates, which eventually permit a robust location of this event. For the sample earthquake reported in fig. 10, both the power of the stacked slowness spectrum and the multichannel coherency be- fore the signal’s onset are very low, thus confirming that the background noise is essen- tially uncorrelated among the different array el- ements. Based on this property, the power of the slowness spectra may be used as a powerful tool for automatic signal discrimination applied to real-time, routine analysis. The above procedures are extended to a fur- ther 25 earthquakes between November 2002 and September 2003. Magnitudes and epicentral distances of these events span the 2.3 ≤ M ≤ 3.5 and 20-140 km ranges, respectively (table I). Figure 11 illustrates the relationship among the measured P- and S-waves ray parameters and the epicentral range, here presented in terms of S-P delay time. For the two farest events (∆ ≈ 140 km), P-wave ray parameters range be- tween 0.1 s/km and 0.15 s/km, corresponding to apparent velocities spanning the 6.5-10 km/s range. Once accounting for the measurement a b c d e f Fig. 9a-f. a) Three-component, array-stacked waveforms for the same event displayed in figs. 7a-e and 8a-e. b, c) Time windows respectively used for P- and S-wave polarization estimates. d, e) P-wave particle motion or- bits over the horizontal and radial-vertical plane, respectively. f) S-wave particle motion over the horizontal plane. P-wave polarization azimuth and incidence angles obtained from the covariance analysis are 35° and 12°, respectively. The S-wave polarization azimuth is 110°. 1051 Performances of the UNDERground SEISmic array for the analysis of seismicity in Central Italy Fig. 10. Time behaviour of the propagation parameters for a Md = 2 event at 12:42 GMT on December 27, 2002, at an epicentral distance of about 80 km. From top to bottom, sample vertical-component seismogram, propaga- tion azimuth, ray parameter, multichannel coherency and power of the frequency-stacked slowness spectra. The P-wave arrival is marked by the sudden increase in both the slowness spectral power and coherency values about 4 s into the event. Table I. List of analysed earthquakes. Event ID Year Month Day OT Ts-Tp (s) Md RSNC Lat (°) RSNC Long (°) 8 2002 11 20 02:31:51 3.1 2.4 42.480 13.240 9 2002 12 1 23:29:23 16.5 3.3 41.660 14.790 10 2002 12 2 20:52:59 18 2.7 41.670 14.860 12 2002 12 23 04:17:40 5.5 2.3 42.688 13.027 14 2002 12 24 07:14:05 3.3 2.7 42.490 13.240 15 2002 12 27 12:42:37 7.3 3 41.940 14.010 1 2003 1 1 20:23:40 7.5 2.6 41.900 13.690 3 2003 1 5 12:41:22 2.7 3 42.500 13.300 4 2003 4 7 15:01:15 2.8 2.5 42.320 13.450 5 2003 4 10 05:50:34 7.5 2.5 41.920 13.870 6 2003 4 11 11:49:10 3 2.5 42.510 13.270 7 2003 5 2 20:14:09 2.8 2.6 42.450 13.260 181 2003 6 30 22:32:20 3.5 2.5 42.300 13.240 186 2003 7 5 05:07:43 3.3 3.1 42.310 13.270 187 2003 7 6 09:13:40 3.4 2.2 42.310 13.230 Fig. 11. P- and S-waves ray parameters as a function of S-P delay times for all the events listed in table I. 1052 Table I (continued). Event ID Year Month Day OT Ts-Tp (s) Md RSNC Lat (°) RSNC Long (°) 188 2003 7 7 12:37:12 3.8 2.2 42.340 13.230 192 2003 7 11 02:05:33 3.9 2.2 42.290 13.270 197 2003 7 16 05:09:26 3.2 2.6 42.330 13.280 238 2003 8 26 08:43:56 3.4 2.4 42.300 13.240 243 2003 8 31 10:33:00 3.6 2.4 42.370 13.310 239 2003 8 27 17:18:13 3.6 2.4 42.520 13.27 244 2003 9 1 00:07:40 3.9 2.4 42.310 13.260 251 2003 9 8 23:03:50 3 2.7 42.290 13.260 261 2003 9 18 21:00:40 3.8 2.4 42.300 13.230 267 2003 9 24 09:03:00 3.2 2.8 42.480 13.230 270 2003 9 27 17:35:25 4.35 3.5 42.220 13.250 errors, these values are compatible with those expected for Pn waves. Events located at small- er epicentral distances depict however similar, or even greater apparent velocities, indicating the steep incidence of Pg waves from sources whose depth is comparable to the epicentral distance (Braun et al., 2004). As previously mentioned by Biondi and Kostov (1989), ray parameter estimates are much more unstable than the measurement of azimuths. The com- Gilberto Saccorotti, Bellina Di Lieto, Fabrizio Tronca, Costantino Fischione, Roberto Scarpa and Rosanna Muscente 1053 Performances of the UNDERground SEISmic array for the analysis of seismicity in Central Italy parison of P- and S-wave ray parameters shows a scattered pattern; however, the estimate of their mutual relationships provides a meaning- ful value (1.78 ± 0.24) for the P-to-S velocity ra- tio (fig. 12). Within the error bounds of individ- ual measurements, a general consistency is also observed among the P- and S-waves wave prop- agation azimuths (fig. 13), as well as among the P-waves polarization and propagation azimuths (fig. 14). The large uncertainties of slowness es- timates observed in figs. 11 and 12 for both P- and S-waves result from the large width of the main peak in the array response pattern, which is in turn controlled by the array aperture. Thus, major resolution improvements are expected following the enlargement of the array present- ly under completion. 6. Location At the beginning, we had hoped to obtain consistent estimates on the locations of UN- DERSEIS data from back-tracing the seismic wave-vectors evaluated from frequency-slow- ness analysis (e.g., Almendros et al., 2001). However, due to the large errors associated with ray parameter measurements, and the intrinsic Fig. 12. Comparison of the ray parameters derived from slowness analyses of P- and S-waves. Fig. 13. Comparison of the propagation azimuths derived from slowness analyses of P- and S-waves. Width and height of the boxes correspond to the un- certainty on slowness estimates. Fig. 14. Comparison of the propagation and polar- ization azimuths derived from analyses of P-waves arrivals. Width and height of the boxes correspond to the uncertainty on slowness and polarization esti- mates. 1054 Gilberto Saccorotti, Bellina Di Lieto, Fabrizio Tronca, Costantino Fischione, Roberto Scarpa and Rosanna Muscente uncertainties in wave types for events located at ranges comparable to the cross-over distance, this procedure didn’t provide satisfactory re- sults. We thus applied a much simpler proce- dure, deriving epicentral locations from propa- gation azimuths and using a factor γ of 7.5 to convert S-P time delays to epicentral ranges. Error bounds on such locations are convenient- ly derived from the uncertainties in azimuths estimates, and considering a variation on γ of ± 0.5. Such uncertainty accounts for errors in S- P time delays estimates, and for the fact that the same S-P delay time may be associated with events having different depths and hence differ- ent epicentral distance. For all the cases, the epicenters given by the CNSN catalogue are in- cluded into the confidence bounds of our esti- mates. Figure 15 reports the discrepancies among UNDERSEIS and CNSN locations as a function of epicentral distance (top) and ray pa- rameter (bottom) for both P- and S-waves. We observe that locations from S-waves depict dis- crepancies which rarely exceed 10 km, and are much more consistent than those retrieved from P-wave data. This is quite surprising: one would expect slowness analysis on the S-wave to provide less robust estimates, caused by path and recording site effects and by inaccurate sensor orientation. On the other side, however, the amplitude of S-waves is generally larger than that of P-waves, thus implying significant SNR improvements. Moreover, we have to con- sider that our slowness spectra were calculated over a cartesian slowness grid of regularly- spaced grid nodes: therefore, slowness meas- urements become progressively less accurate as one moves toward the origin of the grid, i.e. as the apparent velocity increases. Under such conditions, loss of precision is also induced by round-off errors once deriving very small inter- station phase shifts from the eigen-structure of the spatial cross-covariance matrix. Examina- tion of the relationships among location errors and ray parameters shown in fig. 15 confirms these two latter observations. The largest loca- tion errors are in fact associated with P-wave propagating at low ray parameters, thus sug- gesting that wave amplitude and velocity are the main factors controlling the precision of slowness measurements. Keeping in mind these limitations, the above location results may be considered rather satisfactorily, thus indicating that future opera- tions of UNDERSEIS will significantly help to Fig. 15. Location discrepancies as a function of epicentral distance (top) and ray parameter (bottom) from both P- and S-waves data. 1055 Performances of the UNDERground SEISmic array for the analysis of seismicity in Central Italy better delineate the distribution of micro-seis- micity in the Central Apennines. 7. Discussion and conclusions 7.1. Limitations of the array The two key factors controlling the accura- cy of array-based location procedures are: 1) the precision of slowness measurements, and 2) the ability to correctly predict ray-paths through the region encompassing the array and the epi- center. As for the first factor, in the previous section we observed that the largest location er- rors were associated with steeply-inciding P- waves. This is in agreement with theoretical prediction of slowness uncertainties (e.g., Gold- stein and Archuleta, 1991), in which for a given array configuration the robustness of slowness estimates are mostly controlled by SNR and ap- parent velocity. Thus, slowness measurements derived from the faster-propagating, lower-am- plitude P-waves are expected to be less accurate than those retrieved from the corresponding S- waves arrivals. At the time of revising this pa- per (early January 2006), 6 additional elements were added to UNDERSEIS, which may now count on a total aperture of about 800 m. Major technological upgrades are also under de- velopment in order to allow for sampling rates higher than actually available. From such im- provements, we expect a significant increase in measurement precision: on one side, the larger array aperture will improve the resolving capa- bilities toward fast waves, while on the other the higher number of available stations will re- duce the contribution of noise, whose effects are seen on the measurement accuracy scale as . Regarding the second factor mentioned above, the accurate prediction of ray paths would require knowledge of the velocity struc- ture down to scale lengths in the order of a few hundreds of meters, a task which is obviously beyond the purpose of this study. However, the currently missing assessment of source depths may be attained using gradient velocity models. Such procedure, although producing biased es- timates for individual event locations, would preserve the alignments depicted by clusters of N1 adjacent sources; final constraints on absolute locations of these clusters will eventually be gained by joint inversion of array wave-vectors and CNSN arrival-time data associated with the most energetic events. 7.2. Future perspectives In the present study, we were mostly inter- ested in evaluating the location performances of the array: therefore, we selected and analyzed only a subset of the array catalog for which CN- SN locations were available. The present efforts are mostly aimed at developing procedures for the automatic detection of earthquake signals, and a prototype on-line processor based on the algorithms presented above is already under testing. Following these efforts, the most imme- diate goal of UNDERSEIS is to improve the completeness of the seismic catalogue for the Central Apennines. At present, the detection magnitude threshold of the CNSN varies be- tween 1.5 and 2.6 throughout the national terri- tory, and is about 2.2-2.4 for the Central Apen- nines (Cattaneo et al., 2002). However, from the arguments presented above, we do expect UNDERSEIS to contribute significantly to the detection and location of micro-earthquakes un- der poor SNR conditions, thus opening the way to the precise definition of seismicity patterns in the surroundings of the array. In this work, we selected the relevant param- eters for slowness measurements (frequency band and length of the time windows) based on sample estimates of the signal’s frequency content and coherency. These properties, however, vary sig- nificantly from event to event. Any automatic ar- ray processing scheme should therefore include the possibility of performing the analysis only over those time-frequency windows for which the signal has maximum energy and coherency (Cornou et al., 2003), in turn adjusting the length of the window of analysis to the frequency band being analysed. This process could be made faster without loss of resolution by adopting so- phisticated Time-Frequency decompositions of the signal for the calculation of the frequency-de- pendent spatial covariance matrices. For this pur- pose, an extension of the polarization and slow- 1056 Gilberto Saccorotti, Bellina Di Lieto, Fabrizio Tronca, Costantino Fischione, Roberto Scarpa and Rosanna Muscente ness measurements to the wavelet-transformed signals (e.g., Bear and Pavlis, 1997) is currently under development. Such improvements would significantly simplify the computational tasks, thus permitting an easier implementation of auto- matic, on-line detection procedures. Once sys- tematically analysing UNDERSEIS data, a fur- ther factor which should be taken into account re- gards the quantitative definition of the properties of seismic noise, and their variation over time. Al- though we observed that the noise was spatially uncorrelated, there may be periods of particular meterological condition or intense human activi- ty for which coherent seismic signals may propa- gate across the array. Future studies must there- fore address the problem of a systematic analysis of the characteristics of seismic noise throughout extended time spans, in order to fully elucidate the regions of the slowness space which may be affected by the propagation of coherent signals other than those associated with earthquake sources. Acknowledgements Editor A. Rovelli and reviewers T. Braun and G. Milana provided many thoughtful comments which greatly helped to improving the quality of the manuscript. The financial support from Con- sorzio Gran Sasso and Laboratori Nazionali del Gran Sasso made possibile the realization of the seismic array. We acknowledge E. Bellotti and E. Boschi for support and stimulation. This work has been conducted in the framework of the pro- gram agreement between INGV-Osservatorio Vesuviano and Istituto Nazionale di Fisica Nu- cleare. We also acknowledge financial support from MIUR-PRIN2005 project «Analisi e mod- ellistica dei processi sismici e deformativi nel- l’Appennino Centrale». REFERENCES ALMENDROS, J., B.A. CHOUET and P.B. DAWSON (2001): Spatial extent of a hydrothermal system at Kilauea vol- cano, Hawaii, determined from array analyses of shal- low long-period seismicity, 1. Method, J. Geophys. Res., 106, 13565-13580. AMORUSO, A., L. CRESCENTINI and R. SCARPA (1998): Inver- sion of source paraeters from near-and far-field observa- tions: an application to the 1915 Fucino earthquake, Central Appennines, Italy, J. Geophys. Res., 103, 29989- 29999. AMORUSO, A., L. CRESCENTINI, A. MORELLI and R. SCARPA (2002): Slow rupture of an aseismic fault in a seismo- genic region of Central italy, Geophys. Res. Lett., 29, 2.219, doi: 10.129/2002GL016027. BEAR, L.K. and G.L. PAVLIS (1997): Estimation of slowness vectors and their uncertainties using multi-wavelet seismic array processing, Bull. Seismol. Soc. Am., 87, 755-769 BIONDI, B.L. and C. KOSTOV (1989): High-resolution veloc- ity spectra using eigenstructure methods, Geophysics, 54, 832-842. BRATT, S.R. and T.C. BRACE (1985): Locating events with a sparse network of regional arrays, Bull. Seismol. Soc. Am., 78, 780-798. BRAUN, T., J. SCHWEITZER, R.M. AZZARA, D. PICCININI, M. COCCO and E. BOSCHI (2004): Results from the tempo- rary installation of a small aperture seismic array in the Central Apennines and its merits for local event detec- tion and location capabilities, Ann. Geophysics, 47 (5), 1557-1568. CAPON, J. (1969): High-resolution frequency-wavenumber spectrum analysis, Proc. IEEE, 57, 1408-1418. CATTANEO, M., P. AUGLIERA and M. DE MARTIN (2002): Seismic noise measurements for the Italian seismic network, Geophys. Res. Abs., SCA/B-3. CHIARABBA, C., A. AMATO, P. AUGLIERA, S. BAGH, M. CATTA- NEO, L. CHIARALUCE, P. DE GORI, P. DI BARTOLOMEO, A. GOVONI, A. MICHELINI, M. MORETTI, D. PICCININI and M. ROMANELLI (2004): The MAFI project: Mapping Active Faults in Italy using microseismicity data, in Fall Meet- ing 2004, American Geophysical Union, S52A-01. CHIARABBA, C., L. JOVANE and R. DI STEFANO (2005): A new view of Italian seismicity using 20 years of instru- mental recordings, Tectonophysics, 395, 251-268. CHIOU, S.J. and B.A. BOLT (1993): Seismic wave slowness- vector estimation from broad-band array data, Geo- phys. J. Int., 114, 234-248. CORNOU, C., P.Y. BARD and M. DIETRICH (2003): Contribu- tion of dense array analysis to the identification and quantification of basin-edge-induced waves, Part I. Methodology, Bull. Seismol. Soc. Am., 93, 2604-2623. CRESCENTINI, L., A. AMORUSO and R. SCARPA (1999): Con- straints on slow earthquakes dynamics from a swarm in Central Italy, Science, 286, 2132-2134. DE LUCA, G., E. DEL PEZZO, F. DI LUCCIO, L. MARGHERITI, G. MILANA and R. SCARPA (1998): Site response study in Abruzzo (Central Italy): underground array versus surface stations, J. Seismol., 2, 213-226. DE LUCA, G., R. SCARPA, L. FILIPPI, A. GORINI, S. MARCUC- CI, P. MARSAN, G. MILANA and E. ZAMBONELLI (2000): A detailed analysis of two seismic sequences in Abruz- zo, Central Apennines, Italy, J. Seismol., 4, 1-21. GOLDSTEIN, P. and R.J. ARCHULETA (1987): Array analysis of seismic signals, Geophys. Res. Lett., 14, 13-16. GOLDSTEIN, P. and R.J. ARCHULETA (1991): Deterministic frequency-wavenumber methods and direct measure- ments of rupture propagation during earthquakes using a dense array: theory and methods, J. Geophys. Res., 96, 6173-6185. 1057 Performances of the UNDERground SEISmic array for the analysis of seismicity in Central Italy JURKEVICS, A. (1988): Polarization analysis of three-com- ponent array data, Bull. Seismol. Soc. Am., 78, 1725- 1743. KANASEWICH, E.R. (1981): Time Sequence Analysis in Geo- physics (University of Alberta Press, Edmonton), 1-532. ROST, S. and C. THOMAS (2002): Array seismology: meth- ods and applications, Rev. Geophys., 40, 1-27. SACCOROTTI, G. and E. DEL PEZZO (2000): A probabilistic approach to the inversion of data from a seismic array and its application to volcanic signals, Geophys. J. Int., 143, 249-261. SCARPA, R., R. MUSCENTE, F. TRONCA, C. FISCHIONE, P. ROTELLA, M. ABRIL, G. ALGUACIL, W. DE CESARE and M. MARTINI (2004): UNDERSEIS: the UNDERground SEISmic array, Seismol. Res. Lett., 75, 493-504. SCHMIDT, R.O. (1986): Multiple emitter location and signal parameter estimation, IEEE Trans. Antennas Propag., 34, 276-280. SCHWEITZER, J., J. FYEN, S. MYKKELTVEIT and T. KVAERNA (2002): Chapter 9: Seismic Arrays, in IASPEI New Manual of Seismological Observatory Practice (Geo- ForschungsZentrum Potsdam), edited by P. BORMANN, vol. 1, pp. 46. SPUDICH, P. and D. OPPENHEIMER (1986): Dense seismo- graph array observations of earthquake rupture dynam- ics, in Earthquake Source Mechanics, edited by S. DAS, J. BOATWRIGHT and C. SCHOLZ, Am. Geophys. Un., Geophys. Monogr. Ser.., 37, 285-296. WANG, H. and M. KAVEH (1985): Coherent signal-subspace processing for the detection and estimation of angles of arrival of multiple wide-band sources, IEEE Trans. AS- SP, 33, 823-831. (received January 19, 2006; accepted October 20, 2006)