1 / 31

Yan Y. Kagan Dept. Earth and Space Sciences, UCLA, Los Angeles, CA 90095-1567, ykagan@ucla.edu , http://scec.ess.ucla.e

Yan Y. Kagan Dept. Earth and Space Sciences, UCLA, Los Angeles, CA 90095-1567, ykagan@ucla.edu , http://scec.ess.ucla.edu/ykagan.html. Testing spatial distribution forecast of point processes with application to seismology. http://moho.ess.ucla.edu/~kagan/Score .ppt.

xuan
Download Presentation

Yan Y. Kagan Dept. Earth and Space Sciences, UCLA, Los Angeles, CA 90095-1567, ykagan@ucla.edu , http://scec.ess.ucla.e

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Yan Y. Kagan Dept. Earth and Space Sciences, UCLA, Los Angeles, CA 90095-1567, ykagan@ucla.edu, http://scec.ess.ucla.edu/ykagan.html Testing spatial distribution forecast of point processes with application to seismology http://moho.ess.ucla.edu/~kagan/Score.ppt

  2. CMT catalog: Shallow earthquakes, 1976-2005

  3. The table below and accompanying plots are calculated on 2007/ 4/19 at midnight Los Angeles time. The last earthquake with scalar seismic moment M>=10^17.7 Nm (Mw>=5.8) entered in the catalog occurred in the region 0.0 > LAT. > -60.0, -170.0 > LONG. > 110.0 on 2007/ 4/16 at latitude -57.89 and longitude 148.14, Mw = 6.42 ____________________________________________________________________ LONG-TERM FORECAST | SHORT-TERM Probability Focal mechanism | Probability Probability M>5.8 T-axis P-axis M>5.8 ratio eq/day*km^2 Pl Az Pl Az eq/day*km^2 Time- Longitude | | | Rotation Time- dependent/ | Latitude | | | angle dependent independent v v v v degree ……………………………………………………………………………………………………… 154.0 -7.5 1.49E-08 67 24 16 251 58.6 6.17E-11 4.13E-03 154.5 -7.5 2.02E-08 68 18 10 262 60.9 7.17E-11 3.54E-03 155.0 -7.5 2.60E-08 75 18 13 222 28.5 1.53E-09 5.88E-02 155.5 -7.5 3.51E-07 71 20 19 210 8.7 7.66E-07 2.2 156.0 -7.5 6.72E-08 76 20 14 216 20.0 8.49E-07 13. 156.5 -7.5 3.10E-08 76 28 13 231 41.4 5.12E-07 17. 157.0 -7.5 1.90E-08 1 327 16 236 44.7 1.37E-07 7.2 157.5 -7.5 8.92E-09 77 45 13 229 46.6 4.57E-09 0.51 158.0 -7.5 7.42E-09 79 60 11 228 45.3 9.06E-11 1.22E-02 158.5 -7.5 1.05E-08 49 147 4 52 54.2 1.66E-10 1.58E-02 159.0 -7.5 7.64E-09 47 147 4 242 58.8 1.07E-10 1.41E-02 ………………………………………………………………………………………………………

  4. NW Pacific --Here we demonstrate forecast effectiveness: displayed earthquakes (108 events) occurred after smoothed seismicity forecast had been calculated.

  5. SW Pacific --Here we demonstrate forecast effectiveness: displayed earthquakes (170 events) occurred after smoothed seismicity forecast had been calculated.

  6. 1080 events (10 times the actual number of events in 2004-06) have been simulated, using the forecast density. Some points overlap; in low activity area we specified some low-level density to avoid surprises.

  7. Long-term Forecast Efficiency Evaluation • We simulate synthetic catalogs using smoothed seismicity map. • Likelihood function for simulated catalogs and for real earthquakes in the time period of forecast is computed. • If the `real earthquakes’ likelihood value is within 2.5—97.5% of synthetic distribution, the forecast is considered successful. Kagan, Y. Y., and D. D. Jackson, 2000. Probabilistic forecasting of earthquakes, Geophys. J. Int., 143, 438-453.

  8. Cumulative likelihood score distribution Questions: 1. How fast approaches Gaussian; 2. Standard error behavior; 3. Influence of cell size

  9. Earthquake potential in and around China: Estimated from past earthquakes, Yufang Rong and David D. Jackson, GRL, 29(16), 2002.

  10. Kossobokov, 2006. Testing earthquake prediction methods: ``The West Pacific short-term forecast of earthquakes with magnitude MwHRV \ge 5.8", Tectonophysics, 413(1-2), 25-31. See also Kagan & Jackson, TECTO, 2006, pp. 33-38.

  11. Testing alarm-based earthquake predictionsJ. Douglas Zechar and Thomas H. JordanGeophys. J. Int. (2008) 172, 715–724 Results of Molchan trajectory/area skill score analysis for PI (squares) and NSHM (triangles) relative to the RI reference model. Top panel shows complete Molchan trajectories for both predictions and bottom panel shows corresponding area skill score curve. Each plot also shows the α = 1, 5 and 10 per cent critical boundaries.

  12. Calculation of Information Score from EQ density table (log-likelihood) Cell No Long Lat EQ density Cell Probab. EQ Gain EQ size 04-06 77-03 35 107 127.00 7.00 0.001766055 12254 0.601328269 0 68.1973 6 67 40 143.00 40.50 0.001607535 9388 0.419336865 0 62.0759 1 84 30 151.50 45.50 0.001492347 8653 0.358811419 0 57.6279 5 25 73 122.00 24.00 0.001455427 11278 0.456091624 1 56.2022 10 81 32 150.00 44.50 0.001429419 8805 0.349718405 0 55.1979 8 80 32 149.50 44.50 0.001405341 8805 0.343827435 0 54.2681 9 32 96 125.50 12.50 0.001395516 12053 0.467368764 0 53.8887 6 78 33 148.50 44.00 0.001395333 8881 0.344325587 0 53.8816 8

  13. Error diagram (ED), NW Pacific

  14. Error diagram, SW Pacific

  15. Calculation of Information Score from EQ density table (log-likelihood)

  16. Calculation of Information Score from EQ density table (log-likelihood)

  17. Calculation of Information Score from EQ density table (log-likelihood), I=2.3644 Several other log-likelihood scores can be calculated: • Gain for cells where earthquakes occurred (I=2.3625); (2) Gain for actual EQ location (I=2.4014); (3) Gain for simulated events (I=2.3609); (4) Info score for EQ numbers in Error Diagram (I=3.0922) – Kagan, PAGEOPH, 2007. These scores converge if cells are becoming small, or the EQ number increases.

  18. Distribution of log-likelihood score Blue -- SW Red -- NW

  19. Two segments ED

  20. Two segments ED, dependence on slope of the first segment Magenta – standard error; Blue –coefficient of skewness; Green – coefficient of kurtosis.

  21. World seismicity: 1990 – 2000 (PDE)

  22. Error diagram with CMT and PDE concentrations Magenta – CMT 1977-2003; Cyan – PDE 1968-2006; Red – forecast based on CMT 1977-2003; Blue – EQs 04-06; Green – simulation; Black – Poisson rate.

  23. Error diagram with forecast used as template (one simulation) Magenta – CMT 1977-2003; Cyan – PDE 1968-2006; Red – forecast based on CMT 1977-2003; Blue – EQs 04-06; Green – simulation; Black – Poisson rate.

  24. Error diagram with forecast used as template (10 simulations) Magenta – CMT 1977-2003; Cyan – PDE 1968-2006; Red – forecast based on CMT 1977-2003; Blue – EQs 04-06; Green – simulations; Black – Poisson rate.

  25. Molchan, G., and V. Keilis-Borok, 2008.Earthquake Prediction: Probabilistic Aspect,21 pages, http://arxiv.org/abs/0802.1126v1 Typical examples of \phi (strategy) that are used at the research stage are max(n, \tau) and sum (n + \tau).

  26. END Thank you

  27. Abstract We propose a new method to test the effectiveness of spatial point process forecast based on a log-likelihood score for a predicted point density and the information gain for events that actually occurred in the test period. The method largely avoids simulation use and allows to calculate the information score for each event as well as the standard error of each forecast. As the number of predicted events increases, the score distribution approaches the Gaussian law. The degree of its similarity to the Gaussian distribution can be measured by the computed coefficients of skewness and kurtosis. To display the forecasted point density and the point events we use an event concentration diagram or a variant of the Error Diagram (ED).

  28. Abstract (cont.) We demonstrate the application of the method by using our long-term forecast of seismicity in two west Pacific regions. We compare the Error Diagram (ED) for these regions with simplified diagrams based on two-segment approximations. Since in these regions earthquakes are concentrated in narrow subduction belts, the use of the forecast density as a template for the ED is a more convenient display technique. We also show, using simulated event occurrence, that some proposed criteria for measuring forecast effectiveness at EDs, would be strongly biased for small number of events.

  29. Frequency plot for likelihood score distribution

  30. Distribution of log-likelihood score Blue -- SW Red -- NW

  31. Two segments ED

More Related