Off-Network Earthquake Location by Earthquake Early Warning Systems: Methodology and Validation

ABSTRACT Regional source-based earthquake early warning systems perform three consecutive tasks: (1) detection and epicenter location, (2) magnitude determination, and (3) ground-motion prediction. The correctness of the magnitude determination is contingent on that of the epicenter location, and the credibility of the ground-motion prediction depends on those of the epicenter location and the magnitude determination. Thus, robust epicenter location scheme is key for regional earthquake early warning systems. Available source-based systems yield acceptably accurate locations when the earthquakes occur inside the real-time seismic network, but they return erroneous results otherwise. In this study, a real-time algorithm that is intended as a supplement to an existing regional earthquake early warning systems is introduced with the sole objective of ameliorating its off-network location capacity. The new algorithm combines measurements from three or more network stations that are analyzed jointly using an array methodology to give the P-wave slowness vector and S-phase arrival time. Prior to the S-phase picking, the nonarrival of the S phase is used for determining a minimum epicentral distance. This estimate is updated repeatedly with elapsed time until the S phase is picked. Thus, the system timeliness is not compromised by waiting for the S-phase arrival. After the S wave is picked, an epicentral location can be determined using a single array by intersecting the back-azimuth beam with the S-minus-P annulus. When several arrays are assembled, the back azimuth and P and S picks from all arrays are combined to constrain the epicenter. The performance of the array processing for back azimuth and S-wave picking is assessed using a large number of accelerograms, recorded by nine strong motion sensors of the KiK-net seismic network in Japan. The nine stations are treated as three distinct seismic arrays, comprising three stations each. Good agreement is found between array-based and catalog-reported parameters. Finally, the advantage of the new array methodology with respect to alternative schemes for back azimuth and distance is demonstrated.

GPS Derived Seismic Signals for Far Field Earthquake Epicenter Location Estimation

A reliable epicenter estimation method is proposed for Global Positioning System (GPS) derived seismic signal for far-field regional earthquake. The main contribution is the use of time-frequency analysis to estimate the time of arrival (TOA) using multilateration technique. The data from the 2004 Sumatra Andaman earthquake captured from four GPS continuously operating reference stations (GPS CORS) were used in the analysis. To validate the accuracy of the proposed method, the estimated epicenter location was compared with the data released by the United States Geological Survey (USGS). The estimated location shows an error of about 0.0572 degrees in latitude and 0.2848 degrees in longitude. The proposed analysis method could complement existing seismometer measurements, improve in understanding of geo-seismic phenomena, and plan future infrastructure development.

INTERMEDIATE TERM EARTHQUAKE PREDICTION BASED ON INTEREVENT TIMES OF MAINSHOCKS AND ON SEISMIC TRIGGERING

Two models, which contribute to the knowledge on intermediate term earthquake prediction are further examined, improved and applied. The first of these models, called Time and Magnitude Predictable (TIMAP) regional model is based on repeat times of mainshocks generated by tectonic loading on a network of faults which are located in a certain seismic region (faults’ region). The second model, called Decelerating-Accelerating Strain (D-AS) model, is based on triggering of a mainshock by its preshocks. Parameters of the TIMAP model have been specified for the Aegean area and applied by a backward test in 86 circular faults’ regions of this area. The test shows the validity of this time dependent model with 29% false alarms. Data concerning decelerating and accelerating seismic (Benioff) strain, which preceded 46 strong (M≥6.3) recent mainshocks in a variety of global seismotectonic regimes, show that the generation of a mainshock is triggered by quasi-static stress changes due to accelerating preshocks which occur in a broad (critical) region and by static stress changes due to the large number (frequency of occurrence) of small preshocks generated in a narrow (seismogenic) region. Retrospective predictions (postdictions) of these 46 mainshocks by the D–AS model confirms previous results concerning the prediction uncertainties (2σ) of the model in the origin time (± 2.5 years), epicenter location (≤ 150 km) and magnitude (± 0.4) of an ensuing mainshock with a probability ~ 80%. Information is also given on the successful prediction by the D-AS model of: 1) the Cythera strong (M = 6.9) earthquake which occurred on 8 January 2006 in the southwestern part of the Hellenic Arc and 2) of the Rhodes strong (M = 6.4) earthquake which occurred on 15 July 2008 in the Eastern part of this Arc. A backward combined application of both models in the Aegean area shows an uncertainty ≤120km in the epicenter location of an ensuing mainshock.

Appraising the Early-est earthquake monitoring system for tsunami alerting at the Italian Candidate Tsunami Service Provider

