Analysis of surface and seismic sources in dense array data with match field processing and Markov chain Monte Carlo sampling

SUMMARY We introduce a methodology based on array processing to detect and locate weak seismic events in a complex fault zone environment. The method is illustrated using data recorded by a dense array of 1108 vertical component geophones in a 600 m × 600 m area on the Clark branch of the San Jacinto Fault. Because surface and atmospheric sources affect weak ground motion, it is necessary to discriminate them from weak seismic sources at depth. Source epicentral positions and associated apparent velocities are extracted from continuous seismic waveforms using Match Field Processing (MFP). We implement MFP at specific frequencies targeting surface and subsurface sources, using for computational efficiency a forward model of acoustic source in a homogenous medium and Markov Chain Monte Carlo sampling. Surface sources such as Betsy gun shots and a moving vehicle are successfully located. Weak seismic events are also detected outside of the array, and their backazimuth angle is retrieved and found to be consistent with the fault geometry. We also show that the homogeneous acoustic model does not yield satisfying results when extracting microseismic event depth, because of the ambiguity between depth and the apparent velocity based on surface data.