Abstract:Our research program consists of four components, each involving some aspect of multiple-event analysis: (1) highprecision waveform cross-correlation (WCC) for arrival time estimation, (2) robust event clustering, (3) waveform decomposition and source wavelet deconvolution reshaping, (4) double-difference (DD) multiple-event location and tomography. Our research is focusing initially on the development and testing of these seismic analysis methods using "ground-truth" datasets at different scales (local and regional), and will be followed by the application of these methods to "test-bed" seismic datasets, including tests on real-time data streams. We are constructing two regional ground-truth datasets, one for New Mexico and one for California. The former dataset is from the New Mexico Tech short-period network and is dominated by small explosions, whereas the latter is from the UC-Berkeley broadband network and contains relatively large earthquakes (about 260 events of magnitude 4.3 and above). In each case, the ground-truth dataset will be used initially for algorithm assessment purposes and later as a "reference event" dataset for clustering and locating new events. The New Mexico "training" dataset is undergoing final grooming. For the California dataset, we have affirmed the quality of the ground-truth location information by carrying out a DD location analysis using catalog time-difference data. We also present initial results of hierarchical cross-correlation analysis applied to the California dataset. We have completed the modifications to the integer portion of the cross-correlation software, adding the multiple eigentaper approach to provide improved integer correlation statistics. The modifications to the program logic provide the initial step for subsequent adaptation towards multiple frequency-scale correlations. In addition, modifications have been completed to interface the WCC software with our new DD tomography code (discussed below). We are preparing the California dataset for further DD location analysis and for processing with the regional-scale tomography code under development. At the same time, we will experiment with clustering and waveform decomposition and reshaping techniques to expand the ability to use WCC on large datasets. Efforts have also begun on an Earthworm-MatSeis real-time interface, including the generation of CSS database files from Earthworm. We have developed a new DD tomography method that uses high-precision WCC relative arrival time picks as well as differential and absolute catalog arrival-time picks. "Double-difference" refers to the residual between the observed differential arrival-time picks and the predicted differential arrival-times. If two sources have almost identical source mechanisms and the differential signal scattering due to velocity heterogeneity along the ray path is minor, they will produce similar waveforms at the same station, allowing for the use of WCC to determine precise differential times. Compared to absolute travel-time picks, relative arrival-time picks obtained using WCC methods have a precision to the sub-sample level. When two sources are far away from each other so that their waveforms are not similar, we just use the relative catalog arrival-time pick differences. We also use the absolute travel-time picks to constrain the absolute event locations. Our new DD tomography method has been applied to the Hayward fault, California, as an initial test to compare to previously published location results. In comparison to the standard seismic tomography method, the doubledifference tomography method results in a sharper-featured velocity structure that is consistent with the geologic setting. At the same time, it produces absolute and relative event locations accurately, similar to the results obtained by the DD location method developed by Waldhauser and Ellsworth (2000). Our next steps are to apply this algorithm to another test dataset, from Parkfield, California, and then to develop the regional-scale version and apply it to our regional California ground-truth dataset. 438 24th Seismic Research Review – Nuclear Explosion Monitoring: Innovation and Integration

New constraints on seismicity in the Wellington region of New Zealand from relocated earthquake hypo

Triple Junction Kinematics Accounts for the 2016 Mw 7.8 Kaikoura Earthquake Rupture Complexity.

The 2016 Kaikōura Earthquake: Simultaneous Rupture of the Subduction Interface and Overlying Faults

Space Imaging Geodesy Reveals Near Circular, Coseismic Block Rotation During the 2016 Mw 7.8 Kaikōura Earthquake, New Zealand

Tests of the precursory accelerating moment release model using a synthetic seismicity model for Wellington, New Zealand

Analysis of site-response residuals from empirical ground-motion models to account for observed sedimentary basin effects in Wellington, New Zealand

The Magnitude–Frequency Distributions of Earthquakes in Aotearoa New Zealand and on Adjoining Subduction Zones, Using a New Integrated Earthquake Catalog

A New Method for High‐resolution Event Relocation and Application to the Aftershocks of Lushan Earthquake, China

Inverting Geodetic Strain Rates for Slip Deficit Rate in Complex Deforming Zones: An Application to the New Zealand Plate Boundary

Mapping Depth to Bedrock, Shear Stiffness, and Fundamental Site Period at CentrePort, Wellington using Surface Wave Methods: Implications for Local Seismic Site Amplification

Late quaternary activity of the NW Cardrona Fault, Otago, New Zealand

Effective elastic thickness of Zealandia and its implications for lithospheric deformation

Spatiotemporal Evolution of Slow Slip Events at the Offshore Hikurangi Subduction Zone in 2019 Using GNSS, InSAR, and Seafloor Geodetic Data

Structural Heterogeneity of the Midcrust Adjacent to the Central Alpine Fault, New Zealand: Inferences from Seismic Tomography and Seismicity Between Harihari and Ross

Earthquake Relocations Delineate a Discrete Fault Network and Deformation Corridor Throughout Southeast Alaska and Southwest Yukon

th Seismic Research Review – Nuclear Explosion Monitoring : Innovation and Integration

Dominance of transient rheology in far-field postseismic deformation following the 2012 Mw 7.8 Haida Gwaii and 2013 Mw 7.5 Craig earthquakes

Reconstruction of a subsoil model for local seismic response evaluation through experimental and numerical methods: The case of the Wellington CBD, New Zealand

Double-Difference Tomography: the Method and Its Application to the Hayward Fault, California

Segmentation of Shallow Slow Slip Events at the Hikurangi Subduction Zone Explained by Along-Strike Changes in Fault Geometry and Plate Convergence Rates

Aftershocks of the 2012 Off-Coast of Sumatra Earthquake Sequence