HighF 2018
This page documents the High-F activities and decisions for the 2018 verification and validation runs for three groups: Olsen et al. (AWP), Graves (RWG) and Taborda et al. (Hercules).
Contents
Mesh generation rules and parameter constraints
Velocities
- 1. Set Min Vs=500 m/s
- 2. If Vs was lower than 500 m/s and adjusted, then adjust Vp with original Vp/Vs ratio (so that we don’t have the automatic Vs/Vp of 3). We may want to set a minimum value of Vp (Rob to check)
- 3. Then set Max Vp/Vs= 3, if lower Vp to maintain the max of 3 ratio
Lame parameters (mu and lambda)
Use mu and Lambda parameters to fix Vp/Vs issues in the CVM where necessary, as part of the mesh generation. Need to make sure that patches are physical and not only to make the codes run.
- Rob to check if raw model produces lambda of <=zero.
- Note: lambda of zero corresponds to Vp/Vs= sqrt(2)=1.45
- Note: typical Vp speeds are 330 m/s in air, 1450 m/s in water and about 5000 m/s in granite
Anelastic attenuation (Q)
Frequency-independent Q definition
- Qs=100*Vs(km/s)
- Qp=2*Qs
Upper frequency
- 5Hz: based on 500 m/s and 20 m spacing
- Kim and Rob have been low-pass filtering the slip function at 5 Hz
- Ricardo and Naeem don’t filter the source
Source Models
We discussed various sources, both point source (PS) and finite-fault (FF) in the past for our verifications and validation.
- Various models were run in 2016:
- Also ran various smoothed variations of the GP15 FF (need to find a reference for that)
- Rob to follow-up with Wei (last communication was in April 2017).
Proposed models for verification and validation:
- PS: model from En-Jui (moment tensor). Can we replicate En-Jui's results? Would need to use CVMS4.26 (cms5) - not desirable. We could perform a check with CVS4.26.M01 at 5 sec. Need to use the same record En-Jui used. Proposed to use records from Figure 6 in:
- Lee E.-J., P. Chen, and T.H. Jordan (2014). Testing Waveform Predictions of 3D Velocity Models against Two Recent Los Angeles Earthquakes. Seismol. Res. Lett., 85 (6): 1275–1284 (Paper here)
- FF: Summary of agreed-upon FF model (gp.5.3.02)
- FF refined with velocity model Slides from Graves presented on Feb. 21 2018
Site Selection
Simulation results
Description: results for the small region from November 10 2017:
Description: small region simulations for Hercules, AWP (updated since November 2017) and RWG (medium size model to avoid boundary reflexions):
Description: results for the small region from November 10 2017:
Description: small region simulations for Hercules (updated to use simulations including Q), AWP (updated since November 2017) and RWG (medium size model to avoid boundary reflexions):
Latest 2018 results
Description: results for the small and large regions (model dependent) from November 10 2017, computed 20181125 using the latest ts_processor codes:
Summary documentation of extended source simulations submitted at high-f:
Hercules 20171106 notes:
- Simulation of Mw 5.14 Lahabra Earthquake.(2014/03/29 04:09:42.97) with damping and using the extended source model provided by Robert Graves (downsampled version of the SRF for La Habra and based on 3D velocity model.The subfault spacing is 100 m). This simulation is part of High-F project. The simulation is done using the CVM-etree extracted from CVM-S4.26.M01 (28*28*14 km) at the “Small Domain." See readme.txt in their folder for more info.
RWG 20180425 notes:
- 20180425 - RWG simulation results for full-size model region (135 km x 180 km x 60 km). Velocity mesh was created on BlueWaters using UCVM interface that Scott C. created. Model is CVM-S4.26.M01 (at least I believe so, I selected 'cvmsi' as input to single_exe.py). Vp/Vs capped at 3.0, Vsmin=500 ms, Qs=100*Vs(m/s), Qp=2*Qs. Time step is dt=0.001 sec. Waveforms are ground velocity and have been low-pass filtered using zero-phase 4th order butterworth with corner at 5 Hz.
AWP 20181203 no notes:
- no documentation on submission. Latest submission with info was for 20171204 submission, per fabio: rerun of latest anelastic simulation using cap of vp/vs=3 with increased simulation length to 60 seconds
2019-12 summary of plans for 2020
- Rob to test other finite fault models using near-fault data
- Christine to coordinate Jan, sprint for Fabio to allow use of BBP GOF post-processors with high-f results and data and to revive GOF map plotting tools. Need to verify that PSA results are ok (see anomaly on 20181125 results for observations of CE_13066)
- Christine to re-start verification paper development
- AWP and RWG to document what they submitted so far or decide to rerun simulations. Hercules has clear documentation already.
Related Pages