Applying a three-dimensional velocity model, waveform cross correlation, and cluster analysis to locate southern California seismicity from 1981 to 2005

We compute high-precision earthquake locations using southern California pick and waveform data from 1981 to 2005. Our latest results are significantly improved compared to our previous catalog by the following: (1) We locate events with respect to a new crustal P and S velocity model using three-dimensional ray tracing, (2) we examine six more years of waveform data and compute cross-correlation results for many more pairs than our last analysis, and (3) we compute locations within similar event clusters using a new method that applies a robust fitting method to obtain the best locations satisfying all the differential time constraints from the waveform cross correlation. These results build on the relocated catalogs of Hauksson and Shearer (2005) and Shearer et al. (2005) and provide additional insight regarding the fine-scale fault structure in southern California and the relationship between the San Andreas Fault (SAF) and nearby seismicity. In particular, we present results for two regions in which the seismicity near the southern SAF seems to align on dipping faults

Time-dependent seismic tomography of the Coso geothermal area, 1996-2004

Local-earthquake tomographic images were calculated for each of the years 1996 - 2004 using arrival times from the U.S. Navy’s permanent seismometer network at the Coso geothermal area, California. The results show irregular strengthening with time of the wave-speed ratio VP/VS at shallow depths. These changes result predominately from progressive relative increase in VS with respect to VP, and could result from processes associated with geothermal operations such as decrease in fluid pressure and the drying of argillaceous minerals such as illite

A three-dimensional crustal seismic velocity model for southern California from a composite event method

We present a new crustal seismic velocity model for southern California derived from P and S arrival times from local earthquakes and explosions. To reduce the volume of data and ensure a more uniform source distribution, we compute "composite event" picks for 2597 distributed master events that include pick information for other events within spheres of 2 km radius. The approach reduces random picking error and maximizes the number of S wave picks. To constrain absolute event locations and shallow velocity structure, we also use times from controlled sources, including both refraction shots and quarries. We implement the SIMULPS tomography algorithm to obtain three-dimensional (3-D) V_p and V_p/V_s structure and hypocenter locations of the composite events. Our new velocity model in general agrees with previous studies, resolving low-velocity features at shallow depths in the basins and some high-velocity features in the midcrust. Using our velocity model and 3-D ray tracing, we relocate about 450,000 earthquakes from 1981 to 2005. We observe a weak correlation between seismic velocities and earthquake occurrence, with shallow earthquakes mostly occurring in high P velocity regions and midcrustal earthquakes occurring in low P velocity regions. In addition, most seismicity occurs in regions with relatively low V_p/V_s ratios, although aftershock sequences following large earthquakes are often an exception to this pattern

A physics-based earthquake simulator replicates seismic hazard statistics across California

Seismic hazard models are important for society, feeding into building codes and hazard mitigation efforts. These models, however, rest on many uncertain assumptions and are difficult to test observationally because of the long recurrence times of large earthquakes. Physics-based earthquake simulators offer a potentially helpful tool, but they face a vast range of fundamental scientific uncertainties. We compare a physics-based earthquake simulator against the latest seismic hazard model for California. Using only uniform parameters in the simulator, we find strikingly good agreement of the long-term shaking hazard compared with the California model. This ability to replicate statistically based seismic hazard estimates by a physics-based model cross-validates standard methods and provides a new alternative approach needing fewer inputs and assumptions for estimating hazard

Delayed Seismicity Rate Changes Controlled by Static Stress Transfer

