UCVM 14.3.0

Introduction

Unified Community Velocity Model (UCVM) software framework is a collection of software tools designed to provide standard interface to multiple, alternative, California 3D velocity models. One important use of UCVM is in high resolution 3D wave propagation simulations for California. UCVM development is an interdisciplinary research collaboration involving geoscientists and computer scientists. UCVM geoscience research includes identification and assembly of existing California velocity models into a state-wide model and improvements to existing velocity models. UCVM computer science research includes definition of a easy-to-use CVM query interface, integration of regional 3D and geotechnical models, and automated CVM evaluation processing capabilities.

Download

Platform	File	Download	Mirror
Linux	SCEC UCVM 14.3.0 Official Release (391 Mb)	ucvm-14.3.0.tar.gz	N/A
Linux	SCEC UCVM 14.3.0 md5 checksum (< 1Kb)	ucvm-14.3.0.tar.gz.md5	N/A

User Guide

We strongly recommend that you read through the UCVM 14.3.0 user guide for more information on how to install and use UCVM.

Developer's Guide

Learn how to integrate your own velocity models and add other functionality to UCVM by visiting the UCVM 14.3.0 developer guide.

Issue Tracking System

Users can report issues, bugs, and request features using the UCVM Trac system. Requires a SCEC login.

• UCVM Trac System

Known Issues

Illegal Instruction When Running on Titan Head Node

When trying to run UCVM on a head node in Titan, an illegal instruction is thrown. Please note that UCVM does work correctly on the aprun nodes. Because of the fact that the collection of utilities do not run on the head node, doing "make check" at the end of the install will not work. We strongly recommend testing your installation via the following PBS script (please modify for your account and path locations):

#!/bin/bash
#Begin PBS directives
#PBS -A geo015
#PBS -N make_check
#PBS -j oe
#PBS -l walltime=1:00:00,nodes=1
#End PBS directives and begin shell commands
cd /lustre/atlas/proj-shared/geo015/gill/UCVM/test
aprun -n 1 ./accepttest
aprun -n 1 ./unittest
aprun -n 1 ./accept_test.py /lustre/atlas/proj-shared/geo015/gill/ucvm-14.3.0

Using GNU Compilers on Mac OS X

If you have installed the GNU version of GCC on your Mac, you may run into issues installing UCVM. Typically, these versions are installed either with Fink or MacPorts. We are currently investigating this issue and hope for it to be fixed in the next release. In the meantime, please remove the Fink path (/sw) or the MacPorts path from your $PATH so that way UCVM will compile with the default OS X utilities. Please note that you will need the Xcode developer tools installed for this to work.

Supported Model List

UCVM can download and install the following 3D velocity models:

• CVM-S4
• CVM-H 11.9.1
• CVM-S4.26
• CenCal 0.8.0

The following models require manual download, installation, and additional configuration within the UCVM platform.

• SCEC CVM-NCI (requires manual download and installation)
• Magistrale Wasatch Front CVM
• Graves Cape Mendocino
• Lin-Thurber Statewide
• Tape SoCal

Additional Functionality

The following is the coverage region and projection for the included 2D elevation and Vs30 maps. The entire state of California is included, along with portions of Oregon, Nevada, Arizona, and northern Mexico.

Coverage region for UCVM 2D maps (cyan) overlayed upon regions of various California 3D velocity models (CVM-S: red, CVM-H: blue, LinThurber: yellow, Thurber NoCal: green, USGS Bay Area: white)

DEM

A statewide DEM for the proposed coverage region is included within UCVM. Elevation data has been sampled from USGS NED 1 arcsec dataset (~30 m), and bathymetric data from the NOAA ETOPO1 1' dataset (~1.5 km). This DEM is currently sampled at a resolution of approximately 220m but this may be increased. The elevation data is stored as a fixed resolution Etree for the entire 1800 x 900 km region. The elevation at a particular point is smoothed using bilinear interpolation of the surrounding four elevation octants.

UCVM Digital Elevation Model, with elevation set to 0.0 m outside of coverage region. The distortion is caused by differences in the DEM projection (Az-equi) and the plot projection (cylindrical).

Vs30 Maps

Two statewide Vs30 maps for the proposed coverage region are included within UCVM:

• Wills-Wald Vs30 map (default): Vs30 data for the California landmass has been sampled from the Wills (2006) dataset at approx 0.0002197 D resolution, and out-of-state/ocean areas have been sampled from the Wald (2007) dataset at 0.0083333 D resolution.
• Yong-Wald Vs30 map (optional): Vs30 data for the California landmass has been sampled from the Yong (2011) dataset at approx 0.013 D resolution, and out-of-state/ocean areas have been sampled from the Wald (2007) dataset at 0.0083333 D resolution.

These maps are currently sampled at a resolution of approximately 220m but this may be increased. The Vs30 data is stored as a fixed resolution Etree for the entire 1800 x 900 km region. The Vs30 value at a particular point is smoothed using bilinear interpolation of the surrounding four map octants.

UCVM Wills-Wald Vs30 map, with Vs30 set to 0.0 km/s outside of coverage region.

Geotechnical Layer

