1 / 14

Modeling swarms: A path toward determining short-term probabilities

Modeling swarms: A path toward determining short-term probabilities. Andrea Llenos USGS Menlo Park Workshop on Time-Dependent Models in UCERF3 8 June 2011. Outline. Motivation: Why are swarms important for UCERF? Where things stand now Characteristics of swarms

opal
Download Presentation

Modeling swarms: A path toward determining short-term probabilities

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. Modeling swarms:A path toward determining short-term probabilities Andrea Llenos USGS Menlo Park Workshop on Time-Dependent Models in UCERF3 8 June 2011

  2. Outline • Motivation: Why are swarms important for UCERF? • Where things stand now • Characteristics of swarms • Detecting swarms (retrospectively) • What needs to be done • Detecting swarms (prospectively) • Implementation • As ETAS add-on? • As a data assimilation application? Observed seismicity rate Background seismicity rate Aftershock sequences

  3. Time-dependent background rates are needed to account for rate changes due to external (aseismic) processes 2000 Izu Islands swarm (magma/fluids) 2000 Vogtland/Bohemia swarm (fluids) Hainzl and Ogata (2005) Lombardi et al. (2006) 2003-2004 Ubaye swarm (fluid-flow) Daniel et al. (2011)

  4. Salton Trough Time-dependent background rate matches observed seismicity better than stationary ETAS model Transformed Time Llenos and McGuire (2011)

  5. Characteristics of swarms • Increase in seismicity rate above background without clear mainshock • Don’t follow empirical aftershock laws • Bath’s Law • Omori’s Law • These characteristics make them appear anomalous to ETAS Holtkamp and Brudzinski (2011)

  6. Detecting swarms in an earthquake catalog Swarms associated with aseismic transients 2005 Obsidian Buttes, CA (1985-2005, SCEDC) 2005 Kilauea, HI (2001-2007, ANSS) 2002, 2007 Boso, Japan (1992-2007, JMA) Slow slip events on the subduction plate interface off of Boso, Japan observed by cGPS, tiltmeter Shallow aseismic slip on a strike-slip fault in southern CA observed by InSAR and GPS Lohman & McGuire (2007) Slow slip events on southern flank of Kilauea volcano in HI observed by GPS Wolfe et al. (2007) Ozawa et al. (2007)

  7. Data analysis: ETAS model optimization Swarms associated with aseismic transients 2005 Obsidian Buttes, CA (1985-2005, SCEDC) 2005 Kilauea, HI (2001-2007, ANSS) 2002, 2007 Boso, Japan (1992-2007, JMA) • Optimize ETAS model to fit catalog prior to swarm and extrapolate fit through remainder of catalog • Calculate transformed times (~ ETAS predicted number of events in a time interval) • Cumulative number of events vs. transformed time should be linear if seismicity behaving as a point process • Positive deviations occur when more seismicity is being triggered in a time interval than ETAS can explain 2005 Kilauea

  8. Swarms appear as anomalies relative to ETAS 2005 Obsidian Buttes 2002, 2007 Boso, Japan 2005 Kilauea

  9. A path toward determining short-term probabilities • Build off of ETAS-based forecasts • Detect that a swarm is occurring • Has been done retrospectively • Prospectively? • During the swarm • Re-estimate the background rate (and other parameters?) • Re-calculate short-term probabilities • How often? 1x? 2x? Every 5 days? 10 days? • Identify when the swarm is over • Return to pre-swarm background rate? • More sophisticated approaches (e.g., data assimilation)?

  10. Data Assimilation Algorithms • Combines dynamic model with noisy data (e.g. seismicity rates) to estimate the temporal evolution of underlying physical variables (states) • Examples: Kalman filters, particle filters • Applications in navigation, tracking, hydrology Welch & Bishop (2001)

  11. Data Assimilation Example • State-space model based on rate-state equations • States: stressing rate, rate-state state variable g • Algorithm: Extended Kalman Filter • Approach: Optimize ETAS for the catalog, subtract ETAS predicted aftershock rate to obtain time-dependent background rate, use data assimilation algorithm to estimate stressing rate and detect transients that trigger swarms Llenos and McGuire (2011)

  12. A path toward determining short-term probabilities • Build off of ETAS-based forecasts • Detect that a swarm is occurring • Has been done retrospectively • Prospectively? • During the swarm • Re-estimate the background rate (and other parameters?) • Re-calculate short-term probabilities • How often? 1x? 2x? Every 5 days? 10 days? • Identify when the swarm is over • Return to pre-swarm background rate? • More sophisticated approaches (e.g., data assimilation)?

  13. Outline • Why are swarms important for UCERF? • Need time-dependent background rate (mu) to model earthquake rates observed in catalogs accurately • Salton Trough • Ubaye France • CampeiFlagrei • Vogtland Bohemia • Swarms prevalent in Salton Trough, volcanic regions like Long Valley, places where M>6 events have occurred • Characteristics of swarms • Don’t fit empirical models of aftershock clustering, appear anomalous • ETAS parameters change during swarms (primarily stationary background rate) • How to implement this to calculate short-term probabilities? • Where we are now • Detection (retrospective) • How they affect ETAS parameters • Outstanding issues that need to be addressed • Data assimilation?

More Related