by Geoffrey Ely, Patrick Small, Thomas Jordan, Philip Maechling, and Feng Wang
Shallow material properties, S-wave velocity in particular, strongly influence ground motions, so must be accurately characterized for ground-motion simulations. Available near-surface velocity information generally exceeds that which is accommodated by crustal velocity models, such as current versions of the SCEC Community Velocity Model (CVM-S v4.0 by Magistrale, McLaughlin, and Day 1996; Magistrale et al. 2000; Kohler, Magistrale, and Clayton 2003) or the Harvard model (CVM-H v6.3 by Shaw et al. 2015). The elevation-referenced CVM-H voxel model introduces rasterization artifacts in the near-surface due to course sample spacing, and sample depth dependence on local topographic elevation. To address these issues, we propose a method to supplement crustal velocity models, in the upper few hundred meters, with a model derived from available maps of VS30 (the average S-wave velocity down to 30 meters). The method is universally applicable to regions without direct measures of VS30 by using VS30 estimates from topographic slope (Wald and Allen 2007). In our current implementation for Southern California, the geology-based VS30 map of Wills and Clahan (2006) is used within California, and topography-estimated VS30 is used outside of California.
Various formulations for S-wave velocity depth dependence, such as linear spline and polynomial interpolation, were evaluated against the following priorities: (a) capability to represent a wide range of soil and rock velocity profile types; (b) smooth transition to the crustal velocity model; (c) ability to reasonably handle poor spatial correlation of VS30 and crustal velocity data; (d) simplicity and minimal parameterization; and (e) computational efficiency. The favored model includes cubic and square-root depth dependence, with the model extending to a transition depth zT. A transition depth of zT = 350 m is used to ensure adequate sampling of CVM-H (shallower depths may be unsampled by the CVM-H near topographic features). S-wave velocity at the surface is derived from VS30 by a uniform scaling. VP, and in turn density, are inferred from surface VS via the scaling laws of Brocher (2005). VS and VP are independently interpolated between the surface values and those extracted from the crustal velocity model at the transition depth. Density is derived from interpolated VP via the Nafe-Drake law of Brocher. Depth dependence for the interpolation is parameterized with
|z||= zʹ / zT|
|f(z)||= z + b(z - z2)|
|g(z)||= a - az + c(z2 + 2√z - 3z)|
|VS(z)||= f(z)VST + g(z)VS30|
|VP(z)||= f(z)VPT + g(z)P(VS30)|
where zʹ is depth, VST and VPT are S- and P-wave velocities extracted from the crustal velocity model at depth zT, P() is the Brocher P-wave velocity scaling law, and R() is the Nafe-Drake law. The coefficient a controls the ratio of surface velocity to original 30 meter average, b controls overall curvature, and c controls near-surface curvature.
The coefficients a = 1/2, b = 2/3, and c = 3/2 were chosen by trial-and-error fitting Boore and Joyner’s (1997) generic rock profile and CVM-S generic soil profiles, as well as to produce smooth and well-behaved profiles when applied to the CVM-H at the selected CyberShake sites.
Boore, D., J. Gibbs, and M. Rodriguez. 2003. “A Compendium of P- and S- Velocities from Surface-to-Borehole Logging: Summary and Reanalysis of Previously Published Data and Analysis of Unpublished Data.” Open-File Report 03-191. U.S. Geological Survey. http://daveboore.com/data_online.html.
Boore, D., and W. Joyner. 1997. “Site Amplifications for Generic Rock Sites.” Bull. Seism. Soc. Am. 87 (2): 327–41.
Brocher, T. 2005. “Empirical Relations Between Elastic Wavespeeds and Density in the Earth’s Crust.” Bull. Seism. Soc. Am. 95 (6): 2081–92. https://doi.org/10.1785/0120050077.
Kohler, M., H. Magistrale, and R. Clayton. 2003. “Mantle Heterogeneities and the SCEC Reference Three-Dimensional Seismic Velocity Model Version 3.” Bull. Seism. Soc. Am. 93 (2): 757–74. https://doi.org/10.1785/0120020017.
Magistrale, H., S. Day, R. Clayton, and R. Graves. 2000. “The SCEC Southern California Reference Three-Dimensional Seismic Velocity Model Version 2.” Bull. Seism. Soc. Am. 90 (6B): S65–76. https://doi.org/10.1785/0120000510.
Magistrale, H., K. McLaughlin, and S. Day. 1996. “A Geology-Based 3D Velocity Model of the Los Angeles Basin Sediments.” Bull. Seism. Soc. Am. 86 (4): 1161–6.
Shaw, J., A. Plesch, C. Tape, M. Suess, T. Jordan, G. Ely, E. Hauksson, et al. 2015. “Unified Structural Representation of the Southern California Crust and Upper Mantle.” Earth Planet. Sc. Lett. 415 (0): 1–15. https://doi.org/10.1016/j.epsl.2015.01.016.
Wald, D., and T. Allen. 2007. “Topographic Slope as a Proxy for Seismic Site Conditions and Amplification.” Bull. Seism. Soc. Am. 97 (5): 1379–95. https://doi.org/10.1785/0120060267.
Wills, C., and K. Clahan. 2006. “Developing a Map of Geologically Defined Site-Condition Categories for California.” Bull. Seism. Soc. Am. 96 (4A): 1483–1501. https://doi.org/10.1785/0120050179.