1 / 32

Seismological Society of America Meeting 18 April 2006

SCEC Community Modeling Environment (SCEC/CME): SCEC TeraShake Platform: Dynamic Rupture and Wave Propagation Simulations. Seismological Society of America Meeting 18 April 2006. SCEC/CME Project.

earled
Download Presentation

Seismological Society of America Meeting 18 April 2006

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. SCEC Community Modeling Environment (SCEC/CME): SCEC TeraShake Platform: Dynamic Rupture and Wave Propagation Simulations Seismological Society of America Meeting 18 April 2006

  2. SCEC/CME Project Goal:To develop a cyberinfrastructure that can support system-level earthquake science – the SCEC Community Modeling Environment (CME) Support:5-yr project funded by the NSF/ITR program under the CISE and Geosciences Directorates Oct 1, 2001 – Sept 30, 2006 NSF CISE GEO SCEC/ITR Project ISI USGS Information Science Earth Science SDSC IRIS SCEC Institutions www.scec.org/cme

  3. SCEC/CME – All Hands Meeting

  4. SCEC/CME Focus On Seismic Hazard Analysis SCEC/CME System aims to extend and enhance geosciences work already performed in the area of seismic hazard analysis. Metadata for Map: IMT: Peak Acceleration POE: 10% TimeSpan: 50 Years IMR: 1) Boore, Joyner, and Fumal (BJF; 1993, 1994a) with later modifications to differentiate thrust and strike-slip faulting (Boore et al., 1994b), 2) Sadigh et al. (1993) and 3) Campbell and Bozorgnia (1994). ERF: … (and more)

  5. Pathway 4: Ground motion inverse problem Pathway 3: Physics-based earthquake forecasting Pathway 2: Ground motion simulation Other Data Geology Geodesy Unified Structural Representation Invert 4 Faults Motions Stresses Anelastic model Ground Motions FSM RDM AWM SRM 3 2 Earthquake Forecast Model Attenuation Relationship Intensity Measures 1 AWP = Anelastic Wave Propagation SRM = Site Response Model FSM = Fault System Model RDM = Rupture Dynamics Model SCEC/CME Computational Pathways Pathway 1: Standard Seismic Hazard Analysis

  6. SCEC Computational Platform Concept • Definition • A large-scale implementation of computational pathways within a specific computational system (hardware + software + expertise) for producing specific knowledge • Implied components • Validated simulation software and geophysical models • Broadly useful simulation capabilities • Imports from other systems. Exports to other Systems • IT/geoscience collaboration involved in operation • Access to High-performance hardware • May use Workflow management tools

  7. Large Scale Simulation-based Seismic Hazard Computational Platform Development • 2003 • OpenSHA • 2004 • OpenSHA, TeraShake • 2005 • OpenSHA, TeraShake, CyberShake • 2006 • OpenSHA, TeraShake, CyberShake, Earthworks

  8. TeraShake Platform Anelastic Wave Propagation Capabilities Investigate ground motion effects for several different large Southern San Andreas ruptures.

  9. TeraShake Simulation area • 600 km x 300 km x 80 km • dx=200m • Mesh of 1.8 Billion cubes • 0.011 sec time step, 20,000 time steps: 3 min Kinematic source: Cajon Creek to Bombay Beach (or back -60 sec source duration -18,886 point sources, each 6,800 time steps in duration SCEC/CME

  10. Scenario Earthquake Simulations - TeraShake TeraShake 1.3 Cumulative PGV and SA 3.0 Maps – Kim Olsen et al (AWM), Amit Chourasia et al (Viz)

  11. Scenario Earthquake Simulations - TeraShake TeraShake 1.2 and 1.3 Cumulative PGV – Kim Olsen et al (AWM), Amit Chourasia et al (Viz)

  12. Peak Displacements TeraShake N to S Rupture (left) TeraShake S to N Rupture (right)

  13. Particle Velocities Along N50E Profile Largest Peak Motions above ridge between SG and LA Basins

  14. Causes of ‘Hotspot’ in Los Angeles • Buildup of forward directivity pulse • Excitation of guided waves upon striking S.B. basin • Amplification as sedimentary waveguide narrows between Puente Hills and San Gabriel Mtns. • Lateral and vertical focusing effects from 2D horizontally incident plane waves produce amplification patterns similar to those from TS1.3 and TS1.4 • Additional effects possibly from Airy phase of Love wave • Nonlinear effects likely to decrease peak motions

  15. TeraShake Platform Dynamic Rupture Capabilities Incorporate more physics into the TeraShake simulations by introduction of a dynamic rupture-based source description.

  16. TeraShake-2 Dynamic Simulations • RDM run on an inner domain, containing the entire rupture as well as the relevant surrounding environment….So, in addition to rupture physics parameters (yield friction, d0, etc) • Need to include inhomogeneous crust on either side of rupture • This will require validation of rupture simulations • Need absorbing boundaries (PML) • To run the rupture to completion, without need for a large inner domain • Need validation of code coupling scheme • Current plan is to save the rupture and treat it a kinematic source for AWP • This raises issues of rescaling, smoothing, registering • Need to choose a rupture representation format that AWP can ingest

  17. Large Scale Dynamic Rupture Simulations – Terashake 2

  18. TeraShake 2 Simulation Area Spontaneous Rupture Simulation Domain

  19. Olsen Dynamic Rupture Simulations

  20. Olsen Dynamic Rupture Simulation Movies Slip Animation

  21. TeraShake 1 vs TeraShake 2 PGV Map

  22. Planned Application and Development for TeraShake Platform

  23. TeraShake 3 Simulations for 2006 0-3 Hz Deterministic Earthquake Ground Motion Prediction

  24. CMU Hercules Tool chain Highly scalable AWM software tools that uses new technique called In-situ parallel mesh generation David O’Hallaron et al (CMU) Etree Mesh Representation Jacobo Bielak et al (CMU) AWM

  25. Finite Element Dynamic Rupture CodesPeak Surface Velocity Ratio (Topography/No Topography)

  26. Scenario Earthquake Simulations – Puente Hills Peak SA 2.0 magnitude Map Velocity Y Component Animation Puente Hills Simulation Scenario Earthquake (10 Hz) Robert Graves (AWM), Amit Chourasia et al (Viz)

  27. Full 3-D Tomography Model (Preliminary) Po Chen – Inversion-based updates to SCEC CVM3.0 Velocity Model.

  28. SCEC Education and Outreach

  29. Outreach To Geosciences, Computer Sciences, and Public • SC04 • Invited Presentations USC, SDSC Booths • AGU 2004 • Booth, Presentations, Posters • Globus World 2005 • Invited Presentation • SSA 2005 • Booth, Invited Presentations, Posters • GRIDS Workshop 2005 • Invited Presentation • Unavco/IRIS Annual Meeting 2005 • CIG Workshop Presentation, Web Services Workshop, Posters • SCEC Earthquake Spectrum Press Conference 2005 • Earthscope Meeting 2005 • Poster, IT Workshop • GEON All Hands Meeting 2005 • Posters • National Forum for Geoscience Information Technology (FGIT) 2005 • Poster • GSA 2005 • Booth, Poster • SC05 • Invited Presentations USC, SDSC, TeraGrid Booths

  30. Computational Platform Concept New Physics-based Simulations (Validation) Observation Platform Capability Computing Validates Large Scale Scenarios (Full Inner/Outer Scales) PetaShake Platform Data Intensive Capability Computing Applies Seismic Hazard (Application) CyberShake Platform Data Intensive Capacity Computing Community Access (Science Gateways) Computing Requirements SCEC Science SCEC Platforms

  31. Vertical integration SCEC Community Modeling EnvironmentA grid-enabled collaboratory for system-level earthquake science 6 Intelligent services (smart assistants) Integrated system tools (workbench/dashboard) Workflow management Domain applications (webservices/applications) Resource sharing (grids) Hardware (computing, networking, storage) 5 4 3 2 1 Cyberinfrastructure Layering

More Related