Centro di Ricerche Sismologiche


eqkeqk                                                         

              "TRUE WOOD-ANDERSON DIGITAL SYSTEM"

                                                    

         

       IntroduzioneIntroduction

       La nuova leva otticaThe new optical leverage

       Il PSDThe PSD

       AcquisizioneRecording                                                                                                                                               

       TaratureCalibration

       Ingrandimenti staticiStatic magnification

       ConclusioniConclusions

        Riferimenti bibliograficiReferences


jp1


INTRODUZIONEINTRODUCTION

           Nel corso del 2003 sono stati completamente riattivati i due sismografi Wood-Anderson in dotazione alla stazione sismografica di Trieste  operanti dal 1972 fino al  1987. L'originale sorgente luminosa, costituita da una lampadina ad incandescenza, è stata sostituita da un laser a luce rossa visibile e la carta fotosensibile da un "position sensing detector" (PSD). Il PSD  è un dispositivo semiconduttore che fornisce in uscita una tensione proporzionale alla posizione del punto luminoso sulla sua superficie. E' quindi un trasduttore analogico. La sostituzione del vecchio sistema di registrazione è quindi diretta ed immediata. Dopo essere stati filtrati (anti alias) i segnali vengono convertiti in digitale e memorizzati su PC.  Nel corso del 2004 ai due strumenti è stato affiancato un sismometro broadband tre componenti con lo scopo di verificarne l'ingrandimento statico.The original Wood-Anderson instruments in use from 1972 until 1987 in the seismological Trieste department were again put into operation during 2003 was modernized passing the recording system from the analogic to the digital format. The original light source, a bulb lamp, was replaced by a red visible laser source and the photosensitive paper by a "position sensing detector" (PSD). The PSD is a semiconductor device giving as output a voltage  proportional to the spot light position on its surface. Then it is an alogue device. After antialiasing filtering, the data are converted in digital format and recorded in a hard disk. In 2004 a broadband seismometer was placed side by side to the instruments to check their static magnification.
 
 

LA NUOVA LEVA OTTICATHE NEW OPTICAL LEVERAGE jp2

             Ai due strumenti originali è stato tolto il secondo specchio cilindrico e le due lenti poste nella finestra ottica. Togliendo lo specchio cilindrico con il suo supporto si è ottenuta una piccola finestra laterale attraverso la quale si fa passare il raggio laser. Così facendo restano inalterati gli angoli operativi dello specchio mobile e quindi le procedure di taratura dell'apparecchiatura. L'unica differenza con l'originale è che la leva ottica originale, subendo i raggi luminosi due riflessioni, era moltiplicata per quattro e la nuova leva, subendo il raggio una sola riflessione, viene invece moltiplicata per due.The second cylindrical mirror and the two lenses placed in the optical window were removed from the two original instruments. The removal of the cylindrical mirror with its support has produced a small side window through which the laser ray hits the mirror. In this way, operating angles of the mobile mirror and consequently the equipment calibration remain unchanged.  Changing the numbers of the light reflections, now only one, effects the optical lever to change from four time the fisical length of the light path  to only two times.
 
 
 

IL PSDTHE PSD
 

psd
 
 

          Il sensore utilizzato è l' 1L20 della Sitek. E' un sensore monodimensionale con un'area attiva di 20x3 mm sensibile alle radiazioni visibili. La linearità tipica è dello 0.1% entro l'80% della lunghezza della superficie utile. La risoluzione può arrivare ad essere migliore di una parte su un milione e dipende dalle soluzioni ottiche, elettriche e meccaniche adottate. Ha due anodi (Y1 e Y2) ed un catodo (bias). E' sensibile alla posizione del centroide dello spot luminoso incidente sulla sua superficie attiva. La posizione Y del centroide viene calcolata con:The sensor used is the 1L20 of the Sitek. It is a 1D sensor with an active area of 20x3 mm sensitive to visible radiation. The typical linearity is 0.1% within the 80% of the length of the surface. The resolution may be better than one part per million and depends on the optical, electrical and mechanical solutions adopted. It has two anodes (Y1 and Y2) and a cathode (bias). It's sensitive to the centroid position of the spot light incident on its active surface. The centroid position along the Y is calculated by:
 
 

