j. louie, gbcvm workshop 1/14/2008 the model assembler community modeling environment and initial...
Post on 19-Dec-2015
219 views
TRANSCRIPT
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
The Model Assembler Community The Model Assembler Community Modeling Environment and Initial Modeling Environment and Initial Ground-Motion Computations for Ground-Motion Computations for
Reno and Las VegasReno and Las Vegas
J. Louie J. Louie andand L. Preston L. PrestonNevada Seismological LaboratoryNevada Seismological Laboratory
S. LarsenS. LarsenLawrence Livermore National LaboratoryLawrence Livermore National Laboratory
Liz Lenox, Rei Arai, Liz Lenox, Rei Arai, andand Amr Wakwak Amr WakwakDept. Geological Sciences and Engineering, UNRDept. Geological Sciences and Engineering, UNR
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
A code to stitch A code to stitch together existing together existing regional regional geophysical and geophysical and geological data geological data sets.sets.
Generates multi-Generates multi-gigabyte E3D gigabyte E3D input grids.input grids.
ModelModelAssemblerAssembler
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
The ModelAssembler Community The ModelAssembler Community Modeling Environment (MA-CME)Modeling Environment (MA-CME)
Graphical user interface for setting up MA/E3D “.in” files.Graphical user interface for setting up MA/E3D “.in” files. Provides tutorial help to successful setup.Provides tutorial help to successful setup. Data sets, stations, grids all defined with geographic Data sets, stations, grids all defined with geographic
lat/lon lat/lon (have to decide on WGS84, NAD27, etc.).(have to decide on WGS84, NAD27, etc.).
Easy configuration of problems at all scales and Easy configuration of problems at all scales and computational difficulties- 10 Mb to 100 Gb.computational difficulties- 10 Mb to 100 Gb.
After setup, “portal pack” is downloaded to cluster and After setup, “portal pack” is downloaded to cluster and run.run.
Open source: Open source: www.seismo.unr.edu/mawww.seismo.unr.edu/ma Limits:Limits: flat earth; no topography; 1-D variations within flat earth; no topography; 1-D variations within
basins & bedrock; source and stations in same grid.basins & bedrock; source and stations in same grid.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
MA-CME provides a tutorial GUI to E3DMA-CME provides a tutorial GUI to E3DSet up a Grid:Set up a Grid:
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Define Earthquake SourceDefine Earthquake Source
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Arrange Basin-Thickness (& Geotech) InputsArrange Basin-Thickness (& Geotech) Inputs
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Basins in the B&RBasins in the B&R
Four datasets, not well stitched Four datasets, not well stitched
together:together:
1)1) Geologic map with basin Geologic map with basin
depths assumed from bedrock depths assumed from bedrock
proximity, in California.proximity, in California.
2)2) Jachens et al. USGS basin Jachens et al. USGS basin
gravity inversions for Basin & gravity inversions for Basin &
Range; sedimentary plus Range; sedimentary plus
Tertiary volcanic basins.Tertiary volcanic basins.
3)3) Abbott & Louie (2000) Reno Abbott & Louie (2000) Reno
basin gravity study.basin gravity study.
4)4) Langenheim et al. Las Langenheim et al. Las
Vegas model from gravity, Vegas model from gravity,
refraction, a few deep wells.refraction, a few deep wells.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Request Output DataRequest Output Data
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Three Nevada ScenariosThree Nevada ScenariosThree Nevada ScenariosThree Nevada Scenarios
RenoReno
Las VegasLas Vegas
1) Genoa FZ M7.51) Genoa FZ M7.5
3) Furnace Cr FZ M7.53) Furnace Cr FZ M7.5 2) Frenchman M5.02) Frenchman M5.0
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Genoa – Mount Rose Genoa – Mount Rose Scenario:Scenario:
MA-CME adds MA-CME adds Geotech & BasinsGeotech & Basins
Three datasets, not well stitched Three datasets, not well stitched together around Reno or at Calif. together around Reno or at Calif. Border.Border.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Genoa – Mount Genoa – Mount Rose Scenario:Rose Scenario:MA-CME adds MA-CME adds
Geotech & BasinsGeotech & BasinsScott et al. (2004) Vs30 Scott et al. (2004) Vs30 transect.transect.
Pancha measurements at Pancha measurements at ANSS stations. ANSS stations.
““Quadrant” interpolation Quadrant” interpolation between scattered between scattered measurements, respecting measurements, respecting geologic map.geologic map.
Vs30=500 m/s assumed Vs30=500 m/s assumed on sediment, 760 m/s on sediment, 760 m/s (white and cyan) on rock.(white and cyan) on rock.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Genoa Fault – Mount Rose Fault SystemsGenoa Fault – Mount Rose Fault Systems
Dataset mismatches Dataset mismatches plain in model plain in model renderings.renderings.
Extreme directivity Extreme directivity effect.effect.
E3D computation at E3D computation at 0.3 Hz on a 0.3 Hz on a 120x200 grid with 120x200 grid with dh=500 m.dh=500 m.
Olsen and Day Q Olsen and Day Q model included.model included.
QuickTime™ and aPhoto - JPEG decompressor
are needed to see this picture.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Genoa Fault – Genoa Fault – Mount Rose Fault Mount Rose Fault
SystemsSystems
West Reno basin, above end of West Reno basin, above end of rupture, has 5x the PGV of any other rupture, has 5x the PGV of any other basin near fault.basin near fault.
Highest shaking not correlated with Highest shaking not correlated with basin depth- Tahoe artificially deep basin depth- Tahoe artificially deep in this model.in this model.
Dataset boundaries not interfering Dataset boundaries not interfering with these conclusions.with these conclusions.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Frenchman Mountain FaultFrenchman Mountain Fault
Langenheim et Langenheim et al. model from al. model from gravity, gravity, refraction, a few refraction, a few deep wells.deep wells.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Frenchman Mountain FaultFrenchman Mountain Fault
Scott et al. (2006) Scott et al. (2006) Vs30 model from Vs30 model from 1100+ wells, 79 1100+ wells, 79 Vs profiles. Vs profiles.
Thanks to W. Thanks to W. Taylor of UNLV & Taylor of UNLV & G. Wagoner of G. Wagoner of LLNLLLNL
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Frenchman Mountain FaultFrenchman Mountain Fault
Animation combines Animation combines inputs and outputs.inputs and outputs.
E3D computation at E3D computation at 0.3 Hz on a 0.3 Hz on a 144x137 grid with 144x137 grid with dh=308 m.dh=308 m.
Grid resolution is a Grid resolution is a good match to cell good match to cell phone displays.phone displays.
QuickTime™ and aVideo decompressor
are needed to see this picture.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Frenchman Mountain FaultFrenchman Mountain Fault
PGV >5 cm/s PGV >5 cm/s above west-above west-dipping fault at dipping fault at 0.3 Hz.0.3 Hz.
Sub-basins Sub-basins contain the contain the greatest shaking.greatest shaking.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Northern Death Valley Fault ZoneNorthern Death Valley Fault Zone
Jachens Basin Jachens Basin & Range gravity & Range gravity inversions for inversions for sedimentary & sedimentary & volcanic basins.volcanic basins.
Timber Mtn. Timber Mtn. caldera & rifts caldera & rifts up to 8 km up to 8 km deep. deep.
Langenheim et Langenheim et al. model from al. model from gravity, gravity, refraction, a few refraction, a few deep wells deep wells inserted for inserted for LVV.LVV.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Northern Death Valley Fault ZoneNorthern Death Valley Fault Zone
Jachens USGS Jachens USGS geologic map geologic map sets default soil sets default soil and rock Vs30s.and rock Vs30s.
Scott et al. Scott et al. (2006) Vs30 (2006) Vs30 model from model from 1100+ wells, 79 1100+ wells, 79 Vs profiles Vs profiles inserted for inserted for LVV. LVV.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
No. Death Val. - Furnace Cr. Fault ZoneNo. Death Val. - Furnace Cr. Fault Zone
Animation combines Animation combines inputs and outputs.inputs and outputs.
E3D computation at E3D computation at 0.3 Hz on a 281x251 0.3 Hz on a 281x251 grid with dh=1 km.grid with dh=1 km.
Infinite Q.Infinite Q.
The two rupture The two rupture scenarios have very scenarios have very different effects.different effects.
Get the scenario animations on your cell phone from www.seismo.unr.edu/ma/Get the scenario animations on your cell phone from www.seismo.unr.edu/ma/
QuickTime™ and aVideo decompressor
are needed to see this picture.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Northern Death Valley Fault ZoneNorthern Death Valley Fault Zone
For rupture For rupture toward city, toward city, PGV >1 cm/s in PGV >1 cm/s in LVV.LVV.
Clear directivity Clear directivity in this long-in this long-period period simulation.simulation.
Basins between Basins between DV and LVV are DV and LVV are spreading the spreading the directivity effect directivity effect to wider angles to wider angles from the fault from the fault strike.strike.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Northern Death Valley Fault ZoneNorthern Death Valley Fault ZoneCompare PGV of 2679 areas, each 1-kmCompare PGV of 2679 areas, each 1-km22, in LVV against Vs30 and Z2.0., in LVV against Vs30 and Z2.0.
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Northern Death Valley Fault ZoneNorthern Death Valley Fault ZoneFurnace Creek M7.5 Scenario Ruptures, GM in Las Vegas Valley
0.000
0.002
0.004
0.006
0.008
0.010
0.012
0.014
0.016
0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00
Site Vs30, km/s
0.3-Hz E3D PGV, meter/sec
Rupture from NW
Rupture from NW
Rupture from NW
Rupture from NW
Rupture from NW
Rupture from SE
Rupture from SE
Rupture from SE
Rupture from SE
Rupture from SE
Vs30 affecting 1.5-km waves!
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
Northern Death Valley Fault ZoneNorthern Death Valley Fault ZoneFurnace Creek M7.5 Scenario Ruptures, GM in Las Vegas Valley
0.000
0.002
0.004
0.006
0.008
0.010
0.012
0.014
0.016
-0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5
Basin Depth Z2.0, km
0.3-Hz E3D PGV, meter/sec
Rupture from NW
Rupture from NW
Rupture from NW
Rupture from NW
Rupture from NW
Rupture from SE
Rupture from SE
Rupture from SE
Rupture from SE
Rupture from SE
Un
cert
ain
ty in
PG
V
J. Louie, GBCVM WorkshopJ. Louie, GBCVM Workshop1/14/20081/14/2008
ConclusionsConclusions
Students can use MA-CME to set problems up, but Students can use MA-CME to set problems up, but plotting and visualizing are still lots of work.plotting and visualizing are still lots of work.
Scenario Results:Scenario Results: M7.5 Genoa–Mt. Rose–M7.5 Genoa–Mt. Rose– Reno basin at end of rupture gets 5 Reno basin at end of rupture gets 5
times the PGV of any other.times the PGV of any other. M5.0 Frenchman Mtn.–M5.0 Frenchman Mtn.– high PGV accompanied by short high PGV accompanied by short
duration, concentrated by deep basin structure.duration, concentrated by deep basin structure. M7.5 No. Death Valley–M7.5 No. Death Valley– rupture toward city produces up to 15 rupture toward city produces up to 15
times the PGV as rupture away.times the PGV as rupture away.• At long periods, 0.3 Hz, E3D is producing huge uncertainties At long periods, 0.3 Hz, E3D is producing huge uncertainties
in PGV.in PGV.