1 / 33

J. Louie and L. Preston Nevada Seismological Laboratory S. Larsen

Computing Scenario Shaking for Nevada Urban Areas with the Model Assembler Community Modeling Environment (MA-CME). J. Louie and L. Preston Nevada Seismological Laboratory S. Larsen Lawrence Livermore National Laboratory Liz Lenox, Rei Arai, and Amr Wakwak

colman
Download Presentation

J. Louie and L. Preston Nevada Seismological Laboratory S. Larsen

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. Computing Scenario Shaking for Nevada Urban Areas with the Model Assembler Community Modeling Environment (MA-CME) J. Louie and L. Preston Nevada Seismological Laboratory S. Larsen Lawrence Livermore National Laboratory Liz Lenox, Rei Arai, and Amr Wakwak Dept. Geological Sciences and Engineering, UNR

  2. ModelAssembler • A code to stitch together existing regional geophysical and geological data sets. • Generates multi-gigabyte E3D input grids.

  3. The ModelAssembler Community Modeling Environment (MA-CME) • Graphical user interface for setting up MA/E3D “.in” files. • Provides tutorial help to successful setup. • Data sets, stations, grids all defined with geographic lat/lon (have to decide on WGS84, NAD27, etc.). • Easy configuration of problems at all scales and computational difficulties- 10 Mb to 100 Gb. • After setup, “portal pack” is downloaded to cluster and run. • Open source: www.seismo.unr.edu/ma • Limits: flat earth; no topography; 1-D variations within basins & bedrock; source and stations in same grid.

  4. MA-CME provides a tutorial GUI to E3DSet up a Grid:

  5. Define Execution Platform

  6. Define Earthquake Source

  7. Define Rule Base

  8. Arrange Basin-Thickness (& Geotech) Inputs

  9. Request Output Data

  10. Three Nevada Scenarios Reno 2) Genoa FZ M7.5 1) Frenchman M5.0 3) Furnace Cr FZ M7.5 Las Vegas

  11. Frenchman Mountain Fault • 4x4-km N-striking M5.0 scenario • West-dipping oblique-normal • Late Quaternary (<130 ka) • Less than 0.2 mm/yr • By Amr Wakwak (Fulbright PLUS Program)

  12. Frenchman Mountain Fault Scott et al. (2006) Vs30 model from 1100+ wells, 79 Vs profiles. Thanks to W. Taylor of UNLV & G. Wagoner of LLNL

  13. Frenchman Mountain Fault Langenheim et al. model from gravity, refraction, a few deep wells.

  14. Frenchman Mountain Fault Animation combines inputs and outputs. E3D computation at 0.3 Hz on a 144x137 grid with dh=308 m. Grid resolution is a good match to cell phone displays.

  15. Frenchman Mountain Fault PGV >5 cm/s above west-dipping fault at 0.3 Hz. Sub-basins contain the greatest shaking.

  16. Frenchman Mountain Fault Schools closest to the fault experience greatest PGV but shortest durations.

  17. Genoa Fault – Mount Rose Fault Systems Virginia Lake • 80 km long, M7.5 • East-dipping oblique-normal • Latest Quaternary (<15 ka) • 1.0-5.0 mm/yr • Northward rupture scenario from W. Fork Carson Riv. to Virginia Lake • By Rei Arai W. Fork Carson Riv.

  18. Genoa – Mount Rose Scenario:MA-CME adds Geotech & Basins Scott et al. (2004) Vs30 transect. Pancha measurements at ANSS stations. “Quadrant” interpolation between scattered measurements, respecting geologic map. Vs30=500 m/s assumed on sediment, 760 m/s (white and cyan) on rock.

  19. Genoa – Mount Rose Scenario:MA-CME adds Geotech & Basins Three datasets, not well stitched together: • Geologic map with basin depths assumed from bedrock proximity, in California. • Jachens et al. USGS basin gravity inversions for Basin & Range; sedimentary plus Tertiary volcanic basins. • Abbott & Louie (2000) Reno basin gravity study.

  20. Genoa Fault – Mount Rose Fault Systems Dataset mismatches plain in model renderings. Extreme directivity effect. E3D computation at 0.3 Hz on a 120x200 grid with dh=500 m. Olsen and Day Q model included.

  21. Genoa Fault – Mount Rose Fault Systems West Reno basin, above end of rupture, has 5x the PGV of any other basin near fault. Highest shaking not correlated with basin depth- Tahoe artificially deep in this model. Dataset boundaries not interfering with these conclusions.

  22. Genoa Fault – Mount Rose Fault Systems West Reno stations have high PGV and long durations.

  23. Northern Death Valley Fault Zone Ubehebe Crater • 80 km long • Dextral strike-slip • 1.0- 5.0 mm/yr • Latest Quaternary (<15 ka) • Two M7.5 scenarios: • Rupture from SE, epicenter at Furnace Crk. • Rupture from NW, epicenter at Ubehebe Ctr. • By Liz Lenox Furnace Creek

  24. Northern Death Valley Fault Zone Jachens USGS geologic map sets default soil and rock Vs30s. Scott et al. (2006) Vs30 model from 1100+ wells, 79 Vs profiles inserted for LVV.

  25. Northern Death Valley Fault Zone Jachens Basin & Range gravity inversions for sedimentary & volcanic basins. Timber Mtn. caldera & rifts up to 8 km deep. Langenheim et al. model from gravity, refraction, a few deep wells inserted for LVV.

  26. No. Death Val. - Furnace Cr. Fault Zone Animation combines inputs and outputs. E3D computation at 0.3 Hz on a 281x251 grid with dh=1 km. Infinite Q. The two rupture scenarios have very different effects. Get the scenario animations on your cell phone from www.seismo.unr.edu/ma/

  27. Northern Death Valley Fault Zone For rupture away from city, PGV <<1 cm/s in LVV. Clear directivity in this long-period simulation.

  28. Northern Death Valley Fault Zone For rupture toward city, PGV >1 cm/s in LVV. Clear directivity in this long-period simulation. Basins between DV and LVV are spreading the directivity effect to wider angles from the fault strike.

  29. Northern Death Valley Fault Zone Ratio of PGV from rupture toward LVV, over away. Largest amplifications due to directivity are in bedrock. Some margins of LVV also show 1500% amplifications.

  30. Northern Death Valley Fault Zone Compare PGV of 2679 areas, each 1-km2, in LVV against Vs30 and Z2.0.

  31. Vs30 affecting 1.5-km waves! Northern Death Valley Fault Zone

  32. Northern Death Valley Fault Zone Uncertainty in PGV

  33. Conclusions • Students can use MA-CME to set problems up, but plotting and visualizing are still lots of work. • Scenario Results: • M7.5 Genoa–Mt. Rose– Reno basin at end of rupture gets 5 times the PGV of any other. • M5.0 Frenchman Mtn.– high PGV accompanied by short duration, concentrated by deep basin structure. • M7.5 No. Death Valley– rupture toward city produces up to 15 times the PGV as rupture away. • At long periods, 0.3 Hz, E3D is producing huge uncertainties in PGV.

More Related