posmath

            Nel 1996 l'INGV aveva realizzato un sistema simile ma basato su di un sensore CCD che, essendo essenzialmente un trasduttore digitale, richiede elaborazioni aggiuntive rispetto a quelle richieste dal PSD, sia per l'interpolazione della posizione del segnale ottico sia per l'attuazione di filtri anti-alias ( Romeo, 1996)In 1996 the INGV realized a similar system, based on a CCD sensor. The CCD is essentially a digital transducer with additional processing for the interpolation of the optical signal position and more complexity for implementing anti-alias filters (Romeo, 1996) compared to the PSD..

ACQUISIZIONE DEI SEGNALI REGISTRATI DAI DUE TS-220RECORDING
 

       Al segnale in uscita viene applicato un filtro anti alias a 6 poli con 40 Hz di frequenza di taglio. L'acquisitore è a 16 bit con frequenza di campionamento di 100 Hz.  Il salvataggio dei dati avviene tramite trigger STA/LTA. Vengono memorizzati solo gli eventi rilevati dal trigger. La temporizzazione è realizzata tramite un sistema GPS.At the output signal is applied a 6-poles anti alias filter with a 40 Hz frequency cutting. The 16-bit recorder has a 100 Hz sample rate. Data storing is performed via trigger STA / LTA and it is possible a continuous recording. Only triggered data are currently stored.  Timing is achieved through a GPS system.
 

TARATURECALIBRATION
 

       I due strumenti Lhener-Griffith TS-220 sono orientati N-S ed E-O. Sono stati smontati, puliti e completamente ricalibrati.
Nella tabella seguente sono riportate le caratteristiche delle due componenti.
The two instruments Lhener Griffith-TS-220 are oriented NS and EW. They were disassembled, cleaned and completely recalibrated. The following table shows the two components characteristics.
 
 
 

componenteOrientation periodo (s)Period (s) smorzamentoDamping
N-S 0.792 0.787
E-O 0.796 0.818

 

             Originariamente quando l'indice posto sulla base degli strumenti era su uno dei riferimenti laterali riportati sulla base stessa si doveva leggere sulla carta uno spostamento di 50 mm. Ora basta dividere 50 mm per il numero di counts risultanti per stabilire la sensibilità  ed impostarla nel codice di acquisizione dati. In uscita avremo direttamente l'ampiezza in mm delle tracce equivalente a quella che si sarebbe avuta sulla carta fotosensibile.Originally, when the index placed on the basis of the instruments was on one of the marks you had to read the displacement of 50 mm on the paper. Now you need just to divide 50 mm for the resulting number of counts to determine the sensitivity of the instruments and set it in the data acquisition code. In this way in output we will have directly the width in mm of the tracks, equivalent to that which we would have on the photographic paper.
 

  INGRANDIMENTI STATICISTATIC MAGNIFICATION
 

             In letteratura si trovano lavori che hanno tentato di stabilire l'ingrandimento statico del TS-220 dalle sue caratteristiche meccaniche (Uhrhammer et al., 1990)  . Occorre comunque considerare che la costruzione meccanica e la filosofia generale di applicazione dello strumento rendono molto ardua un' analisi di questo tipo.  La valutazione dell'ingrandimento statico dipende soprattutto dall'esatta collocazione del filo di sospensione (diametro 0.02 mm) su di un cilindretto di 2 mm di diametro sul quale è fissato lo specchio mobile. Il filo deve essere esattamente a  22.5°  rispetto al piano dello specchio(Anderson et al., 1925). In fase di calibrazione l'esatto valore di quest'angolo, difficilmente verificabile a causa delle dimensioni in gioco, è fondamentale  nella valutazione della componente dell'accelerazione di gravità applicata. Un errore di, per esempio, 0.1 mm nel fissaggio del filo equivarrebbe ad un errore nell'angolo di circa  6° e ad un errore sulla valutazione dell'ingrandimento di circa  140. Altra causa di errore è la difficoltà nel livellare la base dello strumento in modo da disporre lo specchio mobile a 22.5° rispetto alla superficie esterna del sismometro. Tutto questo porta ad un notevole spread di valori quando si vanno a testare diversi strumenti: Gutenberg nel 1957 ottenne, utilizzando un tavolo vibrante, un valore di 2800 +/- 500.  Il Wood-Anderson reale evita la gestione di questi problemi: lo strumento viene tarato per avere una traccia ampia 50 mm quando l'indice della base è posto su uno dei riferimenti laterali. Il problema si pone quando si vuole simulare lo strumento su registrazioni  effettuate da altri sismometri (es.: broad-band).In literature there are works that have attempted to establish the static magnification of TS-220 from its mechanical characteristics (Uhrhammer et al., 1990). However, the construction mechanics and the general philosophy of the instrument application make it very difficult. The static magnification assessment depends mainly on the exact location of a cylinder of 2 mm in diameter, where the mobile mirror is attached, suspended on a wire of 0.02 mm in diameter. The wire must have exactly an angle of 22.5° with the plan of the mirror (Anderson et al., 1925). The correct calibration of this angle, really difficultly verificable because of the sizes involved, is crucial in assessing the gravity acceleration component applied. An error, for example, of 0.1 mm in the wire fixing procedure means an angle error of about 6° and a magnification evaluation error of about 140. Another cause of errors is the difficulty in leveling the instrument base to put the mobile mirror mobile at a 22.5° angle with the outer surface of the seismometer. All these problems lead to a significant spread of values when different instruments must be tested : in 1957 Gutenberg obtained, using a vibrating table, a value of 2800 +/- 500. The Wood-Anderson usage avoids the management of these problems: the instrument is calibrated to have a track of 50 mm when the base index is placed on one of the side references. Problems arise simulating the instrument with recordings made by other seismometers (eg.: broadband).

            Per l'analisi degli ingrandimenti statici reali dei due wood-anderson, nel corso del 2004 è stato affiancato ad essi un sismografo broadband   Guralp 40-T con periodo massimo di 60 s. L'acquisizione, in continua, avviene tramite digitalizzatore Guralp 24 bit integrato nello strumento. Sui 1252 eventi registrati fino al dicembre 2013   è stato effettuato un confronto tra le ampiezze registrate dagli strumenti originali e quelle ottenute simulandole sui segnali broadband. L'ingrandimento statico  G è definito come:For checking the real static magnification of the two Wood-Anderson instruments, during 2004, a Guralp 40-T broadband seismometer with period extended to 60 s was placed  near our instruments,  acquisition accomplished by an embedded 24 digits digitizer. Until 2013, 1252 events have been recorded and Sandron et al. (2015) compared the amplitude of the real wood-anderson instruments with that ones obtained by their simulation on broadband. The static magnification G is:

     G = A / a  