In this paper we present and discuss the performance of the procedure for earthquake location and characterization implemented in the Italian Candidate Tsunami Service Provider at the Istituto Nazionale di Geofisica e Vulcanologia (INGV) in Rome. Following the ICG/NEAMTWS guidelines, the first tsunami warning messages are based only on seismic information, i.e., epicenter location, hypocenter depth, and magnitude, which are automatically computed by the software Early-est. Early-est is a package for rapid location and seismic/tsunamigenic characterization of earthquakes. The Early-est software package operates using offline-event or continuous-real-time seismic waveform data to perform trace processing and picking, and, at a regular report interval, phase association, event detection, hypocenter location, and event characterization. Early-est also provides mb, Mwp, and Mwpd magnitude estimations. mb magnitudes are preferred for events with Mwp ≲ 5.8, while Mwpd estimations are valid for events with Mwp ≳ 7.2. In this paper we present the earthquake parameters computed by Early-est between the beginning of March 2012 and the end of December 2014 on a global scale for events with magnitude M ≥ 5.5, and we also present the detection timeline. We compare the earthquake parameters automatically computed by Early-est with the same parameters listed in reference catalogs. Such reference catalogs are manually revised/verified by scientists. The goal of this work is to test the accuracy and reliability of the fully automatic locations provided by Early-est. In our analysis, the epicenter location, hypocenter depth and magnitude parameters do not differ significantly from the values in the reference catalogs. Both mb and Mwp magnitudes show differences to the reference catalogs. We thus derived correction functions in order to minimize the differences and correct biases between our values and the ones from the reference catalogs. Correction of the Mwp distance dependency is particularly relevant, since this magnitude refers to the larger and probably tsunamigenic earthquakes. Mwp values at stations with epicentral distance Δ ≲ 30° are significantly overestimated with respect to the CMT-global solutions, whereas Mwp values at stations with epicentral distance Δ ≳ 90° are slightly underestimated. After applying such distance correction the Mwp provided by Early-est differs from CMT-global catalog values of about δ Mwp ≈ 0.0 ∓ 0.2. Early-est continuously acquires time-series data and updates the earthquake source parameters. Our analysis shows that the epicenter coordinates and the magnitude values converge within less than 10 min (5 min in the Mediterranean region) toward the stable values. Our analysis shows that we can compute Mwp magnitudes that do not display short epicentral distance dependency overestimation, and we can provide robust and reliable earthquake source parameters to compile tsunami warning messages within less than 15 min after the event origin time.

Appraising the Early-est earthquake monitoring system for tsunami alerting at the Italian candidate Tsunami Service Provider

In this paper we present the procedure for earthquake location and characterization implemented in the Italian candidate Tsunami Service Provider at INGV in Roma. Following the ICG/NEAMTWS guidelines, the first tsunami warning messages are based only on seismic information, i.e. epicenter location, hypocenter depth and magnitude, which are automatically computed by the software Early-est. Early-est is a package for rapid location and seismic/tsunamigenic characterization of earthquakes. The Early-est software package operates on offline-event or continuous-realtime seismic waveform data to perform trace processing and picking, and, at a regular report interval, phase association, event detection, hypocenter location, and event characterization. In this paper we present the earthquake parameters computed by Early-est from the beginning of 2012 till the end of December 2014 at global scale for events with magnitude M ≥ 5.5, and the detection timeline. The earthquake parameters computed automatically by Early-est are compared with reference manually revised/verified catalogs. From our analysis the epicenter location and hypocenter depth parameters do not differ significantly from the values in the reference catalogs. The epicenter coordinates generally differ less than 20 ∓ 20 km from the reference epicenter coordinates; focal depths are less well constrained and differ generally less than 0 ∓ 30 km. Early-est also provides mb, Mwp and Mwpd magnitude estimations. mb magnitudes are preferred for events with Mwp ≲ 5.8, while Mwpd are valid for events with Mwp ≳ 7.2. The magnitude mb show wide differences with respect to the reference catalogs, we thus apply a linear correction mbcorr = mb · 0.52 + 2.46, such correction results into δmb ≈ 0.0 ∓ 0.2 uncertainty with respect the reference catalogs. As expected the Mwp show distance dependency. Mwp values at stations with epicentral distance Δ ≲ 30° are significantly overestimated with respect the CMT-global solutions, whereas Mwp values at stations with epicentral distance Δ ≳ 90° are slightly underestimated. We thus apply a 3rd degree polynomial distance correction. After applying the distance correction, the Mwp provided by Early-est differs from CMT-global catalog values of about δ Mwp ≈ 0.0 ∓ 0.2. Early-est continuously acquires time series data and updates the earthquake source parameters. Our analysis shows that the epicenter coordinates and the magnitude values converge rather quickly toward the final values. Generally we can provide robust and reliable earthquake source parameters to compile tsunami warning message within less than about 15 min after event origin time.

On the footprints of a major Brazilian Amazon earthquake

Combining historical accounts and seismological studies, three hundred years of dormant information emerged as a source of the largest known seismic event that rocked Brazil since the beginning of our colonization. The probable epicenter location of the 1690 tremor lies on the left bank of the Amazon River, about 45 km downstream from the modern day Manaus. A year later, while passing this area, a missionary met witnesses of the tremor and observed remarkable changes in the topography and vegetation along the margin of the river. By 1692 another priest confirmed this event and the occurrence of large waves in the river, which led to the flooding of the Native Indians' terrains. The tremor spread seismic waves throughout the forest and shook indigenous constructions as far as one thousand kilometers away. A calculation of the seismic parameters shows an estimated magnitude of 7, a maximum intensity of IX MM and a felt area of about 2 million km2. Due to the long recurrence period for this type of tremor, the discovery of one of these events is valuable for seismic global intraplate studies. As for Brazil, it unravels the myth that the country was never hit by severe earthquakes.