A statewide GTL for the proposed coverage region is included within UCVM. This is based on the Vs30-derived GTL method developed by Ely (2010). Interpolation between this GTL and the underlying crustal models is accomplished with the interpolation method described in that same publication. The z range over which interpolation is performed is configurable, as is the interpolation method used. The Ely method takes its input Vs30 value from the Vs30 map included within UCVM (described in the previous section).

UCVM Vs map at depth 0.0 m, showing the Ely GTL overlayed on the 1D background.

Additionally, UCVM is able to support any number of other user-defined GTLs. It will combine them in a manner analogous to how the crustal models are combined. Each GTL may also have its own user-defined interpolation function that is used to blend it with the underlying crustal models. If no interpolation method is specified, linear interpolation is assumed.

These shows how velocity model material properties delivered by a 3D velocity model depend on how the velocity model is queried and whether the model considers elevation, or whether the earth is modeled as a flat region.

Ability to Combine Regional Models

The UCVM API currently combines models by simply tiling them in a user-specified order. The software also provides a post-processing "patching" mechanism to smooth out discontinuities at the interfaces between different regional models. The user identifies a 3D geographic region in the form of an open box where smoothing is to be performed. A patcher is run to extract the material properties along the surfaces of this box at a fixed resolution (eg: 250 m) and it saves them into a new model which may then be imported into and queried by UCVM. Points that fall within the smoothing region are interpolated using an inverse distance weighting function with the 10% closest surface points on the same z-plane. A power factor of 2.0 is used for weighting. Linear interpolation is used along the z-axis to smooth out the edge values.

The following series of images show how this process works:

Fig 1: Unsmoothed map at depth 0m (left). Wire frame outlines of smoothing regions (center). Final smoothed result (right).

Exporting to Etree

Any number of regional models and GTLs may be combined into a single meta-model and then exported to an Etree with variable resolution. The local resolution at a particular point is determined with the following relation:

local rez <= local Vs / (fmax * ppw)

where:
local rez : the local resolution in meters
local Vs  : the local minimum Vs in meters/sec
fmax      : desired max frequency to support in Hz
ppw       : points per wavelength

UCVM Background Material

Numerous resources are available which help explain more about UCVM and the history of its development.

• UCVM Poster From SCEC Annual Meeting 2013
• UCVM Poster From CIG-QUEST Workshop 2013
• UCVM Talk From CIG-QUEST Workshop 2013

Learn more about each CVM as well:

CVM-S

• CVM-S4 Wiki Page
• Collection of Presentations about CVM-S and CVM-H (56Mb ppt file)

CVM-H

• CVM-H Wiki Page
• Harvard Structural Geology & Earth Resources Group
• Collection of Presentations about CVM-S and CVM-H (56Mb ppt file)

Lin-Thurber CVM

• GuoQing Lin's Research

Wasatch Front CVM

• Wasatch Front CVM Home Page

Other

• Qinya Liu Tomography Research

References

1. Ely, G., T. H. Jordan, P. Small, P. J. Maechling (2010), A Vs30-derived Near-surface Seismic Velocity Model Abstract S51A-1907, presented at 2010 Fall Meeting, AGU, San Francisco, Calif., 13-17 Dec. [Ely2010-AGU.pdf]
2. Graves, R. (1994), Rupture History and Strong Motion Modeling of the 1992 Cape Mendocino Earthquake, USGS External Grant Report
3. Lin, G., C. H. Thurber, H. Zhang, E. Hauksson, P. Shearer, F. Waldhauser, T. M. Brocher, and J. Hardebeck (2010), A California statewide three-dimensional seismic velocity model from both absolute and differential Times, Bull. Seism. Soc. Am., 100, in press. supplemental
4. Taborda R., López J., O'Hallaron D., Tu T. and Bielak J. (2007), A review of the current approach to CVM-Etrees, SCEC Annual Meeting, Palm Springs, CA, USA, September 8–12. [1]
5. Thurber, C., H. Zhang, F. Waldhauser, J. Hardebeck, A. Michael, and D. Eberhart-Phillips (2006), Three-dimensional compressional wavespeed model, earthquake relocations, and focal mechanisms for the Parkfield, California, region, Bull. Seism. Soc. Am., 96, S38-S49, 2006. supplemental
6. Thurber, C., T. Brocher, H. Zhang, and V. Langenheim (2007), Three-dimensional P-wave velocity model for the San Francisco Bay region, California, J. Geophys. Res., 112, B07313, doi:10.1029/ 2006JB004682, 2007.
7. Thurber, C., H. Zhang, T. Brocher, and V. Langenheim (2009), Regional three-dimensional seismic velocity model of the crust and uppermost mantle of northern California, J.Geophys. Res., 114, B01304, doi:10.1029/2008JB005766, 2009.
8. Wald, D. J., and T. I. Allen (2007), Topographic slope as a proxy for seismic site conditions and amplification, Bull. Seism. Soc. Am., 97 (5), 1379-1395, doi:10.1785/0120060267.
9. Wills, C. J., and K. B. Clahan (2006), Developing a map of geologically defined site-condition categories for California, Bull. Seism. Soc. Am., 96 (4A), 1483-1501, doi:10.1785/0120050179.
10. Yong, A., Hough, S.E., Iwahashi, J., and A. Braverman (2012), A terrain-based site conditions map of California with implications for the contiguous United States, Bull. Seism. Soc. Am., Vol. 102, No. 1, pp. 114–128, February 2012, doi: 10.1785/0120100262.