dove where A è l'ampiezza di picco della traccia sul sismogramma e a è lo spostamento del terreno. Quindi se Gw e Grispettivamente gli ingrandimenti del
wood-anderson e della sua simulazione sul broadband e Aw e Ab le rispettive ampiezze, ponendo Gb uguale a 2800:
A is the trace peak amplitude and is the ground motion. If  Gw and Gare respectivetly the magnification of wood-anderson and its simulation over broadband signal and Aw and Atheir peak amplitude, setting Gb = 2800 will be:

    Gw = 2800 Aw/Ab G=A/a
 
Nella seguente figura in e b si possono vedere i risultati ottenuti tramite media pesata su una finestra mobile di 50 eventi.  Analizzando l'andamento delle funzioni inaspettatamente si nota che  la risposta risulta  non essere lineare ma seguire una sorta di legge di potenza funzione dell'ampiezza dei segnali (Sandron et al., 2015). Dal grafico N-S  si può constatare che l'ingrandimento statico risulta avvicinarsi a 2800 per ampiezze inferiori  a 0,05 mm, e tendere asintoticamente a 2080 per le ampiezze superiori a 1 mm mentre per la componente E-W la curva risulta spostata verso l'alto di circa 200, tendente quindi a 3000 per le basse ampiezze e 2280 per le ampiezze maggiori. I risultati riguardanti la componente N-S sono più affidabili di quella E-W. La stima dell'ingrandimento di quest'ultima  potrebbe non essere corretta poiché lo specchietto, che era parzialmente staccato dalla massa,  è stato riattaccato cercando di conservare i 22,5° (vedi sopra) cosa per nulla agevole viste le dimensioni in gioco.The following figures show the results achieved by weighted average over a 50 events sliding window.  We can see that the response is not linear but follows a power law, function of the amplitude of the trace.  In the N-S graph the static magnification asympototically approaches 2800 for lower amplitudes and  2080 for the higher ones. In the E-W graph the response is shifted upward of about 200.  Results of the N-S component are more reliable because in the E-W the mirror was partially detached from the wire and it was repaired by us trying to save the 22.5° (see above).

                                                                                              

                                  

                


          

            Fig. 1 -  Ingrandimenti statici e magnitudo corrispondenti dei wood-anderson reali in funzione di quelli simulati su registrazioni broadband delle componenti nord-sud (a, c) e est-ovest (b, d). Da Sandron et al, 2015Fig. 1 - Static magnification of the real wood-anderson instruments and magnitude difference versus their amplitude (simulated) on broadband data. N-S component on left (a, c) and on right the E-W (b, d). See Sandron et al, 2015
 