On 15 June 2010, a M_w5.7 earthquake occurred near Ocotillo, California, in the Yuha Desert. This event was the largest aftershock of the 4 April 2010 M_w7.2 El Mayor-Cucapah (EMC) earthquake in this region. The EMC mainshock and subsequent Ocotillo aftershock provide an opportunity to test the Coulomb failure hypothesis (CFS). We explore the spatiotemporal correlation between seismicity rate changes and regions of positive and negative CFS change imparted by the Ocotillo event. Based on simple CFS calculations we divide the Yuha Desert into three subregions, one triggering zone and two stress shadow zones. We find the nominal triggering zone displays immediate triggering, one stress shadowed region experiences immediate quiescence, and the other nominal stress shadow undergoes an immediate rate increase followed by a delayed shutdown. We quantitatively model the spatiotemporal variation of earthquake rates by combining calculations of CFS change with the rate-state earthquake rate formulation of Dieterich (1994), assuming that each subregion contains a mixture of nucleation sources that experienced a CFS change of differing signs. Our modeling reproduces the observations, including the observed delay in the stress shadow effect in the third region following the Ocotillo aftershock. The delayed shadow effect occurs because of intrinsic differences in the amplitude of the rate response to positive and negative stress changes and the time constants for return to background rates for the two populations. We find that rate-state models of time-dependent earthquake rates are in good agreement with the observed rates and thus explain the complex spatiotemporal patterns of seismicity

Delayed Seismicity Rate Changes Controlled by Static Stress Transfer

On 15 June 2010, a M_w5.7 earthquake occurred near Ocotillo, California, in the Yuha Desert. This event was the largest aftershock of the 4 April 2010 M_w7.2 El Mayor-Cucapah (EMC) earthquake in this region. The EMC mainshock and subsequent Ocotillo aftershock provide an opportunity to test the Coulomb failure hypothesis (CFS). We explore the spatiotemporal correlation between seismicity rate changes and regions of positive and negative CFS change imparted by the Ocotillo event. Based on simple CFS calculations we divide the Yuha Desert into three subregions, one triggering zone and two stress shadow zones. We find the nominal triggering zone displays immediate triggering, one stress shadowed region experiences immediate quiescence, and the other nominal stress shadow undergoes an immediate rate increase followed by a delayed shutdown. We quantitatively model the spatiotemporal variation of earthquake rates by combining calculations of CFS change with the rate-state earthquake rate formulation of Dieterich (1994), assuming that each subregion contains a mixture of nucleation sources that experienced a CFS change of differing signs. Our modeling reproduces the observations, including the observed delay in the stress shadow effect in the third region following the Ocotillo aftershock. The delayed shadow effect occurs because of intrinsic differences in the amplitude of the rate response to positive and negative stress changes and the time constants for return to background rates for the two populations. We find that rate-state models of time-dependent earthquake rates are in good agreement with the observed rates and thus explain the complex spatiotemporal patterns of seismicity

Regional and teleseismic double-difference earthquake relocation using waveform cross-correlation and global bulletin data

We have developed a double-difference algorithm to relocate earthquakes recorded at global seismic networks, using differential arrival times for first and later arriving regional and global phases to invert for the vectors connecting the hypocenters. Differential times are formed from global seismic bulletins and are accurately measured on similar seismograms by time domain waveform cross correlation. We evaluate the performance of this spherical, multiphase double-difference algorithm using three-dimensional regional-scale synthetic data and two sets of earthquake data in different tectonic settings. The first includes 3783 intermediate depth earthquakes that occurred between 1964 and 2000 in the subducting Nazca plate beneath northern Chile, where the relocated seismicity confirms a narrowly spaced double seismic zone previously imaged with temporary local seismic data. Residual statistics and comparison with accurately known locations indicate mean relative location errors at the 90% confidence level of 2.4 km laterally and 1.8 km vertically. Later events typically constrained by cross-correlation data have errors of 1.6 km laterally and 1.4 km vertically. The second data set includes 75 crustal earthquakes in the 1999 Izmit and Düzce, Turkey, aftershock sequences, where the double-difference solutions image orientation and dip of individual fault segments that are consistent with focal mechanisms and near-surface information. Fault complexity likely causes a low level of waveform similarity in this aftershock sequence and thus generates fewer correlated events compared to the Chile earthquakes. Differences between the double-difference locations and corresponding locations in global seismicity catalogs (Earthquake Data Report, EDR; International Seismological Centre, ISC; Engdahl-Hilst-Buland, EHB) are typically greater than 10 km. We evaluate the potential of cross-correlation and double-difference methods to improve hypocenter locations on a global scale