skip to main content

Title: Parsimonious Velocity Inversion Applied to the Los Angeles Basin, CA

The proliferation of dense arrays promises to improve our ability to image geological structures at the scales necessary for accurate assessment of seismic hazard. However, combining the resulting local high‐resolution tomography with existing regional models presents an ongoing challenge. We developed a framework based on the level‐set method that infers where local data provide meaningful constraints beyond those found in regional models ‐ for example the Community Velocity Models (CVMs) of southern California. This technique defines a volume within which updates are made to a reference CVM, with the boundary of the volume being part of the inversion rather than explicitly defined. By penalizing the complexity of the boundary, a minimal update that sufficiently explains the data is achieved. To test this framework, we use data from the Community Seismic Network, a dense permanent urban deployment. We inverted Love wave dispersion and amplification data, from the Mw 6.4 and 7.1 2019 Ridgecrest earthquakes. We invert for an update to CVM‐S4.26 using the Tikhonov Ensemble Sampling scheme, a highly efficient derivative‐free approximate Bayesian method. We find the data are best explained by a deepening of the Los Angeles Basin with its deepest part south of downtown Los Angeles, along with a steeper northeastern basin wall. This result offers new progress toward the parsimonious incorporation of detailed local basin models within regional reference models utilizing an objective framework and highlights the importance of accurate basin models when accounting for the amplification of surface waves in the high‐rise building response band.

more » « less
Award ID(s):
2011079 2011107
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Medium: X
Sponsoring Org:
National Science Foundation
More Like this

    We use broad-band stations of the ‘Los Angeles Syncline Seismic Interferometry Experiment’ (LASSIE) to perform a joint inversion of the Horizontal to Vertical spectral ratios (H/V) and multimode dispersion curves (phase and group velocity) for both Rayleigh and Love waves at each station of a dense line of sensors. The H/V of the autocorrelated signal at a seismic station is proportional to the ratio of the imaginary parts of the Green’s function. The presence of low-frequency peaks (∼0.2 Hz) in H/V allows us to constrain the structure of the basin with high confidence to a depth of 6 km. The velocity models we obtain are broadly consistent with the SCEC CVM-H community model and agree well with known geological features. Because our approach differs substantially from previous modelling of crustal velocities in southern California, this research validates both the utility of the diffuse field H/V measurements for deep structural characterization and the predictive value of the CVM-H community velocity model in the Los Angeles region. We also analyse a lower frequency peak (∼0.03 Hz) in H/V and suggest it could be the signature of the Moho. Finally, we show that the independent comparison of the H and V components with their corresponding theoretical counterparts gives information about the degree of diffusivity of the ambient seismic field.

    more » « less
  2. Abstract

    In this study, we investigate the shallow shear wave velocity structure of the Los Angeles Basin in southern California, using ambient noise correlations between 5 dense arrays and 21 broadband stations from the Southern California Seismic Network (SCSN). We observe clear fundamental mode and first overtone Rayleigh waves in the frequency band 0.25–2.0 Hz, and obtain group velocity maps through tomography. We further derive a 3D shear wave velocity model, covering a large portion of the central LA Basin for the depths shallower than 3 km. We found that the small scale shallow velocity structure heterogeneities are better resolved compared with the SCEC Community velocity models. Our model captures the presence of the Newport‐Inglewood fault by a NW–SE trending high velocity belt. Our model provides more accurate constraints on local ground motion predictions with detailed mapping of structural heterogeneities.

    more » « less

    The near-surface seismic structure (to a depth of about 1000 m), particularly the shear wave velocity (VS), can strongly affect the propagation of seismic waves and, therefore, must be accurately calibrated for ground motion simulations and seismic hazard assessment. The VS of the top (<300 m) crust is often well characterized from borehole studies, geotechnical measurements, and water and oil wells, while the velocities of the material deeper than about 1000 m are typically determined by tomography studies. However, in depth ranges lacking information on shallow lithological stratification, typically rock sites outside the sedimentary basins, the material parameters between these two regions are typically poorly characterized due to resolution limits of seismic tomography. When the alluded geological constraints are not available, models, such as the Southern California Earthquake Center (SCEC) Community Velocity Models (CVMs), default to regional tomographic estimates that do not resolve the uppermost VS values, and therefore deliver unrealistically high shallow VS estimates. The SCEC Unified Community Velocity Model (UCVM) software includes a method to incorporate the near-surface earth structure by applying a generic overlay based on measurements of time-averaged VS in top 30 m (VS30) to taper the upper part of the model to merge with tomography at a depth of 350 m, which can be applied to any of the velocity models accessible through UCVM. However, our 3-D simulations of the 2014 Mw 5.1 La Habra earthquake in the Los Angeles area using the CVM-S4.26.M01 model significantly underpredict low-frequency (<1 Hz) ground motions at sites where the material properties in the top 350 m are significantly modified by the generic overlay (‘taper’). On the other hand, extending the VS30-based taper of the shallow velocities down to a depth of about 1000 m improves the fit between our synthetics and seismic data at those sites, without compromising the fit at well-constrained sites. We explore various tapering depths, demonstrating increasing amplification as the tapering depth increases, and the model with 1000 m tapering depth yields overall favourable results. Effects of varying anelastic attenuation are small compared to effects of velocity tapering and do not significantly bias the estimated tapering depth. Although a uniform tapering depth is adopted in the models, we observe some spatial variabilities that may further improve our method.

    more » « less

    Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

    more » « less
  5. We constrained sedimentary basin structure using a nodal seismic array consisting of ten dense lines that overlie multiple basins in the northern Los Angeles area. The dense array consists of 758 seismic nodes, spaced ~250–300 m apart along linear transects, that recorded ground motions for 30–35 days. We applied the receiver function (RF) technique to 16 teleseismic events to investigate basin structure. Primary basin-converted phases were identified in the RFs. A shear wave velocity model produced in a separate study using the same dataset was incorporated to convert the basin time arrivals to depth. The deepest part of the San Bernardino basin was identified near the Loma Linda fault at a depth of 2.4 km. Basin depths identified at pierce points for separate events reveal lateral changes in basin depth across distances of ~2–3 km near individual stations. A significant change in basin depth was identified within a small distance of ~4 km near the San Jacinto fault. The San Gabriel basin exhibited the largest basin depths of all three basins, with a maximum depth of 4.2 km. The high lateral resolution from the dense array helped to reveal more continuous structures and reduce uncertainties in the RFs interpretation. We discovered a more complex basin structure than previously identified. Our findings show that the basins’ core areas are not the deepest, and significant changes in basin depth were observed near some faults, including the San Jacinto fault, Fontana fault, Red Hill fault and Indian Hill fault.

    more » « less