CONCLUSIONICONCLUSIONS
 

          Alla luce di quanto emerso dalla comparazione delle ampiezze dei segnali reali WA con quelle dei segnali simulati dalle registrazioni broadband appare evidente che nell'utilizzo degli strumenti originali si ha un'incertezza nel calcolo della magnitudo pari a circa +/- 0.13 gradi Richter dovuta alla dipendenza dell'ingrandimento statico dall' ampiezza del segnale. Questo ovviamente non si verifica simulando gli strumenti sulle registrazioni broadband, in questo caso però resta incerta e del tutto arbitraria  la scelta dell'ingrandimento statico da inserire nella funzione di trasferimento.Looking at our results about the comparison of real Wood-Anderson data with their simulation over broadband signals, magnitudes calculated by the original instruments have an uncertainty of about ± 0.129 ±0.129 Richter degrees dues to the dipendence of their static magnification on amplitude of the peak traces. Obvioulsy such a behavior does not  appear in the simulation over broadband data, but in this case the uncertainty of which magnification (2800 or 2080) to use remains. 


                                                                                         DATIDATA
 

       Gli strumenti sono installati presso la sede dell'OGS in Borgo Grotta Gigante (Trieste) (coord: 45.709° N, 13.764° E, alt.: 330 slm).The two instruments are located at the OGS in Borgo Grotta Gigante (Trieste) (coord: 45.709° N, 13.764° E, alt.: 330 m).
  La componente N-S è attiva dal 31/10/2002 e quella E-O dal 13/12/2002.
         I seguenti links danno accesso ai cataloghi dei terremoti  registrati dai wood-anderson e dal broadband posto loro accanto
  2004-2013 2014 2015.         

         Magnitudo dei terremoti registrati della sequenza sismica successiva alla scossa principale del terremoto di Amatrice (RI) del 24/8/2016  link
         Sono riportate sia le magnitudo locali calcolate con i wood-anderson (MAW) sia quelle ottenute simulando gli strumenti sul broadband accanto ad essi con ingrandimento 2800 (ML2800).
The N-S component is working from 31/10/2002 and the E-W one from 13/12/2002. The broadband seismometer is operative from October 2004. In the following links the catalogues of the  earthquakes recorded through them 2004-2013  2014 2015.  
  
           The following link shows the magnitude of the earthquakes after the 08/24/2016  Amatrice mainshock.  Data  recorded by the real woodanderson instruments  (MAW column) and simulated  on the broadband recordings with magnification 2800 (ML2800 column), are reported.   link

 


  RIFERIMENTI BIBLIOGRAFICIREFERENCES
 

- Anderson J. A. and Wood H. O., 1925: Description and Theory of the Torsion Seismometer; BSSA, vol. XV.

- Finetti L., Morelli C., 1972: Earthquake magnitude determination for Trieste WWSS station; Bollettino di Geofisica Teorica ed Applicata, Vol. XIV,   n. 53-54.

- Manuale del sismometro Lehner-Griffith mod. TS-220.

- Richter C. F., 1935: An Instrumental Earthquake Magnitude Scale; BSSA, 25 pp. 1-32.

- Richter C. F., 1958: Elementary Seismology; W. H. Freeman and Company; S. Francisco.

- Uhrhammer R. A. and Collins E. R., 1990: Synthesis of Wood-Anderson Seismograms From Broadband Digital Records; BSSA, vol. 80, pp. 702-716.

- Uhrhammer R. A., Loper S. J. and Romanowicz B., 1996: Determination of Local Magnitude Using BDSN Broadband Records; BSSA, Vol. 86, n. 5,      pp. 1314-1330.

- Romeo G., 1996: Wood Anderson digitale; rapporto interno INGV.
 
  Sandron D., Gentile G. F., Gentili S., Saraò A., Rebez A., Slejko D., 2015: The Wood-Anderson of Trieste (Northeast Italy); One of the Last Operating Torsion Seismometers; Seismological Research Letters, vol. 86.