The Campanian Earthquake Early Warning Project by Aldo Zollo and Claudio Satriano RISCC-Lab, Napoli, Italy www.rissclab.unina.it
Earthquake Mitigation earthquake forecasting seismic risk enigineers/insurers seismic hazard long-term scientists earthquake prediction short-term seismic protection strategies decision makers long-term short-term redrawn from Snieder and van Eck, 1997
Real-Time Systems: Operating Principle Epicenter Seismic waves (~ 3.5 Km/s) Seismic Network Information (speed of light) Target
Real-Time Systems: Operating Principle Epicenter Seismic waves (~ 3.5 Km/s) Seismic Network Information (speed of light) Target T0 time
Real-Time Systems: Operating Principle Epicenter Seismic waves (~ 3.5 Km/s) Seismic Network Information (speed of light) Target P-Arrival T0 TP time
Real-Time Systems: Operating Principle Seismic waves (~ 3.5 Km/s) Information (speed of light) Epicenter Target Seismic Network P-Arrival T0 TP TS time
Real-Time Systems: Operating Principle Epicenter Seismic waves (~ 3.5 Km/s) Seismic Network Information (speed of light) Target P-Arrival Telemetry and computing T0 TP TS time
Real-Time Systems: Operating Principle Epicenter Seismic waves (~ 3.5 Km/s) Seismic Network Information (speed of light) Target P-Arrival Telemetry and computing Seismic Early Warning ~ 10 sec T0 TP TS time
Real-Time Systems: Operating Principle Epicenter Seismic waves (~ 3.5 Km/s) Seismic Network Information (speed of light) Target P-Arrival Telemetry and computing Seismic Early Warning ~ 10 sec Immediate post-event ~ 1 min T0 TP TS time
Seismic Alert Management System t 0 10-20 seconds 100-200 seconds time
Seismic Alert Management System t 0 10-20 seconds 100-200 seconds time Early Warning Data Array: location, magnitude (with uncertainty) Real-time data updating and dissemination of the information (progressive refinement of reliability of estimated parameters) Distributed process of analysis and decision
Seismic Alert Management System t 0 10-20 seconds 100-200 seconds time Early Warning Data Array: location, magnitude (with uncertainty) Real-time data updating and dissemination of the information (progressive refinement of reliability of estimated parameters) Distributed process of analysis and decision Post Event Warning Shake maps Preliminary Scenario: assessment and spatial interpolation of ground motion Detailed Scenario: fracture and wave propagation simulation based on simplified source/ propagation model
Campania region: historical seismicity
Campania region: historical seismicity Irpinia 1980 (Ms=6.9)
Rate of occurrence Instrumental data (Boschi et al, 2003) M>4.0 1 event every 1.5 years M>5.0 1 event every 4 years M>6.0 1 event every 32 years Probability map of moderate to large earthquakes (M 5.5) in Italy for the next 10 years (Cinti et al., G3, 2005)
INGV catalogue (1981-2002), M>2.5
The Irpinia Seismic Network (ISNet) Early Warning Network 28 sites Osiris 24-bit Data Logger 6 channels: 3 accelerometers 3 seismometers (Short Period or Broad Band) Real time data analysis INGV catalogue (1981-2002), M>2.5
The Irpinia Seismic Network (ISNet) Shake Map Network 11 sites Osiris 24-bit Data Logger 6 channels: 3 accelerometers 3 seismometers (Short Period or Broad Band) On request data analysis Early Warning Network 28 sites Osiris 24-bit Data Logger 6 channels: 3 accelerometers 3 seismometers (Short Period or Broad Band) Real time data analysis
The Irpinia Seismic Network (ISNet)
The Irpinia Seismic Network (ISNet) NCC LCC Three levels of data acquisition and analysis Stations (data loggers) Local Control Centers (sub-nets) Network Control Weber et al., 2006
The Irpinia Seismic Network (ISNet) NCC LCC Three levels of data acquisition and analysis Stations (data loggers) Local Control Centers (sub-nets) Network Control Weber et al., 2006 1.5-3.5 s eqk at 4-16 km depth 3-5 s Latency & computation 16-18 s at 60 Km 22-24 s at 80 Km 28-30 s at 100 Km To TP first TS target time
Real-Time Earthquake Location volume defined by stations without arrivals conditional EDT surface B A Satriano et al., 2006 First station detects arrival constraint is Voronoi cells
Real-Time Earthquake Location tt B tt A t now t A conditional EDT surface B A Satriano et al., 2006 Wavefront expands EDT surfaces deform, constraint improves
Real-Time Earthquake Location tt B tt A = t B t A true EDT surface B A Satriano et al., 2006 Second station detects arrival constraint includes EDT surface
Real-Time Earthquake Location Satriano et al., 2006 Third station detects arrival constraint is mainly EDT surfaces
Real-Time Earthquake Location (RTLoc)
Relationship between Magnitude and Log Peak Displacement Zollo and Lancieri, 2006 European Acelerometric Database (Ambraseys et al.2004) 207 events with 4 MW 7.4 (Kokaeli,1999) 376 3-component recordings Epicentre distance 50 km Low-pass filter at 3 Hz
Relationship between Magnitude and Log Peak Displacement Zollo and Lancieri, 2006 P-phase PGA (m/s ) 2 S-phase Time (s)
Relationship between Magnitude and Log Peak Displacement Zollo and Lancieri, 2006 P-phase PGA (m/s ) 2 S-phase Time (s) 0 P2-Sec 0 S1-Sec 0 S2-Sec Log(PGD)t(R=10Km) -1-2 -3-4 -1-2 -3-4 -1-2 -3-4 -5 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5-5 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 Magnitude (Mw) -5 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5
Real-Time Magnitude Estimation 41 30 41 00 40 30 km 0 50 40 00 14 00 14 30 15 00 15 30 16 00 16 30 Simulation study for the calculation of the magnitude in real time Earthquake of M=5.5 Low-pass filtered Ground displacement Synthetics are computed using the discrete wave number method of Bouchon(1981) coded by Coutant (1989) (AXITRA) Mag M5.5 Time from T0 (sec)
Rapid Shake Map Convertito et al., 2006
The End More details can be found on the following articles: A.Zollo, M.Lancieri and S.Nielsen (2006). Earthquake Magnitude Estimation From Peak Amplitudes of Very Early Seismic Signals on Strong Motion Records Geophysical Research Letters C.Satriano, A.Lomax and A.Zollo (2006). Optimal, real-time earthquake location for early Warning In P. Gasparini et al. editors, Seismic Early Warning, Springer-Verlag. E. Weber, G. Iannaccone, A. Zollo, A. Bobbio, L. Cantore, M. Corciulo, V. Convertito, M. Di Crosta, L. Elia, A. Emolo, C. Martino, A. Romeo, C. Satriano (2006). Development and testing of an advanced monitoring infrastructure (ISNet) for seismic early-warning applications in the Campania region of southern Italy In P. Gasparini et al. editors, Seismic Early Warning, Springer-Verlag. I.Iervolino, V. Convertito, M. Giorgio, G Manfredi and A Zollo (2006). Real-time risk analysis in hybrid earthquake early warning systems J.Earthq. Eng
Attenuation with distance Attenuation law for PDt log (P Dt) = f(m, R) = C + B M + A log(r) The observed peak amplitude are corrected for the distance effect by normalizing them to a reference distance of 10Km log ( P D 10Km t ) = log(p D obs t ) Alog ( R 10 )